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.
Chemical kinetic model for hydrocarbon fuels combustion and application to multidimensional finite difference mixing analyses in hypersonic engines and nozzles
A CO-CH4-air reaction system is used to demonstrate the computational singular perturbation (CSP) method for deriving time-resolved simplified chemical kinetics models. CSP provides a programmable algorithm to group the given collection of elementary reactions into reaction groups which are ordered according to their speed. The concept of Importance Index k(m)exp s is introduced: k(m)exp s is defined to be a number between 0 and 1 which measures the importance of the m-th reaction group to the s-th reactant and can readily be computed from data generated by CSP. It is suggested that the robustness of the solutions of the reaction system can be qualitatively assessed by inspecting the Importance Index data.
The purpose of this research project is to develop and validate a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. Hydrazine is used extensively in aerospace propulsion, and although liquid hydrazine is not considered detonable, many fuel handling systems create multiphase mixtures of fuels and fuel vapors during their operation. Therefore, a thorough knowledge of the decomposition chemistry of hydrazine under a variety of conditions can be of value in assessing potential operational hazards in hydrazine fuel systems. To gain such knowledge, a reasonable starting point is the development and validation of a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. A reasonably complete mechanism was published in 1996, however, many of the elementary steps included had outdated rate expressions and a thorough investigation of the behavior of the mechanism under a variety of conditions was not presented. The current work has included substantial revision of the previously published mechanism, along with a more extensive examination of the decomposition behavior of hydrazine. An attempt to validate the mechanism against the limited experimental data available has been made and was moderately successful. Further computational and experimental research into the chemistry of this fuel needs to be completed.
The stoichiometric combustion of propane behind incident shock waves was studied experimentally and analytically over a temperature range from 1700 K to 2600 K and a pressure range from 1.2 to 1.9 atm. Measurements of the concentrations of carbon monoxide (CO) and carbon dioxide (CO2) and the product of the oxygen atom and carbon dioxide concentrations (O)(CO) were made after passage of the incident shock wave. A kinetic mechanism was developed which, when used in a computer program for a flowing, reacting gas behind an incident shock wave predicted experimentally measured results quite well. Ignition delay times from the literature were also predicted quite well. The kinetic mechanism consisted of 59 individual kinetic steps.
The oxidation of stoichiometric mixtures of benzene and toluene behind incident shock waves was studied for a temperature range from 1700 to 2800 K and a pressure range from 1.1 to 1.7 atm. The concentration of CO and CO2 produced were measured as well as the product of the oxygen atom and carbon monoxide concentrations. Comparisons between the benzene experimental data and results calculated by use of a reaction mechanism published in the open literature were carried out. With some additional reactions and changes in rate constants to reflect the pressure-temperature range of the experimental data, a good agreement was achieved between computed and experimental results. A reaction mechanism was developed for toluene oxidation based on analogous rate steps from the benzene mechanism. Measurements of NOx levels in an actual flame device, a jet-stirred combustor, were reproduced successfully by use of the reaction mechanism developed from the shock-tube experiments on toluene. These experimental measurements of NOx levels were reproduced from a computer simulation of a jet-stirred combustor.
The formation of nitrogen oxides NOx during combustion of methane, propane, and a jet fuel, JP-4, was investigated in a jet stirred combustor. The results of the experiments were interpreted using reaction models in which the nitric oxide (NO) forming reactions were coupled to the appropriate hydrocarbon combustion reaction mechanisms. Comparison between the experimental data and the model predictions reveals that the CH + N2 reaction process has a significant effect on NO formation especially in stoichiometric and fuel rich mixtures. Reaction models were assembled that predicted nitric oxide levels that were in reasonable agreement with the jet stirred combustor data and with data obtained from a high pressure (5.9 atm (0.6 MPa)), prevaporized, premixed, flame tube type combustor. The results also suggested that the behavior of hydrocarbon mixtures, like JP-4, may not be significantly different from that of pure hydrocarbons. Application of the propane combustion and nitric oxide formation model to the analysis of NOx emission data reported for various aircraft gas turbines showed the contribution of the various nitric oxide forming processes to the total NOx formed.
Explore the source record for details and available documents.
A unified overview is presented for chemical kinetics and chemical modeling in planetary atmospheres. The recent major advances in the understanding of the chemistry of the terrestrial atmosphere make the study of planets more interesting and relevant. A deeper understanding suggests that the important chemical cycles have a universal character that connects the different planets and ultimately link together the origin and evolution of the solar system. The completeness (or incompleteness) of the data base for chemical kinetics in planetary atmospheres will always be judged by comparison with that for the terrestrial atmosphere. In the latter case, the chemistry of H, O, N, and Cl species is well understood. S chemistry is poorly understood. In the atmospheres of Jovian planets and Titan, the C-H chemistry of simple species (containing 2 or less C atoms) is fairly well understood. The chemistry of higher hydrocarbons and the C-N, P-N chemistry is much less understood. In the atmosphere of Venus, the dominant chemistry is that of chlorine and sulfur, and very little is known about C1-S coupled chemistry. A new frontier for chemical kinetics both in the Earth and planetary atmospheres is the study of heterogeneous reactions. The formation of the ozone hole on Earth, the ubiquitous photochemical haze on Venus and in the Jovian planets and Titan all testify to the importance of heterogeneous reactions. It remains a challenge to connect the gas phase chemistry to the production of aerosols.
A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a tubular combustor is presented. A numerical procedure for the solution of the governing differential equations is described, and models for chemical equilibrium and chemical kinetics calculations are presented. The chemical equilibrium model is used to characterize the hydrocarbon reactions. The chemical kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor consists of a cylindrical duct of varying cross sections with concentric streams of gaseous fuel and air entering the duct at one end. Four sample cases with specified inlet and boundary conditions are considered, and the results are discussed
A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a two-concentric-tube combustor is presented. A numerical procedure for the solution of the governing differential equations is described and models for chemical-equilibrium and chemical-kinetics calculations are presented. The chemical-equilibrium model is used to characterize the hydrocarbon reactions. The chemical-kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor considered consists of two coaxial ducts. Concentric streams of gaseous fuel and air enter the inlet duct at one end; the flow then reverses and flows out through the outer duct. Two sample cases with specified inlet and boundary conditions are considered and the results are discussed.
The aerothermodynamics of thermochemical nonequilibrium carbon dioxide flows is studied. The chemical kinetics models of McKenzie and Park are implemented in separate three-dimensional computational fluid dynamics codes. The codes incorporate a five-species gas model characterized by a translational-rotational and a vibrational temperature. Solutions are obtained for flow over finite length elliptical and circular cylinders. The computed flowfields are then employed to calculate Mach-Zehnder interferograms for comparison with experimental data. The accuracy of the chemical kinetics models is determined through this comparison. Also, the methodology of the three-dimensional thermochemical nonequilibrium code is verified by the reproduction of the experiments.
A two-dimensional/axisymmetric program for the simulation of three-phase reactive flows is presented. The three phases are: a multiple-species gaseous phase, a single-species liquid phase, and a particulate droplet phase. The liquid and gaseous phases are described with continuous Eulerian equations, while the droplets are described in a discrete Lagrangian fashion. The three phases are fully coupled together. An atomization model, a multireaction chemical kinetics model, and a subgrid scale turbulence model completes the description of a general liquid/gas bipropellant combustion process. The code development criteria and qualitative features are highlighted.
A combustion analysis program aimed at upgrading and applying advanced computer programs for gas turbine applications is discussed. 2D and 3D codes, KIVA-II and LeRC-3D, have been used to provide insight into the combustion process and combustor design. The computations performed through these codes show their capability to produce reasonable results, despite such deficiencies in the current models as accurate chemical kinetics modeling of hydrocarbon combustion and turbulence and turbulence combustion interaction modeling.
A two-temperature chemical-kinetic model for air is assessed by comparing theoretical results with existing experimental data obtained in shock-tubes, ballistic ranges, and flight experiments. In the model, named the TTv model, one temperature (T) is assumed to characterize the heavy-particle translational and molecular rotational energies, and another temperature (Tv) to characterize the molecular vibrational, electron translational, and electronic excitation energies. The theoretical results for nonequilibrium air flow in shock tubes are obtained using the computer code STRAP (Shock-Tube Radiation Program), and for flow along the stagnation streamline in the shock layer over spherical bodies using the newly developed code STRAP (Stagnation-Point Radiation Program). Substantial agreement is shown between the theoretical and experimental results for relaxation times and radiative heat fluxes. At very high temperatures the spectral calculations need further improvement. The present agreement provides strong evidence that the two-temperature model characterizes principal features of nonequilibrium air flow. New theoretical results using the model are presented for the radiative heat fluxes at the stagnation point of a 6-m-radius sphere, representing an aeroassisted orbital transfer vehicle, over a range of free-stream conditions. Assumptions, approximations, and limitations of the model are discussed.
A chemical kinetic model is developed to help understand and optimize the production of single-walled carbon nanotubes via the high-pressure carbon monoxide (HiPco) process, which employs iron pentacarbonyl as the catalyst precursor and carbon monoxide as the carbon feedstock. The model separates the HiPco process into three steps, precursor decomposition, catalyst growth and evaporation, and carbon nanotube production resulting from the catalyst-enhanced disproportionation of carbon monoxide, known as the Boudouard reaction: 2 CO(g)-->C(s) + CO2(g). The resulting detailed model contains 971 species and 1948 chemical reactions. A second model with a reduced reaction set containing 14 species and 22 chemical reactions is developed on the basis of the detailed model and reproduces the chemistry of the major species. Results showing the parametric dependence of temperature, total pressure, and initial precursor partial pressures are presented, with comparison between the two models. The reduced model is more amenable to coupled reacting flow-field simulations, presented in the following article.
Opposed Jet Burner (OJB) tools have been used extensively by the authors to measure Flame Strength (FS) extinction limits of laminar H2/N2 air and (recently) hydrocarbon (HC) air Counterflow Diffusion Flames (CFDFs) at one atm. This paper details normalization of FSs of N2- diluted H2 and HC systems to account for effects of fuel composition, temperature, pressure, jet diameter, inflow Reynolds number, and inflow velocity profile (plug, contoured nozzle; and parabolic, straight tube). Normalized results exemplify a sensitive accurate means of validating, globally, reduced chemical kinetic models at approx. 1 atm and the relatively low temperatures approximating the loss of non-premixed idealized flameholding, e.g., in scramjet combustors. Laminar FS is defined locally as maximum air input velocity, U(sub air), that sustains combustion of a counter-jet of g-fuel at extinction. It uniquely characterizes a fuel. And global axial strain rate at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or (sub t)) can be compared directly with computed extinction limits, determined using either a 1-D Navier Stokes stream-function solution, using detailed transport and finite rate chemistry, or (better yet) a detailed 2-D Navier Stokes numerical simulation. The experimental results define an idealized flameholding reactivity scale that shows wide ranging (50 x) normalized FS s for various vaporized-liquid and gaseous HCs, including, in ascending order: JP-10, methane, JP-7, n-heptane, n-butane, propane, ethane, and ethylene. Results from H2 air produce a unique and exceptionally strong flame that agree within approx. 1% of a recent 2-D numerically simulated FS for a 3 mm tube-OJB. Thus we suggest that experimental FS s and/or FS ratios, for various neat and blended HCs w/ and w/o additives, offer accurate global tests of chemical kinetic models at the Ts and Ps of extinction. In conclusion, we argue the FS approach is more direct and fundamental, for assessing, e.g., idealized scramjet flameholding potentials, than measurements of laminar burning velocity or blowout in a Perfectly Stirred Reactor, because the latter characterize premixed combustion in the absence of aerodynamic strain. And FS directly measures a chemical kinetic characteristic of non-premixed combustion at typical flameholding temperatures. It mimics conditions where gfuels are typically injected into a subsonic flameholding recirculation zone that captures air, where the effects of aerodynamic strain and associated multi-component diffusion become important.
The testing and development of existing global and detailed chemical kinetic models for soot formation requires measurements of soot and radical concentrations in flames. A clearer understanding of soot particle inception relies upon the evaluation and refinement of these models in comparison with such measurements. We present measurements of soot formation and hydroxyl (OH) concentration in sequences of flat premixed atmospheric-pressure C2H4/O2/N2 flames and 80-torr C2H4/O2 flames for a unique range of equivalence ratios bracketting the critical equivalence ratio (phi(sub c)) and extending to more heavily sooting conditions. Soot volume fraction and number density profiles are measured using a laser scattering-extinction apparatus capable of resolving a 0.1 percent absorption. Hydroxyl number density profiles are measured using laser-induced fluorescence (LIF) with broadband detection. Temperature profiles are obtained from Rayleigh scattering measurements. The relative volume fraction and number density profiles of the richer sooting flames exhibit the expected trends in soot formation. In near-phi(sub c) visibility sooting flames, particle scattering and extinction are not detected, but an LIF signal due to polycyclic aromatic hydrocarbons (PAH's) can be detected upon excitation with an argon-ion laser. A linear correlation between the argon-ion LIF and the soot volume fraction implies a common mechanistic source for the growth of PAH's and soot particles. The peak OH number density in both the atmospheric and 80-torr flames declines with increasing equivalence ratio, but the profile shape remains unchanged in the transition to sooting, implying that the primary reaction pathways for OH remain unchanged over this transition. Chemical kinetic modeling is demonstrated by comparing predictions using two current reaction mechanisms with the atmospheric flame data. The measured and predicted OH number density profiles show good agreement. The predicted benzene number density profiles correlate with the measured trends in soot formation, although anomalies in the benzene profiles for the richer and cooler sooting flames suggest a need for the inclusion of benzene oxidation reactions.
Chemical kinetic models for the nucleation and growth of clusters and single-walled carbon nanotube (SWNT) growth are developed for numerical simulations of the production of SWNTs. Two models that involve evaporation and condensation of carbon and metal catalysts, a full model involving all carbon clusters up to C80, and a reduced model are discussed. The full model is based on a fullerene model, but nickel and carbon/nickel cluster reactions are added to form SWNTs from soot and fullerenes. The full model has a large number of species--so large that to incorporate them into a flow field computation for simulating laser ablation and arc processes requires that they be simplified. The model is reduced by defining large clusters that represent many various sized clusters. Comparisons are given between these models for cases that may be applicable to arc and laser ablation production. Solutions to the system of chemical rate equations of these models for a ramped temperature profile show that production of various species, including SWNTs, agree to within about 50% for a fast ramp, and within 10% for a slower temperature decay time.