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

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2024

Microstructure effects on high velocity microparticle impacts of copper

Constitutive models can fail to predict high-rate deformation behavior due to their inability to account for microstructural effects. In part, this is because of a dearth of experimental benchmarking data in the high strain-rate, low pressure regime, since many high-rate experiments also probe a region of strong shockwaves, at which point microstructure effects no longer play a primary role. This work uses laser-induced particle impact testing to quantitatively study high velocity impacts of small, rigid alumina microspheres on flat copper substrates with varying amounts of initial cold work in the weak shock regime, but at very high strain rates up to ~10 7 s –1 . Through paired experiments and numerical simulations, this work shows that the initial microstructure condition can have significant influence on dynamical mechanical properties in this range. Specifically, prior work hardening of the copper substrate leads to increased rebounding of the microparticles (i.e., less plastic dissipation in the impact) as well as smaller craters. Each of these experimental measurables can be converted into a strength measure, i.e., the dynamic yield strength or dynamic hardness, respectively, neither of which is well predicted consistently by existing constitutive laws. The general trend of hardening can be captured by such models by incorporating an existing “pre-strain,” suggesting that future calibration of the materials parameters may yield a good fit over a broader range of conditions. Our results emphasize the importance of reporting the microstructural condition in dynamic studies, as well as the necessity of accounting for these factors when formulating and optimizing constitutive models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Megahertz-rate shock-wave distortion cancellation via phase conjugate digital in-line holography

Abstract Holography is a powerful tool for three-dimensional imaging. However, in explosive, supersonic, hypersonic, cavitating, or ionizing environments, shock-waves and density gradients impart phase distortions that obscure objects in the field-of-view. Capturing time-resolved information in these environments also requires ultra-high-speed acquisition. To reduce phase distortions and increase imaging rates, we introduce an ultra-high-speed phase conjugate digital in-line holography (PCDIH) technique. In this concept, a coherent beam passes through the shock-wave distortion, reflects off a phase conjugate mirror, and propagates back through the shock-wave, thereby minimizing imaging distortions from phase delays. By implementing the method using a pulse-burst laser setup at up to 5 million-frames-per-second, time-resolved holograms of ultra-fast events are now possible. This technique is applied for holographic imaging through laser-spark plasma-generated shock-waves and to enable three-dimensional tracking of explosively generated hypersonic fragments. Simulations further advance our understanding of physical processes and experiments demonstrate ultra-high-speed PCDIH techniques for capturing dynamics.

Science & Technology - Other Topics↗
DOE OSTI2025

Integrated simulations of premagnetized and self-magnetizing dynamic screw pinch-driven MagLIF

Magnetically driven implosions such as in magnetized liner inertial fusion (MagLIF) on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that dynamically redistribute liner mass during implosion, limiting fusion fuel compression and confinement, which ultimately degrades performance. Driving the implosion with an initially helical drive field that dynamically shifts the direction of the magnetic field surrounding the liner (i.e., a dynamic screw pinch, DSP) is a method proposed to mitigate MRTI in-flight and improve target performance. In DSPs, the axial drive magnetic field component implodes the liner and diffuses through the shocked, melted liner material into the fuel throughout the implosion. Liners can be designed to enable enough axial magnetic flux to diffuse through the liner material to effectively magnetize the fuel region without the need of an initial axial magnetic field (i.e., from external field coils). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing drive magnetic fields composed of axial and azimuthal components (a helical drive field). These simulated DSP-driven MagLIF targets demonstrate improved fuel conditions and thermonuclear yield compared to a traditional MagLIF target implosion. Synthetic x-ray radiography of the imploding liner material and x-ray emission images of the fuel region at the time of peak neutron yield rate indicate superior implosion morphology for DSP-MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Assessment of Physics Models for Phase Transition Kinetics

The time dependence of phase diagrams and how to model rate dependent transitions remains one of the key unanswered questions in physics. When a material is loaded dynamically through equilibrium phase boundaries, it is the kinetics that determines the real time expression of a phase transition. Here we report the atomic and nanosecond-scale quantification of kinetics of shock-driven phase transition in multiple materials. We uniquely make use of a both a simple shock as well as shock-and-hold loading pathways compress different crystalline solids and induce structural phase transitions below melt. Coupling shock loading with time-resolved synchrotron x-ray diffraction (DXRD), we probe the structural transformations of these solids in the short-lived high pressure and temperature states generated. The novelty and power of using DXRD for the assessment of kinetics of phase transitions lies in the ability to discover and identify new phases and to examine kinetics without prior knowledge of a material's phase diagram. Our results provide a quantified expression and a physics model of kinetics of formation of high-pressure phases under shock loading: transition incubation time, evolution, completion time and crystallization rate.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock compression experiments using the DiPOLE 100-X laser on the high energy density instrument at the European x-ray free electron laser: Quantitative structural analysis of liquid Sn

X-ray free electron laser (XFEL) sources coupled to high-power laser systems offer an avenue to study the structural dynamics of materials at extreme pressures and temperatures. The recent commissioning of the DiPOLE 100-X laser on the high energy density (HED) instrument at the European XFEL represents the state-of-the-art in combining x-ray diffraction with laser compression, allowing for compressed materials to be probed in unprecedented detail. Here, we report quantitative structural measurements of molten Sn compressed to 85(5) GPa and ~3500 K. The capabilities of the HED instrument enable liquid density measurements with an uncertainty of ~1% at conditions which are extremely challenging to reach via static compression methods. We discuss best practices for conducting liquid diffraction dynamic compression experiments and the necessary intensity corrections which allow for accurate quantitative analysis. We also provide a polyimide ablation pressure vs input laser energy for the DiPOLE 100-X drive laser which will serve future users of the HED instrument.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Thermal and mechanical influences on shear band formation and suppression in shocked 1,3,5-trinitroperhydro-1,3,5-triazine (RDX)

High-pressure shear band formation is a critical phenomenon in energetic materials because of its ability to form hotspots and influence mechanical strength. Shear banding is known to occur in a variety of these materials, but the governing dynamics of the mechanisms are not well defined for molecular crystals. Our previous work has found that at high pressures in 1,3,5-trinitroperhydro-1,3,5-triazine (RDX), the initial formation sites for shear bands, called “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth, suppressing the shear banding nucleation and growth. Here, in this work, we assess the influence of a variety of changes to the material state on this phenomenon, including altered initial temperature, lateral strain that confines the system in tension or pressure, and initial molecular vacancies throughout the crystal. Shear band suppression and the nature of the shear band network are assessed as a function of each of these.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydrodynamic instability at impact interfaces and planetary implications

Abstract Impact-induced mixing between bolide and target is fundamental to the geochemical evolution of a growing planet, yet aside from local mixing due to jetting – associated with large angles of incidence between impacting surfaces – mixing during planetary impacts is poorly understood. Here we describe a dynamic instability of the surface between impacting materials, showing that a region of mixing grows between two media having even minimal initial topography. This additional cause of impact-induced mixing is related to Richtmyer-Meshkov instability (RMI), and results from pressure perturbations amplified by shock-wave refraction through the corrugated interface between impactor and target. However, unlike RMI, this new impact-induced instability appears even if the bodies are made of the same material. Hydrocode simulations illustrate the growth of this mixing zone for planetary impacts, and predict results suitable for experimental validation in the laboratory. This form of impact mixing may be relevant to the formation of stony-iron and other meteorites.

58 GEOSCIENCES↗
DOE OSTI2025

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗
DOE OSTI2023

In-Situ X-Ray Imaging High Strain Rate Compression of Laminate Al-Graphene Composite and Mechanical Property Characterization

Owing to the opaque nature of the laminated structures, traditional high-speed optical camera cannot be used to detect the dynamic process of sub-surface deformation. In this article, we report a study of using high speed X-ray imaging to study the high strain rate deformation in laminated Al structures. We used a Kolsky bar apparatus to apply dynamic compression and a high-speed synchrotron X-ray phase contrast imaging (PCI) setup to conduct the in situ X-ray imaging study. The in situ X-ray imaging captures the shock wave propagation in the laminated structures. After shock compression, we characterized the microstructures by using transmission electron microscopy (TEM), which demonstrates an increase of dislocation density. In conclusion, the micro-pillar compression tests show that the yield strength at 0.2% offset of laminated Al-graphene composite has a significant increase of 67%, from 30 to 50 MPa, compared to laminate Al after shock loading.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Equation-of-state measured via x-ray phase contrast imaging for Epon 828/DEA epoxy

Epoxies are a broad class of polymer materials often used as adhesive, structural or binding materials. Epon 828 is an epoxy resin that can be polymerized with a variety of curing agents with the choice of curing agent potentially having an effect on the resulting epoxy polymer’s material properties. In this study, the dynamic behavior of Epon 828 epoxy resin cured with diethanolamine (DEA) is investigated through a series of tamped Richtmyer-Meshkov instability (RMI) experiments measured with x-ray phase-contrast imaging. The measured shock and particle velocities are combined with data in the literature to calibrate Mie-Grüneisen equations-of-state (EOS) for portions and combinations of the collective dataset. The calibrated Mie-Grüneisen EOS are validated against particle velocity profiles extracted from published literature using the Eulerian hydrocode CTH. Here, the Mie-Grüneisen EOS fit to only the tamped RMI experimental data presented here most closely follows the particle velocity profile in the published literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi-frame, ultrafast, x-ray microscope for imaging shockwave dynamics

Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A cavitation and dynamic void growth model for a general class of strain-softening amorphous materials

In this work, we derive a homogenized framework for cavitation and dynamic void growth in a general class of strain-softening materials with particular emphasis on amorphous materials. Analytic solutions are derived by considering an assemblage of spherical shells subject to dynamic hydrostatic tensile pressure. The isochoric response of the matrix material is captured by a free volume theory based amorphous viscoplasticity model. The framework accounts for the coupled rate-dependent, strain-softening behavior typically exhibited in amorphous materials. Validation of the homogenized theory is carried out against numerous direct numerical simulations with excellent agreement. Furthermore, we find that a reduction in cavitation strength and significant acceleration of dynamic void growth rates is induced as a consequence of strain-softening in the amorphous matrix material. Lastly, the model is utilized to interpret some non-intuitive experimental observations of amorphous materials subject to shock compression and subsequent spall failure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On the local measurement of electric currents and magnetic fields using Thomson scattering in Weibel-unstable plasmas

Here, we demonstrate the capability of the Thomson Scattering (TS) diagnostic to measure locally the microscopic electron and ion currents in counter-streaming plasmas unstable to the Weibel or current-filamentation instability. Synthetic TS spectra are calculated with particle distribution functions obtained from particle-in-cell simulations and used to accurately reproduce the simulated currents. We show that this technique allows accurate local measurements of the magnetic field, thus opening the way for the complete experimental characterization of the growth rate, saturation, and nonlinear dynamics of electromagnetic instabilities in plasmas. We illustrate the application of this diagnostic to experimental TS data, which yields local measurements of the magnetic field in Weibel-unstable plasmas and indicates that the magnetic energy density reaches ~1% of the kinetic energy density of the flows, in agreement with previous numerical studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The temperatures of ejecta transporting in vacuum and gases

In this work, we measure continuous thermal radiance from evolving clouds of liquid metal fragments ejected into vacuum, nonreactive, and reactive gas. We implement a model for the thermalization of the ejecta and gas and use this to constrain the absolute temperature of the ejecta cloud. This model enables further analyses of ejecta thermal behavior under a variety of conditions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗