Stationary Framelet Model
Stationary Framelet: Theoretical Foundations
Overview
Professor, what is a flamelet model?
The flamelet model is a model for non-premixed turbulent combustion proposed by Peters in 1984. It views a turbulent diffusion flame as "an ensemble of numerous thin laminar flames (flamelets)" and reduces the internal flame structure to a one-dimensional problem in mixture fraction $Z$.
So you can reduce 3D turbulent combustion to 1D?
Yes. The flamelet assumption is that "the flame thickness is thinner than the smallest turbulent scale (the Kolmogorov scale)". In this case, the internal flame structure becomes locally equivalent to a 1D laminar counterflow flame.
Flamelet Equation
Please explain the flamelet equation.
The flamelet equation for chemical species $Y_i$ can be written as follows.
Here, $\chi$ is the scalar dissipation rate, representing diffusion in $Z$ space.
For a steady flamelet, $\partial/\partial t = 0$, right?
Correct. For a steady flamelet, the time derivative term becomes zero, resulting in an ordinary differential equation with the dissipation rate $\chi_{st}$ (value at the stoichiometric surface) as a parameter. As $\chi_{st}$ increases, the reaction cannot keep up, eventually leading to extinction (quenching dissipation rate $\chi_q$).
S-Curve (S-curve)
What is an S-curve?
Plotting the maximum temperature of the steady flamelet solution against $\chi_{st}$ yields an S-shaped curve. The upper branch represents the burning state, the lower branch the unburned state, and the middle branch is an unstable solution.
What is the typical value of $\chi_q$?
For methane/air, $\chi_q \approx 20-50$ s$^{-1}$; for hydrogen/air, $\chi_q \approx 1000$ s$^{-1}$, which is very large. This means hydrogen is difficult to extinguish.
Flamelet Library
A flamelet library is created by computing multiple steady flamelet solutions with $\chi_{st}$ as a parameter and storing them as a 2D table of $(Z, \chi_{st})$. If the turbulent $\beta$-PDF averaging is applied in advance, it becomes a 3D lookup table of $(\widetilde{Z}, \widetilde{Z''^2}, \widetilde{\chi_{st}})$.
So the core of the flamelet model is the combination of "a 1D laminar flame table + turbulent PDF".
Exactly. The flame structure including detailed chemical reactions is solved offline, and in 3D CFD, only table lookup is required, resulting in very low computational cost.
The Man Who Turned Flames into a "Database" – Norbert Peters' Flamelet Philosophy
The idea proposed by Norbert Peters, the originator of the flamelet model, in 1984 can be summarized in one sentence: "A flame in turbulence is ultimately an ensemble of small laminar flames." If that's the case, then precomputing laminar flames and storing them in a table means that during the main calculation, you only need to reference the mixture fraction. This "precomputation → table lookup" concept dramatically reduced the computational cost of combustion CFD. Peters developed this concept at RWTH Aachen University, creating a lineage that later led to models like FGM (Flamelet-Generated Manifolds) and FPV (Flamelet Progress Variable). Modern gas turbine design cannot be discussed without his legacy.
Computational Methods for Stationary Framelet
Details of Numerical Methods
How do you create a flamelet library?
Library construction has two steps: (1) Calculation of counterflow diffusion flames, (2) Table creation via PDF integration.
Counterflow Diffusion Flame Calculation
How is it calculated specifically?
Solve 1D counterflow diffusion flames using Cantera, FlameMaster, OPPDIF (CHEMKIN-PRO), etc. Gradually increase the strain rate (strain rate) $a$ and track up to the extinction point. The relationship between $a$ and $\chi_{st}$ is given by $\chi_{st} \approx a \cdot f(Z_{st})$.
Typical library construction parameters:
| Parameter | Recommended Value | Remarks |
|---|---|---|
| Grid points in $Z$ direction | 128-256 | Concentrated near $Z_{st}$ |
| $\chi_{st}$ discretization | 30-50 points | Logarithmic spacing |
| Strain rate range | 1 - $a_q$ s$^{-1}$ | Up to extinction |
| Reaction mechanism | GRI-Mech 3.0, etc. | Use detailed mechanism |
Implementation in Fluent
How do you set it up in Fluent?
Configure the following within Fluent's Non-Premixed Combustion model.
1. Select Flamelet Model (switch from Equilibrium Chemistry)
2. Import reaction mechanism in CHEMKIN format
3. Set Number of Flamelets (default 20, recommended 30-50)
4. Set PDF table resolution
5. After calculation starts, transport equations for $\widetilde{Z}$, $\widetilde{Z''^2}$ are solved, and temperature/species are determined via table lookup.
Relationship with FGM (Flamelet Generated Manifold)
How is FGM different from the flamelet model?
FGM is a method proposed by van Oijen et al. (2000) that constructs a low-dimensional manifold from flamelet solutions. In addition to steady flamelets, it introduces a progress variable $C$, enabling representation of transient processes like ignition and extinction.
Comparison of steady flamelet and FGM:
| Characteristic | Steady Flamelet | FGM |
|---|---|---|
| Table dimension | 2-3D ($Z$, $Z''$, $\chi$) | 3-4D ($Z$, $Z''$, $C$, $C''$) |
| Extinction representation | Possible via S-curve | Naturally represented via Progress Variable |
| Auto-ignition | Difficult | Possible |
| Partial premixing | Difficult | Possible |
| Computational cost | Very low | Low (table lookup) |
So FGM is more versatile.
Yes. FGM is becoming mainstream in STAR-CCM+ and OpenFOAM. Fluent has also strengthened its FGM option since R2. However, steady flamelet remains widely used as the simplest and most stable method.
The key to the flamelet model is "creating a high-quality table."
Exactly. The table resolution and the validity of the reaction mechanism directly determine the model's accuracy.
Results Change with the "Shape" of the PDF – Comparison of β-PDF and Clipped Gaussian
An unavoidable aspect of the numerical implementation of flamelet/PDF models is "the choice of the PDF shape for the mixture fraction." The most common is the β-PDF (beta distribution), which is also the default in ANSYS Fluent. However, while β-PDF is mathematically convenient, it can only represent unimodal shapes and has the problem of overly steep behavior at the two boundaries (pure fuel / pure oxidizer). Clipped Gaussian often matches experimental data better, but table generation takes more time. The choice can change peak temperatures by 50–100 K, and the experience of "casually trying different PDF shapes" becomes a difference in practical skill.
Experience the theory firsthand with the interactive simulator for this field
All Simulators