Bifurcation characteristics and flame dynamics of a ducted non-premixed flame with finite rate chemistry

The influence of system parameters such as the flame location, Peclet number and Damköhler number on the bifurcation characteristics and flame dynamics of a ducted non-premixed flame with finite rate chemistry is presented in this paper. In the bifurcation plot with flame location as the bifurcation parameter, subcritical Hopf bifurcation is found for lower values of flame location and supercritical Hopf bifurcation for higher values of flame location, for all the Damköhler numbers used in this study. The flame shapes are captured at eight different phases of a cycle of time series data of acoustic velocity at both the fold and Hopf points for bifurcation with flame location as the parameter. We find that the range of flame height variations at the Hopf point is more than the range of flame height variations obtained at the fold point. We also find that the flame oscillates in the same phase as pressure fluctuation but in a phase different from both velocity and heat release rate fluctuations in the region of hysteresis for bifurcation with flame location. The non-dimensional hysteresis width is plotted as a function of Damköhler number for variation of flame location in the subcritical region. An inverse power law relation is found between the non-dimensional hysteresis width and the Damköhler number. The bifurcation plot with Peclet number as parameter shows a subcritical Hopf bifurcation.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

An elementary model for the validation of flamelet approximations in non-premixed turbulent combustion

The fundamental soundness of three flamelet models for non-premixed turbulent combustion is examined on the basis of their performance in an idealized model problem that merges ideas from the laminar asymptotic theory for non-premixed flames and rigorous homogenization theory for the diffusion of a passive scalar. The overall flame configuration is stabilized by a mean gradient in the passive scalar: large Damköhler number asymptotics results are available for the laminar case to quantify the finite-rate effects that cause the flame to depart from its equilibrium state; the same results can also be used to incorporate higher-order corrections in the approximation of the reactive variables in terms of the passive scalar. The use of such flamelet approximations has been extended well beyond the laminar regime as they lie at the core of practical strategies to simulate non-premixed flames in the turbulent regime: the flamelet representation avoids the problem of turbulence closure for the reactive variables by replacing it by the presumably much simpler closure problem for a passive scalar. It is precisely the validity of this substitution outside the laminar regime that is addressed here in the idealized context of a class of small-scale periodic flows for which extensive rigorous results are available for the passive scalar statistics. Results for this simplified problem are reported here for significant wide ranges of Peclet and Damköhler numbers. Asymptotic convergence is observed in terms of the Damköhler number, with a convergence rate that is found to match the laminar predictions and appears relatively insensitive to the Peclet number. The passive scalar dissipation plays a key role in achieving higher-order corrections for the finite-rate case: replacing its pointwise value by an averaged value is convenient practically and can be rigorously motivated for the class of flows studied here, but while it does achieve an overall improvement over the lower-order equilibrium model, the simplification compromises the higher asymptotic convergence observed with the original finite-rate flamelet model with exact local dissipation.(Some figures in this article are in colour only in the electronic version; see www.iop.org)     

Transient convective burning of interactive fuel droplets in single-layer arrays

The transient convective burning of n-octane droplets interacting within single-layer arrays in a hot gas flow perpendicular to the layer is studied numerically, with considerations of droplet surface regression, deceleration due to the drag of the droplets, internal liquid motion, variable properties, non-uniform liquid temperature and surface tension. Infinite periodic arrays, semi-infinite periodic arrays with one row of droplets (linear array) or two rows of droplets, and finite arrays with nine droplets with centers in a plane are investigated. All arrays are aligned orthogonal to the free stream direction. This paper compares the behavior of semi-infinite periodic arrays and finite arrays with the behavior of previously studied infinite periodic arrays. Furthermore, it identifies the critical values of the initial Damköhler number for bifurcations in flame behavior at various initial droplet spacing for all these arrays. The initial flame shape is either an envelope flame or a wake flame as determined by the initial Damköhler number, the array configuration and the initial droplet spacing. The critical initial Damköhler number separating initial wake flames from initial envelope flames decreases with increasing interaction amongst droplets at intermediate droplet spacing (when the number of rows in the array increases or the initial droplet spacing decreases for a specific number of rows in the array). In the transient process, an initial wake flame has a tendency to develop from a wake flame to an envelope flame, with the moment of wake-to-envelope transition advanced for the increasing interaction amongst droplets at intermediate droplet spacing. For the array with nine droplets with centers in a plane, the droplets at different types of positions have different critical initial Damköhler number and different wake-to-envelope transition time for initial wake flame.     

An experimental study of partially premixed flame sound

The occurrence of oscillating combustion and combustion instability has led to resurgence of interest in the causes, mechanisms, suppression, and control of combustion noise. Noise generated by enclosed flames is of greater practical interest but is more complicated than that by open flames, which itself is not clearly understood. Studies have shown that different modes of combustion, premixed and non-premixed, differ in their sound generation characteristics. However, there is lack of understanding of the region bridging these two combustion modes. This study investigates sound generation by partially premixed flames. Starting from a non-premixed flame, air was gradually added to achieve partial premixing while maintaining the fuel flow rate constant. Methane, ethylene, and ethane partially premixed flames were studied with hydrogen added for flame stabilization. The sound pressure generated by methane partially premixed flames scales with M⁵ compared to M³ for turbulent non-premixed methane flames. Also, the sound pressure generated by partially premixed flames of ethane and ethylene scales as M^4.5. With progressive partial premixing, spectra level increases at all frequencies with a greater increase in the high-frequency region compared to the low-frequency region; flames develop a peak and later a constant level plateau in the low frequency region. The partially premixed flames of methane, ethylene, and ethane generate a similar SPL as a function of equivalence ratio when the fuel volume flow rate is matched. However, when fuel mass flow rate is matched, the ethane and ethylene flames produce a similar SPL, which is lower than that produced by the methane flame.     

Structure and dynamics of edge flames in the near wake of unequal merging shear flows

We examine in this study the structure and dynamic properties of an edge flame formed in the near-wake of two initially separated shear flows, one containing fuel and the other oxidiser. A comprehensive study is carried out within the diffusive-thermal framework where the flow field, computed a-priori, is used for the determination of the combustion field. Our focus is on the effects of three controlling parameters: the Damköhler number controlling the overall flow rate, the oxidiser-to-fuel strain rate ratio of the supply streams that determines the extent of oxidiser entrainment towards the mixing zone, and the Lewis number, assumed equal for the fuel and oxidiser, that depends on the mixture composition. Response curves, representing the edge flame standoff distance as a function of Damköhler number, exhibit two distinct shapes: C-shaped and U-shaped curves characterising the response of low and high Lewis number flames, respectively. Stability considerations show that the upper solution branch of the C-shaped response curve is unstable and hence corresponds to physically unrealistic states, but due to heat conduction toward the cold plate the lower solution branch is always stable. The states forming this solution branch correspond to flame attachment, where the edge flame remains practically attached to the tip of the plate until it is blown off by the flow when the velocity exceeds a critical value. The U-shaped response, on the other hand, consists of equilibrium states that are globally stable. Thus, high Lewis number flames can be always stabilised near the splitter plate, with the edge held stationary or undergoing a back and forth motion, or lifted and stabilised downstream by the flow. Insight into the distinct stabilisation characteristics, exhibited by the different Lewis number cases, is given by examining the relationship between the local flow velocity and the edge propagation speed.     

Numerical investigation of flame–vortex interactions in laminar cross-flow non-premixed flames in the presence of bluff bodies

Flame stabilisation in a combustor having vortices generated by flame holding devices constitutes an interesting fundamental problem. The presence of vortices in many practical combustors ranging from industrial burners to high speed propulsion systems induces vortex–flame interactions and complex stabilisation conditions. The scenario becomes more complex if the flame sustains after separating itself from the flame holder. In a recent study [P.K. Shijin, S.S. Sundaram, V. Raghavan, and V. Babu, Numerical investigation of laminar cross-flow non-premixed flames in the presence of a bluff-body, Combust. Theory Model. 18, 2014, pp. 692–710], the authors reported details of the regimes of flame stabilisation of non-premixed laminar flames established in a cross-flow combustor in the presence of a square cylinder. In that, the separated flame has been shown to be three dimensional and highly unsteady. Such separated flames are investigated further in the present study. Flame–vortex interactions in separated methane–air cross flow flames established behind three bluff bodies, namely a square cylinder, an isosceles triangular cylinder and a half V-gutter, have been analysed in detail. The mixing process in the reactive flow has been explained using streamlines of species velocities of CH₄ and O₂. The time histories of z-vorticity, net heat release rate and temperature are analysed to reveal the close relationship between z-vorticity and net heat release rate spectra. Two distinct fluctuating layers are visible in the proper orthogonal decomposition and discrete Fourier transform of OH mass fraction data. The upper fluctuating layer observed in the OH field correlates well with that of temperature. A detailed investigation of the characteristics of OH transport has also been carried out to show the interactions between factors affecting fluid dynamics and chemical kinetics that cause multiple fluctuating layers in the OH.     

The role of unequal diffusivities in ignition and extinction fronts in strained mixing layers

We have studied flame propagation in a strained mixing layer formed between a fuel stream and an oxidizer stream, which can have different initial temperatures. Allowing the Lewis numbers to deviate from unity, the problem is first formulated within the framework of a thermo-diffusive model and a single irreversible reaction. A compact formulation is then derived in the limit of large activation energy, and solved analytically for high values of the Damköhler number. Simple expressions describing the flame shape and its propagation velocity are obtained. In particular, it is found that the Lewis numbers affect the propagation of the triple flame in a way similar to that obtained in the studies of stretched premixed flames. For example, the flame curvature determined by the transverse enthalpy gradients in the frozen mixing layer leads to flame-front velocities which grow with decreasing values of the Lewis numbers.

The analytical results are complemented by a numerical study which focuses on preferential-diffusion effects on triple flames. The results cover, for different values of the fuel Lewis number, a wide range of values of the Damköhler number leading to propagation speeds which vary from positive values down to large negative values     

Critical conditions at liftoff limit of a laminar n-butane-air jet diffusion flame

The stability mechanism of laminar coflow jet diffusion flames in normal gravity has been studied computationally and experimentally. N-butane, the heaviest alkane in a gaseous state at ambient temperature and pressure, is used as the fuel since the reaction mechanism is similar to that of higher (liquid) hydrocarbons. The critical mean n-butane jet and coflowing air velocities at flame stability limits are measured using a small fuel tube burner (0.8 mm inner diameter). The time-dependent, axisymmetric numerical code with a detailed reaction mechanism (58 species and 540 reactions), molecular diffusive transport, and a radiation model, reveals a flame structure. A fuel-lean peak reactivity spot (i.e., reaction kernel), possessing the hybrid nature of diffusion-premixed flame structure at a constant temperature of ≈1560 K, is formed at the flame base and controls the flame stability. In a near-quiescent environment, the flame base resides below the fuel tube exit plane and thereby premixing is limited. As the coflowing air velocity is increased incrementally under a fixed fuel jet velocity, the flame base moves slightly above (≈1 mm) the burner exit and vigorous premixed combustion becomes prevailing. The local heat-release rate at the reaction kernel nearly doubles due to the increased convective oxygen flux (i.e., a blowing effect). The local Damköhler number, newly defined as a ratio of the square root of the local heat-release rate and the local velocity, decreases gradually first and drops abruptly at a critical threshold value and the flame base lifts off from the burner rim. The calculated coflow air velocity at liftoff is ≈0.38 m/s at the fuel jet velocity of 2 m/s, which is consistent with an extrapolated measured value of 0.41 m/s. This work has determined the critical Damköhler number at the stability limit quantitatively, for the first time, for laminar jet diffusion flames.     

On the sublimit solution branches of the stripe patterns formed in counterflow diffusion flames by diffusional–thermal instability

The nonlinear dynamics of striped diffusion flames, formed in a two-dimensional counterflow by diffusional–thermal instability with Lewis numbers sufficiently less than unity, is investigated numerically by examining various two-dimensional flame-structure solutions bifurcating from the one-dimensional steady solution. The Lewis numbers for fuel and oxidizer are identically set to be 0.3, and an overall single-step Arrhenius-type chemical reaction with a Zel'dovich number of 7 is employed as the chemistry model. Particular attention is focused on the flame-stripe solution branches in the sub-extinction regime and on the hysteresis encountered during the transition between different solution branches. In the numerical simulations, a nonlinear solution with eight stripes is first realized from the one-dimensional solution at a Damköhler number slightly greater than the extinction Damköhler number. The eight-stripe solution survives Damköhler numbers much smaller than the extinction Damköhler number until successive bifurcations, leading to the doubling of the pattern wavelength, occur at the subsequent forward-transition conditions. At the first forward-transition Damköhler number occurs the transition to a four-stripe solution, which in turn transits to a two-stripe solution at the second forward-transition Damköhler number, a value somewhat smaller than the first. However, further transition from a two-stripe solution to a one-stripe solution is not always possible even if a one-stripe solution can be accessed independently for particular initial conditions. The Damköhler-number ranges and shapes for the two-stripe and one-stripe solutions are found to be virtually identical, implying that each stripe could be an independent structure if the distance between stripes is sufficiently large. By increasing the Damköhler number, backward transitions can be observed. In comparison with the forward-transition Damköhler numbers, the corresponding backward-transition Damköhler numbers are always much greater, thereby indicating significant hysteresis between the stripe patterns of strained diffusion flames.     

A numerical study on extinction behaviour of laminar micro-diffusion flames

We conducted a numerical study on the fluid dynamic, thermal and chemical structures of laminar methane–air micro flames established under quiescent atmospheric conditions. The micro flame is defined as a flame on the order of one millimetre or less established at the exit of a vertically-aligned straight tube. The numerical model consists of convective–diffusive heat and mass transport with a one-step, irreversible, exothermic reaction with selected kinetics constants validated for near-extinction analyses. Calculations conducted under the burner rim temperature 300 K and the adiabatic burner wall showed that there is the minimum burner diameter for the micro flame to exist. The Damköhler number (the ratio of the diffusive transport time to the chemical time) was used to explain why a flame with a height of less than a few hundred microns is not able to exist under the adiabatic burner wall condition. We also conducted scaling analysis to explain the difference in extinction characteristics caused by different burner wall conditions. This study also discussed the difference in governing mechanisms between micro flames and microgravity flames, both of which exhibit similar spherical flame shape.     

Understanding cool flames and warm flames

Low-temperature flames such as cool flames, warm flames, double flames, and auto-ignition assisted flames play a critical role in the performance of advanced engines and fuel design. In this paper, an overview of the recent progresses in understanding low-temperature flames and dynamics as well as their impacts on combustion, advanced engines, and fuel development will be presented. Specifically, at first, a brief review of the history of cool flames is made. Then, the recent experimental studies and computational modeling of the flame structures, dynamics, and burning limits of non-premixed and premixed cool flames, warm flames, and double flames are presented. The flammability limit diagram and the temperature-dependent chain-branching reaction pathways, respectively, for hot, warm, and cool flames at elevated temperature and pressure will be discussed and analyzed. After that, the effect of low temperature auto-ignition of auto-igniting mixtures at high ignition Damköhler numbers at engine conditions on the propagation of cool flames, warm flames, and double flames as well as turbulent flames will be discussed. Finally, a new platform using low temperature flames for the development and validation of chemical kinetic models of alternative fuels will be presented. Discussions of future research of the dynamics and control of low temperature flames under engine conditions will be made.     

Response of non-premixed flames to bulk flow perturbations

This paper describes the dynamics of non-premixed flames responding to bulk velocity fluctuations, and compares the dynamics of the flame sheet position and spatially integrated heat release to that of a premixed flame. The space–time dynamics of the non-premixed flame sheet in the fast chemistry limit is described by the stoichiometric mixture fraction surface, extracted from the solution of the

-equation. This procedure has some analogies to premixed flames, where the premixed flame sheet location is extracted from the G = 0 surface of the solution of the G-equation. A key difference between the premixed and non-premixed flame dynamics, however, is the fact that the non-premixed flame sheet dynamics are a function of the disturbance field everywhere, and not just at the reaction sheet, as in the premixed flame problem. A second key difference is that the non-premixed flame does not propagate and so flame wrinkles are convected downstream at the axial flow velocity, while wrinkles in premixed flames convect downstream at a vector sum of the flame speed and axial velocity. With the exception of the flame wrinkle propagation speed, however, we show that that the solutions for the space–time dynamics of the premixed and non-premixed reaction sheets in high velocity axial flows are quite similar. In contrast, there are important differences in their spatially integrated unsteady heat release dynamics. Premixed flame heat release fluctuations are dominated by area fluctuations, while non-premixed flames are dominated by mass burning rate fluctuations. At low Strouhal numbers, the resultant sensitivity of both flames to flow disturbances is the same, but the non-premixed flame response rolls off slower with frequency. Hence, this analysis suggests that non-premixed flames are more sensitive to flow perturbations than premixed flames at O(1) Strouhal numbers.     

Blowout limits of inclined nonpremixed turbulent jet flames

The blowout behavior of inclined nonpremixed turbulent jet flames is investigated by varying the jet inclined angle in the range of -90° to 90° The critical jet velocity at blow-out limit is quantified experimentally for various nozzle diameters, different fuels and inclined angles. Numerical simulations are performed to emphasize the flow field difference for the positive and negative inclined angles. Physical modeling is conducted to incorporate the effect of the inclined angle on blow-out behavior. Major findings include: (1) The negatively inclined jet flames show more intense yellow luminosity with larger sooting zones than the positively inclined jet flames; (2) The blowout limit decreases appreciably with the jet inclined angle for the negatively inclined flames, while for the positively inclined jet flames, this decrease is relatively small; (3) Physical analysis of the flow development of inclined jets is conducted, indicating the centerline velocity along the jet trajectory decreases faster for the flame with smaller inclined angle. And the decrease rate is relatively larger for the negatively inclined jet flames; (4) Based on the analysis of the flow development as well as the characteristic velocity with the inclined angle variation, a model based on the Damköhler number (Da) accounting for the effect of jet inclined angle is developed to characterize the blowout limits of inclined jet flames. The proposed model successfully correlates the experimental data. The present findings provide new data and a basic scaling law for the blowout limit of nonpremixed inclined turbulent jet flames, revealing the effect of the relative angle between the jet momentum and buoyancy.     

An investigation of flame zones and burning velocities of laminar unconfined methane-oxygen premixed flames

Numerical and experimental investigations of unconfined methane-oxygen laminar premixed flames are presented. In a lab-scale burner, premixed flame experiments have been conducted using pure methane and pure oxygen mixtures having different equivalence ratios. Digital photographs of the flames have been captured and the radial temperature profiles at different axial locations have been measured using a thermocouple. Numerical simulations have been carried out with a C2 chemical mechanism having 25 species and 121 reactions and with an optically thin radiation sub-model. The numerical results are validated against the experimental and numerical results for methane-air premixed flames reported in literature. Further, the numerical results are validated against the results from the present methane-oxygen flame experiments. Visible regions in digital flame photographs have been compared with OH isopleths predicted by the numerical model. Parametric studies have been carried out for a range of equivalence ratios, varying from 0.24 to 1.55. The contours of OH, temperature and mass fractions of product species such as CO, CO₂ and H₂O, are presented and discussed for various cases. By using the net methane consumption rate, an estimate of the laminar flame speed has been obtained as a function of equivalence ratio.     

Conditional velocity statistics in the double scalar mixing layer – A mapping closure approach

In this work we use 3D direct numerical simulations (DNS) to investigate the average velocity conditioned on a conserved scalar in a double scalar mixing layer (DSML). The DSML is a canonical multistream flow designed as a model problem for the extensively studied piloted diffusion flames. The conditional mean velocity appears as an unclosed term in advanced Eulerian models of turbulent non-premixed combustion, like the conditional moment closure and transported probability density function (PDF) methods. Here it accounts for inhomogeneous effects that have been found significant in flames with relatively low Damköhler numbers. Today there are only a few simple models available for the conditional mean velocity and these are discussed with reference to the DNS results. We find that both the linear model of Kutznetzov and the Li and Bilger model are unsuitable for multi stream flows, whereas the gradient diffusion model of Pope shows very close agreement with DNS over the whole range of the DSML. The gradient diffusion model relies on a model for the conserved scalar PDF and here we have used a presumed mapping function PDF, that is known to give an excellent representation of the DNS. A new model for the conditional mean velocity is suggested by arguing that the Gaussian reference field represents the velocity field, a statement that is evidenced by a near perfect agreement with DNS. The model still suffers from an inconsistency with the unconditional flux of conserved scalar variance, though, and a strategy for developing fully consistent models is suggested.     

Turbulent scalar flux transport in head-on quenching of turbulent premixed flames: a direct numerical simulations approach to assess models for Reynolds averaged Navier Stokes simulations

The statistical behaviour and the modelling of turbulent scalar flux transport have been analysed using a direct numerical simulation (DNS) database of head-on quenching of statistically planar turbulent premixed flames by an isothermal wall. A range of different values of Damköhler, Karlovitz numbers and Lewis numbers has been considered for this analysis. The magnitudes of the turbulent transport and mean velocity gradient terms in the turbulent scalar flux transport equation remain small in comparison to the pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation when the flame is away from the wall but the magnitudes of all these terms diminish and assume comparable values during flame quenching before vanishing altogether. It has been found that the existing models for the turbulent transport, pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation do not adequately address the respective behaviours extracted from DNS data in the near-wall region during flame quenching. Existing models for transport equation-based closures of turbulent scalar flux have been modified in such a manner that these models provide satisfactory prediction both near to and away from the wall.     

Multiple speeds of flame edge propagation for Lewis numbers above one

Edges of diffusion flames in a counterflow burner are examined numerically for Lewis greater than unity. When the speed of propagation is plotted against Damköhler for a range of Lewis a fold bifurcation is observed. It is shown that there exist stable positively and negatively propagating edges for some Damköhler and Lewis number pairs. It is further shown that changed local conditions can lead to a transition from positive (advancing into the unburnt gasses) to negative (receding) propagation.     

Experimental investigation of 40 kWth methane-assisted and self-sustained pulverized biomass flames

This paper presents an assessment of the effects of methane assistance on pulverized biomass swirl flames specifically regarding flow fields and flame structure. Experiments are carried out using a pilot-scale down-fired cylindrical combustion chamber equipped with a swirl burner and biomass/methane fuel mixtures. Studied conditions have an identical thermal output of 40 kW, with the thermal output share of methane gradually decreasing from 50% to 0% while the biomass share (walnut shells) increased from 50% to 100% (self-sustained condition). A detailed flow field characterization of the respective flames is conducted by in-flame, two-dimensional laser Doppler velocimetry measurements. These measurements are complemented by narrow-band flame imaging conducted at two different wavelengths (OH* and CH* radical band heads). Results show that the methane flames have a significant influence on the ignition and the determination of the flame flow field structure, including higher peak and overall velocities, as well as major changes in the ratio of tangential over axial velocity component. Further on, flame attachment of the self-sustained flame can be permanently achieved by the initial, short-term assistance of a methane flame with comparatively low thermal output. These observations are analysed and discussed, where higher measured velocities and overall changes in the flame structure between the self-sustained and the methane-assisted flames are attributed to important local expansion and momentum changes of the combustion gases introduced by the combustion of methane.     

Flame propagation in a small-scale parallel flow

We consider the propagation of laminar premixed flames in the presence of a parallel flow whose scale is smaller than the laminar flame thickness. The study addresses fundamental aspects with relevance to flame propagation in narrow channels, to the emerging micro-combustion technology, and to the understanding of the effect of small scales in a (turbulent) flow on the flame structure. In part, the study extends the results of a previous analytical study carried out in the thick flame asymptotic limit which has in particular addressed the validity of Damköhler's second hypothesis in the context of laminar steady parallel flows. Several new contributions are made here.

Analytical contributions include the derivation of an explicit formula for the effective speed of a premixed flame U _T in the presence of an oscillatory parallel flow whose scale ? (measured with the laminar flame thickness δ _L ) is small and amplitude A (measured with the laminar flame speed U _L ) is (1). The formula shows a quadratic dependence on both the amplitude and the scale of the flow. The validity of the formula is established analytically in two distinguished limits corresponding to (1) frequencies of oscillations (measured with the natural frequency of the flame U _L /δ _L ), and to higher frequencies of (A/?) (the natural frequency of the flow). The analytical study yields partial support of Damköhler's second hypothesis in that it shows that the flame behaves as a planar flame (to leading order) with an increased propagation speed which depends on both the scale and amplitude of the velocity fluctuation. However our formula for U _T contradicts the formula given by Damköhler in his original paper where U _T has a square root dependence on the scale and amplitude.

Numerical contributions include a significant set of two-dimensional calculations which determine the range of validity of the asymptotic findings. In particular, these account for volumetric heat loss and differential diffusion effects. Good agreement between the numerics and asymptotics is found in all cases, both for steady and oscillatory flows, at least in the expected range of validity of the asymptotics. The effect of the frequency of oscillation is also discussed. Additional related aspects such as the difference in the response of thin and thick flames to the combined effect of heat loss and fluid flow are also addressed. It is found for example that the sensitivity of thick flames to volumetric heat loss is negligibly affected by the parallel flow intensity, in marked contrast to the sensitivity of thin flames. Interestingly, and somewhat surprisingly, thin flames are found to be more resistant to heat loss when a flow is present, even for unit Lewis number; this ceases to be the case, however, when the Lewis number is large enough.