The effect of bulk gas diffusivity on apparent pulverized coal char combustion kinetics
Not Available
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.
Not Available
Explore the source record for details and available documents.
In oxy-fuel combustion, the pure oxygen (O 2 ), diluted with CO 2 is used as oxidant instead air. Hence, the combustion products (CO 2 and H 2 O) are free from pollution by nitrogen oxides. Moreover, high pressures results in the near-liquid density of CO 2 at supercritical state (sCO 2 ). Unfortunately, the effects of sCO 2 on the combustion kinetics are far from being understood. In order to assist in this understanding, in this work we are using quantum chemistry methods. Here we investigate potential energy surfaces of important combustion reactions in the presence of carbon dioxide melocule. All transition states, reactant and product complexes are reported for three reactions: H 2 CO+HO 2 →HCO+H 2 O 2 (R1), 2HO 2 →H 2 O 2 +O 2 (R2), and CO+OH→CO 2 +H (R3). In the reaction R3, covalent binding of CO 2 to OH radical and then CO molecule opens a new pathway, including hydrogen transfer from oxygen to carbon atoms followed by CH bond dissociation. Compared to bimolecular OH+CO mechanism, this pathway reduces the activation barrier by 5 kcal/mol, and is expected to accelerate the reaction. This is the first report of autocatalytic effect in combustion. In case of hydroperoxyl self-reaction 2HO 2 →H 2 O 2 +O 2 the intermediates, containing covalent bonds to CO 2 were found not to be competitive. However, the spectator CO 2 molecule is able to stabilize the cyclic transition state and lower the barrier by 3 kcal/mol. Formation of covalent intermediates was also discovered in H 2 CO+HO 2 →HCO+H 2 O 2 reaction, but these specie lead to substantially higher activation barriers which makes them unlikely to play role in hydrogen transfer kinetics. The van der Waals complexation with carbon dioxide also stabilized transition state and reduces reaction barrier. Lastly, these results indicate that CO 2 environment is likely to have catalytic effect on combustion reactions, which needs to be included in kinetic combustion mechanisms in supercritical CO 2 .
To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.
An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.
Idealized reactors, that are designed to ensure perfect mixing and are used to generate the combustion kinetics for complex hydrocarbon fuels, may depart from the ideal and influence the kinetics model performance. A complex hydrocarbon kinetics model that was established by modeling a jet-stirred combustor (JSC) as a perfectly stirred reactor (PSR), is reevaluated with a simple stochastic process in order to introduce the unmixedness effect quantitatively into the reactor system. It is shown that the comparisons of the predictions and experimental data have improved dramatically with the inclusion of the unmixedness effect in the rich combustion region. The complex hydrocarbon kinetics is therefore verified to be mixing effect free and be applicable to general reacting flow calculations.
A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.
This presentation introduces our research on using droplet burning data to validate combustion kinetic mechanisms of transportation fuels: challenges and opportunities.
The stiff ordinary differential equation solution methods more commonly used for chemical kinetics problems are examined, and the efficiency and accuracy of several of the methods are compared by applying them to two combustion kinetics problems. The local sensitivity analysis methods developed to date are outlined, and their advantages and difficulties are briefly discussed. Examples are used to make comparisons among the different techniques.
Explore the source record for details and available documents.
Acetaldehyde is an important intermediate and a toxic emission in the combustion of fuels, especially for biofuels. To better understand its combustion characteristics, a detailed chemical kinetic model describing the oxidation of acetaldehyde has been developed and comprehensively validated against various types of literature data including laminar flame speeds, oxidation and pyrolysis in shock tubes, chemical structure of premixed flames, and low-temperature oxidation in jet-stirred reactors. To extend the validation range, the chemical structure of a counterflow flame fueled by acetaldehyde at 600 Torr has been measured using vacuum ultra-violet photoionization molecular-beam mass spectrometry. In addition, ignition delay times at 10 atm and 700-1100 K were measured in a rapid compression machine, and a negative temperature coefficient (NTC) behavior was observed. Here, the present kinetic model well reproduces the results of various acetaldehyde combustion experiments covering wide ranges of temperatures (300–2300 K) and pressures (0.02–10 atm), and explains well the observed NTC behavior based on the competition between multiple oxidation pathways for the methyl radicals and their self-recombination forming ethane, a relatively stable species at temperatures below 1000 K.
Not Available
Supercritical carbon dioxide (sCO 2 ) cycles are being investigated for the future of power generation and will contribute to a carbon-neutral future to combat the effects of climate change. These direct-fired closed cycles will produce power without adding significant pollutants to the atmosphere. For such cycles to be efficient, they will need to operate at significantly higher pressures (e.g., 300 atm for Allam Cycle) than existing systems (typically less than 40 atm). There is limited knowledge on combustion at these high pressures and with a high dilution of carbon dioxide. Also, various experimental and computational investigations and model developments have been performed by the University of Central Florida (UCF), Embry-Riddle Aeronautical University (ERAU), and Stanford University (S.U.) to improve the current knowledge base and support the design and development of sCO 2 combustors. This project's technical aspects include chemical kinetics development, fundamental sCO 2 combustion, combustion model and sub-model development, supercritical fluid injection characterization, and heat transfer characterization.
The Ignition Quality Tester (IQT) is a constant volume spray combustion system that can be heated and pressurized to conditions that are similar to a diesel engine at top dead center. With no moving parts and the ability to handle low volatility fuels, the IQT can be a bridge between engines and traditional methods for studying chemical kinetics. By comparing experimental data with model predictions, the IQT has been used to validate skeletal kinetic models of ignition. CFD modeling of the IQT using KIVA-3V was used to predict ignition of n-heptane accurately. Operating the IQT in a regime where chemical kinetics dominates (long ignition delays) allowed NTC behavior to be observed for some isomers of heptane. Experimental results with the low volatility fuel heptamethylnonane also show NTC behavior. At long ignition delays, experimental results can be compared with 0-D detailed chemical mechanisms.
arge-scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. As a first step toward this goal, we developed a procedure for automatically generating the thermochemical data for combustion of an arbitrary fuel. The procedure begins by producing a list of combustion relevant species from a specification of the fuel and combustion conditions of interest. Then, for each element in the list of species, the procedure determines an internal coordinate z-matrix description of its structure, the optimal torsional configuration via Monte Carlo sampling, key rovibrational properties for that optimal geometry (including anharmonic corrections from torsional mappings and/or vibrational perturbation theory), and high level estimates of the electronic and zero-point energies via arbitrarily defined composite methods. This dataset is then converted first to partition functions, then to thermodynamic properties, and finally to NASA polynomial representations of the data. The end product is an automatically generated database of electronic structure results and thermochemical data including representations in a format appropriate for combustion simulations. The utility and functioning of this predictive automated computational thermochemistry (PACT) software package is illustrated through application to the automated generation of thermochemical data for the combustion of n-butane. Butane is chosen for this demonstration as its species list is of reasonably manageable size for debugging level computations, while still presenting most of the key challenges that need to be surmounted in the consideration of larger fuels. Furthermore, its low temperature chemistry is representative of that oc-curring with larger alkanes.
One of the major biofuels used today is biodiesel, composed of fatty acid methyl esters. The fats and oils feedstocks used to make biodiesel are in relatively short supply. To meet increasing global demand, engineered microorganisms have been developed to catalyze conversion of sugar to fatty acid esters that contain unique ..beta..-hydroxy-esters. This study investigated the effect of a hydroxyl group on the combustion characteristics of methyl esters by comparing the chemical behavior of methyl hexanoate (MHx) and methyl 3-hydroxyhexanoate (M3OHHx) - used as surrogates for diesel-boiling-range esters with longer fatty acid chains. The oxidation of these esters was studied experimentally in a flow reactor at 0.84 and 10 bar, 600 to 1,100 K, and stoichiometric conditions; and in a constant volume combustion chamber at 5 and 10 bar, 600 to 900 K, and equivalence ratios of 0.3 and 0.6. MHx was more reactive in the constant volume chamber at temperatures below 800 K at 10 bar and equivalence ratio of 0.6. MHx also exhibited a higher indicated cetane number (16.4) than M3OHHx (8.1). We investigated this reactivity trend using an updated MHx kinetic model and M3OHHx model developed as part of this work. The kinetic models predicted that radicals formed from MHx were consumed by low-temperature chemistry reactions, whereas the M3OHHx reactivity was significantly governed by the ..beta..-radical (on the same carbon as the OH group) chemistry which mainly terminated via the less reactive chain propagation pathway to methyl-3-oxohexanoate + HO2. This study extends our understanding of structural effects on ester reactivity, which will allow for accurate surrogate formulation and performance simulations.
A perturbation analysis modeling approach based on a genetic algorithm was used to identify possible reaction pathways to explain previous experimental observations of the strong acceleration of syngas auto-ignition by trimethylsilanol (TMSO). Organosilicon reactions were taken from an existing chemical kinetic mechanism for tetramethylsilane pyrolysis which also contained TMSO reactions. The experimental targets for optimization were ignition delay times and peak OH radical concentrations in the range 1010–1070 K, 5 atm with dilute (4%) ϕ = 0.1 syngas mixtures in N 2 /Ar containing 0, 200, and 1000 ppm TMSO. Hundreds of thousands of trials resulted in three models which accurately predicted the experimental auto-ignition times and peak OH concentrations. Rate of production analyses indicated that TMSO forms dimethylsilanediol, which then reacts creating two “catalytic loops.” The net effect of the loops is to accelerate the reaction H 2 + O = OH + H at a faster effective rate than H 2 /O 2 kinetics alone, thus enhancing the reactivity of syngas. The formation of the loops may be attributable to the high Si–C bond energies in siloxanols, reducing rates of thermal decomposition and allowing dimethylsilanediol to react with H 2 for a longer duration compared with its alkane analog. Here, the results of the current work provide valuable direction for fundamental studies of these new hypothesized reaction pathways.