Novel insight into the structure and development of detonation.
Structure and development of detonation waves in gaseous media
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.
Structure and development of detonation waves in gaseous media
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.
The interaction of an oblique, overdriven detonation wave with a vorticity disturbance is investigated by a direct two-dimensional numerical simulation using a multi-domain, finite-difference solution of the compressible Euler equations. The results are compared to those of linear theory, which predict that the effect of exothermicity on the interaction is relatively small except possibly near a critical angle where linear theory no longer holds. It is found that the steady-state computational results agree with the results of linear theory. However, for cases with incident angle near the critical angle, moderate disturbance amplitudes, and/or sudden transient encounter with a disturbance, the effects of exothermicity are more pronounced than predicted by linear theory. Finally, it is found that linear theory correctly determines the critical angle.
Paper submitted to 2023 AIAA SciTech Forum. This paper discusses the occurrence of stable, but unsteady (periodic) variations in detonation wave spacing in a Rotating Detonation Combustor.
Utilization of detonation wave to investigate transient processes in gas wave dynamics of rocket propulsion systems
Thermodynamics and propagation of carbon dioxide - air detonation waves at various initial pressures
Oscillation frequency, stability, and structure of detonation waves using von neumann theorem
A rotating detonation engine (RDE) is a realization of pressure-gain combustion, wherein a traveling detonation wave confined in a chamber provides shock-based compression along with chemical heat release. Due to the high wave speeds, such devices can process high mass flow rates in small volumes, leading to compact and unconventional designs. RDEs involve unsteady and multiscale physics, and their operational characteristics are determined by an equilibrium between large- and small-scale processes. While RDEs can provide a significant theoretical gain in efficiency, achieving this improvement requires an understanding of the multiscale coupling. Specifically, unavoidable nonidealities, such as unsteady mixing, secondary combustion, and multiple competing waves associated with practical designs, need to be understood and managed. The secondary combustion processes arise from fuel/air injection and unsteady and incomplete mixing, and can create spurious losses. In addition, a combination of multiple detonation and secondary waves compete and define the dynamical behavior of mixing, heat release distribution, and the overall mode of operation of the device. This review discusses the current understanding of such nonidealities and describes the tools and techniques used to gain insight into the extreme unsteady environment in such combustors.
The interaction of an oblique, over driven detonation wave with a vorticity disturbance is investigated by a direct two-dimensional numerical simulation using a multi-domain, finite-difference solution of the compressible Euler equations. The results are compared to those of linear theory, which predict that the effect of exothermicity on the interaction is relatively small except possibly near a critical angle where linear theory no longer holds. It is found that the steady-state computational results agree with the results of linear theory. However, for cases with incident angle near the critical angle, moderate disturbance amplitudes, and/or sudden transient encounter with a disturbance, the effects of exothermicity are more pronounced than predicted by linear theory. Finally, it is found that linear theory correctly determines the critical angle.
Steady detonation wave velocity in gaseous ozone in tubes of various diameters at initial pressure range and in mixtures with oxygen or argon
Reaction engine operates on the principle of a controlled condensed detonation. In this engine the gas products that are expelled from the engine to produce thrust are generated by the condensed detonation reaction. The engine is constructed of two basic sections consisting of a detonation wave generator section and a condensed detonation reaction section.
Detonation wave study for application to chemical propulsion systems, nuclear propulsion, and pulsed plasma accelerators
A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.
Detonation wave initiation in gaseous media, noting transition process for flames accelerating to detonation and explosions behind reflected shock wave
All liquid propellant rocket instability calculations in current use have limited value in the predictive sense and serve mainly as a correlating framework for the available data sets. The well-known n-t model first introduced by Crocco and Cheng in 1956 is still used as the primary analytical tool of this type. A multitude of attempts to establish practical analytical methods have achieved only limited success. These methods usually produce only stability boundary maps that are of little use in making critical design decisions in new motor development programs. Recent progress in understanding the mechanisms of combustion instability in solid propellant rockets"' provides a firm foundation for a new approach to prediction, diagnosis, and correction of the closely related problems in liquid motor instability. For predictive tools to be useful in the motor design process, they must have the capability to accurately determine: 1) time evolution of the pressure oscillations and limit amplitude, 2) critical triggering pulse amplitude, and 3) unsteady heat transfer rates at injector surfaces and chamber walls. The method described in this paper relates these critical motor characteristics directly to system design parameters. Inclusion of mechanisms such as wave steepening, vorticity production and transport, and unsteady detonation wave phenomena greatly enhance the representation of key features of motor chamber oscillatory behavior. The basic theoretical model is described and preliminary computations are compared to experimental data. A plan to develop the new predictive method into a comprehensive analysis tool is also described.
The objective of the presented work was to develop validated computational fluid dynamics (CFD) based methodologies for predicting propellant detonations and their associated blast environments. Applications of interest were scenarios relevant to rocket propulsion test and launch facilities. All model development was conducted within the framework of the Loci/CHEM CFD tool due to its reliability and robustness in predicting high-speed combusting flow-fields associated with rocket engines and plumes. During the course of the project, verification and validation studies were completed for hydrogen-fueled detonation phenomena such as shock-induced combustion, confined detonation waves, vapor cloud explosions, and deflagration-to-detonation transition (DDT) processes. The DDT validation cases included predicting flame acceleration mechanisms associated with turbulent flame-jets and flow-obstacles. Excellent comparison between test data and model predictions were observed. The proposed CFD methodology was then successfully applied to model a detonation event that occurred during liquid oxygen/gaseous hydrogen rocket diffuser testing at NASA Stennis Space Center.
PBX 9012 is a plastic bonded explosive composed of 90.2 wt % HMX (cyclo-tetramethylene-tetranitramine) and 9.8 % Viton A (vinylidene-hexafluoropropylene copolymer). It has nearly the same composition as LX-07 (90% HMX and 10% Viton A) but with a slightly lower density. Hence it has a higher porosity and is slightly more shock sensitive. Here we describe the SURFplus model calibration for PBX 9012 using data from five shock-to-detonation transition (SDT) gas gun experiments by Burns and Chiquete, detonation front shape data from rate stick experiments and cylinder test experiment reported in Anderson et al., Jackson et al. The SDT experiments provide Pop plot data points (run distance-to-detonation) which characterize shock initiation, Lagrangian velocity time histories from embedded magnetic velocity gauges which give information on the reaction behind the lead shock, and Hugoniot data points for the reactants EOS. The rate stick experiments provide data on the curvature effect which characterize propagating detonation waves. The cylinder test experiments provide data on the detonation release isentrope used to calibrate the products EOS. The HE model uses a Davis reactants EOS with parameters that Burns and Chiquete [2020] fit to Hugoniot data, and Davis products EOS with parameters that Anderson et al. [2020] fit to cylinder test data.
Deflagration to Detonation (DDT) was predicted to occur in porous beds of high-energy solid propellants by solving the unsteady fluid mechanical convective heat transfer from hot gas products, obtained from the rapid burning at high pressures, provides the impetus to develop a narrow combustion zone and a resulting strong shock. A parametric study clearly indicates that DDT occurs only when a combination of the solids loading fraction, the burning rate constants, the propellant chemical energy, and the particle size provide for critical energy and gas release to support a detonation wave. Predictions for the run-up length to detonation as a function of these parameters are presented.