Calculating Spherical Blast Wave Parameters From Height-Of-Burst (HOB) Explosions
Explore the source record for details and available documents.
John Sykes Richardson↗
SEARCH · Engineering Papers
Results for “Blast waves”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Computation of blast wave-obstacle interactions
Numerical simulations of the interaction of a planar blast wave with various obstacles are presented. These obstacles are either ground structures or vehicles flying in the atmosphere. For a structure on the ground, the blast wave encounter is side-on, while for the flying vehicles the encounter is either head-on or oblique. Second-order accurate, finite-difference, and shock-capturing procedures are employed to solve the two-dimensional, axisymmetric, and three-dimensional unsteady Euler equations. Results are presented for the flow field consisting of blast wave striking obstacles that are at rest, moving subsonically and moving supersonically. Comparison of the numerical results with experimental data for a configuration at rest substantiates the validity of this approach and its potential as a flow analysis tool.
Champney, J. M.↗
Characteristics of pressure waves
Air blast characteristics generated by most types of explosions are discussed. Data cover both negative and positive blast load phases and net transverse pressure as a function of time. The effects of partial or total confinement, atmospheric propagation, absorption of energy by ground shock or cratering, and transmission over irregular terrain on blast wave properties were also considered.
Source record↗
An X-ray and optical study of the interaction of the Cygnus Loop supernova remnant with an interstellar cloud
We have used the ROSAT high-resolution imager and optical emission line data to study the structure of a bright optical knot on the southeastern rim of the Cygnus Loop supernova remnant. This knot has been identified as an encounter between the blast wave and a small isolated cloud. The knot appears in projection just behind the blast wave, which is traced by Balmer line filaments which bound the X-ray emission. The knot is a prominent X-ray feature, consisting of a numerb of filaments which are correlated with the optical line emission. These data permit a detailed view of the blast wave interaction. By combining the optical and X-ray data it is possible to trace the blast wave as a continuous surface from its southern edge around the western side of the cloud that then continues on to the north. The southeastern knot is an indentation on the surface of the blast wave. Thus the southeastern knot is not a small cloud that has been overrun by the blast wave, but the tip of a larger cloud. Bright X-ray emission is associated with bright radiative filaments. The location of this emission, upstream of the radiative shocks, implies that the enhanced X-rays come from a reverse shock. The presence of a reverse shock is further evidence that the southeastern knot represents an early stage of a blast wave encountering a large cloud. A consistent picture of the Cygnus Loop as an explosion within a preexisting cavity is emerging. These observations agree with that picture.
Graham, James, R.↗
Two-dimensional Thomson scattering measurements of misaligned electron density and temperature gradients and associated Biermann battery produced fields
We present optical Thomson scattering measurements of electron density and temperature in high Mach number laser-driven blast waves in homogeneous gases. Taylor–Sedov blast waves are launched in nitrogen (N2) or helium (He) at pressures between 0.4 mTorr and 10 Torr by ablating a solid plastic target with a high energy laser pulse (10 J, 1012 W/cm2). Experiments are performed at high repetition rate (1 Hz), which allows one-dimensional and two-dimensional Thomson scattering measurements over an area of several cm2 by automatically translating the scattering volume between shots. Electron temperature and density in the blast wave fronts were seen to increase with increasing background gas pressure. Measured electron density and temperature gradients were used to calculate ∂B/∂t∝∇Te×∇ne. The experimentally measured ∂B/∂t showed agreement with the magnetic field probe (B-dot) measurements, revealing that magnetic fields are generated in the observed blast waves via the Biermann battery effect. The results are compared to numerical three-dimensional collisional magnetohydrodynamic simulations performed with FLASH, and are discussed in the context of spontaneous magnetic field generation via the Biermann battery effect.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Chemical Blast Standard (1 kg)
Chemical explosions create blast waves with large overpressure disturbances. It is important to develop a standard blast model based on data to accurately predict acoustic blast-wave amplitudes near detonations and invert for explosion energy from distant observations of blast-wave signals. However, open data from large, controlled chemical explosions with reliable ground truth can be challenging to find. The lack of access to such data could limit the number of contributions to related research and potentially stifle the rate of discoveries or validation of existing models. Here, to address these data scarcity problem, we have curated and compiled a standardized set of 817 blast-wave waveforms from 19 distinct high-explosive events. The blast-wave waveforms are standardized to a 1 kg trinitrotoluene explosion using scaling laws and corrections for location effects. A brief overview of the dataset is presented along with explosion feature models as well as recommendations for extracting explosion features. The resulting dataset is distributed to an open repository in both Seismic Analysis Code and pandas DataFrame formats containing the waveforms, the scaled distances, and the sample rates.
58 GEOSCIENCES↗
Consequences of impacts of small asteroids and comets with Earth
The fragmentation of a small asteroid in the atmosphere greatly increases its cross sections for aerodynamic braking and energy dissipation. At a typical impact velocity of 22 km/s, the atmosphere absorbs more than half the kinetic energy of stony meteoroids with diameters, D(sub m), less than 220 m and iron meteoroids with D(sub m) less than 80 m. The corresponding diameter for comets with impact velocity 50 km/s is D(sub m) less than 1600 m. Most of the atmospheric energy dissipation occurs in a fraction of a scale height, so large meteors appear to 'explode' or 'flare' at the end of their visible paths. This dissipation of energy in the atmosphere protects the earth from direct impact damage (e.g., craters), but it produces a blast wave that can do considerable damage. The area of destruction around the impact point in which the over-pressure in the blast wave exceeds 4 lb/sq in = 2.8 x 10(exp 5) dynes/cu cm, which is enough to knock over trees and destroy buildings, increases rapidly from zero for chondritic meteoroids less than 56 m in diameter (15 megatons) to about 200 sq km for those 80 m in diameter (48 megatons); the probable diameter of the tunguska impactor of 1908 is about 80 m. Crater formation and earthquakes are not significant in land impacts by stony asteroids less than about 200 m in diameter because of the air protection. A tsunami is probably the most devastating type of damage for asteroids 200 m to 1 km in diameter. An impact by an asteroid this size anywhere in the Atlantic would devastate coastal areas on both sides of the ocean. An asteroid a few kilometers across would produce a tsunami that would reach the foothills of the Appalachian Mountains in the upper half of the East Coast of the United States. Most of Florida is protected from a tsunami by the gradual slope of the ocean off its coast, which causes most of the tsunami energy to be reflected back into the Atlantic. The atmosphere plume produced by asteroids with diameters exceeding about 120 m cannot be contained by the atmosphere, so this bubble of high-temperature gas forms a new layer on top of the atmosphere. The dust entrapped in this hot gas is likely to have optical depths exceeding tau = 10 for asteroids with diameters exceeding about 0.5 to 1 km. The optical flux from asteroids 60 m or more in diameter is enough to ignite pine forests. However, the blast wave from an impacting asteroid goes beyond the radius in which the fire starts. The blast wave tends to blow out the fire, so it is likely that the impact will char the forest, as at Tunguska, but the impact will not produce a sustained fire. Because comets dissipate their energy much higher in the atmosphere than asteroids, they illuminate a much larger region and their blast wave is weaker. So they are much more effective in producing large fires. This suggests that the KT impactor was a comet rather than an asteroid.
Hills, J. G.↗
Moreton Waves
"Moreton waves," named for the observer who popularized them, are a solar phenomenon also known in scientific literature as "Moreton-Ramsey wave," "flare waves," "flare-associated waves," "MHD blast waves," "chromospheric shock fronts" and various other combinations of terms which connote violently propagating impulsive disturbances. It is unclear whether all of the observations to which these terms have been applied pertain to a single physical phenomenon: there has perhaps been some overlap between the observations and the assumed physical properties of the observed occurrence. Moreton waves are ideally observed in the wings of H alpha, and appear as semi-circular fronts propagating at speeds ranging from several hundred to over a thousand km/sec. They form an arc, or "brow shape" which can span up to 180 degrees. Extrapolating the speed and locations of the arc indicates that the phenomenon's origin intersects well with the impulsive phase of the associated H alpha flare (if the flare exhibits an impulsive phase). However, the arc may not form or may not be observable until it is tens of megameters from the flaring region, and subsequently can propagate to distances exceeding 100 megameters. The high speeds and distances of propagation, plus the associated radio and energetic particle observations, provided strong evidence of a coronal, rather than a chromospheric origin. The H alpha manifestation of the wave is assumed to be the "ground track" or "skirt" of a three-dimensional disturbance.
Thompson, B. J.↗
Navier-Stokes simulations of unsteady transonic flow phenomena
Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.
Atwood, C. A.↗
Nonuniform Expansion of the Youngest Galactic Supernova Remnant G1.9+0.3
We report measurements of the X-ray expansion of the youngest Galactic supernova remnant, G1.9+0.3, using Chandra observations in 2007, 2009, and 2011. The measured rates strongly deviate from uniform expansion, decreasing radially by about 60 along the X-ray bright SE-NW axis from 0.84 plus or minus 0.06% yr(exp -1) to 0.52% plus or minus 0.03 yr(exp -1). This corresponds to undecelerated ages of 120-190 yr, confirming the young age of G1.9+0.3 and implying a significant deceleration of the blast wave. The synchrotron-dominated X-ray emission brightens at a rate of 1.9% plus or minus 0.4% yr(exp -1). We identify bright outer and inner rims with the blast wave and reverse shock, respectively. Sharp density gradients in either the ejecta or ambient medium are required to produce the sudden deceleration of the reverse shock or the blast wave implied by the large spread in expansion ages. The blast wave could have been decelerated recently by an encounter with a modest density discontinuity in the ambient medium, such as may be found at a wind termination shock, requiring strong mass loss in the progenitor.
NONUNIFORM EXPANSION↗
Shocks propagate in a 2D dusty plasma with less attenuation than due to gas friction alone
In a dusty plasma, an impulsively generated shock, i.e., blast wave, was observed to decay less than would be expected due to gas friction alone. In the experiment, a single layer of microparticles was levitated in a radio frequency glow-discharge plasma. In this layer, the microparticles were self-organized as a 2D solid-like strongly coupled plasma, which was perturbed by the piston-like mechanical movement of a wire. To excite a blast wave, the wire's motion was abruptly stopped, so that the input of mechanical energy ceased at a known time. It was seen that, as it propagated across the layer, the blast wave's amplitude persisted with little decay. This result extends similar findings, in previous experiments with 3D microparticle clouds, to the case of 2D clouds. In our cloud, out-of-plane displacements were observed, lending support to the possibility that an instability, driven by wakes in the ion flow, provides energy that sustains the blast wave's amplitude despite the presence of gas damping.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Nonequilibrium kinetic freeze-out properties in relativistic heavy ion collisions from energies employed at the RHIC beam energy scan to those available at the LHC
In this paper, we investigate the kinetic freeze-out properties in relativistic heavy ion collisions at different collision energies. We present a study of standard Boltzmann-Gibbs blast-wave (BGBW) fits and Tsallis blast-wave (TBW) fits performed on the transverse momentum spectra of identified hadrons produced in Au + Au collisions at collision energies of $\sqrt{s_{\rm{NN}}}=$7.7 –200 GeV at the Relativistic Heavy Ion Collider (RHIC), and in Pb + Pb collisions at collision energies of $\sqrt{s_{\rm{NN}}}=$ 2.76 and 5.02 TeV at the Large Hadron Collider (LHC). The behavior of strange and multistrange particles is also investigated. We found that the TBW model describes data better than the BGBW one overall, and the contrast is more prominent as the collision energy increases as the degree of nonequilibrium of the produced system is found to increase. From TBW fits, the kinetic freeze-out temperature at the same centrality shows a weak dependence of collision energy between 7.7 and 39 GeV, while it decreases as collision energy continues to increase up to 5.02 TeV. The radial flow is found to be consistent with zero in peripheral collisions at RHIC energies but sizable at LHC energies and central collisions at all RHIC energies. We also observed that the strange hadrons, with higher temperature and similar radial flow, approach equilibrium more quickly from peripheral to central collisions than light hadrons. The dependence of temperature and flow velocity on nonequilibrium parameter (q - 1) is characterized by two second-order polynomials. Both a and dξ from the polynomials fit, related to the influence of the system bulk viscosity, increase toward lower RHIC energies.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
SN 1987A after 18 Years: Mid-Infrared GEMINI and SPITZER Observations of the Remnant
We present high resolution 11.7 and 18.3 micron mid-IR images of SN 1987A obtained on day 6526 since the explosion with the Thermal-Region Camera and Spectrograph (T-ReCS) attached to the Gemini South 8m telescope. The 11.7 micron flux has increased significantly since our last observations on day 6067. The images clearly show that all the emission arises from the equatorial ring (ER). Nearly contemporaneous spectra obtained on day 6184 with the MIPS at 24 micron, on day 6130 with the IRAC in 3.6- 8 micron region, and on day 6190 with the IRS in the 12-37 micron instruments on board the Spitzer Space Telescope's show that the emission consists of thermal emission from silicate dust that condensed out in the red giant wind of the progenitor star. The dust temperature is 1662(sup +18) (sub -12) K, and the emitting dust mass is (2.6(sup +2.0 (sub -1.4)) x 10 (exp -6) M(solar). Lines of [Ne II] 12.82 micron and [Ne III] 15.56 pm are clearly present in the Spitzer spectrum, as well as a weak [Si II] 3 34.8 micron line. We also detect two lines near 26 micron which we tentatively ascribe to [Fe II] 25.99 pm and [0 IV] 25.91 micron. Comparison of the mid-IR Gemini 11.7 micron image with X-ray images obtained by Chandra, UV-optical images obtained by HST, and radio synchrotron images obtained by the ATCA show generally good correlation of the images across all wavelengths. Because of the limited resolution of the mid-IR images we cannot uniquely determine the location. or heating mechanism of the dust giving rise to the emission. The dust could be collisionally heated by the X-ray emitting plasma, providing a unique diagnostic of plasma conditions. Alternatively, the dust could be radiatively heated in the dense UV-optical knots that are overrun by the advancing supernova blast wave. In either case the dust-to-gas mass ratio in the circumstellar medium around the supernova is significantly lower than that in the general interstellar medium of the LMC, suggesting either a low condensation efficiency in the wind of the progenitor star, or the efficient destruction of the dust by the SN blast wave. Overall, we are witnessing the interaction of the SN blast wave with its surrounding medium, creating an environment that is rapidly evolving at all wavelengths. Continuous multiwavelength observations of SN 1987A such as these provide unique snapshots of the very early evolution of supernova remnants.
Bouchet, Patrice↗
An Empirical Non-TNT Approach to Launch Vehicle Explosion Modeling
In an effort to increase crew survivability from catastrophic explosions of Launch Vehicles (LV), a study was conducted to determine the best method for predicting LV explosion environments in the near field. After reviewing such methods as TNT equivalence, Vapor Cloud Explosion (VCE) theory, and Computational Fluid Dynamics (CFD), it was determined that the best approach for this study was to assemble all available empirical data from full scale launch vehicle explosion tests and accidents. Approximately 25 accidents or full-scale tests were found that had some amount of measured blast wave, thermal, or fragment explosion environment characteristics. Blast wave overpressure was found to be much lower in the near field than predicted by most TNT equivalence methods. Additionally, fragments tended to be larger, fewer, and slower than expected if the driving force was from a high explosive type event. In light of these discoveries, a simple model for cryogenic rocket explosions is presented. Predictions from this model encompass all known applicable full scale launch vehicle explosion data. Finally, a brief description of on-going analysis and testing to further refine the launch vehicle explosion environment is discussed.
Blackwood, James M.↗
Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai
The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.
54 ENVIRONMENTAL SCIENCES↗
Introduction to gasdynamics of explosions
Questions concerning the genesis and sustenance of an explosion are investigated, giving attention to the mechanics of explosions, the gasdynamics of explosions, aspects of technological significance, and future prospects. The dynamics of exothermic processes is discussed together with the most prominent effects of explosions. Blast waves are considered, taking into account conservation principles, blast wave transformation, conservative equations in nondimensional form, the equation of state, Eulerian space profiles, Eulerian time profiles, Lagrangian time profiles, boundary conditions and integral relations, and self-similar flow fields.
Oppenheim, A. K.↗
The Impact of SN1987A with its Circumstellar Ring
Following is final report on the study "The Impact of SN1987A with its Circumstellar Ring", which is now complete. In 1994, it was predicted that the blast wave from SN1987A would strike the circumstellar ring in AD1999+/-3, and that the ring would brighten by several hundreds optically as the shock entered the ring. It was also predicted that the emission lines from the shocked ring would have linewidths of a few hundred km/s. Today, we see a "hot spot" on the ring that first appeared in Hubble Space Telescope (HST) images in 1995 and has doubled in brightness between August 1997 and February 1998. Moreover, spectra from STIS show that the emission lines from the hot spot have widths and blue shifts of order 200 km/s, just as we predicted in 1994. Our guess of 1999+/-3 for the impact time was lucky, because we assumed that the hot gas in the bubble between the supernova and the ring had a low density (approximately 10 cm(exp -3)). But it was pointed out that the ROSAT observations of soft X-rays from SN 1987A implied that the intervening gas had a higher density, (approximately 100 cm(exp -3)), which would delay the impact until ca. 2007. Others developed hydrodynamic models to fit the ROSAT X-ray emission spectrum and came to the same conclusion. But these models were oversimplified in that they assumed that the ring was perfectly round. Now we see clearly that the hot spot is a peninsula that protrudes inward from the ring -- the first spot on the ring to be struck by the blast wave. No doubt there are other protrusions on the ring, which we may expect to light up in the next few years until they finally merge to set the entire ring ablaze, probably within 5 - 10 years. It was also predicted that the X-ray emission must be accompanied by optical and ultraviolet emission from atoms in the supernova debris and in the circumstellar gas that cross the reverse shock and the blast wave, respectively, and that the Lyman-alpha and NV emission lines should be bright enough to see with the STIS.
McCray, Richard↗
A sheet description of the emission from middle-aged supernova remnants
The hypothesis that sheetlike shock fronts propagating into large clouds immersed in a fairly uniform intercloud medium can provide a consistent picture accounting for essentially all the observational data for the Cygnus Loop SNR is considered. Filaments are interpreted as line-of-sight projections of portions of the sheet. Variations in the projection of the shock velocity along the line of sight lead naturally to 'turbulent' velocities consistent with observation. An analog model of a large, distorted, bumpy shock front is presented and shown to produce the observed morphology of the remnant's prominent optical filaments. Comparison with the observed morphology supports the idea that regions of optical emission lie on the local outer boundary of the shell and have not been engulfed by the blast wave. Bright X-ray emission arises when the blast wave encounters a preshock medium with an intermediate density higher than average for the adiabatic portion of the remnant, but still lower than that required for the shock to become radiative.
Hester, J. Jeff↗