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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sheath formation around a dielectric droplet in a He atmospheric pressure plasma

Interactions at the interface between atmospheric pressure plasmas and liquids are being investigated to address applications ranging from nanoparticle synthesis to decontamination and fertilizer production. Furthermore, many of these applications involve activation of droplets wherein the droplet is fully immersed in the plasma and synergistically interacts with the plasma. To better understand these interactions, two-dimensional modeling of radio frequency (RF) glow discharges at atmospheric pressure operated in He with an embedded lossy dielectric droplet (tens of microns in size) was performed. The properties of the sheath that forms around the droplet were investigated over the RF cycle. The electric field in the bulk plasma polarizes the dielectric droplet while the electron drift in the external electric field is shadowed by the droplet. The interaction between the bulk and sheath electric fields produces a maximum in E/N (electric field/gas number density) at the equator on one side of the droplet where the bulk and sheath fields are aligned in the same direction and a minimum along the opposite equator. Due to resistive heating, the electron temperature T e is maximum 45° above and below the equator of the droplet where power deposition per electron is the highest. Although the droplet is, on the average, negatively charged, the charge density on the droplet is positive on the poles and negative on the equator, as the electron motion is primarily due to diffusion at the poles but due to drift at the equator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Logarithmic Superdiffusion in Two Dimensional Driven Lattice Gases

The spreading of density fluctuations in two-dimensional driven diffusive systems is marginally anomalous. Mode coupling theory predicts that the diffusivity in the direction of the drive diverges with time as (ln t){sup 2/3} with a prefactor depending on the macroscopic current-density relation and the diffusion tensor of the fluctuating hydrodynamic field equation. Here we present the first numerical verification of this behavior for a particular version of the two-dimensional asymmetric exclusion process. Particles jump strictly asymmetrically along one of the lattice directions and symmetrically along the other, and an anisotropy parameter p governs the ratio between the two rates. Using a novel massively parallel coupling algorithm that strongly reduces the fluctuations in the numerical estimate of the two-point correlation function, we are able to accurately determine the exponent of the logarithmic correction. In addition, the variation of the prefactor with p provides a stringent test of mode coupling theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Simulations of the stratocumulus-topped boundary layer with a third-order closure model

A third order closure model is proposed by Andre et al. (1982), in which the time rate of change terms, the relaxation and rapid effects for the pressure related terms, and the clipping approximation are included along with the quasi-normal closure, to study turbulence in a cloudy layer which is cooled radiatively from above. A spurious oscillation which is strongest near the inversion occurs. An analysis of the problem shows that the oscillation arises from the mean gradient and buoyancy terms of the triple moment equations; these terms are largest near the cloud top. The oscillation is physical, rather than computational. In nature the oscillation is effectively damped, by a mechanism which apparently is not included in our model. In the stably stratified layer just above the mixed layer top, turbulence can excite gravity waves, whose energy is radiated away. Because the closure assumption for the pressure terms does not take into account the transport of wave energy, the model generates spurious oscillations. Damping of the oscillations is possible by introducing diffusion terms into the triple moment equations. With a large enough choice for the diffusion coefficient, the oscillation is effectively eliminated. The results are quite sensitive to the ad hoc eddy coefficient.

Moeng, C. H.↗
NASA NTRS1974

The diffusive idealization of charged-particle transport in random magnetic fields

The uniqueness and accuracy of the equations which describe the transport of charged particles diffusing in a random magnetic field parallel to a relatively large guiding field is examined. With regard to uniqueness, it is found that the same coefficient of diffusion is obtained by three methods that have apparently led to discrepancies in previous work. With regard to accuracy, it is found that two corrections must be added to Fick's law in which the diffusive flux is proportional to the gradient of the density. Explicit expressions are given for a characteristic time and a characteristic length which describe the corrections.

Earl, J. A.↗
NASA NTRS1981

On the transport and propagation of relativistic electrons in galaxies

Analytic solutions to the time-dependent equation describing the transport of relativistic electrons that are being convected and diffused, and may be either gaining or losing energy, are described. The calculations yield a body of solutions from which specific ones can be chosen for the modelling of synchrotron radiation in the Galaxy and in other galaxies. The solutions allow for (1) spatially varying convection velocity; (2) an energy and spatially dependent diffusion coefficient; (3) adiabatic acceleration effects; (4) a term describing spatially dependent synchrotron loss and/or Fermi acceleration; (5) an injection source for the electrons that may be impulsive, steady-state, or variable-injection; (6) a spatial injection or spatially more diffuse source; (7) arbitrary energy dependence of the injection source; and (8) spatial and energy boundaries where the spatial or energy form of the convection, diffusion, and energy supply terms in the basic equation can change character.

Lerche, I.↗
DOE OSTI2013

Estimating anisotropic diffusion of neutrons near the boundary of a pebble bed random system

Due to the arrangement of the pebbles in a Pebble Bed Reactor (PBR) core, if a neutron is located close to a boundary wall, its path length probability distribution function in directions of flight parallel to the wall is significantly different than in other directions. Hence, anisotropic diffusion of neutrons near the boundaries arises. We describe an analysis of neutron transport in a simplified 3-D pebble bed random system, in which we investigate the anisotropic diffusion of neutrons born near one of the system's boundary walls. While this simplified system does not model the actual physical process that takes place near the boundaries of a PBR core, the present work paves the road to a formulation that may enable more accurate diffusion simulations of such problems to be performed in the future. Monte Carlo codes have been developed for (i) deriving realizations of the 3-D random system, and (ii) performing 3-D neutron transport inside the heterogeneous model; numerical results are presented for three different choices of parameters. These numerical results are used to assess the accuracy of estimates for the mean-squared displacement of neutrons obtained with the diffusion approximations of the Atomic Mix Model and of the recently introduced [1] Non-Classical Theory with angular-dependent path length distribution. The Non-Classical Theory makes use of a Generalized Linear Boltzmann Equation in which the locations of the scattering centers in the system are correlated and the distance to collision is not exponentially distributed. We show that the results predicted using the Non-Classical Theory successfully model the anisotropic behavior of the neutrons in the random system, and more closely agree with experiment than the results predicted by the Atomic Mix Model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Discretization Writeup for Grey Flux-Limited Radiation Diffusion

This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Hydrodynamic equations for electrons in graphene obtained from the maximum entropy principle

The maximum entropy principle is applied to the formal derivation of isothermal, Euler-like equations for semiclassical fermions (electrons and holes) in graphene. After proving general mathematical properties of the equations so obtained, their asymptotic form corresponding to significant physical regimes is investigated. In particular, the diffusive regime, the Maxwell-Boltzmann regime (high temperature), the collimation regime and the degenerate gas limit (vanishing temperature) are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Model of line preserving field line motions using Euler potentials

We consider behavior of finite magnetic field lines during reconnection processes. We portray field line motions using Euler potentials representation. Here, we propose a new insight into plasma flow fields related with magnetic reconnection. In this approach reconnection is treated as a breakage of magnetic topology, which results in deviation from the line preserving flow regime. We derive constraints and the general equations for these flows. In our approach the flux preserving flows are treated as a special case of line preserving regime. We also derive a constraint on a non-ideal term in Ohm’s Law within diffusion regions, which relates plasma flow with resistivity, and which must hold for non-reconnective diffusion. We also propose a new method of detecting magnetic reconnection. -- Highlights: •We formulate a line preserving magnetic field flow equation using Euler potentials. •We find constraints on a non-reconnective general resistivity term in Ohm’s Law. •We propose a new method of detecting magnetic reconnection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hierarchical model reduction driven by a proper orthogonal decomposition for parametrized advection-diffusion-reaction problems

This work combines the Hierarchical Model (HiMod) reduction technique with a standard Proper Orthogonal Decomposition (POD) to solve parametrized partial differential equations for the modeling of advection-diffusion-reaction phenomena in elongated domains (e.g., pipes). This combination leads to what we define as HiPOD model reduction, which merges the reliability of HiMod reduction with the computational efficiency of POD. Two HiPOD techniques are presented and assessed by an extensive numerical verification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1952

An active particle diffusion theory of flame quenching for laminar flames / Dorothy M. Simon and Frank E. Belles

An equation for quenching distance based on the destruction of chain carriers by the surface is derived. The equation expresses the quenching distance in terms of the diffusion coefficients and partial pressures of the chain carriers and gas phase molecules, the efficiency of the surface as a chain breaker, the total pressure of the mixture, and a constant which depends on the geometry of the quenching surface. Quenching distances measured by flashback for propane-air flames are shown to be consistent with the mechanism. The derived equation is used with the lean inflammability limit and a rate constant calculated from burning velocity data to estimate quenching distances for propane-air (hydrocarbon lean) flames satisfactorily.

Simon, Dorothy M↗
NASA NTRS2016

New Approaches to Inferences for Steep-Sided Domes on Venus

New mathematical approaches for the relaxation and emplacement of viscous lava domes are presented and applied to steep-sided domes on Venus. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry for two distinct scenarios. In the first scenario, dome relaxation is explored assuming a constant volume of fluid (i.e. lava) has been rapidly emplaced onto the surface. Cooling of lava is represented by a time-variable viscosity and singularities inherent in previous models for dome relaxation have been eliminated. At the onset of relaxation, bulk dynamic viscosities lie in the range between 10(exp 10) - 10(exp 16) Pa s, consistent with basaltic-andesite to rhyolitic compositions. Plausible relaxation times range from 5 to 5000 years, depending on initial lava viscosity. The first scenario, however, is only valid during the final stages of dome relaxation and does not consider the time taken for lava to be extruded onto the surface. In the second scenario, emplacement and growth of a steep-sided dome is considered when the volume of lava on the surface increases over time (i.e. time-variable volume approach). The volumetric flow rate may depend on an arbitrary power of the dome thickness, thus embracing Newtonian as well as other rheologies for describing terrestrial and planetary mass flows. The approach can be used to distinguish between basic flow rate models for fluid emplacement. The formalism results in radial expansion of a dome proportional to t(sup 1/2), consistent with the diffusive nature of the governing equation. The flow at the front is shown to thicken as the front advances for a constant rate of lava supply. Emplacement times are intimately correlated with the bulk rheology. Comparison of the theoretical profiles with the shape of a typical dome on Venus indicates that a Newtonian bulk rheology is most appropriate, consistent with prior studies. However, results here suggest a bulk dynamic viscosity of 10(exp 12) - 10(exp 13) Pa s and emplacement times of approximately 2-16 years. Both scenarios investigated give emplacement times significantly less than prior estimates and compositions consistent with basaltic andesite.

steep-sided domes↗
DOE OSTI2014

The Coal-Seq III Consortium. Advancing the Science of CO 2 Sequestration in Coal Seam and Gas Shale Reservoirs

The Coal-Seq consortium is a government-industry collaborative that was initially launched in 2000 as a U.S. Department of Energy sponsored investigation into CO2 sequestration in deep, unmineable coal seams. The consortium’s objective aimed to advancing industry’s understanding of complex coalbed methane and gas shale reservoir behavior in the presence of multi-component gases via laboratory experiments, theoretical model development and field validation studies. Research from this collaborative effort was utilized to produce modules to enhance reservoir simulation and modeling capabilities to assess the technical and economic potential for CO2 storage and enhanced coalbed methane recovery in coal basins. Coal-Seq Phase 3 expands upon the learnings garnered from Phase 1 & 2, which has led to further investigation into refined model development related to multicomponent equations-of-state, sorption and diffusion behavior, geomechanical and permeability studies, technical and economic feasibility studies for major international coal basins the extension of the work to gas shale reservoirs, and continued global technology exchange. The first research objective assesses changes in coal and shale properties with exposure to CO2 under field replicated conditions. Results indicate that no significant weakening occurs when coal and shale were exposed to CO2, therefore, there was no need to account for mechanical weakening of coal due to the injection of CO2 for modeling. The second major research objective evaluates cleat, Cp, and matrix, Cm, swelling/shrinkage compressibility under field replicated conditions. The experimental studies found that both Cp and Cm vary due to changes in reservoir pressure during injection and depletion under field replicated conditions. Using laboratory data from this study, a compressibility model was developed to predict the pore-volume compressibility, Cp, and the matrix compressibility, Cm, of coal and shale, which was applied to modeling software to enhance model robustness. Research was also conducted to improve algorithms and generalized adsorption models to facilitate realistic simulation of CO2 sequestration in coal seams and shale gas reservoirs. The interaction among water and the adsorbed gases, carbon dioxide (CO2), methane (CH4), and nitrogen (N2) in coalbeds is examined using experimental in situ laboratory techniques to comprehensively model CBM production and CO2 sequestration in coals. An equation of state (EOS) module was developed which is capable of predicting the density of pure components and mixtures involving the wet CBM gases CH4, CO2, and N2 at typical reservoir condition, and is used to inform CO2 injection models. The final research objective examined the effects adsorbed CO2 has on coal strength and permeability. This research studied the weakening or failure of coal by the adsorption of CO2 from empirically derived gas production data to develop models for advanced modeling of permeability changes during CO2 sequestration. The results of this research effort have been used to construct a new and improved model for assessing changes in permeability of coal reservoirs due CO2 injection. The modules developed from these studies and knowledge learned are applied to field validation and basin assessment studies. These data were used to assess the flow and storage of CO2 in a shale reservoir, test newly developed code against large-scale projects, and conduct a basin-oriented review of coal storage potential in the San Juan Basin. The storage potential and flow of CO2 was modeled for shale sequestration of a proprietary Marcellus Shale horizontal gas production well using COMET3 simulation software. Simulation results from five model runs indicate that stored CO2 quantities are linked to the duration of primary production preceding injection. Matrix CO2 saturation is observed to increase in each shale zone after injection with an increase in primary production, and the size of the CO2 plume is also observed to increase in size the longer initial production is sustained. The simulation modules developed around the Coal-Seq experimental work are also incorporated into a pre-existing large-scale numerical simulation model of the Pump Canyon CO2-ECBM pilot in the San Juan Basin. The new model was applied to re-history match the data set to explore the improvements made in permeability prediction against previously published data sets and to validate this module. The assessment of the new data, however, indicates that the impact of the variable Cp is negligible on the overall behavior of the coal for CO2 storage purposes. Applying these new modules, the San Juan Basin and the Marcellus Shale are assessed for their technical ECBM/AGR and CO2 storage potential and the economic potential of these operations. The San Juan Basin was divided into 4 unique geographic zones based on production history, and the Marcellus was divided into nine. Each was assessed based upon each zone’s properties, and simulations were run to assess the potential of full Basin development. Models of a fully developed San Juan Basin suggest the potential for up to 104 Tcf of CO2 storage, and 12.3 Tcf of methane recovery. The Marcellus models suggest 1,248 Tcf of CO2 storage and 924 Tcf of AGR. The economics are deemed favorable where credits cover the cost of CO2 in the San Juan Basin, and in many cases in the Marcellus, but to maximize storage potential, credits need to extend to pay the operator to store CO2.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2014

Calculation of self-diffusion coefficients in iron

On the basis of available P-V-T equation of state of iron, the temperature and pressure dependence of self-diffusion coefficients in iron polymorphs (α, δ, γ and ε phases) have been successfully reproduced in terms of the bulk elastic and expansivity data by means of a thermodynamical model that interconnects point defects parameters with bulk properties. The calculated diffusion parameters, such as self-diffusion coefficient, activation energy and activation volume over a broad temperature range (500-2500 K) and pressure range (0-100 GPa), compare favorably well with experimental or theoretical ones when the uncertainties are considered.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Solutions for Reacting and Nonreacting Viscous Shock Layers with Multicomponent Diffusion and Mass Injection

Numerical solutions are presented for the viscous shocklayer equations where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium. Also the effects of the diffusion model, surface catalyticity, and mass injection on surface transport and flow parameters are considered. The equilibrium calculations for air species using multicomponent: diffusion provide solutions previously unavailable. The viscous shock-layer equations are solved by using an implicit finite-difference scheme. The flow is treated as a mixture of inert and thermally perfect species. Also the flow is assumed to be in vibrational equilibrium. All calculations are for a 45 deg hyperboloid. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented showing the effects of the chemical models; diffusion models; surface catalyticity; and mass injection of air, water, and ablation products on heat transfer; skin friction; shock stand-off distance; wall pressure distribution; and tangential velocity, temperature, and species profiles.

Moss, J. N.↗