Arxiv:2010.13537V2 [Physics.Pop-Ph] 25 May 2021

Noname manuscript No. (will be inserted by the editor) Beirut explosion TNT equivalence from the fireball evolution in the first 170 milliseconds Charles J.Aouad1 · Wissam Chemissany2 · Paolo Mazzali1;3 · Yehia Temsah 4 · Ali Jahami 4 Received: date / Accepted: date Abstract The evolution of the fireball resulting from within which the fire ball is used to measure the shock the August 2020 Beirut explosion is traced using am- wave are in excellent agreement with other experimen- ateur footage taken during the first 400 ms after the tal data. A total TNT equivalence of 200 ± 80 tons at a detonation. 39 frames separated by 16.66 - 33.33 ms distance d ≈ 130 m is found. Finally, the dimensions of are extracted from six different videos located precisely the crater size taken from a hydro-graphic survey con- on the map. Measurements of the time evolution of the ducted 6 days after the explosion are scaled with known radius of the shock wave are traced by the fireball at correlation equations yielding a close range of results. consecutive time sequence until about t ≈ 170 ms and A recent published article by Dewey [12] shows that the a distance d ≈ 128 meters. Pixel scales for the videos Beirut explosion TNT equivalence is an increasing func- are calibrated by de-projecting the existing grains sition of distance. The results of the current paper are los building for which accurate as built drawings are quantitatively in excellent agreement with this finding. available using the length, the width and the height These results present an argument that the actual mass and by defining the line of sight incident angles. In the of ammonium nitrate that participated in the detona- distance range d ≈ 60-128 m from the explosion cen- tion is much less than the quantity that was officially ter, the evolution of the fireball follows a Sedov-Taylor claimed available. model with spherical geometry and an almost instanta- Keywords Explosion physics · Beirut explosion · neous energy release. This model is used to derive the Fireball · Shock wave · Blast · ammonium nitrate energy available to drive the shock front at early times. explosion Additionally, a drag model is fitted to the fireball evolution until its stopping distance at time t ≈ 500 ms at distance d ≈ 145± 5 m. Using the derived TNT equiv- 1 Introduction alent yield, the scaled stopping distance reached by the fireball and the shock wave-fireball detachment epoch On the 4th of August 2020, an explosion occurred in the port of Beirut, Lebanon, after a fire ignited in ware- Charles Aouad house number 12. This tragic event resulted in massive large scale destruction, severe damage to buildings in E-mail: [email protected] an extended radius around the center and loss of lives. 1 Astrophysics Research Institute, Liverpool John Moores arXiv:2010.13537v2 [physics.pop-ph] 25 May 2021 It was claimed by officials that this hangar contained University, IC2, Liverpool Science Park, 146 Brownlow Hill, Liverpool L3 5RF, UK an amount of 2750 tons of ammonium nitrate kept in · 2Institute for Quantum Information and Matter, California the port for around 6 years. Institute of Technology, 1200 E California Blvd, Pasadena, A few attempts have already been made to quantify CA 91125, USA the amplitude of this explosion. Several of these studies · 3 Max-Planck Institut fur Astrophysik, Karl-Schwarzschild- Str. 1, D-85748 Garching, Germany use the time of arrival of the shock wave (from audio · 4 Faculty of engineering, Beirut Arab University, Beirut, and visual inspection of footage) up to distances rang- Lebanon. ing from 500 to about 2000 m [38, 43, 35]. These studies use empirical relations linking the scaled time of arrival 2 Charles J.Aouad1 et al. with the scaled distance to report a TNT equivalence [47, 31]. These energetic burning explosive phenomena range of 0.3-1.1 kt. Pilger et al. [36] use open access are observed in nuclear and in chemical reactions on seismic data to yield a range of 0.5-1 kt TNT. a wide range of magnitudes ranging from small con- Diaz [14] measures the evolution of the fireball until trolled industrial activities to large astrophysical con- about 200 m distance from the center. He uses 26 data texts (e.g. solar flares, supernovae). It is usually ac- points taken from publicly available videos to yield a cepted that a deflagration can transit to a detonation range of 0.5-0.6 kt of TNT. The data of Rigby et al. [38] in suitable conditions [28, 23, 34], although this tran- is used by Dewey [12] to compare the peak hydrostatic sition is still an extensive theoretical area of research. overpressures with experimental measurements for both Once the shock wave reaches the outer boundary of the TNT and ANFO explosions. He finds that Beirut explo- burning fuel, it will be transmitted to the surrounding sion produces overpressures that are weaker than over- medium (whether a fluid or a solid) and will propagate pessures that would have been produced by a TNT ex- isotropically in the form of a blast wave. plosion of the same energy at close distances to the cen- This sudden increase in pressure will cause the ter, while at larger distances it produces slightly larger ignited hot material and the gaseous residues of the re- over-pressures. He concludes that the TNT equivalence action to expand rapidly with high velocities pushing on of the Beirut explosion is an increasing function of dis- the fluid around them [22, 25]. Furthermore, the sudden tance (refer to Fig 3 in [12]). increase in pressures at the vicinity of the detonation In this paper we report the TNT equivalence by will cause an increase in temperature (few thousands measuring the kinematics of the fireball in the close degree kelvin). Finally the transfer of momentum be- proximity of the center (in a distance range of 60-145 tween the shock front and solid particles will also cause m from the explosion center in the first 170 ms). We these solid particles to accelerate spherically away from base our study on both experimental observations of the center [48, 19, 39]. The combined effect will result the fireball evolution generated by chemical explosions in the creation of an optically thick visible fireball. Im- as reported by Gordon, Gross, and Perram [17], and mediately after the explosion, both the fireball and the using a Sedov-Taylor model to derive our results [45, shock front rapidly expand. However, the expansion of 46, 40]. the fireball will decelerate until it halts and reaches a The paper is divided as follows: in section2, we stopping distance while the shock front will detach and present our methodology, in section3, we present our keep expanding depositing energy until it decays into a results and observations and we discuss these results. sonic wave. As the rate of fireball expansion decreases, In section4, we compare our work to the literature and both temperatures and atmospheric over-pressures will in section5 we present some additional thoughts. We also drop making its sharp edge less defined and less draw our conclusions and discuss possible lines of future luminous [see also4, sect303-304] [17,8, 46, 40]. work in section6. 2.2 Mathematical description 2 Methodology The theoretical models for the study of the shock 2.1 Background wave generation and propagation were developed by studying nuclear explosions [8, 40, 45]. Taylor [45] con- Explosions are the swift release of a large amount sidered the solution where the energy is released in- of energy [24, 27,6]. This process is usually caused by stantaneously, in a very small volume (point source) the ignition of a `fuel'. A burning front is then formed and where the mass of the explosives is insignificant: and propagates within the medium burning it as it pro- i.e. nuclear explosion. He then derived the evolution of ceeds. The power of this process: i.e. the energy release the shock wave by solving numerically three differential per unit time, depends not only on the chemical or nu- equations, namely an equation of motion, a continu- clear potential of the fuel, but also on the velocity of ity equation and an equation of gas state taking the propagation of this burning front throughout the ma- boundary conditions given by the Rankine-Hugoniot terial. In a pure deflagration, this burning front propa- equations [37, 18]. Later, he was able to experiment gates subsonically and causes burning by heat transfer the accuracy of his work after the publication of the [20, 33]. photographs of the Trinity nuclear explosion test [44, A more powerful form of burning is a detona- 29]. Although Taylor analysis describes best a nuclear tion, in which case, the burning front travels super- explosion, he later provided evidence, based on exper- sonically creating a shock front ahead of it and can imental data, and showed the boundaries (scaled dis- cause the burning of the fuel by compressive heating tance) where a chemical explosion can resemble a nu- Beirut explosion 3 clear explosion and where his model can still be applied [46]. The Taylor model has the form: 5 2 −1 −1 R = t EK ρ0 : (1) Here ρ0 is the undisturbed gas density, K is a constant that depends on γ the ratio of specific heats of the gas and E the part of the energy that has not been radiated away. This relation can be written as 5 log R = a log t + b (2) 2 10 10 where a is the slope of this linear relation and is ex- pected to be equal to unity if the observation follows the theoretical prediction.

Load more