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.
On the evening of February 24, 1959, routine operations at the S-Site explosives machining building of Los Alamos Scientific Laboratory were interrupted by the unanticipated explosion of approximately 7.5 pounds of high-energy explosives. Tragically, two men working in Building TA-16-260, Bay 2, were killed instantly: Leo Guerin and Ray Means.
Microwave Interferometry (MI) offers the advantage of a continuous time measurement of detonation front velocity from detonation initiation to disassembly, which is an important step to assure the quality of stockpile high explosives. However, the method is currently characterized by areas of poor signal strength, which lead to low confidence measurements. Experiments in inert materials were conducted to determine if reflective hot spots, pockets of plasma that form during detonation, are responsible due to varying hot spot concentrations. Instead, it was found that the copper tube used in a range of standard HE test configurations is the cause of the poor signal reception. Hot spots were represented by microwave reflective aluminum particles. The aluminum was mixed with Titanium Dioxide, a material electrically similar to the insensitive high explosive, triaminotrinitrobenzene (TATB), in volume percent fractions (VPFs) between 0 and 100% aluminum, in increments of 10%. Reflectivity was measured based on input and reflection received from a test apparatus with a layer representing undetonated explosive and another representing an approaching shockwave. The results showed no correlation between VPF and measured reflectivity test cases while enclosed in the standard copper tube. Upon further testing, each sample’s measured reflectivity independent of the copper enclosure did correlate with VPF. This revealed that the test enclosure currently used for MI measurements is causing poor MI signal reception, and new methods must be developed to account for this aberration in MI measurements.
Strong ground motion induces acoustic waves in the atmosphere that can be detected at great distances. These waves provide a record of acceleration at the epicenter of the subterranean event. While this information is valuable for nuclear monitoring purposes, a systematic study of the variation in acoustic parameters with explosive yield and depth has not yet been conducted. Here, we provide a survey of low frequency sound waves generated during the Source Physics Phase 1 experiment, in which six chemical explosions were detonated in granite. We found that pressure amplitudes increase with explosion size but decrease with depth as expected. Pressure amplitude variability increased with signal magnitude. Surprisingly, peak frequency appears to increase with depth. A possible directional signal was identified for one of the events as well. The results presented here may aid the nuclear monitoring community in developing means of determining event depth and yield using acoustic methods. This will complement existing algorithms based on seismic radiation.
The Large Surface Explosion Coupling Experiment (LSECE) project will include two surface energetic events (chemical explosions) conducted as close as possible to ground zero (GZ) of the Source Physics Experiment (SPE) Phase II Dry Alluvium (DAG) site in Yucca Flat at the Nevada National Security Site (NNSS). Seismic and acoustic data from these events will complement DAG buried source data and improve interpretation of speed-of-sound signals from explosions for national security programs. A satellite image of the region near GZ along with some of the recording stations that we will use for the far- eld simulation is shown in Figure 1. The working design of the HE package is shown in the Figure 2 and is roughly a cube of C-4 bricks.
The goal of this study was to use U.S. nuclear explosions with known source parameters (yield, depth, shot point material and/or parameters) to determine moment-derived yield estimates. This work has been accomplished by performing full moment tensor solutions using regional network data from the LLNL network, and other regional broadband stations. As part of this study, we have calculated solutions for 130 U.S. nuclear explosions and 12 additional chemical explosions. We then take several approaches to using moment tensors to estimate yield, considering both the full and isotropic moment tensors, doing straight regression analysis on the whole dataset, then successively refining the calibration with additional information about material and overburden. We have also tried a completely new approach of using the seismic moment to help estimate the radiated seismic energy and tying this to yield through a seismic efficiency. Results appear to be promising, but more work might be required to make this more useful in an operational sense.
Detonation performance experiments and modeling are reported for the explosive PBX 9701, which is composed of 97% 3,3$\acute{-}$diamino-4,4$\acute{-}$azoxyfurazan (DAAF) and 3% vinylidene uoride-chlorotriuoroethylene copolymer (Kel-F 800) binder by weight. PBX 9701 is a newly-developed reducedsensitivity explosive with increased performance relative to the triaminotrinitrobenzene or TATB-based PBX 9502 while still retaining low sensitivity to mechanical insult. The first detonation performance measurements for this formulation are presented, including front-curvature ratesticks and cylinder expansion test data. Prior shock initiation data is also reviewed. These data are used to develop programmed burn and reactive flow calibrations for existing commonly-used performance models which allow engineering calculations with PBX 9701. The calibration process involves several enhancements relative to conventional approaches including the use of an analytical scaling correlation to speed the calibration process and development of a new methodology to improve the consistency between the programmed burn and reactive model calibration process and the resulting calculations. Overall, PBX 9701 is seen to exhibit improved performance relative to insensitive explosives, with a trinitrotoluene (TNT) equivalence of 1.24. The detonation performance properties are found to be well captured by existing models.
Geophysical observations of an explosion scale with the yield of the explosion. Guided by the principal of yield-scaling and observed relationships in the data, we defined and parameterized seismoacoustic forward models in Ford et al. [2014]. This work extends the model with additional data features of maximum velocity and positive phase duration for the seismic and acoustic data types, respectively, and a new optic data type. Parameters for the seismoacoustic models are estimated via Bayesian inference. The parameterized models can then be used to infer the yield of a new explosion.
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.
High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.
Detonation of explosive devices produces extremely hazardous fragments and hot, luminous fireballs. Prior experimental investigations of these post-detonation environments have primarily considered devices containing hundreds of grams of explosives. While relevant to many applications, such large- scale testing also significantly restricts experimental diagnostics and provides limited data for model validation. As an alternative, the current work proposes experiments and simulations of the fragmentation and fireballs from commercial detonators with less than a gram of high explosive. As demonstrated here, reduced experimental hazards and increased optical access significantly expand the viability of advanced imaging and laser diagnostics. Notable developments include the first known validation of MHz-rate optical fragment tracking and the first ever Coherent Anti-Stokes Raman Scattering (CARS) measures of post-detonation fireball temperatures. While certainly not replacing the need for full-scale verification testing, this work demonstrates new opportunities to accelerate developments of diagnostics and predictive models of post-detonation environments.
Dust is ubiquitous in the interstellar medium, serving as both a source of information and headaches to astronomers, depending on who you ask. The out-flows of core collapse supernovae (CCSNe) serve as a factory for dust nucleation and are one of a few known sources of dust in the universe. From observations of local CCSNe such as 1987A and Cassiopeia A, we know that dust is abundant in CCSN ejecta prior to interacting with the ISM. Furthermore, the explosion energy, explosive engine, metallicity, and progenitor mass of the CCSN will all impact the subsequent dust formation history and composition. Observations of the supernova ejecta probe the detailed composition of the ejecta which, in turn, can be used probe the properties of the progenitor star and the process of the explosion
Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.
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 semi-volatile 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 7 and 8. 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. Nitrate and all CAM17 metals have been added to the sampling plan to comply with DTSC comments received in correspondence dated June 4, 2021. 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.
This report documents the development of the Blue Canyon Dome (BCD) testbed, including test site selection, development, instrumentation, and logistical considerations. The BCD testbed was designed for small-scale explosive tests (~5 kg TNT equivalence maximum) for the purpose of comparing diagnostic signals from different types of explosives, the assumption being that different chemical explosives would generate different signatures on geophysical and other monitoring tools. The BCD testbed is located at the Energetic Materials Research and Testing Center near Socorro, New Mexico. Instrumentation includes an electrical resistivity tomography array, geophones, distributed acoustic sensing, gas samplers, distributed temperature sensing, pressure transducers, and high-speed cameras. This SAND report is a reference for BCD testbed development that can be cited in future publications.
Understanding the geophysical response near an underground explosion is crucial for generating insights into the source and emplacement conditions that produce distinct observations in monitoring scenarios occurring at greater distances. Recently, Shot A of the Low Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1) series was conducted at the Nevada National Security Site to provide ground truth for subsurface explosion signal models. This experiment resulted in measuring near-source ground motion at distances ranging from 70 to 1000 m/kt with a 99% success rate, yielding high-fidelity knowledge of the near-field response that can serve as benchmarks for future numerical modeling and experiment planning. However, technical challenges exist in observing near-source phenomena while safeguarding sensitive data acquisition components from the detrimental effects of ground motion in the subsurface. This report outlines tools and techniques to address challenges associated with observing near-source accelerations and within the tunnel drift of the PE1 test bed. Additionally, we describe key systems designed with both modern advancements and legacy guidance to maximize the collection of high-quality ground motion data, which may be applied to constitutive and computational models, leading to new or improved understanding of the near- and far-field signals produced by underground explosions.
A series of large damaging explosions involving fireworks storage occurred near Esparto, Yolo County California on July 1, 2025. Three explosions were located and reported by the University of California Berkeley Seismology Laboratory (UCB/BSL) and the United States Geological Survey (USGS). Analysis of local distance (< 30 km) seismic recordings of the first blast indicates that there were actually three explosions with later blasts delayed by about 3 and 30 seconds. We measured the first arriving P-wave amplitudes on four of these events with good signal-to-noise ratios.
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.
Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae (SNe) resulting from massive stars with solar metallicity and masses from 9.0 to 120 M{sub ⊙}. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using five standard progenitors, and for the Crab SN using a 9.6 M{sub ⊙} progenitor. As a result of using a calibrated central engine, the final kinetic energy of the SN is variable and sensitive to the structure of each pre-SN star. Many progenitors with extended core structures do not explode, but become black holes (BHs), and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the Sun provided that a reasonable contribution from SNe Ia is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star initial mass function has a mean gravitational mass near 1.4 M{sub ⊙}. The average BH mass is about 9 M{sub ⊙} if only the helium core implodes, and 14 M{sub ⊙} if the entire pre-SN star collapses. Only ∼10% of SNe come from stars over 20 M{sub ⊙}, and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored.