Turbulent Premixed Flame Modeling Using the Algebraic Flame Surface Wrinkling Model: A Comparative Study between OpenFOAM and Ansys Fluent

Hereafter, we used the Algebraic Flame Surface Wrinkling (AFSW) model to conduct numerical simulations of the Paul Scherrer Institute (PSI) high-pressure, turbulent premixed Bunsen flame experiments. We implemented the AFSW model in OpenFOAM and in Ansys Fluent, and we compared the outcome of both solvers against the experimental results. We also highlight the differences between both solvers. All the simulations were performed using a two-dimensional axisymmetric model with the standard k−ϵ turbulence model with wall functions. Two different fuel/air mixtures were studied, namely, a 100%CH4 volumetric ratio and a 60%CH4+ 40%H2 volumetric ratio. The thermophysical and transport properties of the mixture were calculated as a function of temperature using the library Cantera (open-source suite of tools for problems involving chemical kinetics, thermodynamics, and transport processes), together with the GRI-Mech 3.0 chemical mechanism. It was found that the outcome of the AFSW model implemented in both solvers was in good agreement with the experimental results, quantitatively and qualitatively speaking. Further assessment of the results showed that, as much as the chemistry, the turbulence model and turbulent boundary/initial conditions significantly impact the flame shape and height.

Effect of Microwave Pulses on the Morphology and Development of Spark-Ignited Flame Kernel

Microwave-assisted spark ignition (MAI) is a promising way to enhance the ignition performance of engines under lean conditions. To understand the effect of microwave-induced flow during MAI, the development and morphology of spark-ignited methane-air flame kernel under various microwave pulse parameters are experimentally studied. Experiments are conducted in a constant volume combustion chamber, and flame development is recorded through a high-speed shadowgraph method. Flame area and deformation index are adopted to evaluate the flame characteristic. Results show that increasing the microwave pulse energy from 0 to 150 mJ exhibits a threshold process for expanding the flame kernel area under 0.2 MPa ambient pressure. When the pulse energy is below the threshold of 90 mJ, the microwave enhancing efficiency is much lower than that beyond the threshold. Increasing microwave pulse repetition frequency (PRF) changes the flow on flame surface and raises the absorption efficiency for microwave energy, and thus helps to improve the MAI performance under higher pressures. Hence, 1 kHz pulses cause more obvious flame deformation than those with higher PRF pulses under 0.2 MPa, while this tendency is reversed as the ambient pressure increases to 0.6 MPa. Besides, microwave pulses of different repetition frequencies lead to different flame kernel morphology, implying the various regimes behind the interaction between a microwave and spark kernel.

A Semi-Analytical Model for Prediction of Wall Quenching Distances of Premixed Flames

In order to predict the variation of the wall quenching distance of a premixed flame under different equivalence ratios and incoming flow velocities, a semi-analytical model applied to both lean single-layer flame and rich double-layer flame has been derived based on the conservation of energy in the quenching zone. In this model the flame surface radiation plays an important role. Factors influencing the radiation have been analyzed, respectively. The model indicates that the factors affecting the quenching distance in premixed flame are more complicated than that in single-wall flames or flames in tube. To fit the empirical coefficient in this model, a methane-air premixed flame quenching distance experiment under both lean and rich conditions has been performed. The comparison between the theoretical prediction and the experiment result shows that this semi-analytical model gives a suitable description of the quenching distance. The relative error of the quenching distance under different equivalence ratios and incoming flow velocities is less than ±15%.

Effects of Body Forces on the Statistics of Flame Surface Density and Its Evolution in Statistically Planar Turbulent Premixed Flames

Body forces such as buoyancy and externally imposed pressure gradients are expected to have a strong influence on turbulent premixed combustion due to the considerable changes in density between the unburned and fully burned gases. The present work utilises Direct Numerical Simulation data of three-dimensional statistically planar turbulent premixed flames to study the influence of body forces on the statistical behaviour of the flame surface density (FSD) and its evolution within the flame brush. The analysis has been carried out for different turbulence intensities and normalised body force values (i.e., Froude numbers). A positive value of the body force signifies an unstable density stratification (i.e., body force is directed from the heavier unburned gas to the lighter burned gas) and vice versa. It is found that for a given set of turbulence parameters, flame wrinkling increases with an increase in body force magnitude in the unstable configuration. Furthermore, higher magnitudes of body force in the unstable density stratification configuration promote a gradient type transport of turbulent scalar and FSD fluxes, and this tendency weakens in the stable density stratification configuration where a counter-gradient type transport is promoted. The statistical behaviours of the different terms in the FSD transport equation and their closures in the context of Reynolds Averaged Navier–Stokes simulations have been analysed in detail. It has been demonstrated that the effects of body force on the FSD and the terms of its transport equation weakens with increasing turbulence intensity as a result of the diminishing relative strength of body force in comparison to the inertial force. The predictions of the existing models have been assessed with respect to the corresponding terms extracted from the explicitly averaged DNS data, and based on this evaluation, suitable modifications have been made to the existing models to incorporate the effects of body force (or Froude number).

Modelling of Conditional Scalar Dissipation Rate in Turbulent Premixed Combustion

In turbulent premixed flames, for the mixing at a molecular level of reactants and products on the flame surface, it is crucial to sustain the combustion. This mixing phenomenon is featured by the scalar dissipation rate, which may be broadly defined as the rate of micro-mixing at small scales. This term, which appears in many turbulent combustion methods, includes the Conditional Moment Closure (CMC) and the Probability Density Function (PDF), requires an accurate model. In this study, a mathematical closure for the conditional mean scalar dissipation rate, <Nc|ζ>, in Reynolds, Averaged Navier–Stokes (RANS) context is proposed and tested against two different Direct Numerical Simulation (DNS) databases having different thermochemical and turbulence conditions. These databases consist of lean turbulent premixed V-flames of the CH4-air mixture and stoichiometric turbulent premixed flames of H2-air. The mathematical model has successfully predicted the peak and the typical profile of <Nc|ζ> with the sample space ζ and its prediction was consistent with an earlier study.

The Effect of Pressure on NOx Entitlement and Reaction Timescales in a Premixed Axial Jet-In-Crossflow

This paper investigates the pressure dependency of a lean premixed jet injected into a lean vitiated crossflow with an experimentally verified detailed chemistry computational fluid dynamics (CFD) model and 53 species considered. Experimental data were taken in an axially staged combustor with an optically accessible test section, allowing the use of particle image velocimetry (PIV) and CH* chemiluminescence techniques as well as point measurement of species concentration, temperature, and pressure. The experimental data cases at one, three, and five atmospheres were selected to describe the flame stabilization dependency on pressure and gain the required knowledge for an extrapolation to engine condition. Simulated exit nitrogen oxide levels were validated with experimental emission data, and a global emission trend for the NO reduction at elevated pressure and constant turbine inlet temperature level was defined. The nitrogen oxide benefit at elevated operating pressure was justified with the significantly smaller flame surface area: the analysis of the simulated spanwise and top-view profiles showed a relatively short receded core flame with nitrogen oxide production in the center at high pressure relative to a longer and larger shear layer flame at atmospheric condition that produced NO toward the inner and outer side of the flame. Decomposition of the Damköhler number revealed the strong influence of the reaction timescales with higher reaction rates at elevated pressure, along with a moderate influence of the turbulent timescales, showing higher turbulence intensity in the lee-side recirculation zone at lower pressure.