SEARCH · Search NASA
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.
Gas gun shock-to-detonation experiments provide data for calibrating HE burn models. At LANL the principal diagnostic is a magnetic gauge package that provides data on the lead shock trajectory, which is used to determine the run-distance to detonation, and the lead shock speed with run distance, which together with jump-off velocity of the velocity gauges determines Hugoniot locus points for the reactants.
ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.
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.
Explosion is the most extreme case of thermal runaway of lithium-ion (Li-ion) batteries. In this study, explosion dynamics of large-format Li-ion cells are investigated experimentally and numerically. Overcharge-to-explosion tests are conducted on 40 Ah Li-ion cells with Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 cathode. Based on the explosion physics, shockwave and detonation models are used to characterize the shock effect of the cell explosion and evaluate the explosion equivalent. Von Neumann peak is observed on the pressure curves, and the shockwave velocity is supersonic at this time; the unwrinkled spherical flame phenomenon observed in the experiment indicates that it is detonation. Additionally, a geometric model is established based on the real testing scenario, and the explosion behavior is numerically studied. The characteristics of the explosion dynamics process are interpreted and the propagation mechanism of the shockwave are revealed; the deflagration to detonation transition (DDT) phenomenon in this process is caused by the formation of “hot spots”, and the explosion of the cells eventually turns into stable combustion. Here this study fills the gap in the research on thermal runaway of Li-ion cells, especially in the extreme cases of fire and explosion, and provide useful guidance for battery safety.
The cylinder expansion experiment is a well-established performance test for condensed explosives that is utilized routinely to determine pressure-energy-volume relationships for detonation products. The modern cylinder test employs optical interferometric techniques to measure velocity, as opposed to older realizations where streak cameras were commonplace. Despite their widespread use, questions sometimes remain as to what kind of data the velocity diagnostics in a cylinder test provide along with how to interpret such measurements and relate them to physical quantities of interest. Here in this study, equations are derived that fully describe the kinematics of the cylinder wall during expansion. These equations are then applied to experimental as well as synthetic data generated via a hydrodynamic simulation in order to verify the mathematical framework developed.
In this study, we present and analyze a near-infrared (NIR) spectrum of the underluminous Type Ia supernova SN 2020qxp/ASASSN-20jq obtained with NIRES at the Keck Observatory, 191 days after B-band maximum. The spectrum is dominated by a number of broad emission features, including the [Fe ιι] at 1.644 μm, which is highly asymmetric with a tilted top and a peak redshifted by ≈2000 km s -1 . In comparison with 2D non-LTE synthetic spectra computed from 3D simulations of off-center delayed-detonation Chandrasekhar-mass (M ch ) white dwarf (WD) models, we find good agreement between the observed lines and the synthetic profiles, and are able to unravel the structure of the progenitor’s envelope. We find that the size and tilt of the [Fe ιι] 1.644 μm profile (in velocity space) is an effective way to determine the location of an off-center delayed-detonation transition (DDT) and the viewing angle, and it requires a WD with a high central density of ~4 × 10 9 g cm -3 . We also tentatively identify a stable Ni feature around 1.9 μm characterized by a “pot-belly” profile that is slightly offset with respect to the kinematic center. In the case of SN 2020qxp/ASASSN-20jq, we estimate that the location of the DDT is ~0.3M WD off center, which gives rise to an asymmetric distribution of the underlying ejecta. We also demonstrate that low-luminosity and high-density WD SN Ia progenitors exhibit a very strong overlap of Ca and 56 Ni in physical space. This results in the formation of a prevalent [Ca ιι] 0.73 μm emission feature that is sensitive to asymmetry effects. Our findings are discussed within the context of alternative scenarios, including off-center C/O detonations in He-triggered sub-M Ch WDs and the direct collision of two WDs. Snapshot programs with Gemini/Keck/Very Large Telescope (VLT)/ELT-class instruments and our spectropolarimetry program are complementary to mid-IR spectra by the James Webb Space Telescope (JWST).
In this work, a high-throughput experimental setup was used to characterize initiation threshold and growth to detonation in the explosives hexanitrostilbene (HNS) and pentaerythritol tetranitrate (PETN). The experiment sequentially launched an array of laser-driven flyers to shock samples arranged in a 96-well microplate geometry, with photonic Doppler velocimetry diagnostics to characterize flyer velocity and particle velocity at the explosive–substrate interface. Vapor-deposited films of HNS and PETN were used to provide numerous samples with various thicknesses, enabling characterization of the evolution of growth to detonation. One-dimensional hydrocode simulations were performed with reactions disabled to illustrate where the experimental data deviate from the predicted inert response. Prompt initiation was observed in 144 μm thick HNS films at flyer velocities near 3000 m/s and in 125 μm thick PETN films at flyer velocities near 2400 m/s. This experimental setup enables rapid quantification of the growth of reactions in explosive materials that can reach detonation at sub-millimeter length scales. These data can subsequently be used for parameterizing reactive burn models in hydrocode simulations, as discussed in Paper II [D. E. Kittell, R. Knepper, and A. S. Tappan, J. Appl. Phys. 131, 154902 (2022)].
Despite the importance of Type Ia supernovae (SNe Ia) throughout astronomy, the precise progenitor systems and explosion mechanisms that drive SNe Ia are still unknown. An explosion scenario that has gained traction recently is double detonation, in which an accreted shell of He detonates and triggers a secondary detonation in the underlying white dwarf. Our research presents a number of high-resolution, multidimensional, full-star simulations of thin-He-shell, sub-Chandrasekhar-mass white dwarf progenitors that undergo a double detonation. This suite of thin-shell progenitors incorporates He shells that are thinner than those in previous multidimensional studies. We confirm the viability of the double detonation across a range of He-shell parameter space, as well as present bulk yields and ejecta profiles for each progenitor. The yields obtained are generally consistent with previous works and indicate the likelihood of producing observables that resemble SNe Ia. The dimensionality of our simulations allow us to examine features of the double detonation more closely, including the details of the off-center secondary ignition and asymmetric ejecta. We find considerable differences in the high-velocity extent of postdetonation products across different lines of sight. The data from this work will be used to generate predicted observables and may further support the viability of the double detonation scenario as an SN Ia channel, as well as show how the properties of the progenitor or viewing angle may influence trends in observable characteristics.
Embedded velocity gauges are standardly used at LANL on gas gun experiments of a shock to-detonation transition. The experiments are intended to be one-dimensional with a planar shock driven by the projectile from a gas gun propagating along the axis of the HE cylinder. As disussed by Gustavsen et al. [1999], the accuracy of the gauges depends on the alignment of the directions of shock propagation, magnetic field and the active gauge element. A source of inaccuracies is from tilt or misalignment angle at impact between the projectile surface and the HE surface; typically a few milliradians. This leads to the shock propagating at a slight angle to the axis of the HE cylinder and can affect the shock rise time of the velocity gauges, the time offset between the tracker gauges for the trajectory of the lead shock, and the gauge positions in the direction of propagation. It is important to note that tilt is characterized by 2 parameters; the magnitude of the angle relative to the HE axis and the direction in the plane transverse to the HE axis. Here we derive the uncertainties from tilt to leading order in the tilt angle.
This report documents the design and data for four related experiments known as the Wedding Cake shots. It is a revised version of LA-UR-17-28352. The basic designs of these four shots are the same. A high-explosive plane-wave lens launches a flyer towards an acceptor material with a smaller diameter. This acceptor is either another piece of high-explosive or a piece of polyethylene. A detonation or shock wave propagates through this material and launches a second flyer. The velocities of both flyers are recorded using Photon Doppler Velocimetry (PDV). The overall shape of the shots is reminiscent of a tiered wedding cake and, consequently, these shots are know as the Wedding Cake shots. Details of the shot assembly are shown in the Wedding Cake Assembly Drawing and component drawings that accompany this report.
Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.
Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.
As the operational time window of experimental rotating detonation engines (RDEs) is expanded and the technology matures toward integration within gas turbines, monitoring techniques must evolve to offer computationally efficient and highly time-resolved diagnostics. In this study, computer vision object detection methodology that seeks to reduce data processing time and calculate wave velocity within drastically reduced time intervals as compared to traditional high-frame-rate RDE images analysis techniques is proposed. The adapted you-only-look-once object detection network is trained to detect individual detonation waves within single down-axis RDE images. The wave location and rotational direction detected within a frame are tracked through a series of high-speed images to calculate the frame-to-frame wave velocity with the time-step resolution of $\mathrm{20 μs}$ across a series of frames. The analysis of the annotation box size and image linearization effects is presented, demonstrating the lowest frame-to-frame velocity total uncertainty of $\mathrm{±3.8\%}$ and the highest classification speed of 9.5 frames per second using linearized images. Linearized images “unwrap” the RDE annulus pixel region to a reduced image size. Here, this new method offers great reductions in data processing times and unsteady detonation behavior insight at intervals more comparable to the timescales of detonation wave interactions via the application of machine learning to experimental RDE data.
Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.
The ability to predict post-ignition explosive response is relevant to explosives safety, particularly weapons assembly/disassembly operations that must operate under a sufficiently conservative safety basis to meet risk mitigation criteria for inadvertent nuclear detonation (IND). In particular, hazard determinations hinge on whether the PBX 9501 main charge can ultimately attain detonation from an initial, non-shock, ignition stimulus. The CISME experiment uses a spherical explosive charge, which is thermally ignited in the center using a fiber-delivered infrared laser pulse. This ignition mechanism - i.e. a localized thermal hotspot - is representative of the ignition-causing insults typical in accident scenarios, including low-velocity impact, drop, skid, and punch events. In fiscal year 2018 CISME experiments were conducted with pristine PBX 9501 spheres, confined and unconfined. Sub-detonative, low-violence reactions were observed. Here we report experiments conducted since 2018. The explosive charge was varied: experiments were conducted with room temperature pour-density HMX powder and prills of PBX 9501, and heated, thermally damaged PBX 9501. Experiments were conducted in a confined configuration, with two 4 mm-thick stainless-steel hemispherical shells bolted around the sphere. Detonation, via DDT (deflagration-to-detonation transition) was observed in the 3-inch diameter powdered HMX and the 6-inch diameter heated and thermally-damaged PBX 9501 tests. No transition to DDT was observed in smaller diameters for those configurations, and no DDT was observed for pour-density prills of PBX 9501 of any diameter. The original test series demonstrated that the response of pristine, undamaged PBX 9501 to a central ignition is relatively benign. This new test series shows that increased porosity, lack of polymer binders and thermal damage can greatly increase the violence of the post-ignition response, up to and including detonation. These results bridge the explosive response of PBX 9501 from sub-detonative through detonation as a function of diameter and provide valuable test cases for modeling efforts that strive to capture the DDT threshold. By exploring explosive response in a more representative 3D geometry, we can also benchmark these new results to more plentiful DDT datasets generated from testing in 1D tube confinement (the historical test configuration)–the configuration upon which almost all DDT models were developed.
In energetic materials, the localization of energy into “hotspots” is known to dictate the initiation of chemical reactions and detonation. Recent all-atom simulations have shown that more energy is localized as internal potential energy (PE) than can be inferred from the kinetic energy (KE) alone. The mechanisms associated with pore collapse and hotspot formation are known to depend on pore geometry and dynamic material response such as plasticity. Therefore, we use molecular dynamics (MD) simulations to characterize shock-induced pore collapse and the subsequent formation of hotspots in 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a highly anisotropic molecular crystal, for various defect shapes, shock strengths, and crystallographic orientations. We find that the localization of energy as PE is consistently larger than the KE in cases with significant plastic deformation. Additionally, an analysis of MD trajectories reveals the underlying molecular- and crystal-level processes that govern the effect of orientation and pore shape on PE localization. We find that the regions of highest PE relate to the areas of maximum plastic deformation, while KE is maximized at the point of impact. Comparisons against octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) reveal less energy localization in TATB, which could be a contributing factor to the latter’s insensitivity.
Understanding the formation and long-timescale evolution of atmospheric plasmas produced by ultrashort, long-wavelength infrared (LWIR) pulses is an important but partially understood problem. Of particular interest are plasmas produced in air with a peak laser intensity ~10 12 W/cm 2 , the so-called clamping intensity observed in LWIR atmospheric guiding experiments where tunneling and multi-photon ionization operative at near-IR or shorter wavelengths are inoperative. We find that avalanche breakdown on the surface of aerosol (dust) particles can act to seed the breakdown of air observed above the 200 GW/cm 2 threshold when a train of 3 ps 10.6 μm laser pulses separated by 18 ps is used. The breakdown first appears at the best focus but propagates backward toward the focusing optic as the plasma density approaches critical density and makes forward propagation impossible. The velocity of the backward propagating breakdown can be as high as 10 9 cm/s, an order of magnitude greater than measured with ns pulse-produced breakdown, and can be explained rather well by the so-called breakdown wave mechanism. Transverse plasma expansion with a similar velocity is assisted by UV photoionization and is observed as a secondary longitudinal breakdown mechanism in roughly 10% of the shots. When a cm-size, TW power beam is propagated, interception of aerosol particles is guaranteed and several (40 cm -3 ) breakdown sites appear, each initially producing a near-critical density plasma. On a 10 ns–1 μs timescale, shockwaves from each site expand radially and coalesce to produce a large hot gas channel. The radial velocity of the expansion agrees well with the prediction of the blast wave theory developed for ultrafast atmospheric detonations.
A remote-operated apparatus for testing the detonation sensitivity of energetic materials is detailed. Using an air ram and rotating disk, the normal force and transverse velocity of the impact plane are controlled independently, enabling the exploration of varying impact conditions over a wide parameter space. A microcontroller local to the apparatus is used to automate apparatus operation and ensure temporal alignment of the impacting ram head with the rotating disk. Calculation of the firing parameters and issuing of operational commands are handled by a remote computer and relayed to the local microcontoller for execution at the hardware level. Finally, impact forces are taken from fast strain measurements obtained from gauges incorporated into the ram head. Infrared imaging of explosive samples provides insight to the peak thermal temperatures experienced at the sample surface during the impact event.