Engineering Papers⌕ Search

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.

At least 145 records · Page 8

Towards the usage of quarry blasts as an active seismic source for subsurface characterization and imaging

SUMMARY Active or man-made seismic sources are often used to delineate subsurface geological structures via seismic imaging techniques. Although conventional controlled seismic sources have provided high-resolution subsurface images, the high cost of data acquisition necessitates improved use of cheaper alternative seismic sources for subsurface imaging, such as traffic noise, quarry blasts and earthquakes. Here, we explore the potential use of quarry blasts as an active seismic source for subsurface imaging. The goal of our study is to demonstrate the utility of quarry blasts as seismic sources for investigating subsurface structure with a field experiment in Mill Creek, southern Oklahoma. We deployed a 7455-m-long receiver line of 72 receivers (105 m spacing) extending away from a quarry, with a recording time of 28 d. We analyse similarities between blasts, wave types, origin time and excitation location of quarry blasts, which is important information for using such signal as an active source. Given a relative homogeneous geological condition in the study area, we could not observe clear reflected waves. However, the use of P- and S-wave velocities led to the delineation of a major basement damage zone in the study area, which further demonstrates the viability of quarry blasts as a useful seismic source for subsurface geological imaging.

58 GEOSCIENCES↗

SOMD LOX/Methane Assessment, Data Analysis Subteam: MOX-E Test Series Analysis Summary

Fundamental explosive testing was conducted by the NASA/MSFC Blast Team (N2L) to determine yields produced from cryogenic liquid mixtures of methane and oxygen for a range of oxygen-to-fuel ratios. The blast wave overpressure and high-speed video data were processed by the NASA led multi-agency data analysis team. This document serves to provide a summary of the optimized data analysis produced and the resulting conclusions regarding the yields achieved during this test program.

Combustion↗

Interstellar Dust: Physical Processes

Dust is formed in stellar environments, and destroyed by sputtering, shattering and vaporization in shock waves due to cloud-cloud collisions and supernova blast waves. Dust is also destroyed during star formation. We review the dust formation and destruction balance. The calculated destruction time-scale is less than or equal to one billion years and the star dust injection time-scale is approx. 2.5 billion years. Hence, the fractions of elemental carbon and silicon locked up in stardust are less than 0.3 and less than 0.15, respectively. An efficient ISM dust formation route is therefore implied. In particular, in dense clouds dust grows; through the processes of coagulation and the accretion of gas phase molecules e.g. H20, CO, CH4. These icy materials may then be photoprocessed to refractory materials in more diffuse regions. The resulting carbonaceous grain mantle may actually be the glue that holds the coagulated grains together.

Jones, A. P.↗

Exhaust Nozzle for a Multitube Detonative Combustion Engine

An improved type of exhaust nozzle has been invented to help optimize the performances of multitube detonative combustion engines. The invention is applicable to both air-breathing and rocket engines used to propel some aircraft and spacecraft, respectively. In a detonative combustion engine, thrust is generated through the expulsion of combustion products from a detonation process in which combustion takes place in a reaction zone coupled to a shock wave. The combustion releases energy to sustain the shock wave, while the shock wave enhances the combustion in the reaction zone. The coupled shockwave/reaction zone, commonly referred to as a detonation, propagates through the reactants at very high speed . typically of the order of several thousands of feet per second (of the order of 1 km/s). The very high speed of the detonation forces combustion to occur very rapidly, thereby contributing to high thermodynamic efficiency. A detonative combustion engine of the type to which the present invention applies includes multiple parallel cylindrical combustion tubes, each closed at the front end and open at the rear end. Each tube is filled with a fuel/oxidizer mixture, and then a detonation wave is initiated at the closed end. The wave propagates rapidly through the fuel/oxidizer mixture, producing very high pressure due to the rapid combustion. The high pressure acting on the closed end of the tube contributes to forward thrust. When the detonation wave reaches the open end of the tube, it produces a blast wave, behind which the high-pressure combustion products are expelled from the tube. The process of filling each combustion tube with a detonable fuel/oxidizer mixture and then producing a detonation repeated rapidly to obtain repeated pulses of thrust. Moreover, the multiple combustion tubes are filled and fired in a repeating sequence. Hence, the pressure at the outlet of each combustion tube varies cyclically. A nozzle of the present invention channels the expansion of the pulsed combustion gases from the multiple combustion tubes into a common exhaust stream, in such a manner as to enhance performance in two ways: (1) It reduces the cyclic variations of pressure at the outlets of the combustion tubes so as to keep the pressure approximately constant near the optimum level needed for filling the tubes, regardless of atmospheric pressure at the altitude of operation; and (2) It maximizes the transfer of momentum from the exhaust gas to the engine, thereby maximizing thrust. The figure depicts a typical engine equipped with a nozzle according to the invention. The nozzle includes an interface section comprising multiple intake ports that couple the outlets of the combustion tubes to a common plenum. Proceeding from its upstream to its downstream end, the interface section tapers to a larger cross-sectional area for flow. This taper fosters expansion of the exhaust gases flowing from the outlets of the combustion tubes and contributes to the desired equalization of exhaust combustion pressure. The cross-sectional area for flow in the common plenum is greater than, or at least equal to, the combined cross-sectional flow areas of the combustor tubes. In the common plenum, the exhaust streams from the individual combustion tubes mix to form a single compound subsonic exhaust stream. Downstream of the common plenum is the throat that tapers to a smaller flow cross section. In this throat, the exhaust gases become compressed to form a compound sonic gas stream. Downstream of the throat is an expansion section, which typically has a bell or a conical shape. (The expansion section can be truncated or even eliminated in the case of an air-breathing engine.) After entering the expansion section, the exhaust gases expand rapidly from compound sonic to compound supersonic speeds and are then vented to the environment. The basic invention admits of numerous variations. For example, the combustion tubes can be arranged around the central axin a symmetrical or asymmetrical pattern other than the one shown in the figure. For another example, the flow cross-sectional area(s) of one or more of the intake ports in the interface section, of the common plenum, the throat, and/or the expansion section can be varied, either symmetrically or asymmetrically, to adjust dynamics of the exhaust stream or to direct the thrust vector away from the central axis.

Bratkovich, Thomas E.↗

Effects of High Fuel Loading and CO2 Dilution on Oxy-Methane Ignition Inside a Shock Tube at High Pressure

Abstract Ignition delay times were measured for methane/O2 mixtures in a high dilution environment of either CO2 or N2 using a shock tube facility. Experiments were performed between 1044 K and 1356 K at pressures near 16 ± 2 atm. Test mixtures had an equivalence ratio of 1.0 with 16.67% CH4, 33.33% O2, and 50% diluent. Ignition delay times were measured using OH* emission and pressure time-histories. Data were compared to the predictions of two literature kinetic mechanisms (ARAMCO MECH 2.0 and GRI Mech 3.0). Most experiments showed inhomogeneous (mild) ignition which was deduced from five time-of-arrival pressure transducers placed along the driven section of the shock tube. Further analysis included determination of blast wave velocities and locations away from the end wall of initial detonations. Blast velocities were 60–80% of CJ-Detonation calculations. A narrow high temperature region within the range was identified as showing homogenous (strong) ignition which showed generally good agreement with model predictions. Model comparisons with mild ignition cases should not be used to further refine kinetic mechanisms, though at these conditions, insight was gained into various ignition behavior. To the best of our knowledge, we present first shock tube data during ignition of high fuel loading CH4/O2 mixtures diluted with CO2 and N2.

Energy & Fuels↗

Acceleration of positrons in supernova shocks

During this project we investigated the acceleration of leptons (electrons and positrons) in collisionless shock waves. In particular, we were interested in how leptons are accelerated in the blast waves existing in the remnants of supernova explosions. Supernova remnants (SNRs) have long been considered as the most likely source of galactic cosmic rays but no definite connection between SNRs and the cosmic rays seen at earth can be made. Only by understanding lepton acceleration in shocks can the rich SNR data base be properly used to understand cosmic ray origins. Our project was directed at the neglected aspects of lepton acceleration. We showed that the efficiency of lepton acceleration depended critically on the lepton injection energy. We showed that, even when infection effects are not important, that proton and lepton distribution functions produced by shocks are quite different in the critical energy range for producing the observed synchrotron emission. We also showed that transrelativistic effects produced proton spectra that were not in agreement with standard results from radio observations, but that the lepton spectra were, in fact, consistent with observations. We performed simulations of relativistic shocks (shocks where the flow speed is a sizable fraction of the speed of light) and discovered some interesting effects. We first demonstrated the power of the Monte Carlo technique by determining the shock jump conditions in relativistic shocks. We then proceeded to determine how relativistic shocks accelerate particles. We found that nonlinear relativistic shocks treat protons and leptons even more differently than nonrelativistic shocks. The transrelativistic effects on the shock structure from the heavy ion component reduces the lepton acceleration to a tiny fraction of the ion acceleration. This effect is dramatic even if high energy leptons (many times thermal energy) are injected, and was totally unexpected. Our results have important consequences for astrophysical environments expected to harbor relativistic flows such as extra-galactic radio sources and accretion onto compact objects.

Ellison, Donald C.↗

Uniform high order spectral methods for one and two dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗

Uniform high-order spectral methods for one- and two-dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one-dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two-dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗

Surface Pressures and Heat Transfer on Unswept Blunt Plates in Helium at High Mach Numbers

Pressure distributions and local convective heat-transfer coefficients on a flat plate at zero angle of attack were measured in helium. Data were obtained with various amounts of leading-edge bluntness at Mach numbers of 12.5 and 14.7. The pressures on a sharp leading-edged plate were not influenced by the leading edge and were predicted by the first-order, hypersonic, weak-interaction theory. Pressures on blunt plates were correlated by introducing the leading-edge Reynolds number as a parameter. Measured heat-transfer coefficients on the sharp plate agreed with predictions obtained form existing exact solutions for hear transfer across the laminar boundary layer. For the blunt plates a comparison of theory with experiment indicated that more knowledge of the flow field between the sock wave and plate surface is necessary before an adequate prediction of convective heat transfer can be made. Shock-wave shapes for the blun plates at a Mach number 12.5 and zero angle of attack were measured. At distances between 2 and 60 leading-edge thicknesses from the shock vertex, the shock-wave shapes were found to be represented by a modified form of the blast-wave analogy.

Marvin, Joseph G.↗

Space-Time Conservation Element and Solution Element Method Being Developed

The engineering research and design requirements of today pose great computer-simulation challenges to engineers and scientists who are called on to analyze phenomena in continuum mechanics. The future will bring even more daunting challenges, when increasingly complex phenomena must be analyzed with increased accuracy. Traditionally used numerical simulation methods have evolved to their present state by repeated incremental extensions to broaden their scope. They are reaching the limits of their applicability and will need to be radically revised, at the very least, to meet future simulation challenges. At the NASA Lewis Research Center, researchers have been developing a new numerical framework for solving conservation laws in continuum mechanics, namely, the Space-Time Conservation Element and Solution Element Method, or the CE/SE method. This method has been built from fundamentals and is not a modification of any previously existing method. It has been designed with generality, simplicity, robustness, and accuracy as cornerstones. The CE/SE method has thus far been applied in the fields of computational fluid dynamics, computational aeroacoustics, and computational electromagnetics. Computer programs based on the CE/SE method have been developed for calculating flows in one, two, and three spatial dimensions. Results have been obtained for numerous problems and phenomena, including various shock-tube problems, ZND detonation waves, an implosion and explosion problem, shocks over a forward-facing step, a blast wave discharging from a nozzle, various acoustic waves, and shock/acoustic-wave interactions. The method can clearly resolve shock/acoustic-wave interactions, wherein the difference of the magnitude between the acoustic wave and shock could be up to six orders. In two-dimensional flows, the reflected shock is as crisp as the leading shock. CE/SE schemes are currently being used for advanced applications to jet and fan noise prediction and to chemically reacting flows.

Chang, Sin-Chung↗

Complex structure of the supernova remnant HB 3

HB 3 is an old, large (84 pc diameter) supernova remnant associated with the W3 H II region/molecular cloud complex. Observations of the imaging proportional counter (IPC) onboard the Einstein X-ray astronomy satellite have been reprocessed to yield a contour map of X-ray brightness and spectra of various regions of this remnant. The measured IPC flux is 2.4 x 10 to the -11th ergs per sq cm per s, giving a 0.2-4 keV luminosity of 1.6 x 10 to the 35th ergs/s for a column densityof 6 x 10 to the 21st per sq cm. The measured X-ray temperatures reveal a decrease from center to limb of the remnant of 1-0.3 keV. HB 3 is in the late adiabatic blast-wave phase of evolution, 30,000 to 50,000 yr old and with an initial blast energy of 3 x 10 to the 50th ergs. The X-ray map is compared with available radio and optical images. In X-rays, HB 3 has two components - a diffuse emission inside the 84 pc radio remnant and a ring of emission at the center of 30 pc in diameter. The diffuse emission is similar to that from other supernova remnants which are moderately obscured (column density, nH approximately 10 to the 22nd per sq cm). Three possibilities for the origin of the ring are explored: (1) a second supernova remnant, (2) a shocked shell in the interstellar medium surrounding HB 3, and (3) reverse-shock heated ejecta. There is no hot neutron star within the remnant.

Leahy, D. A.↗

Investigation of Spherical-Wave-Initiated Flow Fields Around Bodies

Measurements of the velocity flow fields and vortex movements have been made about various simple blunt models undergoing spherical blast waves with a positive overpressure of 4 pounds per square inch. A bullet-optical method was used to determine flow velocities and is applied to velocity fields in which the gradients are largely normal to the free-stream direction. The velocity flow fields are shown at various flow times following passage of the blast front for different models. Vortex movements with time are compared for square-bar models of various aspect ratios. Corner sharpness had no discernible effect on the overall disturbed velocity fields or vortex movements for the square-box models used.

McFarland, Donald R.↗

Review of the PDWA Concept for Combustion Enhancement in a Supersonic Air-Breathing Combustor Environment

This paper reviews the design of the Pulsed Detonation Wave Augmentor (PDWA) concept and the preliminary computational fluid dynamics studies that supported it. The PDWA relies on the rapid generation of detonation waves in a small tube, which are then injected into the supersonic stream of the main combustor. The blast waves thus generated are used to stimulate the mixing and combustion inside the main combustor. The mixing enhancement relies on various forms of the baroclinic interaction, where misaligned pressure and density gradients combine to produce vortical flow. By using unsteady shock waves, the concept also uses the Richtmyer-Meshkov effect to further increase the rate of mixing. By carefully designing the respective configurations of the combustor and the detonation tubes, one can also increase the penetration of the fuel into the supersonic air stream. The unsteady shocks produce lower stagnation pressure losses than steady shocks. Combustion enhancement can also be obtained through the transient shock-heating of the fuel-air interface, and the lowering of the ignition delay in these regions. The numerical simulations identify these processes, and show which configurations give the best results. Engineering considerations are also presented, and discuss the feasibility of the concept. Of primary importance are the enhancements in performance, the design simplicity, the minimization of the power, cost, and weight, and the methods to achieve very rapid cycling.

Canbier, Jean-Luc↗

CFD Applications in Support of the Space Shuttle Risk Assessment

The paper describes a numerical study of a potential accident scenario of the space shuttle, operating at the same flight conditions as flight 51L, the Challenger accident. The interest in performing this simulation is derived by evidence that indicates that the event itself did not exert large enough blast loading on the shuttle to break it apart. Rather, the quasi-steady aerodynamic loading on the damaged, unbalance vehicle caused the break-up. Despite the enormous explosive potential of the shuttle total fuel load (both liquid and solid), the post accident explosives working group estimated the maximum energy involvement to be equivalent to about five hundreds of pounds of TNT. This understanding motivated the simulation described here. To err on the conservative side, we modeled the event as an explosion, and used the maximum energy estimate. We modeled the transient detonation of a 500 lbs spherical charge of TNT, placed at the main engine, and the resulting blast wave propagation about the complete stack. Tracking of peak pressures and impulses at hundreds of locations on the vehicle surface indicate that the blast load was insufficient to break the vehicle, hence demonstrating likely crew survivability through such an event.

Baum, Joseph D.↗

On the Exit Boundary Condition for One-Dimensional Calculations of Pulsed Detonation Engine Performance

In one-dimensional calculations of pulsed detonation engine (PDE) performance, the exit boundary condition is frequently taken to be a constant static pressure. In reality, for an isolated detonation tube, after the detonation wave arrives at the exit plane, there will be a region of high pressure, which will gradually return to ambient pressure as an almost spherical shock wave expands away from the exit, and weakens. Initially, the flow is supersonic, unaffected by external pressure, but later becomes subsonic. Previous authors have accounted for this situation either by assuming the subsonic pressure decay to be a relaxation phenomenon, or by running a two-dimensional calculation first, including a domain external to the detonation tube, and using the resulting exit pressure temporal distribution as the boundary condition for one-dimensional calculations. These calculations show that the increased pressure does affect the PDE performance. In the present work, a simple model of the exit process is used to estimate the pressure decay time. The planar shock wave emerging from the tube is assumed to transform into a spherical shock wave. The initial strength of the spherical shock wave is determined from comparison with experimental results. Its subsequent propagation, and resulting pressure at the tube exit, is given by a numerical blast wave calculation. The model agrees reasonably well with other, limited, results. Finally, the model was used as the exit boundary condition for a one-dimensional calculation of PDE performance to obtain the thrust wall pressure for a hydrogen-air detonation in tubes of length to diameter ratio (L/D) of 4, and 10, as well as for the original, constant pressure boundary condition. The modified boundary condition had no performance impact for values of L/D > 10, and moderate impact for L/D = 4.

Wilson, Jack↗

Numerical Simulation of Bolide Entry with Ground Footprint Prediction

As they decelerate through the atmosphere, meteors deposit mass, momentum and energy into the surrounding air at tremendous rates. Trauma from the entry of such bolides produces strong blast waves that can propagate hundreds of kilometers and cause substantial terrestrial damage even when no ground impact occurs. We present a new simulation technique for airburst blast prediction using a fully-conservative, Cartesian mesh, finite-volume solver and investigate the ability of this method to model far- field propagation over hundreds of kilometers. The work develops mathematical models for the deposition of mass, momentum and energy into the atmosphere and presents verification and validation through canonical problems and the comparison of surface overpressures, and blast arrival times with actual results in the literature for known bolides. The discussion also examines the effects of various approximations to the physics of bolide entry that can substantially decrease the computational expense of these simulations. We present parametric studies to quantify the influence of entry-angle, burst-height and other parameters on the ground footprint of the airburst, and these values are related to predictions from analytic and handbook-methods.

Simulation↗

X-ray structure of the Cygnus loop

An X-ray surface brightness map of the Cygnus Loop in the energy range 0.15-0.85 keV is presented. The X-ray distribution is compared with the optical and radio pictures of the Loop. The location of the center of the supernova shell as derived from the X-ray data is in excellent agreement with the apparent optical and radio centers. We obtain a mean diameter for the X-ray shell of 2.75 (+ or - 0.2) deg. The X-ray spectral data yield a characteristic thermal temperature averaged over the Loop of 2.9 million (+ or - 1.5 million) K. We have calculated anew X-ray parameters for adiabatic blast waves using X-ray production spectra which include line radiation. These parameters are compared with the results of the present experiment. An age of 18,000 years, initial blast energy equal to 3 times 10 to the 50-th power ergs, and ambient interstellar density of 0.25 per cu cm are found.

Rappaport, S.↗

Evolution of a superbubble blastwave in a magnetized medium

Researchers investigate the effects of interstellar magnetic fields on the evolution and structure of interstellar superbubbles, using both analytic and numerical magnetohydrodynamic (MHD) calculations. These cavities of hot gas, surrounded by shells of cold dense material preceded by a shock wave result from the combined action of stellar winds and supernova explosions in OB associations. If the medium in which a superbubble goes off is homogeneous and unmagnetized, the blast wave expands isotropically. As the interstellar gas flows through the shock, it cools significantly and gets strongly compressed such that thermal pressure remains approximately equal to ram pressure. Hence, the swept up material is confined to a very thin shell. However, if the ambient medium is permeated by a uniform magnetic field B sub o approx. 3 mu G (typical value for the interstellar matter (ISM)), the configuration loses its spherical symmetry, and, due to magnetic pressure, the shell of swept up material does not remain thin. Researchers found the following qualitative differences: (1) Except in the immediate vicinity of the magnetic poles, the shell is supported by magnetic pressure. (2) The refraction of field lines at the shock and the thermal pressure gradient along the shell both contribute to accelerating the gas toward the equator. The resulting mass flux considerably decreases the column density at the magnetic poles. (3) Away from the poles, magnetic tension in the shell causes the field lines (particularly the inner boundary) to elongate in the direction of B sub o. In contrast, the shock wave radius increases with increasing theta. (4) The reduced inertia of a parcel in the polar neighborhood makes it easier to decelerate, and accounts for the dimple which appears at the poles in numerical simulations. This dimple also results from the necessity to call on intermediate shocks in order to insure a smooth transition between a purely thermal shock at the poles and a magnetic shock in the rest of the shell. (5) The shock wave propagates faster than in the absence of magnetic field, except near the poles where the reduced mass of the shell allows it to be more efficiently decelerated.

Ferriere, Katia M.↗