SEARCH · Engineering Papers
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.
On February 15, 2013 approximately 17 m asteroid hit Earth, causing shock waves and air blasts over a portion of Russia. A significant amount of material has been recovered from this meteorite fall, officially named Chelyabinsk.
Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.
The detection, location, and monitoring of underground structures are of great importance to national and global security. Tunnels and voids generate seismic signatures detectable at the surface, but using non-invasive seismic data to image near-surface presents several challenges in real-world applications. In this report, we describe the use of a dense surface seismic deployment to generate subsurface models of the Graymont Pleasant Gap mine - a single-layer mine with a complex structure embedded in a high-velocity P-wave limestone bedrock. Our approach consists of three key methods. We use P-wave arrival times from local blast events to perform a tomography inversion with the tomoTD method, constructing a P-wave velocity model of the subsurface. We model the layer above the mine using Rayleigh wave ellipticity and inversion techniques. We leverage ongoing anthropogenic activities to identify and locate noise sources both on the surface and within the subsurface. With this integrated approach we aim to overcome the challenges and enhance our ability to non-invasively characterize underground structures, contributing to improved seismic monitoring techniques.
Objective/Hypothesis We examined a neuroinflammatory response associated with glial activation in the cochlea exposed to blast overpressure and evaluated the potential therapeutic efficacy of specialized pro‐resolving mediators such as neuroprotectin D1, NPD1; (10R, 17S‐dihydroxy‐4Z, 7Z, 11E, 13E, 15Z, 19Z‐docosahexaenoic acid) in a rodent blast‐induced auditory injury model. Study Design Animal Research. Methods A compressed‐air driven shock tube was used to expose anesthetized adult male Long‐Evan rats to shock waves simulating an open‐field blast exposure. Approximately 30 minutes after blast exposure, rats were treated with NPD1 (100 ng/kg body wt.) or vehicle delivered intravenously via tail vein injection. Rats were then euthanized 48 hours after blast exposure. Unexposed rats were included as controls. Tissue sections containing both middle and inner ear were prepared with hematoxylin–eosin staining to elucidate histopathological changes associated with blast exposure. Cochlear tissues were evaluated for relative expression of ionized calcium‐binding adaptor 1 (Iba1), as an indicator of microglial activation by immunohistochemistry and western blot analyses. Results Our animal model resulted in an acute injury mechanism manifested by damage to the tympanic membrane, hemorrhage, infiltration of inflammatory cells, and increased expression of Iba1 protein. Moreover, therapeutic intervention with NPD1 significantly reduced Iba1 expression in the cochlea, suggesting a reduction of a neuroinflammatory response caused by blast overpressure. Conclusions Blast overpressure resulted in an increased expression of proteins involved in gliosis within the auditory system, which were reduced by NPD1. Treatment of NPD1 suggests an effective strategy to reduce or halt auditory microglial cell activation due to primary blast exposure. Level of Evidence NA Laryngoscope , 131:E2018–E2025, 2021
The 15 January 2022 Hunga, Tonga, volcano's explosive eruption produced the most powerful blast recorded in the last century, with an estimated equivalent TNT yield of 100–200 megatons. The blast energy was propagated through the atmosphere as various wave types. The most prominent wave was a long-period (>2000 s) surface-guided Lamb wave with energy comparable to that of the 1883 Krakatoa Lamb wave; both were clearly observed by pressure sensors (barometers) worldwide. Internal gravity, acoustic-gravity, and infrasound waves were captured in great detail by the entire infrasound component of the International Monitoring System (IMS). For instance, infrasound waves (<300 s period) were seen to circumnavigate Earth up to eight times. Atmospheric waves captured by the IMS infrasound network and selected barometers near the source provide insight on Earth's impulse response at planetary scales.
Shock waves are supersonic high-amplitude pressure waves that cause barotrauma when they transfer kinetic energy to the tissues of animals.1–4 Snapping shrimp (Alpheidae) produce shock waves and are exposed to them frequently, so we asked if these animals have evolved mechanisms of physical protection against them. Snapping shrimp generate shock waves by closing their snapping claws rapidly enough to form cavitation bubbles that release energy as an audible ‘‘snap’’ and a shock wave when they collapse.5–8 We tested if snapping shrimp are protected from shock waves by a helmet-like extension of their exoskeleton termed the orbital hood. Using behavioral trials, we found shock wave exposure slowed shelterseeking and caused a loss of motor control in Alpheus heterochaelis from which we had removed orbital hoods but did not significantly affect behavior in shrimp with unaltered orbital hoods. Shock waves thus have the potential to harm snapping shrimp but may not do so under natural conditions because of protection provided to shrimp by their orbital hoods. Using pressure recordings, we discovered the orbital hoods of A. heterochaelis dampen shock waves. Sealing the anterior openings of orbital hoods diminished how much they altered the magnitudes of shock waves, which suggests these helmet-like structures dampen shock waves by trapping and expelling water so that kinetic energy is redirected and released away from the heads of shrimp. Our results indicate orbital hoods mitigate blast-induced neurotrauma in snapping shrimp by dampening shock waves, making them the first biological armor system known to have such a function.
Shock wave diffraction by two dimensional weak disturbances, using aerodynamics blast method
Acoustic metamaterials have proven to be a versatile tool for the precise control and manipulation of sound waves. One of the promising designs of acoustic metamaterials employ the arrays of bubbles and find applications for soundproofing, blast mitigation, and many others. An obvious advantage of bubble-based metamaterials is their ability to be relatively thin while absorbing low-frequency sound waves. The vast majority of theories developed to predict resonant behavior of bubble-based metamaterials capitalize on Minnaert frequency. Here, we propose a novel theoretical approach to characterize bubble-based metamaterials that are based on our previous findings for a single bubble trapped in circular cavity modeled as a thin clamped plate. We obtain analytical expressions for resonant frequencies of bubble metascreens using self-consistent approximation. Two geometry factors, distance between bubble centers and distance between bubble center and interface of acoustic impedance change, are taken into account. We demonstrate the existence of multiple bandgaps and possibility of switching between them via adjustment of geometry parameters and reflector properties.
Blast of gas blows a loop of unexposed film as a wave across a vacuum platen to feed film smoothly into a camera so that 2 successive lengths can be exposed within 50 milliseconds. This technique can be readily applied to multiple aperture cameras as well as to various types of films.
Sources such as atmospheric or buried explosions and shallow earthquakes are known to produce infrasonic pressure waves in the atmosphere. Because of the coupling between neutral particles and electrons at ionospheric altitudes, these acoustic and gravity waves induce variations of the ionospheric electron density. The Global Positioning System (GPS) provides a way of directly measuring the total electron content in the ionosphere and, therefore, of detecting such perturbations in the upper atmosphere. In July and August 1996, three large surface mine blasts (1.5 Kt each) were detonated at the Black Thunder coal mine in eastern Wyoming. As part of a seismic and acoustic monitoring- experiment, we deployed five dual-frequency GPS receivers at distances ranging from 50 to 200 km from the mine and were able to detect the ionospheric perturbation caused by the blasts. The perturbation starts 10 to 15 min after the blast, lasts for about 30 min, and propagates with an apparent horizontal velocity of 1200 meters per second. Its amplitude reaches 3 x 10 (exp 14) el per square meters in the 7-3 min period band, a value close to the ionospheric perturbation caused by the M = 6.7 Northridge earthquake. The small signal-to-noise ratio of the perturbation can be improved by slant-stacking the electron content time-series recorded by the different GPS receivers taking into account the horizontal propagation of the perturbation. The energy of the perturbation is concentrated in the 200 to 300 second period band, a result consistent with previous observations and numerical model predictions. The 300 second band probably corresponds to gravity modes and shorter periods to acoustic modes, respectively. Using a 1-D stratified velocity model of the atmosphere we show that linear acoustic ray tracing fits arrival times at all GPS receivers. We interpret the perturbation as a direct acoustic wave caused by the explosion itself. This study shows that even relatively small subsurface events can produce ionospheric perturbations that are above the detection threshold of the GPS technique. By sensing derivative signals, which can be detected over a relatively broad region, it appears that GPS might be particularly useful for source characterization within the first acoustic quiet zone where infrasound would probably be ineffective. This suggests that dual-frequency GPS monitoring could contribute to Comprehensive Nuclear Test Ban Treaty verification.
Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.
The extinction mechanisms proposed at the Cretaceous-Tertiary geological boundary are reviewed and related to the impact of asteroids or comets in general. For impact energies below 10(exp 4) Megatons (less than 6 x 10(exp 4) years; asteroid diameter less than 650 m), blast, earthquake, and fire may destroy local areas up to 10(exp 5) square m. Tidal waves could flood a kilometer inland over entire ocean basins. The energy range from 105 to 106 Megatons (less than 2 x 10(exp 6) years; asteroid diameter less than 3 km) is transitional. Dust lifted, sulfur released from within impacting asteroids, and soot from fires started by comets can produce climatologically significant optical depths of 10. At energies beyond 10(exp 7) Megatons, blast and earthquake damage is regional (10(exp 6) square cm). Tsunami cresting to 100 m and flooding 20 km inland will sweep the coastal zones of the world's oceans. Fires will be set globally. Light levels may drop so low from the smoke, dust and sulfate that vision is not possible. At energies approaching 10(exp 9) Megatons the ocean surface waters may be acidified by sulfur. The combination of these effects would be devastating.
As a chemical explosion is buried, the mechanism for acoustic wave generation transitions from fully gas-generated at the surface to completely spall-induced at full containment depth. The fully gas-generated and completely spall-induced signals in the acoustic waveform are well described; however, the transition between these two end-members eludes numerical modeling because of the complex phenomena that are involved. The phenomena of crater formation and explosive cloud evolution are simulated using an Eulerian hydrocode that incorporates geomaterials with strength and porosity. Having accurately modeled these phenomena, we can confidently predict the propagation and relative strength of the gas-generated and spall-induced pulses in the recorded acoustic waveform. The numerical predictions agree with observations from the historical Stagecoach experiment as well as modern recordings from the Source Physics Experiment. In particular, the peak pressure p generated by an explosion is initially due to the gas-generated mechanism and decays with scaled depth of burial d s (depth d scaled by the cube-root of explosive yield w 1/3 ) as exp(-d s ) but then transitions near a scaled depth of 6 m/ton 1/3 to the spall-generated mechanism in which the decay is d$_{s}^{-7/4}$. This decay form is related to the strong ground-motion attenuation relationship that affects spall strength. So these results can improve seismoacoustic inverse models for the explosive source that need to account for the gas-generated and spall-induced signals and their effect on peak pressures and other acoustic signal features.
A series of large damaging explosions involving fireworks storage occurred near Esparto, Yolo County California on July 1, 2025. Three explosions were located and reported by the University of California Berkeley Seismology Laboratory (UCB/BSL) and the United States Geological Survey (USGS). Analysis of local distance (< 30 km) seismic recordings of the first blast indicates that there were actually three explosions with later blasts delayed by about 3 and 30 seconds. We measured the first arriving P-wave amplitudes on four of these events with good signal-to-noise ratios.
Following impact, vapor expansion creates an intense airblast that interacts with the ambient atmosphere. The resulting hemi-spherical shock wave leaves a signature on the surface that is dependent on initial atmospheric and surface conditions. Here we propose that the formation of pedestal craters (craters surrounded by an erosion-resistant pedestal) may be a direct consequence of extreme winds and elevated temperatures generated by such an impact-induced atmospheric blast. Pedestal craters, first recognized in Mariner 9 data, are a unique feature on Mars and likely a signature of near-surface volatiles. They are found at high latitudes (small pedestals, Amazonian to Late Hesperian in age) and in thick equatorial mantling deposits (larger pedestals, early Hesperian to Noachian in age). Previously suggested mechanisms for pedestal crater formation (e.g., wind: ejecta curtain vortices or vapor blast; and ejecta dust: armoring) do not provide a complete picture. The clear evidence for near-surface volatiles at high latitudes requires a re-evaluation of these alternative models. The results presented here suggest that a combined atmospheric blast/thermal model provides a plausible formation hypothesis.
A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.
Continuous Reflectometry for Radius versus Time Experiments, or CORRTEX, is a diagnostic that was developed during nuclear testing to measure the distance at which rock walls were disrupted by an explosive device as a function of the time at which disruption occurred. Distances recorded are dynamic measurements of cable length as the shockwaves from an explosive event crush the experiment cables, reducing their usable length. Applications of the CORRTEX diagnostic have been expanded in recent times to include instrumentation of the explosive device itself, to record phenomena such as the speed of the detonation wave in the explosive charge and the crushing or disassembly of material in the charge housing. The diagnostic has likewise been used to monitor explosive performance in rock blasting and oil well hole-clearing activities.
Airburst - In the simulations explored energy from the airburst couples very weakly with the water making tsunami dangerous over a shorter distance than the blast for asteroid sizes up to the maximum expected size that will still airburst (approx.250MT). Future areas of investigation: - Low entry angle airbursts create more cylindrical blasts and might couple more efficiently - Bursts very close to the ground will increase coupling - Inclusion of thermosphere (>80km altitude) may show some plume collapse effects over a large area although with much less pressure center dot Ocean Impact - Asteroid creates large cavity in ocean. Cavity backfills creating central jet. Oscillation between the cavity and jet sends out tsunami wave packet. - For deep ocean impact waves are deep water waves (Phase speed = 2x Group speed) - If the tsunami propagation and inundation calculations are correct for the small (<250MT) asteroids in these simulations where they impact deep ocean basins, the resulting tsunami is not a significant hazard unless particularly close to vulnerable communities. Future work: - Shallow ocean impact. - Effect of continental shelf and beach profiles - Tsunami vs. blast damage radii for impacts close to populated areas - Larger asteroids below presumed threshold of global effects (Ø200 - 800m).