Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical kinetics”

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 181 records · Page 10

AdapChem

AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.

Oluwole, Oluwayemisi O.↗

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

Rate constants for reactions of ClO/x/ of atmospheric interest

Chemical kinetics measurements on 82 gas phase reactions of chlorine containing species are reviewed. Recommended rate constants are given. The principal species of interest are Cl, Cl2, ClO, Cl2O, ClOO, OClO, CINO, HCl and halo derivatives of methane and ethane. Absorption spectra are given for 21 species. In addition the chemical kinetics methods used to obtain these data are discussed with regard to their applicability and reliability.

Watson, R. T.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reaction flows

The use of silane (SiH4) as an effective ignitor and flame stabilizing pilot fuel is well documented. A reliable chemical kinetic mechanism for prediction of its behavior at the conditions encountered in the combustor of a SCRAMJET engine was calculated. The effects of hydrogen addition on hydrocarbon ignition and flame stabilization as a means for reduction of lengthy ignition delays and reaction times were studied. The ranges of applicability of chemical kinetic models of hydrogen-air combustors were also investigated. The CHARNAL computer code was applied to the turbulent reaction rate modeling.

Magnotti, F.↗

Heterogeneous Combustion of Porous Graphite Particles in Microgravity

Recent theoretical investigations on graphite particle combustion have employed several levels of heterogeneous reaction models, ranging from global to elementary models, to describe the oxidation of carbon to gaseous products. Unlike the counterpart homogeneous reaction models, these heterogeneous reaction models are not well developed because of the difficulties associated with decoupling the physical characteristics of the solid (e.g. surface area taking part in combustion) from the chemical kinetic data. This is certainly true for porous graphite particle combustion, where heterogeneous and homogeneous reactions occur within the pores and play an important role in the overall oxidation process. As a result, there are considerable uncertainties of physical phenomena predicted using different heterogeneous kinetic models available in the literature. A good example, discussed later in this paper, is the predicted critical particle size below which the mass burning rate becomes exponentially small. The main goal of this study is to understand the basic mechanism controlling such rapid changes in burning rates, by developing a model where physical contributions are decoupled from chemical rate constants in a consistent manner. Another important goal of the proposed study is to develop a truly intrinsic, detailed heterogeneous reaction model for porous graphite combustion at high-temperatures, and to derive a systematically reduced heterogeneous reaction model in terms of the elementary reaction rate constants of the detailed model. The validation of chemical kinetic models describing the heterogeneous and homogeneous combustion in and around a spherically symmetric porous graphite particle can be considerably simplified by experimental measurements obtained under microgravity conditions. A vital component of this study is to conduct such supporting experiments on particle burning rate and surface temperature using NASA microgravity facilities, in close coordination with the theoretical effort. The basic understanding obtained and models developed as part of this project will be useful for optimal design of coal combustion devices. These models can also be extended to investigate the role of heterogeneous chemistry on pollutant formation pathways in combustion devices. The theoretical approach developed here, with pore diffusion effects decoupled from the chemical effects, can also be extended to understand the heterogeneous combustion of other porous fuels, for example, combustion of magnesium in a CO2 environment for propulsion in the Martian atmosphere.

Chelliah, Harsha K.↗

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

Evaluation of the Uncertainty in JP-7 Kinetics Models Applied to Scramjets

One of the challenges of designing and flying a scramjet-powered vehicle is the difficulty of preflight testing. Ground tests at realistic flight conditions introduce several sources of uncertainty to the flow that must be addressed. For example, the scales of the available facilities limit the size of vehicles that can be tested and so performance metrics for larger flight vehicles must be extrapolated from ground tests at smaller scales. To create the correct flow enthalpy for higher Mach number flows, most tunnels use a heater that introduces vitiates into the flow. At these conditions, the effects of the vitiates on the combustion process is of particular interest to the engine designer, where the ground test results must be extrapolated to flight conditions. In this paper, the uncertainty of the cracked JP-7 chemical kinetics used in the modeling of a hydrocarbon-fueled scramjet was investigated. The factors that were identified as contributing to uncertainty in the combustion process were the level of flow vitiation, the uncertainty of the kinetic model coefficients and the variation of flow properties between ground testing and flight. The method employed was to run simulations of small, unit problems and identify which variables were the principal sources of uncertainty for the mixture temperature. Then using this resulting subset of all the variables, the effects of the uncertainty caused by the chemical kinetics on a representative scramjet flow-path for both vitiated (ground) and nonvitiated (flight) flows were investigated. The simulations showed that only a few of the kinetic rate equations contribute to the uncertainty in the unit problem results, and when applied to the representative scramjet flowpath, the resulting temperature variability was on the order of 100 K. Both the vitiated and clean air results showed very similar levels of uncertainty, and the difference between the mean properties were generally within the range of uncertainty predicted.

Norris, A. T.↗

Analysis of thermo-chemical nonequilibrium models for carbon dioxide flows

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.

Rock, Stacey G.↗

Dynamics of flames near the rich-flammability limit of hydrogen-air mixtures

Flames near the rich-flammability limit of hydrogen-air mixtures are studied using a detailed, time-dependent, one-dimensional Lagrangian model. Results from the numerical simulations indicate that a steady burning velocity is not obtained for very rich hydrogen-air mixtures. As the amount of hydrogen is increased, a damped oscillation is observed in the flame and burning velocities, and then, with further increase in the amount of hydrogen, an undamped oscillation with a complex set of frequencies is observed. Simulations with a simplified one-step irreversible chemical reaction do not show these oscillations, suggesting that chemical kinetics plays a strong role in inducing these oscillations. Further analysis shows that the oscillations are due to a competition for H atoms between chain branching and chain-terminating reactions. Simulations of spherically expanding flames suggest that stretch effects (due to curvature) will cause the oscillations to occur in less rich mixtures than that observed for planar flames. The implications of these oscillations on the rich-flammability limit as well as the role of chemical kinetics in creating a fundamental flammability limit is discussed.

Kailasanath, K.↗

Development of a model and computer code to describe solar grade silicon production processes

Two computer codes were developed for describing flow reactors in which high purity, solar grade silicon is produced via reduction of gaseous silicon halides. The first is the CHEMPART code, an axisymmetric, marching code which treats two phase flows with models describing detailed gas-phase chemical kinetics, particle formation, and particle growth. It can be used to described flow reactors in which reactants, mix, react, and form a particulate phase. Detailed radial gas-phase composition, temperature, velocity, and particle size distribution profiles are computed. Also, deposition of heat, momentum, and mass (either particulate or vapor) on reactor walls is described. The second code is a modified version of the GENMIX boundary layer code which is used to compute rates of heat, momentum, and mass transfer to the reactor walls. This code lacks the detailed chemical kinetics and particle handling features of the CHEMPART code but has the virtue of running much more rapidly than CHEMPART, while treating the phenomena occurring in the boundary layer in more detail.

Gould, R. K.↗

Future fundamental combustion research for aeropropulsion systems

Physical fluid mechanics, heat transfer, and chemical kinetic processes which occur in the combustion chamber of aeropropulsion systems were investigated. With the component requirements becoming more severe for future engines, the current design methodology needs the new tools to obtain the optimum configuration in a reasonable design and development cycle. Research efforts in the last few years were encouraging but to achieve these benefits research is required into the fundamental aerothermodynamic processes of combustion. It is recommended that research continues in the areas of flame stabilization, combustor aerodynamics, heat transfer, multiphase flow and atomization, turbulent reacting flows, and chemical kinetics. Associated with each of these engineering sciences is the need for research into computational methods to accurately describe and predict these complex physical processes. Research needs in each of these areas are highlighted.

Mularz, E. J.↗

Future fundamental combustion research for aeropropulsion systems

Physical fluid mechanics, heat transfer, and chemical kinetic processes which occur in the combustion chamber of aeropropulsion systems were investigated. With the component requirements becoming more severe for future engines, the current design methodology needs the new tools to obtain the optimum configuration in a reasonable design and development cycle. Research efforts in the last few years were encouraging but to achieve these benefits research is required into the fundamental aerothermodynamic processes of combustion. It is recommended that research continues in the areas of flame stabilization, combustor aerodynamics, heat transfer, multiphase flow and atomization, turbulent reacting flows, and chemical kinetics. Associated with each of these engineering sciences is the need for research into computational methods to accurately describe and predict these complex physical processes. Research needs in each of these areas are highlighted.

Mularz, E. J.↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

An analytical study of the hydrogen-air reaction mechanism with application to scramjet combustion

A chemical kinetic mechanism for the combustion of hydrogen has been assembled and optimized by comparing the observed behavior as determined in shock tube and flame studies with that predicted by the mechanism. The reactions contained in the mechanism reflect the current state of knowledge of the chemistry of the hydrogen/air system, and the assigned rate coefficients are consistent with accepted values. It was determined that the mechanism is capable of satisfactorily reproducing the experimental results for a range of conditions relevant to scramjet combustion. Calculations made with the reaction mechanism for representative scramjet combustor conditions at Mach 8, 16, and 25 showed that chemical kinetic effects can be important and that combustor models which use nonequilibrium chemistry should be used in preference to models that assume equilibrium chemistry. For the conditions examined the results also showed the importance of including the HO2 chemistry in the mechanism. For Mach numbers less than 16, the studies suggest that an ignition source will most likely be required to overcome slow ignition chemistry. At Mach 25, the initial temperature and pressure was high enough that ignition was rapid and the presence of an ignition source did not significantly affect reaction rates.

Jachimowski, Casimir J.↗

Opposed-flow Spreading Flames: Effect of Sub-atmospheric Pressure on Spread and Burning Rates

Flame spread over solid fuels is a canonical problem in fire science, due to its direct implications on material flammability and importance in fire development. Flames in a microgravity environment can behave very differently than on Earth, posing additional risks for spaceflight life safety. Sub-atmospheric pressures in ground-based experiments have been proposed to approximately replicate the burning behavior of solid fuels in reduced gravity conditions because of similar effects on heat and mass transfer mechanisms from the flame to the solid. In opposed-flow flame spread, the solid fuel is heated by the flame ahead of its leading edge, and this process is strongly affected by the ambient conditions. In this work, we consider flames spreading over flat thin acrylic samples exposed to a forced flow of 20 cm/s, and pressures between 30 and 100 kPa. When radiation losses and chemical kinetic effects are neglected, experimental and theoretical results suggest an independence of the spread rate on pressure. However, the spread rate shows a decreasing trend with lower pressure in the range 30-60 kPa. The mass burning rate, calculated from the samples weight measured before and after the experiments, shows a similar trend. The flame length, on the other hand, shows a non-monotonic behavior, reaching a maximum at 60 kPa. Additionally, gas emissions measured during the experiments are used to estimate the heat release rate of the spreading flames. The comparison of the heat release rate obtained from the measured emissions and the estimated mass burning rates, suggests that chemical kinetics is not driving the decrease in flame spread observed at low pressures, whereas radiation losses could be responsible in analogy to what has been observed in micro and partial gravity studies. These results could provide more information to guide future Earth-based material flammability testing for spacecraft applications.

fire↗

Opposed Jet Burner Extinction Limits: Simple Mixed Hydrocarbon Scramjet Fuels vs Air

Opposed Jet Burner tools have been used extensively by the authors to measure Flame Strength (FS) of laminar non-premixed H2 air and simple hydrocarbon (HC) air counterflow diffusion flames at 1-atm. FS represents a strain-induced extinction limit based on air jet velocity. This paper follows AIAA-2006-5223, and provides new HC air FSs for global testing of chemical kinetics, and for characterizing idealized flameholding potentials during early scramjet-like combustion. Previous FS data included six HCs, pure and N2-diluted; and three HC-diluted H2 fuels, where FS decayed very nonlinearly as HC was added to H2, due to H-atom scavenging. This study presents FSs on mixtures of (candidate surrogate) HCs, some with very high FS ethylene. Included are four binary gaseous systems at 300 K, and a hot ternary system at approx. 600 K. The binaries are methane + ethylene, ethane + ethylene, methane + ethane, and methane + propylene. The first three also form two ternary systems. The hot ternary includes both 10.8 and 21.3 mole % vaporized n-heptane and full ranges of methane + ethylene. Normalized FS data provide accurate means of (1) validating, globally, chemical kinetics for extinction of non-premixed flames, and (2) estimating (scaling by HC) the loss of incipient flameholding in scramjet combustors. The n-heptane is part of a proposed baseline simulant (10 mole % with 30% methane + 60% ethylene) that mimics the ignition of endothermically cracked JP-7 like kerosene fuel, as suggested by Colket and Spadaccini in 2001 in their shock tube Scramjet Fuels Autoignition Study. Presently, we use FS to gauge idealized flameholding, and define HC surrogates. First, FS was characterized for hot nheptane + methane + ethylene; then a hot 36 mole % methane + 64% ethylene surrogate was defined that mimics FS of the baseline simulant system. A similar hot ethane + ethylene surrogate can also be defined, but it has lower vapor pressure at 300 K, and thus exhibits reduced gaseous capacity. The new FS results refine our earlier idealized reactivity scale that shows wide ranging (50 x) diameter-normalized FSs for various HCs. These range from JP-10 and methane to H2 air, which produces an exceptionally strong flame that agrees within approx. 1% of recent 2-D numerically simulations. Finally, we continue advocating the FS approach as more direct and fundamental, for assessing idealized scramjet flameholding potentials, than measurements of unstrained laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

Ignition of mixtures of SiH sub 4, CH sub 4, O sub 2, and Ar or N sub 2 behind reflected shock waves

Ignition delay times in mixtures of methane, silane, and oxygen diluted with argon and nitrogen were measured behind reflected shock waves generated in the chemical kinetic shock tube at Langley Research Center. The delay times were inferred from the rapid increase in pressure that occurs at ignition, and the ignition of methane was verified from the emission of infrared radiation from carbon dioxide. Pressures of 1.25 atm and temperatures from 1100 K to 1300 K were generated behind the reflected shocks; these levels are representative of those occurring within a supersonic Ramjet combustor. Expressions for the ignition delay time as a function of temperature were obtained from least squares curve fits to the data for overall equivalence ratios of 0.7 and 1.0. The ignition delay times with argon as the diluent were longer than those with nitrogen as the diluent. The infrared wavelength observations at 4.38 microns for carbon dioxide indicated that silane and methane ignited simultaneously (i.e., within the time resolution of the measurement). A combined chemical kinetic mechanism for mixtures of silane, methane, oxygen, and argon or nitrogen was assembled from one mechanism that accurately predicted the ignition of methane and a second mechanism that accurately predicted silane hydrogen ignition. Comparisons between this combined mechanism and experiment indicated that additional reactions, possibly between silyl and methyl fragments, are needed to develop a good silane methane mechanism.

Mclain, A. G.↗