Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explosives”

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 217 records · Page 12

Rapid Analyses of Sparse Seismoacoustic Data Reveals the Timing and Size of the Accurate Energetic Systems Explosion

On 10 October 2025 an explosion occurred at a facility operated by Accurate Energetic Systems in Humphreys County, Tennessee. The incident resulted in 16 fatalities and created a debris field over several square kilometers. To address remaining questions about explosion timing and size, we collected about 20 seismic and 19 acoustic records of the blast from sensors up to hundreds of kilometers away. We then deployed 10 distinct physics-based, reduced order models (ROMs) that used validated geological structure and atmospheric conditions from the time of the event, along with observations of body- and surface-wave energy, as well as acoustic overpressure and phase duration. Each ROM predicted either timing, yield estimates, or both. We binned these estimates and their uncertainties according to each ROMs’ assumptions about confinement (aboveground, buried fully coupled, and buried partially coupled) and combined these estimates with other forensic data to conclude that the event occurred as a single, aboveground explosion on 10 December 2025 12:47:50.8 ±0.1 s with a yield equivalent to 11.8 [2.3,16.5] tons of Trinitrotoluene. Our estimates align with the Bureau of Alcohol, Tobacco, Firearms and Explosives inventory reports of 11–13 tons. This multimethod approach demonstrates the use of remotely observed geophysical data to rapidly aid conventional forensic investigations of accidental explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Detecting Surface Changes from an Underground Explosion in Granite Using Unmanned Aerial System Photogrammetry

Efficient detection and high-fidelity quantification of surface changes resulting from underground activities are important national and global security efforts. In this investigation, a team performed field-based topographic characterization by gathering high-quality photographs at very low altitudes from an unmanned aerial system (UAS)-borne camera platform. The data collection occurred shortly before and after a controlled underground chemical explosion as part of the United States Department of Energy’s Source Physics Experiments (SPE-5) series. The high-resolution overlapping photographs were used to create 3D photogrammetric models of the site, which then served to map changes in the landscape down to 1-cm-scale. Separate models were created for two areas, herein referred to as the test table grid region and the nearfield grid region. The test table grid includes the region within ~40 m from surface ground zero, with photographs collected at a flight altitude of 8.5 m above ground level (AGL). The near-field grid area covered a broader area, 90–130 m from surface ground zero, and collected at a flight altitude of 22 m AGL. The photographs, processed using Agisoft Photoscan® in conjunction with 125 surveyed ground control point targets, yielded a 6-mm pixel-size digital elevation model (DEM) for the test table grid region. This provided the ≤3 cm resolution in the topographic data to map in fine detail a suite of features related to the underground explosion: uplift, subsidence, surface fractures, and morphological change detection. The near-field grid region data collection resulted in a 2-cm pixel-size DEM, enabling mapping of a broader range of features related to the explosion, including: uplift and subsidence, rock fall, and slope sloughing. This study represents one of the first works to constrain, both temporally and spatially, explosion-related surface damage using a UAS photogrammetric platform; these data will help to advance the science of underground explosion detection.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Potential of Argon Detection Capabilities for Nuclear Explosion Monitoring

Abstract Current noble gas detection systems for nuclear explosion monitoring are based on the detection of four radioxenon isotopes—Xe-131m, -133, -133m and -135. The data provided by radioxenon detection could be enhanced by other radionuclide signatures such as Ar-37. Activation of Ca-40 in rock by neutrons produces Ar-37, and monitoring for this additional nuclide could help distinguish detections of nuclear explosions from background sources, such as medical isotope production. This work studies the capabilities of a hypothetical argon detection network. A 10 kt explosion was modeled using MCNP and SCALE to determine the inventory of Ar-37 created in a representative granite rock layer, assuming either 0.1, 1 or 10% of the total inventory was released. The Ar-37 inventory was combined with atmospheric transport data from HYSPLIT compiled in a previous study, along with the detection limits of standard Ar-37 detection systems, to determine how many hypothetical monitoring stations would detect Ar-37 from an explosion. This method was repeated for 365 HYSPLIT data sets to create a year’s worth of hypothetical explosions, releases, and detections. The study quantified the average number of detections per release, the number of stations detecting Ar-37, and the possibility of detecting Ar-37 in coincidence with xenon.

37Ar↗

A comparison of explosively driven shock wave radius versus time scaling approaches

Abstract Explosively driven shock wave radius versus time profiles are frequently used to document energy release and relative explosive performance. Recently, two universal shock wave radius versus time profiles have been presented in the literature, which demonstrate the ability to represent explosively driven shock wave profiles for all explosive sources in any fluid environment. These two universal shock wave profiles are examined here relative to each other and relative to a commonly used nonlinear shock wave profile, which is fit to experimental data for individual explosive materials. The nonlinear profile, originally developed by Dewey, is examined here, and a universal non-dimensional form of the equation is proposed. The universal shock wave profiles are all found to be relatively similar, but with slight variations in a transition region of non-dimensional radii $$0.15\lesssim R^*\lesssim 2$$ 0.15 ≲ R ∗ ≲ 2 . The variations in this region result in different estimations of energy release or blast strength between the curve fits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isolation and characterization of cellulosic fibers from kenaf bast using steam explosion and Fenton oxidation treatment

Degumming is the dominant method to isolate cellulosic fibers in the textile industry. High content of lignin is the principle obstacle in bast-fibers degumming. In this study, in order to remove lignin efficiently, the treatment using steam explosion coupled with Fenton oxidation was conducted on kenaf bast. The influence and mechanism of the combined treatment of steam explosion and Fenton oxidation were studied. Fiber microstructure, composition proportion and distribution, cellulose Segal crystallinity and polymerization as well as residual lignin structure were analyzed. The results showed that the pretreatment of steam explosion could sufficiently extract most of the hemicellulose from kenaf bast. Also, the Fenton oxidation followed by alkaline boiling could effectively remove the lignin from the fiber, especially from the fiber surface. The refined dry fibers produced only remained 8.82% of hemicellulose and 9.36% of lignin. And the whole treatment process didn’t change the cellulose crystalline structure but could gradually degrade hemicellulose, lignin and other amorphous gummy matter as well as disordered areas in cellulose, leading to the increase of fiber Segal crystallinity from 62.02 to 70.12%. Part of cellulose was damaged with the increase of degumming intensity, resulting in the decrease of the polymerization degree of cellulose. Meanwhile, the structure of residual lignin changed during the degumming stages of steam explosion, Fenton oxidation and alkaline boiling. In conclusion, all the results revealed that the treatment of combined steam explosion and Fenton oxidation is significantly effectual for isolating bast-fibers.

36 MATERIALS SCIENCE↗

Explosive Byproduct Gas Transport Through Sorptive Geomedia

Current underground nuclear explosion (UNE) detection strategies rely heavily on atmospheric noble gas sampling of radioxenon. However, discriminating nuclear weapons testing programs from civilian sources is difficult due to highly variable atmospheric radioxenon backgrounds and processes affecting subsurface transport of parent radionuclides. Here, we aim to study the transport of gases produced by subsurface explosions as novel stable signatures for underground nuclear explosion (UNE) monitoring. These gases may be produced in large quantities with distinct molecular ratios, which will be impacted by subsurface transport processes. To demonstrate how ratios of gases produced by explosions can change during transport in geomaterials, we conducted laboratory benchtop experiments on the transport of carbon dioxide (CO 2 ) and hydrogen (H 2 ) gases through variably saturated zeolitic tuff, which is abundant at the historic US testing site. We observed that zeolitic tuff sorbs substantial quantities of CO 2 while allowing H 2 to transport more freely, leading to changes in the molecular ratios of the two gases along the transport pathway. Gas uptake in the dry zeolitic tuff core was 72.3% for CO 2 , compared with 53.4% for xenon and 7.6% for H 2 . The presence of 20% water saturation disrupted the CO 2 sorption process, though to a lesser extent than observed for noble gases, with a 36.7% drop in xenon sorption compared with a 21.9% drop for CO 2 . These results represent the first observations of zeolite sorption altering explosive gas ratios during transport through geomedia relevant to nuclear proliferation monitoring.

54 ENVIRONMENTAL SCIENCES↗

Relating microstructure, temperature, and chemistry to explosive ignition and shock sensitivity

An analysis is put forth that relates explosive properties including microstructure, temperature, and chemistry to explosive ignition and sensitivity. Unlike approaches that focus only on elucidating mesoscopic mechanisms important to ignition, the present analysis seeks a methodology directly applicable to continuum explosive burn models. Because the Scaled Uniform Reactive Front (SURF) burn model has a microstructural basis, it is chosen as the starting point of the present analysis. To build upon the SURF framework, a literal translation is performed of the high-level conceptual notions for which SURF is based, to concrete ignition–combustion parameters, a statistical description of the microstructure, and other thermo-chemical data acquired by non-shock experiments. The analysis requires a void volume fraction distribution acquired from ultra-small angle neutron scattering (USANS). The shock response of PBX 9502 is used to illustrate the theory. While the analysis requires calibration to a shock experiment (the pop plot, for example), the results are a self-consistent set of realistic physical quantities that contribute to the shock initiation process including, initial temperature, chemical kinetics, a statistical description of the microstructure, and the hot spot size, spacing, and temperature. The utility of having a mesoscale based theory is that once the critical hot spot temperature, as a function of shock pressure, is known for a given explosive and initial porosity distribution, the entire set of meso- and macro-scale results, including the pop plot can be calculated for any other porosity. Thus, one can understand changes in sensitivity at different densities (void size distributions), relying only on the assumption that the hot spot temperature curve would not differ significantly if the morphology of the hots spots were similar. Another important utility of the present analysis is to address the question of the role of initial temperature on the observed shock sensitivity of PBX 9502 and other explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗

A digital seismogram archive of nuclear explosion signals, recorded at the Borovoye Geophysical Observatory, Kazakhstan, from 1966 to 1996

Seismologists from Kazakhstan, Russia, and the United States have rescued the Soviet-era archive of nuclear explosion seismograms recorded at Borovoye in northern Kazakhstan during the period 1966–1996. The signals had been stored on about 8000 magnetic tapes, which were held at the recording observatory. After hundreds of man-years of work, these digital waveforms together with significant metadata are now available via the project URL, namely http://www.ldeo.columbia.edu/res/pi/Monitoring/Data/ as a modern open database, of use to diverse communities. Three different sets of recording systems were operated at Borovoye, each using several different seismometers and different gain levels. For some explosions, more than twenty different channels of data are available. A first data release, in 2001, contained numerous glitches and lacked many instrument responses, but could still be used for measuring accurate arrival times and for comparison of the strengths of different types of seismic waves. The project URL also links to our second major data release, for nuclear explosions in Eurasia recorded in Borovoye, in which the data have been deglitched, all instrument responses have been included, and recording systems are described in detail. This second dataset consists of more than 3700 waveforms (digital seismograms) from almost 500 nuclear explosions in Eurasia, many of them recorded at regional distances. It is important as a training set for the development and evaluation of seismological methods of discriminating between earthquakes and underground explosions, and can be used for assessment of three-dimensional models of the Earth’s interior structure.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Correlating Mechanical Sensitivity with Spin Transition in the Explosive Spin Crossover Complex [Fe(Htrz) 3 ] n [ClO 4 ] 2 n

Spin crossover complexes are known to undergo bond length, volume, and enthalpy changes during spin transition. In an explosive spin crossover complex, these changes could affect the mechanical and initiation sensitivity of the explosive and lead to the development of a new class of sensitivity switchable materials. To explore this relationship, the well-known spin crossover compound [Fe(Htrz) 3 ] n [ClO 4 ] 2n (1) was re-evaluated for its explosive properties, and its mechanical impact sensitivity was correlated to spin transition. A variable temperature impact test was developed and used to evaluate the impact sensitivity of 1 in the low spin (LS, S = 0), thermally accessed high spin (HS, S = 2), and mixed LS and HS states. For comparison, the structurally similar Ni compound, [Ni(Htrz) 3 ] n [ClO 4 ] 2n (2), which does not undergo a spin transition at accessible temperatures, was synthesized and characterized, and its explosive properties and variable temperature impact sensitivity measured. These results reveal a correlation between impact sensitivity and spin transition, where 1 exhibits lower impact sensitivity in the LS state and increases in sensitivity upon transition to the HS state. Here, density functional theory was used to predict structural changes that occur upon spin transition that correlate to the change in sensitivity. This demonstrates, for the first time, an explosive spin crossover compound (ExSCO) that exhibits switchable impact sensitivity with a fully reversible internal switching mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization of high-explosive detonations using broadband infrared external cavity quantum cascade laser absorption spectroscopy

Infrared laser absorption spectroscopy provides a powerful tool for probing physical and chemical properties of high explosive detonations. A broadly-tunable swept-wavelength external cavity quantum cascade laser (ECQCL) operating in the mid-wave infrared (MWIR) spectral region is used to measure transmission through explosive fireballs generated from 14 g charges of 4 different explosive types detonated in an enclosed chamber. Analysis of time-resolved transmission and emission at a 2 µs sampling rate shows the evolution of fireball infrared opacity in the first 10 ms after detonation. Broadband high-resolution absorption spectra acquired over the spectral range 2050-2300 cm -1 (4.35-4.88 µm) at a 100 Hz rate are used to measure properties of fireball evolution over longer time scales out to 100 s. Path-integrated concentrations of combustion products CO, CO 2 , H 2 O, and N 2 O are measured and show evolutions over multiple time scales and significant differences between explosive types. Spectral analysis is used to characterize gas temperature and to measure broadband attenuation from absorption and scattering of particulates. In conclusion, analysis of the results provides information on the MWIR optical properties, gaseous detonation/combustion products, and particulates throughout the explosive process including initial detonation, fireball expansion and cooling, and diffusive mixing in the chamber.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Prediction of impact sensitivity, heat of formation and heat of explosion using atomic connectivity

In these proceedings we revisit a large collection of explosives and explosive descriptors with the goal of predicting impact sensitivity using only local atomic environments that can be deciphered from molecular SMILES strings as descriptors without utilizing empirically measured values or computationally expensive electronic structure calculations. From the original database of nearly 500 descriptors, removing empirically measured and electronic structure values decreased the number of descriptors to 135, which we reduced to 18 the most important descriptors using Random Forests. The condensed model predicted impact sensitivity with essentially the same accuracy as the existing, more complex model (R 2 = 0.788 and RSME = 0.312), while remaining applicable to all types of explosives (Peroxides, azides, C-Nitros, Nitroamines, Nitrate Esters, etc.). In addition to impact sensitivity, we proposed similar models to accurately predict values heat of formation (ΔH f ) and heat of explosion (Q), with R 2 = 0.966 and 0.916, respectively. In conclusion, the work in these proceedings allows for prediction of explosive performance and sensitivity with only chemical structure information and an estimate of density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The effect of hardness on polymer-bonded pentaerythritol tetranitrate (PETN) explosive impact sensitivity

Handling sensitivity is an important property to assess when working with explosive samples and can be measured using a variety of tests, including drop-weight impact sensitivity. There exists a longstanding interest in the explosives community on the importance of measurable chemical, physical, and mechanical properties of explosives in impact sensitivity. However, most recent work in this area has explored chemical attributes rather than physical and mechanical properties of explosives. In this study, we measure hardness of explosive samples of pentaerythritol tetranitrate and Sylgard binder (XTX) during the curing process. The samples have been characterized for particle morphology through scanning electron microscopy and handling sensitivity through drop-weight impact testing. The relative importance of states of cure, methods of curing, morphology, and age of material are discussed. The data indicate that although there is a notable difference in morphology and mechanical properties for the samples as the polymer-bonded mixtures cure, the resulting changes to mechanical properties have a minimal effect on the sensitivity of the XTX.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Precursors prior to type IIn supernova explosions are common: Precursor rates, properties, and correlations

There is a growing number of Type IIn supernovae (SNe) which present an outburst prior to their presumably final explosion. These precursors may affect the SN display, and are likely related to poorly charted phenomena in the final stages of stellar evolution. By coadding Palomar Transient Factory (PTF) images taken prior to the explosion, here we present a search for precursors in a sample of 16 Type IIn SNe. We find five SNe IIn that likely have at least one possible precursor event (PTF 10bjb, SN 2010mc, PTF 10weh, SN 2011ht, and PTF 12cxj), three of which are reported here for the first time. For each SN we calculate the control time. We find that precursor events among SNe IIn are common: at the one-sided 99% confidence level, >50% of SNe IIn have at least one pre-explosion outburst that is brighter than 3 × 10{sup 7} L{sub ☉} taking place up to 1/3 yr prior to the SN explosion. The average rate of such precursor events during the year prior to the SN explosion is likely ≳ 1 yr{sup –1}, and fainter precursors are possibly even more common. Ignoring the two weakest precursors in our sample, the precursors rate we find is still on the order of one per year. We also find possible correlations between the integrated luminosity of the precursor and the SN total radiated energy, peak luminosity, and rise time. These correlations are expected if the precursors are mass-ejection events, and the early-time light curve of these SNe is powered by interaction of the SN shock and ejecta with optically thick circumstellar material.

79 ASTRONOMY AND ASTROPHYSICS↗

THE HOST GALAXY OF THE SUPER-LUMINOUS SN 2010gx AND LIMITS ON EXPLOSIVE {sup 56}Ni PRODUCTION

Super-luminous supernovae have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. While these extremely luminous explosions have been observed from z = 0.1 to 1.55, the closest explosions allow more detailed investigations of their host galaxies. We present a detailed analysis of the host galaxy of SN 2010gx (z = 0.23), one of the best studied super-luminous type Ic supernovae. The host is a dwarf galaxy (M{sub g} = -17.42 {+-} 0.17) with a high specific star formation rate. It has a remarkably low metallicity of 12 + log (O/H) = 7.5 {+-} 0.1 dex as determined from the detection of the [O III] {lambda}4363 line. This is the first reliable metallicity determination of a super-luminous stripped-envelope supernova host. We collected deep multi-epoch imaging with Gemini + GMOS between 240 and 560 days after explosion to search for any sign of radioactive {sup 56}Ni, which might provide further insights on the explosion mechanism and the progenitor's nature. We reach griz magnitudes of m{sub AB} {approx} 26, but do not detect SN 2010gx at these epochs. The limit implies that any {sup 56}Ni production was similar to or below that of SN 1998bw (a luminous type Ic SN that produced around 0.4 M{sub Sun} of {sup 56}Ni). The low volumetric rates of these supernovae ({approx}10{sup -4} of the core-collapse population) could be qualitatively matched if the explosion mechanism requires a combination of low-metallicity (below 0.2 Z{sub Sun }), high progenitor mass (>60 M{sub Sun }) and high rotation rate (fastest 10% of rotators).

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity of Seismically Cued Antineutrino Detectors to Nuclear Explosions

Here, we evaluate the sensitivity of large gadolinium-doped water detectors to antineutrinos released by nuclear-fission explosions, using updated signal and background models and taking advantage of the capacity for seismic observations to provide an analysis trigger. Under certain realistic conditions, the antineutrino signature of a 250-kton pure-fission explosion could be identified several hundred kilometers away, in a detector about the size of the largest module currently proposed for a basic physics experiment. In principle, such an observation could provide rapid confirmation that the seismic signal coincided with a fission event, possibly useful for international monitoring of nuclear-weapon tests. We discuss the limited potential for seismically cued antineutrino observations to constrain fission yield, differentiate pure fission from fusion-enhanced weapon tests, indicate that the seismic evidence of an explosion had been intentionally masked, or verify the absence of explosive testing in a targeted area. We conclude that advances in seismic monitoring and neutrino physics have made the detection of explosion-derived antineutrinos more conceivable than previously asserted, but the size and cost of sufficiently sensitive detectors continue to limit applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Magnetorotational Mechanism of the Explosion of Core-Collapse Supernovae

The idea of the magnetorotational explosion mechanism is that the energy of rotation of the neutron star formed in the course of a collapse is transformed into the energy of an expanding shock wave by means of a magnetic field. In the two-dimensional case, the time of this transformation depends weakly on the initial strength of the poloidal magnetic field because of the development of a magnetorotational instability. Differential rotation leads to the twisting and growth of the toroidal magnetic-field component, which becomes much stronger than the poloidal component. As a result, the development of the instability and an exponential growth of all field components occur. The explosion topology depends on the structure of the magnetic field. In the case where the initial configuration of the magnetic field is close to a dipole configuration, the ejection of matter has a jet character, whereas, in the case of a quadrupole configuration, there arises an equatorial ejection. In either case, the energy release is sufficient for explaining the observed average energy of supernova explosion. Neutrinos are emitted as the collapse and the formation of a rapidly rotating neutron star proceeds. In addition, neutrino radiation arises in the process of magnetorotational explosion owing to additional rotational-energy losses. If the mass of a newborn neutron star exceeds the mass limit for a nonrotating neutron star, then subsequent gradual energy losses may later lead to the formation of a black hole. In that case, the energy carried away by a repeated flash of neutrino radiation increases substantially. In order to explain an interval of 4.5 hours between the two observed neutrino signals from SN 1987A, it is necessary to assume a weakening of the magnetorotional instability and a small initial magnetic field (10{sup 9}−10{sup 10} G) in the newly formed rotating neutron star. The existence of a black hole in the SN 1987A remnant could explain the absence of any visible pointlike source at the center of the explosion.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaling Turbulent Combustion Fields in Explosions

We considered the topic of explosions from spherical high-explosive (HE) charges. We studied how the turbulent combustion fields scale. On the basis of theories of dimensional analysis by Bridgman and similarity theories of Sedov and Barenblatt, we found that all fields scaled with the explosion length scale r0. This included the blast wave, the mean and root mean squared (RMS) profiles of thermodynamic variables, combustion variables, velocities, vorticity, and turbulent Reynolds stresses. This was a consequence of the formulation of the problem and our numerical method, which both satisfied the similarity conditions of Sedov. We performed numerical simulations of 1 g charges and 1 kg charges; the solutions were identical (within roundoff error) when plotted in scaled variables. We also explored scaling laws related to three-phase pyrotechnic explosions. We show that although the scaling formally broke down, the fireball still essentially scaled with the explosion length scale r0. However, the discrete Lagrange particles (DLP) (phase 2) and the heterogeneous continuum model (HCM) of the DLP wakes (phase 3) did not scale with r0, and mean and RMS profiles could differ by a factor of 10 in some regions. This was because the DLP particles and wakes introduced an additional scale that broke the similarity conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗