Turbulent premixed hydrogen/air flame-wall interaction with heterogeneous surface reactions

A premixed ${\mathrm{H}}_2$/air flame impinging onto a flat surface in statistically stationary state is studied for both reactive and inert wall cases to gain insights into the effects of the heterogeneous surface reactions. Direct numerical simulation (DNS) results with detailed gas-phase chemistry and surface adsorption and desorption mechanisms indicate differences in the flame front topology and near-wall flame dynamics between these two cases. In the reactive surface case, gas-phase free radicals are inclined to be adsorbed with much reduced near-wall concentration. Consequently, the gas-phase heat release rate (HRR) close to the wall decreases as well because of the low availability of free radicals. However, extra heat released from the reactive surface partially compensates for such difference. Moreover, wall reactions will intensify the turbulent fluctuations of the wall temperature and wall heat flux, while for these two cases the mean temperature profile along the flame propagating direction remains similar.     

Self-organized nanostructures in surface chemical reactions: Mechanisms and mesoscopic modeling

Nanoscale patterns can form in reactive adsorbates on catalytic surfaces as a result of attractive lateral interactions. These structures can be described within a mesoscopic theory that is derived by coarse graining the microscopic master equation thus providing a link between microscopic lattice models and reaction-diffusion equations. Such mesoscopic models allow to systematically investigate mechanisms responsible for the formation of nanoscale nonequilibrium patterns in reactive condensed matter. We have found that stationary and traveling nanostructures may result from the interplay of the attractive lateral interactions and nonequilibrium reactions. Besides reviewing these results, a detailed investigation of a single reactive adsorbate in the presence of attractive lateral interactions and global coupling through the gas phase is presented. Finally, it is outlined how a mesoscopic theory should be constructed for a particular scanning tunneling microscopy experiment [the oxidation of hydrogen on a Pt(111) surface] in order to overcome the failure of a corresponding reaction-diffusion model to quantitatively reproduce the experiments. (c) 2002 American Institute of Physics.     

Self-similar behavior and chemistry tabulation of burnt-gas diluted premixed flamelets including heat-loss

In many combustion systems, the reactive gases feeding the reaction zones are diluted by burnt products, to favor flame stabilization, homogenize the temperature distribution and reduce pollutant emission. The objective of this paper is to discuss a premixed flamelet detailed chemistry tabulation strategy for vitiated and non-adiabatic combustion. Dilution by burnt products is parameterized here with two controlling quantities: the amplitude of the heat-loss in the burnt gases, for instance at walls, and the level of reactant vitiation. The chemical response of premixed flames to variations of these parameters is studied and it is shown that most chemical properties of burnt-gas diluted flames feature self-similar behavior, which can be used to dramatically downsize chemical tables based on canonical flamelets. The self-similar behavior of the flamelets is studied for both molecular diffusion and chemical source budgets in a progress variable composition space. It is found that two different scaling relations are needed to ensure self-similar behavior of both major and radical species.     

Subgrid modeling of intrinsic instabilities in premixed flame propagation

This work is devoted to the investigation and subgrid-scale modeling of intrinsic flame instabilities occurring in the propagation of a deflagration wave. Such instabilities, of hydrodynamic and thermodiffusive origin, are expected to be of particular relevance in recent technological trends such as in the use of hydrogen as a clean energy carrier or as a secondary fuel in hydrogen enriched combustion. A dedicated set of direct numerical simulations is presented and used, in conjunction with coherent literature results, in order to develop scaling arguments for the propagation speed of self-wrinkled flames which are also supported by the outcomes of a weakly non-linear model, namely the Sivashinsky equation. The observed scaling is based on the definition of the number of unstable wavelengths in a reference hydrodynamic lengthscale, in other words the ratio between the neutral or cutoff lengthscale of intrinsic instabilities and the lateral domain of a planar flame. The scalings are then employed to develop an algebraic model for the wrinkling factor in the context of a flame surface density closure approach. An a-priori analysis shows that the model correctly captures the flame wrinkling caused by intrinsic instability at sub grid level. A strategy to include the developed self-wrinkling model in the context of a turbulent combustion model is finally discussed on the basis of the turbulence induced cut-off concept.     

Large eddy simulation of turbulent premixed combustion using tabulated detailed chemistry and presumed probability density function

A method of chemistry tabulation combined with presumed probability density function (PDF) is applied to simulate piloted premixed jet burner flames with high Karlovitz number using large eddy simulation. Thermo-chemistry states are tabulated by the combination of auto-ignition and extended auto-ignition model. To evaluate the predictive capability of the proposed tabulation method to represent the thermo-chemistry states under the condition of different fresh gases temperature, a-priori study is conducted by performing idealised transient one-dimensional premixed flame simulations. Presumed PDF is used to involve the interaction of turbulence and flame with beta PDF to model the reaction progress variable distribution. Two presumed PDF models, Dirichlet distribution and independent beta distribution, respectively, are applied for representing the interaction between two mixture fractions that are associated with three inlet streams. Comparisons of statistical results show that two presumed PDF models for the two mixture fractions are both capable of predicting temperature and major species profiles, however, they are shown to have a significant effect on the predictions for intermediate species. An analysis of the thermo-chemical state-space representation of the sub-grid scale (SGS) combustion model is performed by comparing correlations between the carbon monoxide mass fraction and temperature. The SGS combustion model based on the proposed chemistry tabulation can reasonably capture the peak value and change trend of intermediate species. Aspects regarding model extensions to adequately predict the peak location of intermediate species are discussed.     

Nonlinear dynamic modeling of surface defects in rolling element bearing systems

In this paper an analytical model is proposed to study the nonlinear dynamic behavior of rolling element bearing systems including surface defects. Various surface defects due to local imperfections on raceways and rolling elements are introduced to the proposed model. The contact force of each rolling element described according to nonlinear Hertzian contact deformation and the effect of internal radial clearance has been taken into account. Mathematical expressions were derived for inner race, outer race and rolling element local defects. To overcome the strong nonlinearity of the governing equations of motion, a modified Newmark time integration technique was used to solve the equations of motion numerically. The results were obtained in the form of time series, frequency responses and phase trajectories. The validity of the proposed model verified by comparison of frequency components of the system response with those obtained from experiments. The classical Floquet theory has been applied to the proposed model to investigate the linear stability of the defective bearing rotor systems as the parameters of the system changes. The peak-to-peak frequency response of the system for each case is obtained and the basic routes to periodic, quasi-periodic and chaotic motions for different internal radial clearances are determined. The current study provides a powerful tool for design and health monitoring of machine systems.     

Premixed turbulent combustion modeling using tabulated detailed chemistry and PDF

Tabulated chemistry and presumed probability density function (PDF) approaches are combined to perform RANS modeling of premixed turbulent combustion. The chemistry is tabulated from premixed flamelets with three independent parameters: the equivalence ratio of the mixture, the progress of reaction, and the specific enthalpy, to account for heat losses at walls. Mean quantities are estimated from presumed PDFs. This approach is used to numerically predict a turbulent premixed flame diluted by hot burnt products at an equivalence ratio that differs from the main stream of reactants. The investigated flame, subjected to high velocity fluctuations, has a thickened-wrinkled structure. A recently proposed closure for scalar dissipation rate that includes an estimation of the coupling between flame wrinkling and micromixing is retained. Comparisons of simulations with experimental measurements of mean velocity, temperature, and reactants are performed.     

A kinetic modeling study on the oxidation of primary reference fuel-toluene mixtures including cross reactions between aromatics and aliphatics

A detailed chemical kinetic model for the mixtures of primary reference fuel (PRF: n-heptane and iso-octane) and toluene has been proposed. This model is divided into three parts; a PRF mechanism [T. Ogura, Y. Sakai, A. Miyoshi, M. Koshi, P. Dagaut, Energy Fuels 21 (2007) 3233-3239], toluene sub-mechanism and cross reactions between PRF and toluene. Toluene sub-mechanism includes the low temperature kinetics relevant to engine conditions. A chemical kinetic mechanism proposed by Pitz et al. [W.J. Pitz, R. Seiser, J.W. Bozzelli, et al., in: Chemical Kinetic Characterization of the Combustion of Toluene, Proceedings of the Second Joint Meeting of the U.S. Sections of the Combustion Institute, 2001] was used as a starting model and modified by updating rate coefficients. Theoretical estimations of rate coefficients were performed for toluene and benzyl radical reactions important at low temperatures. Cross reactions between alkane, alkene, and aromatics were also included in order to account for the acceleration by the addition of toluene into iso-octane recently found in the shock tube study of the ignition delay [Y. Sakai, H. Ozawa, T. Ogura, A. Miyoshi, M. Koshi, W.J. Pitz, Effects of Toluene Addition to Primary Reference Fuel at High Temperature, SAE 2007-01-4104, 2007]. Validations of the model were performed with existing shock tube and flow tube data. The model well predicts the ignition characteristics of PRF/toluene mixtures under the wide range of temperatures (500-1700 K) and pressures (2-50 atm). It is found that reactions of benzyl radical with oxygen molecule determine the reactivity of toluene at low temperature. Although the effect of toluene addition to iso-octane is not fully resolved, the reactions of alkene with benzyl radical have the possibility to account for the kinetic interactions between PRF and toluene.     

Fluctuations and stability of stationary non-equilibrium systems in detailed balance

A generalized thermodynamic potential for Markoffian systems with detailed balance and far from thermal equilibrium has been derived in a previous paper. It was shown that the principle of detailed balance is equivalent to a set of conditions fulfilled by this potential (“potential conditions”). The properties of this potential allow us to extend the validity of a number of thermodynamic concepts well known for systems in or near thermal equilibrium to stationary states far from thermal equilibrium. The concept of symmetry breaking phase transitions for these systems is introduced in analogy to thermal equilibrium systems by considering the dependence of the stationary probability density of the system on a set of externally controlled parameters {λ}. A functional of the time dependent probability density of the system is defined in close analogy to the Gibb's definition of entropy. This functional has the properties of a Ljapunov functional of the governing Fokker-Planck equation showing the stability of the stationary probability density. The Langevin equations connected with the Fokker-Planck equation are considered. It is shown that, by means of the potential conditions, generalized “thermodynamic” fluxes and forces may be defined in such a way that the smoothly varying part of the Langevin equations (kinetic equations) constitutes a linear relation between fluxes and forces. The matrix of coefficients is given by the diffusion matrix of the Fokker-Planck equation. The symmetry relations which hold for this matrix due to the potential conditions then lead to the Onsager-Casimir symmetry relations extended to systems with detailed balance near stationary states far from thermal equilibrium. Finally it is shown that under certain additional assumptions the generalized thermodynamic potential may be used as a Ljapunov function of the kinetic equations.     

An experimental and kinetic modeling study of premixed nitroethane flames at low pressure

An experimental and kinetic modeling study is reported on three premixed nitroethane/oxygen/argon flames at low pressure (4.655 kPa) with the equivalence ratios (Φ) of 1.0, 1.5 and 2.0. Over 30 flame species were identified with tunable synchrotron vacuum ultraviolet photoionization mass spectrometry, with their mole fractions quantified as the function of the height above burner. The flame temperature profiles were measured with a Pt–6%Rh/Pt–30%Rh thermocouple. A detailed kinetic mechanism with 115 species and 730 reactions was proposed and validated against experimental results. The computed predictions have shown satisfactory agreement with the experimental results. Basing on the rate-of-production analysis, the reaction pathways that feature the combustion of nitroethane were revealed, including the primary decomposition of C–N bond fission, the oxidation of C2 and C1 hydrocarbons and the formation of nitrogenous species. The presence of NO₂ and NO has been proved to be important for these processes.     

A detailed chemical kinetic reaction mechanism for oxidation of four small alkyl esters in laminar premixed flames

A detailed chemical kinetic reaction mechanism has been developed for a group of four small alkyl ester fuels, consisting of methyl formate, methyl acetate, ethyl formate, and ethyl acetate. This mechanism is validated by comparisons between computed results and recently measured intermediate species mole fractions in fuel-rich, low-pressure, premixed laminar flames. The model development employs a principle of similarity of functional groups in constraining the H atom abstraction and unimolecular decomposition reactions for each of these fuels. As a result, the reaction mechanism and formalism for mechanism development are suitable for extension to larger oxygenated hydrocarbon fuels, together with an improved kinetic understanding of the structure and chemical kinetics of alkyl ester fuels that can be extended to biodiesel fuels. Variations in concentrations of intermediate species levels in these flames are traced to differences in the molecular structure of the fuel molecules.     

Subgrid combustion modeling of 3-D premixed flames in the thin-reaction-zone regime

Large eddy simulation (LES) is used to investigate three-dimensional (3D) lean premixed turbulent methane–air flames in the thin-reaction-zone regime. In this regime, the Kolmogorov scale is smaller than the preheat zone thickness, but larger than the reaction zone thickness. Past numerical studies of similar flames were primarily direct numerical simulation either in two-dimensions or using the artificially thickened flame approach in 3D. For an LES the effect of small (unresolved) scales on the scalar field must be, modeled accurately to capture the correct flame structure. A subgrid combustion model based on the linear-eddy-mixing (LEM) model is used within an LES framework (called LEM–LES hereafter) to capture the 3D flame-structure of the highly stretched premixed flames. A finite-rate, one-step methane–air chemistry with a non-unity Lewis number formulation is used in this study. The simulated flame structure resembles flames experimentally studied in the thin-reaction-zone regime. Even though the preheat zone is broadened by the penetration of small eddies, the chemical reaction zone remains thin and localized. This feature is captured properly in the current LEM–LES approach. The flame structure and other statistics such as the flame area evolution, curvature, and strain-rate statistics computed using the LEM–LES are also in good agreement with the past DNS studies.