Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
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 OSTI2015

Asymptotic diffusion limit of cell temperature discretisation schemes for thermal radiation transport

This paper attempts to unify the asymptotic diffusion limit analysis of thermal radiation transport schemes, for a linear-discontinuous representation of the material temperature reconstructed from cell centred temperature unknowns, in a process known as ‘source tilting’. The asymptotic limits of both Monte Carlo (continuous in space) and deterministic approaches (based on linear-discontinuous finite elements) for solving the transport equation are investigated in slab geometry. The resulting discrete diffusion equations are found to have nonphysical terms that are proportional to any cell-edge discontinuity in the temperature representation. Based on this analysis it is possible to design accurate schemes for representing the material temperature, for coupling thermal radiation transport codes to a cell centred representation of internal energy favoured by ALE (arbitrary Lagrange–Eulerian) hydrodynamics schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2014

Finite-difference schemes for anisotropic diffusion

In fusion plasmas diffusion tensors are extremely anisotropic due to the high temperature and large magnetic field strength. This causes diffusion, heat conduction, and viscous momentum loss, to effectively be aligned with the magnetic field lines. This alignment leads to different values for the respective diffusive coefficients in the magnetic field direction and in the perpendicular direction, to the extent that heat diffusion coefficients can be up to 10{sup 12} times larger in the parallel direction than in the perpendicular direction. This anisotropy puts stringent requirements on the numerical methods used to approximate the MHD-equations since any misalignment of the grid may cause the perpendicular diffusion to be polluted by the numerical error in approximating the parallel diffusion. Currently the common approach is to apply magnetic field-aligned coordinates, an approach that automatically takes care of the directionality of the diffusive coefficients. This approach runs into problems at x-points and at points where there is magnetic re-connection, since this causes local non-alignment. It is therefore useful to consider numerical schemes that are tolerant to the misalignment of the grid with the magnetic field lines, both to improve existing methods and to help open the possibility of applying regular non-aligned grids. To investigate this, in this paper several discretization schemes are developed and applied to the anisotropic heat diffusion equation on a non-aligned grid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2015

Entropy-based artificial viscosity stabilization for non-equilibrium Grey Radiation-Hydrodynamics

The entropy viscosity method is extended to the non-equilibrium Grey Radiation-Hydrodynamic equations. The method employs a viscous regularization to stabilize the numerical solution. The artificial viscosity coefficient is modulated by the entropy production and peaks at shock locations. The added dissipative terms are consistent with the entropy minimum principle. A new functional form of the entropy residual, suitable for the Radiation-Hydrodynamic equations, is derived. We demonstrate that the viscous regularization preserves the equilibrium diffusion limit. The equations are discretized with a standard Continuous Galerkin Finite Element Method and a fully implicit temporal integrator within the MOOSE multiphysics framework. The method of manufactured solutions is employed to demonstrate second-order accuracy in both the equilibrium diffusion and streaming limits. Several typical 1-D radiation-hydrodynamic test cases with shocks (from Mach 1.05 to Mach 50) are presented to establish the ability of the technique to capture and resolve shocks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

An atomistically informed mesoscale model for growth and coarsening during discharge in lithium-oxygen batteries

An atomistically informed mesoscale model is developed for the deposition of a discharge product in a Li-O{sub 2} battery. This mescocale model includes particle growth and coarsening as well as a simplified nucleation model. The model involves LiO{sub 2} formation through reaction of O{sub 2}{sup −} and Li{sup +} in the electrolyte, which deposits on the cathode surface when the LiO{sub 2} concentration reaches supersaturation in the electrolyte. A reaction-diffusion (rate-equation) model is used to describe the processes occurring in the electrolyte and a phase-field model is used to capture microstructural evolution. This model predicts that coarsening, in which large particles grow and small ones disappear, has a substantial effect on the size distribution of the LiO{sub 2} particles during the discharge process. The size evolution during discharge is the result of the interplay between this coarsening process and particle growth. The growth through continued deposition of LiO{sub 2} has the effect of causing large particles to grow ever faster while delaying the dissolution of small particles. The predicted size evolution is consistent with experimental results for a previously reported cathode material based on activated carbon during discharge and when it is at rest, although kinetic factors need to be included. The approach described in this paper synergistically combines models on different length scales with experimental observations and should have applications in studying other related discharge processes, such as Li{sub 2}O{sub 2} deposition, in Li-O{sub 2} batteries and nucleation and growth in Li-S batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗
DOE OSTI2013

Time-independent hybrid enrichment for finite element solution of transient conduction–radiation in diffusive grey media

We investigate the effectiveness of the partition-of-unity finite element method for transient conduction–radiation problems in diffusive grey media. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary diffusion approximation to the radiation in grey media. The coupled equations are integrated in time using a semi-implicit method in the finite element framework. We show that for the considered problems, a combination of hyperbolic and exponential enrichment functions based on an approximation of the boundary layer leads to improved accuracy compared to the conventional finite element method. It is illustrated that this approach can be more efficient than using h adaptivity to increase the accuracy of the finite element method near the boundary walls. The performance of the proposed partition-of-unity method is analyzed on several test examples for transient conduction–radiation problems in two space dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Fractional diffusion on bounded domains

We found that the mathematically correct specification of a fractional differential equation on a bounded domain requires specification of appropriate boundary conditions, or their fractional analogue. In this paper we discuss the application of nonlocal diffusion theory to specify well-posed fractional diffusion equations on bounded domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Raman Spectroscopy of Isotopic Water Diffusion in Ultraviscous, Glassy, and Gel States in Aerosol by Use of Optical Tweezers

The formation of ultraviscous, glassy, and amorphous gel states in aqueous aerosol following the loss of water results in nonequilibrium dynamics due to the extended time scales for diffusive mixing. Existing techniques for measuring water diffusion by isotopic exchange are limited by contact of samples with the substrate, and methods applied to infer diffusion coefficients from mass transport in levitated droplets requires analysis by complex coupled differential equations to derive diffusion coefficients. Here, we present a new technique that combines contactless levitation with aerosol optical tweezers with isotopic exchange (D 2 O/H 2 O) to measure the water diffusion coefficient over a broad range (D w ≈ 10 -12 -10 -17 m 2 s -1 ) in viscous organic liquids (citric acid, sucrose, and shikimic acid) and inorganic gels (magnesium sulfate, MgSO 4 ). For the organic liquids in binary and ternary mixtures, D w depends on relative humidity and follows a simple compositional Vignes relationship. In MgSO 4 droplets, water diffusivity decreases sharply with water activity and is consistent with predictions from percolation theory. These measurements show that, by combining micrometer-sized particle levitation (a contactless measurement with rapid mixing times) with an established probe of water diffusion, D w can be simply and directly quantified for amorphous and glassy states that are inaccessible to existing methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Migration of the separation point on a deforming cylinder

An iterative scheme of solving the Navier-Stokes equations for unsteady two-dimensional flow about a deforming and translating cylinder is given. In the kth iteration, the convected vorticity appears as the source term in an equation of transient diffusion of vorticity. A novel integral transform is used to reduce this transient vorticity equation to a k-dimensional heat equation. The bounded solution of this equation is obtained with a general method of superposition for problems involving a moving boundary. An equation describing the migration of the separation point on a deforming cylinder in unsteady cross flows is derived from the analytically obtained velocity field. The radial expansion of the cylinder surface is shown to hasten the separation time and to increase the separation angle. The results imply that in the steady flow past a body that is moving forward and sinking at constant rates, the locus of points at which the azimuthal component of skin friction changes sign originates on the leeward ray at a point downstream of the front tip.

Lin, S. P.↗