A study of heterogeneous detonations.
Heterogeneous detonations, discussing polydisperse and monodisperse spray detonations and liquid fuel film shock-induced combustion
SEARCH · Search NASA
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.
Heterogeneous detonations, discussing polydisperse and monodisperse spray detonations and liquid fuel film shock-induced combustion
A quasi-two-dimensional, computational fluid dynamic (CFD) simulation of a rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction and other simplifications yield rapidly converging, steady solutions. Viscous effects, and heat transfer effects are modeled using source terms. The effects of potential inlet flow reversals are modeled using boundary conditions. Results from the simulation are compared to measured data from an experimental RDE rig with a converging-diverging nozzle added. The comparison is favorable for the two operating points examined. The utility of the code as a performance optimization tool and a diagnostic tool are discussed.
A Rotating Detonation Rocket Engine (RDRE) model is described which characterizes the device as an essentially infinite number of circumferentially arranged, sequentially firing pulse detonation engine (PDE) tubes. The PDE tubes are in turn treated as lumped-parameter chambers, each of which executes a dynamic (i.e. time dependent) Atkinson cycle with an allowance for a finite Mach number during the refill portion of the cycle. The formulation results in several free parameters which are used to essentially tune the model such that it closely matches higher fidelity computational fluid dynamic RDRE simulations in terms of performance, flow rates, etc. The simplicity of the model allows for straightforward combination with component models (e.g. turbopumps, cooling jackets, etc.) as well as implementation within larger vehicle simulations. Such capability allows for realistic mission analyses and benefits studies which will help to determine the best application for the RDRE concept.
A Rotating Detonation Rocket Engine (RDRE) model is described which characterizes the device as an essentially infinite number of circumferentially arranged, sequentially firing pulse detonation engine (PDE) tubes. The PDE tubes are in turn treated as lumped-parameter chambers, each of which executes a dynamic (i.e. time dependent) Atkinson cycle with an allowance for a finite Mach number during the refill portion of the cycle. The formulation results in several free parameters which are used to essentially tune the model such that it closely matches higher fidelity computational fluid dynamic RDRE simulations in terms of performance, flow rates, etc. The simplicity of the model allows for straightforward combination with component models (e.g. turbopumps, cooling jackets, etc.) as well as implementation within larger vehicle simulations. Such capability allows for realistic mission analyses and benefits studies which will help to determine the best application for the RDRE concept.
Rotating Detonation Rocket Engines (RDRE) have been marketed primarily for their higher specific impulse potential over constant pressure (CP) liquid rocket engines. However, several other performance advantages exist with RDREs over CP engines such as heat transfer advantages for gas expander cycle, increased completeness of combustion at low chamber L*, compact engine design, reduced coolant channel pressure drop potential, and improved injector C* performance. NASA has paved the way for liquid engine system performance enhancement since the Apollo program and continues to do so with metal additive manufacturing (AM), super-alloy materials, and advanced propulsion concepts. A team of propulsion development engineers at NASA are in the process of developing high-performance 7K lbf class RDRE hardware for their potential use in lander, upper stage, and even launch vehicle applications. Clear advantages have been demonstrated with AM including program cost and schedule reductions of up to 50%. It is well known that injector performance is integrally linked to the global performance of a combustion device. This is especially the case for RDREs since detonation stability is heavily dependent on the mixedness of propellants. A major program goal is to rapidly produce ultra-high-performance AM injectors. This paper reviews the available literature on liquid rocket injector design optimization as well as the experimental work conducted to date on injectors tested in RDREs. Major lessons learned are document and suggestions given towards the design of high-performance liquid RDRE injectors. In addition, the integration of metal AM into the design of liquid RDRE injector schemes is discussed. Finally, several candidate AM RDRE injector elements were produced to obtain their diodicity and cold flow characteristics.
The impact of the number of azimuthally propagating waves on the performance of a Rotating Detonation Engine (RDE) is investigated using a simplified two-dimensional computational fluid dynamic simulation. The basic RDE configuration examined has no exit throat. The inlet is assumed lossless and does not allow backflow. The adiabatic, inviscid, and premixed simulation utilizes a particularly simple finite rate reaction mechanism that allows user control over the relative amounts of deflagration and detonation that occur, the presence and extent of a reaction delay associated with fuel and oxidizer mixing, and the number of waves present in the domain. The simplifications and idealizations decouple the effect of wave count from other potential loss mechanisms. Performance is measured using exhaust flow total pressure gain relative to the inlet total pressure. One, two, and three wave solutions are computed under identical boundary conditions and grid resolution. For simulations with no mixing delay and minimal deflagration, the number of waves present has negligible impact on performance. With a mixing delay, performance decreases with increasing wave number. With increased deflagration, performance increases with increasing wave number. With both effects simulated it is found that the two wave solution performs better than either the one or three wave solutions. The causes of these trends are explored. They imply that for practical, RDE’s, the number of waves that are present can impact performance.
The impact of the number of azimuthally propagating waves on the performance of a Rotating Detonation Engine (RDE) is investigated using a simplified two-dimensional computational fluid dynamic simulation. The basic RDE configuration examined has no exit throat. The inlet is assumed lossless and does not allow backflow. The adiabatic, inviscid, and premixed simulation utilizes a particularly simple finite rate reaction mechanism that allows user control over the relative amounts of deflagration and detonation that occur, the presence and extent of a reaction delay associated with fuel and oxidizer mixing, and the number of waves present in the domain. The simplifications and idealizations decouple the effect of wave count from other potential loss mechanisms. Performance is measured using exhaust flow total pressure gain relative to the inlet total pressure. One, two, and three wave solutions are computed under identical boundary conditions and grid resolution. For simulations with no mixing delay and minimal deflagration, the number of waves present has negligible impact on performance. With a mixing delay, performance decreases with increasing wave number. With increased deflagration, performance increases with increasing wave number. With both effects simulated it is found that the two wave solution performs better than either the one or three wave solutions. The causes of these trends are explored. They imply that for practical, RDE’s, the number of waves that are present can impact performance.
Rotating detonation engines (RDEs) represent a novel propulsion technology where incoming reactants are consumed by a self-sustaining detonation wave that travels azimuthally around an annular combustion chamber. In theory, the operating mode of these devices approximates a constant-volume combustion Humphry cycle. This represents a 5-15% increase in thermal efficiency over current rockets and air breathing systems operating on constant-pressure combustion Brayton cycles. The use of liquid fuels in RDEs represent unique challenges in terms of performance and operability. However, the fundamental physics involved in spray breakup, evaporation, and combustion are not well understood inside these devices. This research seeks to provide direct insight into these processes by adding a single liquid fuel injector into an otherwise hydrogen-air RDE to evaluate the breakup and refill dynamics. Visualization is performed by 200-1000 kHz laser induced fluorescence of liquid phase diesel (355 nm) and or kerosene (266 nm).
In combustion simulations, it is usually assumed that highly-energized collision complexes (either radical or stable species) formed in exothermic reactions are always under thermal equilibrium as they evolve through competing networks of reactions. However, in practical flames, non-trivial amounts of reactive radicals such as H, O and OH are present apart from O2. As a result, collisions of these reactive species with the energized collision complexes have been shown in recent studies to induce non-thermal reactivity. These studies have also demonstrated that such non-thermal reactions can be suitably represented in macroscopic kinetics models as chemically termolecular reactions. The present work was focused on identifying and quantifying the effects of including such reactions on the evolution of an initial deflagration front to a developing detonation in H2/CH4-air mixtures under boosted internal combustion (IC) engine conditions. Specifically, fully resolved simulations, with and without non-thermal reactivity, were performed for a constant volume reactor containing stoichiometric H2/CH4-air mixture. It was found that inclusion of termolecular chemistry resulted in a delayed onset of end-gas auto-ignition. Concurrently, the developing detonation intensity was observed to be significantly higher. Chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions responsible for the observed variation in the results.
Soot generated from the detonation of nine explosives was investigated with Raman microspectroscopy and high-resolution transmission electron microscopy (HRTEM). The first order bands, located at approximately 1,580 (G), 1,350 (D1), 1,620 cm –1 (D2), 1,500 (D3), and 1,200 (D4), are used to provide information on the carbon allotropes and defects. The intensity ratios and the full-width at half-maximum provide an indication of the disorder in the material. The Raman spectra of soot acquired from the detonation of nine explosives showed differences which can be used to identify the precursor explosive. Overall, most differences in the spectra were correlated with the presence of amorphous carbon in the soot. In additional to identification of distinct bands in the spectra via deconvolution, principal component analysis was also used to differentiate the different types of soot. HRTEM performed on soot samples with different explosive precursors identified several distinct allotropes of carbon such as nanodiamonds, carbon onions, graphite fibers, and amorphous carbon.
Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.
An update to the numerical shock-fitting scheme from [1] is proposed. The updated scheme locally lowers the order of the scheme in the region of secondary non-fitted waves that interact with the fitted front. In addition, it utilizes a non-uniform mapping from the physical space to the computational space along the streamline direction to offer computational savings over the use of an uniform spacing. The utility of the scheme is first demonstrated on planar ideal gas detonating flows. It is then applied to detonation propagation in a two-dimensional axisymmetric geometry while utilizing an Arrhenius Wescott-Stewart-Davis model for PBX 9502.
Two methods are presented to make measurements of the length of a coaxial cable with submillimeter and submicrosecond accuracy. If the cable is destroyed by the passage of a strong shock, for instance, from the detonation of an explosive, a measure of the phase velocity of the detonation (VOD) can be obtained with high accuracy (≈0.25%). The first method introduces a series of short Gaussian electrical pulses (<300 ps) into the cable that reflect off the cable end and are compared to the time of the original outgoing pulse by using a fast digitizer (100 gigasamples per second). By curve fitting a Gaussian form to the outgoing and reflected pulses, a measurement of the cable length can be made with subsample accuracy. So, as long as the pulses do not overlap, a large number may be present in the cable at any time allowing pulse repetition rates of 40–135 MHz depending on the cable length. Averaging allows accurate VOD measurements to be made that are comparable with the best traditional techniques. The second method replaces the pulser with a wide-spectrum (0.1–10 GHz) amplified random noise source. A similar setup and digitizer are used before autocorrelating short (80–160 ns) subsamples of the total waveform to measure the reflected delay time. That is, a smaller amplitude version of the original waveform appears time shifted by twice the transit time in the cable and this time is estimated via autocorrelation. The advantage of this method is that a new subsample needs only to be selected from a time increment in which new spectral information is present for a new measure of the length to be made. In this manner, many more length estimates can be averaged to establish the VOD than with the Gaussian approach. As a result, a detailed description of the implementation of both methods is presented together with an example of the type of data collected and limitations found.
Fully three-dimensional (3D) microstructure-explicit and void-explicit mesoscale simulations of the shock-to-detonation (SDT) process of pressed granular HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) are performed. The overall size scale of the models is up to 3 × 3 × 15 mm3, with ∼30 000 grains and 206 265 voids. The models account for the heterogeneous material microstructure, constituent distribution, constituent morphology, and voids. Loading conditions considered involve piston velocities in the range of 600–1200 m/s or pressures in the range of 4–8 GPa. The focus is on analyzing the SDT process and the effects of microstructure and voids on the run-to-detonation distance (RDD). Companion two-dimensional (2D) simulations are also carried out to assess the differences between 2D and 3D. Statistically equivalent microstructure sample sets (SEMSSs) are generated and used for both 2D and 3D, allowing the prediction of the statistical and probabilistic Pop plots (PPs). The predictions are in general agreement with trends in available experimental data in the literature. It is found that both the microstructure (heterogeneous grain size, morphology, and size distribution) and voids significantly affect the RDD and the PPs. These effects are systematically delineated and quantified via the use of SEMSSs with different combinations of attributes. A recently developed probabilistic formulation for the PPs is used to characterize the results, allowing uncertainties in the relations between the shock pressure and RDD arising from material heterogeneities to be quantified. The probabilistic formulation is further used to quantify the confidence levels in the ranked order of influences of different combinations of microstructure and voids on the PPs.
A swept-wavelength external cavity quantum cascade laser (ECQCL) is used to perform standoff detection of combustion gases in a plume generated from an outdoor high-explosive (HE) open detonation. The swept-ECQCL system was located at a standoff distance of 830 m from a 41kg charge of LX-14(polymer-bonded high explosive) and was used to measure the infrared transmission/absorption through the post-detonation plume as it propagated through the beam path. The swept-ECQCL was operated continuously to record broadband absorption spectra at 200hz rate over a spectral range from 2050 to 2230 cm -1 (4.48 - 4.88 μm). Fitting of measured spectra was used to determine time-resolved column densities of CO, CO 2 , H 2 O, and N 2 O. Analysis of visible video imagery was used to provide timing correlations and to estimate plume dimensions, from which gas mixing ratios were estimated. Measured emission factors and modified combustion efficiency show good agreement with previously reported values.
A novel set of experiments and reactive flow modeling of pentaerythritol tetranitrate (PETN) is presented. In this work, the specific phenomenon of shock to detonation transition is examined, where an initial, relatively weak shock is propagated into pressed PETN powder at 1.65 g/cm 3 and the subsequent buildup to detonation is observed experimentally. These experiments, in conjunction with reactant and products’ equations of state, are utilized for building reactive flow models.
This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.
Isotope ratio measurements of the actinide elements provide essential information for nuclear detonation forensics and proliferation detection. Resonance Ionization Mass Spectrometry (RIMS) is a high-sensitivity, elementally selective, laser-based form of mass spectrometry that offers the potential to determine the isotopic composition of materials without sample preparation. Due to the elementally selective approach of RIMS, basic research questions of atomic spectroscopy and the probability for producing neutral atoms in the gas phase, must be studied element by element. The studies carried out in this work represent basic research into the application and optimization of RIMS to the analysis of post-detonation debris.