New considerations on causes for combustion instability in liquid propellant rockets
Liquid propellant rocket combustion instability, examining droplet vaporization, wake burning and chemical kinetics
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.
Liquid propellant rocket combustion instability, examining droplet vaporization, wake burning and chemical kinetics
Thermodynamic data are required input for the finite kinetics and equilibrium computer programs needed for modeling the combustion of hydrocarbons in the fields of energy and pollution research. Least squares determined coefficients of the curve-fitted thermodynamic data for 193 species in the C-H-O-N system are presented in card image form and are of suitable format for use by common computer programs.
Propulsion engine combustor design and analysis requires experimentally verified data on the chemical kinetics of fuel. Among the important data is the combustion extinction limit as measured by observed maximum flame strain rate. The extinction limit relates to the ability to maintain a flame in a combustor during operation. Extinction limit data can be obtained for a given fuel by means of a laminar flame experiment using an opposed jet burner (OJB). Laminar extinction limit data can be applied to the turbulent application of a combustor via laminar flamelet modeling. The OJB consists of two axi-symmetric tubes (one for fuel and one for oxidizer), which produce a flat, disk-like counter-flow diffusion flame. This paper presents results of experiments to measure extinction limits for n-heptane and the military specification fuel JP-7, obtained from an OJB. JP-7 is an Air Force-developed fuel that continues to be important in the area of hypersonics. Because of its distinct properties it is currently the hydrocarbon fuel of choice for use in Scramjet engines. This study provides much-desired data for JP-7, for which very little information previously existed. The interest in n-heptane is twofold. First, there has been a significant amount of previous extinction limit study and resulting data with this fuel. Second, n-heptane (C7H16) is a pure substance, and therefore does not vary in composition as does JP-7, which is a mixture of several different hydrocarbons. These two facts allow for a baseline to be established by comparing the new OJB results to those previously taken. Additionally, the data set for n-heptane, which previously existed for mixtures up to 26 mole percent in nitrogen, is completed up to 100% n-heptane. The extinction limit data for the two fuels are compared, and complete experimental results are included.
Earth to orbit launch vehicles are propelled by rocket engines and motors, both liquid and solid. This chapter will discuss liquid engines. The heart of a launch vehicle is its engine. The remainder of the vehicle (with the notable exceptions of the payload and guidance system) is an aero structure to support the propellant tanks which provide the fuel and oxidizer to feed the engine or engines. The basic principle behind a rocket engine is straightforward. The engine is a means to convert potential thermochemical energy of one or more propellants into exhaust jet kinetic energy. Fuel and oxidizer are burned in a combustion chamber where they create hot gases under high pressure. These hot gases are allowed to expand through a nozzle. The molecules of hot gas are first constricted by the throat of the nozzle (de-Laval nozzle) which forces them to accelerate; then as the nozzle flares outwards, they expand and further accelerate. It is the mass of the combustion gases times their velocity, reacting against the walls of the combustion chamber and nozzle, which produce thrust according to Newton's third law: for every action there is an equal and opposite reaction. Solid rocket motors are cheaper to manufacture and offer good values for their cost. Liquid propellant engines offer higher performance, that is, they deliver greater thrust per unit weight of propellant burned. They also have a considerably higher thrust to weigh ratio. Since liquid rocket engines can be tested several times before flight, they have the capability to be more reliable, and their ability to shut down once started provides an extra margin of safety. Liquid propellant engines also can be designed with restart capability to provide orbital maneuvering capability. In some instances, liquid engines also can be designed to be reusable. On the solid side, hybrid solid motors also have been developed with the capability to stop and restart. Solid motors are covered in detail in chapter 11. Liquid rocket engine operational factors can be described in terms of extremes: temperatures ranging from that of liquid hydrogen (-423 F) to 6000 F hot gases; enormous thermal shock (7000 F/sec); large temperature differentials between contiguous components; reactive propellants; extreme acoustic environments; high rotational speeds for turbo machinery and extreme power densities. These factors place great demands on materials selection and each must be dealt with while maintaining an engine of the lightest possible weight. This chapter will describe the design considerations for the materials used in the various components of liquid rocket engines and provide examples of usage and experiences in each.
Combustion synthesis (self propagating, high temperature synthesis-SHS) is a novel technique that is capable of producing many advanced materials. The ignition temperature (Tig) of such combustion synthesis reactions is often coincident with that of the lowest melting point reactant. The resultant liquid metal wets and spreads around the other solid reactant particles of higher melting points, thereby improving the reactant contact and kinetics, followed by formation of the required compounds. This ignition initiates a combustion propagating wave whose narrow reaction front rapidly travels through the reactants. Since this process is highly exothermic, the heat released by combustion often melts the reactant particles ahead of the combustion front and ignites the adjacent reactant layer, resulting in a self-sustaining reaction. Whenever a fluid phase (liquid or gas) is generated by the reaction system, gravity-driven phenomena can occur. Such phenomena include convective flows of fluid by conventional or unstable convection and settling of the higher density phases. A combustion process is often associated with various kinds of fluid flow. For instance, if the SHS reaction is carried out under inert or reactive gas atmospheres, or a volatile, e.g., B2O3, is deliberately introduced as a reactant, convective flows of the gas will occur due to a temperature gradient existing in the atmosphere when a combustion wave is initiated. The increased gas flow will produce a porous (or expanded) SHS product. Owing to the highly exothermic nature of many SHS reactions, liquid phase(s) can also form before, at, or after the combustion front. The huge temperature gradient at the combustion front can induce convective flows (conventional or unstable) of the liquid phase. Each of these types of convective fluid flow can change the combustion behavior of the synthesizing reaction, and, therefore, the resultant product microstructure. In addition, when two or more phases of different density are produced at or ahead of the propagating combustion front settling of the higher density phase will occur resulting in a non-uniform product microstructure and properties.
SiC and Si3N4 materials were tested under various turbine engine combustion environments, chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high speed aircraft. Representative CVD, sintered, and composite materials were evaluated in both furnace and high pressure burner rig exposure. While protective SiO2 scales form in all cases, evidence is presented to support paralinear growth kinetics, i.e. parabolic growth moderated simultaneously by linear volatilization. The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were used to map SiO2 volatility (and SiC recession) over a range of temperature, pressure, and velocity. The functional dependency of material recession (volatility) that emerged followed the form: exp(-QIRT) * P(exp x) * v(exp y). These empirical relations were compared to rates predicted from the thermodynamics of volatile SiO and SiO(sub x)H(sub Y) reaction products and a kinetic model of diffusion through a moving, boundary layer. For typical combustion conditions, recession of 0.2 to 2 micron/h is predicted at 1200- 1400C, far in excess of acceptable long term limits.
High-pressure direct fuel injection plays the most crucial role in energy conversion and improving engine combustion efficiency and emission. The optimization of turbulent and multiphase fuel injection has focused on controlling hydrodynamic parameters such as injection pressure. While the thermodynamic influence is often considered in the flash boiling situation, we inquire into how gasoline-type fuel’s hydro- and thermodynamic properties impact the injection dynamics by fuel-temperature induced cavitation. The turbulent and cavitating flows emanating from the direct-injection nozzle is visualized by ultrafast x-ray imaging with an unprecedented spatiotemporal resolution. The ultrafast liquid-fuel dynamics are dominated by injection pressure as well as fuel temperature through cavitation, an important thermodynamic parameter but often difficult to control in engine combustion. With the most direct and quantitative measurement, we discovered that the near-nozzle fuel-jet dynamics can be perfectly scaled by a single dimensionless parameter, cavitation number, particularly sensitive to the fuel temperature, in a wide operation range. This universal scaling shows that cavitation can be harnessed to elevate the pneumatic-hydraulic to kinetic energy conversion efficiency, critical for promoting fuel atomization and engine combustion performance. This enhancement effect will have even more impact on engine combustion using alternative low-emission fuels with higher saturated vapor pressure.
Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.
Silicon carbide (SiC) and Si3N4 materials were tested in various turbine engine combustion environments chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high-speed aircraft. Representative chemical vapor-deposited (CVD), sintered, and composite materials were evaluated by furnace and high-pressure burner rig exposures. Although protective SiO2 scales formed in all cases, the evidence presented supports a model based on paralinear growth kinetics (i.e., parabolic growth moderated simultaneously by linear volatilization). The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were thus used to map SiO2 volatility (and SiC recession) over a range of temperatures, pressures, and velocities. The functional dependency of material recession (volatility) that emerged followed the form A[exp(-Q / RT)](P(sup x)v(sup y). These empirical relations were compared with rates predicted from the thermodynamics of volatile SiO and SiOxHy reaction products and a kinetic model of diffusion through a moving boundary layer. For typical combustion conditions, recession of 0.2 to 2 micrometers/hr is predicted at 1200 to 1400 C, far in excess of acceptable long-term limits.
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.
The structures of premixed flames propagating in combustible systems, containing uniformly distributed volatile fuel particles, in an oxidizing gas mixtures is analyzed. The experimental results show that steady flame propagation occurs even if the initial equivalence ratio of the combustible mixture based on the gaseous fuel available in the particles, phi(u) is substantially larger than unity. A model is developed to explain these experimental observations. In the model it is presumed that the fuel particles vaporize first to yield a gaseous fuel of known chemical composition which then reacts with oxygen in a one-step overall process. It is shown that the interplay of vaporization kinetics and oxidation process, can result in steady flame propagation in combustible mixtures where the value of phi(u) is substantially larger than unity. This prediction is in agreement with experimental observations.
2,4,dimethyloxetane is an important cyclic ether intermediate that is produced from hydroperoxyalkyl (QOOH) radicals in the low-temperature combustion of n-pentane. However, the reaction mechanisms and rates of consumption pathways remain unclear. In the present work, the pressure- and temperature-dependent kinetics of seven cyclic ether peroxy radicals, which stem from 2,4,dimethyloxetane via H-abstraction and O 2 addition, were determined. The automated kinetic workflow code, KinBot, was used to model the complexity of the chemistry in a stereochemically resolved manner and solve the resulting master equations from 300–1000 K and from 0.01–100 atm. The main conclusions from the calculations include (i) diastereomeric cyclic ether peroxy radicals show significantly different reactivities, (ii) the stereochemistry of the peroxy radical determines which QOOH isomerization steps are possible, (iii) conventional QOOH decomposition pathways, such as cyclic ether formation and HO 2 elimination, compete with ring-opening reactions, which primarily produce OH radicals, the outcome of which is sensitive to stereochemistry. Ring-opening reactions lead to unique products, such as unsaturated, acyclic peroxy radicals, that form direct connections with species present in other chemical kinetics mechanisms through "cross-over" reactions that may complicate the interpretation of experimental results from combustion of n-pentane and, by extension, other alkanes. For example, one cross-over reaction involving 1-hydroperoxy-4-pentanone-2-yl produces 2-(hydroperoxymethyl)-3-butanone-1-yl, which is an iso-pentane-derived ketohydroperoxide (KHP). At atmospheric pressure, the rate of chemical reactions of all seven peroxy radicals compete with that of collisional stabilization, resulting in well-skipping reactions. However, at 100 atm, only one out of seven peroxy radicals undergoes significant well-skipping reactions. Here, the rates produced from the master equation calculations provide the first foundation for the development of detailed sub-mechanisms for cyclic ether intermediates. In addition, analysis of the complex reaction mechanisms of 2,4-dimethyloxetane-derived peroxy radicals provides insights into the effects of stereoisomers on reaction pathways and product yields.
Here, a detailed experimental study of the low-to-intermediate temperature combustion of methylbutene isomers, i.e., branched C 5 alkenes, has been undertaken with multiple experimental facilities. Ignition delay times were measured at equivalence ratios 0.5–2.0, 685–1020 K and up to 45 bar condition from two rapid compression machines and showed slight deviation from an Arrhenius behavior for all three isomers, while their reactivity order differs as temperature changes. Sampled intermediates formed during the oxidation process of mixtures at 900–1150 K and 0.82 bar from a flow reactor and at 730 K and 20 bar from a rapid compression machine were analyzed using gas chromatography techniques. Trends in the formation and consumption of sampled intermediates were modeled using a kinetic model developed in this work for all three isomers. Rate of production and sensitivity analyses emphasize the role of double bond-specific reactions governing the global reactivity of these fuels. Additional studies of the addition reactions of HO 2 radicals to the double bond and to allylic radicals may improve the model performance.
The structure of premixed flames propagating in combustible systems containing uniformly distributed volatile fuel particles in an oxidizing gas mixture is analyzed. This analysis is motivated by experiments conducted at NASA Lewis Research Center on the structure of flames propagating in combustible mixtures of lycopodium particles and air. Several interesting modes of flame propagation were observed in these experiments depending on the number density and the initial size of the fuel particle. The experimental results show that steady flame propagation occurs even if the initial equivalence ratio of the combustible mixture based on the gaseous fuel available in the particles, phi sub u, is substantially larger than unity. A model is developed to explain these experimental observations. In the model, it is presumed that the fuel particles vaporize first to yield a gaseous fuel of known chemical composition which then reacts with oxygen in a one-step overall process. The activation energy of the chemical reaction is presumed to be large. The activation energy characterizing the kinetics of vaporization is also presumed to be large. The equations governing the structure of the flame were integrated numerically. It is shown that the interplay of vaporization kinetics and oxidation process can result in steady flame propagation in combustible mixtures where the value of phi sub u is substantially larger than unity. This prediction is in agreement with experimental observations.
The flow field, spray penetration, and combustion in two-stroke diesel engines are described. Fuel injection begins at 345 degrees after top dead center (ATDC) and n-dodecane is used as the liquid fuel. Arrhenius kinetics is used to calculate the reaction rate term in the quasi-global combustion model. When the temperature, fuel, and oxygen mass fraction are within suitable flammability limits, combustion begins spontaneously. No spark is necessary to ignite a localized high temperature region. Compression is sufficient to increase the gaseous phase temperature to a point where spontaneous chemical reactions occur. Results are described for a swirl angle of 22.5 degrees.
Uni- and bi-molecular reactions involving Criegee intermediates (CIs) have been the focus of many studies due to the role these molecules play in atmospheric chemistry. The reactivity of CIs is known to strongly depend on their structure. The reaction network of the second simplest CI, acetaldehyde oxide (CH 3 CHOO), is investigated in this work in an atmospheric pressure jet-stirred reactor (JSR) during the ozonolysis of trans-2-butene to explore the kinetic pathways relevant to atmospheric chemistry and low-temperature combustion. The mole fraction profiles of reactants, intermediates, and final products are determined by means of molecular-beam mass spectrometry in conjunction with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. A network of CI reactions is identified in the temperature region below 600 K, characterized by CI addition to trans-2-butene, water, formaldehyde, formic acid, and methanol. No sequential additions of the CH 3 CHOO CI are observed, in contrast with the reactivity of the simplest CI (H 2 COO) and the earlier observation of an extensive reaction network with up to four H 2 COO sequential additions (Phys. Chem. Chem. Phys., 2019, 21, 7341–7357). Experimental photoionization efficiency scans recorded at 300 K and 425 K and ab initio threshold energy calculations lead to the identification and quantification of previously elusive intermediates, such as ketohydroperoxide and hydroperoxide species. Specifically, the C 4 H 8 + O 3 adduct is identified as a ketohydroperoxide (KHP, 3-hydroperoxybutan-2-one, CH 3 C([double bond, length as m-dash]O)CH(CH 3 )OOH), while hydroxyacetaldehyde (glycolaldehyde, HCOCH 2 OH) formation is attributed to unimolecular isomerization of the CIs. Other hydroperoxide species such as methyl hydroperoxide (CH 3 OOH), ethyl hydroperoxide (C 2 H 5 OOH), butyl hydroperoxide (OOH), hydroperoxyl acetaldehyde (HOOCH 2 CHO), hydroxyethyl hydroperoxide (CH 3 CH(OH)OOH), but-1-enyl-3-hydroperoxide, and 4-hydroxy-3-methylpentan-2-one (HOCH(CH 3 )CH(CH 3 )C([double bond, length as m-dash]O)CH 3 ) are also identified. Detection of additional oxygenated species such as methanol, ethanol, ketene, and aldehydes suggests multiple active oxidation routes. Finally, these results provide additional evidence that CIs are key intermediates of the ozone-unsaturated hydrocarbon reactions providing critical inputs for improved kinetics models.
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. 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 delay time. NG/NH3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis, demonstrating hydrogen chemistry to be dominant in NG/H2 mixtures and hydrocarbon chemistry to be dominant in NG/NH3 mixtures. 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.