Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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
NASA NTRS1991

Venus ionospheric tail rays - Spatial distributions and interplanetary magnetic field control

The overall properties of Venus ionospheric tail rays (such as density, spatial extent, and distribution) and their relationship to the draped magnetic field configuration behind the planet were investigated using measurements obtained by the Pioneer Venus Orbiter Langmuir probe, a magnetometer, and a plasma-wave detector. The results suggest that tail rays are a normal feature of the steady solar wind interaction with Venus and are not generally associated with a central tail plasma sheet. The statistics of the tail rays occurrence point toward the existence of a distributed terminator ionosphere source, consistent with findings of Brace et al. (1990).

Ong, M.↗
DOE OSTI2020

Low-mass Group Environments Have No Substantial Impact on the Circumgalactic Medium Metallicity

We explore how environment affects the metallicity of the circumgalactic medium (CGM) using 13 low-mass galaxy groups (two to five galaxies) at &lang;z{sub abs}&rang;=0.25 identified near background quasars. Using quasar spectra from the Hubble Space Telescope/COS and from Keck/HIRES or the Very Large Telescope/UVES, we measure column densities of or determine limits on CGM absorption lines. We use a Markov Chain Monte Carlo approach with Cloudy to estimate metallicities of cool (T ~ 10{sup 4} K) CGM gas within groups and compare them to CGM metallicities of 47 isolated galaxies. Both group and isolated CGM metallicities span a wide range (-2 < [Si/H] < 0), where the mean group (-0.54 ± 0.22) and isolated (-0.77 ± 0.14) CGM metallicities are similar. Group and isolated environments have similar distributions of H I column densities as a function of impact parameter. However, contrary to isolated galaxies, we do not find an anticorrelation between H I column density and the nearest group galaxy impact parameter. We additionally divided the groups by member luminosity ratios (i.e., galaxy–galaxy and galaxy–dwarf groups). While there was no significant difference in their mean metallicities, a modest increase in sample size should allow one to statistically identify a higher CGM metallicity in galaxy–dwarf groups compared to galaxy–galaxy groups. We conclude that either environmental effects have not played an important role in the metallicity of the CGM at this stage and expect that this may only occur when galaxies are strongly interacting or merging or that some isolated galaxies have higher CGM metallicities due to past interactions. Thus, environment does not seem to be the cause of the CGM metallicity bimodality.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Collaborative Research Proposal: Time Resolved Optical Emissions Spectroscopy and Laser Induced Fluorescence Spectroscopy of Nanosecond Pulsed Discharges in a Gas-Liquid Water Film Reactor

Electrical discharge plasma formed in contact with liquid water is of interest for a wide range of applications in chemical, biomedical, agricultural, electrical, and materials science and engineering. Such plasma reactors are of very timely importance since they also have significant disinfection capability by inactivating bacteria, viruses and other pathogens. Many types of plasma sources including those driven by AC, DC, RF, microwave, and pulsed electrical power supplies coupled to a wide range of different electrode configurations and reactor designs have been developed and explored which contact the plasma with liquid water. Recent roadmaps have recommended that further work is needed to develop our understanding of the fundamental chemical and physical process which occur at the interface of non-thermal plasma with liquid water solutions in order to advance this large diversity of applications which ultimately depend upon efficient production of key reactive chemical species. There is a wide range of interacting factors that affect the chemical reactions that occur in the plasma, in the liquid phase, and at the interface. These factors ultimately govern the key reactive chemical species formed and used in the various applications and they include a) the reactor design and input parameters, b) the discharge and transport processes, c) the plasma properties, and d) the resulting chemical reactions. For pulsed discharges, the power supply design and output parameters control the applied voltage, frequency, rise time, and width (duration) of the applied pulses. The reactor design involves specification of the gas-liquid contacting methods, electrode gap distance, reactor volume and shape, and gas and liquid flow rates achievable. The gas and liquid compositions as well as the liquid properties such as pH and conductivity are also of key importance in determination of the resulting chemical reactions. In addition, the important plasma properties include plasma gas temperature, electron density, electron energy (distribution), and size of the plasma channels which are all affected by the discharge and transport properties. Many studies have focused on specific aspects of these various processes. Of particular important and relevance to the proposed work is the utilization of nanosecond pulses to generate plasma in gas-liquid systems, liquid bubbles, underwater, and in gases. Recent advances in nanosecond pulses provide significant advantages in utilization with liquid water. For example, fast rise time and short pulses are less sensitive to water conductivity, such short pulses may provide advantages in fast temporal quenching of the plasma, and fundamental analysis of pulse properties, including pulse shape and width, may be facilitated by investigation of single filamentary fast pulses. Many studies have also dealt with the role of the gas composition on formation of reactive oxygen species (ROS) (i.e., hydroxyl radicals – ·OH, hydrogen peroxide – H 2 O 2 , various atomic oxygen species, ozone-O 3 , hydroperoxyl radicals HO 2 ·) and reactive nitrogen species (RNS) (i.e., nitrogen oxides – NO, NO 2 , N 2 O - collectively termed NO X , nitrite-NO 2 - , nitrate-NO 3 - , peroxynitrite-ONOO - ). The hydroxyl radical is the critical species in many chemical oxidation reactions for chemical degradation of toxic compounds in water and gases and for synthesis of some compounds. The mixture of various nitrogen oxide species is important for many biochemical and biological processes involved in biomedical and agricultural applications including disinfection and fertilizer production. The present proposal focuses on determination of the effects of time resolved electron density and hydroxyl radicals on plasma chemical reactions through collaboration with the Princeton Collaborative Low Temperature Plasma Research Facility (PCRF) at the Princeton Plasma Physics Laboratory (PPPL). In order to further investigate the role of the plasma generated electrons and hydroxyl radicals on the overall formation of the key species including hydrogen peroxide, hydroxyl radicals, and nitrogen oxides, the proposed work, thus seeks to determine high resolution time resolved electron density by optical emissions spectroscopy and time resolved hydroxyl radicals using laser induced fluorescence in the nanosecond discharge reactor. The combination of data on electron density and hydroxyl radical concentration will be utilized in the present work to more fully characterize the chemical reaction processes in this system and to advance the design, development, and operation of such chemical reactors for a wide range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Monte Carlo event generator for initial conditions of conserved charges in nuclear geometry

At top collider energies where baryon stopping is negligible, the initial state of heavy-ion collisions is overall charge neutral and predominantly composed of gluons. Nevertheless, there can also be significant local fluctuations of the baryon number, strangeness, and electric charge densities about zero, perturbatively corresponding to the production of quark/antiquark pairs. These previously ignored local charge fluctuations can permit the study of charge diffusion in the quark-gluon plasma (QGP), even at top collider energies. In this paper we present a new model denoted ICCING (initial conserved charges in nuclear geometry) which can reconstruct the initial conditions of conserved charges in the QGP by sampling a (g→q¯q) splitting probability over the initial energy density. We find that the new charge distributions generally differ from the bulk energy density; in particular, the strangeness distribution is significantly more eccentric than standard bulk observables and appears to be associated with the geometry of hot spots in the initial state. The new information provided by these conserved charges opens the door to studying a wealth of new charge- and flavor-dependent correlations in the initial state and ultimately the charge transport parameters of the QGP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2008

Gamma-Ray Burst Afterglows as Probes of Environment and Blast Wave Physics. II. The Distribution of rho and Structure of the Circumburst Medium

We constrain blast wave parameters and the circumburst media ofa subsample of 10 BeppoSAX gamma-ray bursts (GRBs). For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical, and NIR afterglow data. The spectral fits have been done in count space and include the effects ofmetallicity, and are compared with the previously reported optical and X-ray temporal behavior. Using the blast wave model and some assumptions which include on-axis viewing and standard jet structure, constant blast wave energy, and no evolution of the microphysical parameters, we find a mean value ofp for the sample as a whole of 9.... oa -0.003.0" 2 a_ statistical analysis of the distribution demonstrates that the p-values in this sample are inconsistent with a single universal value forp at the 3 _ level or greater, which has significant implications for particle acceleration models. This approach provides us with a measured distribution ofcircumburst density structures rather than considering only the cases of k ----0 (homogeneous) and k - 2 (windlike). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly windlike. The fifth source has a value of 0 < k < 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗
DOE OSTI2025

Programmable Cryogenic Memory in a Ge/GeSi Heterostructure

Programmable memory components that operate optimally at cryogenic temperatures are essential for cryogenic computing architectures that seek to implement computing-in-memory. In this work, we demonstrate highly programmable memory in a Ge/GeSi heterostructure field-effect transistor (HFET). To operate, the HFET is gated to introduce positive carriers within the Ge quantum well, creating a high-conductance state. We show that this device can be set to a low-conductance state by sweeping a negative bias on the device drain, and reset it to its high-conductance state by sweeping a more positive bias on the device gate, thereby creating memory. We then determine that the device can be programmed within a 103 range of conductances using either the SET or the RESET operation. We propose that memory is achieved through charge trapping as carriers tunnel out of the quantum well, and that altering the density and spatial distribution of carriers modulates the device conductance. This mechanism exhibits endurance over 1000 cycles at temperatures ≤ 25 K, suggesting that the carrier traps are located at the oxide-semiconductor interface. As a first demonstration of programmable conductance in a Ge/GeSi HFET, this work highlights the potential of group-IV HFETs to perform as analog cryogenic memory components.

cryogenic memory↗
DOE OSTI2022

Understanding the Role of Order in Y‐Series Non‐Fullerene Solar Cells to Realize High Open‐Circuit Voltages

Abstract Non‐fullerene acceptors (NFAs) as used in state‐of‐the‐art organic solar cells feature highly crystalline layers that go along with low energetic disorder. Here, the crucial role of energetic disorder in blends of the donor polymer PM6 with two Y‐series NFAs, Y6, and N4 is studied. By performing temperature‐dependent charge transport and recombination studies, a consistent picture of the shape of the density of state distributions for free charges in the two blends is developed, allowing an analytical description of the dependence of the open‐circuit voltage V OC on temperature and illumination intensity. Disorder is found to influence the value of the V OC at room temperature, but also its progression with temperature. Here, the PM6:Y6 blend benefits substantially from its narrower state distributions. The analysis also shows that the energy of the equilibrated free charge population is well below the energy of the NFA singlet excitons for both blends and possibly below the energy of the populated charge transfer manifold, indicating a down‐hill driving force for free charge formation. It is concluded that energetic disorder of charge‐separated states has to be considered in the analysis of the photovoltaic properties, even for the more ordered PM6:Y6 blend.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of Ligand Polarity on the Internal Dipoles and Ferroelectric Distortion in BaTiO 3 Nanocubes

Here, surface adsorbates and surrounding matrix species have been demonstrated to affect the properties of nanoscale ferroelectrics and nanoscale ferroelectric composites; potentially counteracting performance losses that can occur in small particle sizes. In this work, the effects of nonpolar oleic acid (OA) and polar tetrafluoroborate (BF 4 - ) ligand capping on the surface of various sizes of BaTiO 3 nanocubes have been investigated with combined neutron diffraction and neutron pair distribution function (PDF), density functional theory (DFT), and ab initio molecular dynamics (AIMD) methods. The low real space PDF region provides an unobstructed view of rhombohedral (split short and long) Ti-O distances in BaTiO 3 nanocubes, mimicking the well-established order-disorder local structure found in bulk BaTiO 3 . Interestingly, the intermediate-range order in nanocubes is found to be orthorhombic, rather than tetragonal. It is concluded that polar ligands adsorbed at BaTiO 3 surfaces stabilize the correlation length scale of local rhombohedral distortions in ferroelectric nanoparticles relative to nonpolar ligands.

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

Elucidation of Design Criteria for V‐based Redox Mediators: Structure‐Function Relationships that Dictate Rates of Heterogeneous Electron Transfer

Redox mediators are attractive solutions for addressing the stringent kinetic stipulations required for efficient energy conversion processes. In this work, we compare the electrochemical properties of four vanadium complexes, namely [V(acac) 3 ], [V 6 O 7 (OMe) 12 ], [ n Bu 4 N] 3 [V 6 O 13 (TRIS NO2 ) 2 ], and [ n Bu 4 N] 5 [V 18 O 46 (NO 3 )] in non-aqueous solutions on glassy carbon electrodes. The goal of this study is to investigate the electron transfer kinetics and diffusivity of these compounds under identical experimental conditions to develop an understanding of structure-function relationships that dictate the physicochemical properties of vanadium oxide assemblies. Complex selection was dictated by two criteria – (1) nuclearity of the transition metal complexes (2) distribution of electron density in the native electronic configuration. In conclusion, our analyses establish that electronic communication between metal centers significantly impacts charge transfer kinetics of these vanadium-based compounds.

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

Explicitly determined sea ice emissivity and emission temperature over the Arctic for surface–sensitive microwave channels

Data assimilation of satellite microwave measurements is one of the important keys to improving weather forecasting over the Arctic region. However, the use of surface-sensitive microwave-sounding channel measurements for data assimilation or retrieval has been limited, especially during winter, due to the poorly constrained sea ice emissivity. In this study, aiming at more use of those channel measurements in the data assimilation, we propose an explicit method for specifying the surface radiative boundary conditions (namely emissivity and emitting layer temperature of snow and ice). These were explicitly determined with a radiative transfer model for snow and ice and with snow/ice physical parameters (i.e., snow/ice depths and vertical distributions of temperature, density, salinity, and grain size) simulated from the thermodynamically-driven snow/ice growth model. Here we conducted 1-D Var experiments in order to examine whether this approach can help to use the surface-sensitive microwave temperature channel measurements over the Arctic sea ice region for data assimilation. Results show that (1) the surface-sensitive microwave channels can be used in the 1-D Var retrieval, and (2) the specification of the radiative boundary condition at the surface using the snow/sea ice emission model can significantly improve the atmospheric temperature retrieval, especially in the lower troposphere (500 hPa to surface). The successful retrieval suggests that useful information can be extracted from surface-sensitive microwave-sounding channel radiances over sea ice surfaces through the explicit determination of snow/ice emissivity and emitting layer temperature.

1-D Var retrieval↗
DOE OSTI2025

Interstitial Atoms and the Frustrated and Allowed Structural Transitions Principle: Tunability in the Electronic Structure of AuCu 3 ‐type Frameworks in Dy 4 T 1− x Ga 12 (T = Ag, Ir)

In this Article, we explore how the chemical pressure (CP) features of an intermetallic phase may provide opportunities to couple perturbations in electron count with the stabilization of the underlying geometrical structure. AuCu 3 ‐type LnGa 3 (Ln = lanthanide or group 3 metal) phases contain octahedral cavities of negative CP held open by overly compressed Ln–Ga interactions, leading to a series of transition metal‐stuffed derivatives. We present new additions to this family with the synthesis and crystal structures of Dy 4 T 1−x Ga 12 with (T, x) = (Ag, 0.29) and (Ir, 0.15), adopting Y 4 PdGa 12 ‐type superstructures of the AuCu 3 ‐type. density functional Ttheory (DFT)‐CP calculations, when adjusted to avoid dipolar CP features, affirm that T atom incorporation provides a mechanism for the relief of packing tensions, while electronic density of states distributions illustrate that the T atoms serve largely as electron or hole donors to the band structure, as needed for them to attain d 10 configurations. The maximum obtainable value for x may be limited by a mismatch between the Fermi energy and pseudogap, in line with the balance of factors envisioned by the frustrated and allowed structural transitions principle. Furthermore, trends in resistivity measurements on T = Ir, Pd, and Ag compounds are interpretable in terms of the varying degrees of disorder arising from x < 1.0.

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

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

High Density of Strong yet Deformable Intermetallic Nanorods Leads to an Excellent Room Temperature Strength-Ductility Combination in a High Entropy Alloy

This paper introduces a new microstructural template for high entropy alloys (HEAs), where the face centered cubic (FCC) complex concentrated solid solution is reinforced with a high density of strong, yet deformable, nanorods of an ordered multi-component intermetallic L12 compound. Thermodynamic modeling has been employed to design this HEA with a large L12 volume fraction. Thermo-mechanical processing by isothermal annealing of the conventionally processed bulk cold-rolled alloy directly at precipitation temperatures, has been applied to produce a high density of uniformly distributed L12 nanorods within refined FCC grains, resulting from concomitant recrystallization and discontinuous precipitation processes. The nanorod morphology of the discontinuous L12 product has been established from three-dimensional atom probe tomography. The refined grains result in a complete coverage of the microstructure with discontinuously precipitated intermetallic nanorods. This nanorod strengthened HEA exhibits an exceptionally high room temperature yield strength of ~1630 MPa, good tensile ductility of ~15%, and an ultimate tensile strength of ~1720 MPa. Furthermore, a single L12 phase alloy, melted based on the precipitate composition in the two-phase FCC + L12 HEA, exhibits very high compressive deformability and strain hardenability, unusual for ordered intermetallic compounds. These results open a new strategy for design of fine-grained microstructures strengthened via ordered intermetallic phases, exploiting the beneficial effects of discontinuous precipitation, for achieving very high room temperature tensile strengths while maintaining good ductility.

Gwalani, Bharat↗
DOE OSTI2024

Operando visualization of porous metal additive manufacturing with foaming agents through high-speed x-ray imaging

Porous metals find extensive applications in soundproofing, filtration, catalysis, and energy-absorbing structures, thanks to their unique internal pore structure and high specific strength. In recent years, there has been an increasing interest in fabricating porous metals using additive manufacturing (AM), leveraging its unique advantages, including improved design freedom, spatial material control, and cost-effective small-batch production. In this study, we conducted pioneering operando visualization of AM porous metal using a laser powder bed fusion (L-PBF) setup combined with a high-speed synchrotron x-ray imaging system. Single track printing experiments using Ti6Al4V (Ti64) combined with titanium hydride (TiH 2 ) and sodium carbonate (Na 2 CO 3 ) as foaming agents, with varying mixing ratios were performed under different processing conditions. Here. the results elucidate the dynamic development of porosity formation. The average pore size is significantly influenced by the particle size of foaming agents when pore coalescence is absent. For all foaming agent content tested in the current study, the number of pores is found to be more sensitive to changes in laser power than in laser scanning speed. Increasing linear energy density (increasing laser power or reducing laser scanning speed) promotes the foaming agent activation thereby porosity formation. However, high linear energy density skews pore distribution towards the surface despite forming deeper melt pools. In addition, the impact of additional factors including foaming agent's laser absorptivity and decomposition kinetics with respect to AM time scales should be carefully considered to avoid ineffective activation of foaming agents during the AM of porous metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanoadaptive strain and functional osseointegration of dental implants in rats

Spatiotemporal implant-bone biomechanics and mechanoadaptive strains in peri-implant tissue are poorly understood. Physical and chemical characteristics of an implant-bone complex (IBC) were correlated in three-dimensional space (along the length and around a dental implant) to gather insights into time related integration of the implant with the cortical portion of a jaw bone in a rat. Rats (N = 9) were divided into three experimental groups with three rats per time point; 3-, 11-, and 24-day. All rats were fed crumbled hard pellets mixed with water (soft-food diet) for the first 3 days followed by a hard-food diet with intact hard-food pellets (groups of 11- and 24-day only). Biomechanics of the IBCs harvested from rats at each time point was evaluated by performing mechanical testing in situ in tandem with X-ray imaging. The effect of physical association (contact area) of a loaded implant with adapting peri-implant tissue, and resulting strain within was mapped by using digital volume correlation (DVC) technique. The IBC stiffness at respective time points was correlated with mechanical strain in peri-implant tissue. Results illustrated that IBC stiffness at 11-day was lower than that observed at 3-day. However, at 24-day, IBC stiffness recovered to that which was observed at 3-day. Furthermore, correlative microscopy and spectroscopy illustrated that the lower IBC stiffness was constituted by softer and less mineralized peri-implant tissue that contained varying expressions of osteoconductive elements. Lower IBC stiffness observed at 11-day was constituted by less mineralized peri-implant tissue with osteoconductive elements that included phosphorus (P) which was co-localized with higher expression of zinc (Zn), and lower expression of calcium (Ca). Higher IBC stiffness at 24-day was constituted by mineralized peri-implant tissue with higher expressions of osteoconductive elements including Ca and P, and lower expressions of Zn. These spatiotemporal correlative maps of peri-implant tissue architecture, heterogeneous distribution of mineral density, and elemental colocalization underscore mechanoadaptive physicochemical properties of peri-implant tissue that facilitate functional osseointegration of an implant. These results provided insights into 1) plausible "prescription" of mechanical loads as an osteoinductive "therapeutic dose" to encourage osteoconductive elements in the peri-implant tissue that would facilitate functional osseointegration of the implant; 2) a "critical temporal window" between 3 and 11 days, and perhaps it is this acute phase during which key candidate regenerative molecules can be harnessed to accelerate osseointegration of an implant under load.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Defect formation and bending properties in graphite under He atom implantation investigated by molecular dynamics method

The formation and evolution of radiation defects and their effects on bending property of a single graphite under low energy helium atom implantation have been studied at atomic scale through molecular dynamics simulations. Various simple interstitial and vacancy defects are observed during the implantation process, which is a function of incident energy and cross section of helium-carbon interaction. Large defects are observed after continuous helium atom implantations. The effects of incident energy on the density, size and distribution of defects are analyzed in detail. Furthermore, the kinetic evolution process of radiation defects at different annealing temperatures have also been simulated, including the migration, aggregation, transformation, coalescence, and self-healing. Here, the large defect clusters, crossing several carbon layers, have been confirmed to be the primary reasons that affect the bending property of graphite. All these results provide a new understanding for further applications of graphite in nuclear reactors.

Atomic implantation↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2023

Pore-scale simulation of multiphase flow and reactive transport processes involved in geologic carbon sequestration

Multiphase flow and reactive transport are two essential physicochemical processes that govern the effectiveness of geological carbon sequestration (GCS). The interaction and feedback among different phases and components during intricate physicochemical processes hold great significance in understanding CO 2 sequestration. Pore-scale simulations can account for multiphase flow and reactive transport processes in porous media and obtain spatial distributions of parameters (density, velocity, concentration, etc.) in the pore space as well as their temporal evolutions. This proves especially valuable considering that experiments can be hindered by constraints in spatial and temporal resolution. The comprehensive insights garnered from pore-scale research can be leveraged for continuum modeling using the representative elementary volume (REV) concept. In this contribution, four sequential mechanisms of CO 2 -brine-rock interaction in three zones delineated by CO 2 saturation are elaborated to elucidate complicated physicochemical processes involved in GCS, which are followed by general descriptions of mathematical equations and pore-scale numerical methods. In addition, as interested and commonly encountered processes, leakage risks during GCS and CO 2 -enhanced oil recovery (CO 2 -EOR) processes are presented. The existing challenges and future directions are discussed for both the performance of the pore-scale models and the current gaps in the field of GCS. Importantly, we expect that this review will prove beneficial for researchers interested in pore-scale simulations, GCS, and related disciplines.

58 GEOSCIENCES↗