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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Lyα heating of inhomogeneous high-redshift intergalactic medium

The intergalactic medium (IGM) prior to the epoch of reionization consists mostly of neutral hydrogen gas. Lyman-α (Lyα) photons produced by early stars resonantly scatter off hydrogen atoms, causing energy exchange between the radiation field and the gas. This interaction results in moderate heating of the gas due to the recoil of the atoms upon scattering, which is of great interest for future studies of the pre-reionization IGM in the H I 21 cm line. We investigate the effect of this Lyα heating in the IGM with linear density, temperature, and velocity perturbations. Perturbations smaller than the diffusion length of photons could be damped due to heat conduction by Lyα photons. The scale at which damping occurs and the strength of this effect depend on various properties of the gas, the flux of Lyα photons, and the way in which photon frequencies are redistributed upon scattering. To find the relevant length scale and the extent to which Lyα heating affects perturbations, we calculate the gas heating rates by numerically solving linearized Boltzmann equations in which scattering is treated by the Fokker-Planck approximation. We find that (1) perturbations add a small correction to the gas heating rate, and (2) the damping of temperature perturbations occurs at scales with comoving wavenumber k ≳ 10{sup 4} Mpc{sup –1}, which are much smaller than the Jeans scale and thus unlikely to substantially affect the observed 21 cm signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

On the evolution of cured voxel in bulk photopolymerization upon focused Gaussian laser exposure

Unconstrained depth photopolymerization is emerging as a promising technique for fabrication of several polymer microstructures such as self propagating waveguides, 3D freeform structures by bulk lithography, and polymer nanoparticles by flash exposure. Experimental observations reveal governing physics beyond Beer Lambert's law and scattering effects. This paper seeks to model unconstrained depth photopolymerization using classical nonlinear Schrödinger equation coupled with transient diffusion phenomenon. The beam propagation part of the proposed model considers scattering effects induced due to spatial variation of the refractive index as a function of the beam intensity. The critical curing energy model is used to further predict profile of polymerized voxel. Profiles of photopolymerized voxel simulated using proposed model are compared with the corresponding experimental results for several cases of exposure dose and duration. The comparison shows close match leading to conclusion that the experimentally observed deviation from Beer Lambert's law is indeed due to combined effect of diffusion of photoinitiator and scattering of light because of change in the refractive index.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

The equilibrium-diffusion limit for radiation hydrodynamics

The equilibrium-diffusion approximation (EDA) is used to describe certain radiation-hydrodynamic (RH) environments. When this is done the RH equations reduce to a simplified set of equations. The EDA can be derived by asymptotically analyzing the full set of RH equations in the equilibrium-diffusion limit. Here, we derive the EDA this way and show that it and the associated set of simplified equations are both first-order accurate with transport corrections occurring at second order. Having established the EDA’s first-order accuracy we then analyze the grey nonequilibrium-diffusion approximation and the grey Eddington approximation and show that they both preserve this first-order accuracy. Further, these approximations preserve the EDA’s first-order accuracy when made in either the comoving-frame (CMF) or the lab-frame (LF). And while analyzing the Eddington approximation, we found that the CMF and LF radiation-source equations are equivalent when neglecting O(β 2 ) terms and compared in the LF. Of course, the radiation pressures are not equivalent. It is expected that simplified physical models and numerical discretizations of the RH equations that do not preserve this first-order accuracy will not retain the correct equilibrium-diffusion solutions. As a practical example, we show that nonequilibrium-diffusion radiative-shock solutions devolve to equilibrium-diffusion solutions when the asymptotic parameter is small.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Final Report of the Project "From the finite element method to the virtual element method"

The Finite Element Method (FEM) is a powerful numerical tool that is being used in a large number of engineering applications. The FEM is constructed on triangular/tetrahedral and quadrilateral/hexahedral meshes. Extending the FEM to general polygonal/polyhedral meshes in straightforward way turns out to be extremely difficult and leads to very complex and computationally expensive schemes. The reason for this failure is that the construction of the basis functions on elements with a very general shape is a non-trivial and complex task. In this project we developed a new family of numerical methods, dubbed the Virtual Element Method (VEM) for the numerical approximation of partial differential equations (PDE) of elliptic type suitable to polygonal and polyhedral unstructured meshes. We successfully formulated, implemented and tested these methods and studied both theoretically and numerically their stability, robustness and accuracy for diffusion problems, convection-reaction-diffusion problems, the Stokes equations and the biharmonic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MFiX: Fractional-Step Method Implementation

A comprehensive, multiphase computational fluid dynamics (CFD) simulation solves several coupled transport equations including continuity, momentum, species, and energy. Chemical reactions further couple these equations through heats of reaction and rates of formation of products and rates of destruction of reactants. A fractional-step method separates changes attributed to chemical reactions from transport phenomena like convection and diffusion. When the governing equations are split into the transport and reacting components, efficient and independent methodologies can be exploited to solve the different systems. Specifically, discretization of field variable transport equations results in large, sparse matrices which are loosely coupled. These systems are solved in succession using iterative techniques that take advantage of the matrix structure. In contrast, chemical reactions tightly couple field variables locally within the domain (e.g., within a single computational cell) resulting in low dimensional but dense, nonlinear systems that are better solved using direct integration techniques.

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

Mathematical model of diffusion-limited evolution of multiple gas bubbles in tissue

Models of gas bubble dynamics employed in probabilistic analyses of decompression sickness incidence in man must be theoretically consistent and simple, if they are to yield useful results without requiring excessive computations. They are generally formulated in terms of ordinary differential equations that describe diffusion-limited gas exchange between a gas bubble and the extravascular tissue surrounding it. In our previous model (Ann. Biomed. Eng. 30: 232-246, 2002), we showed that with appropriate representation of sink pressures to account for gas loss or gain due to heterogeneous blood perfusion in the unstirred diffusion region around the bubble, diffusion-limited bubble growth in a tissue of finite volume can be simulated without postulating a boundary layer across which gas flux is discontinuous. However, interactions between two or more bubbles caused by competition for available gas cannot be considered in this model, because the diffusion region has a fixed volume with zero gas flux at its outer boundary. The present work extends the previous model to accommodate interactions among multiple bubbles by allowing the diffusion region volume of each bubble to vary during bubble evolution. For given decompression and tissue volume, bubble growth is sustained only if the bubble number density is below a certain maximum.

Evaluation Studies↗
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2022

Current challenges in the physics of white dwarf stars

White dwarfs are a class of stars with unique physical properties. They present many challenging problems whose solution requires advanced theories of dense matter, state-of-the-art experimental techniques, and extensive computing efforts. New ground- and space-based observatories will soon provide an increasingly detailed view of white dwarf stars and reveal new phenomena that will challenge our models. This review is an introduction to the nature of white dwarfs, the physical processes that determine their structure and evolution, and the physical conditions they span. Here we discuss a wide variety of currently unsolved or partially resolved problems in their constitutive physics that are broadly related to equations of state, transport processes and opacities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Gravitational encounters and the evolution of galactic nuclei. II. Classical and resonant relaxation

Direct numerical integrations of the Fokker–Planck equation in energy-angular momentum space are carried out for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm, which was described in detail in an earlier paper, includes diffusion coefficients that describe the effects of both random (“classical”) and correlated (“resonant”) encounters. Steady-state solutions are similar to the Bahcall–Wolf solution, n(r)∝r{sup −7/4}, but are modified at small radii due to the higher rate of diffusion in angular momentum, which results in a low-density core. The core radius is a few percent of the influence radius of the SBH. The corresponding phase-space density f(E,L) drops nearly to zero at low energies, implying almost no stars on tightly bound orbits about the SBH. Steady-state rates of stellar disruption are presented, and a simple analytic expression is found that reproduces the numerical feeding rates with good accuracy. The distribution of periapsides of disrupted stars is computed. Time-dependent solutions, f(E,L,t), are also computed, starting from initial conditions similar to those produced by a binary SBH. In these models, feeding rates evolve on two timescales: rapid evolution during which the region evacuated by the massive binary is refilled by angular-momentum diffusion; and slower evolution as diffusion in energy causes the density profile at large radii to attain the Bahcall–Wolf form.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Basic Studies of Non-Diffusive Transport in Plasmas

The project expanded and developed mathematical descriptions, and corresponding numerical modeling, of non-diffusive transport to incorporate new perspectives derived from basic transport experiments performed in the LAPD device at UCLA, and at fusion devices throughout the world. By non-diffusive it is meant that the transport of fundamental macroscopic parameters of a system, such as temperature and density, does not follow the standard diffusive behavior predicted by a classical Fokker-Planck equation. The appearance of non-diffusive behavior is often related to underlying microscopic processes that cause the value of a system parameter, at one spatial position, to be linked to distant events, i.e., non-locality. In the LAPD experiments the underlying process was traced to large amplitude, coherent drift-waves that give rise to chaotic trajectories. Significant advances were made in this project. The results have lead to a new perspective about the fundamentals of edge transport in magnetically confined plasmas; the insight has important consequences for worldwide studies in fusion devices. Progress was also made in advancing the mathematical techniques used to describe fractional diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Diffusion and viscosity coefficients for helium

The first order Boltzmann-Fokker-Planck equation is solved numerically to obtain diffusion and viscosity coefficients for a ternary gas mixture composed of electron, protons, and helium. The coefficients are tabulated for five He/H abundances ranging from 0.01 to 10 and for both He II and He III. Comparison with Burgers's thermal diffusion coefficients reveals a maximum difference of 9-10% for both He II and He III throughout the range of helium abundances considered. The viscosity coefficients are compared to those of Chapman and Cowling and show a maximum difference of only 5-6% for He II but 15-16% for He III. For the astrophysically important gas mixtures, it is concluded that the results of existing studies which employed Burgers's or Chapman and Cowling's coefficients will remain substantially unaltered.

Roussel-Dupre, R.↗
DOE OSTI2019

A Finite-Element Approach for the Partitioning of Carbon in Q&P Steel

This paper reports a Galerkin finite-element analysis of carbon partitioning from martensite into austenite during the quenching and partitioning (Q&P) processes in steels. In contrast to classical or sophisticated diffusion field models, an alternative nonlinear governing equation based on chemical potential and composition is considered. The model is applied to simulate the carbon partitioning of modified 22MnB5 alloys assuming an immobile austenite–martensite phase boundary and, in turn, is compared with experimental measurements of the volume fraction and carbon content of retained austenite. The simulations show outstanding results on the influence of Si and C content on the partitioning parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Enhancement of the droplet nucleation in a dense supersaturated Lennard-Jones vapor

The vapor–liquid nucleation in a dense Lennard-Jones system is studied analytically and numerically. A solution of the nucleation kinetic equations, which includes the elementary processes of condensation/evaporation involving the lightest clusters, is obtained, and the nucleation rate is calculated. Based on the equation of state for the cluster vapor, the pre-exponential factor is obtained. The latter diverges as a spinodal is reached, which results in the nucleation enhancement. The work of critical cluster formation is calculated using the previously developed two-parameter model (TPM) of small clusters. A simple expression for the nucleation rate is deduced and it is shown that the work of cluster formation is reduced for a dense vapor. This results in the nucleation enhancement as well. To verify the TPM, a simulation is performed that mimics a steady-state nucleation experiments in the thermal diffusion cloud chamber. The nucleating vapor with and without a carrier gas is simulated using two different thermostats for the monomers and clusters. The TPM proves to match the simulation results of this work and of other studies.

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

Observation of Weyl exceptional rings in thermal diffusion

A non-Hermitian Weyl equation indispensably requires a three-dimensional (3D) real/synthetic space, and it is thereby perceived that a Weyl exceptional ring (WER) will not be present in thermal diffusion given its purely dissipative nature. Here, we report a recipe for establishing a 3D parameter space to imitate thermal spinor field. Two orthogonal pairs of spatiotemporally modulated advections are employed to serve as two synthetic parameter dimensions, in addition to the inherent dimension corresponding to heat exchanges. We first predict the existence of WER in our hybrid conduction–advection system and experimentally observe the WER thermal signatures verifying our theoretical prediction. When coupling two WERs of opposite topological charges, the system further exhibits surface-like and bulk topological states, manifested as stationary and continuously changing thermal processes, respectively, with good robustness. Our findings reveal the long-ignored topological nature in thermal diffusion and may empower distinct paradigms for general diffusion and dissipation controls.

42 ENGINEERING↗
DOE OSTI2016

A hybridized discontinuous Galerkin method on mapped deforming domains

Here in this paper we present a hybridized discontinuous Galerkin (HDG) discretization for unsteady simulations of convection-dominated flows on mapped deforming domains. Mesh deformation is achieved through an arbitrary Lagrangian-Eulerian transformation with an analytical mapping. We present details of this transformation applied to the HDG system of equations, with focus on the auxiliary gradient equation, viscous stabilization, and output calculation. We discuss conditions under which optimal unsteady output convergence rates can be attained, and we show that both HDG and discontinuous Galerkin (DG) achieve these rates in advection-dominated flows. Results for scalar advection-diffusion and the Euler equations verify the implementation of the mesh motion mapping for both discretizations, show that HDG and DG yield similar results on a given mesh, and demonstrate differences in output convergence rates depending on the choice of HDG viscous stabilization. We note that such similar results bode well for HDG, which has fewer globally-coupled degrees of freedom compared to DG. A simulation of the unsteady compressible Navier-Stokes equations demonstrates again very similar results for HDG and DG and illustrates a pitfall of using steady-state adapted meshes for accurate unsteady simulations.

42 ENGINEERING↗