Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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

The Cosmic Ultraviolet Baryon Survey (CUBS) – III. Physical properties and elemental abundances of Lyman-limit systems at z < 1

ABSTRACT We present a systematic investigation of physical conditions and elemental abundances in four optically thick Lyman-limit systems (LLSs) at z = 0.36–0.6 discovered within the Cosmic Ultraviolet Baryon Survey (CUBS). Because intervening LLSs at z < 1 suppress far-UV (ultraviolet) light from background QSOs, an unbiased search of these absorbers requires a near-UV-selected QSO sample, as achieved by CUBS. CUBS LLSs exhibit multicomponent kinematic structure and a complex mix of multiphase gas, with associated metal transitions from multiple ionization states such as C ii, C iii, N iii, Mg ii, Si ii, Si iii, O ii, O iii, O vi, and Fe ii absorption that span several hundred km s−1 in line-of-sight velocity. Specifically, higher column density components (log N(H i)/cm−2≳ 16) in all four absorbers comprise dynamically cool gas with $\langle T \rangle =(2\pm 1) \times 10^4\,$K and modest non-thermal broadening of $\langle b_\mathrm{nt} \rangle =5\pm 3\,$km s−1. The high quality of the QSO absorption spectra allows us to infer the physical conditions of the gas, using a detailed ionization modelling that takes into account the resolved component structures of H i and metal transitions. The range of inferred gas densities indicates that these absorbers consist of spatially compact clouds with a median line-of-sight thickness of $160^{+140}_{-50}$ pc. While obtaining robust metallicity constraints for the low density, highly ionized phase remains challenging due to the uncertain $N\mathrm{(H\, {\small I})}$, we demonstrate that the cool-phase gas in LLSs has a median metallicity of $\mathrm{[\alpha /H]_{1/2}}=-0.7^{+0.1}_{-0.2}$, with a 16–84 percentile range of [α/H] = (−1.3, −0.1). Furthermore, the wide range of inferred elemental abundance ratios ([C/α], [N/α], and [Fe/α]) indicate a diversity of chemical enrichment histories. Combining the absorption data with deep galaxy survey data characterizing the galaxy environment of these absorbers, we discuss the physical connection between star-forming regions in galaxies and diffuse gas associated with optically thick absorption systems in the z < 1 circumgalactic medium.

Zahedy, Fakhri S. (ORCID:0000000178692551)↗
DOE OSTI2022

Cool outflows in MaNGA: a systematic study and comparison to the warm phase

ABSTRACT This paper investigates the neutral gas phase of galactic winds via the Na i Dλλ5890, 5895Å feature within z ∼ 0.04 MaNGA galaxies and directly compares their incidence and strength to the ionized winds detected within the same parent sample. We find evidence for neutral outflows in 127 galaxies (∼5 per cent of the analysed line-emitting sample). ${\rm Na\, \small {\rm I}\ D}$ winds are preferentially seen in galaxies with dustier central regions and both wind phases are more often found in systems with elevated star formation rate (SFR) surface densities, especially when there has been a recent upturn in the star formation activity according to the $\rm SFR_{5 \, Myr}/SFR_{800 \, Myr}$ parameter. We find the ionized outflow kinematics to be in line with what we measure in the neutral phase. This demonstrates that, despite their small contributions to the total outflow mass budget, there is value to collecting empirical measurements of the ionized wind phase to provide information on bulk motion in the outflow. Depending on dust corrections applied to the ionized gas diagnostics, the neutral phase has ∼1.2−1.8 dex higher mass outflow rates ($\dot{M}_{\rm out}$), on average, compared to the ionized phase. We quantify scaling relations between $\dot{M}_{\rm out}$ and the strengths of the physical wind drivers (SFR, LAGN). Using a radial–azimuthal stacking method, and by considering inclination dependencies, we find results consistent with biconical outflows orthogonal to the disc plane. Our work complements other multiphase outflow studies in the literature that consider smaller samples, more extreme objects, or proceed via stacking of larger samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

CO and [C ii] line emission of molecular clouds: the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic 12 CO, 13 CO, and [C ii] emission maps of molecular cloud (MC) simulations from the SILCC-Zoom project. We present radiation, magnetohydrodynamic zoom-in simulations of individual clouds, both with and without radiative stellar feedback, forming in a turbulent multiphase interstellar medium following on-the-fly the evolution of e.g. H 2 , CO, and C + . We introduce a novel post-processing routine based on cloudy which accounts for higher ionization states of carbon due to stellar radiation in H ii regions. Synthetic emission maps of [C ii] in and around feedback bubbles show that the bubbles are largely devoid of [C ii], as recently found in observations, which we attribute to the further ionization of C+ into C 2+ . For both 12 CO and 13 CO, the cloud-averaged luminosity ratio, $L_\rm {CO}/L_\rm {[C\, \small {II}]}$, can neither be used as a reliable measure of the H 2 mass fraction nor of the evolutionary stage of the clouds. We note a relation between the $I_\rm {CO}/I_\rm {[C\, \small {II}]}$ intensity ratio and the H 2 mass fraction for individual pixels of our synthetic maps. The scatter, however, is too large to reliably infer the H 2 mass fraction. Finally, the assumption of chemical equilibrium overestimates H2 and CO masses by up to 150 and 50 per cent, respectively, and $L_\rm {CO}$ by up to 60 per cent. The masses of H and C + would be underestimated by 65 and 30 per cent, respectively, and $L_\rm {[C\, \small {II}]}$ by up to 35 per cent. Hence, the assumption of chemical equilibrium in MC simulations introduces intrinsic errors of a factor of 2 in chemical abundances, luminosities, and luminosity ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

QUOKKA-based understanding of outflows (QED) – III. Outflow loading and phase structure as a function of galactic environment

ABSTRACT We present results from a suite of 3D high-resolution hydrodynamic simulations of supernova-driven outflows from galactic disc regions with a range of gas surface density, metallicity, and supernova scale height. We use this suite to quantify how outflow properties – particularly the loading factors for mass, metallicity, and energy – vary with these parameters. We find that the winds fall into three broad categories: steady and hot, multiphase and moderately bursty, and cool and highly bursty. The first of these is characterized by efficient metal and energy loading but weak mass loading, the second by moderate loading of mass, metals, and energy, and the third by negligible metal and energy loading but substantial mass loading. The most important factor in determining the kind of wind a galaxy will produce is the ratio of supernova to gas scale heights, with the latter set by a combination of supernova rate, metallicity-dependent cooling rate, and the gravitational potential. These often combine in counterintuitive ways – for example increased cooling causes cold clouds to sink into the galactic midplane more rapidly, lowering the volume-filling factor of dense gas and making the environment more favourable for strong winds. Our findings suggest that the nature of galactic winds is likely highly sensitive to phenomena such as runaway stars occurring at a large height and dense gas and are poorly captured in most simulations, and that metal loading factors for type Ia supernovae may be substantially larger than those for type II, with important implications for galactic chemical evolution.

Vijayan, Aditi (ORCID:0000000277142379)↗
DOE OSTI2016

Understanding Hydraulic Fracturing: A Multi-Scale Problem

Despite the impact that hydraulic fracturing has had on the energy sector, the physical mechanisms that control its efficiency and environmental impacts remain poorly understood in part because the length scales involved range from nano-meters to kilo-meters. We characterize flow and transport in shale formations across and between these scales using integrated computational, theoretical, and experimental efforts. At the field scale, we use discrete fracture network modeling to simulate production at a well site whose fracture network is based on a site characterization of a shale formation. At the core scale, we use triaxial fracture experiments and a finite-element discrete-element fracture propagation model with a coupled fluid solver to study dynamic crack propagation in low permeability shale. We use lattice Boltzmann pore-scale simulations and microfluidic experiments in both synthetic and real micromodels to study pore-scale flow phenomenon such as multiphase flow and mixing. A mechanistic description and integration of these multiple scales is required for accurate predictions of production and the eventual optimization of hydrocarbon extraction from unconventional reservoirs.

03 NATURAL GAS↗
DOE OSTI2019

Effect of Laser-Matter Interaction on Molten Pool Flow and Keyhole Dynamics

Laser-induced keyholing occurs in additive manufacturing and welding processes, but the keyhole dynamics have not been well understood. A multiphase and multiphysics numerical model is used to predict the keyhole shapes recorded in the experimental observations and to predict transient and nonuniform distributions of laser absorption, temperature, and flow velocity in the process. When compared against data from a state-of-the-art dynamic x-ray radiography technique, good agreement is found for the keyhole shapes and fluctuation of the gas-liquid interface, thereby validating the simulation method. A detailed discussion is then given to elucidate the effects of laser absorption on the dynamic behavior of the front and rear keyhole walls. In conclusion, a quantitative comparison of different driving forces on the keyhole is also given to evaluate their significance to the keyhole dynamics.

42 ENGINEERING↗
DOE OSTI2020

Shock compression of vanadium at extremes: Theory and experiment

Here, the equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50 000 K, for densities between ρ = 3 and 15 g / cm 3 , with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Shock melting and the hcp-bcc phase boundary of Mg under dynamic loading

The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enhanced ferroelectric and piezoelectric properties of BCT-BZT at the morphotropic phase boundary driven by the coexistence of phases with different symmetries

The discovery of lead-free piezoelectric materials is crucial for future information and energy storage applications. Enhanced piezoelectric and other physical properties are commonly observed near the morphotropic phase boundary (MPB) composition of ferroelectric solid solutions. The (1 – x)Ba(Zr 0.2 Ti 0.8 )O 3 -x(Ba0.7Ca0.3)TiO 3 (BZT-xBCT) system exhibits a large electromechanical response around its MPB region at x = 0.5. Here, we report experimental and theoretical results of BZT-xBCT over a wide composition range (0.3 ≤ x ≤ 1.0). X-ray diffraction and Raman spectroscopy studies indicate a composition-induced structural phase transition from a rhombohedral (R3m) phase at x ≤ 0.4 to a tetragonal (P4mm) phase at x ≥ 0.6 through a multiphase coexistence region at 0.45 ≤ x ≤ 0.55 involving orthorhombic + tetragonal (Amm2 + P4mm) phases. First-principles calculations elucidate the phase competition in the coexistence region. The critical composition (x = 0.5) displays enhanced dielectric, ferroelectric, and piezoelectric properties, where notably d 33 ~ 320 pC/N. This paper provides clear evidence of Amm2+P4mm crystallographic phases in the MPB region, which is responsible for the improved functional properties.

MPB↗
DOE OSTI2022

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, we theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Quaternary phase diagrams of spinel Li$_y\square_{1–y}$Mn$_x$Ni$_{2–x}$O$_4$ and composite cathode voltages for concentration gradient materials

Core-shell coating structures and concentration gradient materials may enhance Li-ion battery performance by integrating advantages of core and shell components without introducing unfavorable problems associated with general coatings. The fundamental thermodynamic properties of concentration gradient composite materials are complex due to the multicomponent nature of the problem. Here, we systematically study the thermodynamics of ordering and phase separation in the quaternary spinel Li$_y\square_{1–y}$Mn$_x$Ni$_{2–x}$O$_4$ ($\square$ means vacancy) system by density functional theory calculations, together with the coupled cluster expansion method with interactions within and between (Li/$\square$) and (Mn/Ni) sublattices. On the basis of coupled cluster expansion interactions and Monte Carlo simulations, we calculate quaternary phase diagrams as a function of temperature as well as voltage profiles of single ordered phases and multiphase composite structures. The phase diagram and voltage results are in good agreement with available experimental observations. We also predict a stable high-voltage ordered compound LiMnNiO 4 , with a very high delithiation voltage of 4.76 V. For the composite (Mn-rich + Ni-rich) cathode materials, the voltage profiles show combinations of plateaus from each component compound. The computational strategy of combining quaternary phase diagrams with voltage calculations provides a pathway to understand and design concentration gradient materials.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Shock compression of strongly correlated oxides: A liquid-regime equation of state for cerium(IV) oxide

Here, the shock Hugoniot for full-density and porous CeO 2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ = 2.5 to 20 g/cm 3 . The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO 2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO 2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO 2 . In addition, a revised multiphase SESAME-type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Measuring the shock impedance mismatch between high-density carbon and deuterium at the National Ignition Facility

Fine-grained diamond, or high-density carbon (HDC), is being used as an ablator for inertial confinement fusion (ICF) research at the National Ignition Facility (NIF). Accurate equation of state (EOS) knowledge over a wide range of phase space is critical in the design and analysis of integrated ICF experiments. Here in this paper, we report shock and release measurements of the shock impedance mismatch between HDC and liquid deuterium conducted during shock-timing experiments having a first shock in the ablator ranging between 8 and 14 Mbar. Using ultrafast Doppler imaging velocimetry to track the leading shock front, we characterize the shock velocity discontinuity upon the arrival of the shock at the HDC/liquid deuterium interface. Comparing the experimental data with tabular EOS models used to simulate integrated ICF experiments indicates the need for an improved multiphase EOS model for HDC in order to achieve a significant increase in neutron yield in indirect-driven ICF implosions with HDC ablators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Elliptic flow coalescence to identify the $f_0$(980) content

Here we use a simple coalescence model to generate $f_0$(980) particles for three configurations: a $s\bar{s}$ meson, a $u\bar{u}s\bar{s}$ tetraquark, and a $K^+K^–$ molecule. The phase-space information of the coalescing constituents is taken from a multiphase transport (AMPT) simulation of heavy-ion collisions. It is shown that the number of constituent quarks scaling of the elliptic flow anisotropy can be used to discern $s\bar{s}$ from $u\bar{u}s\bar{s}$ and $K^+K^–$ configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗