Investigation of mass and energy coupling between soot particles and gas species in modelling ethylene counterflow diffusion flames

A numerical model is developed aiming at investigating soot formation in ethylene counterflow diffusion flames. The mass and energy coupling between soot solid particles and gas-phase species is investigated in detail. A semi-empirical two-equation model is chosen for predicting soot mass fraction and number density. The model describes particle nucleation, surface growth, and oxidation. A detailed kinetic mechanism is considered for the gas phase and the effect of considering radiation heat losses is also evaluated. Simulations were done for a range of conditions that produce low-to-significant amounts of soot using three strategies: first by changing the strain rate imposed on the flow field, second, by changing the oxygen content in the oxidant stream, and third, by changing the pressure. Additionally, the effect of the transport model chosen was analysed. The results showed that, for the flames studied and within the limits of the present work, the soot and gas radiation terms are of primary importance for numerical simulations. Additionally, it was shown that the soot mass and thermodynamic properties coupling terms are, in general, a second-order effect, with an importance that increases as soot amount increases. As a general recommendation, the radiation terms have always to be considered, whereas full coupling has to be employed only when the soot mass fraction, Y_S, is equal to or larger than 0.008. If a higher precision is required, with errors less than 1%, full coupling should be taken into account for Y_S ≥ 0.002. For lower soot amounts, the coupling through soot mass and thermodynamic properties may be neglected as a first approximation, but an error on the total mass conservation will be present. Additionally, discrepancies from considering different transport models (detailed or simplified) are larger than those found from not fully coupling the phases.     

The effects of strain and curvature on the mass burning rate of premixed laminar flames

The Markstein number characterizes the effect that flame stretch has on the burning velocity. Different expressions for this number are deduced from integral analysis. According to a phenomenological law, the Markstein number can be separated into a part for the curvature of the flame and a part for the straining of the flow. This separation is analysed here. It appears that the Markstein number for curvature and the combined one for both curvature and strain are unique. It is, however, not possible to introduce a separate and unique Markstein number for the flow straining that can be used to describe its influence in different combustion situations. The theoretical and numerical analysis is applied to flat steady counterflow flames as well as to steady, imploding and expanding spherical flames.     

Effects of compression and stretch on the determination of laminar flame speeds using propagating spherical flames

The effects of flow compression and flame stretch on the accurate determination of laminar flame speeds at normal and elevated pressures using propagating spherical flames at constant pressure or constant volume are studied theoretically and numerically. The results show that both the compression-induced flow motion and flame stretch have significant impacts on the accuracy of flame speed determination. For the constant pressure method, a new method to obtain a compression-corrected flame speed (CCFS) for nearly constant pressure spherical bomb experiments is presented. Likewise, for the constant volume method, a technique to obtain a stretch-corrected flame speed (SCFS) at elevated pressures and temperatures is developed. The validity of theoretical results for both constant pressure and constant volume methods is demonstrated by numerical simulations using detailed chemistry for hydrogen/air, methane/air, and propane/air mixtures. It is shown that the present CCFS and SCFS methods not only improve the accuracy of the flame speed measurements significantly but also extend the parameter range of experimental conditions. The results can be used directly in experimental measurements of laminar flame speeds.     

Effect of Lewis number on premixed laminar lean-limit flames stabilized on a bluff body

In this paper, we present a study on the effect of Lewis number, Le, on the stabilization and blow-off of laminar lean limit premixed flames stabilized on a cylindrical bluff body. Numerical simulations and experiments are conducted for propane, methane and two blends of hydrogen with methane as fuel gases, containing 20% and 40% of hydrogen by volume, respectively. It is found that the Le?>?1 flame blows-off via convection from the base of the flame (without formation of a neck) when the conditions for flame anchoring are not fulfilled. Le?≤?1 flames exhibit a necking phenomenon just before lean blow-off. This necking of the flame front is a result of the local reduction in mass burning rates causing flame merging and quenching of the thin flame tube formed. The structure of these flames at the necking location is found to be similar to tubular flames. It is found that extinction stretch rates for tubular flames closely match values at the neck location of bluff-body flames of corresponding mixtures, suggesting that excessive flame stretch is directly responsible for blow-off of the studied Le?≤?1 flames. After quenching of the neck, the upstream part forms a steady and stable residual flame in the wake of the bluff body while the downstream part is convected away.     

Effect of Soret diffusion on lean hydrogen/air flames at normal and elevated pressure and temperature

The influence of Soret diffusion on lean premixed flames propagating in hydrogen/air mixtures is numerically investigated with a detailed chemical and transport models at normal and elevated pressure and temperature. The Soret diffusion influence on the one-dimensional (1D) flame mass burning rate and two-dimensional (2D) flame propagating characteristics is analysed, revealing a strong dependency on flame stretch rate, pressure and temperature. For 1D flames, at normal pressure and temperature, with an increase of Karlovitz number from 0 to 0.4, the mass burning rate is first reduced and then enhanced by Soret diffusion of H₂ while it is reduced by Soret diffusion of H. The influence of Soret diffusion of H₂ is enhanced by pressure and reduced by temperature. On the contrary, the influence of Soret diffusion of H is reduced by pressure and enhanced by temperature. For 2D flames, at normal pressure and temperature, during the early phase of flame evolution, flames with Soret diffusion display more curved flame cells. Pressure enhances this effect, while temperature reduces it. The influence of Soret diffusion of H₂ on the global consumption speed is enhanced at elevated pressure. The influence of Soret diffusion of H on the global consumption speed is enhanced at elevated temperature. The flame evolution is more affected by Soret diffusion in the early phase of propagation than in the long run due to the local enrichment of H₂ caused by flame curvature effects. The present study provides new insights into the Soret diffusion effect on the characteristics of lean hydrogen/air flames at conditions that are relevant to practical applications, e.g. gas engines and turbines.     

A DNS study of the physical mechanisms associated with density ratio influence on turbulent burning velocity in premixed flames

Data obtained in 3D direct numerical simulations of statistically planar, 1D weakly turbulent flames characterised by different density ratios σ are analysed to study the influence of thermal expansion on flame surface area and burning rate. Results show that, on the one hand, the pressure gradient induced within a flame brush owing to heat release in flamelets significantly accelerates the unburned gas that deeply intrudes into the combustion products in the form of an unburned mixture finger, thus causing large-scale oscillations of the burning rate and flame brush thickness. Under the conditions of the present simulations, the contribution of this mechanism to the creation of the flame surface area is substantial and is increased by σ, thus implying an increase in the burning rate by σ. On the other hand, the total flame surface areas simulated at σ = 7.53 and 2.5 are approximately equal. The apparent inconsistency between these results implies the existence of another thermal expansion effect that reduces the influence of σ on the flame surface area and burning rate. Investigation of the issue shows that the flow acceleration by the combustion-induced pressure gradient not only creates the flame surface area by pushing the finger tip into the products, but also mitigates wrinkling of the flame surface (the side surface of the finger) by turbulent eddies. The latter effect is attributed to the high-speed (at σ = 7.53) axial flow of the unburned gas, which is induced by the axial pressure gradient within the flame brush (and the finger). This axial flow acceleration reduces the residence time of a turbulent eddy in an unburned zone of the flame brush (e.g. within the finger). Therefore, the capability of the eddy for wrinkling the flamelet surface (e.g. the side finger surface) is weakened owing to a shorter residence time.