Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Charge sharing in pixelated semiconductor sensors

The charge sharing between neighboring pixels in pixelated sensors can be used to measure particle or x-ray coordinates with accuracy better than the pixel pitch. The accurate model of the charge distribution shape is essential to achieve ultimate coordinate accuracy. The charge sharing is caused by charge carriers' diffusion on the path from the generation point to pixels. This paper is focused on the diffusion of the initially compact charge cloud in the field free region. The diffusion equation solutions are obtained using separation of variable and Fourier synthesis method for different initial conditions and resulting charge distributions are integrated over pixel areas. In conclusion, the look up table containing pre-calculated values for pixel charge fractions is proposed to speed up numerical calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

POLCA8 - modelling of cross section variations inside hexagonal assemblies

This paper presents the POLCA8 approach for modelling non-constant cross section distributions inside hexagonal fuel assemblies. The multigroup diffusion equation is modified to account for intranodal cross section variations. The obtained equation is solved in a node-wise manner based on the Fourier expansion method. As a result of varying cross sections, the solution includes a particular part additionally to the homogeneous one. A method for obtaining the particular solution is derived. Numerical tests on a VVER-1000 core are presented showing the impact of cross-section variations to some key parameters for reactor operation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1975

A cubic spline approximation for problems in fluid mechanics

A cubic spline approximation is presented which is suited for many fluid-mechanics problems. This procedure provides a high degree of accuracy, even with a nonuniform mesh, and leads to an accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several implicit and explicit integration schemes are presented. For two-dimensional flows, a spline-alternating-direction-implicit method is evaluated. The spline procedure is assessed, and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗
NASA NTRS1975

Viscous flow solutions with a cubic spline approximation

A cubic spline approximation is used for the solution of several problems in fluid mechanics. This procedure provides a high degree of accuracy even with a nonuniform mesh, and leads to a more accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several typical integration schemes are presented. For two-dimensional flows a spline-alternating-direction-implicit (SADI) method is evaluated. The spline procedure is assessed and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗
NASA NTRS1976

The absorption of trapped particles by the inner satellites of Jupiter and the radial diffusion coefficient of particle transport

The process of trapped particle absorption by the inner Jovian satellites is considered in detail taking into account both the particle and satellite motions in a magnetic dipole field which is displaced from the center of the planet and tilted with respect to the planetary rotation axis. An expression is derived for computing the sweeping time at a given satellite, defined as the time required for the satellite to sweep up a given fraction of the trapped particles within its sweeping region. By making use of the sweeping time and the radial diffusion equation of particle transport approximate expressions for the diffusion coefficient are derived. Measurements obtained by Pioneer 10 are then used to obtain estimates of the diffusion coefficient at the orbits of Io and Europa. We find that the diffusion coefficient is a function of energy and magnetic latitude for electrons in the energy range 0.7-14 MeV.

Mogro-Campero, A.↗
NASA NTRS1986

Monte Carlo study of a model of diffusion-controlled reactions

The Monte Carlo method is used to perform averages over sink configurations in the present study of diffusion-controlled reactions occurring between solute particles and immobile spherical sinks. In order to determine the average steady state solute concentration profile in a locally perturbed solution for sink volume fractions phi of less than 0.3, the diffusion equation in the monopolar plus dipolar approximation of diffusive couplings between the sinks is solved numerically. The Monte Carlo method is shown to be the most accurate and efficient in the phi = 0.001-0.1 region, where a system of only 25 sinks suffices and the monopolar approximation alone is sufficiently accurate.

Beenakker, C. W. J.↗
NASA NTRS1996

NASA flow fields analysis

The objectives of the present research are to improve design capabilities for low thrust rocket engines through understanding the detailed mixing and combustion processes in a representative combustor. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an on-going experimental program at NASA LERC. Detailed computational modeling involves the solution of both the two- and three-dimensional Navier Stokes equations, coupled with chemical reactions and the species diffusion equations. Computations of interest include both steady state and time-accurate flowfields and are obtained by means of LU approximate factorization in time and flux split upwinding differencing in space. The emphasis in the research is focused on using numerical analysis to understand detailed combustor flowfields, including the shear layer dynamics created between fuel film cooling and the core gas in the vicinity on the nearby combustor wall; the integrity and effectiveness of the coolant film; and three-dimensional fuel and oxidizer jet injection/mixing/combustion characteristics in the primary combustor along with their joint impacts on global engine performance.

Merkle, Charles L.↗
NASA NTRS1987

Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations

The current status of streamline-upwind/Petrov-Galerkin (SUPG) methods for the analysis of flow problems is surveyed in an analytical review. Problem areas addressed include classical Galerkin, upwind, artificial-diffusion, SUPG, discontinuous Galerkin, space-time FEM, and discontinuity-capturing approaches to the scalar advection-diffusion equation; incompressible flows; advective-diffusive systems; and the compressible Euler and Navier-Stokes equations. Graphs and diagrams are provided, and the good stability properties of state-of-the-art SUPG methods are pointed out.

Hughes, Thomas J. R.↗
NASA NTRS1985

Pitch Angle Distributions of Solar Energetic Particles and the Local Scattering Properties of the Interplanetary Magnetic Field

An approximate solution of the Fokker-Planck equation containing pitch angle scattering and adiabatic focusing is discussed. For modest focusing effects the omnidirectional density obeys an ordinary diffusion equation with a modified diffusion coefficient. The anisotropic part of the distribution function is properly normalized and split into an even and an odd part. The even part is determined by the ratio between the scattering mean free path and the focusing length and by the degree of polarization of the magnetic field fluctuations. The odd part is determined by the deviation of the pitch angle scattering from isotropic scattering. The method supplies a powerful tool to obtain the local characteristics of pitch angle scattering. It is insensitive to long lasting solar injections and to moderate radial variations of the mean free path. The method is applied to solar particle events observed on Helio-1 and -2.

Beeck, J.↗
DOE OSTI2023

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗
DOE OSTI2021

AIR Algebraic Multigrid for a Space-Time Hybridizable Discontinuous Galerkin Discretization of Advection(-Diffusion)

This article investigates the efficiency, robustness, and scalability of approximate ideal restriction (AIR) algebraic multigrid as a preconditioner in the all-at-once solution of a space-time hybridizable discontinuous Galerkin discretization of advection-dominated flows. The motivation for this study is that the time-dependent advection-diffusion equation can be seen as a “steady” advection-diffusion problem in $(d+1)$-dimensions and AIR has been shown to be a robust solver for steady advection-dominated problems. Numerical examples demonstrate the effectiveness of AIR as a preconditioner for advection-diffusion problems on fixed and time-dependent domains, using both slab-by-slab and all-at-once space-time discretizations, and in the context of uniform and space-time adaptive mesh refinement. A closer look at the geometric coarsening structure that arises in AIR also explains why AIR can provide robust, scalable, space-time convergence on advective and hyperbolic problems, while most multilevel parallel-in-time schemes struggle with such problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗
NASA NTRS1973

A rigorous cosmic-ray transport equation with no restrictions on particle energy.

A new transport equation for the cosmic-ray omnidirectional intensity is obtained. This equation follows exactly from the coupled pair of differential moment equations we presented earlier. It can be characterized as a nonlocal convection-diffusion equation in which the usual transport coefficients are replaced by time integral operators. The nonlocal equation is shown to reduce to the standard convection-diffusion form if the adiabatic approximation can be applied. In general, the adiabatic approximation does not apply; however, by going to the limits of infinite and zero gyroradius and, in addition, applying the adiabatic approximation, the large- and small-gyroradius transport theories due originally to Jokipii are regained. The validity of these theories as asymptotic limits and as approximate theories in the interplanetary magnetic field is discussed.

Klimas, A. J.↗
NASA NTRS1972

Applicability of a diffusion model to lateral transport in the terrestrial and lunar exospheres.

Kinetic theory is used to determine a series expansion of the vertical flux of particles in an exosphere in terms of time and space derivatives of particle concentration, exobase velocity, and temperature. For sufficiently large scale variations of these parameters in time and space, the series can be truncated to a form that is similar to a diffusion equation. Owing to this analogy, it is possible to unite the mathematical description of molecular diffusion, which governs thermospheric flow, and the corresponding exospheric equation by using effective transport coefficients which change smoothly with altitude through the transition from thermosphere to exosphere. A new definition of the exobase for lateral flow emerges from the analogy of exospheric and thermospheric diffusion, as the altitude where the horizontal mean free path length equals the mean horizontal extent of ballistic trajectories of the transported gas, as opposed to the scale height of the dominant gas which determines the exobase for escape. It is shown that the approximation of exospheric lateral flow as a diffusion process is applicable to global scale problems concerning terrestrial helium and heavier gases, and lunar gases heavier than helium.

Hodges, R. R., Jr.↗
DOE OSTI2022

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1976

Test particle propagation in magnetostatic turbulence. 1. Failure of the diffusion approximation

The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field.

Klimas, A. J.↗
NASA NTRS1985

Diffusion-convection function of cosmic rays

The fundamental properties and some numerical results of the solution of the diffusion equation of an impulsive cosmic-ray point source in an uniform, unbounded and spherically symmetrical moving medium is presented. The diffusion-convection(D-C) function is an elementary composite function of the solution of the D-C equation for the particles injected impulsively from a diffusive point source into the medium. It is the analytic solution derived by the dimensional method for the propagation equation of solar cosmic rays in the heliosphere, i.e. the interplanetary space. Because of the introduction of convection effect of solar wind, a nonhomogeneous term appears in the propagation equation, it is difficult to express its solution in terms of the ordinary special functions. The research made so far has led to a solution containing only the first order approximation of the convection effect.

Zhang, G.↗