Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

CO Laser Absorption Measurements During Syngas Combustion at High Pressure

Abstract Syngas is a desirable fuel for combustion in the Allam-Fetvedt cycle, which involves combustion under supercritical-CO2 conditions. While some work has been conducted in collecting ignition delay times (IDT) at the extreme pressures required by these systems, significant model deficiencies remain. Additionally, considerable barriers in terms of nonideal gas dynamic effects have been shown for these experiments in shock tubes. Further investigation into the fundamental combustion kinetics of H2/CO/CO2 mixtures is required. Time-resolved speciation measurements for target species have been shown to better aid in improving the understanding of underlying chemical kinetics than global ignition delay time measurements. Therefore, laser absorption measurements of CO were measured behind reflected shock waves during combustion of syngas at 5 and 10 bar and temperatures between 1080 and 2100 K. The mixtures investigated utilized H2-to-CO ratios of 1:1 and 1:4, respectively, each at stoichiometric conditions, allowing for discussions of the effect of initial fuel composition. A ratio of fuel to CO2 of 1:2 was also utilized to represent commercially available syngas. The mixtures were diluted in helium and argon (20% He, 76.5% Ar) to minimize thermal effects and to expedite CO thermal relaxation during the experiment. The resulting CO time histories were then compared to modern chemical kinetics mechanisms, and disagreement is seen for this system, which is assumed to be fairly well known. This study elucidates particular chemistry that needs improvement in moving toward a better understanding of syngas combustion at elevated pressures.

Engineering↗
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)↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2022

Limited Dynamic Triggering in the Utah Region, USA

The state of Utah, USA, experiences around 3800 catalogued earthquakes per year, highlighting that the region is seismically active and susceptible to earthquakes. Following the 2002 Denali Fault (M7.9) earthquake in Alaska, the region showed an elevated seismicity rate for 3 weeks following the passage of high amplitude surface waves, suggesting that the region may be particularly susceptible to dynamic triggering. With over 23,396 faults and each fault presenting a potential fault for triggering, we systematically search for dynamic triggering throughout the state of Utah caused by large, global earthquakes with M ≥ 7. Specifically, we analyse earthquake catalogues and all available waveform data to determine statistical increases of seismicity rate following the passage of seismic arrivals. While we find instances of dynamic triggering, our results show that these events occur sparsely in the region. In total, less than 20 per cent of the 273 main shocks that occur from 2000 to the end of 2017 show a statistical indication of dynamic triggering throughout the Utah region, highlighting that dynamic triggering is limited for stresses created by transient signals from global M ≥ 7 earthquakes, with the exception being the Denali Fault (M7.9), Alaska earthquake (i.e. an instance of significant triggering).

58 GEOSCIENCES↗
DOE OSTI2021

In Situ X-ray Diffraction of High Density Polyethylene During Dynamic Compression: Polymer Chain Compression and Phase Transformations

Polymers are widely implemented materials often used for safety and security in our daily lives. However, we do not have a thorough understanding of the phase transformations of these materials when exposed to extreme condition, high strain rate, events. To study polymers in these severe environments, we performed plate impact experiments that impart a one-dimensional shock wave, at high pressures and temperatures, to observe the materials response. Generally, velocimetry is the experimental diagnostic used to capture the movement of high velocity interfaces to then calculate the mechanical variables associated with shock compression. However, when a material undergoes a phase transition, such as between two crystalline phases or melting to decomposition during shock compression, velocimetry does not explicitly indicate what structural change occurred. At the Dynamic Compression Sector at the Advanced Photon Source gas gun driven compression has been paired with x-ray diffraction, providing the ability to study both the shock parameters and structural changes concurrently. We explored high density polyethylene (HDPE) both on its reactants Hugoniot (<10.12 GPa) and within its transition region (24.3 GPa) where HDPE is hypothesized to decompose from reactants to products. With x-ray diffraction we observed both polymer chain compression and crystal structure reorientation from orthorhombic to monoclinic above 6 GPa, as well as either a possible phase change or chemical decomposition at 24.3 GPa. With these experiments we were able to determine the crystalline and amorphous Hugoniot loci at pressures below the polymer transition region.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Internal Digital Image Correlation Technique for High-Strain Rate Dynamic Experiments

Full-field, quantitative visualization techniques, such as digital image correlation (DIC), have unlocked vast opportunities for experimental mechanics. However, DIC has traditionally been a surface measurement technique, and has not been extended to perform measurements on the interior of specimens for dynamic, full-scale laboratory experiments. This limitation restricts the scope of physics which can be investigated through DIC measurements, especially in the context of heterogeneous materials. The focus of this study is to develop a method for performing internal DIC measurements in dynamic experiments. The aim is to demonstrate its feasibility and accuracy across a range of stresses (up to 650 MPa), strain rates (10 3 - 10 6 s -1 ), and high-strain rate loading conditions (e.g., ramped and shock wave loading). Internal DIC is developed based on the concept of applying a speckle pattern at an inner-plane of a transparent specimen. The high-speed imaging configuration is coupled to the traditional dynamic experimental setups, and is focused on the internal speckle pattern. During the experiment, while the sample deforms dynamically, in-plane, two-dimensional deformations are measured via correlation of the internal speckle pattern. In this study, the viability and accuracy of the internal DIC technique is demonstrated for split-Hopkinson (Kolsky) pressure bar (SHPB) and plate impact experiments. The internal DIC experimental technique is successfully demonstrated in both the SHPB and plate impact experiments. In the SHPB setting, the accuracy of the technique is excellent throughout the deformation regime, with measurement noise of approximately 0.2% strain. In the case of plate impact experiments, the technique performs well, with error and measurement noise of 1% strain. The internal DIC technique has been developed and demonstrated to work well for full-scale dynamic high-strain rate and shock laboratory experiments, and the accuracy is quantified. Here, the technique can aid in investigating the physics and mechanics of the dynamic behavior of materials, including local deformation fields around dynamically loaded material heterogeneities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluation of a novel mid-infrared pyrometer for dynamic material experiments

Temperature is a complicated thermodynamic parameter to measure in dynamic compression experiments. Optical pyrometry is a general-purpose “work-horse” technique for measuring temperature from a radiant surface on these experimental platforms. The optical pyrometry channels are commonly held to the visible or Near-Infrared spectrum, which provides high fidelity temperature measurement for shock temperature above ∼1200–1500 K. However, low temperature (T < 1200 K) dynamic material experiments, including low pressure or quasi-isentropic studies, as well as experiments with complex thermodynamic paths, require Mid-Infrared (Mid-IR) for high fidelity measurements. This article outlines the design, testing, and characterization of a novel Mid-IR pyrometer system that can be configured between 2.5 and 5.0 µm, suitable for lower temperature measurements and for increasing the fidelity and precision of higher temperature measurements. Experimental validation was done on two separate gas gun platforms, with two separate impact velocities, achieving temperatures between 450 and 1100 K.

Dutra, E. C. (ORCID:0000000167649195)↗
DOE OSTI2026

ECLEIRS: Exact conservation law embedded identification of reduced states for parameterized nonlinear conservation laws from sparse and noisy data

Multi-query applications such as parameter estimation, uncertainty quantification and design optimization for parameterized partial differential equation (PDE) systems are expensive. While reduced/latent state dynamics approaches for parameterized PDEs offer a viable alternative, these approaches rely on high-quality data and struggle with highly sparse spatiotemporal noisy measurements typically obtained from experiments. Furthermore, there is no guarantee that these models satisfy governing physical conservation laws. In this article, we propose a reduced state dynamics approach, referred to as ECLEIRS, that embeds exact conservation in the solution and flux representation by utilizing a space-time divergence-free neural network formulation. We compare ECLEIRS with other reduced state dynamics approaches, those that do not enforce any physical constraints and those with physics-informed loss functions, for three shock-propagation problems: 1-D advection, 1-D Burgers and 2-D Euler equations. In conclusion, the numerical experiments conducted in this study demonstrate that ECLEIRS provides the most accurate prediction of dynamics for unseen parameters even in the presence of highly sparse and noisy data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Compensating slice emittance growth in high brightness photoinjectors using sacrificial charge

Achieving maximum electron beam brightness in photoinjectors requires detailed control of the 3D bunch shape and precise tuning of the beam focusing. Even in state-of-the-art designs, slice emittance growth due to nonlinear space charge forces and partial nonlaminarity often remains non-negligible. In this work, we introduce a new means to linearize the transverse slice phase space: a sacrificial portion of the bunch’s own charge distribution, formed into a wavebroken shock front by highly nonlinear space charge forces within the gun, whose downstream purpose is to dynamically linearize the desired bunch core. We show that linearization of an appropriately prepared bunch can be achieved via strongly nonlaminar focusing of the sacrificial shock front, while the inner core focuses laminarly. This leads to a natural spatial separation of the two distributions: a dense core surrounded by a diffuse halo of sacrificial charge that can be collimated. Multiobjective genetic algorithm optimizations of the ultracompact x-ray free electron laser injector employ this concept, and we interpret it with an analytic model that agrees well with the simulations. In simulation, we demonstrate a final bunch charge of 100 pC, peak current ∼ 30 A, and a sacrificial charge of 150 pC (250 pC total emitted from cathode) with normalized emittance growth of < 20 nm rad due to space charge. This implies a maximum achievable brightness approximately an order of magnitude greater than existing free electron laser injector designs. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗
DOE OSTI2022

Recovery of release cloud from laser shock-loaded graphite and hydrocarbon targets: in search of diamonds

Abstract This work presents first insights into the dynamics of free-surface release clouds from dynamically compressed polystyrene and pyrolytic graphite at pressures up to 200 GPa, where they transform into diamond or lonsdaleite, respectively. These ejecta clouds are released into either vacuum or various types of catcher systems, and are monitored with high-speed recordings (frame rates up to 10 MHz). Molecular dynamics simulations are used to give insights to the rate of diamond preservation throughout the free expansion and the catcher impact process, highlighting the challenges of diamond retrieval. Raman spectroscopy data show graphitic signatures on a catcher plate confirming that the shock-compressed PS is transformed. First electron microscopy analyses of solid catcher plates yield an outstanding number of different spherical-like objects in the size range between ten(s) up to hundreds of nanometres, which are one type of two potential diamond candidates identified. The origin of some objects can unambiguously be assigned, while the history of others remains speculative.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Shock compression of diamond single crystals to 120 GPa: Refractive index and nonlinear photoelasticity

The optical response of transparent solids at extreme conditions is important for both fundamental science and many applications. Strong transparent solids are of particular interest for use as optical windows in dynamic compression experiments. Due to diamond’s exceptional strength and optical properties, laser-driven shock experiments and plate impact experiments were carried out to examine the diamond optical response for shock wave compression along two different crystal orientations. Using laser interferometry at 532 nm and 1550 nm wavelengths, optical transparency was observed and refractive indices were determined for [100] diamond at stresses up to 119 GPa and for [111] diamond at stresses up to 87 GPa. From these results, the nonlinear photoelastic response for [100] and [111] diamond was determined, revealing significant dependence on both crystal orientation and laser wavelength. To enable [100] diamond as an interferometry window in dynamic compression experiments, the requisite window corrections were determined for 532 nm and 1550 nm wavelengths. Because of diamond’s excellent x-ray transparency, the present findings will be particularly useful for incorporating [100] diamonds as windows in dynamic compression experiments involving x-ray diffraction or other x-ray diagnostics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
DOE OSTI2026

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗