![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
The stability of a thick planar premixed flame, propagating steadily in a direction transverse to that of unidirectional shear flow, is studied. A linear stability analysis is carried out in the asymptotic limit of infinitely large activation energy, yielding a dispersion relation. The relation characterises the coupling between Taylor dispersion (or shear-enhanced diffusion) and the flame thermo-diffusive instabilities, in terms of two main parameters, namely, the reactant Lewis number and the flow Peclet number
. The implications of the dispersion relation are discussed and various flame instabilities are identified and classified in the
-
plane. An important original finding is the demonstration that for values of the Peclet number exceeding a critical value, the classical cellular instability, commonly found for
, exists now for
but is absent when
. In fact, the cellular instability identified for
is shown to occur either through a finite-wavelength stationary bifurcation (also known as type-I
) or through a longwave stationary bifurcation (also known as type-II
). The latter type-II
bifurcation leads in the weakly nonlinear regime to a Kuramoto-Sivashinsky equation, which is determined. As for the oscillatory instability, usually encountered in the absence of Taylor dispersion in
mixtures, it is found to be absent if the Peclet number is large enough. The stability findings, which follow from the dispersion relation derived analytically, are complemented and examined numerically for a finite value of the Zeldovich number. The numerical study involves both computations of the eigenvalues of a linear stability boundary-value problem and numerical simulations of the time-dependent governing partial differential equations. The computations are found to be in good qualitative agreement with the analytical predictions.
1. Introduction
Thick premixed flames propagating in shear flows, such as flames propagating in a narrow-channel Poiseuille flow, are subject to enhanced diffusion in the longitudinal flow direction due to Taylor's dispersion mechanism, as reported in recent investigations [Citation1–4]. Meanwhile, diffusion transport in a direction transverse to the flow direction is only associated with molecular diffusion. Therefore, diffusion appears as effectively anisotropic which in turn significantly affect flame propagation and stability. The thermo-diffusive stability of premixed flames in the presence of Taylor dispersion has been recently studied, analytically and numerically, in the case of flames propagating in the longitudinal (or flow) directions. This has been done assuming adiabatic conditions in [Citation3] and accounting for heat losses in [Citation5]. The equivalent stability problem for flames aligned with the shear flow, propagating in a direction transverse to the flow has not been addressed yet. This problem is important to investigate on account of the anisotropy of diffusion aforementioned, which is expected to markedly alter the previous findings. Another motivation to consider the problem of transverse premixed-flame propagation is due to the interesting findings of a recent investigation on the effect of Taylor dispersion on nonpremixed flames aligned with the direction of a shear flow [Citation6]. This investigation demonstrates that Taylor dispersion can lead unexpectedly to cellular instability of a planar diffusion flame in mixtures with Lewis numbers above unity. This result has been argued in [Citation6] to present a plausible mechanism to explain the formation of so-called diffusion flame streets observed experimentally in non-premixed microcombustors [Citation7–9]. The question arises naturally therefore whether a similar cellular instability can occur for large Lewis numbers in the premixed case; namely, for flames propagating in the transverse direction, noting that such cellular instability for large Lewis numbers does not occur for flames propagating in the longitudinal direction [Citation3, Citation5]. One of the main objectives of this paper is to provide a clear answer to this question; as well as to describe the instabilities encountered and identify the conditions for their occurrence in terms of the parameters, mainly the Lewis number characterising the mixture and the Peclet number characterising the flow. In order to focus the analysis on these objectives, mainly concerned with the coupling between Taylor dispersion and the thermo-diffusive flame instabilities [[Citation10],[pp. 477–479][Citation11]], we shall adopt the thermo-diffusive approximation of constant density and constant transport properties and neglect further the effect of heat losses.
The paper is structured as follows. The problem formulation is given in §2 within a Hele-Shaw or two-dimensional channel configuration[Citation12–14], involving a unidirectional shear flow. In this configuration, the depth-averaged governing equations are written down and account for anisotropic diffusion, with enhanced diffusion in the longitudinal flow direction. A reformulation of the problem suitable for the limit of infinite Zeldovich number, , is provided in §3. This reformulation is used in the linear stability analysis of §4, which is carried out analytically and culminates in the derivation of a dispersion relation. The implications of the dispersion relation are discussed in detail in §5, where various flame instabilities are identified and classified in the parameters space. The stability findings, which follow from the dispersion relation derived analytically in the limit
, are then complemented and examined numerically in §6 for a finite value of the Zeldovich number,
. The numerical study in §6 involves both the computation of the eigenvalues of a linear stability boundary-value problem and the numerical simulations of the time-dependent governing partial differential equations for illustrative cases. Concluding remarks are provided in §7.
2. Problem formulation
Consider a narrow channel as depicted in Figure in which a planar premixed flame is propagating in the negative -direction with laminar flame speed
, in the presence of a shear flow. The shear flow is assumed, for simplicity, to have zero mean such as in the case of a Couette flow. Alternatively, in the case of Poiseuille flow for example, the frame of reference is chosen to be moving with the mean flow speed in the
-direction. The frame is furthermore assumed to move with speed
in the negative
-direction, so that the flame may be considered steady in the unperturbed state. In other words, the flow field adopted is of the form
(1)
(1) where U is the flow (maximum) amplitude and
is the scaled zero-mean shear-flow profile, such that e.g.
for a Couette flow and
for a plane Poiseuille flow.
Figure 1. Planar premixed flame in a channel of width , aligned with the direction (
) of a shear flow and propagating in the negative transverse direction (
).
![Figure 1. Planar premixed flame in a channel of width 2H, aligned with the direction (xˆ) of a shear flow and propagating in the negative transverse direction (yˆ).](/cms/asset/5008f25f-4ab6-4791-81c6-6a94f03b77c9/tctm_a_2254734_f0001_oc.jpg)
To model the chemistry, a single-step irreversible Arrhenius reaction is adopted, with pre-exponential factor B, activation energy E and heat release q per unit mass of fuel consumed. For sufficiently fuel-lean conditions, the fuel burning rate per unit volume can be written as , where ρ is the density assumed to be constant, R is the universal gas constant,
is the fuel mass fraction and T is the temperature. The adiabatic flame temperature
, the Zeldovich number β and the heat release parameter α are defined by
where
and
are the temperature and fuel mass fraction in the unburnt gas and
is the specific heat at constant pressure.
For non-dimensionalization, we shall use the laminar flame thickness as unit length and
as unit time, where
is the thermal diffusivity (assumed constant). Here
refers to the laminar flame speed (for
) which is given by
, where
is the (fuel) Lewis number. The dependent variables
and T are normalised by introducing
(2)
(2) We consider the limit
which leads upon averaging in the
-direction to a problem which is effectively two-dimensional in the first approximation, as shown in [Citation1–4]. The non-dimensional depth-averaged equations read
(3)
(3)
(4)
(4) where
is a parameter proportional to the Peclet number
,
The constant γ is a numerical coefficient, which is determined by the shear-flow profile and is given by
where
, so that
for the Poiseuille flow and
for the Couette flow, introduced above. We note that the parameter p quantifies the enhancement of diffusion by Taylor dispersion and that this enhancement is in the longitudinal x-direction, but not in the transverse y-direction. The boundary conditions in the y-direction are given by
(5)
(5)
(6)
(6) The main focus of this investigation is the determination of steady, one-dimensional solutions of the problem (Equation3
(3)
(3) )–(Equation6
(6)
(6) ) independent of x and t and their stability. The stability analysis is addressed analytically in the asymptotic limit of infinitely large Zeldovich number in §3–§5 and the predictions are then examined numerically for finite value of β in §6.
3. Formulation in the near-equidiffusional flame (NEF) limit
The stability analysis will be carried out asymptotically in the limit using the so-called near-equidiffusional flame (NEF) approximation based on the assumption that the Lewis number deviates little from unity [Citation15, p. 33]. Within this approximation, the reduced Lewis number
is
as
and Equations (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ) can be in terms of the leading-order temperature
and
as
(7)
(7)
(8)
(8) which are applicable outside an infinitely thin reaction sheet, given by
say. The equations are subject to the boundary conditions
(9)
(9)
(10)
(10) and the jump conditions
(11a)
(11a)
(11b)
(11b)
(11c)
(11c) applicable at
. Here we have used the notation
. It is worth pointing out that jump conditions (Equation11a
(11a)
(11a) ), which account for the presence of Taylor dispersion, are derived from an analysis of the structure of the reaction zone. Since the derivation is similar to that reported in [Citation2], its details are not included here.
4. Linear stability analysis for ![](//:0)
![](//:0)
We examine herein the stability of the steady planar flame solution, denoted by an overbar, which satisfies Equations (Equation7(7)
(7) )–(Equation11a
(11a)
(11a) ) with
and
and is given by
(12)
(12) To this solution, we add infinitesimal normal-mode disturbances such that
(13)
(13) where k denotes the real wavenumber and σ a constant (not necessarily real) characterising the growth rate/frequency of the perturbation.
The methodology used to obtain the dispersion relation is classical, see e.g. [Citation10]. We begin by substituting the perturbed solutions (Equation13
(13)
(13) ) into Equations (Equation7
(7)
(7) )–(Equation11a
(11a)
(11a) ). This leads to an eigenvalue boundary problem for the functions
and
which is given by
(14)
(14)
(15)
(15) applicable for
. These equations are subject to the boundary conditions
(16)
(16) and the linearised jump conditions
(17)
(17) applicable at y = 0, i.e.
.
The solution for is given by
(18)
(18) with
, satisfying the solvability condition or the dispersion relation
(19)
(19) Here
(20)
(20) Equation (Equation19
(19)
(19) ) can be solved for
, yielding in general three roots
,
and
, or less. Clearly, the stability of a mode with wavenumber κ is dictated by the root σ whose real part is equal to
. We shall denote this eigenvalue by
(21)
(21) Furthermore,
attains a maximum value at
say, with corresponding complex growth rate
. These values characterise the most unstable mode and are functions of l and λ,
(22)
(22)
5. Implications of the dispersion relation
In this section, we will examine the implications of the dispersion relation (Equation19(19)
(19) ) on flame stability using the notations (Equation20
(20)
(20) )–(Equation22
(22)
(22) ).
5.1. Illustrative cases and terminology
The function assumes various forms depending on the value of l and λ, each one resembling one of the dispersion curves depicted in Figure . In this figure, a first category of dispersion curves can be identified for which the maximum growth rate
is real and occurs for a non-zero value of the wavenumber,
; this is the case of subfigures (a) to (d). In this category, the instability may be termed as a cellular instability since the unstable planar flame solutions are expected to evolve into cellular structures, at least near the instability onset. A second category of dispersion curves can be identified for which the instability may be characterised as being an oscillatory instability, since it corresponds to
having a non-zero imaginary part. This is the case of subfigures (e) to (h). The oscillatory instability appears here as a longwave instability when
as in subfigures (e) and (f) or as a finite-wavelength instability when
as in subfigures (g) and (h). Note that in the case of subfigure (e), a finite wavelength cellular instability can occur instead of the oscillatory longwave instability, if the domain is not large enough. There are also cases such as in subfigure (i) where
and
is real, while
has a non-zero imaginary part except close to
. Solutions for such cases exhibit a non-oscillatory longwave instability near onset, at least in sufficiently large domains; in smaller domains, oscillations corresponding to non-zero wavenumbers are however expected near onset.
Figure 2. Growth rate versus wavenumber κ for selected values of l and λ. Solid lines indicate that
is real (
) and dashed lines indicate that
has a non-zero imaginary part (
). The horizontal and vertical scales are not indicated and are chosen individually for each subfigure for the sake of clarity.
![Figure 2. Growth rate Re(σmax) versus wavenumber κ for selected values of l and λ. Solid lines indicate that σmax is real (Im(σmax)=0) and dashed lines indicate that σmax has a non-zero imaginary part (Im(σmax)≠0). The horizontal and vertical scales are not indicated and are chosen individually for each subfigure for the sake of clarity.](/cms/asset/e56e545d-e497-4ef9-bcaa-767fdb4d3ac7/tctm_a_2254734_f0002_oc.jpg)
5.2. Stability regime diagram
A convenient way of summarising the various stability results is to delimit instability regions in an l-λ plane, which is equivalent to -
plane, as done in Figure . The colour scale in this figure represents the wavenumber
of the most unstable mode (with complex growth rate
). The solid lines in this figure separate the stable (white) region from the unstable ones, whereas the dashed linesFootnote1 for l>0 separate the region of cellular instability (above these lines) from that of non-cellular instability. Note that the region of non-cellular instability below the dashed lines is subdivided into three subregions by the dotted lines. To the right of the vertical dotted line (l>16) below the dashed lines, we have a non-oscillatory longwave instability and to its left
we have an oscillatory instability.
Figure 3. Stability regime diagram of a premixed flame propagating transversely to a shear flow in the l-λ plane, which is equivalent to the -
plane. The figure is based on the dispersion relation (Equation19
(19)
(19) ), obtained in the limit
. The colour scale represents the wavenumber
of the most unstable mode (with complex growth rate
). The types of the bifurcation curves, represented by solid lines, are discussed in §5.3.
![Figure 3. Stability regime diagram of a premixed flame propagating transversely to a shear flow in the l-λ plane, which is equivalent to the Le-Pe plane. The figure is based on the dispersion relation (Equation19(19) 2Γ2(1−Γ)+l(1−Γ+2σ+4λκ2)=0.(19) ), obtained in the limit β→∞. The colour scale represents the wavenumber κm of the most unstable mode (with complex growth rate σm). The types of the bifurcation curves, represented by solid lines, are discussed in §5.3.](/cms/asset/396c9619-28fc-4783-b1d9-7955ca0e1f96/tctm_a_2254734_f0003_oc.jpg)
As is well known [Citation11, pp. 477–479], when (or
), the cellular instability emerges for l<−2 and the oscillatory instability for
, as can be seen in the figure. The most striking observation in the presence of Taylor dispersion (
) is that the cellular instability region comprises positive values of l in addition to negative values. Specifically, cellular instability can now occur, according to our theoretical analysis, for any value of l such that
. It is worth noting in particular that for l>2, the cellular instability can be achieved if the Peclet number is large enough. Such cellular instability appears to be more accessible in real reactive mixtures than the oscillatory instability which requires large values of l, namely
(in the absence of heat loss). Another noteworthy observation pertinent to subunity Lewis number cases is that the cellular instability is hampered by Taylor dispersion when l<0 and suppressed completely for sufficiently strong shear flow,
.
5.3. Bifurcation curves
To characterise the bifurcations from stable to unstable cases in Figure occurring at the marginal condition , we shall adopt the terminology used in [[Citation16], pp. 75–81, [Citation17]]. According to this terminology, the types of bifurcations encountered, as denoted in Figure , are type-I
, type-II
, type-I
and type-III
bifurcations. These are illustrated schematically in Figure where the real growth rate
is plotted for the eigenvalue σ which crosses during the bifurcation from the left-half to the right-half complex plane. In the terminology adopted, the subscript s refers to a stationary or non-oscillatory bifurcation, also known as a zero-eigenvalue bifurcation and corresponds to
. Similarly, the subscript o refers to an oscillatory bifurcation and corresponds to
with
.
Figure 4. A schematic illustration of the four types of bifurcations identified in Figure . Plotted is as a function of wavenumber κ for the eigenvalue σ which crosses during the bifurcation from the left-half to the right-half complex plane. Solid lines indicate that σ is real (
) and dashed lines indicate that σ has a non-zero imaginary part (
).
![Figure 4. A schematic illustration of the four types of bifurcations identified in Figure 3. Plotted is Re{σ(κ)} as a function of wavenumber κ for the eigenvalue σ which crosses during the bifurcation from the left-half to the right-half complex plane. Solid lines indicate that σ is real (Im(σ)=0) and dashed lines indicate that σ has a non-zero imaginary part (Im(σ)≠0).](/cms/asset/d50631d3-f20e-4f64-a56e-c151ca6cb615/tctm_a_2254734_f0004_oc.jpg)
The bifurcation curves in Figure are determined as follows. The bifurcation curve labelled type-II in Figure is obtained from the condition
, which is clear from Figure (a), leading to the explicit relation
(23)
(23) Note that this expression is applicable in the domains
and
; for l>6, the boundary between stability and instability regions corresponds to different types of bifurcation curves, as indicated in Figure . Further details pertaining to this type-II
bifurcation are given in §5.4 mainly dedicated to the derivation of a Kuramoto-Sivashinsky equation.
We turn now to the type-I bifurcation curve identified in Figure . The curve is obtained from the conditions
and
as seen in Figure (b), which determine at the bifurcation the wavenumberFootnote2
and provide the relation
(24)
(24) We note that the type-I
and type-II
boundary curves introduced meet at a tricritical point
(25)
(25) This terminology is borrowed from phase transition theory [Citation18, p. 493] where a tricritical point refers to a point where first-order and second-order phase transition curves meet; such transition curves are mathematically analogous to our type-I
and type-II
bifurcation curves.Footnote3
Furthermore, the type-III bifurcation curve is determined from the requirements
,
at
; see Figure (c). The curve is found to be given by the equation
(26)
(26) which corresponds to a vertical line segment in Figure . The fact that this curve is vertical is straightforward consequence of the dispersion relation (Equation19
(19)
(19) ) being independent of λ when
.
Finally, the type-I bifurcation curve satisfies the conditions
,
at
; see Figure (d). This curve is computed numerically and is found to extend from the point
to the point
.
5.4. Further mathematical implications and weakly nonlinear analysis
The regime diagram of Figure suggests a rich variety of mathematical behaviours in the vicinity of the bifurcation curves and in particular near the tricritical point. Indeed, in the vicinity of the bifurcation curves, weakly nonlinear analyses can be carried out accounting for the coupling between Taylor dispersion and the thermo-diffusive instabilities. For example, in the vicinity of type-II bifurcation curve, a Kuramoto-Sivashinsky (KS) equation can be derived and, in the neighbourhood of type-I
bifurcation curve, a Swift-Hohenberg (SH) equation can be derived [Citation17]. More interestingly, near the tricritical point, sixth-order (in space) partial differential equations can be obtained, describing the flame evolution in the weakly nonlinear regime. These mathematical aspects are discussed elsewhere [Citation19]. Here we shall only address the flame dynamics near the type-II
bifurcation curve, which is classically described by a KS equation. The KS equation in our case takes the form
(27)
(27) where
(28)
(28) Note that the nonlinear term in (Equation27
(27)
(27) ) may be derived using a semi-heuristic kinematic argument as explained in [Citation20], and as done in [Citation5] in the similar problem of flames propagating in a longitudinal direction. The details of the derivation will not be repeated here as they are similar to those in [Citation5].
As for the linear part of (Equation27(27)
(27) ), this can be simply obtained from the dispersion relation (Equation19
(19)
(19) ). To this end, we note that Equation (Equation19
(19)
(19) ) has always a real root
such that
and that the Taylor expansion of
for small values of k is given by
(29)
(29) where
and
are given by
These are the coefficients of
and
in (Equation27
(27)
(27) ), since the linear part of (Equation27
(27)
(27) ) is equivalent to (Equation29
(29)
(29) ) in the case of normal modes
.
It is worth noting that the onset of the (type-II) instability corresponds to the condition
which is equivalent to the condition
used earlier in §5.3. This is so, provided the coefficient
of the fourth-derivative term is positive. Indeed, if
were negative, then the critical wavenumber
at the onset of instability would be necessarily non-zero, which contradicts the requirement
of a type-II
bifurcation. The condition
requires
or p<1, that is equivalently,
or
. Therefore the critical condition given by (Equation28
(28)
(28) ) is only applicable in this range of λ (or p), which determines the black type-II
bifurcation curve in Figure . For
and l>0, cellular flames still exist but the transition from the stable to the unstable region occurs through a type-I
bifurcation, represented by the red curve in Figure .
6. Computational results for finite Zeldovich number β
The theoretical results discussed so far are all based on the dispersion relation (Equation19(19)
(19) ) obtained in the asymptotic limit
. In this section, we shall carry out computations with a finite value of β in order to assess the applicability of theoretical predictions, at least qualitatively, and also to examine the nonlinear evolution of unstable flames. To this end, we consider the problem consisting of governing Equations (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ) and boundary conditions (Equation5
(5)
(5) )–(Equation6
(6)
(6) ), adopting the numerical values
and
, while varying the parameters
and
(or equivalently
).
6.1. Preliminary considerations
Before presenting the finite-β numerical results in the next two subsections, we make a few preliminary remarks which will facilitate the discussion and interpretation of the results. We begin by noting upon examining Equations (Equation3(3)
(3) ) and (Equation4
(4)
(4) ) that the presence of a flow-dependent effective Lewis number in the x-direction, as noted in [Citation2, Citation4], namely,
(30)
(30) where use has been made of the definitions
and
.
An important implication of formula (Equation30(30)
(30) ) is that
as
and
as
. In other words, the effect of a shear flow is such that a weakly diffusing reactant (
) appears effectively as strongly diffusing (in the x-direction,
) and a strongly diffusing reactant (
) appears effectively as weakly diffusing (
) provided the Peclet number is large enough, more precisely when
exceeds the value
. In fact, this value is determined from the fact that
irrespective of
when
, which follows from (Equation30
(30)
(30) ).
6.2. Linear stability analysis based on eigen-boundary value problem
We begin by examining the linear stability of the planar flame solution with the finite value adopted. This base solution satisfies Equations (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ) with
. The scaled burning speed S, which is equal to unity in the limit
, is now computed numerically and plotted as a function of the Lewis number
in Figure . If we denote the dependent variables of the base state as
and
, then the stability of the base solution to small perturbations can be studied by introducing
(31)
(31) into the problem (Equation3
(3)
(3) )–(Equation6
(6)
(6) ). The linear stability is then described by the eigen-boundary value problem
(32)
(32)
(33)
(33) with
(34)
(34)
(35)
(35) It is worth noting the presence in the equations of the longitudinal effective Lewis number
defined in (Equation30
(30)
(30) ), in addition to the Lewis number
. The problem possesses a discrete spectrum of eigenvalues σ, in which the eigenvalue with the maximum real part (growth rate) determines the stability of the base solution. The corresponding computational results are summarised in Figure , where the stability-instability regions are delimited in the
-λ plane. The solid lines in the figure separates the stable (white) region from the unstable (grey) regions. The dash-dotted line represents the condition
, which takes the form
on using
in (Equation30
(30)
(30) ).
Figure 6. Stability regime diagram of a premixed flame propagating transversely to a shear flow in the -λ plane. The figure is based on the computation of the eigenvalues of problem (Equation32
(32)
(32) )–(Equation35
(35)
(35) ) with
and
. The grey region corresponds to unstable flames and the white to stable flames. Compare with Figure , where analogous notations are used.
![Figure 6. Stability regime diagram of a premixed flame propagating transversely to a shear flow in the Le-λ plane. The figure is based on the computation of the eigenvalues of problem (Equation32(32) 1{Le}d2y~Fdy2−Sdy~Fdy−k2{Le}xy~F−{y~F+βy¯Fθ~[1−α(1−θ¯)]2}β22{Le}exp[β(θ¯−1)1+α(θ¯−1)]=σy~F,(32) )–(Equation35(35) y~F=0,dθ~dy=0as y→+∞.(35) ) with β=10 and α=0.85. The grey region corresponds to unstable flames and the white to stable flames. Compare with Figure 3, where analogous notations are used.](/cms/asset/cc41bb6e-ae39-4277-a160-62a8dcd4c921/tctm_a_2254734_f0006_oc.jpg)
Comparing Figure , computed for , with Figure corresponding to
, the following conclusions can be drawn. First, we note the good qualitative agreement between the two figures, notably, regarding the bifurcation curves separating the stability from instability regions. In particular note the presence of a tricritical point at
, identified in Figure . Second, the computations confirm the existence of the cellular instability region for
, an original finding of the asymptotic analysis. More precisely, two cellular instability regions appear in Figure , whose locations are compatible with the conditions
when
and
when
, or equivalently
for any
, as argued in §6.1. Further the flame is found to be stable for all values of λ when
in the case
, while the equivalent stability range predicted by the asymptotic analysis,
, provides the stability range
for
.
6.3. Time-dependent numerical simulations
Finally, we present time-dependent numerical simulations simply to illustrate the occurrence of the cellular instability for and the flame long-time evolution. The simulations are based on the numerical solution of the time-dependent problem (Equation3
(3)
(3) )–(Equation6
(6)
(6) ), with periodic conditions in the x-direction, starting from an initial condition corresponding to a steady, planar premixed flame. The computations are carried out using COMSOL Multiphysics software as described in [Citation4, Citation5]. Two cases are considered corresponding to
and
, whose results are shown in Figures and , respectively. Shown are reaction-rate
fields at selected values of time t. These illustrate the development of the instability and the flame evolution into a cellular structure, which settles into an apparently stable state at large times.
7. Conclusions
In this paper, we have examined the effect of Taylor dispersion, or shear-enhanced diffusion on the stability of a premixed flame, propagating in a direction transverse to that of unidirectional shear flow. A simple thermo-diffusive model is adopted corresponding to a Hele-Shaw channel, whose walls are assumed to be adiabatic and closely spaced. In this configuration, a simple shear flow is prescribed, corresponding for example to a Couette flow. Upon depth-averaging, the problem is governed by two-dimensional transport equations in which diffusion is anisotropic with shear-enhanced diffusion coefficients in the longitudinal flow direction. A linear stability analysis is carried out in the limit , where β is the Zeldovich number. A simple dispersion relation (Equation19
(19)
(19) ) is derived analytically involving the Lewis number
and the Peclet number
as parameters. A stability regime diagram (Figure ) is constructed in the
-
plane, which categorises the various instabilities and bifurcations encountered.
A remarkable finding is the demonstration that the classical cellular instability, which is usually expected to occur in mixtures, exists now only for
mixtures when the Peclet number exceeds a critical value. The necessary conditions needed to observe cellular flames in
mixtures include that the geometry must be slender such as in the case of a Hele-Shaw channel, that the Peclet number
must exceed a critical value above
and that the flame must be aligned with the flow direction. Furthermore, the cellular instability identified for
is shown to occur either through a finite-wavelength stationary bifurcation (also known as type-I
) or through a longwave stationary bifurcation (also known as type-II
). In the weakly nonlinear regime, a Kuramoto-Sivashinsky equation in the vicinity of type-II
bifurcation is derived. Moreover, it is found that the oscillatory instability, usually encountered in
mixtures persists under the influence of Taylor dispersion if the Peclet number is below a critical value and disappears above this value.
The stability results aforementioned, which follow from the dispersion relation obtained in the limit , are complemented by numerical computations carried out for a finite value of the Zeldovich number. The computations involve the determination of the eigenvalues of a linear stability boundary-value problem and numerical simulations of the time-dependent governing partial differential equations. The numerical results are found to be in good qualitative agreement with the analytical predictions. In particular, a stability regime diagram (Figure ) computed for
is found to be consistent with the stability regime diagram of Figure , corresponding to
.
To close this paper, we note that the stability results obtained for the case of transverse flame propagation differ markedly from those obtained for flames propagating in the longitudinal direction [Citation3]. For example, the appearance of cellular instability in mixtures when the Peclet number is above a critical value, does not occur in the case of longitudinal propagation. In fact, in the latter case, only stable flames are encountered when
and
. A natural follow-up of the current work is to extend the stability analysis to flames propagating in an arbitrary direction relative to that of the shear flow. It is also worthwhile to investigate the influence of heat losses on our predictions, following [Citation5], since such influence is significant in practice. Interestingly, the cellular instability identified herein for
mixtures has been recently shown to also take place for diffusion flames aligned with the direction of the shear flow [Citation6]. This unexpected instability has been proposed in [Citation6] as a plausible mechanism for the formation of diffusion flame streets observed in experiments [Citation7–9]. It is desirable to have similar experiments in the premixed case to test the presence of the cellular instability identified herein for
mixtures in the presence of a strong shear flow.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
Notes
1 The dashed curves meet at a cusp located at . The crossing from the non-cellular to the cellular region as λ is increased involves a discontinuous jump in
for l<16 and a continuous transition for l>16. Similar transitions have been observed in diffusion flames which have been described with more detail in [Citation6].
2 An explicit relation for can be obtained by substituting the expression for λ in (Equation24
(24)
(24) ) into the equation
, which follows from the condition
. From this relation, we can deduce that
as
.
3 As in phase transition theory [Citation18, p. 496], the curves of type-I (Equation24
(24)
(24) ) and type-II
(Equation23
(23)
(23) ) are continuous and have continuous first derivatives but discontinuous second and higher-order derivatives at the tricritical point.
References
- P. Pearce and J. Daou, Taylor dispersion and thermal expansion effects on flame propagation in a narrow channel, J. Fluid Mech. 754 (2014), pp. 161–183.
- J. Daou, P. Pearce, and F. Al-Malki, Taylor dispersion in premixed combustion: Questions from turbulent combustion answered for laminar flames, Phys. Rev. Fluids 3 (2018), p. 023201.
- J. Daou, Effect of Taylor dispersion on the thermo-diffusive instabilities of flames in a Hele–Shaw burner, Combust. Theory Model 25 (2021), pp. 765–783.
- P. Rajamanickam and J. Daou, A thick reaction zone model for premixed flames in two-dimensional channels, Combust. Theory Model 27 (2023), pp. 487–507.
- J. Daou, A. Kelly, and J. Landel, Flame stability under flow-induced anisotropic diffusion and heat loss, Combust. Flame 248 (2023), p. 112588.
- P. Rajamanickam, A. Kelly, and J. Daou, Stability of diffusion flames under shear flow: Taylor dispersion and the formation of flame streets, Combust. Flame 257 (2023), p. 113003.
- C. Miesse, R.I. Masel, M. Short, and M.A. Shannon, Diffusion flame instabilities in a 0.75 mm non-premixed microburner, Proc. Combust. Inst. 30 (2005), pp. 2499–2507.
- C. Miesse, R.I. Masel, M. Short, and M.A. Shannon, Experimental observations of methane–oxygen diffusion flame structure in a sub-millimetre microburner, Combust. Theory Model 9 (2005), pp. 77–92.
- B. Xu and Y. Ju, Studies on non-premixed flame streets in a mesoscale channel, Proc. Combust. Inst.32 (2009), pp. 1375–1382.
- G.I. Sivashinsky, Nonlinear analysis of hydrodynamic instability in laminar flames–I. Derivation of basic equations, Acta Astronaut 4 (1977), pp. 1177–1206.
- P. Clavin and G. Searby, Combustion waves and fronts in flows: Flames, shocks, detonations, ablation fronts and explosion of stars, Cambridge University Press, Cambridge, 2016.
- G. Joulin and G.I. Sivashinsky, Influence of momentum and heat losses on the large-scale stability of quasi-2D premixed flames, Combust. Sci. Technol. 98 (1994), pp. 11–23.
- D. Fernández-Galisteo, V.N. Kurdyumov, and P.D. Ronney, Analysis of premixed flame propagation between two closely-spaced parallel plates, Combust. Flame 190 (2018), pp. 133–145.
- E. Al Sarraf, C. Almarcha, J. Quinard, B. Radisson, B. Denet, and P. Garcia-Ybarra, Darrieus–Landau instability and markstein numbers of premixed flames in a Hele–Shaw cell, Proc. Combust. Inst. 37 (2019), pp. 1783–1789.
- J.D. Buckmaster and G.S.S. Ludford, Lectures on mathematical combustion, SIAM, Philadelphia, Pennsylvania, 1983.
- M. Cross and H. Greenside, Pattern formation and dynamics in nonequilibrium systems, Cambridge University Press, Cambridge, 2009.
- M.C. Cross and P.C. Hohenberg, Pattern formation outside of equilibrium, Rev. Mod. Phys. 65 (1993), p. 851.
- L.D. Landau and E.M. Lifshitz, Statistical physics: Volume 5, Elsevier, Oxford, 2013.
- P. Rajamanickam and J. Daou, Tricritical point as a crossover between type- Is and type- IIs bifurcations, Prog. Scale Model. Int. J. (2023).
- G.I. Sivashinsky, Instabilities, pattern formation, and turbulence in flames, Annu. Rev. Fluid Mech.15 (1983), pp. 179–199.