Scaling Laws for Thermodiffusively Unstable Lean Premixed Turbulent Hydrogen-Air Flames

This study evaluates and adapts two existing models for thermodiffusive instabilities in lean premixed turbulent hydrogen-air flames using 91 DNS cases, revealing that while both models are physically equivalent and reduce to a Karlovitz-number-dependent form under typical conditions, explicit consideration of specific stability parameters is necessary for accurate scaling in ultra-low flame speed regimes.

The Gist

The "Super-Fast" Hydrogen Fire: A Tale of Two Rules

Imagine you are trying to predict how fast a fire will burn. Now, imagine that fire isn't just burning wood, but a mixture of hydrogen and air. This isn't your average campfire; it's a lean, high-tech flame used in things like jet engines and gas turbines.

The problem is that hydrogen flames are tricky. They are prone to a weird glitch called thermodiffusive instability. Think of it like a crowd of people trying to walk through a narrow door. If the people (hydrogen molecules) move much faster than the air around them, they start jostling, creating ripples and bumps on the edge of the fire. These ripples make the fire front "wrinkle," and just like a crumpled piece of paper has more surface area than a flat one, a wrinkled fire burns much, much faster than a smooth one.

Scientists have been trying to write a "rulebook" (a mathematical model) to predict exactly how much faster this fire will burn based on how turbulent the air is. Recently, two different teams came up with two different rulebooks.

The Two Competing Rulebooks

1. The "Instability" Rulebook (The $\omega^2$ Model): This team looked at the physics of why the fire wrinkles. They used a complex number (let's call it the "Wobble Factor") to measure how unstable the flame is. Their rule says: "If the flame wobbles a lot, multiply the speed by this factor."
2. The "Ratio" Rulebook (The $Ze/Pe$ Model): This team took a different approach. They looked at the ratio of two other numbers: how reactive the fuel is versus how easily it spreads. Their rule says: "If this ratio is high, the fire speeds up."

For a long time, scientists didn't know which rulebook was better, or if they were actually saying the same thing in different languages. They had tested these rules on simple, calm fires in a box, but they weren't sure if the rules would work for real-world jet flames, which are messy, turbulent, and change shape as they fly.

The Great Experiment: 91 Simulations

The authors of this paper decided to settle the debate. They didn't just run one test; they ran 91 massive computer simulations (like running a video game 91 times with different settings).

They tested:

• Different pressures: From normal air pressure to the crushing pressure inside a jet engine.
• Different temperatures: From cool air to scorching hot exhaust.
• Different turbulence: From a gentle breeze to a violent storm.
• Different shapes: From simple box flames to complex jet flames shooting out of a nozzle.

The Big Discovery: Two Worlds, One Truth

After crunching the numbers, the authors found something fascinating. The world of hydrogen flames actually has two distinct regimes (two different "worlds"), and the rulebooks behave differently in each:

1. The "Normal" World (Low Pressure)

This covers most gas turbines and burners.

• The Result: Both rulebooks turned out to be identical.
• The Analogy: Imagine two people giving you directions to a store. One says, "Walk 5 blocks, then turn left." The other says, "Walk 5 blocks, then turn left." They used different words to describe the steps, but they gave the exact same instructions.
• The Takeaway: In this regime, you don't need to worry about which complex formula you use. They both simplify to a single, easy rule: The faster the turbulence, the faster the fire burns, and the amount of "wobble" doesn't change the basic math.

2. The "Extreme" World (High Pressure)

This covers the most extreme conditions, like inside an engine during a specific part of the cycle.

• The Result: Here, the two rulebooks look different on the surface, but they are actually deeply connected.
• The Analogy: Imagine one person says, "The car is fast because of the engine," and another says, "The car is fast because of the tires." At first, they seem to be arguing about different parts. But the scientists realized that in this extreme world, the engine and the tires are so tightly linked that you can't have one without the other. They are two sides of the same coin.
• The Takeaway: In these extreme conditions, you do need to include the specific "Wobble Factor" or the "Ratio" to get the math right. However, the authors proved that these two factors are physically related. If you understand one, you mathematically understand the other.

Why Does This Matter?

Think of designing a clean-energy engine as trying to build a race car.

• If your math for how fast the engine burns fuel is wrong, the car might stall, or worse, explode.
• For a long time, engineers had to guess which "rulebook" to use for hydrogen engines.
• This paper is like a unified manual. It tells engineers: "Don't panic. Whether you are in a normal engine or an extreme one, these two models are actually saying the same thing. You can trust them to predict how fast the fire will burn."

The Bottom Line

The scientists took two competing theories, tested them against a mountain of data (91 simulations!), and found that they are physically equivalent.

• In normal conditions, both models simplify to the same easy-to-use formula.
• In extreme conditions, they look different but are mathematically linked, meaning they describe the same physical reality.

This is a huge step forward for designing sustainable, hydrogen-powered engines. It means we can finally build better, safer, and more efficient engines because we finally have a reliable map for how hydrogen flames behave in the wild.

Technical Summary

1. Problem Statement

Lean premixed hydrogen-air flames are highly susceptible to thermodiffusive (TD) instabilities. Due to the high diffusivity of hydrogen (Lewis number $Le < 1$), these instabilities amplify flame front perturbations, leading to cellular structures that make the flame locally hotter, thinner, and faster. This significantly enhances local reactivity, often causing the local flame speed to exceed the unstretched laminar burning velocity by several times.

While turbulence-chemistry interaction models exist, there is a lack of a first-principles model for the stretch factor ($I_0$) that explicitly accounts for TD responses under varying local turbulence and thermodynamic conditions. Two recent empirical models have been proposed to address this:

1. The $\omega^2$-model (Howarth et al.): Uses an instability parameter derived from linear stability theory.
2. The $Ze/Pe$-model (Rieth et al.): Uses the ratio of the Zel'dovich number ($Ze$) to the Peclet number ($Pe$).

The Gap: These models have been validated primarily in simple, statistically homogeneous, isotropic turbulence ("flame-in-a-box"). It remains unclear how they perform in technically relevant, complex configurations (e.g., turbulent jet flames with shear and inhomogeneities) and whether the two formulations are physically equivalent or applicable across different pressure regimes.

2. Methodology

The authors conducted a systematic evaluation and consolidation of the two models using a comprehensive dataset of 91 Direct Numerical Simulation (DNS) cases.

• Dataset Scope:

  • Configurations: Forced and decaying homogeneous isotropic turbulence, temporally evolving shear layers, and piloted slot jet flames.
  • Conditions: Pressures from 1 to 40 atm, equivalence ratios ($\phi$) from 0.2 to 0.57, unburnt temperatures ($T_u$) from 300 to 750 K, and Karlovitz numbers ($Ka$) ranging from 1.3 to 2580.
  • Key Distinction: The study includes both "low-pressure" regimes (typical of gas turbines) and "high-pressure" regimes (characterized by ultra-low flame speeds, typical of internal combustion engines with exhaust gas recirculation).

• Model Evaluation:

  • The authors fitted the DNS data to both the original $\omega^2$-model and a reformulated $Ze/Pe$-model.
  • The $Ze/Pe$-model was modified to satisfy the laminar limit ($I_0 \to I_0^*$ as $Ka \to 0$) and to use the refined Karlovitz number ($Ka^*$) based on the stretched laminar flame speed.
  • Jet Flame Validation: Jet flame data were excluded from the model fitting process to serve as an independent test of predictive capability for complex, shear-driven flows.

3. Key Contributions

1. Systematic Consolidation: The first comprehensive comparison of the $\omega^2$ and $Ze/Pe$ models across a wide spectrum of thermodynamic conditions and flame geometries.
2. Identification of Two Distinct Regimes: The analysis revealed that a single scaling law cannot describe the entire dataset. Instead, the behavior splits into two distinct regimes based on a critical pressure ($\Pi_c$):
   • Low-Pressure Regime ($p \le \Pi_c$): Typical of gas turbines.
   • High-Pressure Regime ($p > \Pi_c$): Characterized by ultra-lean mixtures and high pressures.
3. Physical Equivalence Proof: Through dimensional analysis and Taylor series expansions, the authors demonstrated that both models are governed by the same underlying functional relationship involving the non-dimensional group $Ze(1-Le)$.
4. Adapted Scaling Relation: Proposed a unified framework where the models reduce to a simpler form in the low-pressure regime but require explicit dependence on $\omega^2$ or $Ze/Pe$ in the high-pressure regime.

4. Key Results

A. Performance in Low-Pressure Regime

• Convergence: In the low-pressure regime, both models reduce to an identical functional form:
  $$I_0 \approx (1 + q Ka^{*m}) I_0^*$$
• Independence: The explicit dependence on $\omega^2$ or $Ze/Pe$ becomes negligible because the prefactors vary by less than 20%. The stretch factor depends primarily on the Karlovitz number ($Ka^*$) and the laminar stretch factor ($I_0^*$).
• Accuracy: Both models predict the DNS data well (Mean Absolute Percentage Error, MAPE $\approx$ 8–10%) and perform reasonably well on the independent jet flame data.

B. Performance in High-Pressure Regime

• Divergence: In the high-pressure regime, the models do not converge to a simple form. The explicit consideration of the instability parameter ($\omega^2$) or the ratio ($Ze/Pe$) is required for accurate scaling.
• Accuracy: The modified models still capture the trends well (MAPE $\approx$ 6–8% for high-pressure cases), but the prefactors are no longer constant.
• Jet Flames: The models show predictive capability for jet flames but exhibit reduced accuracy near the flame tip, likely due to flame-flame interactions and expansion effects not captured by the local Karlovitz number.

C. Physical Unification

• The authors derived an analytical link showing that the scaling coefficients in both models are related through the physical relation:
  $$C Ze Le^{-1/3} \approx \exp(0.022 \omega^2)$$
• This confirms that despite different mathematical formulations, both models describe the same physical mechanism: the interplay between the Zel'dovich number and the Lewis number.

5. Significance

• Model-Based Design: This work provides the necessary validation and reconciliation of turbulence-chemistry interaction models required for the design of sustainable hydrogen combustion systems (e.g., gas turbines and engines).
• Bridging Theory and Application: By testing models developed in idealized turbulence against complex jet flames, the study validates the transferability of these scaling laws to real-world engineering applications.
• Regime-Specific Modeling: The identification of distinct low- and high-pressure regimes prevents the misuse of simplified scaling laws in conditions where they fail (specifically high-pressure, ultra-lean conditions).
• Future Development: The unification of the models into a single physical framework ($Ze(1-Le)$) offers a robust foundation for developing next-generation low-order combustion models that can accurately predict burning rates in thermodiffusively unstable hydrogen flames.