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

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Shock compression of poly(methyl methacrylate) PMMA in the 1000 GPa regime: Z machine experiments

Hydrocarbon polymers are used in a wide variety of practical applications. In the field of dynamic compression at extreme pressures, these polymers are used at several high energy density (HED) experimental facilities. One of the most common polymers is poly(methyl methacrylate) or PMMA, also called Plexiglass® or Lucite®. Here, we present high-fidelity, hundreds of GPa range experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for PMMA to more than threefold compression up to 650 GPa and re-shock Hugoniot states up to 1020 GPa in an off-Hugoniot regime, where experimental data are even sparser. These data can be used to put additional constraints on tabular equation of state (EOS) models. The present results provide clear evidence for the need to re-examine the existing tabular EOS models for PMMA above ~120 GPa as well as perhaps revisit EOSs of similar hydrocarbon polymers commonly used in HED experiments investigating dynamic compression, hydrodynamics, or inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Behavior of soda-lime silicate glass under laser-driven shock compression up to 315 GPa

Shock experiments give a unique insight into the behavior of matter subjected to extremely high pressures and temperatures. Understanding the behavior of materials under such extreme conditions is key to modeling material failure and deformation dynamics under impact. While studies on pure silica are extensive, the shock behavior of other commercial silicates that contain additional oxides has not been systematically investigated. To better understand the role of composition in the dynamic behavior of silicates, we performed laser-driven dynamic compression experiments on soda-lime glass (SLG) up to 315 GPa. Using the accurate pulse shaping offered by the long pulse laser system at the Matter in Extreme Conditions end-station at the Linac Coherent Light Source, SLG was shock compressed along the Hugoniot to multiple pressure-temperature points. Velocity Interferometer System for Any Reflector was used to measure the velocity and determine the pressure inside the SLG. The U s –u p relationship obtained agrees well with the previous parallel plate impact studies. Within the error bars, no transformation to the crystalline phase was observed up to 70 GPa, which is in contrast to the behavior of pure silica under shock compression. Our studies show that the glass composition strongly influences the shock compression behavior of the silicate glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proton acceleration in an overdense hydrogen plasma by intense CO 2 laser pulses with nonlinear propagation effects in the underdense pre-plasma

We report on proton acceleration from intense CO 2 laser-irradiated hydrogen plasmas at near-critical densities, with the density gradient steepened by Nd:YAG laser ablation-driven hydrodynamic shocks. While the experimental results, such as the quasi-monoenergetic proton spectra and their scaling with respect to the laser energy, are generally in agreement with the simulations, certain laser shots produced significantly higher proton energies than anticipated during the experiment. The increased proton energy may be linked to nonlinear propagation effects in the steepened plasma density ramp before the critical surface, including relativistic self-focusing and, for the case of temporally-structured laser pulses observed in the experiment, focusing of the trailing pulse through the plasma channel formed by the leading pulse 25 ps ahead. The occurrence of channel focusing in the underdense hydrogen plasma is supported by a subsequent pump-probe experiment with a dark-field imaging technique, where the formation of ion channels was observed after the passage of an intense CO 2 laser pulse

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An assessment of spall failure modes in laser powder bed fusion fabricated stainless steel 316L with low-volume intentional porosity

Here, this paper reports on spall failure and damage modes in Laser Powder Bed Fusion fabricated Stainless Steel 316L (SS316L) with intentional levels of low-volume (1–5 vol. %) porosity and pore sizes of 200, 350, and 500 μm. The fabricated specimens were subjected to uniaxial-strain plate-impact loading at ~4.5 GPa, to initiate incipient spall failure. Analysis of velocimetry profiles measured using multi-probe photon-Doppler velocimetry coupled with post-mortem analysis of soft-recovered samples reveals local suppression of spall failure (termed as spall-dominated) as a function of porosity, as the failure mechanism transitions from spall-centered tensile stress dominated to a pore-centered microstructure-dominated damage mode involving void/crack nucleation and growth at pre-existing pores. The critical porosity level where the suppression of spall failure is first observed, as well as the spall location, is dependent on both the volume fraction and the size of the initially fabricated pores. In samples of 500 μm pore size, the suppression of spall failure is observed with as little as 1 vol. % porosity, while samples with smaller pores (200 μm) still experience spall-centered tensile stress dominated failure with higher levels (5 vol. %) of porosity. In the case of pore-centered microstructure-dominated failure, spall damage can occur but the spall plane is shifted toward the rear free surface, or more generally in areas further away from the region with pores. Highly heterogeneous deformation twinning, shear banding, grain rotation, and cracking are observed in the vicinity of pre-existing pores and expected spall failure sites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the importance of three-dimensional modeling for high-energy-density physics experiments

Laser-driven cylindrical implosion experiments enable direct measurements of hydrodynamic instability growth in convergent geometries, providing a wealth of validation data in the high-energy-density regime. These experiments are designed to be nearly axially invariant, allowing for modeling with complementary two-dimensional slices of the cylinder. Two distinct hydrodynamics codes are employed to model a subset of these experiments, and the results are shown to be in very good agreement with each other and the available experimental data. While this 2D modeling approach adequately captures most of the physics of the implosion and ensuing instability growth, there are crucial aspects from the three-dimensional nature of the experiments that are missed in 2D. The first fully 3D simulations of these experiments are presented, and small but significant differences are found to arise from both the axial and azimuthal non-uniformity in the laser drive. Recent experimental results confirming the drive asymmetry are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The colliding planar shocks platform to study warm dense matter at the National Ignition Facility

We have developed an experimental platform at the National Ignition Facility that employs colliding planar shocks to produce warm dense matter with uniform conditions and enable high-precision equation of state measurements. The platform uses simultaneous x-ray Thomson scattering and x-ray radiography to measure the density, electron temperature, and ionization state in warm dense matter. The experimental platform is designed to create a large volume of uniform plasma (approximately 700×700×150μm 3 ) at pressures approaching 100 Mbar and minimize the distribution of plasma conditions in the x-ray scattering volume, significantly improving the precision of the measurements. Here, in this study, we present the experimental design of the platform and compare hydrodynamic simulations to x-ray radiography data from initial experiments studying hydrocarbons, producing uniform densities within ±25% of the average probed condition. We show that the platform creates a homogeneous plasma that can be characterized using x-ray Thomson scattering. Thus, the new platform enables accurate measurements of plasma conditions necessary to test models for the equation of state and ionization potential depression in the warm dense matter regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

High precision control of laser energy for laser–matter interaction studies

Precise, highly reproducible control of the laser energy is required for high confidence laser–matter interaction research such as in dynamic compression science and high energy density physics. The energy must be adjustable without affecting the pulse shape (time varying intensity) or beam smoothness. We have developed a convenient two-stage energy tuning method for a nominal 100 J, 351 nm (UV) laser. The energy is adjusted in 10 J (10%) increments by operating the laser at full energy and inserting a beam splitter in the laser output. As the splitter is located after the final frequency tripling optics, the UV pulse shape is unchanged. The energy is varied by substituting a splitter of different reflectivity. For finer 3 J (3%) increments, the infrared pulse is attenuated inside the laser before the final amplifier. This requires modest tuning to preserve the pulse shape. The demonstrated variation in shot-to-shot reproducibility is less than +/–2.5 J (5% of the full energy), irrespective of the laser output energy. These approaches can be adapted to most ~100 J class lasers. We describe these techniques and show two examples where they have elucidated the underlying physics in laser shock compression experiments. One used only the beam splitters to establish the pressure for melting in iron. The other combined both techniques to finely increment the peak stress (~2 GPa steps) in germanium to precisely determine the onset and completion of melting—including the melting kinetics. Furthermore, these unambiguous results would not be possible without the developments described here.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hugoniot properties of porous stainless steel: Insights from molecular dynamics simulations

Here, we present a systematic study of Hugoniot properties of porous 316L stainless steel using both a simple interpolation scheme and direct shock simulations in order to probe pore collapse kinetics as well as final thermodynamic states. Both methods indicate that equilibrated Hugoniot properties depend on pore density only and not on the pore distribution or size. We then create a simple porous equation of state model that is shown to be accurate for a range of validation data. This allows us to extend our simulations to make direct comparison to experimental data that have generally significantly larger system sizes and durations. In addition, our direct shock simulations indicate that the relaxation time after hotspot formation is system size dependent and can reach nanosecond timescales for the largest pores investigated in our study, thereby possibly having a measurable effect on fast dynamic loading experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Laser material interactions in tamped materials on picosecond time scales in aluminum

Here, a 100 ps laser is used to probe the pressure generation, depth of the non-solid ablator, and the non-linear optical effects through tamper materials. Samples consisted of an aluminum ablator with tampers of sapphire and coverslip glass. In general, the sapphire tamped sample achieves higher pressures at lower laser intensities as compared to the coverslip glass tamped sample. Attempts to model the details of this set of experimental data with standard available radiation coupled hydrodynamic codes make clear that more physics is needed in these simulations to accurately predict the impact of the tamper material on the pressure generation and the depth of non-solid aluminum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗