A simplified model for coupling Darrieus-Landau and diffusive-thermal instabilities

This paper proposes a simplified phenomenological model that unifies the long-wave Darrieus-Landau and short-wave diffusive-thermal instabilities in premixed flames by introducing a cubic coupling term, revealing a distinguished crossover regime governed by a hydro-diffusive number that explains the emergence of chaotic, multi-scale flame front structures.

The Gist

The Big Picture: A Tug-of-War on a Burning Rope

Imagine a flame not as a smooth, steady sheet of fire, but as a wobbly, dancing rope. Scientists have long known that this "rope" wants to wiggle and break apart due to two different reasons.

1. The "Puff" Instability (Darrieus–Landau): When the flame burns, the hot gas expands and puffs out. This expansion makes the flame want to stretch into big, smooth, rolling waves (like a gentle ocean swell). If left alone, it would form large, lazy bumps.
2. The "Leak" Instability (Diffusive-Thermal): Heat and fuel molecules don't always move at the same speed. Sometimes, heat escapes faster than the fuel can keep up. This creates tiny, jagged ripples (like static on a TV screen) that want to tear the flame apart into fine, chaotic fuzz.

The Problem: For decades, scientists studied these two wiggles separately. They had one set of math for the big waves and a completely different set for the tiny fuzz. They assumed you could just add them together.

The Discovery: This paper argues that you can't just add them. When the flame is right on the edge of being unstable, these two forces fight each other in a very specific way. The author, Prabakaran Rajamanickam, proposes a new, simpler model to describe this fight.

---

The New Ingredient: The "Hydro-Diffusive Area"

The paper introduces a new concept called the Hydro-Diffusive Area.

• The Analogy: Imagine you are trying to smooth out a crumpled piece of paper.
  • The "Puff" force tries to make the paper billow out into big folds.
  • The "Leak" force tries to crumple it into tiny, sharp creases.
  • The Hydro-Diffusive Area is like the size of the table you are working on. It represents the specific "zone" where the big billows and the tiny creases crash into each other.

In the past, scientists ignored this zone. They thought the big waves and tiny creases happened in separate universes. This paper says: "No, they are interacting right here, and that interaction creates a cubic (three-sided) force that acts as a brake."

The Two Regimes: Smooth vs. Chaotic

The author shows that depending on how strong the "Leak" force is, the flame behaves in two very different ways:

1. The "Gentle Giant" Mode (Positive Markstein Number)

• What it is: The "Leak" force is weak. The flame is mostly dominated by the "Puff."
• The Result: The flame forms large, smooth, rolling waves with sharp points (cusps) at the top, like a series of gentle hills. It's chaotic, but in a predictable, rhythmic way.
• The Math: This is the classic behavior scientists already knew about (the Michelson–Sivashinsky equation).

2. The "Crossover" Mode (The New Discovery)

• What it is: The "Leak" force gets stronger, right on the edge of taking over. The two forces are now equal partners in the fight.
• The Result: This is where the magic happens. The flame doesn't just get bigger or smaller; it changes its texture.
  • Instead of just big hills, you get fine, cellular structures (like a honeycomb or a cracked mud surface).
  • The flame grows faster than before.
  • The movement becomes chaotic. The big "cusps" (the sharp points) try to form, but the tiny "fuzz" keeps destroying them, only for the big cusps to reform again. It's a constant, frantic cycle of building and breaking.
• The Analogy: Imagine a crowd of people trying to march in a line (the big waves). Suddenly, a group of kids starts running around their legs, tripping them and making them stumble (the tiny fuzz). The result isn't just a messy line; it's a chaotic, high-energy scrum where the formation is constantly being broken and reformed.

Why Does This Matter?

1. It Explains the "Fine Print": In real experiments, flames often show these tiny, complex cellular structures that old models couldn't explain well. This new model explains why they appear: it's the result of the two instabilities fighting on equal footing.
2. It's a Simpler Map: Instead of needing a massive, complicated computer simulation of every gas molecule to understand the flame, this paper offers a "simplified map" (a single equation) that captures the essence of the chaos.
3. The "Brake" That Never Lets Go: The paper highlights a special "braking" force (the cubic term) that kicks in when the flame is most unstable. Even if the usual stabilizing forces vanish, this new "Hydro-Diffusive Area" force keeps the flame from exploding into total disorder, keeping it in a state of "controlled chaos."

Summary in One Sentence

This paper proposes that when a flame is on the verge of instability, the big, smooth waves and the tiny, jagged ripples don't just happen side-by-side; they collide in a specific "zone" that creates a new, chaotic dance of fine-scale structures, which can be described by a simpler, unified mathematical model.

Technical Summary

1. Problem Statement

Premixed flame instability is traditionally analyzed through two separate, parallel frameworks:

• Darrieus–Landau (DL) Instability: A long-wavelength hydrodynamic instability driven by thermal expansion (density change). It is dominant when the Markstein number ($M$) is positive and of order unity.
• Diffusive-Thermal (DT) Instability: A short-wavelength instability arising from the disparity between thermal and molecular diffusion (preferential diffusion), occurring when the Lewis number is less than unity (implying $M < 0$).

The Gap: Existing theories treat these instabilities in isolation. While the Matalon–Matkowsky–Clavin–Joulin theory explains DL stabilization at short wavelengths, and Sivashinsky's work explains DT instability, a unified framework capturing the coupling between these mechanisms—particularly near the stability threshold where $M$ is small—has been lacking. Previous attempts often treated them as additive or relied on complex full conservation equations that are difficult to solve analytically.

2. Methodology

The author proposes a minimal phenomenological model based on a modified linear dispersion relation and asymptotic analysis.

• Modified Dispersion Relation: The paper identifies a missing leading-order coupling term between the long-wave DL mode and short-wave DT mode.
  • Standard DL: $\sigma_{DL} = \varepsilon|k| - Mk^2$
  • Standard DT: $\sigma_{DT} = -Mk^2 - 4k^4$ (for $M<0$)
  • Proposed Coupled Relation: $\sigma = \varepsilon|k| - Mk^2 - N|k|^3$
  • Here, $\varepsilon$ represents thermal expansion, $M$ is the Markstein number, and $N$ is a new parameter called the hydro-diffusive number.
• New Physical Parameter ($N$): Defined as $N = A/\delta_L^2$, where $A$ is the "hydro-diffusive area" (the characteristic area where hydrodynamic and diffusive processes interact) and $\delta_L$ is the flame thickness. The term $-N|k|^3$ acts as a nonlocal stabilizing mechanism.
• Asymptotic Regimes: The study analyzes two distinct regimes based on the scaling of $M$ relative to the thermal expansion parameter $\varepsilon$:
  1. Classical Regime ($M \gg \sqrt{\varepsilon N}$): Recovers the standard Michelson–Sivashinsky (MS) equation.
  2. Crossover Regime ($M \sim \sqrt{\varepsilon N}$): A distinguished limit where both instabilities participate at equal order.
• Numerical Simulation: The derived evolution equations are solved numerically in periodic domains to observe flame front dynamics, propagation speeds, and chaotic behavior.

3. Key Contributions

A. Identification of the Cubic Coupling Term

The paper's primary theoretical contribution is the derivation of the $-N|k|^3$ term in the dispersion relation.

• This term represents the lowest-order nonlocal correction consistent with symmetry.
• It becomes the dominant stabilizing mechanism when the classical quadratic stabilization ($Mk^2$) vanishes (i.e., as $M \to 0$).
• Unlike the quartic term ($-k^4$) used in previous DT models, the cubic term emerges naturally from the coupling of hydrodynamic and diffusive transport.

B. The Generalized Evolution Equation

In the crossover regime, the author derives a generalized evolution equation for the flame amplitude $f(x,t)$:
$$f_t + \frac{1}{2}f_x^2 = M f_{xx} + \varepsilon H(f_x) + N H(f_{xxx})$$
where $H$ denotes the Hilbert transform.

• This equation bridges the gap between the Michelson–Sivashinsky (MS) equation (pure DL) and the Kuramoto–Sivashinsky (KS) equation (pure DT).
• The term $N H(f_{xxx})$ introduces a nonlocal stabilizing effect that persists even when Markstein stabilization ($M$) is negligible.

C. Distinct Scaling Laws

The paper establishes that the crossover regime exhibits unique scaling laws compared to the classical DL regime:

• Wavenumber: $k \sim \sqrt{\varepsilon}$ (finer cellular structures than the classical $k \sim \varepsilon$).
• Growth Rate: $\sigma \sim \varepsilon^{3/2}$ (faster growth than the classical $\sigma \sim \varepsilon^2$).
• Amplitude: $f \sim \sqrt{\varepsilon}$ (smaller flame deformations).

4. Results

Numerical Observations

Numerical solutions reveal a distinctive chaotic regime in sufficiently large domains:

• Competition of Scales: The flame front exhibits a persistent competition between large-scale DL cusp structures and small-scale DT wrinkles.
• Chaotic Dynamics:
  • For $M < 0$ (Destabilizing DT): The system behaves like the Kuramoto–Sivashinsky equation, showing sustained chaotic activity with fine cellular structures.
  • For $M > 0$ (Stabilizing DT): In large domains, the system still exhibits chaos. The characteristic DL cusps are intermittently destroyed by short-wavelength wrinkles and then reform, creating complex recurrent cycles.
• Phase Portraits: The dynamics of the propagation speed $U$ versus its derivative show complex attractors, confirming the system acts as an effectively infinite-dimensional dynamical system.

Comparison with Previous Models

The paper compares its scaling ($M \sim \sqrt{\varepsilon}$) with Sivashinsky's earlier proposal ($M \sim \varepsilon^{2/3}$).

• The current model suggests the cubic stabilizing term becomes active earlier as the Markstein number decreases toward the instability threshold.
• Consequently, the cubic model predicts larger characteristic dimensions and slower growth rates compared to models relying on the quartic ($k^4$) term, which only dominates in deeper asymptotic limits.

5. Significance

• Unified Framework: The model provides a tractable, unified description of flame front instability that captures the essential coupled dynamics of DL and DT mechanisms without requiring the full complexity of Navier-Stokes conservation equations.
• Explanation of Fine-Scale Structures: It offers a theoretical explanation for the fine-scale cellular structures and accelerated growth rates observed in experiments near the stability threshold, which pure DL or pure DT theories cannot fully explain.
• New Physical Parameter: The introduction of the hydro-diffusive area ($A$) provides a concrete physical interpretation for the coupling strength, distinct from the Lewis-number-dependent Markstein number.
• Mathematical Insight: The paper highlights that the inclusion of the cubic term ($-( -\Delta)^{3/2} f$) fundamentally alters the singularity structure of the solutions, precluding the simple "pole decomposition" (finite number of cusps) characteristic of the classical MS equation. This suggests the dynamics are inherently richer and more complex.

In conclusion, this work advances the understanding of premixed flame instability by demonstrating that the interaction between hydrodynamic and diffusive-thermal effects is not merely additive but involves a critical, nonlocal cubic coupling that governs the transition between stability regimes.