Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK WAVE INTERACTION”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparative analysis of internal energy excitation and dissociation of nitrogen predicted by independently developed ab initio potential energy surfaces

In this article we present a comparative atomic level study analyzing the vibrational excitation and dissociation of molecular nitrogen due to N 2 ⁢( 1 Σ g ⁢ + )+N⁢( 4 S u ) and N 2 ⁢( 1 Σ g ⁢ + )+N 2 ⁢( 1 Σ g ⁢+ ) interactions governed by independently developed potential energy surfaces at the University of Minnesota and NASA Ames Research Center. Further, vibrational excitation was studied for N 2 +N 2 interactions from T=10000 to 25000 K and for N 2 +N from T=5000 to 30000 K. Nonequilibrium dissociation is studied from T=10000 to 30000 K under the quasi-steady-state condition for N 2 +N 2 and N 2 +N interactions. Finally, an inviscid Mach 20 dissociating nitrogen flow over a cylinder with a Knudsen number of 0.015 is carried out to study the impact of molecular interactions predicted by independently developed potential energy surfaces on a canonical hypersonic flow.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Interactions of laser-driven tin ejecta microjets over phase transition boundaries

Ejecta microjets offer an experimental methodology to study high-speed particle laden-flow interactions, as microjets consist of millions of particulates traveling at velocities of several kilometers per second and are easily generated by most common shock drives. Previous experiments on the OMEGA Extended Performance laser found that collisions between two counter-propagating laser-driven tin ejecta microjets varied as a function of drive pressure; jets generated near shock pressures of 10 GPa passed through each other without interacting, whereas jets generated at shock pressures of over 100 GPa interacted strongly, forming a cloud around the center interaction point. In this paper, we present a more systematic scan of tin ejecta microjet collisions over intermediate pressure regimes to identify how and at what shock pressure interaction behavior onsets. Radiographs of interacting microjets at five different laser drive energies qualitatively demonstrate that interaction behavior onsets slowly as a function of laser drive energy. Quantitative mass and density metrics from each radiograph provide trends on jet characteristics and collisional mass dispersion. It is observed that jetting mass, jet densities, and mass dispersion increase with increasing drive pressures and that the increased jet density at the higher drive energies may account for the increased mass dispersion. This work provides an important step in the understanding of tin ejecta microjet collisions and paves the way for future studies on the physics dominating high-speed particle-laden flow interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fast imaging of laser ablation of cerium: Dependence of plasma plume dynamics and structure on native oxide growth

We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Laser Ablation Plasmas and Spectroscopy for Nuclear Applications

The development of measurement methodologies to detect and monitor nuclear-relevant materials remains a consistent and significant interest across the nuclear energy, nonproliferation, safeguards, and forensics communities. Optical spectroscopy of laser-produced plasmas is becoming an increasingly popular diagnostic technique to measure radiological and nuclear materials in the field without sample preparation, where current capabilities encompass the standoff, isotopically resolved and phase-identifiable (e.g., UO and UO[Formula: see text]) detection of elements across the periodic table. These methods rely on the process of laser ablation (LA), where a high-powered pulsed laser is used to excite a sample (solid, liquid, or gas) into a luminous microplasma that rapidly undergoes de-excitation through the emission of electromagnetic radiation, which serves as a spectroscopic fingerprint for that sample. This review focuses on LA plasmas and spectroscopy for nuclear applications, covering topics from the wide-area environmental sampling and atmospheric sensing of radionuclides to recent implementations of multivariate machine learning methods that work to enable the real-time analysis of spectrochemical measurements with an emphasis on fundamental research and development activities over the past two decades. Background on the physical breakdown mechanisms and interactions of matter with nanosecond and ultrafast laser pulses that lead to the generation of laser-produced microplasmas is provided, followed by a description of the transient spatiotemporal plasma conditions that control the behavior of spectroscopic signatures recorded by analytical methods in atomic and molecular spectroscopy. High-temperature chemical and thermodynamic processes governing reactive LA plasmas are also examined alongside investigations into the condensation pathways of the plasma, which are believed to serve as chemical surrogates for fallout particles formed in nuclear fireballs. Laser-supported absorption waves and laser-induced shockwaves that accompany LA plasmas are also discussed, which could provide insights into atmospheric ionization phenomena from strong shocks following nuclear detonations. Furthermore, the standoff detection of trace radioactive aerosols and fission gases is reviewed in the context of monitoring atmospheric radiation plumes and off-gas streams of molten salt reactors. Finally, concluding remarks will present future outlooks on the role of LA plasma spectroscopy in the nuclear community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Mechanics of pore array collapse and interaction in shock-compressed polymethyl methacrylate (PMMA)

Recent studies on dynamic pore collapse have revealed significant development of shear localization, which can lead to material failure in porous structures and hot spot generation in energetic materials. These findings have dramatically improved the understanding of failure mechanisms during pore collapse but also prompt further investigation of realistic porous materials. In particular, porous media consist of many pores and porous networks. Even in low-porosity materials, pores can form in close proximity during the manufacturing process, leading to the critical question of pore–pore interaction during collapse under dynamic loading conditions. This study investigates, via plate impact experiments coupled with high-speed internal digital image correlation and shadowgraphy techniques, the collapse of two pores in shock-compressed PMMA at stresses between 0.4 and 1 GPa. The results of these experiments provide new insights into shear localization in pore collapse, in addition to distinct interactions between pores. Shadowgraphy measurements reveal novel, direct visualization of shear band development and crack evolution from pore surfaces. Spacing between adiabatic shear bands is measured over a range of impact stresses and is predicted accurately by the Grady–Kipp model. Pore interactions are found to effect a transition in the impact stress threshold at which different failure mechanisms initiate and are also found to possibly influence preferential sites for shear cracking. Throughout the study, numerical and theoretical models are leveraged to understand shear localization behavior. The role of baroclinicity and wave interactions between the pores is used to elucidate interaction mechanisms between pores.

Lawlor, Barry P. [California Institute of Technolo↗
DOE OSTI2024

Numerical simulations of laser-driven experiments of ion acceleration in stochastic magnetic fields

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ΔE of ~ 264 keV for 242 MeV calcium ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamics of a relativistic jet through magnetized media

The merger of two neutron stars (NSs) produces the emission of gravitational waves, the formation of a compact object surrounded by a dense and magnetized environment. If the binary undergoes delayed collapse a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB), may be launched. The interaction of the jet with the environment has been shown to play a major role in shaping the structure of the outflow that eventually powers the gamma-ray emission. In this paper, we present a set of 2.5D RMHD simulations that follow the evolution of a relativistic non-magnetized jet through a medium with different magnetization levels, as produced after the merger of two NSs. We find that the predominant consequence of a magnetized ambient medium is that of suppressing instabilities within the jet and preventing the formation of a series of collimation shocks. One implication of this is that internal shocks lose efficiency, causing bursts with low-luminosity prompt emission. On the other hand, the jet-head velocity and the induced magnetization within the jet are fairly independent of the magnetization of the ambient medium. Future numerical studies with a larger domain are necessary to obtain light curves and spectra in order to better understand the role of magnetized media.

(magnetohydrodynamics) MHD↗
DOE OSTI2024

Evolving Outer Heliosphere: Tracking Solar Wind Transients from 1 au to the VLISM with IBEX and Voyager 1

Abstract Interstellar Boundary Explorer (IBEX) observations of energetic neutral atom (ENA) fluxes from the heliosphere have greatly enriched our understanding of the interaction of the solar wind (SW) with the local interstellar medium (LISM). However, there has been recent controversy surrounding the inability of most ENA models to produce as high an intensity of ∼0.5–6 keV ENAs as IBEX observes at 1 au, especially as a function of time. In our previous study (E. J. Zirnstein et al.), we introduced a new model that utilizes a data-driven magnetohydrodynamic simulation of the SW–LISM interaction to propagate pickup ions through the heliosheath (HS) after they are nonadiabatically heated at the heliospheric termination shock. E. J. Zirnstein et al. only simulated and analyzed IBEX observations from the direction of Voyager 2. In this study, we expand our model to include fluxes from the direction of Voyager 1, as well as in the low-latitude part (middle) of the ribbon (10° below the ecliptic plane). We show that the model results at Voyager 1 are consistent with E. J. Zirnstein et al.’s results at Voyager 2 in terms of a secondary ENA source contribution of ≲20% from both directions. Our results in the middle of the ribbon also reproduce the data, when including a time-dependent secondary ENA source. Finally, we demonstrate with our simulation that three large pressure waves likely merged in the VLISM and were observed by Voyager 1 as “pf2,” while at least one of the wave’s effects in the HS was observed by IBEX as a brief enhancement in ENA flux in early 2016.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗
DOE OSTI2025

Surface and sub-surface porosity effects on the initial free surface expansion in shocked single crystal aluminum

Defects and roughness at a material’s surface can allow for the formation of a material jet when a shockwave reaches a surface. The formation of a jet can lead to significant material damage of the shocked material, as well as the downstream material that the jet impacts. Here, we utilize planar shocks in single crystal aluminum to assess initial jetting formation for surface notches and sub-surface porosity. In the case of surface notches, we find that changing the width of the notch can cause the shock breakout response to transition from a standard jetting mechanism to an atomic atomization at high velocities. This transition occurs due to the mechanisms associated with lateral relaxations within the notch and how this interacts with the shock focusing event that induces jetting. For sub-surface pores, the pore collapse induces a re-shock, localized near the pore, that is stronger than the initial shock, but unsupported. Changing the pore dimensions and location relative to the surface influences the strength of the shock that emanates from its collapse when it reaches the surface, leading to a variety of different breakout conditions including jetting and localized material failure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of pre-existing heterogeneities in materials under shock and spall

There has been a challenge for many decades to understand how heterogeneities influence the behavior of materials under shock loading, eventually leading to spall formation and failure. Experimental, analytical, and computational techniques have matured to the point where systematic studies of materials with complex microstructures under shock loading and the associated failure mechanisms are feasible. This is enabled by more accurate diagnostics as well as characterization methods. As interest in complex materials grows, understanding and predicting the role of heterogeneities in determining the dynamic behavior becomes crucial. Early computational studies, hydrocodes, in particular, historically preclude any irregularities in the form of defects and impurities in the material microstructure for the sake of simplification and to retain the hydrodynamic conservation equations. Contemporary computational methods, notably molecular dynamics simulations, can overcome this limitation by incorporating inhomogeneities albeit at a much lower length and time scale. This review discusses literature that has focused on investigating the role of various imperfections in the shock and spall behavior, emphasizing mainly heterogeneities such as second-phase particles, inclusions, and voids under both shock compression and release. Pre-existing defects are found in most engineering materials, ranging from thermodynamically necessary vacancies, to interstitial and dislocation, to microstructural features such as inclusions, second phase particles, voids, grain boundaries, and triple junctions. This literature review explores the interaction of these heterogeneities under shock loading during compression and release. Systematic characterization of material heterogeneities before and after shock loading, along with direct measurements of Hugoniot elastic limit and spall strength, allows for more generalized theories to be formulated. Further, continuous improvement toward time-resolved, in situ experimental data strengthens the ability to elucidate upon results gathered from simulations and analytical models, thus improving the overall ability to understand and predict how materials behave under dynamic loading.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Global Hybrid Simulations of Interaction Between Interplanetary Rotational Discontinuity and Bow Shock/Magnetosphere: Can Ion-Scale Magnetic Reconnection be Driven by Rotational Discontinuity Downstream of Quasi-Parallel Shock?

Ion-scale magnetic reconnection has been observed downstream of the terrestrial quasi-parallel (Q-∥) shock. Whether it is driven by interplanetary discontinuities or turbulent Q-∥ shock, however, is unclear. Using three-dimensional global hybrid simulation, we investigate the generation of magnetic reconnection downstream of the Q-∥ shock, while an interplanetary rotational discontinuity (RD) is launched to the bow shock. Cases with various solar wind Alfvén Mach numbers, M A = 3.0 to 8, and propagation directions of the RD are presented. The propagation direction n is assumed to be in the GSE xz plane and pointing earthward, with n = (-sin(θ 12 /2),0,-cos(θ 12 /2)), where θ 12 is the angle between the upstream (B1) and downstream (B2) magnetic fields across the transmitted RD. It is found that magnetic reconnection occurs inside the RD downstream of the Q-∥ shock, forming flux ropes extending along the dawn-dusk direction about tens of ion inertial lengths. Large-amplitude low-frequency waves originated from the Q-∥ shock lead to the bending and squeezing of the field lines around the RD, which play an important role in triggering reconnection inside the RD. As the RD impacts the dayside magnetopause, magnetopause reconnection takes place between the field lines behind the RD and geomagnetic field lines. Nevertheless, no reconnection is found downstream of the Q-∥ shock itself or outside the RD in the magnetosheath. The existent and structure of reconnection in the magnetosheath are found to strongly depend on the parameters M A and n. Our simulation shows that ion-scale magnetic reconnection is driven by an external driver in the form of the compression of an RD around the bow shock and in the magnetosheath, rather than caused by the turbulent Q-∥ shock alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗