Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock-induced compression”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2025

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗
DOE OSTI2025

Exploring scenarios for enhanced fuel compression and performance on the National Ignition Facility with machine-learning-aided design techniques

Recent fusion experiments on the National Ignition Facility (NIF) have achieved ignition, producing multi-MJ fusion yields for input laser energies of roughly 2 MJ [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Building on the success of the target designs that have achieved ignition, we explore new implosion scenarios predicted to generate significantly more compression of the dense DT ice layer and correspondingly higher yields while preserving many of the key physics characteristics of present-day ignition designs. Our main result is a novel 3-shock implosion scheme that effectively minimizes the shock-induced entropy in the dense, accelerating DT shell and maximizes the resulting fuel compression subject to a fixed leading shock strength consistent with present-day ignition experiments, which is necessary to melt the crystalline high-density carbon ablator. Compared to the first NIF experiment to fulfill Lawson's ignition criterion, shot N210808 [Abu-Shawareb et al., Phys. Rev. Lett. 129, 075001 (2022)], our design exhibits a 40% increase in simulated peak areal density (ρR) and a 5× increase in 1D fusion yield using a 4% lighter ablator and identical DT payloads. We also present a complete integrated 2D hohlraum design and laser pulse specifications capable of generating the desired 3-shock drive and maintaining control of the low-mode capsule implosion symmetry, where the increase in simulated 2D yield relative to N210808 is > 10×. This new implosion regime was discovered with help from a machine-learning-enabled capsule design optimization framework. We outline the workflow this automated tool uses to identify improved design candidates by running several rounds of capsule simulations, constructing a surrogate model mapping input variations to key physics output quantities, and querying the resulting statistical model to propose adjustments to the x-ray drive and capsule to reach a set of physics objectives prescribed by the designer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

A shock-induced phase change in orthoclase.

New shock compression data to 340 kb for single-crystal orthoclase (along (001)), demonstrate the onset of a shock-induced phase change at about 115 kb. Along the Hugoniot a mixed-phase region extends to about 300 kb, above which the data are believed to correspond to the properties of a high-pressure phase having the hollandite structure (zero pressure density of 3.84 g/cu cm) reported by Ringwood et al. If the hollandite value for the zero pressure density is used, the zero pressure bulk modulus of this phase is approximately 1.8 plus or minus 0.2 Mb.

Ahrens, T. J.↗
DOE OSTI2026

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)↗
NASA NTRS1975

Control of shock-wave boundary-layer interactions by bleed in supersonic mixed compression inlets

An experimental investigation has been conducted to determine the effect of bleed region geometry and bleed rate on shock wave-boundary layer interactions in an axisymmetric, mixed-compression inlet at a Mach number of 2.5. The full realizable reduction in transformed form factor is obtained by bleeding off about half the incident boundary layer mass flow. Bleeding upstream or downstream of the shock-induced pressure rise is preferable to bleeding across the shock-induced pressure rise. Slanted holes are more effective than normal holes. Two different bleed hole sizes were tested without detectable difference in performance.

Fukuda, M. K.↗
DOE OSTI2025

Chemical decomposition in shock-compressed 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals: Time-resolved Raman spectroscopy

Understanding the evolution of chemical decomposition—at the molecular level—in shock-compressed insensitive high explosive (IHE) single crystals is an important need. Toward this goal, time-resolved Raman spectra were measured in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals shock compressed to stresses above 20 GPa. At 22 GPa and higher stresses, the Raman peak intensities were significantly reduced, showing the onset and extent of chemical changes. At 33 GPa, no Raman peaks were observed, suggesting complete decomposition of FOX-7. The 22 GPa onset from Raman data is in marked contrast to the much higher (32 GPa) onset determined previously from wave profile measurements—showing the significantly greater sensitivity of Raman measurements for gaining insight into shock-induced molecular-level changes in IHE single crystals. The close match between the complete loss of Raman peaks at 33 GPa and the 32 GPa reaction threshold determined from continuum data makes a good case that the two results are related and they reflect significant energy release due to decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS1974

Numerical simulation of high Reynolds number transonic flows

A code has been developed for simulating high Reynolds number transonic flow fields of arbitrary configuration. An explicit finite-difference method with time splitting is used to solve the time-dependent equations for compressible turbulent flow. A nonorthogonal computational mesh of arbitrary configuration facilitates the description of the flow field. The code is applied to simulate the flow over an 18 percent thick circular-arc biconvex airfoil at zero angle of attack and free-stream Mach number of 0.775. A simple mixing-length model is used to describe the turbulence and chord Reynolds numbers of 1, 2, 4, and 10 million are considered. The solution describes in sufficient detail both the shock-induced and trailing-edge separation regions, and provides the profile and friction drag.

Deiwert, G. S.↗
DOE OSTI2026

Shock-induced phase transitions and stacking fault formation in additively manufactured eutectic high-entropy alloy Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2

Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 additively manufactured using laser powder bed fusion (LPBF) is among the toughest as-built alloys reported and is a promising candidate for use in extreme environments. However, its behavior under multi-megabar pressure regimes remains unexplored. We used femtosecond in situ X-ray diffraction to investigate the shock response of LPBF Ni40Co20Fe10Cr10Al18W2 under laser-driven shock compression and release. Our results reveal that the initial dual-phase face-centered + body-centered structure transforms to a single face-centered phase over a wide pressure range of 84±10 to 277±55 GPa, followed by a transition to a single body-centered phase at 431±48 GPa. We establish the Hugoniot equation-of-state of LPBF Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 and compare it to the benchmark alloy AlCoCrFeNi 2.1 , demonstrating the effects of W-doping and increased Al content. High stacking fault probabilities, close to those measured in Au and Ag, are observed upon compression. A portion of the stacking faults are annihilated upon release to ambient pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Effect of a collapsing gas bubble on the shock-to-detonation transition in liquid nitromethane

We studied the shock-induced collapse of butane gas bubbles in the homogeneous explosive nitromethane (NM) to investigate the effects of hot spot formation on the detonation process. A butane bubble was injected into a sample of NM, and a shock wave from a flat plate impactor compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and bubble collapse processes. A multiband optical fiber pyrometer measured the time-resolved thermal radiance, and we used the results and emissivity values extracted from spectral fits to estimate temperatures. We measured the characteristics of the shock-to-detonation transition in NM with and without a bubble. All experiments were performed at shock pressures near 8 GPa, where neat NM can detonate. A single bubble in this system was shown to sensitize NM, leading to a reduced run-to-detonation time. We used hydrodynamic modeling to predict shock wave propagation, the extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. We used a temperature-dependent Arrhenius burn model for simulations, and it yielded much better results than reactive burn models that depend only on pressure and density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗
DOE OSTI2025

Nonseeded linewise temperature measurements by resonantly ionized photoemission thermometry in a Mach 4 Ludwieg tube

A one-dimensional (1D) thermometry using oxygen-tagging resonantly ionized photoelectron thermometry (O 2 RIPT) was employed to investigate thermal gradients within a Mach 4 Ludwieg tube. The Ludwieg tube is pulsed with a test duration of approximately 100 ms, providing a cold supersonic flow at Mach 4 ideal for studying aerothermal effects. This study focused on measuring freestream temperatures, capturing shock-induced heating behind a detached bow shock from a blunt cylinder, and resolving sharp temperature variations across a bow shock generated by a cylinder. The O 2 RIPT technique produced strong emission signals extending approximately 4 cm long, demonstrating its capability for precise temperature measurements in high-speed wind tunnel environments. The results confirm that O 2 RIPT is well-suited for applications in large-scale aerodynamic testing facilities, particularly in regions with strong compression effects, enabling the resolution of sharp thermal gradients. This method presents a promising solution for thermometry in dynamic flow conditions relevant to various experimental ground-test facilities.

McCord, Walker (ORCID:0000000179747550)↗
DOE OSTI2025

Numerical Treatment of Shock-induced Nuclear Burning in Double Detonation Type Ia Supernovae

We present a benchmark problem to assess the treatment of shock-induced nuclear burning in the context of double detonation Type Ia supernovae. In a stratified white dwarf model, we implement a shock-detection criterion that suppresses burning in zones characterized by compression and significant pressure gradients, controlled by a tunable parameter, f shock . One-dimensional simulations, using the open-source Castro suite, were conducted across three treatments—burning fully enabled, and burning suppressed with f shock = 2/3 and f shock = 1—across three spatial resolutions (5.0, 2.5, and 0.3125 km). At the finest resolution, the burning-enabled and f shock = 1 models converge, while the f shock = 2/3 front continues to show slight offset behavior. Since most simulations are carried out at much lower resolutions, our tests support the idea that burning in shocks should always be disabled in practice. We also observe that the behavior of lower-resolution simulations remains extremely sensitive to the choice of f shock .

GPU computing↗
NASA NTRS1972

Shock melting and vaporization of metals.

The effect of initial porosity on shock induction of melting and vaporization is investigated for Ba, Sr, Li, Fe, Al, U, and Th. For the less compressible of these metals, it is found that for a given strong shock-generation system (explosive in contact, or flyer-plate impact) an optimum initial specific volume exists such that the total entropy production, and hence the amount of metal liquid or vapor, is a maximum. Initial volumes from 1.4 to 2.0 times crystal volumes, depending on the metal sample and shock-inducing system, will result in optimum post-shock entropies.

Ahrens, T. J.↗
NASA NTRS1973

Prediction of space shuttle fluctuating pressure environments, including rocket plume effects

Preliminary estimates of space shuttle fluctuating pressure environments have been made based on prediction techniques developed by Wyle Laboratories. Particular emphasis has been given to the transonic speed regime during launch of a parallel-burn space shuttle configuration. A baseline configuration consisting of a lightweight orbiter and monolithic SRB, together with a typical flight trajectory, have been used as models for the predictions. Critical fluctuating pressure environments are predicted at transonic Mach numbers. Comparisons between predicted environments and wind tunnel test results, in general, showed good agreement. Predicted one-third octave band spectra for the above environments were generally one of three types: (1) attached turbulent boundary layer spectra (typically high frequencies); (2) homogeneous separated flow and shock-free interference flow spectra (typically intermediate frequencies); and (3) shock-oscillation and shock-induced interference flow spectra (typically low frequencies). Predictions of plume induced separated flow environments were made. Only the SRB plumes are important, with fluctuating levels comparable to compression-corner induced separated flow shock oscillation.

Plotkin, K. J.↗
NASA NTRS1973

Shock remagnetization associated with meteorite impact at planetary surfaces

The stages of a cratering event are considered together with four distinct magnetization volumes, taking into account the jetted material, the ejecta, the impact melt, and shock-compressed material which remains in situ. Shock classifications are discussed along with mineralogical aspects of shock remagnetization, the behavior of iron sulfides, defect cubic spinels, ilmenite-hematite series, the Fe-Ti-O system, alloys phases, and silicates. Remagnetization mechanisms considered include a first order reversible crystallographic transition in bcc iron-nickel, a second order Curie point transition in fcc iron-nickel, shock-induced anisotropy, and shock melting of iron containing silicates. Experimental results are reported together with a hypothetical crater model.

Wasilewski, P. J.↗