SEARCH · Search NASA

Results for “Detonation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Modelling the Fluctuations of Reactive Shock Waves in Heterogeneous Solid Explosives as Stochastic Processes

Although current phenomenological burn models are useful for describing the average or bulk reactive flow behaviour of heterogeneous explosives, one fundamental weakness inherent to these models is the loss of detailed microstructural information at the scale of the calculation. In order to include the effects of the microstructure, and in particular the underlying material heterogeneities that influence the build-up to detonation, a new paradigm is put forth for modelling sub-grid, reaction-induced fluctuations (i.e. “hot spots”) at the continuum level. This modelling approach assumes that the reaction rate is stochastic, rather than deterministic, and it uses Langevin-type equations with a mathematical framework built upon Itô calculus and Lambourn's CIM model for shocked heterogeneous explosives. This approach follows directly from our previous letter, Ref. [1], and is inspired by the probability density function (pdf) methods used in turbulent reactive flows. Here, the stochastic burn model is derived, implemented, and exercised far beyond what has been shown in previous work. New hydrocode simulation results demonstrate the role of stochastic fluctuations during shock initiation; these fluctuations are approximated by collections of discrete particles, that evolve with drift (i.e. deterministic) and diffusion (i.e. stochastic) coefficients. Furthermore, the particle values are propagated and averaged to calculate the heat release, yield strength, and material impedance in each computational cell. Hydrocode simulation results further show how the fluctuating hot spot energy may or may not be transmitted to the wave front, and result in a detonation wave-like structure. The fundamental stochastic nature of this model permits simulations to have varying outcomes with the same initial conditions; this allows for go/no-go estimation (e.g., marginal or failed detonations), which might possibly be calibrated using the statistical distributions from real materials. Mesoscale calculations of shocked heterogeneous explosives also show that these fluctuations are physically justified (i.e., Ref. [2]), and it is hypothesized that the pdf functions provide a link between the meso(grain) and continuum scales for practical engineering calculations. Hence, our approach is a paradigm shift for building efficient, reduced order continuum burn models, that represent the sub-grid stochastic behaviour of shocked heterogeneous solid explosives.

42 ENGINEERING↗

DOE OSTI2020

Shock Initiation Microscopy with High Time and Space Resolution

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Implications of Underground Nuclear Explosion Cavity Evolution for Radioxenon Isotopic Composition

Abstract Isotopic ratios of radioxenons sampled in the atmosphere or subsurface can be used to verify the occurrence of an underground nuclear explosion (UNE). Differences in the half-lives of radioactive xenon precursors and their decay-chain networks produce different time-dependent concentration profiles of xenon isotopes allowing isotopic ratios to be used for tracking UNE histories including estimating the time of detonation. In this study, we explore the potential effects of post-detonation cavity processes: precipitation of iodine precursors, gas seepage, and prompt venting on radioxenon isotopic evolution which influences UNE histories. Simplified analytical models and closed-form solutions yielding a potentially idealized radioactive decay/ingrowth chain in a closed and well-mixed system typically have limited application by not including the partitioning of the radionuclide inventory between a gas phase and rock melt created by the detonation and by ignoring gas transport from the cavity to host rock or ground surface. In reality, either subsurface transport or prompt release that is principally responsible for gas signatures violates the closed-system (or batch-mode) assumption. A closed-form solution representing time-dependent source-term activities is extended by considering the cavity partitioning process, slow seepage, and/or prompt release of gases from the cavity and applied to realistic systems.

Sun, Yunwei (ORCID:0000000194801439)↗

DOE OSTI2024

Ensuring Σ s Y s = 1 in transport of species mass fractions

When transporting species mass fractions in reacting flow simulations, there are physical constraints that must be met. Unfortunately, nonlinear transport schemes such as weighted essentially non-oscillatory (WENO) schemes do not ensure that the sum of mass fractions equals 1. In detonation simulations, errors in the sum of mass fractions are observed to worsen over time when the standard WENO scheme is used. To prevent this, typically one species is forced to absorb all or most of the error in the sum of mass fractions. Here, this work presents an alternative method for correcting the WENO-interpolated mass fractions such that inert species do not change. The method is demonstrated for both argon and nitrogen-diluted hydrogen detonations in one dimension. Using the modified WENO interpolations, the error in the sum is reduced significantly. In addition, the new approach is better able to capture the physical instability expected for nitrogen-diluted detonations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Dirty bomb source term characterization and downwind dispersion: Review of experimental evidence

Dirty bombs are considered one of the easiest forms of radiological terrorism, a form of terrorism based on the deliberate use of radiological material to cause adverse effects in a target population. One U.S. Government official has even described a dirty bomb attack as “all but inevitable”. While people in the vicinity of the blast may experience acute radiation effects, people downwind may unknowingly be contaminated by the radioactive airborne particulate and face increased long-term cancer risk. The likelihood of increased cancer risk depends on the radionuclide used and its specific activity, its aerosolization potential, the particle sizes generated in the blast, and where a person is with respect to the detonation. Different studies have reported that plausible radionuclides for dirty bomb include 60 Co, 90 Sr, 137 Cs, 192 Ir, 241 Am based on their availability in commercial sources as well as safeguards, the amount needed for adverse health effects, previous mishandling of radionuclides and malicious uses. In order to have increased long-term cancer risk, the radionuclide would have to deposit inside the body by entering the respiratory tract and then possibly migrate to other organs or bones (ground shine is not considered in this paper because areas affected by the event will likely become inaccessible). This implies that the particles will have to be smaller than 10 μm to be inhaled. Experiments involving the detonation of dirty bombs have shown that particles or droplets smaller than 10 μm are generated, independently from the initial radionuclide or its state (e.g., powder, solution). Atmospheric tests have shown that in unobstructed terrain, the radionuclide laden cloud can travel kilometers downwind even for relatively small amounts of explosives. Furthermore, buildings in the path of the cloud can change the dose rate. For instance, in one experiment with a single building, the dose rate was 1–2 orders of magnitude lower behind the obstacle compared to its front face. For people walking around, the amount of particulate deposited on them and inhaled will depend on their path relative to the cloud, resulting in the counterintuitive result that the closer people may actually not be the ones more at risk because they could simply miss the bulk of the cloud in their wandering. In summary, the long-term cancer risk for people caught in a dirty bomb cloud away from the detonation requires considering where and when the people are, which radionuclide was used, and the layout of the obstacles (e.g., buildings, vegetation) in the path of the cloud.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2025

Overcoming Variability: A Reproducible Approach to SERS Detection of Nanodiamonds

Detonation nanodiamonds (DNDs) are formed at specific pressures and temperatures during explosions. Different explosives produce varied yields of DNDs within their detonation soot, with Composition B producing the highest yield. Raman spectroscopy (RS) is often used for the characterization of sp 2 - and sp 3 -hybridized carbon allotropes in carbonaceous materials because of distinct disorder and graphitic bands. Bulk diamond also gives a distinct Raman peak at 1332 cm –1 . Furthermore, as bulk diamond decreases in size to nanometer-sized species, the peak red-shifts and broadens, becoming increasingly difficult to detect with RS using visible excitations. Therefore, surface-enhanced Raman spectroscopy (SERS) was used to enhance the diamond peak of DNDs, enabling better detection and faster examination of DNDs within detonation soot. Previous literature of the SERS of DNDs delivered inconsistent results in spectral signatures and SERS substrates. Herein, refining of the methodology for the acquisition of SERS spectra of DNDs was achieved. Before any SERS experiments, the DNDs were first characterized with normal Raman (NR) and scanning electron microscopy. Two routes for SERS enhancement were evaluated: colloidal noble metal nanoparticles and evaporated silver films. Silver films produced the best signal enhancement of DNDs with the best signal-to-noise and peak enhancements observed at 20–30 nm thick silver films at 5% (∼300 μW) laser power. Consistent, reproducible SERS spectra were acquired of small aggregates of DNDs down to ∼500 nm. NR and SERS mapping analysis of DNDs before and after evaporation of silver films revealed the improvements in the detection capabilities of SERS compared with NR.

Carbon↗

DOE OSTI2022

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

The reactive flow evolution of the polymer-bonded explosive PBX 9502: Experiments and model validation in extreme pressure regimes

The shock-to-detonation transition properties of the triaminotrinitrobenzene based PBX 9502 high explosive (HE) are experimentally and computationally explored in extremely high input pressure conditions. These include both slightly sub-Chapman–Jouguet and overdriven input pressure conditions, namely, ∼25 and ∼31 GPa, respectively. Our experiments capture the transient buildup of a shock-induced reaction via measurement of HE and polymethyl methacrylate window interface particle velocity profiles for a variety of sample thicknesses for this insensitive HE. These observations necessitate extremely thin explosive samples, and the high rates of reaction provide a considerable challenge to optical diagnostics. Samples at these thicknesses also provide an opportunity for evaluation of potential micro-structure effects on the resulting shock-to-detonation-transition measurements. To address this, the thin samples are also characterized via x-ray micro-computed tomography. Finally, a pair of previously established continuum-level detonation performance modeling approaches for PBX 9502 were used to analyze the experiments. The employed model variants crucially differ in their definition of each model’s empirical reaction rate functional form, utilization of shock state quantities, and local flow variable dependencies. As a result, the present experiments provide a novel platform to evaluate the quantitative and qualitative consequences stemming from these modeling choices in a challenging initiation scenario, largely beyond the chosen calibration range of either model. This new experimental information will provide a platform for both improved physics and model parameterizations for this well-studied explosive.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Linear spectropolarimetry of 35 Type Ia supernovae with VLT/FORS: an analysis of the Si$\small {II}$ line polarization

Spectropolarimetry enables us to measure the geometry and chemical structure of the ejecta in supernova explosions, which is fundamental for the understanding of their explosion mechanism(s) and progenitor systems. We collected archival data of 35 Type Ia supernovae (SNe Ia), observed with Focal Reducer and Low-Dispersion Spectrograph (FORS) on the Very Large Telescope at 127 epochs in total. We examined the polarization of the Si ιι λ6355 Å line ($p_{\rm Si\, \small {II}}$) as a function of time, which is seen to peak at a range of various polarization degrees and epochs relative to maximum brightness. We reproduced the $\Delta m_{15}\!-\!p_{\rm Si\, \small {II}}$ relationship identified in a previous study, and show that subluminous and transitional objects display polarization values below the $\Delta m_{15}\!-\!p_{\rm Si\, \small {II}}$ relationship for normal SNe Ia. We found a statistically significant linear relationship between the polarization of the Si ii λ6355 Å line before maximum brightness and the Si ιι line velocity and suggest that this, along with the $\Delta m_{15}\!-\!p_{\rm Si\, \small {II}}$ relationship, may be explained in the context of a delayed-detonation model. In contrast, we compared our observations to numerical predictions in the $\Delta m_{15}\!-\!v_{\rm Si\, \small {II}}$ plane and found a dichotomy in the polarization properties between Chandrasekhar and sub-Chandrasekhar mass explosions, which supports the possibility of two distinct explosion mechanisms. A subsample of SNe displays evolution of loops in the q–u plane that suggests a more complex Si structure with depth. This insight, which could not be gleaned from total flux spectra, presents a new constraint on explosion models. To conclude, we compared our statistical sample of the Si ιι polarization to quantitative predictions of the polarization levels for the double-detonation, delayed-detonation, and violent-merger models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Investigating Fracture Network Deformation Using Noble Gas Release

We investigate deformation mechanics of fracture networks in unsaturated fractured rocks from subsurface conventional detonation using dynamic noble gas measurements and changes in air permeability. We dynamically measured the noble gas isotopic composition and helium exhalation of downhole gas before and after a large subsurface conventional detonation. These noble gas measurements were combined with measurements of the subsurface permeability field from 64 discrete sampling intervals before and after the detonation and subsurface mapping of fractures in borehole walls before well completion. We saw no observable increase in radiogenic noble gas release from either an isotopic composition or a helium exhalation point of view. Large increases in permeability were observed in 13 of 64 discrete sampling intervals. Of the sampling intervals which saw large increases in flow, only two locations did not have preexisting fractures mapped at the site. Given the lack of noble gas release and a clear increase in permeability, we infer that most of the strain accommodation of the fractured media occurred along previously existing fractures, rather than the creation of new fractures, even for a high strain rate event. These results have significant implications for how we conceptualize the deformation of rocks with fracture networks above the percolation threshold, with application to a variety of geologic and geological engineering problems.

Gardner, W. Payton↗

DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

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.

Geochemistry & Geophysics↗

DOE OSTI2020

Project 57 Air Monitoring Report: (Jan 1 - Dec 31 2018)

During the late 1950s, the Atomic Energy Commission (AEC) (now the U.S. Department of Energy [DOE]) conducted a series of safety experiments to determine if a nuclear device subjected to a large conventional explosives detonation would result in a nuclear yield. The AEC obtained temporary use of a large portion of western Emigrant Valley from the U.S. Department of Defense (DOD) for one of these experiments: Project 57. Following the Project 57 safety experiment, the AEC fenced the contaminated area and returned control of the surrounding land to the DOD. Emigrant Valley is part of the Nevada Test and Training Range (NTTR). For safety and security reasons, access to the NTTR is controlled using both physical (i.e., fences) and administrative (e.g., signs and postings) controls. Therefore, the public cannot access the Project 57 site and there are no known human receptors that routinely access the site. Project 57 was detonated on April 24, 1957, in Emigrant Valley approximately 13 mi (21 km) northeast of the north end of Yucca Flat. This test was undertaken to develop: (1) a means of estimating immediate distribution and long-term redistribution of plutonium dispersed during a nonnuclear detonation; (2) biomedical evaluation techniques for use in likely plutonium-laden environments; (3) methods of decontaminating ground areas, pavements, and building materials; and (4) alpha survey instruments and field monitoring procedures to promptly estimate contaminant deposition. Although the test did not result in the fission of nuclear materials, it did disseminate plutonium across the ground surface. Various radiological surveys have been performed in the area since Project 57 was conducted. The AEC delineated the area containing radioactive material based on radioactivity surveys performed shortly after the test was conducted. A designated Contamination Area (CA), marked by a fence, was established later. The distribution of americium-241 (Am-241) in the area was determined again in a 1997 flyover, which showed Am-241 ranging from as much as 70,000 counts per second (cps) at ground zero to background (<70 cps). This survey also documented Am-241 radioactivity on the ground surface beyond the east side CA fence at levels of up to 150 cps. In 2007, the DOE expanded the CA by posting “Contamination Area” signs 200 ft to 400 ft (60 m to 120 m) beyond the original fence, which formed a new CA boundary.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Near-Field Imaging of Shallow Chemical Explosions in Granite Using Change Detection Methods with Surface and Borehole Seismic Data

Explosions detonated in geologic media damage it in various ways via processes that include vaporization, fracturing, crushing of interstitial pores, etc. Seismic waves interact with the altered media in ways that could be important to the discrimination, characterization, and location of the explosions. As part of the Source Physics Experiment, we acquired multiple pre- and post-explosion near-field seismic datasets and analyzed changes to seismic P-wave velocity. Our results indicate that the first explosion detonated in an intact media can cause fracturing and, consequently, a decrease in P-wave velocity. After the first explosion, subsequent detonations in the pre-damaged media have limited discernible effects. We hypothesize this is due to the stress-relief provided by a now pre-existing network of fractures into which gasses produced by the explosion migrate. We also see an overall increase in velocity of the damaged region over time, either due to a slow healing process or closing of the fractures by subsequent explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2020

Soil Sample Plan in Support of the Site 300 Explosives Waste Treatment Facility Ecological Risk Assessment Pursuant to the Hazardous Waste Facility Permit, Part V, Special Condition #14 a-f. August 2020

The Department of Toxic Substances Control issued the Final Hazardous Waste Facility Permit to Lawrence Livermore National Laboratory, Site 300 in 2017. As part of the continual process to ensure waste treatment emissions are not impacting ecological receptors, LLNL analyzed soil samples in the vicinity of the Explosive Waste Treatment Facility (EWTF) Burn Units (i.e., Burn Cage and Burn Pan) and the Detonation Pad. Soil samples were analyzed by EPA Solid Waste 846 Methods for furans, explosives, metals, and semivolatile compounds. In all samples, furans, explosives, and semi-volatile organic compounds were not detected. Metals were detected at concentrations below background levels. This soil sample plan will sample the same four downwind areas of the Burn Units and the same downwind areas of the Detonation, as well as including four new sample locations. The four new locations will be located to the northeast and southwest of the Burn Units and Detonation (as shown in Figure 3). The constituents to be analyzed are in Tables 2 and 3. Perchlorate was not identified as a constituent of concern in the 2007 Soil Sampling Plan. However, for this plan, all soil samples will be analyzed for perchlorate. After completion of the sampling event, a soil sample report will be submitted to DTSC in accordance with requirements specified in the Hazardous Waste Facility Permit, Part V, Special Conditions 16 a-f. This proposed soil sampling plan will be fully implemented no later than one year after receiving DTSC approval, in accordance with the Hazardous Waste Facility Permit, Part V, Special Condition 15.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Scaling Law for Cloud-Rise Velocity vs. Scaled Height of Burst (SHOB)

The shockwave produced by a nuclear detonation leaves a low-density region in the vicinity of air zero (i.e., the point of detonation). As the shockwave degenerates into an acoustic wave, the buoyant force takes over and causes the low-density region to begin to rise at speeds that can reach as high as ~200-300 mph. The objective of this project was to find out if this initial cloud rise followed any scaling law, by stimulating a cloud rise for a 1 kiloton (kt) detonation and determining if a correlation existed between the scaled-height-of-burst (SHOB) and the scaled-cloud-rise velocity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Verification and Validation of High Explosive Reactive Burn Models Implemented in LANL's EAP and LAP Code Base

Reactive burn models represent a significant leap in high explosive (HE) modeling capability. The first generation of engineering models of HE detonation are called programmed burn models and they are largely based on the distance between a prescribed detonation point and each zone in a simulation. There have been many advancements to programmed burn models over the years and when the assumptions upon which they are based are met, a properly tuned programmed burn model can be highly accurate but if any of their assumptions is not met, as is the case for corner turning or weakly initiated HE burn, they will give the wrong answer. Reactive burn models represent an entirely new way of modeling HE burn. They use the local conditions of a zone – e.g. temperature, pressure or density – as calculated by a hydrocode to determine if and when the zone is going to detonate and if so, how rapidly. This difference opens up an entirely new set of capabilities for HE modeling. It makes it possible to accurately and predictively model phenomena like the effect of confinement and the formation of dead zones. Reactive burn models have seen sustained development effort at LANL for at least the last decade but several recent developments make it timely to transition reactive burn models from a research topic to a production tool. The main goal of this milestone is to facilitate and accelerate the adoption of reactive burn as a commonly available modeling option, with recommendations on the resolution that will be required and uncertainties associated with their modeling choices. To achieve this, we have performed verification, validation, and uncertainty quantification (UQ) assessments of AWSD and SURF/SURFplus in xRage and FLAG on a variety of different problems.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2021

Shocked Transport: Experiments to Study Dynamic Particle Behavior under Varying Volume Fraction Conditions

Shocks interacting with densely-packed porous media is common in volcanic eruptions or in engineering applications such as ordinance. Large-scale, exploratory experiments have been performed by Frost and collaborators examining the behavior of a bed of particles subjected to a detonation wave. The detonation replaces the central explosive with a high-density, high-temperature gas that propagates outward. A shock wave travels through the particle bed, followed by a contact discontinuity between the shocked air and the explosive products. In the early times, the detonation wave causes many of the particles to collide and compact with each other. As the particles disperse, the experiments demonstrate a series of large-scale jets or fingers that form in the advancing particle front.

42 ENGINEERING↗