On nonlinear model equations for the response of premixed flames to acoustic like accelerations

The response of a dynamical flame model to imposed acoustic accelerations is studied analytically and numerically. Through linear stability analyses, two analytical approximations for the primary and the parametric stability boundaries are found. The approximation for the primary instability boundary is accurate for any periodic accelerations, in the limit of large acoustic frequencies. The critical acoustic amplitude u _a for Landau–Darrieus instability suppression is identified and found to depend only on the density contrast and the shape of the periodic acoustic stimuli. The proposed model evolution equation is next integrated numerically with various imposed acoustic accelerations; the primary and parametric flame responses are identified. It is shown analytically and numerically that in the presence of a fully developed, yet weakened by acoustics, Landau–Darrieus (or primary) instability the wrinkle amplitude and the mean flame speed oscillate at the same frequency as the acoustic stimuli; the threshold for suppression of primary instability by acoustic forcing is determined exactly. Increasing the acoustic amplitude allows the flame to respond parametrically to the acoustics. This response is characterised by troughs and crests interchanging their roles while the mean flame speed again oscillates with the same frequency as the acoustic stimuli and at twice that of wrinkle amplitude oscillations.     

Impact of pressure fluctuations on the dynamics of laminar premixed flames

Thermo-acoustic instabilities are problematic in the design of continuous-combustion propulsion systems such as gas turbine engines, rocket motors, jet engine afterburners, and ramjets. Conceptually, the coupling between acoustics and flame dynamics can be divided into two categories: flame area fluctuations and changes in the local flame speed. The latter can be caused by the thermodynamic fluctuations that accompany an acoustic wave. This coupling is the focus of the present work. In this paper, we are concerned with the dynamics of laminar premixed flames involving large hydrocarbon species. Through high-fidelity numerical simulations, we investigate the flame response for a wide range of fuels and acoustic frequencies. The combustion of hydrogen and methane is considered for verification purposes and as baseline cases for comparison with two large hydrocarbon fuels, n-heptane and n-dodecane. We extract the phase and gain of the unsteady heat release response, which are directly related to the Rayleigh criterion and thus the stability of the system. For all fuels, we observe a local peak in the heat release gain. At high frequencies, we find that the fluctuations of the different species mass fractions decrease with the inverse of the acoustic frequency, leading to chemistry being “frozen” in the high-frequency limit. This allows us to predict the flame behavior directly from the steady-state solution.     

Memory effects of local flame dynamics in turbulent premixed flames

A premixed and thermo-diffusively unstable turbulent hydrogen-air flame-in-a-box case is simulated in conjunction with the flame particle tracking (FPT) method. The flame is located in the flamelet regime. The focus lies on the assessment of memory effects in local flame dynamics. By tracking flame particles on an iso-surface of the flame during flame-turbulence interaction, the time history of flame speed and flame stretch can be recorded for each point on the flame iso-surface in a Lagrangian reference frame. The results reveal a time delay between the local flame speed and flame stretch signal, showing that previous values of flame stretch affect currently observed values of flame speed. Furthermore, by choosing flame particles whose trajectories are dominated by single frequencies, the time delay can be quantified. While plotting instantaneous values of flame speed and flame stretch results in a large scattering for turbulent flames, a quasi-linear correlation can be achieved by shifting the time signal of flame stretch according to the time delay. The time delay itself depends on the local flow time scale, which is expressed as a local Damköhler number. There is, however, an important difference between consumption and displacement speed. While most analyses in the literature are limited to the flame displacement speed, the flame consumption speed is evaluated for each flame particle in this work as well, which shows a strong correlation with the local equivalence ratio even at unsteady conditions. As the flame particles move toward regions with more negative flame stretch, the consumption speed decreases as the flame locally extinguishes. At the same time, the diffusive component of the displacement speed increases, as the tangential component of the diffusive flux increases in regions with strong negative flame curvature.     

Ignition and transient dynamics of sub-limit premixed flames in microchannels

We examine, via two-dimensional numerical simulation of a model system, some unsteady transient ignition scenarios and sustained oscillatory combustion modes that can occur in a single-pass, conductive channel, premixed microburner. These issues are relevant to the problem of ignition, evolution to stable combustion and the operational modes of microcombustors. First, we describe an unsteady ignition sequence that may occur when a single-pass microburner with initially cold walls has its exit walls heated and maintained at a fixed temperature. In particular, we demonstrate that as the heat from the exit walls propagates down the microburner walls, a reaction wave is driven rapidly down the channel towards the inlet via a sequence of oscillatory ignition and quenching transients. This scenario has been observed experimentally during the ignition of a single-pass microburner. Secondly, we show how an initial axial wall temperature gradient can lead to a variety of sustained combustion modes within the channel, including stable stationary flames, regimes of periodic motion involving quenching and re-ignition, regimes of regular oscillatory combustion, and regimes consisting of a combination of re-ignition/quenching events and regular oscillatory motions, all of which have been observed experimentally.     

Explosive dynamics of bluff-body-stabilized lean premixed hydrogen flames at blow-off

Two-dimensional direct numerical simulation (DNS) databases of bluff-body-stabilized lean hydrogen flames representative of complicated reactive–diffusive system are analysed using the combined approach of computational singular perturbation (CSP) and tangential stretching rate (TSR) to investigate chemical characteristics in blow-off dynamics. To assess the diagnostic approaches in flame and blow-off dynamics, Damköhler number and TSR variables are applied and compared. Four cases are considered in this study showing different flame dynamics such as the steadily stable mode, local extinction by asymmetric vortex shedding, convective blow-off and lean blow-out. DNS data points in positive explosive eigenvalue conditions were subdivided into four different combinations in TSR and extended TSR space and categorized in four distinct characteristic regions, such as kinetically explosive or dissipative and transport-enhanced or dissipative dynamics. The TSR analysis clearly captures the local extinction point in the complicated vortex shedding and allows an improved understanding of the distinct chemistry-transport interactions occurring in convective blow-off and lean blow-out events.     

Minimum ignition energy and propagation dynamics of laminar premixed cool flames

Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.     

Pressure gradient tailoring effects on vorticity dynamics in the near-wake of bluff-body premixed flames

We investigate the role of mean streamwise pressure gradients in the development of a bluff-body-stabilized premixed flame in the near-wake of the bluff body. To this end, a triangular prism flame holder is situated in three different channel geometries: a nominal case with straight walls, a nozzle with a stronger mean pressure gradient, and a diffuser with a comparatively weaker mean pressure gradient. All geometries are implemented using embedded boundaries, and adaptive mesh refinement is used to locally resolve all relevant thermal (i.e., flame) and fluid-mechanical (i.e., vorticity) scales. A premixed propane flame, modeled using a 66-step skeletal mechanism, interacts with vorticity in the boundary layer of the triangular bluff body in the presence of each mean pressure gradient. Analysis of flame-related enstrophy budget terms reveals key differences in the behavior of baroclinic torque between cases, the specifics of which are tied to larger variations in the mean flow structure, recirculation zone structure, and confinement effects. Our results show that the baroclinic torque changes significantly among the configurations, with the nozzle exhibiting the largest baroclinic torque production. However, these differences are shown to be only a secondary consequence of the background pressure gradient, with the primary consequence being the change in the recirculation zone length resulting from the different channel configurations. These results are relevant for flame stabilization with bluff bodies, where clear understanding of the sensitivities to global mean pressure gradient is important to engineering design.     

Nonperturbative renormalization group approach to turbulence

We suggest a new, renormalization group (RG) based, nonperturbative method for treating the intermittency problem of fully developed turbulence which also includes the effects of a finite boundary of the turbulent flow. The key idea is not to try to construct an elimination procedure based on some assumed statistical distribution, but to make an ansatz for possible RG transformations and to pose constraints upon those, which guarantee the invariance of the nonlinear term in the Navier-Stokes equation, the invariance of the energy dissipation, and other basic properties of the velocity field. The role of length scales is taken to be inverse to that in the theory of critical phenomena; thus possible intermittency corrections are connected with the outer length scale. Depending on the specific type of flow, we find different sets of admissible transformations with distinct scaling behaviour: for the often considered infinite, isotropic, and homogeneous system K41 scaling is enforced, but for the more realistic plane Couette geometry no restrictions on intermittency exponents were obtained so far. Received: 28 December 1997 / Accepted: 6 August 1998     

Effects of gas compressibility on the dynamics of premixed flames in long narrow adiabatic channels

We examine the dynamics of premixed flames in long, narrow, adiabatic channels focusing, in particular, on the effects of gas compressibility on the propagation. Recognising the importance of the boundary conditions, we examine and compare three cases: flame propagation in channels open at both ends, where the pressure must adjust to the ambient pressure at both ends and the expanding gas is allowed to leave the channel freely, and flame propagation in channels that remain closed at one of the two ends, where the burned/unburned gas remains trapped between the flame and one of the two walls. Earlier studies have shown that a flame accelerates when travelling down a narrow channel as a result of the combined effects of wall friction and thermal expansion. In the present work we show that compressibility effects enhance the transition to fast accelerating flames in channels open at both ends and in channels closed at the ignition end. In both situations, the accelerating flames could reach values that, depending on the effective Mach number, are as large as fifty times the laminar flame speed. In contrast, when the channel is closed at the far end, the acceleration is limited and the propagation speed is damped as the flame approaches the far boundary. Moreover, we show that, in channels closed at their ignition end, the flame in sufficiently long channels evolves into a steadily propagating compression-driven flame. The propagation speed of these flames depends exponentially on the constant-volume equilibrium temperature, which is higher than the (constant pressure) adiabatic flame temperature, and is therefore larger than for ordinary isobaric flames. Fast propagating compression waves cannot emerge in channels that remain open at their ignition end because of the reduced pressure forced by the open boundary.     

Phenomenological model of chain-branching premixed flames

In this work, we introduce a global kinetic model that includes fuel, oxygen, products and two radical species involved in the reversible chain-branching, chain-propagation and chain-termination reactions. The model naturally extends the Zeldovich–Liñán and Zeldovich–Barenblatt–Dold models and can be used to describe both premixed and diffusion flames. Here it is applied to the problem of the deflagration wave propagation in the hydrocarbon fuel/air mixture with arbitrary equivalence ratio under the simplifying thermal-diffusive approximation. The conservation equations are solved numerically in order to obtain the velocity and structure of the combustion wave. It is demonstrated that the peak values of the adiabatic flame temperature and deflagration velocity are shifted towards the rich mixture composition if the reverse reactions of product decomposition are taken into account. The dependence of the flame speed and temperature on parameters of the system is analysed. The prospects of further investigation are discussed.     

Structure and propagation of two-dimensional,partially premixed,laminar flames in diesel engine conditions

We investigate the influence of inflow velocity (V_in) and scalar dissipation rate (χ) on the flame structure and stabilisation mechanism of steady, laminar partially premixed n-dodecane edge flames stabilised on a convective mixing layer. Numerical simulations were performed for three different χ profiles and several V_in (V_in = 0.2 to 2.5m/s). The ambient thermochemical conditions were the same as the Engine Combustion Network’s (ECN) Spray A flame, which in turn represents conditions in a typical heavy duty diesel engine. The results of a combustion mode analysis of the simulations indicate that the flame structure and stabilisation mechanism depend on V_in and χ. For low V_in the flame is attached. Increasing V_in causes the high-temperature chemistry (HTC) flame to lift-off, while the low-temperature chemistry (LTC) flame is still attached. A unique speed S_R associated with this transition is defined as the velocity at which the lifted height has the maximum sensitivity to changes in V_in. This transition velocity is negatively correlated with χ. Near $V_{in}=S_R$ a tetrabrachial flame structure is observed consisting of a triple flame, stabilised by flame propagation into the products of an upstream LTC branch. Further increasing the inlet velocity changes the flame structure to a pentabrachial one, where an additional HTC ignition branch is observed upstream of the triple flame and ignition begins to contribute to the flame stabilisation. At large V_in, the LTC is eventually lifted, and the speed at which this transition occurs is insensitive to χ. Further increasing V_in increases the contribution of ignition to flame stabilisation until the flame is completely ignition stabilised. Flow divergence caused by the LTC branch reduces the χ at the HTC branches making the HTC more resilient to χ. The results are discussed in the context of identification of possible stabilisation modes in turbulent flames.     

The response of turbulent premixed flames to normal acoustic excitation

To model the thermo-acoustic excitation of flames in practical combustion systems, it is necessary to know how a turbulent flame front responds to an incident acoustic wave. This will depend partly on the way in which the burning velocity responds to the wave. In this investigation, the response of CH₄/air and CH₄/H₂/air mixtures has been observed in a novel flame stabilisation configuration, in which the premixture of fuel and air is made to decelerate under controlled conditions in a wide-angle diffuser. Control is provided by an annular wall-jet of air and by turbulence generators at the inlet. Ignition from the outlet of the diffuser allows an approximately flat flame to propagate downwards and stabilise at a height that depends on the turbulent burning velocity. When the flow is excited acoustically, the ensemble-averaged height oscillates. The fluctuations in flow velocity and flame height are monitored by phase-locked particle image velocimetry and OH-planar laser induced fluorescence, respectively. The flame stabilised against a lower incident velocity as the acoustic amplitude increased. In addition, at the lowest frequency of 52 Hz, the fluctuations in turbulent burning velocity (as represented by the displacement speed) were out-of-phase with the acoustic velocity. Thus, the rate of displacement of the flame front relative to the flow slowed as the flow accelerated, and so the flame movement was bigger than it would have been if the burning velocity had not responded to the acoustic fluctuation. With an increase in frequency to 119 Hz, the relative flame movement became even larger, although the phase-difference was reduced, so the effect on burning velocity was less dramatic. The addition of hydrogen to the methane, so as to maintain the laminar burning velocity at a lower equivalence ratio, suppressed the response at low amplitude, but at a higher amplitude, the effect was reversed.     

Nonlinear kinematic response of premixed flames to harmonic velocity disturbances

This paper analyzes the nonlinear dynamics of premixed flames responding to harmonic velocity disturbances. These nonlinear dynamics were studied by solving a constant flame speed front tracking equation for the flame’s response to harmonically oscillating velocity disturbances. The solution to these equations is used to quantify the transfer function relating the ratio of the normalized flame area to velocity fluctuations, G = (A′/A_o)/(u′/u_o), upon the amplitude of velocity oscillations, ε = u′/u_o. Due to nonlinearities, the amplitude of this transfer function relative to its linear value decreases with increasing amplitude of velocity oscillation, u′/u_o. In contrast, the transfer function phase exhibits almost no amplitude dependence. The velocity amplitude where transfer function nonlinearities become significant depends strongly upon three parameters: a Strouhal number, St = ωL_f/u_o (where L_f is the flame length), the ratio of the flame length to width, β = L_f/R, and the flame shape in the absence of perturbations (i.e., conical, inverted wedge, etc.). In the linear case, the transfer function, G, depends only upon an algebraic combination of the first two parameters, given by St₂ = St (1 + β²)/β². In general, however, G exhibits a distinct dependence upon both parameters St and β. In particular, we show that the nonlinear response of G is an intrinsically dynamic phenomenon; i.e., its quasi-steady response (St 1) is purely linear. As such, nonlinearity is enhanced with increasing Strouhal numbers. In contrast, nonlinearity is suppressed at large β values; as such, the response of a long flame remains quite similar to its linear value, even at large ε values where the flame front exhibits substantial corrugation and cusping. Finally, we show that the response of conical flames remains much more linear at comparable disturbance amplitudes than for “V” or wedge-shaped flames. These predictions are shown to be consistent with available experimental data.     

PDF calculations of piloted premixed jet flames

A series of piloted premixed jet flames with strong finite-rate chemistry effects is studied using the joint velocity-turbulence frequency-composition PDF method. The numerical accuracy of the calculations is demonstrated, and the calculations are compared to experimental data. It is found that all calculations show good agreement with the measurements of mean and rms mixture fraction fields, while the reaction progress is overpredicted to varying degrees depending on the jet velocity. In the calculations of the flame with the lowest jet velocity, the species and temperature show reasonable agreement with the measurements, with the exception of a small region near the centerline where products and temperature are overpredicted and fuel and oxidizer are underpredicted. In the calculations of the flame with the highest jet velocity, however, the overprediction of products and temperature and underprediction of fuel and oxidizer is far more severe. An extensive set of sensitivity studies on inlet boundary conditions, turbulence model constants, mixing models and constants, radiation treatment, and chemical mechanisms is conducted to show that any parameter variation offers little improvement from the base case. To shed light on these discrepancies, diagnostic calculations are performed in which the chemical reactions are artificially slowed. These diagnostic calculations serve to validate the experimental data and to quantify the amount by which the base case calculations overpredict reaction progress. Improved calculations of this flame are achieved only through artificially slowing down the chemical reaction by a factor of about 10. The mixing model behavior in this combustion regime is identified as a likely cause for the observed discrepancy in reaction progress.     

New approach to nonlinear dynamics of fullerenes and fullerites

A new type of nonlinear (anharmonic) excitations—bushes of vibrational modes—in physical systems with point or space symmetry is discussed. All infrared-active and Raman-active bushes for C₆₀ fullerene are found by means of special group-theoretical methods.