Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hugoniot”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Shock effects in olivine and implications for Hugoniot data

New observations of shock deformation in single-crystal olivine are presented for the range in peak pressures from about 16 to 75 GPa. A nonequilibrium model is suggested for the behavior of olivine under shock which explains the densification associated with a mixed-phase region and which is consistent with all available observations, yet which bears little similarity to the equilibrium behavior of olivine at high pressure.

Jeanloz, R.↗
NASA NTRS1985

Fast and optimal solution to the Rankine-Hugoniot problem

A new, definitive, reliable and fast iterative method is described for determining the geometrical properties of a shock (i.e., theta sub Bn, yields N, V sub s and M sub A), the conservation constants and the self-consistent asymptotic magnetofluid variables, that uses the three dimensional magnetic field and plasma observations. The method is well conditioned and reliable at all theta sub Bn angles regardless of the shock strength or geometry. Explicit proof of uniqueness of the shock geometry solution by either analytical or graphical methods is given. The method is applied to synthetic and real shocks, including a bow shock event and the results are then compared with those determined by preaveraging methods and other iterative schemes. A complete analysis of the confidence region and error bounds of the solution is also presented.

Vinas, A. F.↗
NASA NTRS1986

Fast and optimal solution to the 'Rankine-Hugoniot problem'

A new, definitive, reliable and fast iterative method is described for determining the geometrical properties of a shock (i.e., theta sub Bn, yields N, V sub s and M sub A), the conservation constants and the self-consistent asymptotic magnetofluid variables, that uses the three dimensional magnetic field and plasma observations. The method is well conditioned and reliable at all theta sub Bn angles regardless of the shock strength or geometry. Explicit proof of uniqueness of the shock geometry solution by either analytical or graphical methods is given. The method is applied to synthetic and real shocks, including a bow shock event and the results are then compared with those determined by preaveraging methods and other iterative schemes. A complete analysis of the confidence region and error bounds of the solution is also presented.

Vinas, A. F.↗
NASA NTRS1991

Effect of void-size distribution on the Hugoniot state at low shock pressures

In most theoretical and experimental investigations into the shock response of underdense solid media, the influence of the medium's mesostructure on the resulting pressure and degree of compaction has not been taken into account. In typical cases examined, shock pressures are well in excess of 1 GPa and this approach is clearly justified. However, at low pressures, calculations show that the distribution of void sizes can affect the final state achieved upon shocking the medium from a given initial porosity. This paper analyzes the response of porous aluminum to low pressure shocking and demonstrates a dependence of the final shocked state on the distribution of void sizes.

Griffiths, David J.↗
NASA NTRS1994

An improved solution to the 'Rankine-Hugoniot' problem

This paper presents an extension of the nonlinear least squares fitting technique of Vinas and Scudder (1986) (VS), which finds the physical and geometrical properties of nondissipational magnetohydrodynamic (MHD) shocks. The new method incorporates plasma temperature observations in the form of normal momentum flux and energy density flux conservation as well as plasma density, velocity, and magnetic field data. The new technique is capable of using known standard deviations in the individual measurement points to properly weight the fitting procedure. The new fitting code is validated through the analysis of synthetic shocks with known physical and geometrical properties. Finally, it is compared to the original VS method and the preaveraged velocity coplanarity technique.

Szabo, A.↗
DOE OSTI2024

Sound speed measurements in shock compressed cemented tungsten carbide: Evolution of elastic moduli with damage at pressures to 100 GPa

The motivation of the present study is to gain insights into the evolution of elastic properties of cemented tungsten carbides (WC) shock compressed to 100 GPa. Seven plate impact experiments—two front surface impact and five release wave overtake—are conducted to make simultaneous measurements of Hugoniot states and longitudinal sound speeds in shocked WC with 3.7wt.% cobalt binder. The sound speeds along with estimates for bulk sound speeds, obtained using the Birch–Murnaghan EoS, are analyzed to determine the elastic moduli—longitudinal, bulk, and shear—as a function of Hugoniot stress. The longitudinal and bulk sound speeds at Hugoniot states of interest are found to increase linearly with longitudinal stress. Consistent with the increase in sound speeds, the longitudinal and bulk moduli also increase with Hugoniot stress; however, the increase in longitudinal modulus is modest when compared to predictions of theoretical models that account for pressure and temperature dependence of elastic moduli, but with no damage. The shear moduli remain nearly constant at ~318 GPa over the range of Hugoniot states investigated. These values are, however, much lower than those predicted by the Steinberg–Guinan model with no damage. Poisson’s ratio decreases initially from its ambient value of 0.208 to ~0.199 for Hugoniot stress ≤10 GPa indicating consolidation of the WC microstructure with low initial stress; however, with an increase in Hugoniot stress to ~100 GPa, Poisson’s ratio increases to ~0.317, indicating degradation of shear moduli with increasing stress. The product of density and Grüneisen parameter (ρΓ), after an initial spike, remains nearly constant for volumetric strains ≥0.07. The maximum average temperature rise is estimated to be ~286°C at the highest Hugoniot stress employed in the study.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fluid-like shock compression of dilute polymer nanocomposites

The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.

Moore, Nathan W. [Sandia National Laboratories (SN↗
DOE OSTI2022

Vibrational response and temperature of shock-compressed Pt: In situ extended x-ray absorption fine structure measurements to 325 GPa

Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Peak states of molybdenum single crystals shock compressed to high stresses

In this work, to determine crystal anisotropy effects at high stresses, peak states behind the plastic shock waves were examined in BCC single crystals. Using plate impact experiments, molybdenum (Mo) single crystals were shock compressed up to 190 GPa elastic impact stress along [100], [110], and [111] orientations. Laser interferometry was used to measure wave velocities and particle velocity profiles at the Mo–LiF window interface. These data were analyzed to obtain in-material quantities in the peak states. The Hugoniots for [100] and [110] orientations were comparable, but the Hugoniot for the [111] orientation was different from the other two orientations. Also, these Mo single crystal Hugoniots display differences from the polycrystalline Mo Hugoniots. Although none of the differences can be considered large, the present results demonstrate that, unlike FCC metal single crystals (Cu, Al), some anisotropy is preserved in Mo single crystal Hugoniots even at high stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗
NASA NTRS1973

Shock compression of a recrystallized anorthositic rock from Apollo 15

Hugoniot measurements on 15,418, a recrystallized and brecciated gabbroic anorthosite, yield a value of the Hugoniot elastic limit (HEL) varying from 45 to 70 kbar as the final shock pressure is varied from 70 to 280 kbar. Above the HEL and to 150 kbar, the pressure-density Hugoniot is closely described by a hydrostatic equation of state constructed from ultrasonic data for single-crystal plagioclase and pyroxene. Above 150 kbar, the Hugoniot states indicate that a series of one or more shock-induced phase changes are occurring in the plagioclase and pyroxene. From Hugoniot data for both the single-crystal minerals and the Frederick diabase, we infer that the shock-induced high-pressure phases in 15,418 probably consists of a 3.71 g/cu cm density, high-pressure structure for plagioclase and a 4.70 g/cu cm perovskite-type structure for pyroxene.

Ahrens, T. J.↗