Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion equations”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Correlation functions for the fractional generalized Langevin equation in the presence of internal and external noise

We study generalized fractional Langevin equations in the presence of a harmonic potential. General expressions for the mean velocity and particle displacement, the mean squared displacement, position and velocity correlation functions, as well as normalized displacement correlation function are derived. We report exact results for the cases of internal and external friction, that is, when the driving noise is either internal and thus the fluctuation-dissipation relation is fulfilled or when the noise is external. The asymptotic behavior of the generalized stochastic oscillator is investigated, and the case of high viscous damping (overdamped limit) is considered. Additional behaviors of the normalized displacement correlation functions different from those for the regular damped harmonic oscillator are observed. In addition, the cases of a constant external force and the force free case are obtained. The validity of the generalized Einstein relation for this process is discussed. The considered fractional generalized Langevin equation may be used to model anomalous diffusive processes including single file-type diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Stochastic models for surface diffusion of molecules

We derive a stochastic model for the surface diffusion of molecules, starting from the classical equations of motion for an N-atom molecule on a surface. The equation of motion becomes a generalized Langevin equation for the center of mass of the molecule, with a non-Markovian friction kernel. In the Markov approximation, a standard Langevin equation is recovered, and the effect of the molecular vibrations on the diffusion is seen to lead to an increase in the friction for center of mass motion. This effective friction has a simple form that depends on the curvature of the lowest energy diffusion path in the 3N-dimensional coordinate space. We also find that so long as the intramolecular forces are sufficiently strong, memory effects are usually not significant and the Markov approximation can be employed, resulting in a simple one-dimensional model that can account for the effect of the dynamics of the molecular vibrations on the diffusive motion.

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

The B{sup −3/2} diffusion in magnetized plasma

Fluid equations with transport effects have been derived from summing up discrete particle effects in the configuration space. The diffusion coefficient takes a form different from, but equivalent to, the Kubo formula. By taking into account the decoherence along the field lines, the diffusive damping in the perpendicular direction, and the finite Larmor radius effect in the cyclotron resonance, it is shown that the diffusion of magnetized 3d plasma has a regime of B{sup −3/2} scaling which lies between the Bohm diffusion of B{sup −1} and the classical diffusion of B{sup −2}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Synthesis of Jicamarca data during the great storm of March 8, 1970.

Continuous measurement of electron temperature and electron concentration during the great magnetic storm of Mar. 8, 1970, when the F2 peak rose as high as 800 km. Continuous measurements of the vertical component of E x B drift were made from 1130 to 1730 LT. A synthesis of the observed electron concentrations was made using a model that solved the time-dependent electron-ion continuity equation for O(+) and the molecular ions; the latter chemistry was included to improve the solutions at lower heights. It was found, using the observed values of vertical drift, that most of the observed features of the F2-region ionosphere were synthesized by the model. An interpretation of the synthesis that emphasizes the importance of horizontal diffusion at the equator is given. Observations are presented for Mar. 7, 8, and 9, 1970, including the Huancayo magnetograms, which display a close correlation with the drift measurements.

Woodman, R. F.↗
DOE OSTI2020

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the polynomial order used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Verification of the DIF3D Software to Support Fast Reactor Analysis (Rev. 3)

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2020

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗
NASA NTRS1991

Single Doppler radar and lidar studies of the dynamics of mixed layers and overlying stable regions

Some examples of profiles of PBL quantities retrieved from Doppler radar and lidar observations are shown, and the extent to which the deduced quantities agree with the traditional eddy viscosity concept is examined. It is found in one case that the concept of 'top down bottom up' diffusion introduced by Wyngaard (1983) is more pertinent than that of the eddy viscosity. In another case it is found that, in stable layers with sufficiently large Richardson number, the momentum fluxes could be counter-gradient. In this latter case the use of a negative turbulence diffusion coefficient would result in an ill-posed differential equation, thus confusing the 'diffusion problem'.

Xu, Mei↗
NASA NTRS1990

Injection and diffusive transport of suprathermal through energetic solar flare protons (35 keV to 20 MeV)

Consideration is given to the injection and interplanetary propagation of low-energy protons caused by the solar particle event of July 20, 1981, in which flare protons in the range from 35 keV to 20 MeV were observed by instruments on ISEE 3. The observed time-intensity and time-anisotropy profiles were fitted over the entire energy range using a model based on the spherically symmetric Fokker-Plank equation, including convection, diffusion, and adiabatic deceleration. The results are used to discuss the behavior of the radial interplanetary diffusion coefficient and the scattering mean free path for protons. Also, evidence is found for diffusive coronal shock acceleration of protons during the event.

Beeck, J.↗
NASA NTRS1975

Wave theory of turbulence in compressible media (acoustic theory of turbulence)

The generation and the transmission of sound in turbulent flows are treated as one of the several aspects of wave propagation in turbulence. Fluid fluctuations are decomposed into orthogonal Fourier components, with five interacting modes of wave propagation: two vorticity modes, one entropy mode, and two acoustic modes. Wave interactions, governed by the inhomogeneous and nonlinear terms of the perturbed Navier-Stokes equations, are modeled by random functions which give the rates of change of wave amplitudes equal to the averaged interaction terms. The statistical framework adopted is a quantum-like formulation in terms of complex distribution functions. The spatial probability distributions are given by the squares of the absolute values of the complex characteristic functions. This formulation results in nonlinear diffusion-type transport equations for the probability densities of the five modes of wave propagation.

Kentzer, C. P.↗
DOE OSTI2014

Projection of two-dimensional diffusion in narrow asymmetric channels onto the longitudinal direction

Diffusive transport of particles is a ubiquitous feature of physical, chemical and biological systems. Typical structures like pores, tubes or fibers, are quasi one-dimensional, such that we need to solve 2+1 or 3+1 dimensional differential equations to describe correctly transport along them. The so-called Fick-Jacobs approach dramatically simplifies the problem if one assumes that a solute distribution in any cross-section of the channel is uniform at equilibrium. That study focuses on the mapping of the diffusion equation in a two-dimensional narrow asymmetric channel of varying cross section onto the longitudinal coordinate. We present a generalization to the case of an asymmetric channel using the projection method introduced earlier by Kalinay and Percus. We derive an expansion of the effective diffusion coefficient, which represents corrections to the Fick-Jacobs equation and contains the well-known previous results as special cases. Finally, we study numerically some specific two-dimensional asymmetric channel configurations to test and show the broader applicability of this effective diffusion coefficient formula.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

Analytical Modeling of Operating Characteristics of Premixing-prevaporizing Fuel-air Mixing Passages

A model for predicting the distribution of liquid fuel droplets and fuel vapor in premixing-prevaporizing fuel-air mixing passages of the direct injection type is described. This model consists of three computer programs: a calculation of the two dimensional or axisymmetric air flow field neglecting the effects of fuel; a calculation of the three dimensional fuel droplet trajectories and evaporation rates in a known, moving air flow; and a calculation of fuel vapor diffusing into a moving three dimensional air flow with source terms dependent on the droplet evaporation rates. The air flow calculation can treat compressible swirling flows in arbitrary ducts with arbitrary distributions of temperature and velocity as initial conditions. The fuel droplets are treated as initial conditions. The fuel droplets are treated as individual particle classes each satisfying Newton's law, a heat transfer, and a mass transfer equation. The vapor diffusion calculation treats three dimensional, gas phase, turbulent diffusion processes with the turbulence level determined by the air flow calculations and the source terms determined by the droplet evaporation rates.

Anderson, O. L.↗
NASA NTRS1984

Pitch angle diffusion in the Jovian magnetodisc

Using a linear approximation to the field variation and a mathematical technique for distorting the orbit integration into the complex gyrophase plane, a formula is derived for the nonadiabatic change in magnetic moment. This change varies as the cosine of the particle gyrophase at the equator, so that if a particle traverses the equator several times with different values of this gyrophase, it suffers randomly correlated values for the nonadiabatic change in magnetic moment. The algebraic formula for this nonadiabatic change is found to be in excellent agreement with previous numerical computations. The phase space density, averaged over several bounce periods, is found to satisfy a diffusion equation in pitch angle. The diffusion coefficient is dependent on particle energy, pitch angle, and the field line along which the particles are moving.

Birmingham, T. J.↗
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2023

Merging of the superbanana plateau and $\sqrt{v}$ transport regimes in nearly quasisymmetric stellarators

We report alpha particle confinement is one of the most demanding issues for stellarators. It now seems clear that it is possible to design optimized stellarators that confine the background plasma at near tokamak radial transport levels. Moreover, adequate collisionless alpha particle confinement is possible in the core of a highly optimized stellarator. Here, the collisional confinement of barely trapped alphas in an optimized stellarator is considered by accounting for the resonance due to the reversal in direction of the drift within a flux surface and investigating the sensitive role of magnetic shear in keeping this resonance close to the passing boundary in some nearly quasisymmetric stellarator configurations. The treatment relies on a narrow collisional boundary layer formulation that combines the responses of both these resonant pitch angle alphas and the remaining barely trapped alphas. A novel merged regime treatment leads to explicit expressions for the energy diffusivity for both superbanana plateau (or resonant plateau) and $\sqrt{v}$ transport in the large aspect ratio limit for a slowing down tail alpha distribution function, where v is the effective pitch angle scattering collision frequency of the trapped alphas off the background ions. Depending on the details of the optimization scheme and the sign of the magnetic shear, modest magnetic shear can be used to reduce superbanana (or resonant) plateau transport to below the $\sqrt{v}$ transport level. In addition, a quasilinear equation retaining spatial diffusion is derived for a general alpha distribution function that allows the radial alpha transport to modify the distribution so it is no longer isotropic in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Statistical CT noise reduction with multiscale decomposition and penalized weighted least squares in the projection domain

Purposes: The suppression of noise in x-ray computed tomography (CT) imaging is of clinical relevance for diagnostic image quality and the potential for radiation dose saving. Toward this purpose, statistical noise reduction methods in either the image or projection domain have been proposed, which employ a multiscale decomposition to enhance the performance of noise suppression while maintaining image sharpness. Recognizing the advantages of noise suppression in the projection domain, the authors propose a projection domain multiscale penalized weighted least squares (PWLS) method, in which the angular sampling rate is explicitly taken into consideration to account for the possible variation of interview sampling rate in advanced clinical or preclinical applications. Methods: The projection domain multiscale PWLS method is derived by converting an isotropic diffusion partial differential equation in the image domain into the projection domain, wherein a multiscale decomposition is carried out. With adoption of the Markov random field or soft thresholding objective function, the projection domain multiscale PWLS method deals with noise at each scale. To compensate for the degradation in image sharpness caused by the projection domain multiscale PWLS method, an edge enhancement is carried out following the noise reduction. The performance of the proposed method is experimentally evaluated and verified using the projection data simulated by computer and acquired by a CT scanner. Results: The preliminary results show that the proposed projection domain multiscale PWLS method outperforms the projection domain single-scale PWLS method and the image domain multiscale anisotropic diffusion method in noise reduction. In addition, the proposed method can preserve image sharpness very well while the occurrence of 'salt-and-pepper' noise and mosaic artifacts can be avoided. Conclusions: Since the interview sampling rate is taken into account in the projection domain multiscale decomposition, the proposed method is anticipated to be useful in advanced clinical and preclinical applications where the interview sampling rate varies.

60 APPLIED LIFE SCIENCES↗