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 109 records · Page 6

Advanced Diagnostics For Megahertz Imaging Of Mixing, Fuel Spray, And Combustion Processes For Rotating Detonation Combustors

Recent advancements in megahertz rate, high-power, burst-mode laser technology are leveraged to perform and explore imaging measurements that spatially and temporally resolve the mixing, combustion, and detonation flow field in two laboratory-scale rotating detonation combustors (RDCs). In a non-premixed annular RDC, multiple imaging diagnostics are explored to investigate gaseous and liquid injector behavior, the detonation wave structure, and the propellant refill. In one instance, OH planar laser-induced fluorescence (OH-PLIF) imaging is performed up to a 2 MHz repetition rate to track the combustion products and reaction zone locations. In the same annular RDC, a single liquid fuel jet is injected, and laser-based 355-nm imaging of the fuel spray is performed up to a 1 MHz repetition rate. For this configuration, the annular RDC is used as a detonation driver to impose periodic detonation waves to interact with the fuel spray. Moreover, in this annular RDC, a range of tracer-based laser imaging measurements are explored to time-resolve the unsteady oxidizer air recovery and refill process. In a non-premixed linear RDC, planar imaging measurements of the fuel mixing are performed up to a 200 kHz repetition rate using PLIF of a tracer in the fuel supply. The fuel mixing imaging helps explain the origin of the observed pre and post wave burning, detonation structure, and enables quantifying injector recovery timescales. This paper will provide a high-level broad survey of diagnostics applied in these RDCs, lessons learned, and interesting observations.

Combustion↗

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗

On the theory of the propagation of detonation in gaseous systems

The existing theory of detonation is critically examined. It is shown that the considerations with which the steady value of the velocity of detonation is chosen are not convincing. In connection with the problem of the process of the chemical reaction in a detonation wave, the objections raised against the conceptions of Le Chatelier and Vieille of the 19th century with regard to the ignition of the gas by the shock wave are refuted. On the basis of this concept, it is possible to give a rigorous foundation for the existing method of computing the detonation velocity. The distributions of the temperature, the pressure, and the velocity in the detonation wave front as the chemical reaction proceeds, are considered. On the assumption of the absence of losses, the pure compression of the gas in the shock wave at the start of the chemical reaction develops a temperature that is near the temperature of combustion of the given mixture at constant pressure.

THEORIES - DETONATION WAVES↗

Mass flow rate effect on a rotating detonation combustor with an axial air injection

Changing the flow rate of reactants being injected into a rotating detonation combustor (RDC) results in interesting behavior of the system. Prior studies have found that an increase in mass flow rate gradually increases the detonation wave speed before splitting the wave into multiple fronts. Here, the focus of this study is in understanding the physics of such behavior through a combination of experiments and numerical simulations. For this purpose, the axial air inlet-based RDC geometry is used. In the experiments, the wave velocity increased similar to prior studies. Corresponding full-scale simulations show that increase in mass flow rate by increasing pressure of the feed plenums shortens the recovery time of the injectors. This causes a more uniform fuel–air mixture to form prior to the arrival of the detonation wave. As a result, a more ideal detonation is observed, which leads to an increased wave velocity. Details of the detonation structure and variation with mass flow rate are analyzed. The presence and variation of the different deflagration zones in such non-premixed RDCs are also discussed.

42 ENGINEERING↗

Design and Operability of a Variably Premixed Rotating Detonation Engine for the Evaluation of Mixing Effects

The design of a variably premixed rotating detonation engine (RDE) is presented with preliminary experimental results. Premixed operation decouples the mixing process from the detonation cycle, enabling direct comparison with premixed simulations. A porous medium, formed into a ring of 6.45 mm radial thickness, serves as the premixture injector, arresting potential flash-back events. Two chamber geometries are tested, the first aimed at eliminating recirculation zones and comprised of a 7.62 mm wide channel, the width of which the porous injector occupies 85%. The second chamber geometry features a 10.7 mm channel width and a 3.08 mm wide backward-facing step on the internal diameter of the premixture injector. Transverse optical access is incorporated into the design using a transparent outer body. Tests were conducted at mixing conditions ranging from non-premixed to fully premixed with gaseous hydrogen and air reactants in the narrow-channel configuration. Engine operating modes and detonation wave speeds were characterized using aft-end high-speed chemiluminescence imaging, while detonation wave topology and heat release in the refill zone are captured by transverse imaging. A two-wave counterpropagating mode was most commonly observed in the narrow-channel configuration with wave speeds ranging from 49-67% of the theoretical Chapman-Jouguet (CJ) velocity. Fully premixed operation in the backward-facing step configuration rendered wave speeds up to 86% of CJ velocity, consistent with previous non-premixed results using the same channel width and a backward-facing step. These results demonstrate successful premixed RDE operation with high wave speeds and highlight the role of product recirculation zones in stabilizing the detonation cycle.

detonation↗

Canonical solutions for unsteady flow fields

The initial value problem of one-dimensional gas-dynamics involving discontinuous, nonuniform initial data is discussed. Canonical solutions which are valid in a small x, t region aroung a discontinuity, and which include the first order effects of nonuniformities in the data, are derived explicitly. The theory is derived by considering a group of elementary piston problems. Solutions with a shock or with a centered expansion wave are worked out individually in order to relate initial flow properties and their gradients to the speed and acceleration of the discontinuity waves. They are then combined to represent the solution of a general initial value problem by regarding the piston path as a contact line. In addition, problems with chemical reaction are discussed in terms of elementary piston problems which involve strong detonation waves, Chapman-Jouguet detonation waves, and deflagration waves.

Liu, G. C.↗

High-Resolution Genuinely Multidimensional Solution of Conservation Laws by the Space-Time Conservation Element and Solution Element Method

In this overview paper, we review the basic principles of the method of space-time conservation element and solution element for solving the conservation laws in one and two spatial dimensions. The present method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. In contrast to the modern upwind schemes, the approach here does not use the Riemann solver and the reconstruction procedure as the building blocks. The drawbacks of the upwind approach, such as the difficulty of rationally extending the 1D scalar approach to systems of equations and particularly to multiple dimensions is here contrasted with the uniformity and ease of generalization of the Conservation Element and Solution Element (CE/SE) 1D scalar schemes to systems of equations and to multiple spatial dimensions. The assured compatibility with the simplest type of unstructured meshes, and the uniquely simple nonreflecting boundary conditions of the present method are also discussed. The present approach has yielded high-resolution shocks, rarefaction waves, acoustic waves, vortices, ZND detonation waves, and shock/acoustic waves/vortices interactions. Moreover, since no directional splitting is employed, numerical resolution of two-dimensional calculations is comparable to that of the one-dimensional calculations. Some sample applications displaying the strengths and broad applicability of the CE/SE method are reviewed.

Himansu, Ananda↗

The Use of Steady and Pulsed Detonations for Propulsion Systems

Objectives of the ODWE concept studies are: demonstrate the feasibility of the oblique detonation wave engine (ODWE) for hypersonic propulsion; demonstrate the existance and stability of an oblique detonation wave in hypersonic wind tunnels; develop engineering codes which predict the performance characteristics of the ODWE including specific impulse and thrust coefficients for various operating conditions; develop multi-dimensional computer codes which can model all aspects of the ODWE including fuel injection, mixing, ignition, combustion and expansion with fully detailed chemical kinetics and turbulence models; and validate the codes with experimental data use the simulations to predict the ODWE performance for conditions not easily obtained in wind tunnels.

Adelman, Henry G.↗

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. Previous analyses of this concept were based on inviscid assumptions. The present results indicate that viscous effects are of primary importance; in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure, and turbulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. In the unstable case, a detonation wave is formed, which propagates upstream and unstarts the ram accelerator. In the stable case, a solution that converges to steady-state is obtained, in which the combustion wave remains stationary with respect to the ram accelerator projectile. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. In addition to these studies, two numerical techniques were tested. These two techniques are vector extrapolation to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices that arise in chemically reacting flows.

Yungster, Shaye↗

Detonation duct gas generator demonstration program

The feasibility of the generation of detonation waves moving periodically across high speed channel flow is experimentally demonstrated. Such waves are essential to the concept of compressing requirements and increasing the engine pressure compressor with the objective of reducing conventional compressor requirements and increasing the engine thermodynamic efficiency through isochoric energy addition. By generating transient transverse waves, rather than standing waves, shock wave losses are reduced by an order of magnitude. The ultimate objective is to use such detonation ducts downstream of a low pressure gas turbine compressor to produce a high overall pressure ratio thermodynamic cycle. A 4 foot long, 1 inch x 12 inch cross-section, detonation duct was operated in a blow-down mode using compressed air reservoirs. Liquid or vapor propane was injected through injectors or solenoid valves located in the plenum or the duct itself. Detonation waves were generated when the mixture was ignited by a row of spark plugs in the duct wall. Problems with fuel injection and mixing limited the air speeds to about Mach 0.5, frequencies to below 10 Hz, and measured pressure ratios of about 5 to 6. The feasibility of the gas dynamic compression was demonstrated and the critical problem areas were identified.

Wortman, Andrew↗

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↗

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

Computational Study of Near-limit Propagation of Detonation in Hydrogen-air Mixtures

A computational investigation of the near-limit propagation of detonation in lean and rich hydrogen-air mixtures is presented. The calculations were carried out over an equivalence ratio range of 0.4 to 5.0, pressures ranging from 0.2 bar to 1.0 bar and ambient initial temperature. The computations involved solution of the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing (TVD) scheme, and a point implicit, first-order-accurate, time marching algorithm. The hydrogen-air combustion was modeled with a 9-species, 19-step reaction mechanism. A multi-level, dynamically adaptive grid was utilized in order to resolve the structure of the detonation. The results of the computations indicate that when hydrogen concentrations are reduced below certain levels, the detonation wave switches from a high-frequency, low amplitude oscillation mode to a low frequency mode exhibiting large fluctuations in the detonation wave speed; that is, a 'galloping' propagation mode is established.

Yungster, S.↗

Numerical and boundary condition effects on the prediction of detonation engine behavior using detailed numerical simulations

High-fidelity numerical simulations of an experimental rotating detonation engine with discrete fuel/air injection were conducted. A series of configurations with different feed-plenum pressures but with constant equivalence ratio were studied. Detailed chemical kinetics for the hydrogen/air system is used. A resolution study for the full rotating detonation engine (RDE) system simulation is also conducted. Two kinds of boundary conditions, a total pressure boundary and a constant mass flow rate boundary, are used to assess the effects of the inlet boundary. As mass flow rate is increased, the total pressure boundary causes more error in the axial pressure distribution while the constant mass flow rate gives a better solution for all cases ran. The simulations confirm experimental findings, and reproduce qualitative as well as some of the quantitative trends. These results demonstrate that a) fuel-air mixing is highly non-uniform within the detonation chamber, leading to variations in local equivalence ratio, b) the fuel and oxidizer injectors experience significant backflow as the detonation wave passes over, but recover at different rates which further augments the inefficiencies in mixing, and c) parasitic combustion in the mixing region makes the detonation wave weak by extending the reaction zone across the wave.

33 ADVANCED PROPULSION SYSTEMS↗

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗

Heat Transfer Measurements of NASA Liquid Kerosene/Oxygen Rotating Detonation Rocket Engine

The RDRE has been identified as a viable high performance propulsion system with numerous advantages over the state-of-the-art (SOA) liquid rocket engine. NASA has been investigating this combustion device for applications ranging from lander, upper stage, thruster, and hypersonics. All of the activities funded to date have been focused on closing major technology gaps preventing the RDRE from being used more broadly by industry. One of those gaps include the prediction and management of heat transfer to the walls from the extreme combustion environment. This work overviews the heat transfer measurements made using new and existing hardware. A liquid oxygen / liquid RP-1 RDRE was tested utilizing a calorimeter outer body and outer body nozzle extension. An actively cooled inner body with axial running cooling channels was also used. A bimetallic GRCop-42 / Monel K500 injector was developed and demonstrated to be a viable technology for RDRE environments. Trends in bulk heat loads are discussed along with a direct comparison to constant pressure theory predictions of wall heat transfer. The experimental data obtained in this investigation showed similar heat fluxes to constant pressure theory indicating deflagration may have dominated the flow field. Only a single detonation wave was observed in all tests which imparted a significant dynamic load (vibration) on the test article making measurement of combustor performances extremely challenging. Both accelerometer and load cell data corroborate the extreme G-forces measured. The single wave mode yielded relatively low performances compared with other hot fire test data sets available. This leads to the conclusion that a single detonation wave in this geometry is not sufficient for high performance. Wave multiplicity, or rather, a specific number of waves may yield higher performances.

Thomas Teasley↗

Propagation velocities of laser-produced plasmas from copper wire targets and water droplets

Experiments were performed to determine the plasma propagation velocities resulting from KrF laser irradiation of copper wire target (75 microns diameter) and water droplets (75 microns diameter) at irradiance levels ranging from 25 to 150 GW/sq cm. Plasma propagation velocities were measured using a streak camera system oriented orthogonally to the high-energy laser propagation axis. Plasma velocities were studied as a function of position in the focused beam. Results show that both the shape of the plasma formation and material removal from the copper wire are different and depend on whether the targets are focused or slightly defocused (approximately = 0.5 mm movement in the beam axis). Plasma formation and its position relative to the target is an important factor in determining the practical focal point during high-energy laser interaction with materials. At irradiance of 100 GW/sq cm, the air plasma has two weak-velocity components which propagate toward and away from the incident laser while a strong-velocity component propagates away from the laser beam as a detonation wave. Comparison of the measured breakdown velocities (in the range of 2.22-2.27 x 10(exp 5) m/s) for air and the value calculated by the nonlinear breakdown wave theory at irradiance of 100 GW/sq cm showed a quantitative agreement within approximately 50% while the linear theory and Gaussian pulse theory failed. The detonation wave velocities of plasma generated from water droplets and copper wire targets for different focused cases were measured and analyzed theoretically. The propagation velocities of laser-induced plasma liquid droplets obtained by previous research are compared with current work.

Song, Kyo-Dong↗