Transition Model (γ-Reθ)
Transition Model (γ-Reθ): Theoretical Foundations
Overview
Professor, how is the transition from laminar to turbulent flow handled in CFD? With k-epsilon or SST k-omega, it's fully turbulent everywhere, right?
Exactly. Standard RANS models assume fully turbulent flow, so they overpredict drag in laminar regions and cannot reproduce the transition location. That's why the $\gamma$-$Re_\theta$ transition model (Langtry-Menter, 2009) was developed. It adds two transport equations for the intermittency $\gamma$ and the transition momentum thickness Reynolds number $\widetilde{Re}_{\theta t}$ to the SST k-omega model.
What is intermittency?
$\gamma$ takes a value from 0 to 1, where $\gamma = 0$ represents fully laminar flow and $\gamma = 1$ represents fully turbulent flow. In the transition region, $0 < \gamma < 1$, representing a state where laminar and turbulent flows are locally mixed.
Governing Equations
Please tell me the specific equations.
The two additional transport equations are as follows.
Transport equation for intermittency $\gamma$:
Transport equation for transition momentum thickness Reynolds number $\widetilde{Re}_{\theta t}$:
How does $\gamma$ affect the turbulence model?
$\gamma$ modifies the production term in the $k$ equation of the SST k-omega model.
In laminar regions ($\gamma = 0$), turbulence production becomes zero, keeping $k$ at a low value. This suppresses eddy viscosity, resulting in the reproduction of a laminar-like velocity profile.
Physical Mechanisms of Transition
What specific types of transition are there?
The main transition mechanisms are as follows.
| Transition Type | Physical Mechanism | Typical Scenario |
|---|---|---|
| Natural Transition | Growth of Tollmien-Schlichting waves | Airfoil surfaces with low turbulence intensity |
| Bypass Transition | High freestream turbulence directly disturbs the boundary layer | Gas turbine blade rows |
| Separation-Induced Transition | Transition occurs within a laminar separation bubble | Low-Re airfoils, compressor blades |
| Crossflow Instability Transition | Instability of the crossflow velocity profile in a 3D boundary layer | Leading edge of swept wings |
The $\gamma$-$Re_\theta$ model mainly covers natural and bypass transitions. It also handles separation-induced transition to some extent, but additional correlation functions are needed for crossflow instability transition.
"The Moment It Becomes Turbulent" — A 100-Year Question in Transition Research
The "Transition" from laminar to turbulent flow is one of the oldest research themes among the unsolved problems in fluid mechanics. Starting from Reynolds' experiments in the 19th century, followed by the discovery of Tollmien–Schlichting waves (1929) and the recognition of bypass transition (1960s), there is still no universal model to predict "under what conditions and from where turbulence begins." The γ-Reθ model is one of the most robust engineering solutions within these constraints, a practical compromise that "forces transition into the framework of turbulence models" using four equations.
Computational Methods for Transition Model (γ-Reθ)
Transition Correlation Functions
How is it determined when transition begins?
The core of the $\gamma$-$Re_\theta$ model is the empirical correlation function. It determines the critical Reynolds number $Re_{\theta c}$ from the freestream turbulence intensity $Tu$ and the pressure gradient parameter $\lambda_\theta$.
This correlation is based on experimental correlations like Abu-Ghannam-Shaw (1980) and Mayle (1991). The specific functional form is published in Menter et al.'s papers.
| Parameter | Effect |
|---|---|
| Freestream Turbulence Intensity $Tu$ | Higher values move transition upstream ($Re_{\theta c}$ decreases) |
| Pressure Gradient $\lambda_\theta$ | Favorable pressure gradient delays transition, adverse pressure gradient promotes it |
Role of the $\widetilde{Re}_{\theta t}$ Equation
What is the purpose of the $\widetilde{Re}_{\theta t}$ equation? Isn't $\gamma$ alone enough?
Good question. The transition correlation function requires the freestream $Tu$ and $\lambda_\theta$ as inputs. However, in CFD, it is technically difficult to refer to freestream values near the wall (non-local information is needed).
The transport equation for $\widetilde{Re}_{\theta t}$ is a mechanism to solve this non-locality problem. It calculates the correlation function value in the freestream and diffuses it to the near-wall region via the transport equation. This allows access to freestream Tu information even near the wall.
So $\widetilde{Re}_{\theta t}$ is like a "postman" delivering freestream information to the wall?
Exactly that image.
Coupling with SST k-omega
In terms of implementation, is it added to the SST k-omega model?
Exactly. Transition SST ($\gamma$-$Re_\theta$ SST) solves the following four equations simultaneously.
1. $k$ equation (SST k-omega, but production term modified by $\gamma$)
2. $\omega$ equation (SST k-omega, unchanged)
3. $\gamma$ equation (intermittency)
4. $\widetilde{Re}_{\theta t}$ equation (transition Re number)
Mesh Requirements
Does the transition model require a finer mesh than standard RANS?
Higher resolution than standard RANS is needed to correctly capture the transition region.
| Parameter | SST k-omega (Wall Function) | Transition SST |
|---|---|---|
| Wall $y^+$ | 30~100 | 1 or less |
| Layers within boundary layer | 8~15 | 20~30 |
| Streamwise (chordwise) | 100 cells/chord | 200~300 cells/chord |
| Spanwise | -- | Finer if solving 3D structure of transition front |
$y^+ \approx 1$ is mandatory because transition is a subtle phenomenon near the wall, right?
γ-Rθ Transition Model—"RANS Prediction of Transition" Made Practical by Menter (2006)
The Menter-Langtry γ-Reθ transition model (2006, 2009), which predicts the "transition" from laminar to turbulent flow using RANS, is the most widely used transition model in engineering. It adds two transport equations for the transition onset momentum thickness Reynolds number Reθt and the intermittency coefficient γ to the k-ω SST, handling three modes—natural transition, bypass transition, and separation-induced transition—within a single framework. It is implemented as standard in Fluent and OpenFOAM and is used for transition prediction in aircraft wings, gas turbine blades, and wind turbine blades. However, it has high sensitivity to the settings of turbulence intensity (Tu) and length scale (Λ), and the handling of inlet conditions greatly affects accuracy.
Experience the theory firsthand with the interactive simulator for this field
All Simulators