Computations of turbulent lean premixed combustion using conditional moment closure

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k–? model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC.     

Kinetic parameters,collision rates,energy exchanges and transport coefficients of non-thermal electrons in premixed flames at sub-breakdown electric field strengths

The effects of an electric field on the collision rates, energy exchanges and transport properties of electrons in premixed flames are investigated via solutions to the Boltzmann kinetic equation. The case of high electric field strength, which results in high-energy, non-thermal electrons, is analysed in detail at sub-breakdown conditions. The rates of inelastic collisions and the energy exchange between electrons and neutrals in the reaction zone of the flame are characterised quantitatively. The analysis includes attachment, ionisation, impact dissociation, and vibrational and electronic excitation processes. Our results suggest that Townsend breakdown occurs for E/N = 140 Td. Vibrational excitation is the dominant process up to breakdown, despite important rates of electronic excitation of CO, CO₂ and N₂ as well as impact dissociation of O₂ being apparent from 50 Td onwards. Ohmic heating in the reaction zone is found to be negligible (less than 2% of peak heat release rate) up to breakdown field strengths for realistic electron densities equal to 10¹⁰ cm^?3. The observed trends are largely independent of equivalence ratio. In the non-thermal regime, electron transport coefficients are insensitive to mixture composition and approximately constant across the flame, but are highly dependent on the electric field strength. In the thermal limit, kinetic parameters and transport coefficients vary substantially across the flame due to the spatially inhomogeneous concentration of water vapour. A practical approach for identifying the plasma regime (thermal versus non-thermal) in studies of electric field effects on flames is proposed.     

Measurement of the structure coefficient of refractive index fluctuations in a turbulent premixed butane-air flame by means of a laser-based interferometer technique

With a view to measuring the structure coefficient of refractive index fluctuations in a turbulent premixed butane-air flame, a thin laser beam is sent into the flame perpendicular to the flow direction. The laser beam generally undergoes fluctuations of direction, phase, and amplitude. Only the random deflections of the laser beam may be taken into account. After having traversed the flame, the perturbed laser beam enters into an interferometric system. Materials and experimental procedure are described. In the unperturbed interference pattern, the zones only sensitive to fluctuations of the angle-of-arrival of the laser beam are detected. From the random displacements of the central bright fringe, the structure coefficient of refractive index fluctuations in the flame is measured. To prove that the method of measurement is satisfactory, the result obtained is applied for computing the power spectral density of the angle-of-arrival of the laser beam from the formula of correlations of the laser beam deflection angles which we have demonstrated in previous works. This computed power spectral density is compared to that measured from the effective position of the detector. A good agreement is observed between the two results.     

对冲火焰流场的二维分布对层流火焰传播速度测量的影响

在φ=0.6～0.75的条件下,利用PⅣ测得对冲火焰的二维速度场,分析使用偏移中心线一定距离提取的轴线速度分布曲线所得到的层流火焰传播速度SOL与沿中心线提取的速度分布曲线得到的SOL之间的误差大小.实验结果表明,在弱火焰时对冲火焰速度场具有明显的二维效应.φ=0.6,偏离中心轴线5 mm时,相对误差达到33%;而φ=...     

Stagnation point heat transfer from laminar,high temperature methane flames

Combustion of methane-rich fuels frequently provides forced convective heating in industry, and the ability to predict the rate of heat transfer from such flames to solid surfaces is often desirable. Mathematical modelling of stagnation point heat flux has been achieved by numerical solution of the boundary layer equations, and by an analytical equation modified to include the effects of chemical reaction in the free stream flow and to allow for the enhancement in heat flux caused by the diffusion and exothermic recombination of reactive species in the boundary layer surrounding the heat receiving body. Predictions from these models have been compared with experimental data obtained in high temperature methane flames of various equivalence ratios. Within the equilibrium region of these flames, predictions from the modified analytical equation based on total Lewis numbers equal to and greater than one form a tight envelope around the experimental results, and hence provide a relatively simple method of predicting heat flux. Numerical solutions tend to slightly underestimate predictions from the analytical equation and experimental data, although agreement with the alternative prediction method increases with the surface temperature of the heat receiving body     

Measurement of unstable burning velocities of iso-octane–air mixtures at high pressure and the derivation of laminar burning velocities

A new technique is reported for measuring burning velocities at high pressures in the final stages of two inwardly propagating flame kernels in an explosion bomb. The flames were initiated at diametrically opposite spark electrodes, close to the wall, in quiescent mixtures. Measurements of pressure and flame kernel propagation speeds by high-speed photography showed the burning velocities to be elevated above the corresponding laminar burning velocities as a result of the developing flame instabilities. The enhancement increased with increase in pressure and decreased with increase in Markstein number. When the Markstein number was negative, instabilities could be appreciable, as could the enhancement. For the iso-octane–air mixtures investigated, where the mixtures had well-characterised Markstein numbers or critical Peclet numbers at the relevant pressures and temperatures, it was possible to explain the enhancement quantitatively by the spherical explosion flame instability theory of Bechtold and Matalon, provided the critical Peclet number was that observed experimentally, and allowance was made for the changing pressure. With this theoretical procedure, it was possible to derive values of laminar burning velocity from the measured values of burning velocity over a wide range of equivalence ratios, pressures, and temperatures. The values became less reliable at the higher temperatures and pressures as the data on Markstein and critical Peclet numbers became less certain. It was found that with iso-octane as the fuel the laminar burning velocity decreased during isentropic compression.     

Structural response of different Lewis number premixed flames interacting with a toroidal vortex

Simultaneous measurements of temperature, CH* and OH* chemiluminescent species are carried out to explore the impact of stretch rate and curvature on the structure of premixed flames. The configuration of an initially flat premixed flame interacting with a toroidal vortex is selected for the present study and reasons for this choice are discussed. Lewis number effects are assessed by comparing methane and propane flames. It is emphasized that the flame structure experiences very strong variations. In particular, the flame is shrunk (broadened) in the initial (final) period of the interaction with the vortex where strain rate (curvature) contribution of the stretch rate is predominant. By further analysing independently the thickness of the preheat and reaction zones, it is shown that for propane flames, not only the former but also the latter is significantly altered in zones where the flame curvature is negative. Changes in the reaction zone properties are further emphasized using CH* and OH* radicals. It is demonstrated that higher thermal diffusivity plays a significant role around curved regions, in which the enhanced diffusion of heat leads to a strong increase of CH* compared to OH* intensity. As an overall conclusion, this study suggests that it would be interesting to reassess the internal flame structure at lower and moderate Karlovitz numbers since changes might appear for a moderate vortex intensity with typical size much larger than the flame thickness.     

Experimental study of a confined premixed metal combustor: Metal flame stabilization dynamics and nitrogen oxides production

This work presents a study of a magnesium/air combustion process in the context of innovative zero carbon dioxide (CO₂) energy carriers for reducing global warming effects. In order to analyze more deeply the confined combustion of magnesium under fluctuating overpressure conditions (0 to 24 hPa) and the generated gaseous by-products, magnesium/air flames have been realized in a combustion chamber with a conical bluff-body as flame holder and different contraction ratios diaphragms at the exit duct. Sieved magnesium samples with two size-fractions were tested: 20–50?µm and 50–70?µm. The gaseous emissions of nitrogen oxides (NO_x) and dioxygen (O₂) were analyzed with on-line infrared, ultraviolet and paramagnetic analyzers. A flame pulsating behavior was clearly observed from light emission intensity (monitored by a photodiode) and pressure fluctuations (monitored by a pressure sensor); the frequencies obtained ranged between 3 and 10?Hz. The frequency of the pulsation exhibited strong dependence on the geometric configuration of the chamber: a contraction diaphragm divided by two the frequency level of the fluctuations in the studied range of maximum overpressure. Such fluctuations may probably be the consequence of periodic perturbations of the recirculation zone behind the bluff-body. These periodic perturbations are themselves caused by strong periodic overpressure variations due to stiff contraction downstream responding to gas velocity fluctuations. This feed-back-loop mechanism was considered in this study. NO_x emissions produced through the thermal pathway were analyzed for equivalence ratios ranging from 0.29 to 1. The representation of NO_x versus equivalence ratio exhibited a parabolic shape with a maximum for an equivalence ratio of 0.4. Moreover, NO_x emissions of this metal combustor have shown a similar order of magnitude than current internal combustion engines.