Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “detonation wave”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Observation of breakdown wave mechanism in avalanche ionization produced atmospheric plasma generated by a picosecond CO 2 laser

Understanding the formation and long-timescale evolution of atmospheric plasmas produced by ultrashort, long-wavelength infrared (LWIR) pulses is an important but partially understood problem. Of particular interest are plasmas produced in air with a peak laser intensity ~10 12 W/cm 2 , the so-called clamping intensity observed in LWIR atmospheric guiding experiments where tunneling and multi-photon ionization operative at near-IR or shorter wavelengths are inoperative. We find that avalanche breakdown on the surface of aerosol (dust) particles can act to seed the breakdown of air observed above the 200 GW/cm 2 threshold when a train of 3 ps 10.6 μm laser pulses separated by 18 ps is used. The breakdown first appears at the best focus but propagates backward toward the focusing optic as the plasma density approaches critical density and makes forward propagation impossible. The velocity of the backward propagating breakdown can be as high as 10 9 cm/s, an order of magnitude greater than measured with ns pulse-produced breakdown, and can be explained rather well by the so-called breakdown wave mechanism. Transverse plasma expansion with a similar velocity is assisted by UV photoionization and is observed as a secondary longitudinal breakdown mechanism in roughly 10% of the shots. When a cm-size, TW power beam is propagated, interception of aerosol particles is guaranteed and several (40 cm -3 ) breakdown sites appear, each initially producing a near-critical density plasma. On a 10 ns–1 μs timescale, shockwaves from each site expand radially and coalesce to produce a large hot gas channel. The radial velocity of the expansion agrees well with the prediction of the blast wave theory developed for ultrafast atmospheric detonations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗
NASA NTRS2003

Pulsed Ejector Wave Propogation Test Program

The development of, and initial test data from, a nondetonating Pulse Detonation Engine (PDE) simulator tested in the NASA Glenn 1 x 1 foot Supersonic Wind Tunnel (SWT) is presented in this paper. The concept is a pulsed ejector driven by the simulated exhaust of a PDE. This pro- gram is applicable to a PDE entombed in a ramjet flowpath, i.e., a PDE combined-cycle propulsion system. The ejector primary flow is a pulsed, uiiderexpanded, supersonic nozzle simulating the supersonic waves ema- nating from a PDE, while the ejector secondary flow is the 1 x 1 foot SWT test section operated at subsonic Mach numbers. The objective is not to study the detonation details, but the wave physics including t,he start- ing vortices, the extent of propagation of the wave front, the reflection of the wave from the secondary flowpath walls, and the timing of these events of a pulsed ejector, and correlate these with Computational Fluid Dynamics (CFD) code predictions. Pulsed ejectors have been shown to result in a 3 to 1 improvement in LID (length-to-diameter) and a near 2 to 1 improvement in thrust augmentation over a steady ejector. This program will also explore the extent of upstream interactions between an inlet and large, periodically applied, backpressures to the inlet as would be present due to combustion tube detonations in a PDE. These interactions could result in inlet unstart or buzz for a supersonic mixed compression inlet. The design of the present experiment entailed the use of an 2-t diagram characteristics code to study the nozzle filling and purging timescales as well as a series of CFD analyses conducted using the WIND code. The WIND code is a general purpose CFD code for solution of the Reynolds averaged Navier-Stokes equations and can be applied to both steady state and time-accurate calculations. The first, proof-of-concept, test entry (spring 2001) pressure distributions shown here indicate the simulation concept was successful and therefore the experimental approach is sound.

Fernandez, Rene↗
NASA NTRS1976

Numerical analysis of a two-dimensional nonsteady detonations

In the present work a system of two-dimensional nonsteady hydrodynamic and chemical kinetic equations was numerically integrated for an exothermic system. Assumed two-step reaction model simulates practically an oxyhydrogen mixture. The calculation starts from a plane Chapman-Jouguet detonation as an initial condition. Two-dimensional disturbances are generated by artificially placing nonuniformities ahead of the detonation front. Regardless of the difference of the given initial disturbances, a fixed number of triple shock waves were produced for a fixed combination of mixture model and geometry when the transition period was over. This shows that for a given detonation tube geometry any exothermic system has its own characteristic multidimensional structure. The obtained number of triple shock waves contained in the detonation front was in agreement with existing experimental observations under the same condition.

Taki, S.↗
DOE OSTI2020

Experimental and modeling analysis of detonation in circular arcs of the conventional high explosive PBX 9501

We examine the diffraction dynamics of a two-dimensional (2D) detonation in a circular arc of the conventional HMX-based, high performance, solid explosive PBX 9501, for which the detonation reaction zone length scale is estimated to be of the order of 100–150 µm. In this configuration, a steady propagating detonation will develop, sweeping around the arc with constant angular speed. We report on results from three PBX 9501 arc experiments, exploring the variation in linear speed on the inner and outer arc surfaces for the steady wave along with the structure of the curved detonation front, as a function of varying inner surface radius and arc thickness. Comparisons of the properties of the motion of the steady wave for each arc configuration are then made with a spatially-distributed PBX 9501 reactive burn model, calibrated to detonation performance properties in a 2D planar slab geometry. We show that geometry-induced curvature of the detonation near the inner arc surface has a significant effect on the detonation motion even for conventional high explosives. We also examine the detonation driving zone structure for each arc case, and thus the subsonic regions of the flow that determine the influence of the arc geometry on the detonation propagation. In addition, streamline paths and reaction progress isolines are calculated. Overall, we conclude that a common approximation for modeling conventional high explosive detonation, wherein the shock-normal detonation speed is assumed equal to the Chapman–Jouguet speed, can lead to significant errors in describing the speed at which the detonation propagates.

42 ENGINEERING↗
DOE OSTI2020

Type Ia supernovae deflagration-to-detonation transition explosions powered by the Zel’dovich reactivity gradient mechanism

ABSTRACT Our aim in this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyse weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for the deflagration-to-detonation transition (DDT) to occur. We observe the formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel’dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated, aided either by reactivity gradients or by collisions of carbon detonations. The observed evolutionary time-scales are found to be sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to the DDT might be of consequence for the chemical composition of the outer ejecta regions of Type Ia supernovae from the single degenerate channel, and offer the potential for validation of the proposed model.

Brooker, E.↗
DOE OSTI2020

Shock-induced kinetics and cellular structures of liquid nitromethane detonation

Using a combination of high-speed diagnostics, optical pyrometry, velocimetry and video photography, we examine the spatiotemporal reaction kinetics of a prototypical high-energy explosive, liquid nitromethane (NM), as we drive this high explosive into detonation. The detonations were initiated by powerful shock waves whose durations (4 ns) were shorter than the characteristic time associated with the reaction (7 ns). Simple optical spectroscopy alone cannot characterize the kinetics with high time resolution because the reactants and products are flowing at a high velocity of ~6 km/s (6 µm/ns). Additionally, every detection volume behind the shock front contains material initiated over a range of times. High spatial resolution can be obtained by probing interfaces where the shock enters or breaks out of the NM. In addition, we obtain cellular patterns imprinted on the luminous shock front by the two-stage explosion in NM. These spatiotemporal patterns arise naturally as a result of the asymmetry produced by the moving shock front. In this way the shock front serves as a thin moving intrinsic optical gauge that reports and characterizes the two-stage NM explosion behind the front.

42 ENGINEERING↗
NASA NTRS1994

Computer program for calculation of complex chemical equilibrium compositions and applications. Part 1: Analysis

This report presents the latest in a number of versions of chemical equilibrium and applications programs developed at the NASA Lewis Research Center over more than 40 years. These programs have changed over the years to include additional features and improved calculation techniques and to take advantage of constantly improving computer capabilities. The minimization-of-free-energy approach to chemical equilibrium calculations has been used in all versions of the program since 1967. The two principal purposes of this report are presented in two parts. The first purpose, which is accomplished here in part 1, is to present in detail a number of topics of general interest in complex equilibrium calculations. These topics include mathematical analyses and techniques for obtaining chemical equilibrium; formulas for obtaining thermodynamic and transport mixture properties and thermodynamic derivatives; criteria for inclusion of condensed phases; calculations at a triple point; inclusion of ionized species; and various applications, such as constant-pressure or constant-volume combustion, rocket performance based on either a finite- or infinite-chamber-area model, shock wave calculations, and Chapman-Jouguet detonations. The second purpose of this report, to facilitate the use of the computer code, is accomplished in part 2, entitled 'Users Manual and Program Description'. Various aspects of the computer code are discussed, and a number of examples are given to illustrate its versatility.

Gordon, Sanford↗
DOE OSTI2024

Calculation of Velocities from Explosive Shot Test Fast-Frame Imagery

A cylinder of the experimental HE is detonated inside of an aquarium and sequential fast-frame images are taken of the resulting shock wave. We wish to use the sequence of images to extract quantitative data regarding the detonation velocity and velocity of the shock wave which can then be used to tune the parameters of the equation of state for the experimental HE material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2022

NASAs Succeeds in Testing of Advanced Rotating Detonation Rocket Engine for New US Space Flight Capability

June-August 2022, NASA engineers at Marshall Space Flight Center in Huntsville Alabama have successfully fired two regeneratively cooled advanced rotating detonation rocket engines (RDRE’s). The engines have accumulated 17 starts at over 600 seconds of total duration. Multiple firings were achieved of greater than 110 seconds each with detonation modes. A single full throttle test produced over 4000 lbf for 15 seconds with detonation modes. The mean pressure at a single point on the injector face was 620 psia. These tests completed the project’s main objective: demonstrating that additive GRCop-alloy hardware could survive long durations while subjected to the detonative events. 4-5 co-rotating detonations were observed during most tests with a single test showing 2-3 waves. Several other milestones were also achieved including successful demonstration of active throttling with detonation modes, successful ignition without a predetonator, and the use of novel additive manufacturing techniques. The primary collaborator was IN Space, LLC (West Lafayette, IN) through an STMD announcement for collaborative opportunity (ACO). Hot fire testing was conducted at Marshall Space Flight Centers heritage east test area at Test Stand 115 in collaboration with Marshall ET10.

Thomas Teasley↗
DOE OSTI2020

Simulations of Sweeping Wave Propagation in a Boron Carbide Plate

In February 2020, experiments were conducted at Los Alamos National Laboratory measuring the propagation of a sweeping wave through a plate of boron carbide which was induced by detonating an adjacent charge of high explosive. Six such shots were red and in each case progress of the wave along the exposed surface of the plate was tracked using photon doppler velocimetry (PDV). In this report we discuss calculations performed using the Lagrangian hydrocode FLAG to model data taken from these experiments. We discuss the models chosen for the different experimental components, the choice of mesh (focusing on, both, two- and three-dimensional geometries), and compare the results from these calculations to experimentally-extracted PDV traces. Notably, we find that the PDV traces for this experiment are most faithfully reproduced using the full three-dimensional geometry, even when the mesh resolution for the three-dimensional calculations was significantly lower than the resolution for the corresponding two-dimensional calculations. We attribute this to the geometry of the experimental setup. Moreover, we find only negligible differences for the PDV traces obtained using the sesame table equation of state for boron carbide and the Gruneisen equation of state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scaling Law for Cloud-Rise Velocity vs. Scaled Height of Burst (SHOB)

The shockwave produced by a nuclear detonation leaves a low-density region in the vicinity of air zero (i.e., the point of detonation). As the shockwave degenerates into an acoustic wave, the buoyant force takes over and causes the low-density region to begin to rise at speeds that can reach as high as ~200-300 mph. The objective of this project was to find out if this initial cloud rise followed any scaling law, by stimulating a cloud rise for a 1 kiloton (kt) detonation and determining if a correlation existed between the scaled-height-of-burst (SHOB) and the scaled-cloud-rise velocity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1971

The development of a sonic boom simulator with detonable gases

A sonic boom pressure profile was simulated in the far-field by detonation of a methane-oxygen mixture contained in a slender, shaped Mylar envelope. Ideal N-waves were synthesized with peak overpressures from two to five psf and durations of 30 to 75 milliseconds. The detonation of the gas mixture was initiated by a single Primacord strand running the length of balloon. The N-wave producing balloon was synthesized as a composite structure, utilizing experimental pressure profiles obtained from the detonations of slender, axisymmetric balloons with elementary, non-cylindrical shapes.

Strugielski, R. T.↗
NASA NTRS1976

Blast wave in a nozzle for propulsive applications

The reported investigation has been conducted in connection with studies concerning the development of a propulsion system based on the use of a detonating fluid propellant. Measurements have been made of the pressure and shock wave velocity in a conical nozzle at various ambient pressures and at an ambient temperature of 25 C. In the experiments a small amount of explosive was placed at the end wall of a conical aluminum nozzle and detonated by a microdetonator inside the nozzle. Differences regarding the characteristics of conventional chemical propulsion and detonation propulsion are illustrated with the aid of a graph. One- and two-dimensional numerical flow calculations were performed and compared with the experimental data.

Varsi, G.↗
NASA NTRS1982

Dynamic effects of combustion

The dynamic effects of combustion are due to the evolution of exothermic energy and its deposition in the compressible medium where the process takes place. The paper examines the dynamics of combustion phenomena, including ignition, turbulent flame propagation (inflammation), explosion, and detonation, with emphasis on their exothermic characteristics. Ignition and explosion are treated as problems of nonlinear mechanics, and their dynamic behavior is described in terms of phase space models and cinematographic laser shear interferograms. The results of a numerical random vortex model of turbulent flame propagation are confirmed in a combustion tunnel experiment, where it was observed that a fresh mixture of burnt and unburnt gases can sustain combustion with a relatively small expenditure of overall mass flow, due to the increasing specific volume of burnt gases inside the flame front. An isentropic pressure wave is found to precede the accelerating flame in the process of detonation, and components of this presssure wave are shown to propagate at local sonic velocities.

Oppenheim, A. K.↗
DOE OSTI2020

New Yield Estimates for Nuclear Detonations Over Water

One of the primary methods used by the United States to estimate the yields of atmospheric nuclear tests performed in the 1940’s, 50’s, and 60’s was to film each detonation with a high-speed camera, measure the blast radius as a function of time, and then use Taylor’s spherical blast wave equation to calculate the source energy (i.e., the yield). For surface detonations, the yield in Taylor’s equation was assumed to be twice as large as the actual yield to account for the energy reflected from the surface, assuming the ground was a perfect reflector. Recent simulations of surface detonations over water (referred to as barge shots) have shown that Taylor’s equation must be modified further to account for material entrainment and heat losses. We have computed this lost energy using a high-order radiation-hydrodynamics code to determine the reduction in blast radius vs. time. By analyzing simulation data in the same manner as film data, we have determined that a yield correction factor of 1.27± 0.04 should be applied to the barge shots.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗