Flame dynamics and Markstein numbers in Hele-Shaw cells and porous media under Darcy's law

This study presents a hydrodynamic model demonstrating that premixed flames in Hele-Shaw cells and porous media under Darcy's law exhibit unique dynamics characterized by distinct Markstein numbers for curvature, tangential strain, and gravity, where viscosity-driven tangential velocity discontinuities fundamentally alter flame stability and instability mechanisms compared to conventional combustion.

The Gist

Imagine a fire burning. Usually, we think of fire spreading through open air, like a campfire or a candle flame. In those cases, the air moves freely, swirling and twisting around the flame. But what happens if you squeeze that fire into a very narrow space, like the tiny gap between two glass plates (a Hele-Shaw cell) or force it through a sponge-like material (porous media)?

This paper explores exactly that scenario. It turns out that when fire is squeezed into these tight spaces, it behaves less like a wild, swirling flame and more like a fluid flowing through a crowded hallway. The rules change completely.

Here is a breakdown of the paper's key discoveries using simple analogies:

1. The "Crowded Hallway" Effect (Darcy's Law)

In open air, fire follows the rules of standard fluid dynamics (Navier-Stokes), where air can swirl and spin freely. But in a narrow gap or a sponge, the walls or the material create so much friction that the air can't swirl. It has to move in a straight line, like a crowd of people trying to walk through a narrow corridor.

The scientists used Darcy's Law to describe this. Think of it as the "traffic rule" for fire in a jammed hallway: the air moves slowly and straight, and the friction from the walls is the most important thing, not the air's ability to spin.

2. The "Two-Headed" Flame (The Markstein Numbers)

In normal fire, scientists use a single number (called a Markstein number) to predict how a flame reacts when it gets curved or stretched. It's like having one "thermostat" that controls how the fire behaves.

The Big Discovery: Under these crowded conditions, that single thermostat breaks. The fire now needs two different thermostats:

• Thermostat A (Curvature): How the flame reacts when it bends (like a curved road).
• Thermostat B (Tangential Strain): How the flame reacts when the air along the flame's surface is being pulled or stretched.

Why the difference? In open air, the air moves smoothly across the flame. In a crowded hallway, the air on one side of the flame can suddenly move at a different speed than the air on the other side, creating a "speed bump" or a discontinuity. This means the fire has to react differently to bending versus being stretched sideways.

3. The "Gravity Switch"

The paper also found a third thermostat that only exists in these crowded conditions: Gravity.
In normal fire, gravity doesn't really change how the flame stretches. But in these narrow gaps, gravity can actually pull the air sideways along the flame, stretching it in a way that wouldn't happen otherwise. It's like having a third hand that can tug on the fire depending on which way is "up."

4. The "Traffic Refraction" (Streamline Bending)

Imagine a line of cars (air) driving toward a bend in the road (the flame).

• Normal Fire: The cars turn slightly as they cross the flame.
• Crowded Fire: Because of the friction and the "speed bump" mentioned earlier, the cars turn much more sharply. The paper calls this "augmented streamline refraction." It's like the cars are forced to take a much tighter turn than usual. This sharp turning makes the fire much more unstable and prone to breaking apart into little cells.

5. The "Counter-Flow" Surprise

The researchers looked at a specific setup where two streams of air blow toward each other, with a flame in the middle.

• Normal Fire: The stability of the flame depends on how much the air expands when it burns (density).
• Crowded Fire: The stability depends on how "thick" or "sticky" the air is (viscosity).
  It's a complete switch in the rules. A flame that would be stable in open air might get snuffed out instantly in a narrow gap because the "stickiness" of the air changes the game.

6. The "Wobbly Flame" (Instabilities)

Finally, the paper looks at what happens when the flame starts to wiggle or become unstable.

• Tight Squeeze (Strong Confinement): The flame behaves according to a specific mathematical equation (Michelson–Sivashinsky) that predicts long, wavy ripples. The friction makes these ripples grow faster and bigger than they would in open air.
• Loose Squeeze (Moderate Confinement): The flame behaves differently, following a different set of rules (Ginzburg–Landau), where it might form specific patterns or spots rather than long waves.

The Bottom Line

This paper tells us that confinement changes the personality of fire.
When you squeeze a flame into a narrow space or a porous material, it stops behaving like a free-spirited dancer and starts acting like a rigid, friction-bound flow. It needs new rules to predict how it will bend, stretch, and react to gravity. The most surprising part is that the "stickiness" of the air becomes more important than the "heaviness" of the air, and the flame develops a split personality, reacting differently to curves than to stretches.

This is crucial for understanding how fires behave in things like building insulation, underground gas reservoirs, or specialized industrial burners where space is tight.

Technical Summary

1. Problem Statement

The study addresses the fundamental differences in the hydrodynamics of premixed flames propagating in strongly confined environments, specifically Hele-Shaw cells (narrow gaps) and porous media.

• Context: In these systems, momentum losses due to friction with solid matrices or channel walls are dominant. Consequently, the flow is governed by Darcy's law rather than the Navier-Stokes equations.
• Gap in Knowledge: While previous studies have analyzed flame stability under Darcy's law, the specific role of confinement on Markstein numbers (which characterize the response of burning rates to curvature and stretch) remained unclear. Literature often applies Markstein number expressions derived for unconstrained flames to confined settings without verification.
• Objective: To derive a hydrodynamic model for confined flames, explicitly determine the Markstein numbers under Darcy's law, and analyze how confinement alters flame dynamics, stability, and extinction limits compared to conventional flames.

2. Methodology

The authors employ a hydrodynamic model treating the flame as a discontinuity surface (a thin flame approximation) within the framework of Darcy's law.

• Governing Equations:
  • The flow satisfies Darcy's law: $\mu \mathbf{v} = -\nabla p + \rho \mathbf{g}$.
  • Mass conservation ($\nabla \cdot \mathbf{v} = 0$) and pressure continuity across the flame front are enforced.
  • The flame front is defined by a level set $G(x,t)=0$.
• Asymptotic Analysis:
  • A multi-scale asymptotic analysis is performed using a one-step Arrhenius chemistry model with large activation energy.
  • This analysis generalizes previous work (e.g., Clavin-Williams) to non-unit Lewis numbers and accounts for the specific transport properties (viscosity $\mu$, thermal diffusivity $\lambda$, density $\rho$) in the unburnt and burnt gases.
• Key Parameters:
  • Density ratio: $r = \rho_u / \rho_b$.
  • Viscosity/Permeability ratio: $m = (\mu_u / \mu_b)(\kappa_b / \kappa_u)$. In Hele-Shaw cells, permeability is constant, so $m$ is the viscosity ratio.
  • Markstein Numbers: Defined as coefficients linking burning rate to flame stretch and curvature.

3. Key Contributions and Theoretical Findings

A. Discovery of Distinct Markstein Numbers

The most significant theoretical finding is that under Darcy's law, the burning rate depends on three distinct Markstein numbers, unlike the single Markstein number ($M_c$) typically found in Navier-Stokes flames.

1. Curvature Markstein Number ($M_c$): Relates to flame curvature.
2. Tangential Flow Markstein Number ($M_t$): Relates to tangential flow strain.
3. Gravity Markstein Number ($M_g$): A unique term arising from gravity-induced straining along the flame front.

The Burning Rate Formula:
The derived burning rate $\dot{m}$ is given by:
$$\dot{m} = 1 - M_c \mathcal{K} + (M_c - M_t)\nabla_t \cdot \mathbf{v}_t - M_g \nabla_t \cdot \mathbf{g}_t$$
Where $\mathcal{K}$ is kinematic stretch and $\nabla_t \cdot \mathbf{v}_t$ is the surface divergence of tangential velocity.

Why $M_c \neq M_t$?
In conventional flames (Navier-Stokes), tangential velocity is continuous across the flame, leading to $M_c = M_t$. Under Darcy's law, tangential velocity discontinuities are permitted due to viscosity variations ($\mu_u \neq \mu_b$). This discontinuity decouples the curvature and tangential strain effects, making $M_c$ and $M_t$ generally unequal.

B. Augmented Streamline Refraction

The paper demonstrates that streamlines refract (deflect) significantly more across a flame front in Darcy flows than in classical flows.

• Classical: Refraction ratio $\tan \theta_b / \tan \theta_u = 1/r$.
• Darcy: Refraction ratio $\tan \theta_b / \tan \theta_u = m/r$.
• Since $m \approx 1/3$ (typical viscosity ratio), streamlines refract approximately three times more under Darcy's law. This "augmented refraction" is a primary driver for enhanced hydrodynamic instabilities.

C. Application to Canonical Configurations

1. Radially Symmetric Flames: In purely radial flows, tangential strain is zero ($\nabla_t \cdot \mathbf{v}_t = 0$). Consequently, $M_c$ and $M_t$ do not appear separately, and the burning rate reduces to the classical form $\dot{m} = 1 - M_c \mathcal{K}$. The Darcy-specific effects vanish here.
2. Planar Counterflow Flames:
   • The strain rate jump across the flame is governed by the viscosity ratio $m$ rather than the density ratio $r$.
   • Extinction: The extinction strain rate is $A_{ext} = 1/M_t$. Because $M_t$ depends on viscosity variations, strained flames in confined media are more prone to extinction (lower $A_{ext}$) when $M_t < M_c$ (which occurs for certain Lewis numbers).
3. Linear Stability Analysis:
   • The dispersion relation for planar flames incorporates contributions from Darrieus-Landau (DL), Saffman-Taylor (ST), and Rayleigh-Taylor (RT) instabilities.
   • The growth rate depends on $M_t$ and $M_g$, but not on $M_c$ in the linear limit (due to the vanishing of the $(1-\mathbf{v}\cdot\mathbf{n})$ term).
   • Viscosity variations ($\mu \neq 1$) are shown to be the root cause of the flow's sensitivity to uniform imposed flows.

4. Nonlinear Dynamics and Confinement Regimes

The study analyzes the transition from strong confinement (Darcy limit) to weak confinement (Euler limit) using a unified dispersion relation.

• Strong Confinement (Low Permeability / Narrow Gap):
  • Dynamics follow a Michelson-Sivashinsky (MS) equation with modified coefficients.
  • The DL instability is amplified compared to unconstrained flames because the growth rate scales as $s \sim |k|(r-1)/(1+m)$, and since $1+m < 2$, the instability is stronger.
• Moderate Confinement:
  • Dynamics transition to Ginzburg-Landau (GL) type behavior.
  • The instability onset occurs at a non-zero critical wavenumber ($k_c \neq 0$), characteristic of finite-wavenumber instabilities.
• Weak Confinement (Unconstrained):
  • Recovers the classical MS equation and standard DL instability scaling.

5. Significance and Implications

• Theoretical Revision: The paper establishes that standard Markstein number formulations are insufficient for confined combustion. The inequality $M_c \neq M_t$ and the existence of $M_g$ are fundamental consequences of Darcy's law.
• Predictive Capability: The derived formulas allow for accurate prediction of flame extinction limits and instability growth rates in porous burners and Hele-Shaw experiments, which are critical for safety and efficiency in industrial applications (e.g., porous media combustion, micro-combustors).
• Physical Insight: The identification of "augmented streamline refraction" provides a clear physical mechanism explaining why confined flames are often more unstable than their free-space counterparts.
• Gravity Effects: The introduction of $M_g$ highlights that buoyancy effects in vertical confined channels can induce straining that significantly alters flame propagation, a factor previously overlooked in Darcy-based models.

In summary, this work provides a rigorous mathematical framework that bridges the gap between classical combustion theory and confined flow physics, revealing that confinement fundamentally alters the coupling between flame geometry, flow strain, and propagation speed.