Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The structure of mass-loading shocks

A new two-fluid model which describes mass loading in the solar wind (e.g., the interaction of the solar wind with a cometary coma or the local interstellar medium) is presented. The self-consistent back-reaction of the mass-loaded ions is included through their effective scattering in low-frequency MHD turbulence and the invocation of a diffusive approximation. Such an approximation has the advantage of introducing self-consistent dissipation coefficients into the governing equations, thereby facilitating the investigation of the internal structure of shocks in mass-loading environments. To illustrate the utility of the new model, we consider the structure of cometary shocks in the hypersonic one-dimensional limit, finding that the incoming solar wind is slowed by both mass loading and the development of a large cometary ion pressure gradient. The shock is broadened and smoothed by the cometary ions with a thickness of the order of the cometary ion diffusion scale.

Zank, G. P.↗
NASA NTRS1997

A Three-Wave Model of the Stratosphere with Coupled Dynamics, Radiation and Photochemistry

A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral dynamical model is used to simulate the distributions of trace gases in the present day atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential height are solved in the semi-spectral dynamical model. The residual mean circulation is derived from these dynamical variables and used to advect the chemical species in the 2- D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients for chemical tracers are expressed in terms of the amplitude, frequency and growth rate of dynamical waves; local chemical loss rates; and a time constant parameterizing small scale mixing. The contributions to eddy flux are from the time varying wave amplitude (transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of the high truncation in the dynamical module (only three longest waves are resolved), the model has simulated many observed characteristics of stratospheric dynamics and distribution of chemical species including ozone. Compared with the values commonly used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this model are smaller, especially in the subtropics. It is also found that the chemical eddy diffusion has only a small effects in determining the distribution of most slow species, including ozone in the stratosphere.

Shia, Run-Lie↗
DOE OSTI2021

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗
NASA NTRS1989

Stoichiometric effects on the optical properties of LiInSe(2)

The diffuse reflectance of LiInSe(2) between 0.67 and 3.54 eV have been measured and the Kubelka-Munk theory was used to obtain the absorption coefficient from the data. The band gap in the samples is located at approximately 1.6 eV. The In and Se content of the samples was determined from atomic absorption measurements. Absorption peaks at 0.890 and 0.896 eV have been observed which are correlated, respectively, to the selenium and indium deficiencies in the samples.

Smith, Cecily J.↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗
DOE OSTI2022

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

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

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗
NASA NTRS1980

The effect of convection on the propagation of relativistic galactic electrons

The paper presents exact analytical solutions which describe steady-state transport of relativistic electrons subject to diffusion, convection, and radiation losses. The diffusion coefficient is spatially and energy dependent; the bulk convection velocity is considered to be constant and directed outward from the galactic plane. The electron spectral behavior with respect to galactic height and energy indicates that convection provides a 'break' in the steady-state energy spectrum without recourse to distribution of sources. Radio observations of NGC 891 and NGC 4631 galaxies can be used to estimate their diffusion coefficient and its energy dependence. Application of analytical solutions to our Galaxy indicates that the diffusion coefficient and outflow velocity deduced from relativistic electron measurements and radio observations of their synchrotron radiation agree with values based on the transport of the cosmic ray nucleonic component.

Lerche, I.↗
NASA NTRS1959

Theoretical and Experimental Investigation of Heat Conduction in Air, Including Effects of Oxygen Dissociation

Solutions are presented for the conduction of beat through a semi-infinite gas medium having a uniform initial temperature and a constant boundary temperature. The coefficients of thermal conductivity and diffusivity are treated as variables, and the solutions are extended to the case of air at temperatures where oxygen dissociation occurs. These solutions are used together with shock-tube measurements to evaluate the integral of thermal conductivity for air as a function of temperature.

Hansen, C. Frederick↗
DOE OSTI2023

Inner Belt Electron Decay Timescales: A Comparison of Van Allen Probes and DREAM3D Losses Following the June 2015 Storm

Abstract NASA's Van Allen Probes observed significant, long‐lived fluxes of inner belt electrons up to ∼1 MeV after geomagnetic storms in March and June 2015. Reanalysis of Magnetic Electron Ion Spectrometer (MagEIS) data with improved background correction showed a clearer picture of the relativistic electron population that persisted through 2016 and into 2017 above the Fennell et al. (2015, https://doi.org/10.1002/2014gl062874 ) limit. The intensity and duration of these enhancements allow estimation of decay timescales for comparison with simulated decay rates and theoretical lifetimes. We compare decay timescales from these data and DREAM3D simulations based on them using geomagnetic activity‐dependent pitch angle diffusion coefficients derived from plasmapause‐indexed wave data (Malaspina et al., 2016, https://doi.org/10.1002/2016gl069982 , 2018, https://doi.org/10.1029/2018gl078564 ) and phase space densities derived from MagEIS observations. Simulated decay rates match observed decay rates more closely than the theoretical lifetime due to significantly nonequilibrium pitch angle distributions in simulation and data. We conclude that nonequilibrium effects, rather than a missing diffusion or loss process, account for observed short decay rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

MODELLING FIBER ORIENTATION DURING ADDITIVE MANUFACTURING COMPRESSION MOLDING PROCESSES

The production of high-performance thermoplastic composites reinforced with short carbon fibers can be achieved by a novel “additive manufacturing-compression molding” technique. An advantage of such combination is two-fold: controlled fiber orientation in additive manufacturing and less void content by compression molding. In this study, a computational fluid dynamics model has been developed to predict the behavior of printed layers during fiber-reinforced thermoplastic extrusion and subsequent compression molding. The fiber orientation was modelled with a simple quadratic closure model. The interaction between the fibers is included using a rotary diffusion coefficient which becomes significant in concentrated regimes. Finally, the second order orientation tensor is coupled with the momentum equation as an anisotropic part of the stress term. The effect of processing parameters on the behavior of printed layers was investigated to determine the favorable printing scenarios. The developed numerical model enables design of high-performance composites with tunable mechanical properties.

Seta, Berin↗
NASA NTRS2002

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗