Blast enhancement from metalized explosives

Shock Waves (2021) 31:203–230 https://doi.org/10.1007/s00193-021-00994-z ORIGINAL ARTICLE Blast enhancement from metalized explosives Q. Pontalier1 · J. Loiseau2 · S. Goroshin1 · F. Zhang3 · D. L. Frost1 Received: 27 May 2020 / Revised: 10 December 2020 / Accepted: 13 December 2020 / Published online: 5 May 2021 © The Author(s) 2021 Abstract Experiments are carried out to determine the effects of particle size and mass loading on the free-field blast wave from spherical, constant volume metalized explosive charges. The charges are comprised of gelled nitromethane with uniformly embedded aluminum, magnesium, or glass particles. Particle sizes are varied over an order of magnitude with particle mass fractions up to 50%. Peak blast overpressures are directly measured within the fireball with piezoelectric pressure gauges and outside the fireball are inferred by tracking the velocity of the blast wave and using the Rankine–Hugoniot relation. With the addition of inert particles, the peak blast overpressure is initially mitigated, but then recovers in the far field. For charges with reactive particles, the particles react promptly with oxidizers in the detonation products and release energy as early as within the first few hundred microseconds in all cases. The particle energy release enhances the peak blast overpressures in the far field by up to twice the values for a constant volume charge of the baseline homogenous explosive. By plotting the peak blast overpressure decay as a function of energy-scaled distance, it is inferred that at least half of the particle energy release contributes to the blast overpressure in the far field of higher mass loadings, and nearly all of the particle energy for a particle mass fraction of 10%. For aluminum, the blast augmentation is not a systematic function of particle size. This observation implies that conventional models for particle combustion that depend on particle surface area are not appropriate for describing the rapid aluminum reaction that occurs in the extreme conditions within the detonation products, which influences the blast wave propagation. Keywords Blast wave enhancement · Al/Mg particles · Reaction timescales · Particle–blast interaction · Mass fraction · Particle size 1 Introduction Communicated by D. Ranjan. Supplementary Information The online version contains supplementary material available at https://doi.org/10.1007/s00193-021-00994-z. * Q. Pontalier quentin.pontalier@mail.mcgill.ca D. L. Frost david.frost@mcgill.ca 1 McGill University, Macdonald Engineering Building, 817 Sherbrooke Street West, Montreal, QC H3A 0C3, Canada 2 Chemistry and Chemical Engineering Department, Royal Military College, 17 General Crerar Crescent, Kingston, ON K7K 7B4, Canada 3 Defence Research and Development Canada (DRDC) - Suffield, PO Box 4000, Stn Main, Medicine Hat, AB T1A 8K6, Canada The addition of metallic particles to explosive formulations has been a common technique for enhancing the performance of conventional high explosives for over a century, going back to the study in 1899 by Escales [1] and the patent by Roth [2] in 1900. Blast enhancement studies with metalized explosives are based upon the fact that the volumetric energy release from metal reactions is generally larger than the base explosive itself, and the assumption that particle reaction may occur sufficiently rapidly to enhance the blast wave strength. It is assumed that for a constant volume charge, the energy release from the metal reaction more than compensates for the reduction in explosive mass with the particle addition, so that the overall energy output and, in turn, the blast strength are augmented. Various reactive metals, such as aluminum, boron, magnesium, titanium, tungsten, and zirconium [3], have been added to explosives to enhance their performance. Aluminum is 13 Vol.:(0123456789) 204 the most common metal used, due to the high heat of formation of condensed aluminum oxide, low metal density, and relatively low cost. Energy release from the oxidation of the added metal particles can increase the blast overpressure and impulse, and hence such metalized explosives are often referred to as enhanced blast explosives (EBX). Comprehensive reviews of such explosives have been prepared by Frost and Zhang [4], Trzciński and Maiz [5], and Türker [6]. The so-called combined effects aluminized explosives (see [7] for a specific review) are formulated to enhance their metal pushing performance as well as blast effects [8]. The addition of particles to an explosive influences the pressure history behind the blast wave produced in several ways. For inert particles, as they are accelerated by the expanding combustion products, and subsequently decelerate, the particles exchange momentum with the surrounding gas and perturb the local pressure field non-uniformly as the particles typically form particle jets or filaments [9]. This effect leads to a reduction in peak blast overpressure in the near field as the particles accelerate and to a corresponding augmentation farther from the charge as the particles decelerate and return energy to the flow [10–13]. In the case of reactive particles, the contribution of the energy release from the particle reaction to the local pressure field is more complex. The coupling between the particle reaction and the blast wave pressure depends on the heat transfer rate between the particles and the surrounding gas [14], the oxidizing species present [15], the reaction rate of the particles, and the location of the energy deposition [16]. Other phenomena that may influence the strength of the blast wave from a metalized explosive include afterburning of the explosive detonation products [17, 18] or afterburning of the reactive gases generated by particle combustion. Secondary fragmentation of partially reacted particles due to relative flow, the interaction of reflected waves, or impact with nearby structure may also contribute to the blast overpressure augmentation. Even after more than a century of research, the physical mechanisms associated with the particle reaction and the augmentation of the blast pressure remain poorly understood. One of the main outstanding issues for aluminized explosives is quantifying the fraction of the energy release from the aluminum reaction that occurs on the timescale of the detonation propagation (i.e., on the order of microseconds), relative to the fraction that occurs on the timescale of the blast wave propagation (i.e., on the order of milliseconds). From research into the propagation of detonations in metalized explosives, or the so-called detonics regime, there is a growing body of evidence that metal particles can react on microsecond timescales behind the detonation wave and contribute to the ability to accelerate a metal casing [19–28]. The early signs of metal reaction seem to appear regardless of the explosive formulations used. In particular, metallic particles were observed to react at microsecond timescales 13 Q. Pontalier et al. in RDX [19–21], RDX–ammonium perchlorate (AP) [22], HMX [23, 24], nitroguanidine (NG), and bis(2,2,2-trinitroethyl)nitramine (BTNEN) [24]. Early-time contribution of Al reaction was also observed in hydrogen-free compounds such as benzotrifuroxan (BTN) [25], carbon-free substances such as ammonium nitrate (AN) [26], or liquids such as nitromethane (NM) [27–29]. For example, Loiseau et al. [28] measured the ability of mixtures of aluminum with gelled NM to accelerate a flyer plate. They found, for aluminum particles from 3.5 to 55 μm in diameter, that the detonation products delivered more energy to the flyer plate than the baseline explosive within 1–3 μs and concluded that the aluminum reaction begins very close to the detonation sonic plane. In mixtures of NM and packed beds of aluminum particles, Kato et al. [29] observed a detonation pressure increase due to Al reaction about 2.5 μs behind the leading shock for 8-μm Al particles. Additional evidence that aluminum can react promptly is given by recent experimental results that a mixture consisting of fine aluminum particles mixed with water (with 10% added hydrogen peroxide) can support the propagation of a detonation wave [30, 31]. This implies that aluminum powder must be capable of reacting with water on a microsecond timescale to ensure that the propagation of the detonation front is self-sustained. Since the oxidizers present in the detonation products of common CHNO explosives consist primarily of water, with smaller amounts of carbon dioxide and carbon monoxide present, it is reasonable to assume that aluminum can react with the water in the products at the extreme conditions behind the detonation front or during the early stages of product expansion. Despite the evidence that aluminum can react on a microsecond timescale within detonation products, previous models for the propagation of blast waves from aluminized explosives have largely ignored this observation. Within the blast and shock wave community, it is often assumed that no particle reaction occurs during the detonics stage, and the particles then react with the combustion products and surrounding air on millisecond timescales, with the reaction rate assumed based on the wealth of single-particle combustion experimental results available at relatively low pressure and quiescent conditions; for example, see the particle reaction time correlations of Beckstead [32]. Particle sizes used in aluminized explosives have ranged over four orders of magnitude, from hundreds of nanometers to millimeters. Although nanometric aluminum particles were expected to react more rapidly than larger particles and contribute significantly to the blast strength, the predicted improvement in blast performance with adding nanoaluminum powder to an explosive has largely not been realized. Contributing factors to this include the difficulties in deagglomerating and dispersing nanopowder uniformly within an explosive, and the lower active aluminum content of the powder relative to larger particles due to the preexisting oxide layer which represents a Blast enhancement from metalized explosives larger fraction of the overall mass. For particles in the micrometric range (e.g., 1–200 μm), previous studies showed little dependence of the particle size on the blast enhancement [15, 33, 34]. This is not consistent with most existing particle reaction models in which the reaction rate is determined by the amount of particle surface area available for reaction. We will return to this issue of the lack of dependence of particle size on blast performance and the implications for the particle reaction mechanism later in the Sect. 5. Another parameter that influences the blast performance of aluminized explosives is the particle mass loading. As metal particles are mixed into an explosive, the explosive is diluted, reducing the mass fraction of explosive in the charge. However, the energy release from the reaction of the metal particles, which is larger on a volumetric basis than that of the explosive itself, partially compensates for the dilution effect. If the particle mass fraction is maximized, two other factors come into play. First, the resulting heterogeneous explosive will be highly fuel-rich and there will be insufficient oxidizer in the detonation products to fully react the metal. Secondly, the heat loss to the metal particles will reduce the temperature of the products and may inhibit or prevent the reaction of the particles altogether. In an earlier extensive series of experiments with packed beds of metal particles saturated with the liquid explosive-sensitized nitromethane (the solid mass fraction varied from 76 to 80%, for aluminum particles larger than 8 μm ), we observed a range of different particle reaction behaviors and blast effects, depending on the particle and charge size [4]. In many cases, the particles did not fully react, making it difficult to systematically determine the influence of particle size and mass loading on the blast enhancement from the metal particles. Similar effects were also observed by Trzciński and Maiz [5] for fully dense nitromethane–particle mixtures. In addition to the particle mass loading, another parameter that strongly influences the blast performance from metalized explosives is the scale of the charge, which governs the residence time of the particles within the hot detonation products before expansion cooling of the products may quench further particle reaction. From experimental results with kg-scale charges containing packed beds of magnesium particles saturated with nitromethane, a critical charge size for particle reaction was clearly observed, which depends on particle size [35]. The particle reaction behavior and blast wave parameters will not be preserved with a change in scale. The importance of considering the scale in interpreting results regarding the effect of particle size on blast enhancement in aluminized explosives has been pointed out by Peuker et al. [15] who carried out experiments with subgram quantities of aluminized explosive. The trade-off between dilution and enhanced energy density effects with the addition of particles suggests a critical mass fraction should exist that optimizes the explosive 205 blast performance. The optimum mass fraction will depend on the type of reactive particle used as well as the nature of the explosive and is best determined experimentally. To find the optimum particle mass loading with the nitromethane–particle system, it is necessary to gel the NM to enable the suspension of the particles. This technique was used by Loiseau et al. [28] and Trzciński et al. [36], who carried out experimental studies in which nitromethane was gelled with the addition of 4 wt% PMMA to allow a reduction in the particle loading. Trzciński et al. [36] investigated the effect of adding different fractions (15–60 wt%) of an Al–Mg alloy on the detonation properties and blast wave properties in a confined area from 200 g cylindrical charges. The highest peak overpressure and impulse values at 2 m from the charge were found for explosives with a particle mass fraction of 30%. The mean particle size was 63 μm. In the present study, gelled nitromethane is also used as the base explosive, with various mass fractions of aluminum, magnesium, or glass particles added. Of particular interest is the effect of varying the particle size as well as the mass fraction on the blast wave properties. The motivation of the work is to determine if the rapid (sub-millisecond timescale) reaction of aluminum in detonation products that has been observed in earlier detonics experiments can be shown experimentally to directly contribute to blast wave performance for charges with a prototypical scale (kg quantities) and a particle mass loading in a range that ensures uniform reaction of the particles. In the following sections, the experiments will first be described (Sect. 2), followed by a description of the pressure gauge analysis and video data reduction procedure (Sect. 3). Presentation and discussion of the experimental results are given, respectively, in Sects. 4 and 5 before the conclusions (Sect. 6). Additional results are given in the Appendix. 2 Experimental overview The baseline explosive charge utilized a thin-walled G40 light bulb for a casing (1 mm thickness and nominal diameter of 12.7 cm), with the filament removed and filled with nitromethane (NM). The NM was gelled to the consistency of liquid honey to suspend the particles by adding 4 wt% polymethyl methacrylate (PMMA), with 0.4 wt% of 3M K1 glass microballoons (GMBs) added to sensitize the mixture. This explosive mixture was previously used in another study [28]. Various sized spherical aluminum powders (H-10, 13 μm dia.; H-50, 54 μm dia.; H-95, 114 μm dia.; Valimet, Stockton, CA) were added at mass fractions ranging from 10 to 50%. Since the charge volume was constant at about 1 L, the fraction of gelled NM, as well as the Al mass fraction, varied between trials. Tests were also carried out with spherical glass particles (+ 325/− 170 mesh; 66 μm 13 206 dia.) to serve as a non-reactive particle case for comparison, as well as a limited number of trials using spherical Mg powder (GRAN 17, 240 μm dia.; Reade, PA). The charges were detonated by inserting a detonator (Reynolds RP-83) into a booster charge of 85.7–92.7 g (depending on the trial, see trial matrix in Table 1) of C-4 formed by hand, placed within a plastic sphere and held in place at the center of the charge with a glass tube through which the detonator was inserted. The charge was placed on the top of a wooden post, with a height of burst of 1.52 m (HOB ≈ 12 times the charge diameter). For each trial, the blast overpressure was recorded with three fast-response (1.5 μs rise time) piezoelectric pressure transducers mounted in lollipop-style gauges (dia. 30 cm). The transducers were located at 1.0 m (PCB 113B26/PCB 113A31 depending of the trial), 1.5 m (PCB 113A31), and 2.0 m (PCB 113A31) from the charge, with the furthest gauge at almost 90◦ as shown schematically in Fig. 1. The gauge at 2.0 m is also visible behind the charge in Fig. 2. The distance of this gauge from the charge was chosen such that the ground-reflected wave arrived at the transducer location after the end of the positive phase of the blast wave signature so that the blast overpressure profile was not influenced by the reflected wave. The motion of the blast front was visualized using two Phantom v1610 high-speed video cameras recording at 16,000 frames/s (camera 1) and 35,000 frames/s (camera 2), respectively. The camera 1 operating at the slower speed had a higher resolution (1280 × 800 pixels) and was used Q. Pontalier et al. for the analysis of the shock wave motion. This camera was offset to the left relative to the charge to provide the maximum field of view, as shown in Fig. 1 ( C′ E = 1.67 m). The distance from the camera to the point C′ shown in Fig. 1 was CC′ = 34.50 m. The resulting absolute Euclidean distance from the camera to the charge was CE = 34.54 m. To enhance the visualization of the blast front, a zebra board (width of stripes is about 15 cm) was placed behind the charge at a distance of C′ Z = 5.00 m, resulting in a distance of CZ = 39.50 m between the camera and the backdrop. The camera was focused on the plane of the charge center. Using the two black fiducial markers on the backdrop (shown in Fig. 2), spaced 3.95 m apart, for scaling, and with the camera resolution, the physical size of one pixel in the plane of the backdrop is estimated to be 3.7 mm. The corresponding size of a pixel in the� plane of the charge is calculated to be 3.2 mm (= 3.7 × CC ). This value multiplied by the blast CZ radius measured in pixels gives the value Rapp . For all the trials, to increase again the field of view, the camera was rotated slightly to the left so that the charge center was just out of view to the right. Since the rotation angle is relatively small (≈ 0.4◦), the value of Rapp is not expected to vary significantly and, consequently, no additional correction has been applied. Camera 2 was set to a resolution of 1152 × 400 pixels and was placed adjacent to camera 1 slightly to the right as shown in Fig. 1. This camera was used exclusively for the qualitative analysis of the fireball structures in Sect. 4.1. Contrary to camera 1, the charge is Table 1 Trial matrix The baseline corresponds to the gelled NM (95.6 wt% NM + 0.4 wt% GMB + 4 wt% PMMA). P0 and T0 correspond, respectively, to the atmospheric pressure and temperature 13 Blast enhancement from metalized explosives 207 12 trials conducted with added aluminum powder (3 with H-10 Al, 5 with H-50 Al, 4 with H-95 Al), 4 with magnesium powder (GRAN 17 Mg), and 3 trials conducted with added spherical glass beads. A summary of the trials conducted is presented in Table 1, and the videos of selected trials are available online as supplementary materials. 3 Analysis methodology In this section, the different methods for extracting the blast properties are described. Since the blast properties may depend on the actual distance from the charge, to avoid confusion with the terminology, we denote the near field as the region up to 2 m from the charge, the mid-field region 2–3 m from the charge, and the far-field region beyond 3 m from the charge. 3.1 Blast overpressure profiles from pressure gauges Fig. 1 Experimental overview (top view). The cameras were focused on the charge plane (dashed blue line). When the blast is viewed through a certain angle (𝛼 ), the apparent blast radius ( Rapp) overestimates the actual blast radius (R). The optical axis of the camera 1 (CZ) is quasi-perpendicular to the charge plane. The pressure gauges denoted G1 (∢g1 = 54.1◦), G2 (∢g2 = 53.8◦), and G3 (∢g3 = 88.4◦) are located, respectively, EG1 = 1.0 m , EG2 = 1.5 m , and EG3 = 2.0 m from the charge. Camera 2 was only used to visualize the fireball structures (Sect. 4.1) Fig. 2 Photograph of the experimental setup with zebra board backdrop (front view). The scale shown is in the plane of the zebra board in the field of view of camera 2. In total, 20 trials were conducted with, respectively, 1 trial with the baseline charge (95.6 wt% gelled NM + 4 wt% PMMA + 0.4 wt% GMB), The pressure–time histories were recorded by three transducers placed at 1.0 m, 1.5 m, and 2.0 m from the charge. The transducers were located close enough to the charge that in most cases they were engulfed by the fireball, with the farthest transducer (at 2 m from the charge) typically near the maximum extent of the fireball. The non-uniform density field within the fireball, which is manifested by the conical jets that are visible projecting out of the fireball interface, leads to fluctuations in the pressure signals. Misalignment of the face of the lollipop gauge with the direction of travel of the blast wave can also generate expansion waves from the edge of the lollipop, which influences the pressure recorded shortly after the pressure reaches a maximum. An example of the overpressure history (i.e., the pressure above atmospheric pressure, ΔP = P − P0 ) recorded at 2.0 m from the charge is shown in Fig. 3 for the trial with the charge containing 30 wt% of H-50 Al. The peak blast overpressure (ΔPs) corresponds to the maximum overpressure value. Values of all peak overpressures at 2.0 m from the charge are included in the Appendix (Tables 5 and 6). The blast time of arrival corresponds to the time separating the beginning of the oscillatory signal generated by the burst of the detonator (exactly correlated with time zero in Fig. 3) and the time of the maximum overpressure. Note that the previous definitions may introduce error in the estimation of the peak blast overpressure and blast arrival times, especially if the blast front is not associated with the maximum recorded overpressure. This scenario may happen if the gauges are misaligned with the direction of propagation of the blast wave. Both the secondary shock and the ground-reflected wave are evident on the pressure signal. The timing of the arrival of the ground-reflected shock can 13 208 Q. Pontalier et al. 11 400 Blast (from video with parallax correction) Blast (from pressure gauges) Jet (from video) Bow shock (from video) Fireball (from video) 30% H-50 (pressure gauge data) 10 350 9 Ground-reflected shock 8 NF 250 Arrival time (ms) Overpressure (kPa) 300 200 Ps Secondary shock 150 100 MF FF 7 6 5 4 3 50 2 0 ta -50 1 1 2 3 4 5 Time (ms) Fig. 3 Overpressure profile captured at 2.0 m from the charge in the case of a charge containing H-50 aluminum at a mass fraction of 30% (trial #6) also be determined from the high-speed videos. The gauge located at 2.0 m is closest to the zebra board, although the reflected wave from the zebra board arrives at the gauge after a delay of 13.83 ms, which is not visible on the signal shown. The pressure profiles recorded within the fireball (at 1.0 m and 1.5 m from the charge) differ substantially from that of the (modified) Friedlander form [37], making it difficult to determine unambiguous values for the blast duration and peak overpressure. As a consequence, the trial-by-trial comparison of peak overpressure data is not considered at these distances. However, the blast times of arrival at 1.0 m are analyzed globally to evaluate the shortest aluminum reaction timescales. Due to the oscillatory nature of the pressure signal and the difficulty in choosing an unambiguous positive-phase blast duration, the calculated impulse values are not analyzed further. 3.2 Video analysis The blast front trajectories were determined through analysis of the high-speed videography images. To enhance the visibility of the blast wave, the function find edges was applied to the images using the open-source image processing software ImageJ. Since the vertical motion of the blast wave is perturbed by the flow of detonation products out of the tube holding the detonator, the blast front is tracked along a horizontal line emanating from the charge center. As the location of the blast front extracted from the images ( Rapp ) overestimates the actual location (R) due to the parallax effect, a correction is implemented before calculating the blast velocity. 13 0 0 1 2 3 4 5 6 Distance (m) Fig. 4 Trajectories of the blast wave with the parallax correction for the case of a charge with 30 wt% of H-50 (trial #6) aluminum particles. Also shown are the blast arrival times obtained from the three pressure transducers as well as the trajectory of the bow shock that is driven by a jet of particles that is projected ahead of the blast front. The approximate trajectory of the maximum extent of the fireball is also plotted. The abbreviations NF, MF, and FF correspond, respectively, to near field, mid field, and far field The formula used to calculate the actual distance from the blast front to the charge center is given as follows: R = Rapp × √ CC� R2app + CE2 − 2 × Rapp × C� E (1) where R < Rapp is generally excepted when R = Rapp = C� E. The blast trajectory with taking into account the parallax is plotted in Fig. 4. Note that the blast wave is not visible until it has propagated about 2 m due to the saturation of the video images from the luminosity of the fireball. Note also that data points are gathered by a group of 2 or 3, at regularly spaced intervals. These points are extracted when the blast front approaches the intersection of two consecutive stripes (shown in the background of Fig. 2), where the contrast is optimum. To determine the blast velocity, the blast trajectory is fitted by the monotonic fit of Kleine et al. [38] and then differentiated. The parallax correction is relatively small, with a typical error in the velocity of about 2% (4% when converted into pressure), with the error highest when the blast wave is furthest from the charge. An approximate estimation of the trajectory of the maximum extent of the fireball is also plotted in Fig. 4. As mentioned previously, the first two gauges are engulfed in the fireball and the last gauge is typically very close to the contact surface. Blast enhancement from metalized explosives The peak blast overpressure can be extracted by using the blast wave Mach number together with the Rankine–Hugoniot relation, provided the local temperature and pressure are known. For each trial, the local temperature (ranging from 21 to 26 °C) and atmospheric pressure (ranging from 92.71 to 93.83 kPa at an altitude of 700 m) were recorded (see Table 1) and used in the Rankine–Hugoniot relations. Further details of the method of determining the peak blast overpressure in this way can be found in an earlier publication [10]. Since the direct visualization of the blast front on the tracking line is not generally possible for distances from the charge less than 2.0 m, the peak blast overpressure is extrapolated in the near field by inferring the blast curvature, as shown in Fig. 5, for a direct comparison with gauge data. The results of this comparison are gathered in Table 7, in the Appendix. When extracting the blast wave trajectories, care must be taken to distinguish the actual blast front from the localized pressure disturbances (or bow shocks) created by the jets that are formed during the explosive event, as illustrated in Fig. 6. Due to the white color of these jet-like structures, they likely are comprised of agglomerated aluminum oxide particles or perhaps partially reacted aluminum in the case of large initial solid mass loading. The bow shocks may be either ahead of or behind the blast front. Since the two features are not propagating at the same speed, one propagating in ambient air and the other in the post-shocked medium, they ultimately merge after some time. In Fig. 4, trajectories of the jet and the bow shock visible ahead of the blast front are plotted. Initially separated in space, the blast front and the bow shock tend to converge while propagating outward. On the contrary, the bow shocks tend to separate from the Fig. 5 Processed image with the find edges function showing the extrapolation method of the blast front in the near-to-mid field, on the blast front tracking line 209 Fig. 6 Single frame from the video record of a test for a charge with 30 wt% of H-50 (trial #6) aluminum particles illustrating shock tracking against zebra board particle jets, due to aerodynamic drag with the surrounding atmosphere. The effect of the localized bow shocks on the peak blast overpressure has been discussed in previous publications [10, 11]. 4 Results 4.1 Visualization of fireball expansion and blast propagation To illustrate the dynamics of the fireball expansion and afterburning and the blast wave propagation, single frames from the video records of camera 2 have been extracted from a series of trials. Figure 7 shows the images for the detonation of the baseline gelled nitromethane explosive. After about 0.5 ms, the fireball is highly luminous, but the processed image indicates that small-scale perturbations and jet-like structures are present on the fireball interface, likely due to action of the Rayleigh–Taylor (RT) instability. After 1 ms, the expansion has cooled the fireball gases and the interface between the detonation products and surrounding atmosphere appears highly convoluted. The fireball continues to expand and cool until the secondary and the ground-reflected waves propagate up through the fireball, enhancing the mixing and afterburning of the products. The top interface of the highly luminous region in the image at 4 ms corresponds to the reflected wave that is propagating upward at this time. Figure 8 shows the effect of adding 30 wt% of inert glass particles to the gelled nitromethane. The image taken at 0.314 ms shows the luminous fireball, and the processed image indicates that conical jet-like structures are present, perhaps caused by entrainment of the glass particles into the RT perturbations. Fine streaks of particles are also present in this image near the fireball interface, which may be from 13 210 Fig. 7 Single frames from the video record of the detonation of a charge containing the baseline gelled nitromethane explosive (trial #0). The images on the left are the original video images, and the images on the right are the same images, but after an edge Laplacian 5 × 5 filter has been applied to highlight the gradients in the image. Times for the image are (top to bottom): 0 ms, 0.514 ms, 1.000 ms, 2.000 ms, and 4.000 ms, respectively. The arrows indicate the locations of the blast wave Fig. 8 Single frames from the video record of the detonation of a charge containing 30 wt% glass particles (trial #18). Similar to Fig. 7, the images on the left are the original video images and the same images are on the right, after an edge Laplacian 5 × 5 filter has been applied to highlight the gradients in the image. Times for the image are (top to bottom): 0 ms, 0.314 ms, 1.000 ms, 2.000 ms, and 4.000 ms, respectively. The arrows indicate the locations of the blast wave 13 Q. Pontalier et al. Blast enhancement from metalized explosives 211 ∑imax ∑jmax L= i=1 j=1 Pi,j 3500 Baseline (#0) 30% H-50 (#6) 30% Glass (#18) 3000 2500 Averaged pixel intensity the motion of the fine glass particles from the shattering of the glass light bulb casing. At times of 1 ms and later, the fireball has expanded and cooled. The secondary and the ground-reflected waves do not induce significant afterburning of the soot from the NM products, due to the absorption of heat by the added inert glass. Figure 9 illustrates the fireball afterburning and blast wave propagation in the case of a charge with 30 wt% H-50 Al particles. Prominent conical particle jets form and project ahead of the primary blast front, which is only visible emerging from the fireball interface mixing region after about 2 ms. The fireball intensity does not diminish during the recording time of the video (10 ms), as the afterburning of the soot from the detonation products as well as energy release from the aluminum reaction continues to contribute to the fireball luminosity. The reader can refer to the videos available online for a dynamic visualization of Figs. 7, 8, and 9. To illustrate the variation of luminosity between these three trials (#0, #6, #18), the spatially averaged pixel intensity ( L ) is defined as: 2000 1500 1000 500 0 -500 -1000 0 5 10 15 20 Time (ms) Fig. 10 Comparison of spatially-averaged pixel intensity ( L ) as a function of time between the baseline explosive (trial #0) and with mixtures with 30 wt% H-50 (trial #6) and 30 wt% glass beads (trial #18) (2) N with imax = 1280 , the number of pixels by width; jmax = 800 , the number of pixels by height; N = imax × jmax = 1280 × 800 = 1,024,000, the total number of pixels in a frame; and Pi,j , the intensity value of the pixel of coordinate (i, j) from the original 16-bit cine (Phantom format) files, is plotted as a function of time, in Fig. 10. The cine file format is preferred to the RGB format to avoid correction of Fig. 9 Single frames from the video record of the detonation of a charge containing 30 wt% H-50 Al particles (trial #6). As shown in Figs. 7 and 8, the images on the left are the original video images and the images on the right have been processed with an edge Laplacian 5 × 5 filter to highlight the gradients in the image. Times for the image are (top to bottom): 0 ms, 0.657 ms, 1.000 ms, 2.029 ms, and 4.000 ms, respectively. The arrows indicate the locations of the blast wave 13 212 Q. Pontalier et al. the pixel intensity value according to pixel color. For each trial, the intensity of the first image is subtracted from all the images, so the L value is equal to 0 before detonation of the charge. The detonation time corresponds to the start of the sharp rise of the luminosity that saturates the camera. The integrated intensity decreases sharply after the detonation but increases again during the afterburning phase which occurs after the primary blast wave reflects off the ground, after about 2 ms. The value of L captured at a mass fraction of 30% H-50 Al is significantly larger than that of the baseline. Conversely, adding the same mass fraction of glass particles reduces the luminosity considerably. Values of L can be even negative due to the soot formed that causes the video image to be darker than the initial image before the detonation. Note also the absence of an afterburning phase for the inert particles. The analysis of the spatially integrated luminescence presented here remains purely qualitative due to the numerous camera settings that influence the pixel intensity, such as the lens aperture, exposure time, and extreme dynamic range (EDR) setting, which is a feature of the operation of Phantom cameras that reduces the specific intensity of pixels that are approaching saturation. 4.2 Fireball expansion and blast propagation (pressure gauges and videography) The relative position of the maximum extent of the luminous fireball and the shock wave for the three previous trials is shown in Fig. 11. Note that the maximum extent of the fireball is approximate and cannot be considered a precisely defined measurement since in cases where some of 11 Blast - Baseline (#0) Fireball - Baseline (#0) Blast - 30% H50 (#6) Fireball - 30% H50 (#6) Blast - 30% Glass (#18) Fireball - 30% Glass (#18) 10 9 Arrival time (ms) 8 7 6 5 4 4.3 Blast time of arrival in the near‑to‑mid field (pressure gauges) 3 2 1 0 the pixels on the camera sensor are saturated, some charge leakage (or pixel bleed) between adjacent pixels can occur. The blast arrival times captured with the pressure gauges are also reported. In general, the blast time of arrival does not match perfectly with the blast time of arrival obtained with the videography analysis in the mid field. The discrepancy is explained by the fact that the pressure gauge at 2.0 m from the charge is not located within the field of view of camera 1 (see Fig. 1) and implies that the blast propagation is not isotropic. Nevertheless, the blast times of arrival obtained from videography measurements are offset to match the arrival time at the pressure gauge in the mid field (Fig. 11). As mentioned before, for the reactive case (30% H-50), the first two pressure gauges are engulfed in the fireball and the last gauge located typically very close to the fireball surface. Analyzing the blast times of arrival is of interest to highlight the contribution of the particle reaction to the blast wave propagation. The energy deposited into the gas by the reaction of the particles may couple with the blast front and increase its velocity. Conversely, the absorption of the explosive energy by inert particles may reduce the energy available to support the blast, hence decreasing the blast velocity. Figure 11 indicates that the blast times of arrival are smaller for the trial with aluminum particles compared to the baseline from 0.426 ms (at 1.0 m) and more significantly from 0.830 ms (at 1.5 m). This suggests that the aluminum reaction contributes to increasing the blast velocity. Conversely, the blast time of arrival is higher when glass particles are added to the explosive mixture, indicating that the addition of inert particles contributes to decreasing the blast velocity. A contraction of the fireball size after a few milliseconds due to rarefaction waves is visible for the inert case (30% glass) and for the baseline case, but not observed for the trial with aluminum. The comparison between the trajectories of the fireballs indicates that the fireballs generated by the aluminized explosives grow to a larger size than the fireballs generated with the baseline explosive. Conversely, the size of the fireball generated by the explosive diluted with inert particles tends to be smaller than that of the baseline explosive. 0 1 2 3 4 5 Distance (m) Fig. 11 Blast front and fireball trajectories plotted for the baseline explosive (trial #0) and with mixtures with 30 wt% H-50 (trial #6) and 30 wt% glass beads (trial #18) 13 Figure 12a summarizes the times of arrival, for a given type of powder, at a distance of 2.0 m from the charge. Values of times of arrival for the different trials are normalized by the time of arrival for the baseline charge (i.e., ta ∕ta-baseline ). The performance of the baseline is equal to 1 and is represented by the dashed blue line. When the explosive is diluted with glass particles, the blast times of arrival are systematically delayed compared to the baseline. The arrival times increase Blast enhancement from metalized explosives 213 1.25 1.15 10% 20% 30% 40% 50% 2 Normalized peak overpressure Normalized time of arrival 1.2 1.1 1.05 1 0.95 0.9 0.85 1.6 1.4 1.2 1 0.8 0.8 0.75 1.8 10% 20% 30% 40% 50% H-10 H-50 H-95 Mg Glass 0.6 H-10 H-50 H-95 Materials Materials (a) (b) Mg Glass Fig. 12 a Effect of the mass fraction on the normalized blast times of arrival (ta ∕ta-baseline) calculated from pressure gauge data at 2.0 m from the charge and sorted by type of powders. The baseline value (gelled NM without particles) is represented by the horizontal dashed blue line. b Same plot but for the normalized peak blast overpressure (ΔPs ∕ΔPs-baseline) monotonically with the mass fraction by up to 23% of the baseline value at a mass fraction of 50%. Conversely, the addition of reactive (Al or Mg) particles in the gelled NM systematically reduces the blast times of arrival, for all trials tested, compared to the baseline case. The reduction in the arrival times is generally a monotonic function of the mass fraction. The addition of 20 wt% of magnesium particles reduces the blast time of arrival by 13% compared to the pure gelled NM. The addition of 40–50 wt% of aluminum particles typically reduces the blast time of arrival by 17–23% compared to the baseline explosive depending on the powder considered. Figure 12a also indicates that the particle size of the aluminum particles has very little influence on the blast arrival times. Note the variation of the normalized arrival times between the repeated trials (H-50 particles at a mass fraction of 40%). This discrepancy suggests that the small variations of the normalized arrival times, when comparing the different aluminum particle sizes, are on the order of the stochastic variation of the experiments. The blast times of arrival for each trial, recorded by the pressure gauges, are presented in Tables 3 and 4 in the Appendix. particles, compared to the baseline. Anomalously, the charge with 20 wt% magnesium particles outperforms all the other charges at enhancing the peak blast overpressures with an augmentation of 89% compared to the baseline, although data were not obtained for higher Mg mass loadings due to problems with the pressure sensors. For charges with aluminum particles, the highest performance was obtained at mass loadings of 30 wt% (H-10, H-50) or 40 wt% (H-95), with an augmentation of 47–65% compared to the baseline value. This suggests that the peak blast overpressure enhancement does not change monotonically with the mass fraction, contrary to the blast arrival times. This lack of correlation between the arrival times and peak overpressure may be due to the error associated with the extraction method mentioned in Sect. 3. Conversely, the addition of 50 wt% of glass particles in the gelled NM mitigates the peak overpressure by around 27%. As observed for the blast arrival times, the normalized peak blast overpressures show little dependence on the aluminum particle size excepted for a mass fraction of 40%. At this mass fraction, the normalized peak blast overpressures with H-95 particles outperform those obtained with H-10 particles by 44%. Nevertheless, this discrepancy is cautionary due to the stochastic variation between repeated trials (see the results for H-50 particles). All the peak overpressure values are presented in Tables 5 and 6 in the Appendix. 4.4 Peak blast overpressures in the mid field (pressure gauge) Similar to Fig. 12a, the normalized peak blast overpressure values at the same distance are reported in Fig. 12b. The graph shows that the peak blast overpressures are systematically augmented for all charges tested with reactive materials and mitigated for all the charges tested with inert glass 13 214 Q. Pontalier et al. 250 1.3 Baseline (#0) 10% Glass (#17) 30% Glass (#18) 50% Glass (#19) Normalized peak overpressure Peak overpressure (kPa) 200 1.2 150 100 50 10% Glass (#17) 30% Glass (#18) 50% Glass (#19) 1.1 1 0.9 0.8 0.7 0.6 0 1.5 2 2.5 3 3.5 4 4.5 5 Distance (m) (a) 0.5 1.5 2 2.5 3 3.5 4 4.5 5 Distance (m) (b) Fig. 13 a Effect of particle mass fraction on peak blast wave overpressure decays for glass particles. The baseline charge (#0) contains only gelled NM with no particles. b Peak overpressure efficiencies relative to the baseline explosive (dashed blue line) 4.5 Effect of particle material and mass loading on blast decay in the mid‑to‑far field (videography) Figure 13a shows the decay of the peak blast overpressure with distance, in the mid-to-far field, for trials with glass particles, with various mass fractions of the additive powder. Note that the relative errors resulting mainly from the tracking of the blast front on the videos are typically of the order of the size of the data markers and are consequently not shown. To highlight the performance, the normalized peak overpressures are also shown in Fig. 13b. As shown in Fig. 12, the performance of the baseline trial is equal to 1 and is represented by a blue dashed line. The peak blast overpressure decays for all the other materials are reported in the Appendix in Figs. 19 and 20. Adding the glass particles reduces the strength of the blast wave in the mid field, with the reduction in pressure increasing as more glass powder is added (by up to 37.5% at around 2 m). In the far field, the peak blast overpressure recovers almost to the same value as the base explosive charge for 50 wt% addition of glass and even exceeds the baseline case for 10 wt% and 30 wt% addition, by up to 22.8% even though the charges with added glass have a reduced amount of explosive. This behavior is consistent with the earlier results from Pontalier et al. [10–13], who attributed the near-field reduction in blast pressure to the energy transfer from the detonation products to the kinetic energy and heating of the inert particles, and the later time 13 blast pressure recovery to the pressure perturbations induced by the particle–flow interactions. With the addition of reactive particles, Fig. 14 shows that in all cases, the particle reaction enhances the strength of the blast wave over the baseline explosive, with the enhancement greatest in the far field. In the mid field, i.e., just outside the maximum fireball radius, the particle reaction counteracts the pressure reduction observed for inert particle addition, such that the peak blast overpressure is higher than that of the baseline case by between 5% (10 wt% H-95) and 27% (40 wt% H-10) at around 2.2 m. Further from the charge, the combination of the particle reaction and the motion of the particles within the products sustains the peak blast overpressure at a higher level than for the homogeneous liquid explosive, so that the peak blast enhancement with the addition of reactive particles can attain a value of up to almost twice the baseline value (for 20 wt% Mg), at 4.5 m. The effect of the mass fraction of added particles, for a given type of particle, on the blast overpressure decay is shown in Figs. 13 and 14. For the addition of glass particles, the attenuation of the blast wave overpressure increases monotonically with the mass fraction, as expected. Conversely, with added H-10 aluminum particles, the blast enhancement increases monotonically with the mass fraction. In the case of H-50 and H-95 particles, the blast enhancement depends on the location from the charge and is not necessarily a monotonic function of the mass fraction. In the far field, the maximum peak overpressure is attained at a mass fraction of 30–40% for H-50 and 50% for H-95. However, this assertion is subject to the scatter observed Blast enhancement from metalized explosives 215 2 2 1.8 Normalized peak overpressure Normalized peak overpressure 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) 1.6 1.4 1.2 1 1.5 1.8 40% H-50 (#7) 40% H-50 (#8) 30% H-50 (#6) 20% H-50 (#5) 10% H-50 (#4) 1.6 1.4 1.2 1 2 2.5 3 3.5 4 4.5 5 1.5 2 2.5 Distance (m) (a) H-10 Normalized peak overpressure Normalized peak overpressure 4 4.5 5 4 4.5 5 2 50% H-95 (#12) 40% H-95 (#11) 30% H-95 (#10) 10% H-95 (#9) 1.6 1.4 1.2 1 1.5 3.5 (b) H-50 2 1.8 3 Distance (m) 1.8 40% Mg (#16) 30% Mg (#15) 20% Mg (#14) 10% Mg (#13) 1.6 1.4 1.2 1 2 2.5 3 3.5 4 4.5 5 Distance (m) (c) H-95 1.5 2 2.5 3 3.5 Distance (m) (d) Mg Fig. 14 Effect of the mass fraction on the peak overpressure efficiencies relative to the baseline explosive (dashed blue line) plotted for a H-10 Al, b H-50 Al, c H-95 Al, and d GRAN 17 Mg between the two repeated trials (40% H-50, trials #7 and #8). The results for Mg indicate that the blast enhancement is greatest for a particle mass fraction of 20%, at least in the far field. Note that for aluminum particles, the blast enhancement increases more rapidly when the mass fraction is high (30–50%) compared to low mass fractions (10–20%). For aluminum particles, at a specific mass fraction of 30%, Fig. 15 shows that the peak blast overpressure exhibits no systematic dependence on particle size. The effects of particle size on blast overpressure at other mass fractions are shown in Fig. 20 in the Appendix. The results for aluminum particles again do not show a systematic variation with particle size. In the mid field, the smallest aluminum particles yield slightly higher normalized peak overpressures for all mass fractions, whereas in the far field, the larger particles generate larger normalized overpressures at a mass fraction of 10% (after 3.5 m). However, at a mass fraction of 40%, H-10 Al particles yield the highest peak overpressure values in the far field. Figure 20 also shows that with the addition of magnesium particles the normalized peak overpressures in the mid field are comparable to those for aluminum particles, but are lower than the aluminum values in the far field for mass fractions of 10 wt%, 30 wt%, and 40 wt%, which is consistent with the 13 216 Q. Pontalier et al. Table 2 Blast times of arrival for all trials at 1.0 m from the charge indicating the early contribution of the energy release of Al/Mg particles to the blast wave Normalized peak overpressure 2 30% H-10 (#2) 30% H-50 (#6) 30% H-95 (#10) 30% Mg (#15) 30% Glass (#18) 1.5 1 0.5 1.5 2 2.5 3 3.5 4 4.5 5 Distance (m) Fig. 15 Effect of material and particle size on the peak overpressure efficiencies relative to the baseline explosive (dashed blue line), for a particle mass fraction of 30% smaller specific energy release from magnesium relative to aluminum. The comparison between the peak overpressure values at 2.0 meters from the charge, obtained with the pressure gauge or the videography analysis, is shown in the Appendix, in Table 7. The peak blast overpressure discrepancies between the two methods are not systematic. This observation suggests that the strength of the blast wave is not uniform in every direction. 5 Discussion The results presented in this study show strong evidence of prompt reaction of the aluminum or magnesium particles embedded in the gelled NM mixture. The reactions are sufficiently exothermic to counteract the heat and momentum transfer from the detonation products to the particles and contribute to enhancing the luminosity (see Fig. 10) and size (see Fig. 11) of the fireball, as well as the velocity and overpressure of the blast wave (see Fig. 12). Conversely, adding inert glass particles serves to absorb some of the explosive energy, which reduces the luminosity and size of the fireball as well as the velocity and overpressure of the blast wave, at least in the mid field. Farther from the charge, some of the kinetic energy of the particles is transferred back to the surrounding gas which leads to a recovery of the blast wave overpressure, in some cases to values above that of the baseline explosive alone (see Fig. 13). This blast recovery phenomenon, attributed to the pressure disturbances created by the particle–flow interactions, was described in a previous 13 work with explosive charges containing a variety of different embedded particles [10]. The contribution of these disturbances may also depend on the particle size. Small particles tend to adapt more quickly to the flow in comparison with larger particles and hence transfer their energy to the flow more rapidly [11]. The influence of the energy release from the particles on the blast wave propagation depends on the location of the energy release. If the energy release occurs immediately behind the blast wave, the local increase in gas pressure will promptly enhance the blast overpressure. Even if the particle energy release occurs some distance behind the blast wave, the pressure disturbances generated can catch up to the blast wave in the mid or far field since the local sound speed behind the blast wave is increased due to the shock heating as well as the particle reaction. In general, with the addition of reactive particles, the blast wave peak overpressure increases linearly with distance from the charge, when normalized relative to the baseline explosive, due to the effect noted above (see Fig. 14). Effectively, the particle reaction continues to contribute to augmenting the blast wave overpressure in the far field due to the continual arrival of pressure disturbances at the blast wave front. In the near field within the fireball, the blast wave is not directly visible, but the motion of the blast wave can be tracked using signals from the pressure gauges that are engulfed by the fireball. For example, the times for the arrival of the blast wave at the first gauge, located 1.0 m from Blast enhancement from metalized explosives the charge, are reported in Table 2 and are less than that for the baseline explosive, which is consistent with the higher blast wave speed noted above. From Fig. 11, it is evident that the combustion products fireball extends beyond this pressure gauge for the baseline explosive charge and eventually grows to a radius of about 1.5 m. Hence, it is expected that when aluminum particles are added to the explosive, by the time that the blast wave reaches this first pressure gauge, the aluminum particles will have been exposed primarily to the oxidizers within the detonation products, which predominantly consist of water for charges with NM. The minimum observed arrival time is 353 μs, and hence the particles must have started to react before this time. Since the nearest pressure transducer is located 1 m from the charge, with the current diagnostics used, it is not possible to resolve the rate of energy release from the particles at earlier times than this. It is of interest to consider the nature of this early-time reaction of aluminum with the water in the detonation products. This reaction produces molecular hydrogen gas which will contribute to inflating the size of the fireball. This hot hydrogen gas is highly reactive, and it will react promptly when it mixes with the oxygen in the surrounding atmosphere at the contact surface between the detonation products and air. The rate of reaction of the hydrogen will be limited by the mixing rate with the surrounding atmosphere which in turn is influenced by the rate of growth of the perturbations that develop at the contact surface due to the Rayleigh–Taylor instability as the interface decelerates. In Fig. 11, the growth of the radius of the fireball with added aluminum particles begins to separate from that of the baseline explosive as early as 200 μs. This may be due to the production of hydrogen and subsequent reaction with the surrounding air. After this time, the fireball continues to expand over several milliseconds. The reactions within the fireball, due to the particles as well as from residual soot from the baseline oxygen-deficient explosive, are also enhanced by the mixing induced by the passage of the secondary and the ground-reflected blast waves after about 2 ms, which leads to an increase in fireball luminosity (see Fig. 10). The oxygen balance of an explosive provides a measure of the availability of oxidizers for reacting with aluminum in the explosive decomposition products. Alternatively, equilibrium calculations can be carried out to estimate the amount of aluminum that could theoretically react within the detonation reaction zone as well as during adiabatic expansion of the products. Calculations were carried out using the equilibrium code CHEETAH 2.0 [39] to determine the maximum amount of aluminum that could potentially completely react with the products of the detonation of the gelled NM at the CJ detonation state as well as if the products were allowed to adiabatically expand to 1 atm. The result of 217 these calculations with the BKWC library indicated that a maximum of 26 wt% of Al mixed with the gelled NM could potentially completely react by the end of the detonation zone assuming that equilibrium conditions are attained. Furthermore, if the products are allowed to expand adiabatically to atmospheric pressure, additional reaction of Al is possible and in this case an initial mixture containing 36 wt% would result in complete reaction of the Al after expansion. If more species are allowed in the detonation products (JCZS library), the percentages go up to 39% and 63%, respectively. These values provide an upper bound for the amount of aluminum that could react with the detonation products without additional mixing with oxygen in the surrounding air. Due to the finite rates of aluminum oxidation and mixing, the actual amount of Al that could react in the products will be less. An important observation from the present experimental results is that the blast enhancement with added aluminum particles does not systematically depend on particle size. Despite varying the particle size by one order of magnitude (or the surface area by two orders of magnitudes), the normalized blast overpressure was not significantly different for the various particle sizes, and in the far field the results did not vary systematically with particle size at the various mass loadings considered. In a previous work, Peuker et al. [15] found negligible particle size effects for 3–40-μm aluminum particles on blast enhancement from small-scale aluminized explosive charges. These results suggest that conventional scaling laws for aluminum particle combustion based on surface area correlations that are derived from single-particle experimental results at quiescent atmospheric conditions must be modified to describe the particle reaction within the extreme conditions within the detonation products at early times. In the detonation reaction zone for condensed explosives, pressures are on the order of tens of gigapascals [14], with temperatures ranging from 3000 to 5000 K [40] and slip velocities above 1 km/s [41]. Under these extreme conditions, the particle reaction mechanism and transport properties are completely unknown. Within the detonation products, the density of the flow is initially of the same order as the density of the metallic particles such that extremely high heat transfer rates may be possible. In supercritical conditions, aluminum oxide dissolves in water [42], and hence the preexisting oxide coating on aluminum particles provides no barrier against the dissolution of the particle and subsequent anaerobic reaction with the detonation products. Note that the reaction of the aluminum with the oxidizers in the detonation products (primarily H2 O) is not coincident with the energy release from the reactions. The formation of some aluminum sub-oxides (e.g., AlO, Al2 O) is weakly exothermic, but the primary energy release is associated with the chemical condensation of the stable aluminum oxide Al2 O3, 13 218 Q. Pontalier et al. which occurs only after a degree of expansion of the products has taken place. In contrast with the lack of systematic dependence of blast performance on the particle size, the ratio of particle mass to the explosive mass does influence the blast enhancement. As the particle mass loading is increased for a fixed volume, the explosive is diluted. Initially, the energy release from the added particles overcomes this dilution effect, as the maximum energy release for complete oxidation of the particles is several times that of the energy release of the explosive itself (31.5 kJ/g for Al [4], 24.7 kJ/g for Mg [43] vs 6.4 kJ/g for NM [44]). Eventually, with a sufficiently high solid mass fraction, the heat and momentum losses to the particles will weaken the blast strength. This suggests that there exists an optimal particle loading for maximum blast strength, which will depend on both the particle and explosive properties. In the present experiments, the normalized peak overpressure is maximized for mass fractions between 30 and 50% for aluminum particles (see Fig. 14). The smallest mass loading (10%) consistently gave the smallest blast overpressures, and the overpressure augmentation increases with distance from the charge, possible due to the energy release from hydrogen gas as it mixes and rapidly burns with the surrounding air. 5.1 Energy‑based blast overpressure scaling To compare the blast wave overpressure decay with distance for explosive charges of different sizes, it is conventional to scale the distance with the mass of explosive (Hopkinson scaling [45, 46]), or the explosive energy release, which depends on the explosive mass. For a multi-component explosive, it appears to be more meaningful to scale the distance with a parameter that depends on the total energy release, or the sum of the energy released by each component, which in general will be different, rather than the total mass of the charge [4]. Using Bach and Lee’s convention [47], we define a scaled distance R0 (or explosive characteristic length), defined by the following expression: ( R0 = Etotal 𝜌0 c20 kj ) 1 j+1 (3) where 𝜌0 and c20 are, respectively, the density and the speed of sound of the blast propagating medium (atmospheric air in this case); j and kj are geometrical factor set, respectively, at 2 and 4𝜋 for a 3D spherical blast wave. If it is assumed that the ambient air is an ideal gas, the previous equation can be rewritten as follows: 13 R0 = ( Etotal 4𝜋𝛾P0 )1 3 (4) where 𝛾 is the ratio of specific heats, taken to be 1.405 for air (in the mid-to-far field) and P0 is the value of the atmospheric pressure, which changes from trial to trial. To account for the multiple sources of energy release, Etotal is defined as follows: Etotal = mNM × qNM + mC4 × qC4 + 𝜀 × mme × qme . (5) Here, m (g) is the mass of each component. For mNM , the mass of explosive is calculated without the 4.4 wt% inert additives (GMB + PMMA). Note that in the previous definition, the total energy Etotal is assumed to be released instantaneously after the detonation. The quantity q is the energy release per unit mass of each component (kJ/g), and the last term refers to the energy release from the metal combustion, with the subscript “me” (or metal) corresponding to either Al or Mg, depending on the trial. The parameter 𝜀, which may vary between 0 and 1, corresponds to the fraction of the energy released from the metal combustion that is assumed to contribute to supporting the blast wave propagation. Note that this 𝜀 value does not represent the amount of Al or Mg that reacts. Indeed, even if the timescale of particle energy release is short, it is not instantaneous, as assumed in (4). This finite timescale for the energy release implies that the percentage of metal reacting is generally higher than the 𝜀 value. In the limit of long times for the energy release of the metal particles, 𝜀 will tend to zero, whereas for instantaneous and complete energy release 𝜀 will approach unity. Note that since the charges in the present investigation have a constant volume, the relative contribution of the base explosive and the particle reaction to the blast overpressure will change as the particle material or mass loading is changed. Values of energy release for complete combustion of the various components, i.e., qNM , qC4 , qAl , and qMg , are set, respectively, to 6.4 kJ/g [44], 6.7 kJ/g [48], 31.5 kJ/g [4], and 24.7 kJ/g [43]. Note that the values of heats of combustion for magnesium and aluminum particles are those for complete reactions with molecular oxygen. However, these values are equivalent to the total heat of combustion assuming that the particles first react primarily with water in the detonation products to form molecular hydrogen gas, followed by the combustion of the hydrogen with the oxygen present in air, i.e., { 2Al(s) + 3H2 O(g) → Al2 O3 (s) + 3H2 (g) qAl (31.5 kJ/g) 3H2 (g) + 1.5O2 (g) → 3H2 O(g) (6) and Blast enhancement from metalized explosives 219 { Mg(s) + H2 O(g) → MgO(s) + H2 (g) qMg (24.7 kJ/g) . (7) H2 (g) + 0.5O2 (g) → H2 O(g) If the two reactions given in each case above are added together, as expected the resulting reaction gives the direct oxidation of the metal with molecular oxygen to form the metal oxide. In the case of aluminum, complete oxidation of aluminum to Al2 O3 releases an energy of 31.5 kJ/g. The two intermediate reactions given above release roughly equal amounts of energy. In particular, the first reaction of aluminum with water releases 15.1 kJ/g, with the remainder of the energy released by the reaction of hydrogen gas. It should be reemphasized that the parameter 𝜀 introduced above corresponds to the fraction of the energy release that is assumed to contribute to the blast wave propagation, and cannot be directly related to the fraction of aluminum that reacts. For example, a value of 𝜀 of 0.5 could correspond roughly to either half of the aluminum present fully reacting, or all of the aluminum present reacting only with water to form hydrogen. Water vapor is not the only oxidizer present in the detonation products of NM [49] but is the most prevalent. According to CHEETAH 2.0 (BKWC library), the detonation product gases of NM at the CJ detonation state are composed, by order of prevalence, of water ( H2 O ), carbon monoxide (CO), nitrogen ( N2 ), hydrogen ( H2 ), and carbon dioxide ( CO2 ). The quantity of water vapor exceeds the quantity of carbon monoxide by more than one order of magnitude. The Al/Mg particles could alternatively react with all the other oxidizers than water. In particular, previous publications have shown the possibility of metallic particles to react, at microsecond timescales, in hydrogen-free explosives, through other reaction pathways including reaction with CO2 as an oxidizer [25]. This implies that early-time metal reaction can contribute to the blast strength through different reaction pathways and that for different explosive compositions. Nevertheless, the reaction timescales and pathways may influence the overall blast wave dynamics as well as the peak blast overpressure values. The q values presented in (6) and (7) are the most energetic ones and do not change if the primary oxidizer is substituted. For instance, the Al reaction with CO2 and the subsequent combustion of the produced CO with the O2 of the air give the same value of energy release (qAl = 31.5 kJ/g). If the peak blast overpressure results for charges containing inert glass particles are plotted versus scaled distance, the results are shown in Fig. 16. Figure 16a shows the overpressure decay in dimensional coordinates (same as the plot in Fig. 13a). In Fig. 16b, the results are shown versus non-dimensional distance R/ R0 [with 𝜀 = 0 in the definition of Etotal in (5)]. As expected, the curves collapse reasonably well, with small discrepancies with the baseline results apparent in the mid and far field, suggesting that scaling with the energy released by the explosive (both NM and C4) for a liquid explosive mixture diluted by inert particles is reasonable. In the scaled coordinates, the small reduction in overpressure with the addition of particles in the near field 250 250 Baseline (#0) 10% Glass (#17) 30% Glass (#18) 50% Glass (#19) 200 Peak overpressure (kPa) Peak overpressure (kPa) 200 150 100 150 100 50 50 0 1.5 Baseline (#0) 10% Glass (#17) 30% Glass (#18) 50% Glass (#19) 0 2 2.5 3 3.5 Distance (m) (a) 4 4.5 5 1 1.5 2 2.5 3 3.5 Scaled distance (R/R 0) (b) Fig. 16 Peak blast overpressure decays with distance for glass particles with a dimensional radius, and b with scaled distance, R∕R0 , where R0 accounts for the explosive energy only [NM+C4, 𝜀 = 0 in (5)] 13 220 Q. Pontalier et al. 300 300 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) Baseline (#0) 250 Peak overpressure (kPa) Peak overpressure (kPa) 250 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) Baseline (#0) 200 150 100 50 200 150 100 50 0 1.5 2 2.5 3 3.5 4 4.5 0 5 1 Distance (m) 1.5 2 (a) H-10 300 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) Baseline (#0) 200 150 100 50 1 1.5 2 2.5 3 150 100 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) Baseline (#0) 200 150 100 50 1 1.5 2 2.5 Scaled distance (R/R0) (e) H-10, ε = 0.5 1 1.5 2 (d) H-10, ε = 0.25 300 250 0 2.5 Scaled distance (R/R0) (c) H-10, ε = 1 Peak overpressure (kPa) 200 50 Scaled distance (R/R0) 13 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) Baseline (#0) 250 Peak overpressure (kPa) Peak overpressure (kPa) 250 0 3 (b) H-10, ε = 0 300 0 2.5 Scaled distance (R/R0) 3 3 Blast enhancement from metalized explosives ◂Fig. 17 Peak blast overpressure decays with distance for H-10 alu- minum particles with a dimensional radius, and b with scaled distance, R∕R0, where R0 accounts for the explosive energy only [NM + C4, 𝜀 = 0 in (5)]. c Peak blast overpressure decays with distance for H-10 aluminum particles assuming a full contribution of the Al chemical energy (𝜀 = 1). d Peak blast overpressure decays with distance for H-10 aluminum particles assuming that 25% (𝜀 = 0.25) of the Al chemical energy contributes to the blast front. e Peak blast overpressure decays with distance for H-10 aluminum particles assuming that 50% (𝜀 = 0.5) of the Al chemical energy contributes to the blast front is accounted for by the energy lost to heating and accelerating the particles, whereas the small overpressure excess in the far field, as noted earlier, is due to particle–gas interactions [10, 11]. This figure shows that the reduction in the blast performance for explosives with added inert particles is primarily due to the dilution effect. In Fig. 17a, the blast overpressure decay for the addition of H-10 Al particles is shown with respect to dimensional distance. In the remainder of the plots, the results are shown plotted versus non-dimensional distance, in which different values of 𝜀 were assumed in the definition of scaled distance, to determine if the overpressure results collapse. The working assumption is that if the overpressure results collapse in the scaled coordinates, at a given distance, then the value of 𝜀 chosen represents the effective contribution of the total energy release by the particle reaction to the blast front. In Fig. 17b, c, the scaled distance is chosen assuming, respectively, no aluminum energy release (𝜀 = 0) and that the full aluminum energy release contributes to the blast propagation ( 𝜀 = 1). In both cases, the curves generally fail to collapse, indicating that these assumptions are not consistent with the experimental results. The sole exception to that is the assumption of complete energy release for a solid mass fraction of 10%, which results in the overpressure curve nearly overlapping the baseline explosive case in the far field (Fig. 17c). This suggests, at this small mass fraction, that the energy release associated with complete 221 oxidation of the aluminum powder contributes to the peak blast overpressure by the time the blast reaches the far field. The complete reaction of the particles in the detonation products is theoretically possible at this mass fraction, as demonstrated earlier with CHEETAH calculations. Figure 17d shows the blast overpressure results plotted with a scaling assuming that 25% of the particle energy release contributes to the blast (i.e., 𝜀 = 0.25). The results collapse in the near field, at a distance near the edge of the fireball interface. At this distance, the particles will have reacted primarily with the detonation products, so one interpretation of the results is that one-half of aluminum particles fully react with water within the products and the energy release associated with this reaction contributes to the nearfield blast overpressure. In Fig. 17e, the results are plotted assuming that 50% of the particle energy release contributes to the blast strength. The results collapse in the far field, suggesting that additional energy is released, possibly due to reaction of the hydrogen produced, that contributes to strengthening the blast wave as it moves from the near field to the far-field region. A similar methodology is used to plot the results for H-50 and H-95 aluminum and for magnesium particles in Fig. 18, with values chosen for 𝜀 ranging from 0.15 to 0.5 for the aluminum particles and 0.3 to 0.65 for the magnesium particles. In the far field, for aluminum particles, with 𝜀 = 0.5, the curves for mass loadings of 40–50% approach that of the baseline explosive, suggesting that the fraction of the particle energy release that has contributed to the blast wave overpressure by the far-field approaches 50%. For smaller particle mass loadings (10–30%), the scaled overpressure values actually exceed the baseline values, suggesting that at these mass loadings even more than half the particle energy release contributes to the blast wave. Note again that the energy that contributes to the blast does not reflect the real mass of reactive particles that react. 13 222 Q. Pontalier et al. 300 300 40% H-50 (#7) 40% H-50 (#8) 30% H-50 (#6) 20% H-50 (#5) 10% H-50 (#4) Baseline (#0) 250 Peak overpressure (kPa) Peak overpressure (kPa) 250 200 150 100 50 0 40% H-50 (#7) 40% H-50 (#8) 30% H-50 (#6) 20% H-50 (#5) 10% H-50 (#4) Baseline (#0) 200 150 100 50 1 1.5 2 2.5 0 3 1 1.5 R/R0 (a) H50, ε = 0.15 Peak overpressure (kPa) Peak overpressure (kPa) 50% H-95 (#12) 40% H-95 (#11) 30% H-95 (#10) 10% H-95 (#9) Baseline (#0) 250 200 150 100 50 200 150 100 50 1 1.5 2 2.5 0 3 1 1.5 R/R0 2.5 3 (d) H-95, ε = 0.5 300 300 40% Mg (#16) 30% Mg (#15) 20% Mg (#14) 10% Mg (#13) Baseline (#0) 40% Mg (#16) 30% Mg (#15) 20% Mg (#14) 10% Mg (#13) Baseline (#0) 250 Peak overpressure (kPa) 250 Peak overpressure (kPa) 2 R/R0 (c) H-95, ε = 0.15 200 150 100 50 200 150 100 50 1 1.5 2 R/R0 (e) Mg, ε = 0.3 13 3 300 50% H-95 (#12) 40% H-95 (#11) 30% H-95 (#10) 10% H-95 (#9) Baseline (#0) 250 0 2.5 (b) H-50, ε = 0.5 300 0 2 R/R0 2.5 3 0 1 1.5 2 R/R0 (f) Mg, ε = 0.65 2.5 3 Blast enhancement from metalized explosives ◂Fig. 18 Same graphs as in Fig. 17 considering that a assuming that 15% (𝜀 = 0.15) of the H-50 Al chemical energy contributes to the blast front, b assuming that 50% (𝜀 = 0.5) of the H-50 Al chemical energy contributes to the blast front, c assuming that 15% (𝜀 = 0.15) of the H-95 Al chemical energy contributes to the blast front, d assuming that 50% (𝜀 = 0.5) of the H-95 Al chemical energy contributes to the blast front, e assuming that 30% (𝜀 = 0.3) of the Mg chemical energy contributes to the blast front, f assuming that 65% (𝜀 = 0.65) of the Mg chemical energy contributes to the blast front 6 Conclusions The effect of adding both inert (glass) and reactive (aluminum and magnesium) particles to gelled nitromethane on the performance of the blast wave generated by the detonation of the resulting heterogeneous explosive in constant volume charges was determined experimentally. The results were compared with the baseline case of the homogeneous explosive itself. Blast arrival times and peak blast overpressures in the near-to-mid field were analyzed from pressure gauge data. High-speed videography was used to infer the peak overpressure decays in the mid-to-far field. The main conclusions of this study are: 1. The addition of inert glass particles mitigates the peak blast overpressure by up to 38% in the mid field (i.e., just outside the maximum extent of the fireball), with the mitigation increasing with particle mass loading. When the peak blast overpressure is plotted as a function of non-dimensional distance accounting for the total energy release of the explosive, the blast decay results largely collapse, suggesting that the blast performance depends mainly on the explosive energy available. The addition of particles slightly attenuates the blast in the near field (due to heat and momentum transfer to the particles) and augments the blast in the far field (due to particle–flow interactions), consistent with previous results [10, 11]. 2. Based on the reduced arrival time of the blast wave at the nearest pressure gauge with Al or Mg particles, relative to that of the baseline explosive, the energy release from the particles influences the blast wave propagation within about 300 μs of charge detonation in all cases. The peak blast overpressure was found to be increased by up to almost twice the baseline value in the far field. The best blast performance with aluminum particles was observed at a particle mass loading between 30% and 50% depending on the distance from the charge. 3. Although the Al particle size was varied over an order of magnitude, the blast enhancement effect was rela- 223 tively insensitive to particle size. This suggests that the use of correlations for particle reaction time that depend monotonically on particle size is not sufficient to predict the blast performance. New particle reaction models are needed that are also consistent with the early metal reactions observed in metal acceleration experiments with metalized explosives. 4. By plotting the blast overpressure results versus an energy-scaled distance, and comparing the resulting overpressure decay with the baseline explosive, it is possible to gain insight into the degree to which the particle reaction contributes to the blast wave overpressure when the blast wave reaches various distances. The results for the metalized explosives containing aluminum or magnesium particles suggest that at a particle mass fraction of 10%, most of the energy release from the particle reaction has contributed to the blast propagation by the far field, which here is defined as distances greater than 3 m from 1 L volume charges. At particle mass fractions of 30–50%, more than half of the particle reaction energy contributes to the blast wave overpressure by the far field. The most probably reaction pathway for the aluminum in nitromethane consists of early submillisecond reaction of the aluminum with the water in the detonation products, releasing hydrogen, followed by combustion of the hydrogen with the surrounding air in the highly turbulent mixing zone near the fireball interface. Nevertheless, it is hypothesized that the possibility of the contribution to the blast of the early-time metal energy release would also occur in hydrogen-free compounds through other reaction pathways. The present results demonstrate that sub-millisecond energy release from the reaction of the metal particles with the water in the detonation products contributes to augmenting the blast strength. Additional experimental diagnostics very near the charge are needed to further resolve the rate of particle energy release at early times. Some insight into the partition of the particle energy release history between the detonics (microseconds), early expansion (hundreds of microseconds), and later fireball expansion (millisecond) stages of the phenomena may be obtained by comparing hydrocode predictions of the blast decay with experimental results, assuming a parametric variation in the particle energy release rate, and efforts in this direction are currently underway [50]. 13 224 Q. Pontalier et al. Table 3 Blast times of arrival for all trials at 1.0 m, 1.5 m, and 2.0 m from the charge, captured with pressure transducers At every distance, the arrival times are normalized by the baseline value (ta ∕ta-baseline) for comparison. Values of arrival times for two specific trials (trials #15 and #16), at 2.0 m, are not available due to record problems. At each distance, and for every type of particles, the most efficient trial at reducing (resp. augmenting) the blast arrival is in red (resp. blue) Table 4 Same as Table 3 but reordered for comparison at a specific mass fraction 13 Blast enhancement from metalized explosives 225 Table 5 Absolute peak overpressure (ΔPs) values extracted from the pressure gauge at 2.0 from the charge for all the trials tested Normalized peak overpressure values by the peak overpressure value of the baseline referred as normalized ΔPs ∕ΔPs-baseline are also calculated. For every type of particles, the most efficient trial at enhancing (resp. mitigating) the peak blast overpressure is in red (resp. blue). The peak overpressure obtained for the trial #4 (10% H-50) is not considered as reliable Table 6 Same as Table 5 but reordered for direct comparison at a specific particle mass fraction 13 226 Q. Pontalier et al. Table 7 Comparison between peak overpressures at 2.0 from the charge obtained through videography or with pressure gauge measurements extp: peak overpressure values extrapolated by inferring the blast curvature as shown in Fig. 5, sd: overpressure profile presenting a slight deviation from the modified Friedlander form, id: overpressure profile presenting an important deviation from the modified Friedlander form, R.E.: relative error between the two methods Appendix Tables 3 and 4 present the times of arrival captured with the pressure gauges at all distances from the charge (1.0 m, 1.5 m, and 2.0 m). Table 3 is sorted by particle size, whereas Table 4 is sorted by mass fraction. The same tables, but for the peak blast overpressures obtained at 2.0 m from the charge, are shown in Tables 5 and 6. These peak blast overpressures are compared with the peak overpressures obtained with the videography analysis at a distance of 2.0 from the charge in Table 7. Note that most of the time the blast front was not visible at this distance, at the same height as the charge. Consequently, the blast front position at 2.0 m from the charge was extrapolated by inferring the blast curvature as shown in Fig. 5. Discarding the non-reliable peak overpressure recorded for trial #4, the maximum discrepancy observed between the gauge measurements and 13 the video analysis method is 42.2%. However, the discrepancies between the two methods are not similar for all trials tested, suggesting that they are of a stochastic nature and not a systematic error resulting from the methodology employed. The discrepancies are explained by the fact that the blast is not tracked in the same directions around the charge (pressure gauges are not in the camera field of view as shown in Fig. 1). In other words, the blast propagation is not isotropic. Figure 19 exhibits the peak blast overpressure decays obtained with the videography in the mid-to-far field for trials with aluminum and magnesium particles. The peak overpressure decays are also plotted in Fig. 20a, c, e at three specific mass fractions (10%, 30%, and 40%) to illustrate the effect of the material and particle sizes. The normalized plots are also shown in Fig. 20b, d, f to highlight the performance. The trial videos corresponding to Figs. 7, 8, and 9 are available online as supplementary materials. Blast enhancement from metalized explosives 227 300 300 40% H-10 (#3) 30% H-10 (#2) 10% H-10 (#1) Baseline (#0) 250 Peak overpressure (kPa) Peak overpressure (kPa) 250 200 150 100 50 0 1.5 40% H-50 (#7) 40% H-50 (#8) 30% H-50 (#6) 20% H-50 (#5) 10% H-50 (#4) Baseline (#0) 200 150 100 50 2 2.5 3 3.5 4 4.5 0 1.5 5 2 2.5 Distance (m) (a) H-10 4 4.5 5 300 50% H-95 (#12) 40% H-95 (#11) 30% H-95 (#10) 10% H-95 (#9) Baseline (#0) 40% Mg (#16) 30% Mg (#15) 20% Mg (#14) 10% Mg (#13) Baseline (#0) 250 Peak overpressure (kPa) 250 Peak overpressure (kPa) 3.5 (b) H-50 300 200 150 100 50 0 1.5 3 Distance (m) 200 150 100 50 2 2.5 3 3.5 Distance (m) (c) H-95 4 4.5 5 0 1.5 2 2.5 3 3.5 4 4.5 5 Distance (m) (d) Mg Fig. 19 Peak blast overpressure decays plotted in the mid-to-far field for trials with a H-10 Al, b H-50 Al, c H-95 Al, and d Gran 17 Mg 13 228 Q. Pontalier et al. 300 2 Peak overpressure (kPa) 250 Normalized peak overpressure 10% H-10 (#1) 10% H-50 (#4) 10% H-95 (#9) 10% Mg (#13) Baseline (#0) 10% Glass (#17) 200 150 100 50 0 1.5 2 2.5 3 3.5 4 4.5 10% H-10 (#1) 10% H-50 (#4) 10% H-95 (#9) 10% Mg (#13) 10% Glass (#17) 1.5 1 0.5 1.5 5 2 2.5 Distance (m) (a) 10 wt% 300 Normalized peak overpressure Peak overpressure (kPa) 200 150 100 50 2 2.5 3 3.5 4 4.5 5 4 4.5 5 4 4.5 5 1 2 2.5 3 3.5 Distance (m) (d) 30 wt% 2 Normalized peak overpressure 40% H-10 (#3) 40% H-50 (#8) 40% H-50 (#7) 40% H-95 (#11) 40% Mg (#16) Baseline (#0) 250 Peak overpressure (kPa) 4.5 1.5 0.5 1.5 5 300 200 150 100 50 2 2.5 3 3.5 Distance (m) (e) 40 wt% 13 4 30% H-10 (#2) 30% H-50 (#6) 30% H-95 (#10) 30% Mg (#15) 30% Glass (#18) Distance (m) (c) 30 wt% 0 1.5 3.5 2 30% H-10 (#2) 30% H-50 (#6) 30% H-95 (#10) 30% Mg (#15) Baseline (#0) 30% Glass (#18) 250 0 1.5 3 Distance (m) (b) 10 wt% 4 4.5 5 40% H-10 (#3) 40% H-50 (#8) 40% H-50 (#7) 40% H-95 (#11) 40% Mg (#16) 1.5 1 0.5 1.5 2 2.5 3 3.5 Distance (m) (f) 40 wt% Blast enhancement from metalized explosives 229 ◂Fig. 20 Effect of particle material and particle sizes on the peak blast overpressure decays with dimensional distance, for particle mass fractions of 10 wt%, 30 wt%, and 40 wt% Acknowledgements The authors thank the field operation team at DRDC Suffield for assistance with setting up the experimental test site and S. Trebble for carrying out the high-speed video recording. The authors would like to thank A. Longbottom and A. Milne at Fluid Gravity Engineering for their valuable comments on the pressure gauge analysis. The authors would also like to thank the four anonymous reviewers for their many constructive comments. Compliance with ethical standards Conflict of Interest The authors declare that they have no conflict of interest. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 11. 12. 13. 14. 15. 16. 17. 18. 19. References 20. 1. Escales, R.: UK24377 (1899) 2. Roth, G.: Performance of explosives (1900) 3. Yen, N.H., Wang, L.Y.: Reactive metals in explosives. Propellants Explos. Pyrotech. 37(2), 143–155 (2012). https://doi.org/10.1002/ prep.200900050 4. Frost, D.L., Zhang, F.: The nature of heterogeneous blast explosives. 19th International Symposium on Military Aspects of Blast and Shock (MABS 19), Calgary, Canada (2006) 5. Trzciński, W.A., Maiz, L.: Thermobaric and enhanced blast explosives—properties and testing methods. Propellants Explos. Pyrotech. 40(5), 632–644 (2015). https://doi.org/10.1002/ prep.201400281 6. Türker, L.: Thermobaric and enhanced blast explosives (TBX and EBX). Def. Technol. 12(6), 423–445 (2016). https://doi.org/ 10.1016/j.dt.2016.09.002 7. Vadhe, P.P., Pawar, R.B., Sinha, R.K., Asthana, S.N., Subhananda Rao, A.: Cast aluminized explosives (review). Combust. Explos. Shock Waves 44(4), 461–477 (2008). https:// doi. org/ 10. 1007/ s10573-008-0073-2 8. Capellos, C., Baker, E.L., Nicolich, S., Balas, W., Pincay, J., Stiel, L.I.: Eigenvalue detonation of combined effects aluminized explosives. AIP Conf. Proc. 955(2007), 357–360. https://doi.org/10. 1063/1.2833055 9. Frost, D.L.: Heterogeneous/particle-laden blast waves. Shock Waves 28(3), 439–449 (2018). https:// doi. org/ 10. 1007/ s00193-018-0825-1 10. Pontalier, Q., Loiseau, J., Goroshin, S., Frost, D.L.: Experimental investigation of blast mitigation and particle-blast interaction during the explosive dispersal of particles and liquids. 21. 22. 23. 24. 25. Shock Waves 28(3), 489–511 (2018). https:// doi. org/ 10. 1007/ s00193-018-0821-5 Pontalier, Q., Lhoumeau, M., Milne, A., Longbottom, A., Frost, D.L.: Numerical investigation of particle-blast interaction during explosive dispersal of liquids and particles. Shock Waves 28(3), 513–531 (2018). https://doi.org/10.1007/s00193-018-0820-6 Pontalier, Q., Lhoumeau, M., Frost, D.L.: Blast wave mitigation in granular materials. AIP Conf. Proc. 1979(1), 110014 (2018). https://doi.org/10.1063/1.5044933 Loiseau, J., Pontalier, Q., Milne, A., Goroshin, S., Frost, D.L.: Terminal velocity of liquids and granular materials accelerated by a high explosive. Shock Waves 28(3), 473–487 (2018). https:// doi.org/10.1007/s00193-018-0822-4 Ripley, R .C., Zhang, F., Lien, F.-S.: Acceleration and heating of metal particles in condensed matter detonation. Proc. R. Soc. Lond. A 468(February), 1564–1590 (2012). https://doi.org/10. 1098/rspa.2011.0595 Peuker, J.M., Krier, H., Glumac, N.: Particle size and gas environment effects on blast and overpressure enhancement in aluminized explosives. Proc. Combust. Inst. 34(2), 2205–2212 (2013). https:// doi.org/10.1016/j.proci.2012.05.069 Zarei, Z., Frost, D.L.: Simplified modeling of blast waves from metalized heterogeneous explosives. Shock Waves 21(5), 425–438 (2011). https://doi.org/10.1007/s00193-011-0316-0 Balakrishnan, K., Menon, S.: On the role of ambient reactive particles in the mixing and afterburn behind explosive blast waves. Combust. Sci. Technol. 182(2), 186–214 (2010). https://doi.org/ 10.1080/00102200903341579 McNesby, K.L., Homan, B.E., Ritter, J.J., Quine, Z., Ehlers, R.Z., McAndrew, B.A.: Afterburn ignition delay and shock augmentation in fuel rich solid explosives. Propellants Explos. Pyrotech. 35(1), 57–65 (2010). https://doi.org/10.1002/prep.200800084 Duan, J., Wu, C., An, F.J., Liao, S.S.: A quasi-analytical model of aluminized explosive products under strong constraint. Adv. Mech. Eng. 9(11), 1–10 (2017). https://doi.org/10.1177/16878 14017739511 Yuan, X., Wu, C., An, F., Liao, S., Zhou, M., Shi, Z., Xue, D.: A non-isentropic model of aluminized explosives involved with the reaction degree of aluminum powder for post-detonation burning behavior. Phys. Fluids 32, 023306 (2020). https://doi.org/10. 1063/1.5139692 Li, X., Pei, H., Zhang, X., Zheng, X.: Effect of aluminum particle size on the performance of aluminized explosives. Propellants Explos. Pyrotech. 45(5), 807–813 (2020). https://doi.org/10.1002/ prep.201900308 Lefrancois, A., Baudin, G., Le Gallic, C., Boyce, P., Coudoing, J-P.: Nanometric aluminum powder influence on the detonation efficiency of explosives. 12th International Detonation Symposium, San Diego, U.S.A. (2002) Manner, V.W., Pemberton, S.J., Gunderson, J.A., Herrera, T.J., Lloyd, J.M., Salazar, P.J., Rae, P., Tappan, B.C.: The role of aluminum in the detonation and post-detonation expansion of selected cast HMX-based explosives. Propellants Explos., Pyrotech. 37(2), 198–206 (2012). https://doi.org/10.1002/prep.20110 0138 Makhov, M.N., Gogulya, M.F., Dolgoborodov, A.Yu., Brazhnikov, M.A., Arkhipov, V.I., Pepekin, V.I.: Acceleration ability and heat of explosive decomposition of aluminized explosives. Combust. Explos. Shock Waves 40(4), 458–466 (2004). https://doi.org/10. 1023/B:CESW.0000033569.77449.d9 Tappan, B.C., Hill, L.G., Manner, V.W., Pemberton, S.J., Lieber, M.A., Johnson, C.E., Sanders, V.E.: Reactions of powdered aluminum with explosives that selectively form carbon dioxide or water as oxidizers. Int. J. Energ. Mater. Chem. Propuls. 15(4), 339–350 (2016). https:// doi. org/ 10. 1615/ IntJE nergeticMa teria lsChemProp.2017011503 13 230 26. Robbins, D.L., Anderson, E.K., Anderson, M.U., Jackson, S.I., Short, M.: Cylinder test characterization of an ammonium nitrate and aluminum powder explosive. 15th International Symposium on Detonation, vol. 803, pp. 826–835 (2014) 27. Baudin, G., Lefrancois, A., Bergues, D., Bigot, J., Champion, Y.: Combustion of nanophase aluminum in the detonation products of nitromethane and TNT. 11th International Detonation Symposium, Colorado, U.S.A. (1998) 28. Loiseau, J., Goroshin, S., Frost, D.L., Higgins, A.J., Zhang, F.: Ability of metalized gelled nitromethane to accelerate a flyer plate. 16th International Detonation Symposium, Cambridge, U.S.A. (2018) 29. Kato, Y., Murata, K., Itoh, S.: Detonation characteristics of packed beds of aluminum saturated with nitromethane. 13th International Detonation Symposium, IDS 2006, pp. 187–195 (2006) 30. Schmitt, M.M., Bowden, P.R., Tappan, B.C., Henneke, D.: Steadystate shock-driven reactions in mixtures of nano-sized aluminum and dilute hydrogen peroxide. J. Energ. Mater. 36(3), 266–277 (2018). https://doi.org/10.1080/07370652.2017.1376233 31. Loiseau, J., Chase, G., Goroshin, S., Frost, D.L.: Detonation and metal acceleration of aluminum–water mixtures. 16th International Detonation Symposium, Cambridge, U.S.A. (2018) 32. Beckstead, M.W.: Correlating aluminum burning times. Combust. Explos. Shock Waves 41(5), 533–546 (2005). https://doi.org/10. 1007/s10573-005-0067-2 33. Brousseau, P., Dorsett, H.E., Cliff, M.D., Anderson, C.J.: Detonation properties of explosives containing nanometric aluminum powder. 12th Internation Detonation Symposium, pp. 1–10, San Diego, U.S.A. (2002) 34. Maiz, L., Trzciński, W.A.: Detonation and combustion of new heterogeneous composite explosives containing aluminum particles. 11th International Armament Conference on Scientific Aspects of Armament and Safety Technology, Ryn, Poland (2016) 35. Frost, D.L., Zhang, F., Murray, S., McCahan, S.: Critical conditions for ignition of metal particles in a condensed explosive. 12th International Detonation Symposium, pp. 1–10 (2002) 36. Trzciński, W.A., Paszula, J., Grys, S.: Detonation and blast wave characteristics of nitromethane mixed with particles of an aluminium–magnesium alloy. Propellants Explos. Pyrotech. 35(2), 85–92 (2010). https://doi.org/10.1002/prep.200900041 37. Dewey, J.M.: The shape of the blast wave: studies of the Friedlander equation. 21st International Symposium on Military Aspects of Blast and Shock (MABS 21), pp. 1–9, Jerusalem, Israel (2010) 38. Kleine, H., Dewey, J.M., Ohashi, K., Mizukaki, T., Takayama, K.: Studies of the TNT equivalence of silver azide charges. 13 Q. Pontalier et al. 39. 40. 41. 42. 43. 44. 45. 46. 47. 48. 49. 50. Shock Waves 13(2), 123–138 (2003). https:// doi. org/ 10. 1007/ s00193-003-0204-3 Fried, L.E., Howard, W.M., Souers, P.C.: Cheetah 2.0 user’s manual. Technical report, Lawrence Livermore National Laboratory, Livermore, CA, USA (1998) Keshavarz, M.H.: Detonation temperature of high explosives from structural parameters. J. Hazard. Mater. 137(3), 1303–1308 (2006). https://doi.org/10.1016/j.jhazmat.2006.04.057 Tanguay, V., Goroshin, S., Higgins, A.J., Zhang, F.: Aluminum particle combustion in high-speed detonation products. Combust. Sci. Technol. 181(4), 670–693 (2009). https://doi.org/10.1080/ 00102200802643430 Weingärtner, H., Franck, E.U.: Supercritical water as a solvent. Angew. Chem. Int. Ed. 44(18), 2672–2692 (2005). https://doi.org/ 10.1002/anie.200462468 Holley, C.E., Huber, E.J.: The heats of combustion of magnesium and aluminum. J. Am. Chem. Soc. 73(12), 5577–5579 (1951). https://doi.org/10.1021/ja01156a020 Persson, P.-A., Holmberg, R., Lee, J.: Rock Blasting and Explosives Engineering. CRC Press, Boca Raton (1993) Baker, W.E.: Explosions in Air. University of Texas Press, Austin (1973) Kinney, G.F., Graham, K.J.: Explosive Shocks in Air, 2nd edn. Springer, Berlin (1985) Bach, G.G., Lee, J.H.S.: An analytical solution for blast waves. AIAA J. 8(2), 271–275 (1970). https://doi.org/10.2514/3.5655 Dobratz, B.M., Crawford, P.C.: LLNL explosives handbook— properties of chemical explosives and explosive simulants. Technical report, Lawrence Livermore National Laboratory (LLNL), Livermore, California (1985) Leadbetter, J., Ripley, R.C., Zhang, F., Frost, D.L.: Multiple energy scaling of blast waves from heterogeneous explosives. 21st International Colloquium on the Dynamics of Explosions and Reactive Systems (ICDERS), Poitiers, France (2007) Pontalier, Q., Longbottom, A., Cargill, S., Grapes, C., Loiseau, J., Goroshin, S., Frost, D.L., Zhang, F.: Numerical investigation of blast enhancement with metalized explosives. 32nd International Symposium on Shock Waves, pp. 1–13, Singapore (2019). https:// doi.org/10.3850/978-981-11-2730-4_0323-cd Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.