Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high explosive experiments”

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 55 records · Page 3

The Multiphysics on Advanced Platforms Project

In 2015, the Lawrence Livermore National Laboratory started development of next-generation multiphysics simulation capabilities for the National Nuclear Security Administration under the Advanced Technologies Development and Mitigation (ATDM) element of the Advanced Simulation and Computing program in collaboration with the Exascale Computing Project (ECP). A key driver for this effort across the NNSA tri-lab was the emergence of advanced high performance computing (HPC) architectures based on heterogeneous compute capabilities, including GPU based systems, as part of the national drive toward exascale computing platforms at multiple Department of Energy (DOE) facilities. Developing a multiphysics code capable of meeting the various simulation needs of the NNSA as defined by the current generation of integrated codes (or ICs), initially developed as part of the Accelerated Strategic Computing Initiative (ASCI) program beginning in 1996, and able to scale to the current 100 petaflop class pre-exascale systems, as well the forthcoming exaflop class computers, is a daunting challenge. To accomplish this ambitious goal, LLNL has embraced two key themes: use of high-order numerical methods and a modular approach to code development. The LLNL next generation effort is organized under the Multi-Physics on Advanced Platforms Project (MAPP). A foundational component of MAPP is the Axom computer science (CS) toolkit which provides infrastructure for the development of modular, performance portable, multi-physics application codes. MARBL is a next-generation application code built on the Axom base to address the modeling needs of the high energy density physics (HEDP) community for simulating high-explosive, magnetic or laser driven experiments such as inertial confinement fusion (ICF), pulsed-power magneto-hydrodynamics (MHD), equation of state (EOS) and material strength studies as part of the NNSA’s stockpile stewardship program (SSP).

97 MATHEMATICS AND COMPUTING↗

Magnetic Flux Compression Concept for Nuclear Pulse Propulsion and Power

The desire for fast, efficient interplanetary transport requires propulsion systems having short acceleration times and very high specific impulse attributes. Unfortunately, most highly efficient propulsion systems which are within the capabilities of present day technologies are either very heavy or yield very low impulse such that the acceleration time to final velocity is too long to be of lasting interest, One exception, the nuclear thermal thruster, could achieve the desired acceleration but it would require inordinately large mass ratios to reach the range of desired final velocities. An alternative approach, among several competing concepts that are beyond our modern technical capabilities, is a pulsed thermonuclear device utilizing microfusion detonations. In this paper, we examine the feasibility of an innovative magnetic flux compression concept for utilizing microfusion detonations, assuming that such low yield nuclear bursts can be realized in practice. In this concept, a magnetic field is compressed between an expanding detonation driven diamagnetic plasma and a stationary structure formed from a high temperature superconductor (HTSC). In general, we are interested in accomplishing two important functions: (1) collimation of a hot diamagnetic plasma for direct thrust production; and (2) pulse power generation for dense plasma ignition. For the purposes of this research, it is assumed that rnicrofusion detonation technology may become available within a few decades, and that this approach could capitalize on recent advances in inertial confinement fusion ICF) technologies including magnetized target concepts and antimatter initiated nuclear detonations. The charged particle expansion velocity in these detonations can be on the order of 10 (exp 6)- 10 (exp 7) meters per second, and, if effectively collimated by a magnetic nozzle, can yield the Isp and the acceleration levels needed for practical interplanetary spaceflight. The ability to ignite pure fusion micro-bursts with reasonable levels of input energy is an equally challenging scientific problem. It remains to be seen, however, whether an effective ignition driver can be developed which meets the requirements for practical spaceflight application (namely high power density, compactness, low weight, and low cost). In this paper, system level performance and design issues are examined including generator performance, magnetic flux compression processes, magnetic diffusion processes, high temperature superconductor (HTSC) material properties, plasmadynamic processes, detonation plasma expansion processes, magnetohydrodynamic instabilities, magnetic nozzle performance, and thrust production performance. Representative generator performance calculations based on a simplified skin layer formulation are presented as well as the results of exploratory small-scale laboratory experiments on magnetic flux diffusion in HTSC materials. In addition, planned follow-on scientific feasibility experiments are described which utilize high explosive detonations and high energy gas discharges to simulate the plasma conditions associated with thermonuclear micro-detonations.

Litchford, Ronald J.↗

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Numerical calculation of shock-induced initiation of detonations

Results of numerical calculations of the impact of steel cylinders and spheres on the plastic bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor acceptor cell method to compute mixed cell fluxes. The parameters in the Forest Fire burn model are developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.

Cort, G. E.↗

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

The L = 6.6 Oosik barium plasma injection experiment and magnetic storm of March 7, 1972

A high-explosive shaped charge vaporizing a hollow conical liner of Ba metal and producing a fast field-aligned jet of plasma was detonated at high altitude, during a quiescent phase of a magnetic storm initiated by an ssc 10 hrs prior to the experiment, in an attempt to trace out and observe the dynamics of an auroral field line in the magnetosphere. Observations offer evidence for an upward Birkeland current sheet at the poleward edge of the auroral spiral of 8 x 10 to the minus second power A/m. Unusual features of the substorm leave open the possibility that the plasma injection may have triggered it.

Wescott, E. M.↗

L = 1.24 conjugate magnetic field line tracing experiments with barium shaped charges

Description of three experiments involving the injection of barium ions into magnetic flux tubes with the aid of high-explosive shaped charges with hollow conical liners of barium metal. In these experiments (called Alco, Bubia, and Loro, respectively) the explosive charges were detonated at altitudes above 500 km, producing jets of barium plasma with initial velocities ranging from 8 to 20 km/sec. The most interesting result of the experiments, which were carried out near L = 1.24 and were successful in tracing an entire field line some 7000 km to the conjugate ionosphere, was the observation that for Loro and probably Alco the direction and rate of drift of the two ends of the field line were identical, whereas those for Bubia differed significantly. In the case of the Bubia event significant differences in magnitude and direction at the conjugate points are noted, which lead to the conclusion that the field line could not have been equipotential during the interval of observation.

Wescott, E. M.↗

Development of spinning-disk solid sample delivery system for high-repetition rate x-ray free electron laser experiments

X-ray free-electron lasers (XFELs) deliver intense x-ray pulses that destroy the sample in a single shot by a Coulomb explosion. Experiments using XFEL pulse trains or the new generation of high-repetition rate XFELs require rapid sample replacement beyond those provided by the systems now used at low repletion-rate XFELs. We describe the development and characterization of a system based on a spinning disk to continuously deliver a solid sample into an XFEL interaction point at very high speeds. We tested our system at the Linac Coherent Light Source and European XFEL hard x-ray nano-focus instruments, employing it to deliver a 25 μm copper foil sample, which can be used as a gain medium for stimulated x-ray emission for the proposed x-ray laser oscillator.

47 OTHER INSTRUMENTATION↗

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↗

Detection of the Large Surface Explosion Coupling Experiment by a Sparse Network of Balloon-Borne Infrasound Sensors

In recent years, high-altitude infrasound sensing has become more prolific, demonstrating an enormous value especially when utilized over regions inaccessible to traditional ground-based sensing. Similar to ground-based infrasound detectors, airborne sensors take advantage of the fact that impulsive atmospheric events such as explosions can generate low frequency acoustic waves, also known as infrasound. Due to negligible attenuation, infrasonic waves can travel over long distances, and provide important clues about their source. Here, we report infrasound detections of the Apollo detonation that was carried on 29 October 2020 as part of the Large Surface Explosion Coupling Experiment in Nevada, USA. Infrasound sensors attached to solar hot air balloons floating in the stratosphere detected the signals generated by the explosion at distances 170–210 km. Three distinct arrival phases seen in the signals are indicative of multipathing caused by the small-scale perturbations in the atmosphere. We also found that the local acoustic environment at these altitudes is more complex than previously thought.

47 OTHER INSTRUMENTATION↗

Extracting the Breakout Distance from the ECOT Trajectories: Gaussian Process Regression Approach

Enhanced Corner Turning (ECOT) experiments provide an important metric of performance of high explosive (HE) formulations. The breakout distance is a single scalar value that characterizes the corner turning efficiency of an HE. Extracting the breakout distance from the raw ECOT results, whether experimental or simulated, is a conceptually straightforward procedure which, however, is non-unique, especially in the presence of noise. More specifically, this procedure involves numerical smoothing and selecting particular values for parameters of this smoothing introduces human bias. In this work, we propose to use the Gaussian process regression to analyze ECOT results. This analysis involves the effective smoothing of the data, thus allowing for accurate extraction of the breakout distance. Most importantly, the parameters of this smoothing can be inferred from the ECOT data itself, rendering the approach effectively parameter-free and thus diminishing the human bias. An additional benefit of the Gaussian process regression, being a statistical inference method, is that not just the value of the breakout distance, but also its confidence interval can be extracted from the data. This report introduces the Gaussian process regression, as applied to ECOT, and demonstrates its usefulness by extracting the breakout distances for a selection of experimental and simulated data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A15 Nb 3 Si: a ‘high’ T c superconductor synthesized at a pressure of one megabar and metastable at ambient conditions

A15 Nb 3 Si is, until now, the only 'high' temperature superconductor produced at high pressure (~110 GPa) that has been successfully brought back to room pressure conditions in a metastable condition. Based on the current great interest in trying to create metastable-at-room-pressure high temperature superconductors produced at high pressure, we have restudied explosively compressed A15 Nb 3 Si and its production from tetragonal Nb 3 Si. First, diamond anvil cell pressure measurements up to 88 GPa were performed on explosively compressed A15 Nb 3 Si material to trace T c as a function of pressure. T c is suppressed to ~5.2 K at 88 GPa. Then, using these T c (P) data for A15 Nb 3 Si, pressures up to 92 GPa were applied at room temperature (which increased to 120 GPa at 5 K) on tetragonal Nb 3 Si. Measurements of the resistivity gave no indication of any A15 structure production, i.e. no indications of the superconductivity characteristic of A15 Nb 3 Si. This is in contrast to the explosive compression (up to P ~ 110 GPa) of tetragonal Nb 3 Si, which produced 50%–70% A15 material, T c = 18 K at ambient pressure, in a 1981 Los Alamos National Laboratory experiment. This implies that the accompanying high temperature (1000 °C) caused by explosive compression is necessary to successfully drive the reaction kinetics of the tetragonal → A15 Nb 3 Si structural transformation. Our theoretical calculations show that A15 Nb 3 Si has an enthalpy vs the tetragonal structure that is 70 meV atom –1 smaller at 100 GPa, while at ambient pressure the tetragonal phase enthalpy is lower than that of the A15 phase by 90 meV atom –1 . Furthermore, the fact that 'annealing' the A15 explosively compressed material at room temperature for 39 years has no effect shows that slow kinetics can stabilize high pressure metastable phases at ambient conditions over long times even for large driving forces of 90 meV atom –1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cookoff of Powdered and Pressed Explosives using a Micromechanics Pressurization Model

Cookoff experiments of powdered and pressed TATB-based plastic bonded explosives (PBXs) have been modeled using a pressure-dependent universal cookoff model (UCM) in combination with a micromechanics pressurization (MMP) model described in a companion paper. The MMP model is based on the accumulation of decomposition gases at nucleation sites that load the surrounding TATB crystals and binder. This is the first cookoff model to use an analytical mechanics solution for compressibility and thermal expansion to describe internal pressurization caused by both temperature and decomposition occurring within closed-pore explosives. This approach produces more accurate predictions of ignition time and pressurization within high-density explosives than simple equation-of-state models. The current paper gives details of the reaction chemistry, model parameters, predicted uncertainty, and validation using experiments from multiple laboratories with errors less than 6 %. The UCM/MMP model framework gives more accurate thermal ignition predictions for high density explosives that are initially impermeable to decomposition gases.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Abbreviated Technical Report: Experimentally Interrogating Detonation Chemistry on Sub-Nanosecond to Nanosecond Timescales

Direct experimental measurement of chemical reactions during high explosive detonation remains challenging. Theory and modeling have long preceded experiment in the fundamental physical and chemical kinetic properties of detonation, and experimentation at the relevant timescales are needed to both validate models and provide fundamental understanding of detonation. In this LDRD-ER project, two approaches, x-ray diffraction and core-level x-ray Raman, were developed and used to further experimental capabilities to address this gap. We further developed dynamic x-ray diffraction to directly detect nanodiamond formation during detonation, providing experimental data towards resolving longstanding controversy in the scientific literature, and although the full kinetics have not yet been fully mapped, diamond diffraction appears on the same timescales as detonation soot formation. In the second research thrust, we have developed core-level x-ray Raman for use with high explosives. This technique provides information analogous to x-ray absorption spectroscopy and electron energy loss spectroscopy, but uses inelastic scattering of hard x-rays that can interrogate chemistry around light elements much deeper into the material. The low cross section and requisite high solid angle collection have hindered its use for ultra-fast spectroscopy. We developed and tested a high-q spectrometer which will substantially increase cross section and signal-to-noise, showing this method will also not dramatically alter, compared to x-ray absorption, the most discriminating spectral features of C, N, and O from various high explosives and expected detonation products. We have also used x-ray Raman combined with OCEAN electronic structure calculations to explore dynamic photodegradation mechanisms in PETN and CL-20 explosives. This provides a pathway towards implementing capability to dynamically explore chemistry at an x-ray free electron laser.

36 MATERIALS SCIENCE↗

Nonthermal electron and ion acceleration by magnetic reconnection in large laser-driven plasmas

Magnetic reconnection is a fundamental plasma process that is thought to play a key role in the production of nonthermal particles associated with explosive phenomena in space physics and astrophysics. Experiments at high-energy-density facilities are starting to probe the microphysics of reconnection at high Lundquist numbers and large system sizes. We have performed particle-in-cell (PIC) simulations to explore particle acceleration for parameters relevant to laser-driven reconnection experiments. We study particle acceleration in large system sizes that may be produced soon with the most energetic laser drivers available, such as at the National Ignition Facility. In these conditions, we show the possibility of reaching the multi-plasmoid regime, where plasmoid acceleration becomes dominant. Our results show the transition from X point to plasmoid-dominated acceleration associated with the merging and contraction of plasmoids that further extend the maximum energy of the power-law tail of the particle distribution for electrons. We also find for the first time a system-size-dependent emergence of nonthermal ion acceleration in driven reconnection, where the magnetization of ions at sufficiently large sizes allows them to be contained by the magnetic field and energized by direct X point acceleration. For feasible experimental conditions, electrons and ions can attain energies of ϵ max , e / k B T e > 100 and ϵ max , i / k B T i > 1000 . Using PIC simulations with binary Monte Carlo Coulomb collisions, we study the impact of collisionality on plasmoid formation and particle acceleration. The implications of these results for understanding the role reconnection plays in accelerating particles in space physics and astrophysics are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New Design Yields Robust Large-Area Framing Camera

Manuscript to be submitted to Review of Scientific Instruments to determine interest in possible publication. Website: https://aip.scitation.org/rsi/authors/manuscript High-speed imaging is a core diagnostic capability for the Nevada National Security Site (NNSS). Data from experiments that measure the behavior of materials under extreme conditions are essential for the stockpile mission. For example, during subcritical experiments, materials that are shocked and compressed under high-explosive drive change behavior rapidly over timescales measured in millionths or even billionths of a second. Characterizing and understanding this behavior requires radiographic and/or optical imaging to track motion on these timescales. Framing cameras provide a discrete series of images over a short period of time in a manner that mimics a high-speed movie camera, but with individual shutter times that must be extremely short in order to capture freeze-frame images of rapidly evolving phenomena such as high-explosive shock-driven ejecta, dynamic compression of metals, and high-velocity fluid flow.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗