Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK WAVES”

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

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2021

Filamentary surface plasma discharge flow length and time scales

Nanosecond surface dielectric barrier discharges (ns-SDBDs) are a class of plasma actuators that utilize a high-voltage pulse of nanosecond duration between two surface-mounted electrodes to create an electrical breakdown of air, along with rapid heating. These actuators usually produce multiple filaments when operated at high pulse frequencies, and the rapid heating leads to the formation of shock waves and complex flow fields. In this work we replicate a single filament of the ns-SDBDs and characterize the induced flow using velocity measurements from particle image velocimetry and density measurements from background-oriented schlieren. The discharge is produced by a high voltage electrical pulse between two copper electrodes on an acrylic base. A hot gas kernel characterizes the flow field formed close to the electrodes that expands and cools over time and a vortex ring that propagates away from the surface while entraining cold ambient fluid. The gas density deficit inside the kernel displays a power-law decay over time. Based on the observations, we develop a simplified theoretical model based on vortex-driven cooling and perform a scaling analysis to obtain the induced flow length and time scales. The results show that the cooling process's time scales correspond to a circulation-based time scale of the vortex ring, and the length scale of the kernel corresponds to the vortex ring radius. Furthermore, these findings can guide the choice of optimal filament spacing and pulse frequencies in the design, deployment, and operation of ns-SDBDs for flow control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of the density field around a supersonic conical projectile using quantitative schlieren

Quantitative schlieren imaging is a flow measurement technique that is capable of measuring density fields throughout refractive flowfields. The technique was applied here to measure the density field surrounding supersonic conical projectiles in free flight. Shock waves attached to a supersonic conical projectile offer a simple geometry with the well-established Taylor–Maccoll analytical flow solution to which these experimental measurements were compared. The schlieren images recorded a projection of the index of refraction field surrounding the 10° half-angle cones which was converted to density first through an Abel inversion and then the Gladstone–Dale law. Three Abel inversion methods—two-point, three-point, and arbitrary ray axisymmetric projection (ARAP)—were applied to deconvolute the three-dimensional flow within the constrained axisymmetric flow field. The resulting reconstructed density profiles were compared to the Taylor–Maccoll solution, parameterized by cone geometry and Mach number. The experimental density fields demonstrated strong agreement with the theoretical profiles. Experimental consistency was confirmed across various projectile speeds, demonstrating quantitative schlieren’s capability to accurately reconstruct the density of the flow field, even within the resolution constraints imposed by high-speed imaging. In conclusion, an assessment of experimental uncertainties in the density reconstruction was performed.

Abel Inversion↗
DOE OSTI2021

The effects of laser pulse length and collisional ionization on the acceleration of titanium ions

Abstract The interaction of a relativistic laser pulse ( > 10 18 W cm −2 µ m 2 ) with foil targets can accelerate ions to energies of tens of MeV u −1 with optimized laser and target parameters. We report the results on simulations of the interaction of a 6.0 × 10 20 W cm −2 laser pulse incident on ultrathin (10–500 nm) titanium foils to investigate the roles of laser pulse duration and ionization mechanisms in the acceleration of titanium ions. While holding peak intensity constant, two laser pulse durations were investigated, 140 and 650 fs. The optimum thickness is dependent on pulse duration, as it requires the concurrence of target transparency with the incidence of the peak laser intensity. The collisional processes do not play a significant role in Ti ion beam generation from the 140 fs laser pulse duration at the optimum thickness (30 nm). However, for the 650 fs laser optimum, collisions improve the conversion efficiency of highly energetic ( > 10 MeV u −1 ), high charge titanium (Ti 21 − 22 + ) by a factor of 20, and the titanium ion cutoff energy by ∼15%. This improvement is due to the fact that collisional ionization increases the electron density of the plasma, which delays the time of relativistic transparency, causing collisions to decrease the optimum foil thickness from 150 to 100 nm. Additionally, collisional ionization increases the charge-to-mass ratio of the titanium, and injects more electrons into the accelerating sheath field. At the optimum thickness, target normal sheath acceleration is the dominant mechanism of acceleration, with additional contributions from radiation pressure and shock wave acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Staged Z-pinch experiment at the double-EAGLE pulsed power facility

The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.

Physics - Plasma physics↗
DOE OSTI2020

Quantum Monte Carlo determination of the principal Hugoniot of deuterium

Here, we present coupled electron-ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock-wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the potential energy surface by wave-function optimization within variational Monte Carlo and projection quantum Monte Carlo methods. Compared to a previous study, our calculations also include low pressure-temperature (P,T) conditions resulting in close agreement with experimental data, while our revised results at higher (P,T) conditions still predict a more compressible Hugoniot than experimentally observed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Comparing theoretical predictions of radiation-free velocities of edge dislocations to molecular dynamics simulations

Transonic defect motion is of interest for high strain-rate plastic deformation as well as for crack propagation. Ever since Eshelby's 1949 prediction [J. D. Eshelby, Proc. Phys. Soc. A 62, 307 (1949)] in the isotropic limit of a “radiation-free” transonic velocity $v$ RF = $\sqrt{2}$c T , where shock waves are absent, there has been speculation about the significance of radiation-free velocities (if they truly exist) for defect mobility. Here, in this work, we argue that they do not play any significant role in dislocation dynamics in metals, based on comparing theoretical predictions of radiation-free velocities for transonic edge dislocations with molecular dynamics simulations for two face-centered cubic metals: Ag, where theory predicts radiation-free states, and Cu, where it does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Time dependent signatures of core-collapse supernova neutrinos at HALO

In this work, we calculate the response of a lead-based detector, such as the Helium and Lead Observatory (HALO) or its planned upgrade HALO-1kt to a Galactic core-collapse supernova. We pay particular attention to the time dependence of the reaction rates. All reaction rates decrease as the neutrino luminosity exponentially drops during the cooling period, but the ratio of one-neutron (1n) to two-neutron (2n) event rates in HALO is independent of this overall decrease. Nevertheless, we find that this ratio still changes with time due to the changing character of neutrino flavor transformations with the evolving conditions in the supernova. In the case of inverted hierarchy, this is caused by the fact that the spectral splits become less and less sharp with the decreasing luminosity. In the case of normal hierarchy, it is caused by the passage of the shock wave through the Mikheyev-Smirnov-Wolfenstein resonance region. However, in both cases, we find that the change in the ratio of 1n to 2n event rates is limited to a few percent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Surface stress tensor and junction conditions on a rotating null horizon

The general form of the surface stress tensor of an infinitesimally thin shell located on a rotating null horizon is derived, when different interior and exterior geometries are joined there. Even though the induced metric on the surface must be the same approached from either side, the first derivatives of the metric need not be. Such discontinuities lead to a Dirac δ-distribution in the Einstein tensor localized on the horizon. For a general stationary axisymmetric geometry the surface stress tensor can be expressed in terms of two geometric invariants that characterize the surface, namely the discontinuities [κ] and [J] of the surface gravity κ and angular momentum density J. The Komar energy and angular momentum are given in coordinates adapted to the Killing symmetries, and the surface contributions to each determined in terms of [κ] and [J] . Guided by these, a simple modification of the original Israel junction conditions is verified directly from the Einstein tensor density to give the correct finite result for the surface stress, when the normal n to the surface is allowed to tend continuously to a null vector. The relation to Israel’s original junction conditions, which fail on null surfaces, is given. The modified junction conditions are suitable to the matching of a rotating “black hole” exterior to any interior geometry joined at the Kerr null horizon surface, even when the surface normal is itself discontinuous and the Barrabès-Israel formalism is also inapplicable. This joining on a rotating null horizon is purely of the matter shell type and does not contain a propagating gravitational shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Correspondence between Color Glass Condensate and High-Twist Formalism

The color glass condensate (CGC) effective theory and the collinear factorization at high twist (HT) are two well-known frameworks describing perturbative QCD multiple scatterings in nuclear media. It has long been recognized that these two formalisms have their own domain of validity in different kinematic regions. Taking direct photon production in proton-nucleus collisions as an example, we clarify for the first time the relation between CGC and HT at the level of a physical observable. We show that the CGC formalism beyond shock-wave approximation, and with the Landau-Pomeranchuk-Migdal interference effect is consistent with the HT formalism in the transition region where they overlap. Such a unified picture paves the way for mapping out the phase diagram of parton density in nuclear medium from dilute to dense region.

Parton distribution functions↗
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 OSTI2024

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2020

A computational study of the thermodynamic conditions leading to autoignition in nanosecond pulsed discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results towards resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model. This can afford the use of very detailed chemical kinetic information. The full non-equilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

ignition, plasma, kinetics↗
DOE OSTI2025

Characterization of ablator dynamics initiated by picket-pulse conditions

Direct-drive laser-fusion targets may include low-cost plastic ablator material. Here, this work describes a tabletop experimental system designed to enable the investigation of the interaction of such ablator material as a function of the laser wavelength and the associated excitation mechanism. This experimental system offers time-resolved imaging of the material response with subpicosecond temporal resolution and ≈1-μm spatial resolution. The system involves synchronized femtosecond and picosecond lasers to initiate and characterize the dynamics of polystyrene planar targets irradiated by conditions similar to a “picket” prepulse from a direct-drive pulse shape. Ultrafast imaging of the initial excitation process followed by the resulting critical plasma formation, shock-wave propagation, and blowoff plasma expansion is used to compare the impact of excitation wavelength, 355 nm or 266 nm, on the ensuing dynamics. This work aims to investigate fundamental laser material interactions to help advance physics models that can inform the development of laser driver and target designs.

42 ENGINEERING↗
DOE OSTI2023

Ultrafast x-ray detection of low-spin iron in molten silicate under deep planetary interior conditions

The spin state of Fe can alter the key physical properties of silicate melts, affecting the early differentiation and the dynamic stability of the melts in the deep rocky planets. The low-spin state of Fe can increase the affinity of Fe for the melt over the solid phases and the electrical conductivity of melt at high pressures. However, the spin state of Fe has never been measured in dense silicate melts due to experimental challenges. We report detection of dominantly low-spin Fe in dynamically compressed olivine melt at 150 to 256 gigapascals and 3000 to 6000 kelvin using laser-driven shock wave compression combined with femtosecond x-ray diffraction and x-ray emission spectroscopy using an x-ray free electron laser. The observation of dominantly low-spin Fe supports gravitationally stable melt in the deep mantle and generation of a dynamo from the silicate melt portion of rocky planets.

58 GEOSCIENCES↗
DOE OSTI2021

Stiff and strong, lightweight bi-material sandwich plate-lattices with enhanced energy absorption

Plate-based lattices are predicted to reach theoretical Hashin–Shtrikman and Suquet upper bounds on stiffness and strength. However, simultaneously attaining high energy absorption in these plate-lattices still remains elusive, which is critical for many structural applications such as shock wave absorber and protective devices. In this article, we present bi-material isotropic cubic + octet sandwich plate-lattices composed of carbon fiber-reinforced polymer (stiff) skins and elastomeric (soft) core. This bi-material configuration enhances their energy absorption capability while retaining stretching-dominated behavior. We investigate their mechanical properties through an analytical model and finite element simulations. Our results show that they achieve enhanced energy absorption approximately 2–2.8 times higher than their homogeneous counterparts while marginally compromising their stiffness and strength. When compared to previously reported materials, these materials achieve superior strength-energy absorption characteristics, making them an excellent candidate for stiff and strong, lightweight energy absorbing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Empirical Acoustic Source Model for Chemical Explosions in Air

Chemical explosions generate pressure disturbances in air that radiate as nonlinear shock waves near the source and transition into acoustic waves with distance. Because low-frequency acoustic waves generally travel large distances without significant loss of energy, they are often used for explosion monitoring and yield estimation. However, quantitative relationships between acoustic energy and explosion yields are required for accurate yield estimation. Here, we develop an empirical acoustic source model for chemical explosions from experimental data. The empirical model returns the acoustic pressure waveform for the detonation of 1 kg of trinitrotoluene, which is conventionally used to represent the explosive release of 4.184 MJ of explosion energy. The full-waveform model can be used to predict acoustic signals for an arbitrary yield of a high-explosive detonation based on the standard scaling law and to estimate acoustic energies in a specific frequency range. We evaluate the accuracy of the acoustic source model independently by estimating the yield of other explosive events that are not included in the model development. Finally, statistical characteristics of the model and their implications for the uncertainty quantification of estimated yields are discussed.

58 GEOSCIENCES↗