Search NASA⌕ Search

SEARCH · Search NASA

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
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 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 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↗
DOE OSTI2021

Model investigation of the longitudinal broadening of the transverse momentum two-particle correlator

Here, the multiphase transport model is used to investigate the longitudinal broadening of the transverse momentum two-particle correlator C 2 (Δη,Δφ), and its utility to extract the specific shear viscosity, η/s, of the quark-gluon plasma formed in ultrarelativistic heavy ion collisions. The results from these model studies indicate that the longitudinal broadening of C 2 (Δη,Δφ) is sensitive to the value of η/s. However, reliable extraction of the longitudinal broadening of the correlator requires the suppression of possible self-correlations associated with the definition of the collision centrality.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Elliptic flow splittings in the Polyakov–Nambu–Jona-Lasinio transport model

To incorporate the effect of gluons on the evolution dynamics of the quark matter produced in relativistic heavy-ion collisions, we extend the three-flavor Nambu–Jona-Lasinio (NJL) transport model to include the contribution from the Polyakov loops. Imbedding the resulting pNJL partonic transport model in an extended multiphase transport (extended AMPT) model, we then study the elliptic flow splittings between particles and their antiparticles in relativistic heavy-ion collisions at beam energy scan energies. We find that a weak quark vector interaction in the partonic phase is able to describe the elliptic flow splitting between protons and antiprotons in heavy-ion collisions at √sNN=7.7 to 39 GeV. Knowledge of the quark vector interaction is useful for understanding the equation of state of quark matter at large baryon chemical potentials and thus the location of the critical point in the QCD phase diagram.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Utilization of event shape in search of the chiral magnetic effect in heavy-ion collisions

The search for the chiral magnetic effect (CME) has been a subject of great interest in the field of high-energy heavy-ion collision physics, and various observables have been proposed to probe the CME. Experimental observables are often contaminated with background contributions arising from collective motions (specifically elliptic flow) of the collision system. We present a method study of event-shape engineering (ESE) that projects the CME-sensitive γ 112 correlator and its variations (γ 132 and γ 123 ) to a class of events with minimal flow. We discuss the realization of the zero-flow mode, the sensitivity on the CME signal, and the corresponding statistical significance for Au + Au, Ru + Ru, and Zr + Zr collisions at √s NN = 200 GeV with a multiphase transport model, as well as a new event generator, event-by-event anomalous-viscous fluid dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗