A new constitutive model for underground chemical explosions in partially saturated porous rock
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Discrimination of underground explosions from naturally occurring earthquakes and other anthropogenic sources is one of the fundamental challenges of nuclear explosion monitoring. In an operational setting, the number of events that can be thoroughly investigated by analysts is limited by available resources. The capability to rapidly screen out events that can be robustly identified as not being explosions is, therefore, of great potential benefit. Nevertheless, possible mis-classification of explosions as earthquakes currently limits the use of screening methods for verification of test-ban treaties. Moment tensors provide a physics-based classification tool for the characterization of different seismic sources and have enabled the advent of new techniques for discriminating between earthquakes and explosions. Following normalization and projection of their six-degree vectors onto the hypersphere, existing screening approaches use spherically symmetric metrics to determine whether any new moment tensor may have been an explosion. Here, we show that populations of moment tensors for both earthquakes and explosions are anisotropically distributed on the hypersphere. Distributions possessing elliptical symmetry, such as the scaled von Mises–Fisher distribution, therefore provide a better description of these populations than the existing spherically symmetric models. We describe a method that uses these elliptical distributions in combination with a Bayesian classifier to achieve successful classification rates of 99 per cent for explosions and 98 per cent for earthquakes using existing catalogues of events from the western United States. The 1983 May 5 Crowdie underground nuclear test and 2018 July 20 DAG-1 deep-borehole chemical explosion are the only two explosions out of 140 that are incorrectly classified. Application of the method to the 2006–2017 nuclear tests in the Democratic People’s Republic of Korea yields 100 per cent identification rates and we provide a simple routine MTid for general usage. The approach provides a means to rapidly assess the likelihood of an event being an explosion and can be built into monitoring workflows that rely on simultaneously assessing multiple different discrimination metrics.
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.
A research effort was conducted to investigate the effects of an unintentional internal detonation of energetic materials on a small transportation container called the “Reinforced Ammo Can” (RAC), which is currently approved for the transport of energetic materials up to 2 grams net explosive weight (NEW) onsite at Lawrence Livermore National Laboratory (LLNL). The project used experimental testing of the containers based on simulated results conducted by a Missouri S & T explosives engineering graduate student while on an internship working in explosives safety at LLNL. The proposed experimental testing consisted of 20 central detonation tests in separate Reinforced Ammo Can containers to verify their high explosive (HE) containment capability. RAC containers are required for transporting secondary explosives on site at LLNL Site 200 in quantities between 300 mg and 2 grams of Trinitrotoluene (TNT) equivalent NEW. Currently, the LLNL Environmental Safety and Health (ES&H) Manual states that the RAC containers are permissible for on-site transportation of secondary explosives up to 2 grams NEW by approved HE handlers with the approval of an Explosives Safety Expert (ESE). Simulated results indicated that tensile failure of the closure mechanism on the container is the most likely failure mode, and showed potential occurrence at the 2 gram NEW threshold. Experimental validation tests were conducted indicating failure/rupture of the containers lid with venting of the detonation products at the 2 gram NEW test threshold with an explosive mass safety factor of 1.3. Failure consisted of impact opening the container latch mechanism resulting in potential ejection of the lid depending on the orientation of the charge within the confinement. This analysis solely is an evaluation of the response of the containment if a detonation were to occur, it does NOT consider the probability of such an event.
Two balls were dropped approximately 150 m down an open borehole as part of an investigation into the aftershocks from a 51 ton TNT-equivalent chemical explosion. The ball drop sources, a concrete-filled basketball and a bowling ball, were intended to mimic falling debris within an underground explosion borehole and cavity. Falling debris or sloughing of the borehole and cavity walls is thought to be one of the sources of explosion aftershocks. The signals from the ball drop sources were recorded by the same seismic network that monitored the aftershocks of the chemical explosion, so the ball drop recordings were used as templates to search through the continuous stream of aftershocks and detect correlated signals as a proportion of all aftershocks. The search resulted in only four similar event recordings out of the thousands of detected aftershocks. A synthetic test showed that the small, high-frequency ball drop templates were limited to producing highly correlated detections for events occurring in a small region around the impact location at 150 m depth. The high-frequency templates were necessary since the small ball drop signals were only above the noise at frequencies greater than 10 Hz. Unfortunately, the impact location is where the borehole was filled with sand in preparation for the subsequent chemical explosion in the test series, making it 150 m above the explosion created cavity centered at 300 m depth. This puts the ball drop templates 150 m away from the cavity, where many more events from falling debris would be expected. Further synthetic tests confirmed events emanating from the cavity 150 m below the impact site would have correlated with the lower frequency template and been detected had the template been usable below 10 Hz. Therefore, we conclude that the four detected events were most probably falling debris within the open borehole that impacted the floor at 150 m depth just as the ball drop sources had, and that this type of explosion aftershock mechanism comprises less than 1% of all explosion aftershocks.
The Source Physics Experiment (SPE) series is a long-term NNSA research and development effort designed to improve U.S. arms control and nuclear nonproliferation verification and monitoring capabilities. The findings from the SPE will advance the United States’ nuclear explosion monitoring capabilities, particularly with respect to detection, discrimination and determination of yields associated with small nuclear explosions that can be lost amid the noisy seismo-acoustic background from other sources. The data generated from the SPE, a series of well-designed and recorded chemical explosions, will contribute to the development and validation of first-principles explosive source generated seismo-acoustic modeling codes. These codes will then facilitate the update of semi-empirical methods, currently based on historic test site data, such that key explosion observables can be reproduced, thus improving confidence in nuclear test monitoring in new areas and/or under novel emplacement conditions. The overall SPE project is comprised of both the development of the new explosion simulation codes and the chemical explosion test series. The chemical explosion test series will generate the empirical data required to both develop and validate the new simulation codes.
Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters—such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. The post-bounce temperature and density evolution of zero-age main-sequence 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio isotopes that are of particular interest to the next generation of gamma-ray astronomical observations: 43 K, 47 Ca, 44 Sc, 47 Sc, 48 V, 48 Cr, 51 Cr, 52 Mn, 59 Fe, 56 Co, 57 Co, and 57 Ni. Furthermore, these nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43 K, 47 Ca, 44 Sc, 47 Sc, and 59 Fe, those that are dependent on the progenitor structure are 48 V, 51 Cr, and 57 Co, and those that probe neither are 48 Cr, 52 Mn, 57 Ni, and 56 Co. We discuss the prospects of observing these radionuclides in supernova remnants.
SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.
The Dry Alluvium Geology (DAG) project was Phase II of the Source Physics Experiment and consisted of a series of four chemical explosive tests conducted in the same source hole on the Nevada National Security Site. This hole is located at 37.1146°N and -116.0693°W, with a surface elevation of 1,285.2 meters (m) (4,216.5 feet [ft]) above sea level. The first test (DAG-1) was conducted on July 20, 2018, at 16:51:52.67838 Coordinated Universal Time (UTC). The explosive source for DAG-1 was nitromethane initiated by a small plastic-bonded explosive (PBX) charge, detonated at the depth of 385.0 m (1,263.2 ft) below ground surface. DAG-1 had a trinitrotoluene (TNT) equivalent yield of 0.908 metric tons (2,002 pounds [lbs]). DAG-2 was conducted on December 19, 2018, at 18:45:56.92115 UTC. This test was the largest in the series, with a TNT equivalent yield of 50.997 metric tons (112,429 lbs). The explosive source for DAG-2 was nitromethane initiated by a small PBX charge, detonated at the depth of 299.8 m (983.6 ft) below ground surface. DAG-3 was conducted on April 27, 2019, at 15:49:01.84183 UTC. The explosive source for this test was nitromethane initiated by a small PBX charge, detonated at the depth of 149.9 m (492.0 ft) below ground surface. DAG-3 had a TNT equivalent yield of 0.908 metric tons (2,002 lbs). The final DAG test (DAG-4) was conducted on June 22, 2019, at 21:06:19.87632 UTC. The explosive source for DAG-4 was nitromethane initiated by a small PBX charge, detonated at the depth of 51.6 m (169.3 ft) below ground surface. DAG-4 had a TNT equivalent yield of 10.357 metric tons (22,833 lbs). The four tests were recorded by an extensive set of instrumentation that included sensors both at near-field (less than 200 m) and far-field (200 m or greater) distances. The near-field instruments consisted of three-component (3C) accelerometers installed at various depths ranging from 51.6 to 385 m (169.3 to 1,263.1 ft) below ground surface in boreholes positioned around the source hole, and arrays of single-component and 3C accelerometers on the surface. The far-field network comprised a variety of seismic and acoustic sensors, including short-period geophones, broadband seismometers, and 3C accelerometers at distances of 200 m to 400 kilometers. In addition, the DAG-2, DAG-3, and DAG-4 explosions were recorded by a temporary array of 496 geophones arranged in a densely spaced grid pattern known as “Large N.” This report coincides with the release of these data for analysts and organizations that are not participants in this program. This report describes the four DAG tests and the various types of near-field, far-field, and other data that are available. Assembled data sets are accessible through: Incorporated Research Institutions for Seismology, Data Management Center 1408 NE 45th Street, Suite 201, Seattle, Washington 98105 USA. www.iris.washington.edu
Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.
Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities, especially development of convection, will determine the morphology of the explosion. We have performed 3D simulations using spherical-polar coordinates covering a reduced angular extent (90° computational domain), and with angular resolutions of 2°, 1°, 1/2°, and 1/4°, to study the development of turbulent-like flows in core-collapse supernova explosions on a time scale of order 100 ms. We have employed the multi-physics Chimera code that includes detailed nuclear physics and spectral neutrino transport. Coarse resolution models do not develop an inertial range, presumably due to the bottleneck effect, such that the energy is prevented from cascading down to small scales and tends to accumulate at large scales. High-resolution models instead, start to recover the k −5/3 scaling of Kolmogorov's theory. Stochasticity and few simulation samples limit our ability to predict the development of explosions. Over the simulated time period, our models show no clear trend in improving (or diminishing) conditions for explosion as the angular resolution is increased. However, we find that disordered flow provides an effective pressure, as characterized by the Reynolds stress, behind the shock. This contribution, ~40%–50% of the thermal pressure, aids shock revival and thus the development of the explosion. Finally, we show that the kinetic energy power spectrum of reduced angular extent and full 4π models are consistent, thus indicating that a 90° computational domain is an adequate configuration to study the character of turbulence in CCSNe.
ABSTRACT We present observations of SN 2020fqv, a Virgo-cluster type II core-collapse supernova (CCSN) with a high temporal resolution light curve from the Transiting Exoplanet Survey Satellite (TESS) covering the time of explosion; ultraviolet (UV) spectroscopy from the Hubble Space Telescope (HST) starting 3.3 d post-explosion; ground-based spectroscopic observations starting 1.1 d post-explosion; along with extensive photometric observations. Massive stars have complicated mass-loss histories leading up to their death as CCSNe, creating circumstellar medium (CSM) with which the SNe interact. Observations during the first few days post-explosion can provide important information about the mass-loss rate during the late stages of stellar evolution. Model fits to the quasi-bolometric light curve of SN 2020fqv reveal 0.23 M⊙ of CSM confined within 1450 R⊙ (1014 cm) from its progenitor star. Early spectra (<4 d post-explosion), both from HST and ground-based observatories, show emission features from high-ionization metal species from the outer, optically thin part of this CSM. We find that the CSM is consistent with an eruption caused by the injection of ∼5 × 1046 erg into the stellar envelope ∼300 d pre-explosion, potentially from a nuclear burning instability at the onset of oxygen burning. Light-curve fitting, nebular spectroscopy, and pre-explosion HST imaging consistently point to a red supergiant (RSG) progenitor with $M_{\rm ZAMS}\approx 13.5\!-\!15 \, \mathrm{M}_{\odot }$, typical for SN II progenitor stars. This finding demonstrates that a typical RSG, like the progenitor of SN 2020fqv, has a complicated mass-loss history immediately before core collapse.
Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 and 100 Myr can be used to probe the origin of the Solar system. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: 26 Al, 36 Cl, 41 Ca, 53 Mn, 60 Fe, 92 Nb, 97 Tc, 98 Tc, 107 Pd, 126 Sn, 129 I, 135 Cs, 146 Sm, 182 Hf, and 205 Pb. Here, we probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25 M ⊙ , explosion energies between 3.4 × 10 50 and 1.8 × 10 52 erg and compact remnant masses between 1.5 and 4.89 M ⊙ . We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as 92 Nb and 97 Tc, are the most affected by the remnant mass, 92 Nb varying by five orders of magnitude. Isotopes synthesized primarily in explosive C-burning and explosive He-burning, such as 60 Fe, are most affected by explosion energies. 60 Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15 M ⊙ model. The final yield of each examined SLR is used to compare to literature models.
Chemical explosions create blast waves with large overpressure disturbances. It is important to develop a standard blast model based on data to accurately predict acoustic blast-wave amplitudes near detonations and invert for explosion energy from distant observations of blast-wave signals. However, open data from large, controlled chemical explosions with reliable ground truth can be challenging to find. The lack of access to such data could limit the number of contributions to related research and potentially stifle the rate of discoveries or validation of existing models. Here, to address these data scarcity problem, we have curated and compiled a standardized set of 817 blast-wave waveforms from 19 distinct high-explosive events. The blast-wave waveforms are standardized to a 1 kg trinitrotoluene explosion using scaling laws and corrections for location effects. A brief overview of the dataset is presented along with explosion feature models as well as recommendations for extracting explosion features. The resulting dataset is distributed to an open repository in both Seismic Analysis Code and pandas DataFrame formats containing the waveforms, the scaled distances, and the sample rates.
The Deflagration-to-Detonation Transition (DDT) in one-dimensional porous explosive, where combustion in an explosive transitions to detonation, can be described by the following model. This simplified model proceeds in five steps, as follows: 1) Ignition of the explosive, surface burning. 2) Convective burning, with the flame front penetrating through the porous network of the explosive. This proceeds until the pressure grows high enough to result in choked flow in the pores restricting the convective burn. 3) The choked flow results in the formation of a high-density compact of explosive. This compact is driven into undisturbed material by the pressure of the burning explosive. 4) The compression of the undisturbed porous explosive by the compact leads to the ignition of a compressive burn. This builds in pressure until a supported shock forms. 5) The shock builds in pressure until detonation occurs. See Figures 1 and 2 for an overview of the proceeding steps. It has been assumed in the past that the same mechanism which drives DDT in the one dimensional case will apply to both two and three dimensions. However, this has never been experimentally verified and it is possible that an entirely new mechanism drives DDT higher dimensions. In order to test this a series of two dimensional DDT tests were performed to help provide empirical evidence of the mechanism in higher dimensions.
In this paper, we describe studies of shock initiation and shock-to-detonation transitions in energetic materials using a tabletop shock compression microscope with nanosecond time resolution and micrometer spatial resolution. Planar input shocks with durations of 4–20 ns are produced using 0–4.5 km/s laser-launched flyer plates. Emphasis is on measurements of temperature, velocities, pressure, and microstructure using photon Doppler velocimetry (PDV), optical pyrometry and high-speed videography. Techniques are discussed for fabricating disposable shock target arrays of tiny plastic-bonded explosives (PBX), liquid and powder explosives, and single-crystal explosives for high-throughput studies. Optical temperature measurements of shocked triaminotrinitrobenzene (TATB) are discussed. Since TATB is yellow, we developed methods to correct for the blue absorption to obtain more accurate temperatures. Hot spots in shocked polymer-encased HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) crystals are observed in real-time, showing a hot spot produced in a collapsing void that ignites a deflagration. Despite the small dimensions of our explosive charges (typically 1 mm diameter and 250 μm length), we produced reproducible detonation states in solid and liquid explosives using short-duration shocks near the von Neumann spike (VNS) pressure. We show the VNS pressure is associated with a transition to high-efficiency gas production from the explosives. In studies of NM, prior to detonation, we see reaction originating at hot spots which coalesce to form a superdetonation.