Search NASA⌕ Search

SEARCH · Search NASA

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

An efficient transport solver for tokamak plasmas

A simple approach to efficiently solve a coupled set of 1-D diffusion-type transport equations with a stiff transport model for tokamak plasmas is presented based on the 4th order accurate Interpolated Differential Operator scheme along with a nonlinear iteration method derived from a root-finding algorithm. Here, numerical tests using the Trapped Gyro-Landau-Fluid model show that the presented high order method provides an accurate transport solution using a small number of grid points with robust nonlinear convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Stability analysis of nonlinear two-grid method for multigroup neutron diffusion problems

We present theoretical analysis of a nonlinear acceleration method for solving multigroup neutron diffusion problems. This method is formulated with two energy grids that are defined by (i) fine-energy groups structure and (ii) coarse grid with just a single energy group. The coarse-grid equations are derived by averaging of the multigroup diffusion equations over energy. The method uses a nonlinear prolongation operator. We perform stability analysis of iteration algorithms for inhomogeneous (fixed-source) and eigenvalue neutron diffusion problems. To apply Fourier analysis the equations of the method are linearized about solutions of infinite-medium problems. The developed analysis enables us to predict convergence properties of this two-grid method in different types of problems. Numerical results of problems in 2D Cartesian geometry are presented to confirm theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Nonclassical Particle Transport in the 1-D Diffusive Limit

A nonclassical linear Boltzmann equation has been recently proposed to address nonexponential attenuation of the particle flux in certain inhomogeneous random media applications. In particular, this effect arises in Pebble Bed reactor cores, in which the locations of the pebbles are spatially correlated. In this paper, we investigate nonclassical particle transport taking place in a 1-D random periodic diffusive system. We provide computational results that validate the theoretical predictions, demonstrating for the first time that the solution of the nonclassical particle transport equation is well-approximated by the solution of the nonclassical diffusion equation. For simplicity, we consider the following assumptions: (i) transport occurs in rod geometry, in which particles can only move in the directions μ = ±1; (ii) transport is monoenergetic; and (iii) scattering is isotropic. This paper provides numerical results that demonstrate the validity of the nonclassical diffusion approximation to the nonclassical transport equation in certain 1-D diffusive systems. To our knowledge, this is the first time computational results that validate the predictions of the asymptotic theory for the nonclassical model are presented. This result provides a more solid foundation in which to improve this theory for relevant nuclear applications. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

A Study of the von Neumann Stability Analysis of a Semi-implicit Numerical Method Applied Over the Radial Impurity Transport Equation in Tokamak Plasma

The radial impurity transport equation for tokamak plasma is a form of diffusion–convection–reaction equation. The impurity transport equation is solved to determine the distribution of impurity (non-fuel) ion species with different ionization states perpendicular to magnetic surfaces of tokamak plasma. The equation for each charge (ionization) state Z is a non-linear, second-order in space, first-order in time, parabolic partial differential equation coupled to the previous Z − 1 and the next Z + 1 charge states of the impurity species through its reaction term. The number of differential equations to be solved simultaneously is hence determined by the number of ionization states of the impurity species studied. The solution to the set of these coupled equations can be obtained using a semi-implicit numerical method applied on it. The present study describes the application of von Neumann stability analysis over the semi-implicit numerical method applied over the radial impurity transport equation and determines a generic stability criterion for the method. The stability analysis is further illustrated using the geometry of Aditya tokamak installed at the Institute for Plasma Research Gandhinagar, India as an example. The impurity species considered is oxygen (Atomic number = 8). This leads to a set of eight coupled equations for charge states Z = 1 to 8 over which von Neumann analysis is illustrated in present study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1998

Multi-Component Diffusion with Application To Computational Aerothermodynamics

The accuracy and complexity of solving multicomponent gaseous diffusion using the detailed multicomponent equations, the Stefan-Maxwell equations, and two commonly used approximate equations have been examined in a two part study. Part I examined the equations in a basic study with specified inputs in which the results are applicable for many applications. Part II addressed the application of the equations in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) computational code for high-speed entries in Earth's atmosphere. The results showed that the presented iterative scheme for solving the Stefan-Maxwell equations is an accurate and effective method as compared with solutions of the detailed equations. In general, good accuracy with the approximate equations cannot be guaranteed for a species or all species in a multi-component mixture. 'Corrected' forms of the approximate equations that ensured the diffusion mass fluxes sum to zero, as required, were more accurate than the uncorrected forms. Good accuracy, as compared with the Stefan- Maxwell results, were obtained with the 'corrected' approximate equations in defining the heating rates for the three Earth entries considered in Part II.

Sutton, Kenneth↗
NASA NTRS1977

The relative diffusion of strong magnetic fields and tenuous gases

General equations are given for treating the diffusion of a magnetic field in a fluid with electrical conductivity and a resistive diffusion coefficient. Consequences of nonnegligible gas motion are investigated by considering the one-dimensional problem where a tenuous isothermal gas with pressure expressed as a function of x-coordinate and time is enmeshed in a magnetic field in the y-direction. Approximate equations are obtained for the case of weak fluctuations in gas or field pressure, and the properties of a nonlinear diffusion equation are illustrated with the example of an initial thin slab of gas confined by a uniform field on either side. The advance of a gas front into a nonvanishing gas density is illustrated by using the progressive wave solution. Application of the results to the evolution of streamers and filaments observed in the solar corona is briefly discussed

Parker, E. N.↗
DOE OSTI2016

Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene, and carbon tetrachloride

Mutual diffusion coefficients of all 20 binary liquid mixtures that can be formed out of methanol, ethanol, acetone, benzene, cyclohexane, toluene, and carbon tetrachloride without a miscibility gap are studied at ambient conditions of temperature and pressure in the entire composition range. The considered mixtures show a varying mixing behavior from almost ideal to strongly non-ideal. Predictive molecular dynamics simulations employing the Green-Kubo formalism are carried out. Radial distribution functions are analyzed to gain an understanding of the liquid structure influencing the diffusion processes. It is shown that cluster formation in mixtures containing one alcoholic component has a significant impact on the diffusion process. The estimation of the thermodynamic factor from experimental vapor-liquid equilibrium data is investigated, considering three excess Gibbs energy models, i.e., Wilson, NRTL, and UNIQUAC. It is found that the Wilson model yields the thermodynamic factor that best suits the simulation results for the prediction of the Fick diffusion coefficient. Four semi-empirical methods for the prediction of the self-diffusion coefficients and nine predictive equations for the Fick diffusion coefficient are assessed and it is found that methods based on local composition models are more reliable. Finally, the shear viscosity and thermal conductivity are predicted and in most cases favorably compared with experimental literature values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Using Directional Diffusion Coefficients for Nonlinear Diffusion Acceleration of the First Order S{sub N} Equations in Near-Void Regions

The common definition of the diffusion coefficient as the inverse of three times the transport cross section is not compatible with voids. Morel introduced a non-local tensor diffusion coefficient that remains finite in voids. It can be obtained by solving an auxiliary transport problem without scattering or fission. Larsen and Trahan successfully applied this diffusion coefficient for enhancing the accuracy of diffusion solutions of very high temperature reactor (VHTR) problems that feature large, optically thin channels in the z-direction. It was demonstrated that a significant reduction of error can be achieved in particular in the optically thin region. Along the same line of thought, non-local diffusion tensors have been applied to modeling the TREAT reactor confirming the findings of Larsen and Trahan. Previous work of the authors have introduced a flexible Nonlinear-Diffusion Acceleration (NDA) method for the first order S{sub N} equations discretized with the discontinuous finite element method (DFEM). This NDA method uses a scalar diffusion coefficient in the low-order system that is obtained as the flux weighted average of the inverse transport cross section. Hence, it suffers from very large and potentially unbounded diffusion coefficients in the low order problem. However, it was noted that the choice of the diffusion coefficient does not influence consistency of the method at convergence and hence the diffusion coefficient is essentially a free parameter. The choice of the diffusion coefficient does, however, affect the convergence behavior of the nonlinear diffusion iterations. Within this work we use Morel's non-local diffusion coefficient in the aforementioned NDA formulation in lieu of the flux weighted inverse of three times the transport cross section. The goal of this paper is to demonstrate that significant enhancement of the spectral properties of NDA can be achieved in near void regions. For testing the spectral properties of the NDA with non-local diffusion coefficients, the periodical horizontal interface problem is used. This problem consists of alternating stripes of optically thin and thick materials both of which feature scattering ratios close to unity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Oxygen Reduction on Stainless Steel in Concentrated Chloride Media

In this work, a rotating disk electrode was used to measure the cathodic kinetics on stainless steel as a function of diffusion layer thickness (6 to 60 pm) and chloride concentration (0.6 to 5.3 M NaCl). It was found that while the cathodic kinetics followed the Levich equation for large diffusion layer thicknesses, the Levich equation over predicts the mass-transfer limited current density for diffusion layer thicknesses less than 20 Am. Also, an unusual transitory response between the activation and mass-transfer controlled regions was observed for small diffusion layer thicknesses that was more apparent in lower concentration solutions. The presence and reduction of an oxide film and a transition in the oxygen reduction mechanism were identified as possible reasons for this response. The implications of these results on atmospheric corrosion kinetics under thin electrolyte layers is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A score-based diffusion model approach for adaptive learning of stochastic partial differential equation solutions

In this paper, we propose a novel framework for adaptively learning the time-evolving solutions of stochastic partial differential equations (SPDEs) using score-based diffusion models within a recursive Bayesian inference setting. SPDEs play a central role in modeling complex physical systems under uncertainty, but their numerical solutions often suffer from model errors and reduced accuracy due to incomplete physical knowledge and environmental variability. To address these challenges, we encode the governing physics into the score function of a diffusion model using simulation data and incorporate observational information via a likelihood-based correction in a reverse-time stochastic differential equation. This enables adaptive learning through iterative refinement of the solution as new data becomes available. To improve computational efficiency in high-dimensional settings, we introduce the ensemble score filter, a training-free approximation of the score function designed for real-time inference. Numerical experiments on benchmark SPDEs demonstrate the accuracy and robustness of the proposed method under sparse and noisy observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Topological Hall effect in diffusive ferromagnetic thin films with spin-flip scattering

In this paper, we study the topological Hall (TH) effect in a diffusive ferromagnetic metal thin film by solving a Boltzmann transport equation in the presence of spin-flip scattering. A generalized spin-diffusion equation is derived which contains an additional source term associated with the gradient of the emergent magnetic field that arises from skyrmions. Because of the source term, spin accumulation may build up in the vicinity of the skyrmions. This gives rise to a spin-polarized diffusion current that in general suppresses the bulk TH current. Only when the spin-diffusion length is much smaller than the skyrmion size does the TH resistivity approach the value derived by Bruno et al. [Phys. Rev. Lett. 93, 096806 (2004)]. Finally, we derive a general expression of the TH resistivity that applies to thin-film geometries with spin-flip scattering, and show that the corrections to the TH resistivity become large when the size of room temperature skyrmions is further reduced to tens of nanometers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2011

A Simple, Analytical Model of Collisionless Magnetic Reconnection in a Pair Plasma

A set of conservation equations is utilized to derive balance equations in the reconnection diffusion region of a symmetric pair plasma. The reconnection electric field is assumed to have the function to maintain the current density in the diffusion region, and to impart thermal energy to the plasma by means of quasi-viscous dissipation. Using these assumptions it is possible to derive a simple set of equations for diffusion region parameters in dependence on inflow conditions and on plasma compressibility. These equations are solved by means of a simple, iterative, procedure. The solutions show expected features such as dominance of enthalpy flux in the reconnection outflow, as well as combination of adiabatic and quasi-viscous heating. Furthermore, the model predicts a maximum reconnection electric field of E(sup *)=0.4, normalized to the parameters at the inflow edge of the diffusion region.

Hesse, Michael↗
DOE OSTI2025

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗
DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗
NASA NTRS1975

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗
NASA NTRS1973

Vertical distribution of vibrational energy of molecular nitrogen in a stable auroral red arc and its effect on ionospheric electron densities

Previous solutions of the problem of the distribution of vibrationally excited molecular nitrogen in the thermosphere have either assumed a Boltzmann distribution and considered diffusion as one of the loss processes or solved for the energy level populations and neglected diffusion. Both of the previous approaches are combined by solving the time dependent continuity equations, including the diffusion process, for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to a stable auroral red arc. The primary source of molecular nitrogen excitation was subexcitation, and inelastic collisions between thermal electrons and molecular nitrogen. The reaction rates for this process were calculated from published cross section calculations. The loss processes for vibrational energy were electron and atomic oxygen quenching and vibrational energy exchange. The coupled sets of nonlinear, partial differential equations were solved numerically by employing finite difference equations.

Newton, G. P.↗