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.
Conference paper for the ASME 2022 Turbomachinery Technical Conference & ExpositionJune 13-17, 2022, Rotterdam, The Netherlands<br>
Coal dust explosions can be hazardous; however, they can also generate a significant rise in stagnation pressure if adequately harnessed. Rotating detonation combustors seek to take advantage of the stagnation pressure rise phenomenon in a more sustained and controlled manner via confinement to a physical annulus, leading to increased overall thermodynamic efficiency. Here this investigation presents an analysis of detonations fueled by Carbon Black, a solid particulate consisting of virtually pure carbon molecules and lean Hydrogen-Air mixtures. It is realized that with the addition of Carbon Black, an increase of lean mixture detonability and detonation velocities extending the operating limit over that of a pure hydrogen-air mixture is experienced. For all testing conditions, the total equivalence ratio is held at φ = 1, while the fuel mixture's carbon mass fraction is increased from 0 to 0.7 while the hydrogen is decreased. Detonation wave velocities are extracted from high-speed imaging through applying a Discrete Fourier Transform algorithm to determine changes to the wave speed as Carbon Black particles are introduced. As a result, due to the addition of Carbon Black as an auxiliary fuel source, detonations were formed instead of deflagrations in operating conditions where one would expect deflagrations at the same hydrogen-air equivalence ratios without Carbon Black addition. The detonation formation provides evidence that the coal particles are reacting within the detonation wave in a large enough capacity to support a detonation wave within the annulus. Furthermore, the wave speed is shown to increase with the additional of carbon particles. At a constant global equivalence ratio, the detonation wave velocities were found to decrease with hydrogen's incremental replacement with coal particles. Whereby, through a theoretical comparison of the heat of combustion as computed from the experimentally derived detonation wave velocities, a linear relationship of the two was shown to exist. Therefore, the heat of combustion has the potential to describe an operational limit to sustaining a detonation wave.
Explosives are often used in industry, geology, mining, and other applications, but it is not always clear what remains after a detonation or the fate and transport of any residual material. The goal of this study was to determine to what extent intact molecules of high explosive (HE) compounds are detectable and quantifiable from post-detonation dust and particulates in a field experiment with varied topography. We focused on HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane), which is less studied in field detonation literature, as the primary explosive material and RDX (1,3,5-Trinitroperhydro-1,3,5-triazine) as the secondary material. The experiment was conducted at Site 300, Lawrence Livermore National Laboratory’s Experimental Test Site, in California, USA. Two 20.4 kg and one 40.8 kg above ground explosions (primarily comprised of LX-14, an HMX-based polymer-bonded high explosive) were detonated on an open-air firing area on separate days. The complex terrain of the firing area (e.g., buildings, berm, low-height obstacles) was advantageous to study HE deposition in relation to plume dynamics. Three types of samples were collected up to 100 m away from each shot: surface swipes of aluminum plates, surface swipes of fixed objects, and filters from air samples. We used atmospheric flow tube-mass spectrometry (AFT-MS) to quantify picogram levels of molecular residue of HE material in the post-detonation dust. An aliquot of sample extract in methanol (e.g., 1 µL of 0.5 mL) was placed onto a resistive material and then thermally desorbed into the AFT-MS. We successfully detected and quantified both HMX and RDX in many of the samples. Based on mass (pg) detected and solution dilution, we back-calculated the mass collected on the swipe or filter (ng per sample). The aerial distribution of molecular residue was consistent with the path of the plume, which was strongly determined by wind speed and direction at the time of each shot. The quantity of material detected appeared to correlate more with distance from the shot and the wind conditions than with shot size. This study demonstrates that the picogram detection levels of AFT-MS are well-suited for quantification of analytes (e.g., HMX and RDX) in environmental samples.
The impact of periodic detonation wave impact on a liquid fuel jet is investigated in a linear detonation combustor. The linear detonation combustor operated with gaseous natural gas and oxygen generates sustained, self-excited detonation waves that propagate along its length at approximately 8 kHz, representing a wave propagation frequency in typical rocket rotating detonation engines. The effect of the detonation wave on the dynamic injection and break-up of a single diesel jet injected into the combustion chamber at varying injection pressures is evaluated with chemiluminescence, fuel planar laser induced fluorescence and Mie scattering measurements at 100 kHz. The detonation wave significantly impacts the liquid jet trajectory with its deflection in both windward and leeward direction as the adverse pressure gradient across it changes between wave passages. The maximum recovery height of the liquid jet is observed to be consistent across all operating conditions, and dependent on the detonation wave strength in the chamber.
Thermochemical models of detonation are widely used to estimate energy delivery, but they are based on the assumption that the carbon-rich condensates (soot) formed during detonation are very similar to bulk carbon. Here, we present an analytic equation of state (EOS) based on experimental detonation data for soot formed during the detonation of triaminotrinitrobenzene (TATB)-based high explosives. X-ray photoelectron spectra of several detonation soots are used to determine the elemental nitrogen abundance, with surprisingly high values for TATB. The proposed TATB soot EOS is highly compressible at low pressures and shares some features of glassy carbon, exhibiting graphite- and diamond-like behavior as a function of pressure. We demonstrate the influence of formed soot on detonation performance, including a lowering of the detonation velocity at typical charge densities, and a more compressive product Hugoniot at overdriven conditions. The soot model improves the accuracy of thermochemical calculations for TATB-based explosives across a wide range of states. Detonation velocity predictions for HMX (cyclotetramethylene-tetranitramine)-TATB blends with 80% or more TATB content, as well predictions for 1,3-diamino-2,4,6-trinitrobenzene (DATB) and 3-nitro-1,2,4-triazol-5-one (NTO), which share some features with TATB, are also improved.
Although inherently unstable, existing research in rotating detonation combustion supports its application in notionally steady processes resulting in greater availability compared to conventional, constant pressure combustion. Further improvements rely on a more in-depth understanding of system losses and identifying conditions which optimize device performance. Within this study, the presence and proportion of ideal and nonideal combustion regimes are compared across a variety of process conditions and wave modes. Here, large-scale data analysis seeks to summarize proportional heat release associated with commensal, parasitic, and detonative combustion averaged across individual traces of OH* chemiluminescent data acquired at the detonation plane. Means of regime partitioning based on the anatomy of the time-resolved OH* signal are proposed to ensure consistent analysis throughout the current and future studies concerning combustion regimes. Of particular interest is the possible influence of wave on the nonideal combustion in relative proportion to the desired detonation. Results showed improved percent detonation with increasing significance for the following trends: decreasing equivalence ratio, increasing wave count, decreasing wave velocity, and increasing detonation time. Increased wave number, brought on by decreased equivalence ratios and wave velocities, is thought to decrease fill region surface area, and therefore, decrease nonideal contact burning. Additional analysis is performed to consider possible trend variation due to the presence of stable galloping waves, which were found to have minimal influence on relative percent detonation values. The outcome of this study suggests operational states, which correspond to increased wave quantities for increased proportions of reactants consumed by the targeted detonative combustion regime.
Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.
A variably premixed rotating detonation engine using gaseous hydrogen and air reactants is introduced to enable investigation of key cycle processes while varying the homogeneity of the reactant inlet conditions. Two chamber configurations are investigated, the first with reactants filling the entire span of a straight annular channel and the second with a slightly larger channel and a backward-facing step. The first configuration permits both premixed and non-premixed fuel injection, enabling mixing quality modulation. The second configuration is operated only at fully premixed conditions. Operating modes and detonation wave speeds are characterized using exhaust-plume imaging, while the chamber heat release field is captured by transverse imaging through a transparent outer body. Tests using the first configuration were characterized by unstable counterpropagating modes with low detonation wave speeds regardless of the state of premixing, while the second configuration rendered single-wave behavior with wave speeds up to 86% of the Chapman–Jouguet velocity. Comparisons with a simple computational fluid dynamics model of the second configuration indicate that reactant preheating significantly influences the detonation wave topology, highlighting the potential utility of the test platform for isolating key physics associated with the effects of reactant premixing, preheating, and chamber geometry on rotating detonation engine operation.
Disturbances caused by the detonation wave in a rotating detonation combustor (RDC) propagate upstream through the inlet, and can potentially affect and couple to upstream components, such as turbomachinery or isolators. These disturbances can potentially also affect the operation of the RDC itself. By drawing from the analogy of a detonation wave bounded by an inert gas, the pressure disturbances observed upstream of the inlet are explained as the consequence of the passage of an upstream propagating oblique shock. In this study, the pressure rise in the plenum from the oblique shock is measured in an axial air inlet RDC. The speed of the upstream propagating wave is estimated to be moderately above the acoustic speed of the oxidizer in the plenum. The wave propagates into the plenum despite local regions of choking in the inlet. It is estimated that the time it takes a fluid particle to transit from the plenum to the detonation channel through the inlet is much larger than the rotational time of the detonation wave. This implies that a fluid particle experiences multiple shocks prior to entering the detonation channel. The oblique shock propagating upstream through the inlet area change is modeled by leveraging an analogy with a quasi-1D shock wave moving in a variable area duct with mean (incoming) flow. Due to flow expansion along the area change, fluid particles are found to experience stronger shocks in the inlet than in the plenum, thereby creating different thermodynamic states within the fill region as the oxidizer emerges from the inlet.
This paper presents results from numerical simulations of a non-premixed hydrogen-air rotating detonation combustor with radial injection. The fuel and air mass flow rate are varied in order to hold a unity global equivalence ratio. The calculations show that multiple detonation waves co-exist when the mass flow rate is increased. Conditional statistics of the detonation structure and combustion processes suggest similarities across co-existing waves. Quantification of the injection response to the rotation of a detonation indicates that at higher flow rate the refill time is short enough to allow for a quick and well mixed composition prior to the new front passage. Details of the combustion characteristics are analyzed. Finally, the results elucidate the correlation between initial injection conditions and detonation multiplicity on the overall physics within the combustor.
Injector response to self-excited and sustained detonation wave propagation in an optically-accessible nonpremixed natural gas-oxygen linear detonation combustor is investigated. 100 kHz measurements of planar laser induced fluorescence (PLIF) of trace quantities of acetone injected into the fuel supply and simultaneous OH∗ chemiluminescence measurements were used to correlate the effect on injection dynamics on detonation wave structure and propagation. Measurements over a range of equivalence ratios (0.61-1.48) revealed wave propagation through regions of high fuel stratification in the axial and transverse direction of reactant injection. Consequently, the propagation of the detonation wave along the channel is through a flow-field with significant variation in fuel mass fraction, reactivity, and sound speeds. At lean conditions, wave propagation frequency is generally lower (∼6 kHz) resulting in longer time between wave passages for mixing and pre-heat resulting in steeper detonation wave-fronts anchored near the injector exit plane. Conversely, at rich equivalence ratios the wave frequency is higher (∼10 kHz) resulting in shorter intra-cycle period for reactant fill and mixing resulting waves lifted from the injection plane that are more compact. Higher chemiluminescence intensity in the post-wave region with deflagrative combustion is observed in these cases. Phase averaged representation of the acetone-PLIF measurements indicate that the fraction of time between wave passages available for reactant fill, mixing and pre-heat is consistent across all equivalence ratios.
The investigation, by use of two-dimensional numerical hydrodynamics simulations, of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is continued. Previously we found that the deflagration is insufficient to unbind the star. Expansion shuts off the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. The burning was allowed to develop into a detonation in these nonspherical models. The detonation grows toward spherical symmetry at late times. At these densities (rho approx. 10(exp 7) to 10(exp 8) g cm(exp -3)), either Ni-56 or nuclei of the Si-Ca group are the dominant products of the burning. The bulk yields are sensitive to the density of the star when the transition to detonation occurs. The relevance of the abundances, velocities, mixing, and total energy release to the theory and interpretation of Type Ia supernovae is discussed.
As rotating detonation engines (RDEs) progress in maturity, the importance of monitoring advancements toward development of active control becomes more critical. Experimental RDE data processing at time scales which satisfy real-time diagnostics will likely require the use of machine learning. This study aims to develop and deploy a novel real-time monitoring technique capable of determining detonation wave number, direction, frequency, and individual wave speeds throughout experimental RDE operational windows. To do so, the diagnostic integrates image classification by a convolutional neural network (CNN) and ionization current signal analysis. Wave mode identification through single-image CNN classification bypasses the need to evaluate sequential images and offers instantaneous identification of the wave mode present in the RDE annulus. Here, real-time processing speeds are achieved due to low data volumes required by the methodology, namely one short-exposure image and a short window of sensor data to generate each diagnostic output. The diagnostic acquires live data using a modified experimental setup alongside Pylon and PyDAQmx libraries within a python data acquisition environment. Lab-deployed diagnostic results are presented across varying wave modes, operating conditions, and data quality, currently executed at 3–4 Hz with a variety of iteration speed optimization options to be considered as future work. These speeds exceed that of conventional techniques and offer a proven structure for real-time RDE monitoring. The demonstrated ability to analyze detonation wave presence and behavior during RDE operation will certainly play a vital role in the development of RDE active control, necessary for RDE technology maturation toward industrial integration.
Presentation from ASME Turbo Expo 2021
Limit cycle oscillation (LCO) detonation wave behaviors are presented and analyzed for test times exceeding 20 s in a water-cooled rotating detonation engine (RDE). LCO detonation waves exhibit cyclic acceleration and deceleration, resulting in oscillating wave spacing at unique process conditions. In previous RDE studies, similar behaviors have been studied as microsecond-scale instabilities leading to ascending or descending modal transitions. In the current work, however, LCO waves are considered a persistent wave mode, occupying unique portions of the operational envelope adjacent to those of their equally spaced counterparts. These occurrences of LCO waves are repeatable, enduring behaviors. A method to generate shifted contour surfaces specifically intended to extract and analyze wave spacing variation through time, termed limit cycle oscillation visualization (LCOV) surfaces, is presented. LCOV surfaces transform data into the reference frame of a primary traveling wave and are used to analyze quasi-steady, short-timescale, and transitional LCO modes. Results are leveraged to understand the relationship between fill height, wave strength, local wave acceleration, and subsequent LCO wave spacing for individual wave sets. In conclusion, quasi-steady LCO waves display wave spacing oscillations between equal spacing values associated with ±1 wave across runs exceeding 18 s.
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.
Within this study, the presence and proportion of ideal and non-ideal combustion regimes is compared across a variety of process conditions and wave modes. Large-scale data analysis seeks to summarize proportional heat release associated with commensal, parasitic, and detonative combustion averaged across individual traces of OH* chemiluminescent data acquired at the detonation plane. Means of regime partitioning based on the anatomy of the time-resolved OH* signal are proposed to ensure consistent analysis throughout the current and future studies concerning combustion regimes.
Presentation presenting results of the NETL water-cooled RDE.