Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPLOSIONS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Simulation Study of Near-Surface Coupling of Nuclear Devices vs. Equivalent High-Explosive Charges

A computational study was performed to examine the differences in near-surface ground-waves and air-blast waves generated by high-explosive energy sources and those generated by much higher energy - density low - yield nuclear sources. The study examined the effect of explosive-source emplacement (i.e., height-of-burst, HOB, or depth-of-burial, DOB) over a range from depths of -35m to heights of 20m, for explosions with an explosive yield of 1-kt . The chemical explosive was modeled by a JWL equation-of-state model for a ~14m diameter sphere of ANFO (~1,200,000kg – 1 k t equivalent yield ), and the high-energy-density source was modeled as a one tonne (1000 kg) plasma of ‘Iron-gas’ (utilizing LLNL’s tabular equation-of-state database, LEOS) in a 2m diameter sphere, with a total internal-energy content equivalent to 1 k t . A consistent equivalent-yield coupling-factor approach was developed to compare the behavior of the two sources. The results indicate that the equivalent-yield coupling-factor for air-blasts from 1 k t ANFO explosions varies monotonically and continuously from a nearly perfec t reflected wave off of the ground surface for a HOB ≈ 20m, to a coupling factor of nearly zero at DOB ≈ -25m. The nuclear air - blast coupling curve, on the other hand, remained nearly equal to a perfectly reflected wave all the way down to HOB’s very near zero, and then quickly dropped to a value near zero for explosions with a DOB ≈ -10m. The near - surface ground - wave traveling horizontally out from the explosive source region to distances of 100’s of meters exhibited equivalent - yield coupling - factors t hat varied nearly linearly with HOB/DOB for the simulated ANFO explosive source, going from a value near zero at HOB ≈ 5m to nearly one at DOB ≈ -25m. The nuclear-source generated near-surface ground wave coupling-factor remained near zero for almost all HOB’s greater than zero, and then appeared to vary nearly - linearly with depth-of-burial until it reached a value of one at a DOB between 15m and 20m. These simulations confirm the expected result that the variation of coupling to the ground, or the air, change s much more rapidly with emplacement location for a high-energy-density (i.e., nuclear-like) explosive source than it does for relatively low - energy - density chemical explosive sources. The Energy Partitioning, Energy Coupling (EPEC) platform at LLNL utilizes laser energy from one quad (i.e. 4-laser beams) of the 192 - beam NIF Laser bank to deliver ~10kJ of energy to 1mg of silver in a hohlraum creating an effective small-explosive ‘source’ with an energy density comparable to those in low-yield nuclear devices. Such experiments have the potential to provide direct experimental confirmation of the simulation results obtained in this study, at a physical scale (and time-scale) which is a factor of 1000 smaller than the spatial- or temporal-scales typically encountered when dealing with nuclear explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

LLNL Macroscopic Anisotropic Explosives Research at INL National Security Test Range - Test Results

A select team of 23 engineers, scientists, and explosives specialists from LLNL, LANL, INL, and Marine Raiders from Marine Special Operations Command (MARSOC) and U.S. Special Operations Command (SOCOM) assembled during the second week of November at the INL National Security Test Range near Idaho Falls to investigate and demonstrate fundamental principles of explosives anisotropy. Today's explosives are isotropic in their detonation performance. That is, no matter what direction a detonation runs through bulk explosive, the performance is the same; whereas, anisotropic explosives exhibit different performance, depending on which direction the detonation wave moves through the explosive. The ANISO Team worked in subfreezing temperatures on the Snake River Plain, carrying out 55 experimental explosives shots in four days that lead to a clear understanding of the performance and behavior of an assembly of small, linerless, C4 shaped charges. These shots clearly demonstrated, for the first time, on a macroscopic scale, the principle of anisotropy in measured progression of the detonation through the explosive assembly. The outputs of nine piezo timing pins in the explosive assembly clearly showed detonation progressing through the assembly faster than nominal detonation velocity and moving slower than nominal detonation velocity in the opposite direction. Basic data from these experiments will be used to design and construct explosives assemblies that will be shot in the LLNL High Explosives Applications Facility's (HEAF). These experimental tests will provide refined basic data that will then be used by modelers to develop high explosives models. Computer simulations using these models will then be run to predict performance and design inhomogeneous, anisotropic bulk explosive charges that will be tested at LLNL.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2025

Time-Dependent Satellite Explosion Probabilities for Long-Term Orbital Debris Environment Modeling

On-orbit accidental explosions are a significant contributor to the growth of the orbital debris population. Certain types of satellites have been seen to exhibit different behaviors and timelines of explosions. Recent assessments of on-orbit explosions for three categories of satellites – Russian Sistema Obespecheniya Zapuska (SOZ, Proton 4th stage attitude and ullage motors) units, spacecraft, and rocket bodies – suggest that on-orbit explosions can be modeled using a continuous, time-dependent probability of explosion. This paper discusses the methodology for developing a time-dependent probability of explosion as a function of the satellite’s orbital lifetime based on historical explosions. The SOZ units are found to be best modeled by a modified Gaussian probability distribution, with the peak probability of explosion occurring around 10.4 years on-orbit. Spacecraft and rocket bodies both exhibit exponentially decaying explosion probabilities, with rocket bodies showing fast-, medium-, and slow-time modes of decay. Total cumulative probabilities of explosion are approximately 57% for SOZ units, 4% for spacecraft, and 2% for rocket bodies. Examples of implementation in a simulation of the orbital debris environment using NASA's long-term evolutionary model are given, showing the overall explosion behavior over two centuries is significantly different from that using a previously implemented constant-probability type-dependent explosion rate model and provides a better representation of the historical explosion record.

Alyssa P Manis↗
DOE OSTI2022

Effects of Model Uncertainties in Underground Chemical Explosions on Far-field Results

We used the CTH shock physics code to simulate the explosion of an 18-t chemical explosive at a depth of 250 m. We used the CTH in the two-dimensional axisymmetric (cylindrical) geometry (2DC) and most simulations included fully tamped explosions in wet tuff. Our study focused on parametric studies of three of the traditional strength models available in CTH, namely, geologic-yield, elastic perfectly-plastic von Mises, and Johnson-Cook strength (flow stress) models. We processed CTH results through a code that generates Reduced Displacement Potential (RDP) histories for each simulation. Since RDP is the solution of the linear wave equation in spherical coordinates, it is mainly valid at far-enough distance from the explosion the elastic radius. Among various parameters examined, we found the yield strength to have the greatest effect on the resulting RDP, where the peak RDP reduces almost linearly in log-log space as the yield strength increases. Moreover, an underground chemical explosion results in a cavity whose final diameter is inversely proportional to the material yield strength, i.e., as the material's yield strength increases the resulting final cavity radius decreases. Additionally, we found the choice of explosive material (COMP-C4 versus COMP-B) has minor effects on the peak RDP, where denser COMP-C4 shows higher peak RDP than the less dense COMP-B by a factor of ~1.1. In addition to wet tuff, we studied explosions in dry tuff, salt, and basalt, for a single strength model and yield strength value. We found wet tuff has the highest peak RDP value, followed by dry tuff, salt, and basalt. 2DC simulations of explosions in 11 m radius spherical, hemispherical, and cylindrical cavities showed the RDP signals have much lower magnitude than tamped explosions, where the cavity explosions mimicked nearly decoupled explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Explosion behavior investigation and safety assessment of large-format lithium-ion pouch cells

We report large-format lithium-ion (Li-ion) batteries with high energy density for electric vehicles are prone to thermal runaway (or even explosion) under abusive conditions. In this study, overcharge induced explosion behaviors of large-format Li-ion pouch cells with Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 cathode at different current rates (C-rates) (0.5C, 1C, 2C) were investigated. The explosion characteristics of the cells were elucidated by discussing the evolution of the cell voltage, the surface temperature and the shock wave pressure. Generally, the whole overcharge process could be divided into four stages according to the evolution of several key parameters and the overcharge behaviors; the overcharge C-rate has a great influence on cells’ thermal behaviors. The experimental results showed that the thermal runaway process of Li-ion cells caused by overcharging consisted of two kinds of explosions, physical explosion and chemical explosion. The existence of observable negative pressure zone in the pressure curves indicated that the Li-ion cells are not a self-supplying oxygen system during the explosion. Further, the explosion dynamics parameters were matched. An explosion TNT-equivalent conversion strategy that depended on the pressure of the shock wave was utilized to evaluate the released energy and its hazards. In addition, with respect to the overcharge of Li-ion pouch cells, a safety assessment method and a safety management method were proposed based on the explosion behaviors. From the perspective of battery safety, this study is of great significance for the safety design of Li-ion cells and can provide guidance for engineers to optimize the safety function of battery packs.

25 ENERGY STORAGE↗
DOE OSTI2025

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
DOE OSTI2023

Investigation on the explosion dynamics of large-format lithium-ion pouch cells

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.

42 ENGINEERING↗
DOE OSTI2021

Exploring the Effects of Emplacement Conditions on Explosion P/S Ratios across Local to Regional Distances

High-frequency (~> 2 Hz) seismic P/S amplitude ratios are well-established as a discriminant to distinguish between natural earthquakes and underground explosions at regional distances (~200–1500 km). As research shifts toward identifying lower-yield events, work has begun to investigate the potential of this discriminant for use at local distances (<200 km), in which initial results raise questions about its effectiveness. Here, we utilize data from several chemical explosion experiment series at the Nevada National Security Site in southern Nevada in the United States to study explosion Pg/Lg ratios across the range of local to regional distances. Additionall, the experiments are conducted over differing emplacement conditions, with contrasting geologies and a variety of yields and depths of burial, including surface explosions. We first establish the similarities of Pg/Lg ratios from chemical explosions to those from historic nuclear tests and conclude that, as previous data have suggested, chemical explosion ratios are good proxies for nuclear tests. We then examine Pg/Lg ratios from the new experiment series as functions of distance, yield, depth of burial, and scaled depth of burial (SDOB). At far-local and regional distances, we observe consistently higher ratios from hard-rock explosions compared to ones in a weaker dry alluvium medium, consistent with prior regional distance results. No other trends with yield, depth of burial, or SDOB are strongly evident. Scatter in the observed ratios is very high, particularly at the shortest event-to-station distances, suggesting that small-scale path effects play a significant role. On average, the local distance explosion Pg/Lg ratios show remarkable consistency across all the variations in emplacement. Explosion source models will need to reproduce these results.

58 GEOSCIENCES↗
DOE OSTI2021

Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes

Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10 −13 to 6.3 × 10 −12 m 2 , which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3 −3 × 10 4 bubbles per mm 3 ), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.

Valdivia, Pedro↗
DOE OSTI2018

Categorizing An Explosive Device by its Nuclear Signature (Rev. 1)

This report documents a study to eliminate ambiguity in identifying a very low-yield nuclear detonation in the early stages of the response to such an event. The ambiguity arises from the very similar explosive and radiation effects produced by three categories of explosive devices: An explosively driven radiological dispersal device (RDD) or an accidental explosion in which radiologic material is scattered by chemical explosives. A failed nuclear device in which nuclear material, such as uranium or plutonium, has been scattered by conventional explosives. A nuclear device that, in addition to scattering nuclear material, produced a very low-yield from nuclear fission. The study found that the properties of fission products provide a ready solution to the problem. The radiation dose rate from fission products that were produced in a nearly instantaneous “fission event” decreases extremely rapidly over the first eight hours after that event, whereas nuclear materials and common radiological materials produce an essentially constant dose rate. This contrast provides an unequivocal fission product signature. The measurements required to utilize this signature can be made by first responders with radiation detectors already in the field, on time scales of an hour, although later confirmation using more capable detectors, including spectroscopic capability, is advised. The process of categorizing an explosives-related incident, both pre- and post-detonation, can be represented as a Decision Tree. The Decision Tree illustrates the use of already-fielded radiation detectors to categorize incidents as likely to involve only conventional explosives, an explosively driven RDD, or a nuclear device. The identification of an explosion as a fission event is extremely important to evacuation planners. It has been shown that delayed evacuation after sheltering in place for up to 24 hours after a fission event has the potential to reduce radiation exposure to the affected public and responders by a factor of ten. Traditional evidence collection after a nuclear detonation is likely to be similarly delayed, as the value and urgency of collection efforts must be weighed against radiation and other hazardous exposure of the collectors. The evidence that survives the immediate effects of a nuclear detonation will, in many cases, still be available for collection once radiation levels have diminished to acceptable levels. Certain classes of evidence, notably information stored on digital media, may significantly degrade from exposure to post-detonation radiation fields. Technical nuclear forensics collections will be performed in all events with nuclear yield. Modification of the high-yield collection strategy may be necessary for a very low-yield device.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS2001

Onset of the Magnetic Explosion in Solar Flames and Coronal Mass Ejections

We present observations of the magnetic field configuration and its transformation in six solar eruptive events that show good agreement with the standard bipolar model for eruptive flares. The observations are X-ray images from the Yohkoh soft X-ray telescope (SXT) and magnetograms from Kitt Peak National Solar Observatory, interpreted together with the 1-8 Angstrom X-ray flux observed by Geostationary Operational Environmental Satellites (GOES). The observations yield the following interpretations: (1) Each event is a magnetic explosion that occurs in an initially closed single bipole in which the core field is sheared and twisted in the shape of a sigmoid, having an oppositely curved elbow on each end. The arms of the opposite elbows are sheared past each other so that they overlap and are crossed low above the neutral line in the middle of the bipole. The elbows and arms seen in the SXT images are illuminated strands of the sigmoidal core field, which is a continuum of sheared/twisted field that fills these strands as well as the space between and around them; (2) Although four of the explosions are ejective (appearing to blow open the bipole) and two are confined (appearing to be arrested within the closed bipole), all six begin the same way. In the SXT images, the explosion begins with brightening and expansion of the two elbows together with the appearance of short bright sheared loops low over the neutral line under the crossed arms and, rising up from the crossed arms, long strands connecting the far ends of the elbows; and (3) All six events are single-bipole events in that during the onset and early development of the explosion they show no evidence for reconnection between the exploding bipole and any surrounding magnetic fields. We conclude that in each of our events the magnetic explosion was unleashed by runaway tether-cutting via implosive/explosive reconnection in the middle of the sigmoid, as in the standard model. The similarity of the onsets of the two confined explosions to the onsets of the four ejective explosions and their agreement with the model indicate that runaway reconnection inside a sheared core field can begin whether or not a separate system of overlying fields, or the structure of the bipole itself, allows the explosion to be ejective. Because this internal reconnection apparently begins at the very start of the sigmoid eruption and grows in step with the explosion, we infer that this reconnection is essential for the onset and growth of the magnetic explosion in eruptive flares and coronal mass ejections.

Moore, Ronald L.↗
DOE OSTI2017

Predicting Large-scale Effects During Cookoff of Plastic-Bonded Explosives (PBX 9501 PBX 9502 and LX-14)

In this study, we have made reasonable cookoff predictions of large-scale explosive systems by using pressure-dependent kinetics determined from small-scale experiments. Scale-up is determined by properly accounting for pressure generated from gaseous decomposition products and the volume that these reactive gases occupy, e.g. trapped within the explosive, the system, or vented. The pressure effect on the decomposition rates has been determined for different explosives by using both vented and sealed experiments at low densities. Low-density explosives are usually permeable to decomposition gases and can be used in both vented and sealed configurations to determine pressure-dependent reaction rates. In contrast, explosives that are near the theoretical maximum density (TMD) are not as permeable to decomposition gases, and pressure-dependent kinetics are difficult to determine. Ignition in explosives at high densities can be predicted by using pressure-dependent rates determined from the low-density experiments as long as gas volume changes associated with bulk thermal expansion are also considered. In the current work, cookoff of the plastic-bonded explosives PBX 9501 and PBX 9502 is reviewed and new experimental work on LX-14 is presented. Reactive gases are formed inside these heated explosives causing large internal pressures. The pressure is released differently for each of these explosives. For PBX 9501, permeability is increased and internal pressure is relieved as the nitroplasticizer melts and decomposes. Internal pressure in PBX 9502 is relieved as the material is damaged by cracks and spalling. For LX-14, internal pressure is not relieved until the explosive thermally ignites. The current paper is an extension of work presented at the 26th ICDERS symposium [1].

42 ENGINEERING↗
DOE OSTI2025

Additive manufacturing of high explosives with inert dilution for wave shaping and detonation velocity grading

The performance of a high explosive charge is highly influenced by the detonation velocity and any inherent internal structure. Additively manufactured (AM) high explosive charges allow for a precisely engineered internal structure and the selective placement of discrete or graded volumes of a tailored material, which are unachievable by conventional means. The dilution of a solid explosive by an inert additive for detonation velocity reduction is herein evaluated for a method of material extrusion called a direct-ink-write AM, where an 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane-based ultraviolet -curable energetic-binder paste explosive was diluted by the inert density mock 5-iodo-2′-deoxyuridine for use in the manufacture of multi-material explosive charges. The explosive ink was characterized for printing by rheological and particle size measurements, and the detonation velocities for 0–40 wt. % inert diluent were found to have a continuous linear reduction in detonation velocity down to 14.3% difference. A two-component 3D plane wave explosive lens with 8.038 km/s “fast” and 7.491 km/s “slow” diluted ink in an internal cone demonstrated a reduced time of breakout at the output face from 361 to 53 ns. A linearly graded “fast–slow–fast” cylinder was successfully printed and fired, showing engineered control along the length of the charge as the detonation velocity dropped from 7.851 to 6.700 km/s and then recovered to 7.241 km/s. The extreme control over the internal composition of a high explosive charge via a dual-material extruder/mixer challenges conventional manufacturing approaches and opens new avenues for engineered detonation wavefronts, high explosive charge performance, and new applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗
DOE OSTI2022

Legacy surface change analyses of the 1993 Rock Valley earthquake sequence for direct comparison to planned NA-22 underground conventional high-explosive experiments (Source Physics Experiment 3 (SPE3) - RV/DC Task 1.3 FY22 Final Report)

Recent work under two previous phases of the NNSA NA-22 Source Physics Experiment (SPE) have shown that underground chemical high-explosive experiments can produce detectable surface changes that differ in spatial extent and vertical magnitude depending on the geologic media at the site (Schultz-Fellenz et al., 2018; 2020; Crawford et al., 2021). Neither of these two prior phases of SPE identified natural earthquake-related surface effects in the same region that occurred at a similar depth as the explosive experiments for direct comparison. While earthquakes can also produce surface changes that are detectable using remote sensing data analyses, it is expected that the pattern and spatial extent of surface changes would vary between earthquakes and explosions. However, no direct-observed surface-change signature comparison between earthquakes and explosions has ever been performed. The SPE Phase 3 Rock Valley/Direct Comparison (RV/DC) program presents a unique opportunity to investigate and characterize co-occurrence of both earthquakes and explosions. In this project, we worked to address five tasks in a workflow, as follows: 1. Identify and obtain existing high-resolution legacy satellite and aerial imagery as close in time before and after the 1993 Rock Valley earthquake sequence to temporally constrain the analyses. 2. Transform these pre-earthquake and post-earthquake datasets into digital elevation models (DEMs) using geospatial analysis software packages (e.g., ArcGIS, Agisoft Metashape, and Google Earth Engine). 3. Perform DEM differencing analyses to assess and quantify earthquake-related changes from the 1993 sequence, and develop map products that visualize these analyses. 4. Use the analyses from (3) to: (a) assess spatial distribution and magnitude of surface changes due to the 1993 Rock Valley earthquake sequence, and (b) determine parameters of forthcoming, planned explosion-related surface change data collection from sensors mounted on unmanned aerial vehicles (UAVs) (e.g., spatial extent of collection, design and density of survey control, sensors to deploy, forward speed and line spacing of UAV flight lines, flight altitude). 5. Develop a summary report on the analyses, including how the analyses define parameters and identify focus areas for any future surface change analytical field campaigns related to the explosive experiment. Analyzing these legacy data and identifying whether they can detect any surface changes related to the earthquake sequence facilitates opportunities for direct signature comparison of surface change from explosions at one location, which has never previously been performed. Comparing the surface change signatures from a co-located and depth-equivalent earthquake and an explosion could help to advance remote sensing event discrimination techniques. This report summarizes the work completed toward this ambitious goal.

42 ENGINEERING↗
NASA NTRS2000

Onset of the Magnetic Explosion in Solar Flares and Coronal Mass Ejections

We present observations of the magnetic field configuration and its transformation in six solar eruptive events that show good agreement with the Hirayama-Shibata model for eruptive flares. The observations arc X-ray images from the Yohkoh Soft X-ray Telescope (SXT) and magnetograms from Kitt Peak National Solar Observatory, interpreted together with the 1-8 Angstrom X-ray flux observed by GOES. The observations show: 1. Each event is a magnetic explosion that occurs in an initially closed magnetic bipole in which the core field is sheared and twisted in the shape of a sigmoid, having an oppositely curved elbow on each end. The arms of the opposite elbows are sheared past each other so that they overlap and are crossed low above the neutral line in the middle of the bipole. 2. Although four of the explosions arc ejective (blow open the bipole) and two are confined (arc arrested within the closed bipole), all six begin the same way. In the SXT iniages, the explosion begins with brightening and expansion of the two elbows together with the appearance of short bright sheared loops low over the neutral line under the crossed arms and, rising up from the crossed arms, long strands connecting the far ends of the elbows. 3. All six events arc single-bipole events in that during the onset and early development of the explosion they show no evidence for reconnection between the exploding bipole and any surrounding magnetic fields. We conclude that in each of our events the magnetic explosion was unleashed by runaway tether-cutting via implosive/explosive rcconnection in the middle of the sigmoid, as in the Hirayama-Shibata model. The similarity of the onsets of the two confined explosions to the onsets of the four ejective explosions and their agreement with the model indicate that runaway reconnection inside a sheared core field can begin whether or not the overlying fields (the envelope field and contiguous fields that press down on the envelope) allow the explosion to be ejective. We consider the implications ofthis for the unleashing of magnetic explosions in field configurations composed of to or more impacted bipoles.

Moore, R. L.↗
NASA NTRS2015

An Empirical Non-TNT Approach to Launch Vehicle Explosion Modeling

In an effort to increase crew survivability from catastrophic explosions of Launch Vehicles (LV), a study was conducted to determine the best method for predicting LV explosion environments in the near field. After reviewing such methods as TNT equivalence, Vapor Cloud Explosion (VCE) theory, and Computational Fluid Dynamics (CFD), it was determined that the best approach for this study was to assemble all available empirical data from full scale launch vehicle explosion tests and accidents. Approximately 25 accidents or full-scale tests were found that had some amount of measured blast wave, thermal, or fragment explosion environment characteristics. Blast wave overpressure was found to be much lower in the near field than predicted by most TNT equivalence methods. Additionally, fragments tended to be larger, fewer, and slower than expected if the driving force was from a high explosive type event. In light of these discoveries, a simple model for cryogenic rocket explosions is presented. Predictions from this model encompass all known applicable full scale launch vehicle explosion data. Finally, a brief description of on-going analysis and testing to further refine the launch vehicle explosion environment is discussed.

Blackwood, James M.↗
NASA NTRSDate not supplied

An Interim Set of TNT Curves for LOX/LNG Explosions

Space launch companies are actively developing, or in some cases have already developed, new vehicles that use large quantities of liquid oxygen (LOX) and liquid natural gas (LNG) propellants. This propellant mixture currently lacks sound/verified explosive safety standards for use in explosive siting and flight safety analysis. In the U.S., multiple government and commercial organizations have conducted limited testing, or are preparing to conduct tests, intended to provide a sound technical basis for LOX/LNG explosive standards applicable to space launch vehicle ground and flight operations. The limited testing and analysis conducted so far indicates the potential to produce extremely energetic explosions due to the miscibility of LOX/LNG (methane), a unique feature relative to previously used propellant mixtures such as hydrogen/oxygen or kerosene/oxygen. MOX is a homogenous mixture of methane and oxygen that is possible because of methane’s 100% miscibility in LOX. Preliminary testing by N2L, Inc. described in this paper shows that it is possible to form MOX with a variety of mixing methods, and that MOX is a sensitive high explosive capable of producing overpressures greater than comparable masses of C-4. The limited amount and quality of large-scale LOX/LNG explosive test data, and the potential formation of high-explosive MOX, create significant unknowns in the determination of hazard areas (cleared of the public) for ground operations (such as a static fire test) and flights of launch vehicles with LOX/LNG propellant. This paper will review the methods used by NASA to develop and validate the LOX/LH2 blast model using large-scale explosive test programs such as Project PYRO, the Hydrogen-Oxygen Vertical Impact (HOVI), and the Large-Scale Hydrogen-Oxygen Explosion (LSHOE) tests. This paper will describe how the same process was used to prepare an interim LOX/LNG blast model that applies to various LV failure scenarios and conservatively accounts for potential MOX formation. This paper also summarizes past and future testing and modeling efforts funded by a consortium of NASA, the US Federal Aviation Administration (FAA), and the US Space Force (USSF). These test programs should be concluded within the next 3 years and are intended to provide empirical data for model verification and validation.

Liquid Oxygen↗