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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bridging material models across scales: An integrated approach to equation of state and molecular dynamics modeling of copper

New uncertainty-aware equation of state (EOS) and electrical conductivity (EC) models for copper have been developed. The multiphase EOS/EC models are fit to experimental solid/liquid EC isobar measurements as well as density-functional theory molecular dynamics (DFT-MD) EC calculations in both expanded and compressed regimes (0.1–16 g/ cm 3 ⁠). The liquid and solid EOS phases were fit to available experimental data along with additional DFT-MD data over the same range as the EC. Leveraging the DFT-MD data, a corresponding machine-learned interatomic potential (MLIAP) for copper was trained using genetic-algorithm optimization. The copper MLIAP was constrained by EOS shock points at high compressions. The final EOS bounded MLIAP proves to be stable over a large density range (approximately 0.1–20 g/ cm 3 ) with good agreement to an isothermal compression curve, shock Hugoniot, and liquid speed of sound measurements at high pressures (100s of GPa).

Acoustic measurements and instrumentation↗
DOE OSTI2026

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗
DOE OSTI2020

The role of turbulent fluctuations in aerosol activation and cloud formation

Aerosol indirect effects are one of the leading contributors to cloud radiative properties relevant to climate. Aerosol particles become cloud droplets when the ambient relative humidity (saturation ratio) exceeds a critical value, which depends on the particle size and chemical composition. In the traditional formulation of this problem, only average, uniform saturation ratios are considered. Using experiments and theory, we examine the effects of fluctuations, produced by turbulence. Our measurements, from a multiphase, turbulent cloud chamber, show a clear transition from a regime in which the mean saturation ratio dominates to one in which the fluctuations determine cloud properties. The laboratory measurements demonstrate cloud formation in mean-subsaturated conditions (i.e., relative humidity <100%) in the fluctuation-dominant activation regime. The theoretical framework developed to interpret these measurements predicts a transition from a mean- to a fluctuation-dominated regime, based on the relative values of the mean and standard deviation of the environmental saturation ratio and the critical saturation ratio at which aerosol particles activate or become droplets. The theory is similar to the concept of stochastic condensation and can be used in the context of the atmosphere to explore the conditions under which droplet activation is driven by fluctuations as opposed to mean supersaturation. Overall, it provides a basis for future development of cloud droplet activation parameterizations that go beyond the internally homogeneous parcel calculations that have been used in the past.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Metastable–solid phase diagrams derived from polymorphic solidification kinetics

Nonequilibrium processes during solidification can lead to kinetic stabilization of metastable crystal phases. A general framework for predicting the solidification conditions that lead to metastable-phase growth is developed and applied to a model face-centered cubic (fcc) metal that undergoes phase transitions to the body-centered cubic (bcc) as well as the hexagonal close-packed phases at high temperatures and pressures. Large-scale molecular dynamics simulations of ultrarapid freezing show that bcc nucleates and grows well outside of the region of its thermodynamic stability. An extensive study of crystal–liquid equilibria confirms that at any given pressure, there is a multitude of metastable solid phases that can coexist with the liquid phase. We define for every crystal phase, a solid cluster in liquid (SCL) basin, which contains all solid clusters of that phase coexisting with the liquid. Additionally, a rigorous methodology is developed that allows for practical calculations of nucleation rates into arbitrary SCL basins from the undercooled melt. It is demonstrated that at large undercoolings, phase selections made during the nucleation stage can be undone by kinetic instabilities amid the growth stage. On these bases, a solidification–kinetic phase diagram is drawn for the model fcc system that delimits the conditions for macroscopic grains of metastable bcc phase to grow from the melt. We conclude with a study of unconventional interfacial kinetics at special interfaces, which can bring about heterogeneous multiphase crystal growth. A first-order interfacial phase transformation accompanied by a growth-mode transition is examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nonequilibrium statistical thermodynamics of multicomponent interfaces

Nonequilibrium interfacial thermodynamics has important implications for crucial biological, physical, and industrial-scale transport processes. Here, we discuss a theory of local equilibrium for multiphase multicomponent interfaces that builds upon the “sharp” interface concept first introduced by Gibbs, allowing for a description of nonequilibrium interfacial processes such as those arising in evaporation, condensation, adsorption, etc. By requiring that the thermodynamics be insensitive to the precise location of the dividing surface, one can identify conditions for local equilibrium and develop methods for measuring the values of intensive variables at the interface. We then use extensive, high-precision nonequilibrium molecular dynamics (NEMD) simulations to verify the theory and establish the validity of the local equilibrium hypothesis. In particular, we demonstrate that equilibrium equations of state are also valid out of equilibrium, and can be used to determine interfacial temperature and chemical potential(s) that are consistent with nonequilibrium generalizations of the Clapeyron and Gibbs adsorption equations. We also show, for example, that, far from equilibrium, temperature or chemical potential differences need not be uniform across an interface and may instead exhibit pronounced discontinuities. However, even in these circumstances, we demonstrate that the local equilibrium hypothesis and its implications remain valid. These results provide a thermodynamic foundation and computational tools for studying or revisiting a wide variety of interfacial transport phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Biphasic decay of intact SHIV genomes following initiation of antiretroviral therapy complicates analysis of interventions targeting the reservoir

The latent reservoir for HIV-1 in resting CD4 + T cells persists despite antiretroviral therapy (ART) and precludes cure. Reservoir-targeting interventions are evaluated in ART-treated macaques infected with simian immunodeficiency virus (SIV) or simian-human immunodeficiency virus (SHIV). Efficacy is determined by reservoir measurements before and after the intervention. However, most proviruses persisting in the setting of ART are defective. In addition, intact HIV-1 and SIV genomes undergo complex, multiphasic decay observable when new infection events are blocked by ART. Intervention-induced elimination of latently infected cells must be distinguished from natural decay. Here, we address these issues for SHIV. We describe an intact proviral DNA assay that allows digital counting of SHIV genomes lacking common fatal defects. We show that intact SHIV genomes in circulating CD4 + T cells undergo biphasic decay during the first year of ART, with a rapid first phase (t 1/2 = 30.1 d) and a slower second phase (t 1/2 = 8.1 mo) that is still more rapid that the slow decay observed in people with HIV-1 on long-term ART (t 1/2 = 3.7 y). In SHIV models, most interventions are tested during 2nd phase decay. Natural 2nd phase decay must be considered in evaluating interventions as most infected cells present at this time do not become part of the stable reservoir. In addition, for interventions tested during 2nd phase decay, a caveat is that the intervention may not be equally effective in people with HIV on long-term ART whose reservoirs are dominated by latently infected cells with a slower decay rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The kinetics of SARS-CoV-2 infection based on a human challenge study

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Model Development for Thermal-Hydrology Simulations of a Full-Scale Heater Experiment in Opalinus Clay

Disposal of commercial spent nuclear fuel in a geologic repository is studied. In situ heater experiments in underground research laboratories provide a realistic representation of subsurface behavior under disposal conditions. Here, this study describes process model development and modeling analysis for a full-scale heater experiment in opalinus clay host rock. The results of thermal-hydrology simulation, solving coupled nonisothermal multiphase flow, and comparison with experimental data are presented. The modeling results closely match the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Corrosion rates measured under potentiostatic conditions to represent performance of a cermet nuclear waste form

Electrochemical tests were conducted to assess the corrosion performance of a developmental multiphase cermet nuclear waste form in different conditions. Corrosion currents were monitored as the surface equilibrated in an acidic brine solution during potentiostatic tests at different imposed potentials. Electrochemical impedance spectroscopy was used to quantify differences in the electrical properties of surfaces stabilised at different potentials. The passivating oxides that will stabilise the surface under disposal conditions are formed and maintained during testing. The measured corrosion currents were two orders of magnitude lower for passivated surfaces than for polished surfaces over a wide range of redox conditions. The electrochemical results were correlated with microscopic analyses of the corroded surfaces and physical models of the surface passivation and degradation processes are proposed. The measured corrosion rates are suitable for the long-term performance modeling of cermet nuclear waste forms.

Gattu, V. K.↗
DOE OSTI2020

Controlled plasma–droplet interactions: a quantitative study of OH transfer in plasma–liquid interaction

Plasmas in contact with liquids are a rich source of OH radicals and have been extensively studied in the last decade to leverage the ability to generate chemically reactive species in gas phase plasmas to decompose organics. Multiphase transfer of OH radicals is highly transport limited and to overcome transport limits, the plasma activation of aerosols, small liquid droplets, interspersed in the plasma has been proposed. In this work, we report a combined experimental and modeling study of a controlled plasma–droplet interaction experiment using a diffuse RF glow discharge in He + 0.2% H 2 O with detailed plasma diagnostics, ex situ analysis of the plasma-induced chemistry in the droplet containing formate, droplet trajectory and size measurements. This enables a quantitative study of the reactivity transfer of OH from the gas phase plasma to the liquid phase and how its diffusion limitations impact formate decomposition in the water droplet. For a droplet with a diameter of 36 μm, we observed 50% reduction in formate concentration in the droplets after plasma treatment for droplet residence times in the plasma of ~10 ms. These short droplet residence times in the plasma allow in some cases for droplet size reductions of ~5% in spite gas temperatures of 360 K. A one-dimensional reaction–diffusion model was used to calculate the OH transport and formate oxidation inside the droplet and was able to predict the conversion of formate by plasma in a droplet without any fitting parameters. The model further shows that formate conversion is dominated by near-interfacial reactions with OH radicals and is limited by diffusion of formate in the droplet. Here, the results show that a controlled plasma–micro-droplet reactor as reported in this study might be an excellent tool for detailed quantitative plasma–liquid interaction studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Droplet bioprinting of acellular and cell-laden structures at high-resolutions

Advances in digital light projection(DLP) based (bio) printers have made printing of intricate structures at high resolution possible using a wide range of photosensitive bioinks. A typical setup of a DLP bioprinter includes a vat or reservoir filled with liquid bioink, which presents challenges in terms of cost associated with bioink synthesis, high waste, and gravity-induced cell settling, contaminations, or variation in bioink viscosity during the printing process. Here, we report a vat-free, low-volume, waste-free droplet bioprinting method capable of rapidly printing 3D soft structures at high resolution using model bioinks and model cells. A multiphase many-body dissipative particle dynamics model was developed to simulate the dynamic process of droplet-based DLP printing and elucidate the roles of surface wettability and bioink viscosity. Process variables such as light intensity, photo-initiator concentration, and bioink formulations were optimized to print 3D soft structures (∼0.4–3 kPa) with a typical layer thickness of 50 µm, an XY resolution of 38 ± 1.5 μm and Z resolution of 237 ± 5.4 µm. To demonstrate its versatility, droplet bioprinting was used to print a range of acellular 3D structures such as a lattice cube, a Mayan pyramid, a heart-shaped structure, and a microfluidic chip with endothelialized channels. Droplet bioprinting, performed using model C3H/10T1/2 cells, exhibited high viability (90%) and cell spreading. Additionally, microfluidic devices with internal channel networks lined with endothelial cells showed robust monolayer formation while osteoblast-laden constructs showed mineral deposition upon osteogenic induction. Overall, droplet bioprinting could be a low-cost, no-waste, easy-to-use, method to make customized bioprinted constructs for a range of biomedical applications.

DLP↗
DOE OSTI2020

The age–chemical abundance structure of the Galactic disc – II. α -dichotomy and thick disc formation

ABSTRACT We extend our previous work on the age–chemical abundance structure of the Galactic outer disc to the inner disc (4 < r < 8 kpc) based on the SDSS/APOGEE survey. Different from the outer disc, the inner disc stars exhibit a clear bimodal distribution in the [Mg/Fe]–[Fe/H] plane. While a number of scenarios have been proposed in the literature, it remains challenging to recover this bimodal distribution with theoretical models. To this end, we present a chemical evolution model embedding a complex multiphase inner disc formation scenario that matches the observed bimodal [Mg/Fe]–[Fe/H] distribution. In this scenario, the formation of the inner disc is dominated by two main starburst episodes $6\,$Gyr apart with secular, low-level star formation activity in between. In our model, the first starburst occurs at early cosmic times ($t\sim 1\,$ Gyr) and the second one $6\,$ Gyr later at a cosmic time of $t\sim 7\,$ Gyr. Both these starburst episodes are associated with gas accretion events in our model, and are quenched rapidly. The first starburst leads to the formation of the high-α sequence, and the second starburst leads to the formation of the metal-poor low-α sequence. The metal-rich low-α stars, instead, form during the secular evolution phase between the two bursts. Our model shows that the α-dichotomy originates from the rapid suppression of star formation after the first starburst. The two starburst episodes are likely to be responsible for the formation of the geometric thick disc (z >1 kpc), with the old inner thick disc and the young outer thick disc forming during the first and the second starbursts, respectively.

Lian, Jianhui↗
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 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↗