Modeling of Low-temperature oxidation and combustion of droplets

Key features of radiation extinguishing of spherical hot flame around a single droplet with its subsequent low-temperature oxidation and combustion under microgravity conditions—a phenomenon discovered in experiments onboard the International Space Station—have been reproduced using the mathematical model of droplet combustion and detailed kinetic mechanism of n-heptane oxidation and combustion. It has been demonstrated that experimentally observed repeated temperature flashes were blue flame flashes, and their duration was determined by the hydrogen peroxide decomposition time. In addition to this phenomenon, calculations predict the existence of new modes of low-temperature oxidation and combustion of droplets without the hot flame stage. In such modes, the basic reaction is concentrated very close to the droplet surface, and fuel vapor reacts in it only partially.     

A comprehensive kinetic mechanism for CO,CH2O,and CH3OH combustion

New experimental profiles of stable species concentrations are reported for formaldehyde oxidation in a variable pressure flow reactor at initial temperatures of 850–950 K and at constant pressures ranging from 1.5 to 6.0 atm. These data, along with other data published in the literature and a previous comprehensive chemical kinetic model for methanol oxidation, are used to hierarchically develop an updated mechanism for CO/H₂O/H₂/O₂, CH₂O, and CH₃OH oxidation. Important modifications include recent revisions for the hydrogen–oxygen submechanism (Li et al., Int J Chem Kinet 2004, 36, 565), an updated submechanism for methanol reactions, and kinetic and thermochemical parameter modifications based upon recently published information. New rate constant correlations are recommended for CO + OH = CO₂ + H ( R23 ) and HCO + M = H + CO + M ( R24 ), motivated by a new identification of the temperatures over which these rate constants most affect laminar flame speed predictions (Zhao et al., Int J Chem Kinet 2005, 37, 282). The new weighted least‐squares fit of literature experimental data for ( R23 ) yields k₂₃ = 2.23 × 10⁵T^1.89exp(583/T) cm³/mol/s and reflects significantly lower rate constant values at low and intermediate temperatures in comparison to another recently recommended correlation and theoretical predictions. The weighted least‐squares fit of literature results for ( R24 ) yields k₂₄ = 4.75 × 10¹¹T^0.66exp(?7485/T) cm³/mol/s, which predicts values within uncertainties of both prior and new (Friedrichs et al., Phys Chem Chem Phys 2002, 4, 5778; DeSain et al., Chem Phys Lett 2001, 347, 79) measurements. Use of either of the data correlations reported in Friedrichs et al. (2002) and DeSain et al. (2001) for this reaction significantly degrades laminar flame speed predictions for oxygenated fuels as well as for other hydrocarbons. The present C₁/O₂ mechanism compares favorably against a wide range of experimental conditions for laminar premixed flame speed, shock tube ignition delay, and flow reactor species time history data at each level of hierarchical development. Very good agreement of the model predictions with all of the experimental measurements is demonstrated. © 2007 Wiley Periodicals, Inc. 39: 109–136, 2007     

Plasma Assisted Low Temperature Combustion

This paper presents recent kinetic and flame studies in plasma assisted low temperature combustion. First, the kinetic pathways of plasma chemistry to enhance low temperature fuel oxidation are discussed. The impacts of plasma chemistry on fuel oxidation pathways at low temperature conditions, substantially enhancing ignition and flame stabilization, are analyzed base on the ignition and extinction S-curve. Secondly, plasma assisted low temperature ignition, direct ignition to flame transition, diffusion cool flames, and premixed cool flames are demonstrated experimentally by using dimethyl ether and n-heptane as fuels. The results show that non-equilibrium plasma is an effective way to accelerate low temperature ignition and fuel oxidation, thus enabling the establishment of stable cool flames at atmospheric pressure. Finally, the experiments from both a non-equilibrium plasma reactor and a photolysis reactor are discussed, in which the direct measurements of intermediate species during the low temperature oxidations of methane/methanol and ethylene are performed, allowing the investigation of modified kinetic pathways by plasma-combustion chemistry interactions. Finally, the validity of kinetic mechanisms for plasma assisted low temperature combustion is investigated. Technical challenges for future research in plasma assisted low temperature combustion are then summarized.     

Validation of the extended simultaneous kinetics and ringdown model by measurements of the reaction NH2 + NO

The determination of rate constants for fast chemical reactions from nonexponential cavity ringdown profiles requires a consideration of the interfering laser bandwidth effect that arises if the line width of the ringdown probe laser exceeds the absorption line width of the detected species. The deconvolution of the kinetics and the bandwidth effect can be accomplished with the extended simultaneous kinetics and ringdown (eSKaR) model presented by Guo et al. (Guo, et al. Phys. Chem. Chem. Phys. 2003, 5, 4622). We present a detailed validation of this eSKaR model by a corresponding investigation of the well-known rate constant for the reaction NH2 + NO. Line profiles of the pulsed ringdown probe laser and the NH2 absorption line were determined from forward modeling of experimental ringdown profiles and verified by narrow-bandwidth laser absorption measurements. In addition, the rate constant for the title reaction was evaluated using the eSKaR model and also by means of a conventional pump-probe approach with variable time delays between the photolysis (pump) and ringdown (probe) laser pulses. The resulting room temperature rate constant for the NH2 + NO reaction, k1= (8.5 +/- 1.0) x 10(12) cm(3) mol(-1) s(-1), and the room temperature pressure broadening coefficient of NH2, = 2.27 GHz/bar, measured on the A2A1<-- X2B1 transition at wavelengths around lambda = 597 nm, were found to be in excellent agreement with the available literature data.     

Recombination of ozone via the chaperon mechanism

The recombination of ozone via the chaperon mechanism, i.e., ArO+O2 --> Ar+O3 and ArO2+O --> Ar+O3, is studied by means of classical trajectories and a pairwise additive Ar-O3 potential energy surface. The recombination rate coefficient has a strong temperature dependence, which approximately can be described by T(-n) with n approximately 3. It is negligible for temperatures above 700 K or so, but it becomes important for low temperatures. The calculations unambiguously affirm the conclusions of Hippler et al. [J. Chem. Phys. 93, 6560 (1990)] and Luther et al. [Phys. Chem. Chem. Phys. 7, 2764 (2005)] that the chaperon mechanism makes a sizable contribution to the recombination of O3 at room temperature and below. The dependence of the chaperon recombination rate coefficient on the isotopomer, studied for two different isotope combinations, is only in rough qualitative agreement with the experimental data. The oxygen atom isotope exchange reaction involving ArO and ArO2 van der Waals complexes is also investigated; the weak binding of O or O2 to Ar has only a small effect.     

Reaction mechanism of oxidation, hydroxylation, and epoxidation by hypofluorous acid: a theoretical study of unusual H-bond-assisted catalysis

The oxidation of organic molecules by hypofluorous acid (HOF) was studied extensively and systematically by Rozen et al. Therefore, it seems appropriate to refer to the process as Rozen oxidation. An entire set of model molecules was selected for quantum chemical investigation of the oxidation mechanism: a C=C double bond in ethylene, sulfur and selenium in dimethyl derivatives, nitrogen and phosphorus in trimethyl derivatives, as well as methyl azides. In the gas phase, van der Waals complexes between HOF and the previously mentioned species easily are formed, but these complexes are reluctant to undergo oxidation. The addition of another HOF molecule connected with the formation of a cyclic complex (i.e., substrate and two molecules of HOF) seems to be decisive for the oxidation process. The attempt to substitute the second HOF molecule with H2O demonstrated the superiority of HOF. Complexes of this kind decompose along the reaction path smoothly (i.e., with a low activation energy) to the respective oxidation product. A potential role of the hydroxyl cation (HO+) in the oxidation step is mentioned. Besides an oxidation product, one HOF molecule is released (an essential feature of catalysis), and furthermore, hydrogen fluoride is formed. It was suggested by Sertchook et al. (J. Phys. Chem. A 2006, 110, 8275) that the interaction between the substrate to be oxidized and HOF is catalytically influenced by the HF molecule. The mechanism suggested here is more feasible and, particularly at the early stages of the oxidation process, decisive. Also, the role of acetonitrile, used as a solvent by Rozen et al., is discussed in terms of a continuum model. Moreover, passing from potential energies to Gibbs energies is considered.     

Iodine oxidation by hydrogen peroxide and Bray-Liebhafsky oscillating reaction: effect of the temperature

This work presents a new experimental kinetic study at 39° and 50° of the iodine oxidation by hydrogen peroxide. The results allow us to obtain the temperature effect on the rate constants previously proposed at 25° for our model of the Bray-Liebhafsky oscillating reaction (G. Schmitz, Phys. Chem. Chem. Phys. 2010, 12, 6605.). The values calculated with the model are in good agreement with many experimental results obtained under very different experimental conditions. Numerical simulations of the oscillations observed formerly by different authors are presented, including the evolutions of the iodine, hydrogen peroxide, iodide ions and oxygen concentrations. Special attention is paid to the perturbing effects of oxygen and of the iodine loss to the gas phase.     

Application of scaling and kinetic equations to helium cluster size distributions: Homogeneous nucleation of a nearly ideal gas

A previously published model of homogeneous nucleation [Villarica et al., J. Chem. Phys. 98, 4610 (1993)] based on the Smoluchowski [Phys. Z. 17, 557 (1916)] equations is used to simulate the experimentally measured size distributions of 4He clusters produced in free jet expansions. The model includes only binary collisions and does not consider evaporative effects, so that binary reactive collisions are rate limiting for formation of all cluster sizes despite the need for stabilization of nascent clusters. The model represents these data very well, accounting in some cases for nearly four orders of magnitude in variation in abundance over cluster sizes ranging up to nearly 100 atoms. The success of the model may be due to particularities of 4He clusters, i.e., their very low coalescence exothermicity, and to the low temperature of 6.7 K at which the data were collected.     

异辛烷/正庚烷/乙醇三组分燃料着火的化学动力学模型

提出一个包括异辛烷、正庚烷和乙醇的三组分燃料的着火动力学模型, 该机理包括50 个组分和193 个反应. 通过路径分析和灵敏度分析, 给出了基础燃料在高低温条件下的不同反应路径和影响氧化过程的重要基元反应. 该机理预测的单组分(异辛烷、正庚烷、乙醇)燃料、双组分基础燃料和三组分燃料的点火延迟时间与实验值有很高一致性. 本文机理包含较少的组分数与反应数, 因而可适用汽油掺烧乙醇的多维计算流体动力学(CFD)数值模拟.     

Chemical kinetic study of the effect of a biofuel additive on jet-A1 combustion

The kinetics of oxidation of kerosene Jet A-1 and a kerosene/rapeseed oil methyl ester (RME) mixture (80/20, mol/mol) (biokerosene) was studied experimentally in a jet-stirred reactor at 10 atm and constant residence time, over the temperature range 740-1200 K, and for variable equivalence ratios (0.5-1.5). Concentration profiles of the reactants, stable intermediates, and final products were obtained by probe sampling followed by on-line and off-line gas chromatography analyses. The oxidation of these fuels in these conditions was modeled using a detailed kinetic reaction mechanism consisting of 2027 reversible reactions and 263 species. The surrogate biokerosene model fuel used here consisted of a mixture of n-hexadecane, n-propylcyclohexane, n-propylbenzene, and n-decane, where the long-chain methyl ester fraction was simply represented by n-hexadecane. The proposed kinetic reaction mechanism used in the modeling yielded a good representation of the kinetics of oxidation of kerosene and biokerosene under jet-stirred reactor conditions and of kerosene in a premixed flame. The data and the model showed the biokerosene (Jet A-1/RME mixture) has a slightly higher reactivity than Jet A-1, whereas no major modification of the product distribution was observed besides the formation of small unsaturated methyl esters produced from RME's oxidation. The model predicts no difference in the ignition delays of kerosene and biokerosene. Using the proposed kinetic scheme, the formation of potential soot precursors was studied with particular attention.     

Experimental and modeling study of the oxidation of benzene

This paper describes an experimental and modeling study of the oxidation of benzene. The low‐temperature oxidation was studied in a continuous flow stirred tank reactor with carbon‐containing products analyzed by gas chromatography. The following experimental conditions were used: 923 K, 1 atm, fuel equivalence ratios from 1.9 to 3.6, concentrations of benzene from 4 to 4.5%, and residence times ranging from 1 to 10 s corresponding to benzene conversion yields from 6 to 45%. The ignition delays of benzene–oxygen–argon mixtures with fuel equivalence ratios from 1 to 3 were measured behind shock waves. Reflected shock waves permitted to obtain temperatures from 1230 to 1970 K and pressures from 6.5 to 9.5 atm. A detailed mechanism has been proposed and allows us to reproduce satisfactorily our experimental results, as well as some data of the literature obtained in other conditions, such as in a plug flow reactor or in a laminar premixed flame. The main reaction paths have been determined for the four series of measurements by sensitivity and flux analyses. © 2003 Wiley Periodicals, Inc. Int J Chem Kinet 35: 503–524, 2003     

Structures and correlation functions of multicomponent and polydisperse hard-sphere mixtures from a density functional theory

The structures of nonuniform binary hard-sphere mixtures and the correlation functions of uniform ternary hard-sphere mixtures were studied using a modified fundamental-measure theory based on the weight functions of Rosenfeld [Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)] and Boublik-Mansoori-Carnahan-Starling-Leland equation of state [Boublik, J. Chem. Phys. 53, 471 (1970); Mansoori et al., J. Chem. Phys. 54, 1523 (1971)]. The theoretical predictions agreed very well with the molecular simulations for the overall density profiles, the local compositions, and the radial distribution functions of uniform as well as inhomogeneous hard-sphere mixtures. The density functional theory was further extended to represent the structure of a polydisperse hard-sphere fluid near a hard wall. Excellent agreement was also achieved between theory and Monte Carlo simulations. The density functional theory predicted oscillatory size segregations near a hard wall for a polydisperse hard-sphere fluid of a uniform size distribution.     

Isotope branching and tunneling in O(3P)+HD-->OH+D; OD+H reactions

The O((3)P)+HD and O((3)P)+D(2) reactions are studied using quantum scattering calculations and chemically accurate potential energy surfaces developed for the O((3)P)+H(2) system by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)]. Cross sections and rate coefficients for OH and OD products are calculated using accurate quantum methods as well as the J-shifting approximation. The J-shifting approach is found to work remarkably well for both O+HD and O+D(2) collisions. The reactions are dominated by tunneling at low temperatures and for the O+HD reaction the hydrogen atom transfer leading to the OH product dominates at low temperatures. Our result for the OH/OD branching ratio is in close agreement with previous calculations over a wide range of temperatures. The computed OH/OD branching ratios are also in close agreement with experimental results of Robie et al. [Chem. Phys. Lett. 134, 579 (1987)] at temperatures above 400 K but the theoretical results do not reproduce the rapid rise in the experimental values of the branching ratio for temperatures lower than 350 K. We believe that new measurements could resolve the long-standing discrepancy between experiment and theory for this benchmark reaction.     

On the possible catalytic role of a single water molecule in the acetone + OH gas phase reaction: a theoretical pseudo-second-order kinetics study

In this work, we have revisited the mechanism of the acetone + OH radical reaction assisted by a single water molecule simulating atmospheric conditions. Density functional methods are employed in conjunction with CCSD(T) and large basis sets to explore the potential energy surface of this radical-molecule reaction. Computational kinetics calculations in a pseudo-second order mechanism have been performed, taking into account average atmospheric water concentrations and temperatures. We have used this method recently to study the single-water molecule-assisted H-abstraction in acetaldehyde (Iuga et al. in J Phys Chem Lett 1:3112, 2010) and in glyoxal (Iuga et al. in Chem Phys Lett 501:11, 2010) by OH radicals, and we showed that the initial water complexation step is essential in the rate constant calculation. In both cases, the amount of complex formed is only about 0.01% of the total organic molecule concentration, and as a consequence, water does not accelerate the reaction. In the acetone reaction with OH radicals under atmospheric conditions, we also find that the water–acetone complex concentration is much too small to be relevant, and thus, the rate constant of the water-assisted mechanism is orders of magnitude smaller than the water-free corresponding value.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Origin of bulklike structure and bond length disorder of Pt37 and Pt6Ru31 clusters on carbon: comparison of theory and experiment

We describe a theoretical analysis of the structures of self-organizing nanoparticles formed by Pt and Ru-Pt on carbon support. The calculations provide insights into the nature of these metal particle systems-ones of current interest for use as the electrocatalytic materials of direct oxidation fuel cells-and clarify complex behaviors noted in earlier experimental studies. With clusters deposited via metallo-organic Pt or PtRu(5) complexes, previous experiments [Nashner et al. J. Am. Chem. Soc. 1997, 119, 7760; Nashner et al. J. Am. Chem. Soc. 1998, 120, 8093; Frenkel et al. J. Phys. Chem. B 2001, 105, 12689] showed that the Pt and Pt-Ru based clusters are formed with fcc(111)-stacked cuboctahedral geometry and essentially bulklike metal-metal bond lengths, even for the smallest (few atom) nanoparticles for which the average coordination number is much smaller than that in the bulk, and that Pt in bimetallic [PtRu(5)] clusters segregates to the ambient surface of the supported nanoparticles. We explain these observations and characterize the cluster structures and bond length distributions using density functional theory calculations with graphite as a model for the support. The present study reveals the origin of the observed metal-metal bond length disorder, distinctively different for each system, and demonstrates the profound consequences that result from the cluster/carbon-support interactions and their key role in the structure and electronic properties of supported metallic nanoparticles.     

Electrochemical Oxidation of Carbon‐Containing Fuels and Their Dynamics in Low‐Temperature Fuel Cells

Fuel cells can convert the energy that is chemically stored in a compound into electrical energy with high efficiency. Hydrogen could be the first choice for chemical energy storage, but its utilization is limited due to storage and transport difficulties. Carbon‐containing fuels store chemical energy with significantly higher energy density, which makes them excellent energy carriers. The electro‐oxidation of carbon‐containing fuels without prior reforming is a more challenging and complex process than anodic hydrogen oxidation. The current understanding of the direct electro‐oxidation of carbon‐containing fuels in low‐temperature fuel cells is reviewed. Furthermore, this review covers various aspects of electro‐oxidation for carbon‐containing fuels in non‐steady‐state reaction conditions. Such dynamic investigations open possibilities to elucidate detailed reaction kinetics, to sense fuel concentration, or to diagnose the fuel‐cell state during operation. Motivated by the challenge to decrease the consumption of fossil fuel, the production routes of the fuels from renewable resources also are reviewed.     

High-temperature oxidation chemistry of n-butanol--experiments in low-pressure premixed flames and detailed kinetic modeling

An automated reaction mechanism generator is used to develop a predictive, comprehensive reaction mechanism for the high-temperature oxidation chemistry of n-butanol. This new kinetic model is an advancement of an earlier model, which had been extensively tested against earlier experimental data (Harper et al., Combust. Flame, 2011, 158, 16-41). In this study, the model's predictive capabilities are improved by targeting isomer-resolved quantitative mole fraction profiles of flame species in low-pressure flames. To this end, a total of three burner-stabilized premixed flames are isomer-selectively analyzed by flame-sampling molecular-beam time-of-flight mass spectrometry using photoionization by tunable vacuum-ultraviolet synchrotron radiation. For most species, the newly developed chemical kinetic model is capable of accurately reproducing the experimental trends in these flames. The results clearly indicate that n-butanol is mainly consumed by H-atom abstraction with H, O, and OH, forming predominantly the α-C(4)H(9)O radical (CH(3)CH(2)CH(2)˙CHOH). Fission of C-C bonds in n-butanol is only predicted to be significant in a similar, but hotter flame studied by O?wald et al. (Combust. Flame, 2011, 158, 2-15). The water-elimination reaction to 1-butene is found to be of no importance under the premixed conditions studied here. The initially formed isomeric C(4)H(9)O radicals are predicted to further oxidize by reacting with H and O(2) or to decompose to smaller fragments via β-scission. Enols are detected experimentally, with their importance being overpredicted by the model.     

Examining the robustness of first-principles calculations for metal hydride reaction thermodynamics by detection of metastable reaction pathways

First principles calculations have played a useful role in screening mixtures of complex metal hydrides to find systems suitable for H(2) storage applications. Standard methods for this task efficiently identify the lowest energy reaction mechanisms among all possible reactions involving collections of materials for which DFT calculations have been performed. The resulting mechanism can potentially differ from physical reality due to inaccuracies in the DFT functionals used, or due to other approximations made in estimating reaction free energies. We introduce an efficient method to probe the robustness of DFT-based predictions that relies on identifying reactions that are metastable relative to the lowest energy reaction path predicted with DFT. An important conclusion of our calculations is that in many examples DFT cannot unambiguously predict a single reaction mechanism for a well defined metal hydride mixture because two or more mechanisms have reaction energies that differ by a small amount. Our approach is illustrated by analyzing a series of single step reactions identified in our recent work that examined reactions with a large database of solids [Kim et al., Phys. Chem. Chem. Phys. 2011, 13, 7218].