Ignition and chemical kinetics in hypersonic ramjets
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
High throughput ignition screening would be useful for evaluating fuels before scaled production. One method for doing this involves generating a map of ignition delay times over different conditions for each fuel blend. Constant condition ignition delay time maps can be used to estimate ignition for engines in prescribed thermodynamic trajectories. A Miniaturized Ignition Screening Rapid Compression Machine (MISR) has been developed and built at the University of Illinois Chicago (UIC) intending to generate those constant condition ignition delay time maps for low temperature chemistry relevant conditions. Some design features include a bellows-actuated driver section, which is lighter compared to other RCM designs, an adjustable brake section, and full optical access. The setup was designed to achieve reaction conditions (T > 600K and P < 100 bar with repetition rates of about 8 seconds per experiment). These conditions were meant to be relevant to low-to-intermediate temperature auto-ignition and for generation of data for developing relevant chemical kinetics models or ignition maps. Because of the compression characteristics and the heat loss characteristics of the MISR, the apparatus currently tests higher temperatures over shorter test times than what would be desirable, however there are some pathways to improved compression characteristics as detailed in this report. A detailed description of the MISR design and its operating mechanism is discussed. The team faced several challenges and setbacks in the construction and commissioning of the MISR, and those are highlighted, along with the modifications that were performed to get the MISR operational. To demonstrate the performance characteristics of the MISR, limited ignition delay experiments were conducted for ethanol fuel at different equivalence ratios (Φ), and those results were usable, but again test a higher temperature and shorter ignition delay time than initially intended.
Abstract Ignition delay times from undiluted mixtures of natural gas (NG)/H2/Air and NG/NH3/Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000–1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data is compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH3 fuel addition made no noticeable difference in ignition time. NG/NH3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. The reaction CH3 + O2 = CH3O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.
Kinetics in hydrogen air flow systems - calculation of autoignition delays for hypersonic ramjets
Methanol is a potentially attractive fuel for marine and off−road engines owing to its availability at bunkering and global distribution locations. Although methanol is well−distributed worldwide, its fuel chemistry and ignition properties make it poorly suited as a direct drop−in replacement for diesel fuel in compression−ignition engines. However, industrial processes are regularly used to convert methanol, via catalytic dehydration, to dimethyl ether (DME) over nonprecious metal catalysts. This chemical conversion can occur at relatively low pressures, temperatures, and catalyst space velocities, highlighting a potential opportunity to generate DME via onboard catalytic dehydration of methanol. DME’s fuel kinetic and ignition properties for compression ignition are much more favorable than those of methanol or even diesel fuel, but DME is more challenging than diesel fuel or methanol to pump, store, and deliver through conventional diesel fueling injection hardware. Thus, a potential opportunity exists to use the ignition and kinetic properties of DME, with the transportation and delivery advantages of methanol, in a methanol−fueled mixing−controlled compression−ignition engine. The present work explores performance, combustion behavior, and emissions reduction opportunities for methanol mixing−controlled combustion, enabled by a HCCI of DME that represents a small fraction of the total fuel energy that can be generated onboard via catalytic dehydration of methanol.
Flame structure and ignition and extinction processes associated with a strained fuel strip are explored numerically using detailed transport and complex kinetics for a propane-air reaction. Ignition modes are identified that are similar to those predicted by one-step activation energy asymptotics, i.e., modes in which diffusion flames can ignite as independent or dependent interfaces and modes in which single premixed or partially premixed flames ignite and burn. These ignition modes have been found to be dependent on critical combinations of strain rate, fuel strip thickness and initial reactant temperatures. The formation of NO/NO2 is found to be strongly dependent on strain rate and the local molecular mixing of reactants which occurs as a consequence of strain. Extinction in this configuration is seen to occur due to fuel consumption by adjacent flames, although viscosity is seen to have the effect of delaying extinction by reducing the effective strain rate experienced by the flames. Response of the flames to oscillatory strain rates is seen to be strongly dependent on the amplitude and frequency of the oscillation.
Ignition delay times from undiluted mixtures of natural gas (NG)/H 2 /Air and NG/NH 3 /Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000 and 1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data are compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H 2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH 3 fuel addition made no noticeable difference in ignition time. NG/NH 3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. Here, the reaction CH 3 + O 2 = CH 3 O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.
A pressure dependent cookoff model for PBX 9501 was developed. This cookoff model was implemented in the finite element (FE) heat transfer code Aria and was used to simulate a set of cookoff experiments. This set is called the Large Scale Annular Cookoff (LSAC) experiments. Three dimensional and axisymmetric heat conduction FE models were implemented to simulate these experiments. The predictive ability of the PBX 9501 pressure dependent cookoff model was evaluated using Latin Hypercube sampling (LHS). Predictions of the thermal times to ignition and temperatures are compared with the experiments.
An anaytical study of hydrogen air kinetics was performed. Calculations were made over a range of pressure from 0.2 to 4.0 atm, temperatures from 850 to 2000 K, and mixture equivalence ratios from 0.2 to 2.0. The finite rate chemistry model included 60 reactions in 20 species of the H2-O2-N2 system. The calculations also included an assessment of how small amounts of the chemicals H2O, NOx, H2O2, and O3 in the initial mixture affect ignition and reaction times, and how the variation of the third body efficiency of H2O relative of N2 in certain key reactions may affect reaction time. The results indicate that for mixture equivalence ratios between 0.5 and 1.7, ignition times are nearly constant; however, the presence of H2O and NO can have significant effects on ignition times, depending on the mixture temperature. Reaction time is dominantly influenced by pressure but is nearly independent of initial temperature, equivalence ratio, and the addition of chemicals. Effects of kinetics on reaction at supersonic combustor conditions are discussed.
A data-driven chemical kinetic mechanism for the military version of Jet A, F-24, is developed for numerical simulations of the ignition process. The main purpose of this study is to obtain a practical F-24 mechanism across wide temperature and equivalence ratio ranges, with a particular focus on the negative temperature coefficient and low temperature regions. The new mechanism (ARLMech-HC-F24) is based on the HyChem model of a similar fuel and optimized using a micro-genetic algorithm against an experimental ignition delay data set of the target fuel. The development and optimization processes include reaction selection, population creation, shuffled tournament implementation based on a merit function, and child-individual creation for the next generation. Several techniques and parameters are proposed to generate an accurate mechanism through an efficient process. The newly introduced data-driven mechanism based on these techniques shows better merit value convergence and represents the ignition behavior more accurately than that without the techniques. This practical mechanism is suitable for the numerical simulations of the F-24 or Jet A ignition problem, and the suggested strategies can be employed in similar problems of rate coefficients estimation.
Ignition delay times (IDTs) of 2-methyl decane (C 11 H 24 -2) are measured in a high-pressure shock tube and in a rapid compression machine at equivalence ratios in the range 0.5–2.0 at 90 % dilution, at temperatures in the range 600–1430 K and at pressures of 15 and 30 bar. To clarify the effect of the branched methyl group on fuel oxidation, IDTs of n-undecane (nC 11 H 24 ) are also measured at similar conditions to those measured for C 11 H 24 -2. A new chemical kinetic mechanism, using C3MechV4.0.1 as the core chemistry, is developed and validated against the new experimental data. The thermodynamic properties of the fuel (RH), alkyl (Ṙ), alkyl peroxy (RȮ 2 ), hydroperoxy-alkyl (Q̇OOH), and peroxy hydroperoxy alkyl (Ȯ 2 QOOH) radicals are updated in both the C 11 H 24 -2 and nC 11 H24 models using THERM25. A reaction path flux analysis for C 11 H 24 -2 at different temperatures was conducted. Compared to nC 11 H 24 , C 11 H 24 -2 shows slower reactivity. At low and intermediate temperatures, the chain propagation pathway Q̇OOH ↔ C 11 cyclic ether + ȮH is favored for C 11 H 24 -2, while the chain branching pathway Q̇OOH ↔ Ȯ 2 QOOH ↔ C 11 carbonyl hydroperoxide (KHP) + ȮH is suppressed, leading to lower reactivity compared to nC 11 H 24 . At high temperatures, the presence of the branched methyl group inhibits the direct decomposition of the fuel, resulting in reduced C 2 H 4 formation, which in turn suppresses the reactivity of the fuel.
Fuel blending effects on chemically-dominated fuel properties, such as gasoline anti-knock quality, are influenced by fundamental chemical kinetic interactions between the blending agent and the base fuel. Historically, quantification of such interactions has focused on changes to the radical pool, including $\dot{O}$H and HO 2 , while intermolecular interactions pertaining to carbonated, non-fuel-specific intermediates are typically overlooked. In this regard, this work aims to derive new insight into intermolecular chemical kinetic interactions that are intrinsic to fuel blending effects via a case study on blends of 0–30% ethanol (by volume) into FGF-LLNL (a multi-component gasoline surrogate for FACE-F research gasoline) using a rapid compression machine at a diluted/stoichiometric fuel loading, compressed pressure of 40 bar and low- to intermediate-temperature regimes that are representative of boosted SI engine operation. Ethanol blending effects on the intermediate temperature heat release (ITHR) of FGF-LLNL are characterized using experimental measurements, where ethanol is found to promote the extent of ITHR and suppress the transition from ITHR to main ignition. Chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model. Sensitivity analyses on ITHR characteristics further corroborate the ethanol blending effects, and highlight the significant dependence of ITHR on both fuel-specific and non-fuel-specific reactions. An approach allowing comprehensive characterization of the complex intermolecular chemical kinetic interactions between constitutes in a fuel blend is then proposed. Application of the approach to FGF-LLNL/E0–E30 reveals that ethanol perturbs the heat release and autoignition characteristics of FGF-LLNL not only by directly changing the $\dot{O}$H and HO 2 radical pools via fuel-specific reactions, but also through intermolecular interactions where participating intermediates can be produced and consumed by various sub-chemistries. Disabling the intermolecular interactions in carbonated species between ethanol and FGF-LLNL sub-chemistries leads to somewhat slower ITHR evolution and lower ignition reactivity. The role of the individual intermolecular interaction is also characterized using the proposed approach. Lastly, implications of intermolecular interactions for future studies that aim to improve model performance are highlighted, where it is found that further investigations on the core C0-C4 chemistry are needed for developing highly accurate chemistry models for complex fuel blends.
The sputter deposition of alternating layers of Ni(V) and Al forms a reactive multilayer known to undergo self-propagating formation reactions when ignited. The sequential deposition process leads to nanometer-scale premixing of reactants at each included interface, which ultimately affects multilayer exothermicity. This work performs the direct measurement of a disordered face-centered cubic (FCC) solid solution premixed phase at the interfaces of Ni(V)/Al multilayers via scanning transmission electron microscopy. The crystallinity of the observed phase differs from previously reported a priori predictions of an amorphous interlayer. The disordered FCC phase retains its symmetry after annealing for 16 h at 135 ± 5 °C, but the lattice parameter shifts consistently with an Al-rich composition. The existence of a crystalline premix in Ni(V)/Al is attributed to the electronic contribution to the entropy of crystallization. The importance of electronic entropy to the phase formation of energetic materials motivates its inclusion when constructing digital twins for atomistic kinetics and ignition sensitivity.
Recent interest in nonequilibrium plasma discharges as sources of ignition for the automotive industry has not yet been accompanied by the availability of dedicated models to perform this task in computational fluid dynamics (CFD) engine simulations. The need for a low-temperature plasma (LTP) ignition model has motivated much work in simulating these discharges from first principles. Most ignition models assume that an equilibrium plasma comprises the bulk of discharge kernels. LTP discharges, however, exhibit highly nonequilibrium behavior. In this work, a method to determine a consistent initialization of LTP discharge kernels for use in engine CFD codes like CONVERGE is proposed. The method utilizes first principles discharge simulations. Such an LTP kernel is introduced in a flammable mixture of air and fuel, and the subsequent plasma expansion and ignition simulation is carried out using a reacting flow solver with detailed chemistry. Finally, the proposed numerical approach is shown to produce results that agree with experimental observations regarding the ignitability of methane-air and ethylene-air mixtures by LTP discharges.
The PBX 9501 self-ignited thermally damaged explosive (SITX) test was carried out by Holmes et al. A 6-inch diameter sphere made of the PBX 9501 was heated inside a steel confinement vessel until it underwent a thermal explosion. This test provided measured temperature data that was used to define a boundary condition to carry out post-experiment finite element (FE) thermal simulations and comparisons of the FE thermal response with the experimental time to thermal ignition and internal temperature measurements at the center of the PBX 9501 sphere. Finite element models of this experiment were constructed using the heat transfer code Aria. A cookoff model for PBX 9501 was used to simulate the thermal decomposition of PBX 9501 and to predict the time to thermal ignition. An inverse heat conduction scheme was designed and applied to estimate a 1D temperature boundary condition.
In this work, a new Inverse Heat Conduction (IHC) algorithm is applied to estimate the surface temperatures at twelve locations on the inner and outer cylindrical boundaries of an HMX-based Plastic Bonded Explosive (PBX) annulus. This IHC algorithm was developed in references using a set of Direct Heat Conduction (DHC) solutions and a temperature correction method. The DHC solutions were calculated using a Galerkin based finite element (FE) method. This HMX based PBX annulus was used in the Large Scale Annular Cookoff (LSAC) experiment, Shot 5. The reason Shot 5 was chosen as a prototype mathematical model for this study is that the temperature was measured at eighteen locations in the midplane of the HMX-based PBX annulus. In addition, this annulus underwent an experimental thermal ignition and a deflagration that caused a thermal explosion and the disassembly of the experiment. The objective of this study is to describe how the application of the temperature correction algorithm produced the convergence of the DHC solutions to the measured temperatures at twelve internal locations in the midplane of the HMX-based PBX annulus.
An approach to the systematic optimization of a large-scale dynamic model is proposed which consists in parameterization of simulation results as response surfaces. The optimization procedure is carried out using a second-order orthogonal design. The approach proposed here is demonstrated by an example involving the shock-initiated ignition of methane.
An overview is given of an effort that focused on using CFD analysis to complement design and configuration definition of Lean-Direct Injection (LDI) combustion concepts for NASA's Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for second and third generation LDI configurations at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray modeling, ignition and kinetics modeling. Emissions (EINOx) characteristics were predicted for CST cruise conditions, and compared with emissions data from experimental measurements to evaluate the fidelity of the CFD modeling approach to predict emissions changes in response to changes in supersonic cycle conditions.