SEARCH · Engineering Papers

Results for “HUGONIOT EQUATION OF STATE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Hypervelocity impact heating of porous aluminum

Estimates of the thermal energy and other properties of the postshock condition are obtained for porous aluminum targets impacted by iron or aluminum projectiles, in the shock pressure range between about 0.4 and 8 Mbar. The starting point is the determination of a distinct Hugoniot equation for each value of porosity, from available experimental data. Rankine-Hugoniot equations and a Mie-Gruneisen equation of state supply the relations necessary for finding the thermodynamic properties of adiabatically relaxed materials. The results are of interest in such fields as ablation studies, meteoritics, and lunar topography.

Rehfuss, D. E.↗

DOE OSTI2020

Shock Hugoniot of forged and additively manufactured 304L stainless steel

The purpose of this research was to measure the equation of state for additively manufactured (AM) and forged 304L stainless steel. An understanding of the dynamic behavior of AM metals is integral to their timely adoption into industrial and governmental applications. The Hugoniot of the AM 304L was compared to that of the forged 304L. This comparison enabled us to determine the sensitivity of the equation of state to processing parameters and the resulting microstructure. Our results show that there is a measurable difference in the Hugoniot between the AM 304L and forged 304L. The shock wave velocities for the AM 304L were found to be slower than those for the forged 304L at similar particle speeds perhaps because of the presence of porosity in the AM 304L. While our measurements of density showed no distinct trend toward a measurable decrease in density in the AM 304L, there was a fair amount of scatter in these results.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Hugoniot properties of porous stainless steel: Insights from molecular dynamics simulations

Here, we present a systematic study of Hugoniot properties of porous 316L stainless steel using both a simple interpolation scheme and direct shock simulations in order to probe pore collapse kinetics as well as final thermodynamic states. Both methods indicate that equilibrated Hugoniot properties depend on pore density only and not on the pore distribution or size. We then create a simple porous equation of state model that is shown to be accurate for a range of validation data. This allows us to extend our simulations to make direct comparison to experimental data that have generally significantly larger system sizes and durations. In addition, our direct shock simulations indicate that the relaxation time after hotspot formation is system size dependent and can reach nanosecond timescales for the largest pores investigated in our study, thereby possibly having a measurable effect on fast dynamic loading experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Shock Hugoniot of Forged and Additively Manufactured 304L Stainless Steel

The purpose of this research was to measure the equation of state for additively manufactured (AM) and forged 304L stainless steel using a novel experimental technique. An understanding of the dynamic behavior of AM metals is integral to their timely adoption into various applications. The Hugoniot of the AM 304L was compared to that of the forged 304L at particle velocities where the material retains a two-wave structure. This comparison enabled us to determine the sensitivity of the equation of state to microstructure as varied due to processing. Our results showed that there was a measurable difference in the measured shock velocity between the AM and forged 304L. The shock wave velocities for the AM 304L were found to be ~3% slower than those for the forged 304L at similar particle velocities. To understand these differences, properties such as densities, sound speeds, and texture were measured and compared between the forged and AM materials. Our results showed that no measurable difference was found in these properties. Additionally, it is possible that differing elastic wave amplitudes may influence shock velocity.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Shock Hugoniot of forged and additively manufactured 304L stainless steel

36 MATERIALS SCIENCE↗

DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Understanding the effects of radiative preheat and self-emission from shock heating on equation of state measurement at 100s of Mbar using spherically converging shock waves in a NIF hohlraum

Over the last six years many experiments have been done at the National Ignition Facility to measure the Hugoniot of materials, such as CH plastic at extreme pressures, up to 800 Mbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. Here, the shock front conditions are characterized using x-ray radiography. In this paper we examine the role of radiation in heating the unshocked material in front of the shock to understand the impact it has on equation of state measurements and how it drives the measured data off the theoretical Hugoniot curve. In particular, the two main sources of radiation heating are the preheating of the unshocked material by the high-energy kilo-electron-volt x-rays in the hohlraum and the heating of the material in front of the shock, as the shocked material becomes hot enough to radiate significantly. Using our model, we estimate that preheating can reach 4 eV in unshocked material, and that radiation heating can begin to drive data off the Hugoniot significantly, as pressures reach above 400 Mb.

Equations of state↗

NASA NTRS1971

Equation of state for porous metals under strong shock compression.

An equation of state is developed for metals in the high-pressure high-temperature region that can be extended in a thermodynamically consistent manner to the vapor phase. The constants are obtained from observed thermodynamic properties and are specified for Al, Fe, Cu, Ni, Pb, and W. This equation of state was checked against experimental data by computing the Hugoniot curves for normal and porous samples for cases where experimental data are available. Also, the velocity of sound and release isentropes were computed for various shock conditions.

Naumann, R. J.↗

DOE OSTI2020

Nonideal mixing effects in warm dense matter studied with first-principles computer simulations

Here, we study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO 3 . First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2 × 10 5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, which emerge because of the ionization of L and K shell electrons, are well represented, and the maximum compression ratio of the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly important role. However, the results imply that the equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Shock initiation of the HMX-based explosive PBX 9012: Experiments, uncertainty analysis, and unreacted equation-of-state

In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.

36 MATERIALS SCIENCE↗

DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Equation of state of hot, dense magnesium derived with first-principles computer simulations

Using two first-principles computer simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium in the regime of warm dense matter, with densities ranging from 0.43 to 86.11 g cm –3 and temperatures from 20,000 K to 5×10 8 K. These conditions are relevant for the interiors of giant planets and stars as well as for shock compression measurements and inertial confinement fusion experiments. Here, we study ionization mechanisms and the electronic structure of magnesium as a function of density and temperature. We show that the L shell electrons, 2s and 2p energy bands, merge at high densities. This results in gradual ionization of the L-shell with increasing density and temperature. In this regard, Mg differs from MgO, which is also reflected in the shape of its principal shock Hugoniot curve. For Mg, we predict a single broad pressure-temperature region, where the shock compression ratio is approximately 4.9. Mg thus differs from Si and Al plasmas that exhibit two well-separated compression maxima on the Hugoniot curve for L and K shell ionizations. Finally, we study multiple shocks and effects of preheat and precompression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗