Computational and experimental study of oxygen-enhanced axisymmetric laminar methane flames

Three axisymmetric laminar coflow diffusion flames, one of which is a nitrogen-diluted methane/air flame (the ‘base case’) and the other two of which consist of nitrogen-diluted methane vs. pure oxygen, are examined both computationally and experimentally. Computationally, the local rectangular refinement method is used to solve the fully coupled nonlinear conservation equations on solution-adaptive grids. The model includes C2 chemistry (GRI 2.11 and GRI 3.0 chemical mechanisms), detailed transport, and optically thin radiation. Because two of the flames are attached to the burner, thermal boundary conditions at the burner surface are constructed from smoothed functional fits to temperature measurements. Experimentally, Raman scattering is used to measure temperature and major species concentrations as functions of the radial coordinate at various axial positions. As compared to the base case flame, which is lifted, the two oxygen-enhanced flames are shorter, hotter, and attached to the burner. Computational and experimental flame lengths show excellent agreement, as do the maximum centreline temperatures. For each flame, radial profiles of temperature and major species also show excellent agreement between computations and experiments, when plotted at fixed values of a dimensionless axial coordinate. Computational results indicate peak NO levels in the oxygen-enhanced flames to be very high. The majority of the NO in these flames is shown to be produced via the thermal route, whereas prompt NO dominates for the base case flame.     

Detailed numerical simulation of laminar flames by a parallel multiblock algorithm using loosely coupled computers

A parallel algorithm for the detailed multidimensional numerical simulation of laminar flames able to work efficiently with loosely coupled computers is described. The governing equations have been discretized using the finite volume technique over staggered grids. A SIMPLE-like method has been employed to solve the velocity–pressure fields while the species equations have been calculated in a segregated manner using an operator-splitting technique. The domain decomposition method is used to optimize the domain's discretization and to parallelize the code. The main attributes and limitations, together with the computational features (computational effort, parallel performance, memory requirements, etc), are shown, taking into account different degrees of chemical modelling and two benchmark problems: a premixed methane/air laminar flat flame and a confined co-flow non-premixed methane/air laminar flame. In order to assess the validity of the numerical solutions, a post-processing procedure, based on the generalized Richardson extrapolation for h-refinement studies and on the grid convergence index, has been used.     

An adaptive multilevel local defect correction technique with application to combustion

The standard local defect correction (LDC) method has been extended to include multilevel adaptive gridding, domain decomposition and regridding. The domain decomposition algorithm provides a natural route for parallelization by employing many small tensor-product grids, rather than a single large unstructured grid; this algorithm can greatly reduce memory usage. The above properties are illustrated by successfully applying the new algorithm to a simple heat transfer problem with an analytical solution, and by subsequently solving the more complex problem of an axisymmetric laminar Bunsen flame with one-step chemistry. The simulation data show excellent agreement with results previously published in the literature.     

Parallelization strategies for an implicit Newton-based reactive flow solver

The solution of reactive flows using fully implicit methods on distributed memory machines is investigated in detail. Three different parallel implementations of Newton's method are described and tested on the solution of two-dimensional laminar axisymmetric coflow diffusion flames. Each implementation has different computational requirements, both in the amount of communication among the processes and in the computational overhead due to the calculation of physical quantities at the interfaces between subdomains. An effective trade-off is established between communications and calculations so that the most communication-intensive implementation results in computational speedup only if the network is sufficiently fast.

Benchmark results are presented for a variety of chemical mechanisms, grid decomposition techniques, and hardware. Parallelization efficiencies of about 80% and speedups of 20–100 are reported for most test cases. The method developed here is well suited for complex chemistry problems with very large mechanisms; in particular, the numerical solution of a laminar axisymmetric JP-8/air coflow diffusion flame with a 222-species mechanism is made possible using this approach.     

How fast can we burn, 2.0

We review the state of the art in measurements and simulations of the behavior of premixed laminar and turbulent flames, subject to differential diffusion, stretch and curvature. The first part of the paper reviews the behavior of premixed laminar flames subject to flow stretch, and how it affects the accuracy of measurements of unstrained laminar flame speeds in stretched and spherically propagating flames. We then examine how flow field stretch and differential diffusion interact with flame propagation, promoting or suppressing the onset of thermodiffusive instabilities. Secondly, we survey the methodology for and results of measurements of turbulent flame speeds in the light of theory, and identify issues of consistency in the definition of mean flame speeds, and their corresponding mean areas. Data for methane at a single operating condition are compared for a range of turbulent conditions, showing that fundamental issues that have yet to be resolved for Bunsen and spherically propagating flames. Finally, we consider how the laminar flame scale response of flames to flow perturbations interacting with differential diffusion leads to very different outcomes to the overall sensitivity of the burning rate to turbulence, according to numerical simulations (DNS). The paper concludes with opportunities for future measurements and model development, including the perennial recommendation for robust archival databases of experimental and DNS results for future testing of models.     

A computational framework for predicting laminar reactive flows with soot formation

Numerical modeling is an attractive option for cost-effective development of new high-efficiency, soot-free combustion devices. However, the inherent complexities of hydrocarbon combustion require that combustion models rely heavily on engineering approximations to remain computationally tractable. More efficient numerical algorithms for reacting flows are needed so that more realistic physics models can be used to provide quantitative soot predictions. A new, highly-scalable combustion modeling tool has been developed specifically for use on large multiprocessor computer architectures. The tool is capable of capturing complex processes such as detailed chemistry, molecular transport, radiation, and soot formation/destruction in laminar diffusion flames. The proposed algorithm represents the current state of the art in combustion modeling, making use of a second-order accurate finite-volume scheme and a parallel adaptive mesh refinement (AMR) algorithm on body-fitted, multiblock meshes. Radiation is modeled using the discrete ordinates method (DOM) to solve the radiative transfer equation and the statistical narrow-band correlated-k (SNBCK) method to quantify gas band absorption. At present, a semi-empirical model is used to predict the nucleation, growth, and oxidation of soot particles. The framework is applied to two laminar coflow diffusion flames which were previously studied numerically and experimentally. Both a weakly-sooting methane–air flame and a heavily-sooting ethylene–air flame are considered for validation purposes. Numerical predictions for these flames are verified with published experimental results and the parallel performance of the algorithm analyzed. The effects of grid resolution and gas-phase reaction mechanism on the overall flame solutions were also assessed. Reasonable agreement with experimental measurements was obtained for both flames for predictions of flame height, temperature and soot volume fraction. Overall, the algorithm displayed excellent strong scaling performance by achieving a parallel efficiency of 70% on 384 processors. The proposed algorithm proved to be a robust, highly-scalable solution method for sooting laminar flames.     

On polyhedral structures of lean methane/hydrogen Bunsen flames: Combined experimental and numerical analysis

In premixed flame propagation of lean hydrogen or hydrogen-enriched blends, both hydrodynamic and thermo-diffusive instabilities are governing the flame front shape and affect its propagation velocity. As a result, different types of cellular patterns can occur along the flame front in a laminar scenario. In this context, an interesting phenomenon is the formation of polyhedral flames which can be observed in a Bunsen burner. It is the objective of this work to systematically characterize the polyhedral structures of premixed methane/hydrogen Bunsen flames in a combined experimental and numerical study. A series of lean flames with hydrogen content varying between 20 and 85% at two equivalence ratios is investigated. The experiments encompass chemiluminescence imaging together with Planar Laser-induced Fluorescence (PLIF) measurements of the OH radical. Characteristic cell sizes are quantified from the experiments and related to the characteristic length scales obtained from a linear stability analysis. In the experiments, it is observed that the cell sizes at the base of the polyhedral Bunsen flames decrease almost linearly with hydrogen addition and only a weak dependence on the equivalence ratio is noted. These trends are well reflected in the numerical results and the length scale comparison further shows that the wavelength with the maximum growth rate predicted by the linear stability analysis is comparable to the cell size obtained from the experiment. The correlation between the experimental findings and the linear stability analysis is discussed from multiple perspectives considering the governing time and length scales, furthermore drawing relations to previous studies on cellular flames.     

Emission spectroscopy of flame fronts in aluminum suspensions

Spatially resolved emission spectra from Bunsen-type flames stabilized in aluminum suspensions in air and oxygen–argon/helium mixtures were obtained using a mechanical-optical scanning system. A low resolution (1.5 nm) spectrometer was used to acquire the broad spectra over the 350–1000 nm range, and a high-resolution (0.04 nm) instrument was used for observation of AlO molecular bands and non-ionized atomic aluminum. The temperature of condensed phase emitters in the flame was derived using polychromatic fitting of the continuum spectra to Planck’s law. AlO temperature was found by fitting of the theoretically calculated shape of the band to experimental data. Peak temperatures of the condensed emitters were found to be approximately 3250 K in aluminum-air flames and approximately 3350 K for oxygen–argon/helium flames. Temperatures derived from AlO spectra coincide with the temperature of the condensed emitters with measurement accuracy and are only 100–200 °C lower than the computed equilibrium flame temperatures. The radial distribution of the temperature profile of the continuous emitters was found via Abel deconvolution and recovered the double-front structure of the Bunsen flame cone, with the outer flame being attributed to a diffusion flame of the fuel-rich products with ambient air. The observation of atomic aluminum lines seen in emission from the outer flame edge and partial self-absorption from the inner flame confirms the structure associated with the double-front structure. The implications of these results for the regime of particle combustion in a dust flame are discussed.     

Characteristics of microjet methane diffusion flames

Characteristics of microjet methane diffusion flames stabilized on top of the vertically oriented, stainless-steel tubes with an inner diameter ranging from 186 to 778 μ m are investigated experimentally, theoretically and numerically. Of particular interest are the flame shape, flame length and quenching limit, as they may be related to the minimum size and power of the devices in which such flames would be used for future micro-power generation. Experimental measurements of the flame shape, flame length and quenching velocity are compared with theoretical predictions as well as detailed numerical simulations. Comparisons of the theoretical predictions with measured results show that only Roper's model can satisfactorily predict the flame height and quenching velocity of microjet methane flames. Detailed numerical simulations, using skeletal chemical kinetic mechanism, of the flames stabilized at the tip of d = 186, 324 and 529 μ m tubes are performed to investigate the flame structures and the effects of burner materials on the standoff distance near extinction limit. The computed flame shape and flame length for the d = 186 μm flame are in excellent agreement with experimental results. Numerical predictions of the flame structures strongly suggest that the flame burns in a diffusion mode near the extinction limit. The calculated OH mass fraction isopleths indicate that different tube materials have a minor effect on the standoff distance, but influence the quenching gap between the flame and the tube.     

Large scale effects in weakly turbulent premixed flames

In this study we numerically investigate large scale premixed flames in weakly turbulent flow fields. A large scale flame is classified as such based on a reference hydrodynamic lengthscale being larger than a neutral (cutoff) lengthscale for which the hydrodynamic or Darrieus–Landau (DL) instability is balanced by stabilizing diffusive effects. As a result, DL instability can develop for large scale flames and is inhibited otherwise. Direct numerical simulations of both large scale and small scale three-dimensional, weakly turbulent flames are performed at constant Karlovitz and turbulent Reynolds number, using two paradigmatic configurations, namely a statistically planar flame and a slot Bunsen flame. As expected from linear stability analysis, DL instability induces its characteristic cusp-like corrugation only on large scale flames. We therefore observe significant morphological and topological differences as well as DL-enhanced turbulent flame speeds in large scale flames. Furthermore, we investigate issues related to reaction rate modeling in the context of flame surface density closure. Thicker flame brushes are observed for large scale flames resulting in smaller flame surface densities and overall larger wrinkling factors.     

Methane flames in a jet impinging onto a wall

Traditionally, research has focused on positive stretch in the stagnation flow and negative stretch along the Bunsen flame. Only a very limited amount of research has been devoted to studying the behavior of a conical Bunsen flame established in a stagnation flow, which is significantly affected by the combined effects of the curvature stretch and the aerodynamic straining. This investigation is aimed at studying the characteristics of laminar conical premixed flames in an impinging jet flow experimentally and theoretically. First, we analyze the transport processes of a nonreactive impinging jet flow numerically. For lower burner-to-plate distance, the potential core becomes concave at the top. Hence, a conical Bunsen flame established in such a flow field may suffer positive flow stretch. The predicted flame shapes using a simple model incorporated with the numerical results agree well with the experimental observations. Flame shapes exhibit double-solution characteristics in a certain range of methane concentrations. Experimentally, by following different paths of adjusting methane concentration (decreasing from rich to lean or increasing from lean to rich), two different flame configurations (planar or conical flame) may exist at the same flow conditions, namely burner-to-plate distance, inlet velocity, and methane concentration. At the higher (or lower) critical methane concentration, the transition from a flat flame to a conical flame (or from a conical flame to a flat flame) occurs. The calculation of stretch and measurement of flame temperature for the low inlet velocity, 0.8 m/s, show that the stretch of a conical flame established in a stagnation flow is negative (dominated by the flame curvature). However, it is important to emphasize that at high velocity, e.g., U_in = 1.6 m/s, a negatively stretched flame tip can suffer positive flow stretch. This significant finding has been verified in the experiment since the conical flame tip is higher than the positively stretched flat flame.     

On the existence of steady-state gaseous microgravity spherical diffusion flames in the presence of radiation heat loss

Whether steady-state gaseous microgravity spherical diffusion exist in the presence of radiation heat loss is an important fundamental question and has important implications for spacecraft fire safety. In this work, experiments aboard the International Space Station and a transient numerical model are used to investigate the existence of steady-state microgravity spherical diffusion flames. Gaseous spherical diffusion flames stabilized on a porous spherical burner are employed in normal (i.e., fuel flowing into an ambient oxidizer) and inverse (i.e., oxidizer flowing into an ambient fuel) flame configurations. The fuel is ethylene and the oxidizer oxygen, both diluted with nitrogen. The flow rate of the reactant gas from the burner is held constant. It is found that steady-state gaseous microgravity spherical diffusion flames can exist in the presence of radiation heat loss, provided that the steady-state flame size is less than the flame size for radiative extinction, and the flame develops fast enough that radiation heat loss does not drop the flame temperature below the critical temperature for radiative extinction (1130 K). A simple model is provided that allows for the identification of initial conditions that can lead to steady-state spherical diffusion flames. In the spherical, infinite domain configuration, the characteristic time for the diffusion-controlled system to effectively reach steady-state is found to be on the order of 100,000 s. Despite a narrow range of attainable conditions, flames that exhibit steady-state behavior are observed aboard the ISS for up to 870 s, even with the constraint of a finite boundary. Steady-state flames are simulated using the numerical model for over 100,000 s.     

Measurements of three-dimensional mean flame surface area ratio in turbulent premixed Bunsen flames

A data processing scheme with particular emphasis on proper flame contour smoothing is developed and applied to measure the three-dimensional mean flame surface area ratio in turbulent premixed flames. The scheme is based on the two-sheet imaging technique such that the mean flame surface area ratio is an average within a window covering a finite section of the turbulent flame brush. This is in contrast to the crossed-plane tomograph technique which applies only to a line. Two sets of Bunsen flames have been investigated in this work with the turbulent Reynolds number up to 4000 and the Damköhler number ranging from less than unity to close to 10. The results show that three-dimensional effects are substantial. The measured three-dimensional mean flame surface area ratio correlates well with a formula similar to the Zimont model for turbulent burning velocity but with different model constants. Also, the mean flame surface area ratio displays a weak dependency on turbulence intensity but a strong positive dependency on the turbulence integral length scale.     

Modulation of a methane Bunsen flame by upstream perturbations

In this paper the effects of an upstream spatially periodic modulation acting on a turbulent Bunsen flame are investigated using direct numerical simulations of the Navier-Stokes equations coupled with the flamelet generated manifold (FGM) method to parameterise the chemistry. The premixed Bunsen flame is spatially agitated with a set of coherent large-scale structures of specific wave-number, K. The response of the premixed flame to the external modulation is characterised in terms of time-averaged properties, e.g. the average flame height ?H? and the flame surface wrinkling ?W?. Results show that the flame response is notably selective to the size of the length scales used for agitation. For example, both flame quantities ?H? and ?W? present an optimal response, in comparison with an unmodulated flame, when the modulation scale is set to relatively low wave-numbers, 4π/L ? K ? 6π/L, where L is a characteristic scale. At the agitation scales where the optimal response is observed, the average flame height, ?H?, takes a clearly defined minimal value while the surface wrinkling, ?W?, presents an increase by more than a factor of 2 in comparison with the unmodulated reference case. Combined, these two response quantities indicate that there is an optimal scale for flame agitation and intensification of combustion rates in turbulent Bunsen flames.     

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

多介质流界面高精度自适应欧拉算法

采用欧拉网格自适应算法捕捉多介质流界面,获得了高精度界面特征,对不同物质引入不同位标函数跟踪界面运动,将位标函数方程与流体动力学方程非耦合求解,在笛卡尔坐标系中运用二阶精度有限体积算法,在保持流场守恒条件下,采用多层网格级对笛卡尔网格嵌套细化,实现了多介质流物质界面的高精度自适应跟踪.方法逻辑简单,大大节省了CPU时间,且能够对局部参数急剧变化的流场(如激波)进行自适应跟踪.     

Euler及N-S方程的二维四边形/叉树混合网格自适应算法

描述了一种基于直角叉树网格的Euler和N-S方程自适应算法.由于考虑了粘性的作用,提出并使用了四边形叉树混合网格的方法,在几何表面附近生成贴体的四边形网格,外流场使用直角叉树网格.采用中心有限体积法,对Euler及N-S方程进行数值求解,对N-S方程的计算中加入了B-L代数湍流模型.在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离等特性的捕捉和分辨能力.采用上述方法,数值分析了单段和多段翼型的绕流问题.     

Estimation of laminar flame speeds using axisymmetric bunsen flames: Molecular transport effects

Intricacies associated with the estimation of laminar flame speed using the axisymmetric Bunsen flame technique were assessed, through parametric direct numerical simulations. The study involved methane-air mixtures at atmospheric pressure and temperature, and both the flame cone angle and flame surface area methods were utilized to estimate the laminar flame speeds based on conditions used in recent relevant experimental studies. The results provided insight into the details of the flame structure and allowed for the assessment of various non-idealities and the attendant uncertainties associated with the estimation of laminar flame speeds. Additionally, molecular transport effects were investigated by altering the fuel diffusivity, in order to evaluate its impact on the flame structure and propagation under the presence of negative stretch. The modification of fuel diffusivity was found to affect the burning rate as stretch varies. Under fuel rich conditions, decreasing the fuel diffusivity was found to have an opposite effect on the heat release and thus the burning rate, when compared to positively stretched flames that have been investigated recently in a similar manner. The reported results are expected to provide guidance in flame propagation experiments using the convenient Bunsen flame method at near-atmospheric or elevated pressures, as well as insight into the effects of negative stretch that has, compared to positive, attracted less attention in past studies.     

Numerical simulation of flames as gas-dynamic discontinuities

The dynamics of thin premixed flames is computationally studied within the context of a hydrodynamic theory. A level-set method is used to track down the flame, which is treated as a free-boundary interface. The flow field is described by the incompressible Navier–Stokes equations, with different densities for the burnt and unburnt gases, supplemented by singular source terms that properly account for thermal expansion effects. The numerical scheme has been tested on several benchmark problems and was shown to be stable and accurate. In particular, the propagation of a planar flame front and the dynamics of hydrodynamically unstable flames were successfully simulated. This includes recovering the planar front in narrow domains, the Darrieus–Landau linear growth rate for long waves of small amplitude, and the nonlinear development of cusp-like structures predicted by the Michelson–Sivashinsky equation for a small density change. The stationary flame of a Bunsen burner with uniform and parabolic outlet flows were also simulated, showing in particular a careful mapping of the flow field. Finally, the evolution of a hydrodynamically unstable flame was studied for finite amplitude disturbances and realistic values of thermal expansion. These results, which constitute one of the main objectives of this study, elucidate the effect of thermal expansion on flame dynamics.     

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.