Flame dynamics of azimuthal forced spinning and standing modes in an annular combustor

Azimuthal forcing has been applied to flames in a laboratory scale annular combustor in order to accurately control the azimuthal mode of excitation. A new forcing configuration permitted not only the pressure amplitude, but also the spin ratio and mode orientation to be accurately controlled, in order to generate standing modes and for the first time strong spinning modes in both a clockwise (CW) and anti-clockwise (ACW) direction. The phase averaged heat release dynamics of these modes was compared and a number of differences observed depending on the direction of pressure wave propagation, demonstrating characteristic ACW and CW heat release patterns. A new spin compensating averaging method was then introduced to analyse the flame dynamics, and it was shown that through the application of this method the dynamics of standing wave oscillations could be decomposed to recover the characteristic ACW and CW heat release responses. The global heat release response was also assessed during strongly spinning modes, and the magnitude of the response was shown to depend strongly on the direction of propagation, demonstrating the importance of the local swirl direction on the global heat release response, with important implications for the modelling of such flows.     

Flame sheet dynamics of bluff-body stabilized flames during longitudinal acoustic forcing

Bluff-body stabilized flames are susceptible to combustion instabilities due to interactions between acoustics, vortical disturbances, and the flame. In order to elucidate these flow-flame interactions during an instability, an experimental and computational investigation of the flame-sheet dynamics of a harmonically excited flame was performed. It is shown that the flame dynamics are controlled by three key processes: excitation of shear layer instabilities by the axially oscillating flow, anchoring of the flame at the bluff body, and the kinematic response of the flame to this forcing. The near-field flame features are controlled by flame anchoring and the far-field by kinematic restoration. In the near-field, the flame response grows with downstream distance due to flame anchoring, which prevents significant flame movement near the attachment point. Theory predicts that this results in linear flame response characteristics as a function of perturbation amplitude, and a monotonic growth in magnitude of the flame-sheet fluctuations near the stabilization point, consistent with the experimental data. Farther downstream, the flame response reaches a maximum and then decays due to the dissipation of the vortical disturbances and action of flame propagation normal to itself, which acts to smooth out the wrinkles generated by the harmonic flow forcing. This behavior is strongly non-linear, resulting in significant variation in far-field flame-sheet response with perturbation amplitude.     

Impact of swirl and bluff-body on the transfer function of premixed flames

The frequency response of three lean methane/air flames submitted to flowrate perturbations is analyzed for flames featuring the same equivalence ratio and thermal power, but a different stabilization mechanism. The first flame is stabilized by a central bluff body without swirl, the second one by the same bluff body with the addition of swirl and the last one only by swirl without central insert. In the two last cases, the swirl level is roughly the same. These three flames feature different shapes and heat release distributions, but their Flame Transfer Function (FTF) feature about the same phase lag at low frequencies. The gain of the FTF also shows the same behavior for the flame stabilized by the central insert without swirl and the one fully aerodynamically stabilized by swirl. Shedding of vortical structures from the injector nozzle that grow and rollup the flame tip controls the FTF of these flames. The flame stabilized by the swirler-plus-bluff-body system features a peculiar response with a large drop of the FTF gain around a frequency at which large swirl number oscillations are observed. Velocity measurements in cold flow conditions reveal a strong reduction of the size of the vortical structures shed from the injector lip at this forcing condition. The flame stabilized aerodynamically only by swirl and the one stabilized by the bluff body without swirl do not exhibit any FTF gain drop at low frequencies. In the former case, large swirl number oscillations are still identified, but large vortical structures shed from the nozzle also persist at the same forcing frequency in the cold flow response. These different flame responses are found to be intimately related to the dynamics of the internal recirculation region, which response strongly differs depending upon the injector used to stabilize the flame.     

Hydrodynamic and chemical scaling for blow-off dynamics of lean premixed flames stabilized on a meso-scale bluff-body

Direct numerical simulations were conducted to investigate the effect of two parameters, density ratio and laminar flame speed, on the conditions of the onset of local extinction and blow-off of lean premixed flames, stabilized on a meso-scale bluff-body in hydrogen-air and syngas–air mixtures. A total of six simulation cases were considered as isolated comparison of the two parametric effects of the fluid dynamic instability and flame time scale. For all cases under study, the general flame development towards the blow-off limit showed a sequence of five distinct modes, with possible cyclic patterns among the different modes for a range of velocity conditions. The onset of local extinction was observed during the asymmetric vortex shedding and vortex street mode. As the density ratio is decreased, the flow inunder reviewstability is promoted through the increased sinuous mode, and such behavior was properly scaled by the Strouhal number. Although the blow-off velocity is altered by the fluid dynamic effects, the condition for the onset of local extinction and blow-off was mainly dictated by the competition between flow residence time associated with the lateral flame motion and ignition delay of the local mixtures. Time scale analysis supported the validity of the findings across all the cases investigated.     

Unsteady propagation of premixed methane/propane flames in a mesoscale disk burner of variable-gaps

Unsteady flame propagation of air-premixed methane and propane flames was investigated in a new mesoscale disk burner, of which disk-gap could be precisely varied. To begin with, the quenching disk-gaps on the flammability limits were measured. In most cases, with the slight increase of the disk-gap, cellular flame structures could be generated. The initiation of such cellular structures could be explained by the thermally induced hydrodynamic instability, and it could be enhanced if the Lewis number was sufficiently small. When the disk-gap was sufficiently larger than a critical value (approximately 1.5 times the quenching distance), the cellular flame structure changed into a smooth one in the azimuthal direction. With a further increase of the disk-gap, the flame propagation velocities approached to constant values. These values were comparable to the laminar burning velocities except for the propane-rich conditions, in which much larger propagation velocities were observed. The flame stretch effects (coupled with Le-effects) within a narrow space were suspected as the reason. The structural transition of the premixed flame could be investigated successfully through various disk-gaps, from the smallest quenching-scale to the ordinary large scales via whole mesoscales including Hele–Shaw scales.     

Leading point behavior during the ignition of an annular combustor with liquid n-heptane injectors

Experimental and numerical investigations of ignition in combustors with multiple burners have recently emerged and have provided new insights on the last phase of ignition in gas turbine-like annular geometries where the flame propagates from burner to burner. Previous comparisons between calculations and experiments of light-round in a laboratory scale annular combustion chamber have demonstrated the ability of large-eddy simulation to predict such processes for perfectly premixed conditions and, more recently, for n-heptane spray injection. The present analysis focuses on two additional operating points with liquid n-heptane sprays and the turbulent flame propagation in the two-phase mixture is examined through the behavior of its leading points. The validation of the light-round process is characterized in terms of ignition delays. The detailed analysis of the propagation through the definition of a leading point enables to highlight some key phenomena responsible for the flame behavior, such as the influence of the liquid droplet spray and its vaporization in the chamber. Calculations indicate that the volumetric expansion due to the chemical reaction at the flame induces a strong azimuthal flow in the fresh stream at a distance of several sectors ahead of the flame, which modifies conditions in this region. This creates heterogeneities in the gas composition and wakes on the downstream side of the swirling jets formed by the injectors, with notable effects on the motion of the leading point and on the absolute flame velocity.     

Nonlinear hydrodynamic and thermoacoustic oscillations of a bluff-body stabilised turbulent premixed flame

Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.     

Measurements of flame orientation and scalar dissipation in turbulent partially premixed methane flames

Three turbulent flames were studied using a new experimental facility developed at Sandia National Laboratories. Line imaging of Raman and Rayleigh scattering and CO laser-induced fluorescence (LIF) yielded information on all major species, temperature, mixture fraction, and a 1D surrogate measure of scalar dissipation. Simultaneously, crossed planar OH LIF imaging provided information on the instantaneous flame orientation, allowing estimation of the full 3D (flame-normal) scalar dissipation rate. The three flames studied were methane–air piloted jet flames (Sandia flames C, D, and E), which cover a range in Reynolds number from 13,400 to 33,600. The statistics of the instantaneous flame orientation are examined in the different flames, with the purpose of studying the prevailing kinematics of isoscalar contours. The 1D and 3D results for scalar dissipation rate are examined in detail, both in the form of conditional averages and in the form of probability density functions. The effect of overall strain and Reynolds number on flame suppression and eventual extinction is also investigated, by examining the doubly conditional statistics of temperature in the form of S-shaped curves. This latter analysis reveals that double conditioning of temperature on both mixture fraction and scalar dissipation does not collapse the data from these flames onto the same curve at low scalar dissipation rates, as might be expected from simple flamelet concepts.     

突扩燃烧室湍流预混反应流的数值模拟

对突扩燃烧室这一典型工程燃烧装置内的湍流预混反应流进行了数值模拟。时平均控制方程组的封闭采用k-ε湍流输运模型和EBU-Arhenius湍流反应模型。模拟结果给出了突扩燃烧室内湍流预混反应流的气体时均流场、组分浓度场与温度场的分布。通过数值模拟结果与实验的比较对EBU-Arhenius模型进行了讨论与评价。     

Reaction zone structures and mixing characteristics of partially premixed swirling CH4/air flames in a gas turbine model combustor

The mixing, reaction progress, and flame front structures of partially premixed flames have been investigated in a gas turbine model combustor using different laser techniques comprising laser Doppler velocimetry for the characterization of the flow field, Raman scattering for simultaneous multi-species and temperature measurements, and planar laser-induced fluorescence of CH for the visualization of the reaction zones. Swirling CH₄/air flames with Re numbers between 7500 and 60,000 have been studied to identify the influence of the turbulent flow field on the thermochemical state of the flames and the structures of the CH layers. Turbulence intensities and length scales, as well as the classification of these flames in regime diagrams of turbulent combustion, are addressed. The results indicate that the flames exhibit more characteristics of a diffusion flame (with connected flame zones) than of a uniformly premixed flame.     

Effect of inlet flow turbulence on the combustion instability in a premixed backward-facing step combustor

This paper reports the effect of inlet flow turbulence intensity on the combustion instability characteristics in a backward facing step combustor. The inlet turbulence intensity is varied by a turbulence generator. Unsteady pressure measurements and OH* chemiluminescence images are recorded over a wide range of operating conditions at different inlet turbulence intensities. The study shows an early onset of instability at low turbulence level, i.e., higher turbulence postpones the onset of instability to higher Reynolds number Re and/or higher equivalence ratio Φ. The early onset of instability in the Re and Φ parameter spaces is due to the change in system parameters such as flame speed and size of the recirculation zone downstream of the step at different turbulence levels. Further, the onset is characterized as subcritical bifurcation. At low Re, the hysteresis zone width is small for low turbulence levels and it is large at higher turbulence levels; and at higher Re, the hysteresis width remains constant at all turbulence levels. Investigation of instability characteristics reveals that there are momentary slippages from limit cycle orbit into brief silent regimes in an intermittent manner. The frequency of occurrence of the momentary silent regimes increases with reduction in turbulence, indicating that higher turbulence helps in maintaining the system in a stable limit cycle orbit. High-speed chemiluminescence imaging reveals the necessity of the vortex rollup in the recirculation zone to grow up to the top wall by dilatation from the heat release for the onset of instability. Considerations of the effect of turbulence on both the flame speed and the recirculation zone size together explain all the observed bifurcation trends. These results suggest that inlet flow turbulence should not just be considered as background noise. The turbulence effects on both the flame and flow should be considered in predicting the instability characteristics.     

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.     

Modelling of lifted turbulent diffusion flames in a channel mixing layer by the flame hole dynamics

The partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics as an effort to develop a prediction model for the turbulent flame lift off. The essence of the flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of the quenching holes initially created by the local quenching events. The numerical simulation for the flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for a constant-density fuel–air channel mixing layer to obtain the background turbulent flow and mixing fields, from which a time series of two-dimensional scalar-dissipation-rate array is extracted. Subsequently, a Lagrangian simulation of the flame hole random walk mapping, projected to the scalar dissipation rate array, yields a temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. In particular, the probability of encountering the reacting state, while conditioned with the instantaneous scalar dissipation rate, is examined to reveal that the conditional probability has a sharp transition across the crossover scalar dissipation rate, at which the flame edge changes its direction of propagation. This statistical characteristic implies that the flame edge propagation instead of the local quenching event is the main mechanism controlling the partial quenching events in turbulent flames. In addition, the conditional probability can be approximated by a heavyside function across the crossover scalar dissipation rate.     

Effects of Lewis number,density ratio and gravity on burning velocity and conditional statistics in stagnating turbulent premixed flames

DNS is performed to analyse the effects of Lewis number (Le), density ratio and gravity in stagnating turbulent premixed flames. The results show good agreement with those of Lee and Huh (Combustion and Flame, Vol. 159, 2012, pp. 1576–1591) with respect to the turbulent burning velocity, S_T, in terms of turbulent diffusivity, flamelet thickness, mean curvature and displacement speed at the leading edge. In all four stagnating flames studied, a mean tangential strain rate resulting in a mean flamelet thickness smaller than the unstretched laminar flame thickness leads to an increase in S_T. A flame cusp of positive curvature involves a superadiabatic burned gas temperature due to diffusive–thermal instability for an Le less than unity. Wrinkling tends to be suppressed at a larger density ratio, not enhanced by hydrodynamic instability, in the stagnating flow configuration. Turbulence is produced, resulting in highly anisotropic turbulence with heavier unburned gas accelerating through a flame brush by Rayleigh–Taylor instability. Results are also provided on brush thickness, flame surface density and conditional velocities in burned and unburned gas and on flame surfaces to represent the internal brush structures for all four test flames.     

Computational investigations of the coupling between transient flame dynamics and thermo-acoustic instability in a self-excited resonance combustor

Combustion instability due to thermo-acoustic interactions is a critical combustion problem that requires a thorough understanding because of its adverse impact on stable and reliable operation of combustors in high-speed propulsion devices like gas turbines and rockets. This work conducts computational investigations of the coupling between the transient flame dynamics such as the ignition delay and local extinction and the thermo-acoustic instability developed in a self-excited resonance combustor to gain deep insights into the mechanisms of thermo-acoustic instability. A 2D modelling framework that employs different flamelet models (the steady flamelet model and the flamelet/progress variable approach) is developed to enable the examination of the effect of the transient flame dynamics caused by the strong coupling of the turbulent mixing and finite-rate chemical kinetics on the occurrence of thermo-acoustic instability. The models are validated by using the available experimental data for the pressure signal. Parametric studies are performed to examine the effect of the occurrence of the transient flame dynamics, the effect of artificial amplification of the Damköhler number, and the effect of neglecting mixture fraction fluctuations on the predictions of the thermo-acoustic instability. The parametric studies reveal that the occurrence of transient flame dynamics has a strong influence on the onset of the thermo-acoustic instability. Further analysis is then conducted to localise the effect of a particular flame dynamic event, the ignition delay, on the thermo-acoustic instability. The reverse effect of the occurrence of the thermo-acoustic instability on the transient flame dynamics in the combustor is also investigated by examining the temporal evolution of the local flame events in conjunction with the pressure wave propagation. The above observed two-way coupling between the transient flame dynamics (the ignition delay) and the thermo-acoustic instability provides a plausible mechanism of the self-excited and sustained thermo-acoustic instability observed in the combustor despite the fact that the results are obtained from 2D simulations. The same analysis is expected to be extensible to fully 3D simulations.     

The differential diffusion effect of the intermediate species on the stability of premixed flames propagating in microchannels

The propagation of premixed flames in adiabatic and non-catalytic planar microchannels subject to an assisted or opposed Poiseuille flow is considered. The diffusive–thermal model and the well-known two-step chain-branching kinetics are used in order to investigate the role of the differential diffusion of the intermediate species on the spatial and temporal flame stability. This numerical study successfully compares steady-state and time-dependent computations to the linear stability analysis of the problem. Results show that for fuel Lewis numbers less than unity, Le_F < 1, and at sufficiently large values of the opposed Poiseuille flow rate, symmetry-breaking bifurcation arises. It is seen that small values of the radical Lewis number, Le_Z, stabilise the flame to symmetric shape solutions, but result in earlier flashback. For very lean flames, the effect of the radical on the flame stabilisation becomes less important due to the small radical concentration typically found in the reaction zone. Cellular flame structures were also identified in this regime. For Le_F > 1, flames propagating in adiabatic channels suffer from oscillatory instabilities. The Poiseuille flow stabilises the flame and the effect of Le_Z is opposite to that found for Le_F < 1. Small values of Le_Z further destabilise the flame to oscillating or pulsating instabilities.     

Direct evaluation of the subgrid scale scalar flux in turbulent premixed flames with conditioned dual-plane stereo PIV

With the dual-plane stereo PIV technique the instantaneous three-dimensional resolved rate-of-strain tensor is directly measured in turbulent premixed flames. Simultaneously, also the instantaneous subgrid scale (SGS) scalar flux is measured with fine resolution, where for the latter term the conditioned particle image velocimetry (CPIV) technique is applied. The subgrid resolution reaches 118 μm, allowing a 9 × 9 resolution of a subgrid filter with width Δ = 1 mm. This combined measurement approach allows the a-priori comparison of models for the SGS scalar flux term with direct measurements which is important for large eddy simulation methods in turbulent premixed flames. Two different flame conditions of a premixed V-shaped turbulent flame are investigated where the turbulence intensity is varied by a factor of nearly three. The instantaneous radial and axial SGS fluxes are compared with the following three models: gradient model with Smagorinsky approach for the turbulent viscosity, Clark model, and extended gradient model with an anisotropy term. None of these models shows a good correlation with the directly measured flux. The anisotropy term alone (being nearly similar to the Clark model) shows, however, a right trend behaviour. An analysis of the data indicates a significant dependency of the experimentally determined SGS flux on the Favre averaged reaction progress (spatially averaged over the SGS area). A relatively simple closure for the SGS flux, which describes the dilatation due to the gasdynamic expansion, and which is a function proportional to , shows a rather good correlation with direct measurement for some of the components. A successful SGS scalar flux model for premixed turbulent flames most likely needs to include at least two different effects.     

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.     

Numerical computation of an immiscible two-phase flow in comparison with experiment

The flow of fluids with interfaces or free surfaces is of great interest in both basic and applied research. Despite enormous efforts to develop numerical methods for solving such flow problems, no undisputed standards have emerged yet. We compare here three different approaches for the problem of a liquid film falling down along a vertical wall under the action of gravity. Due to an instability, waves form on the film surface and a nontrivial flow pattern emerges. This provides a test case for two-phase flow simulation methods under conditions of practical relevance. The methods compared are a moving boundary method and a volume of fluid method implemented in commercial software packages and a new phase field method. A set of experimental data is used as a reference for the comparison. All three methods correctly reproduce the main features of the experiment; however, none are free of quantitative deviations from the data. The text was submitted by the authors in English.