Turbulence Model Selection

Flow & Geometry
Simulation Goals
Configure inputs and press ▶ RUN SELECTOR.
Y-Axis Grouping
Highlight Type
Cursor Inspector (Click map to pin)
Hover map to see models at Re. Click to pin.
Regime Boundaries
Pipe laminar→turbulentRe ≈ 2300
Flat plate transition onsetRe_x ≈ 3×10⁵
Flat plate fully turbulentRe_x ≈ 5×10⁵
Cylinder drag crisisRe ≈ 2–5×10⁵
LES affordable (research)Re_L ≲ 10⁵
LES requires HPCRe_L ~ 10⁶
WMLES required (wall BL)Re_L ≳ 5×10⁶
DDES preferred range10⁵–10⁸
Click any model card to expand full profile. Includes exact equations, failure modes, validation history, and OpenFOAM tips.
Model Re Ma/Compress. APG Rotation Tu Mesh Notes
How to Read
Sensitivity Scale
0 — NoneModel is insensitive to this parameter
1 — LowMinor effect, manageable
2 — ModerateNoticeable; review setup
3 — HighSignificant; requires care
4 — CriticalDominant; results strongly affected
Abbreviations & Full Forms
ReReynolds Number (ρUL/μ)
MaMach Number (U/c)
APGAdverse Pressure Gradient
FPGFavourable Pressure Grad.
ZPGZero Pressure Gradient
TuTurbulence Intensity (%)
y⁺Wall distance (u_τ·Δy/ν)
BLBoundary Layer
TKE / kTurbulent Kinetic Energy
ωSpecific dissipation rate
εTurbulent dissipation rate
RANSReynolds-Averaged N-S
LESLarge Eddy Simulation
DESDetached Eddy Simulation
DDESDelayed DES
IDDESImproved DDES
RSMReynolds Stress Model
SGSSubgrid-Scale (model)
SASScale Adaptive Simulation
SASpalart-Allmaras
SSTShear Stress Transport
WMLESWall-Modelled LES
Model Type Cost Separation Heat Xfr Transition Acoustics Rotating y⁺ Range Anisotropy Comp. Ready Min Cells (Est.) OF Class
Head-to-Head: k-ω SST vs Competitors
Canonical Flow SST Result Quality Competitor Model & Quality Winner Engineering Note
Flat Plate (ZPG) Excellent (captures log-law perfectly) SA: Excellent (designed for it) Tie SA is 20% cheaper computationally.
RAE2822 Airfoil (Transonic APG) Excellent (predicts shock location well) k-ε: Poor (shock too far aft, no sep) SST SST limits shear stress in APG via Bradshaw's assumption.
Backward Facing Step (Sep.) Good (reattachment slightly late) v2f / LRR: Excellent (exact reattachment) RSM / v2f SST overpredicts TKE in shear layer. RSM handles anisotropy.
Pipe Flow (Internal ZPG) Good (requires fine mesh near wall) k-ε Realizable: Excellent + cheap k-ε Real. k-ε with wall functions is vastly cheaper for bulk pipe flows.
Tandem Cylinder (Acoustics) Poor (Steady RANS damps shedding) SST-DDES: Excellent (resolves wake) DDES RANS cannot resolve broadband noise. Scale-resolving required.
Turbine Cascade (Curvature) Moderate (blind to streamline curvature) SARC / SSG: Good (captures curvature) SARC / RSM Standard eddy viscosity assumes turbulence is blind to rotation.
1. RANS Closure & Boussinesq
Boussinesq Approximation
τ_ij = 2μ_t S_ij − (2/3)ρk δ_ij S_ij = 0.5·(∂U_i/∂x_j + ∂U_j/∂x_i)
k-ω SST Exact Transport Equations
∂(ρk)/∂t + ∂(ρU_j k)/∂x_j = P_k - β*ρkω + ∂/∂x_j[(μ + σ_k μ_t) ∂k/∂x_j] ∂(ρω)/∂t + ∂(ρU_j ω)/∂x_j = (γ/ν_t)P_k - βρω² + ∂/∂x_j[(μ + σ_ω μ_t) ∂ω/∂x_j] + 2(1-F1)ρσ_{ω2}(1/ω) ∂k/∂x_j ∂ω/∂x_j
RSM (Reynolds Stress Model)
∂(ρR_ij)/∂t + C_ij = P_ij + D_ij + Φ_ij - ε_ij

Φ_ij (Pressure-Strain) is the critical closure. LRR uses a linear model; SSG uses a quadratic model (better for swirl). Realizability (Lumley triangle) ensures physical positive-definite turbulent normal stresses.

2. DES / DDES Formulation
DES Length Scale
d̃ = min(d_wall, C_DES · Δ) Δ = max(Δx, Δy, Δz)
DDES Shielding Function
d̃ = d_wall − F_d · max(0, d_wall − C_DES·Δ) F_d = 1 − tanh((8·r_d)³) r_d = (ν_t + ν) / (κ² · d_wall² · √|S_ij|²)
IDDES & SAS & WMLES
  • IDDES: Blends DDES shielding with WMLES formulation. Introduces elevating function to fix log-layer mismatch.
  • SAS-SST: Uses von Kármán length scale L_vK = κS / |∇²U|. Break RANS into LES structures without explicit grid dependence.
  • WMLES: Algebraic wall model imposes wall shear stress τ_w as a BC for the resolved outer layer.
3. Wall Treatment & y⁺ Physics
Spalding Unified Wall Law
y⁺ = u⁺ + e^{-κB} [ e^{κu⁺} - 1 - κu⁺ - 0.5(κu⁺)² - (1/6)(κu⁺)³ ]

Valid from y⁺=0 to y⁺=300 continuously. Used in OpenFOAM's nutUSpaldingWallFunction.

OpenFOAM Wall Function Mapping
Modely+ < 1y+ > 30
k-εnutLowReWallFnnutkWallFunction
k-ω SSTnutLowReWallFnnutkWallFunction
SAnutLowReWallFnnutUSpalding...
y⁺ Physics

y⁺ = u_τ·Δy / ν. First cell centroid dictates valid assumptions. Buffer layer (y⁺ 5–30) violates both viscous and log-law assumptions.

4. LES Cost & Scales
Cost Scaling vs Re
ModelGrid (N)TimestepsTotal Cost
RANSRe^0Re^0 (steady)O(1)
WMLESRe^1.0Re^0.5O(Re^1.5)
WRLESRe^1.8Re^0.6O(Re^2.4)
DNSRe^2.25Re^0.75O(Re^3.0)
Turbulence Scales
  • Integral (L): Macro scale, contains most energy.
  • Taylor (λ): Inertial subrange. λ/L ~ Re_L^(-1/2)
  • Kolmogorov (η): Dissipation scale. η = (ν³/ε)^(1/4). η/L ~ Re_L^(-3/4). Resolved by DNS.
5. Transition Physics
Transition Mechanisms
  • Natural: TS waves → secondary instability → breakdown. Needs very low freestream Tu (< 0.1%).
  • Bypass: High freestream Tu (> 1%) penetrates BL, bypassing TS waves. Typical in turbomachinery.
  • Separation-induced: Laminar BL separates in APG, transitions in shear layer, reattaches as turbulent (LSB).
γ-Reθ Formulation

Uses empirical correlations (e.g. Mayle) linking local Tu to transition onset Re_θt. γ (intermittency) equation multiplies TKE production, switching it on smoothly.

6. Compressible Turbulence
Favre Averaging
ũ_i = \bar{ρ u_i} / \bar{ρ} u_i = ũ_i + u_i''

Density-weighted averaging eliminates complex density fluctuation terms in N-S equations.

Dilatation Dissipation (Sarkar)
ε_c = α M_t² ε_s M_t = √(2k) / c

As Mach increases, compressibility drains energy via shocks/dilatation. Reduces spreading rate of supersonic jets (comp. shear layer anomaly). Corrected via Sarkar/Wilcox modifiers.

Curated validation resources with direct links to databases, landmark papers (DOI), benchmark datasets, and experimental data. Each canonical flow maps to specific model recommendations.
🗄 Official Validation Databases
📄 Canonical Test Cases — Model Map
Flat Plate — ZPG Boundary Layer
ExperimentWieghardt (1951), Klebanoff (1955)
Best RANSk-ω SST, SA
Backward Facing Step (BFS)
ExperimentDriver & Seegmiller (1985), Re_h=37,400
Best RANSv2f, SST (reattachment ±5%)
Data sourceNASA TMR BFS
RAE 2822 Airfoil (Transonic)
ExperimentCook, McDonald & Firmin (1979), ARC R&M 3744
Best RANSk-ω SST (Case 9: Ma=0.73)
Data sourceNASA TMR RAE2822
2D NACA 0012 Airfoil
ExperimentGregory & O'Reilly (1970)
Best RANSSA, SST (pre-stall)
Periodic Hill (Wavy Channel)
LES referenceBreuer et al. (2009), Re_b=10595
Best approachDDES, LES (SST overestimates sep)
Data sourceERCOFTAC C69
ERCOFTAC T3 Series — Transition
ExperimentRoach & Brierley (1992)
Best RANSγ-Reθ, γ SST
NASA Common Research Model (CRM)
ExperimentVassberg et al. (2008), AIAA 2008-6919
Best RANSk-ω SST, SA-RC
📚 Landmark Papers — Full DOI References
📖 Textbooks & Review Articles
✅ Per-Model Validation Evidence Summary
Model Best Validated Case Error vs Experiment Reference Experiment Known Failure Flow Peer-Reviewed Source
k-ω SSTRAE2822 Case 9 (transonic)Cp within 5%, skin friction ±8%Cook et al. 1979, ARC R&M 3744Massive post-stall separationMenter (1994)
k-ε StandardPlane mixing layer, pipe flowVelocity profile ±3%Rodi (1972) mixing layerRound jet (15% error in spreading rate)Launder & Spalding (1974)
k-ε RealizableRound jet, swirling flowSpreading rate correctWygnanski & Fiedler (1969)APG separation onsetShih et al. (1995)
Spalart-AllmarasNACA 0012, flat plate BLCf within 2% (ZPG)Gregory & O'Reilly (1970)Separated free shear, heat transferSpalart & Allmaras (1992)
v2fImpinging jet heat transferNu distribution ±5%Baughn & Shimizu (1989)Elliptic f solver divergence on coarse meshDurbin (1991)
γ-Reθ SSTT3A/T3B/T3C flat platesTransition location ±5–10%Roach & Brierley (1992), ERCOFTACTu decay between inlet and LELangtry & Menter (2009)
LRR RSMSquare duct secondary flowVelocity components ±8%Gessner & Jones (1965)Divergence from poor initialisationLaunder et al. (1975)
SST-DDESTandem cylinderCL frequency ±3%, PSD shapeJenkins et al. (2005) NASA TM-2005-213649Grey area MSD in thin BLSpalart et al. (2006)
WALE LESChannel flow Re_τ=395u+ within 1%, TKE ±5%Moser, Kim & Mansour (1999) DNSHPC required for Re_τ>1000Nicoud & Ducros (1999)
Smagorinsky LESDecaying isotropic turbulenceE(k) spectrum shape matchesComte-Bellot & Corrsin (1971)Over-dissipative in laminar near-wall regionsSmagorinsky (1963)
IDDESChannel flow Re_τ=395, 2000Log-layer mismatch <5%DNS Moser et al. (1999)Highly mesh-design dependentShur et al. (2008)
📐 Grid Independence & V&V Best Practice
Richardson Extrapolation (Roache GCI)

The Grid Convergence Index (GCI) is the AIAA/ASME standard for reporting discretization error:

GCI = Fs × |ε| / (r^p - 1) ε = (f_2 - f_1) / f_1 (relative change) r = h_coarse / h_fine (refinement ratio ≥ 1.3) p = ln(|ε_32/ε_21|) / ln(r) (observed order) Fs = 1.25 (safety factor) Target: GCI_fine < 5% for engineering, <1% for research.
Three-Mesh GCI Protocol
  • Coarse: baseline (~100k cells for 3D RANS)
  • Medium: ×2–3 refinement by volume
  • Fine: ×2–3 additional refinement
  • Report GCI for key QoI (Cd, Nu, Cp)
  • For LES: report SGS contribution fraction <20%
y+ Verification
  • Post-run: check y+ histogram in ParaView / CFD-Post
  • For wall-resolved: 95% cells must have y+ < 1
  • For wall functions: 70% cells must have 30 < y+ < 300
  • Buffer layer contamination (5 < y+ < 30) invalidates both approaches
⚠ Common Validation Pitfalls
Inlet Condition Mismatch

If freestream turbulence intensity (Tu) and length scale (L) at the leading edge do not match the experimental conditions, all transition and separation predictions will be systematically off. The k and ω fields decay with distance — run a separate decay study to match Tu at the measurement plane, not just at the inlet boundary.

Reference Area / Length Confusion

Force coefficients (Cd, Cl) are non-dimensionalised by a reference area and dynamic pressure. Validate that your reference area, reference velocity, and density match the experimental setup exactly. A 5% mismatch in reference chord length = 5% systematic error in all force coefficients.

Incorrect Experimental Re

Many classic experiments (e.g. RAE2822) were conducted at fixed chord Reynolds number, not a body Reynolds number. Ensure your CFD inlet conditions reproduce the same Re_c, including the correct viscosity at the experimental temperature.

Statistical Convergence for URANS/LES

LES and DDES require sufficient flow-through times before statistics are sampled. Minimum: 5–10 flow-through times for wash-out, then sample for at least 5 additional flow-through times. For turbomachinery: minimum 10 blade passages. Insufficient averaging produces spurious results.

Compressibility Mismatch

At Ma > 0.3, use a compressible solver (rhoPimpleFoam, rhoCentralFoam). Incompressible simulations at Ma = 0.5 introduce >10% error in pressure coefficient even before shocks appear.

Wall Function Buffer Zone

The buffer layer (5 < y+ < 30) is in-between the viscous sublayer and log-law. Both wall function assumptions fail here. If y+ values cluster in this range, the result is neither resolved nor accurately modelled. Avoid by design: target y+ < 1 or y+ > 30, not both on the same mesh.

Disclaimer

This tool compiles information from the above peer-reviewed sources and publicly available databases. All recommendations should be validated against case-specific benchmark data before production use. Turbulence model selection is inherently problem-dependent — no single model is universally optimal. The author assumes no liability for engineering decisions made solely on the basis of this tool.

Curated YouTube lectures, CFD blog posts, open-source tutorials, and community tools for turbulence modeling practitioners.
▶ YouTube — Lectures & Tutorials
📝 Blogs, Forums & Online Resources
🔧 Tools & Calculators
📊 Decision Framework — Quick Reference
Model Selection Flowchart (Text)
  • Attached BL, external aero, cruise: → Spalart-Allmaras or k-ω SST. If Ma > 0.3, enable compressibility correction.
  • Mild separation, APG, force coefficient: → k-ω SST. If transition important → γ-Reθ SST.
  • Pipe flow, HVAC, internal duct: → k-ε Realizable with wall functions (y+ 30–100).
  • Swirling flow, cyclone, combustion: → k-ε RNG or RSM (LRR/SSG).
  • Heat transfer, conjugate: → Low-Re k-ε LS, v2f, or k-ω SST with y+ < 1.
  • Laminar–turbulent transition: → γ-Reθ SST (4-eq) or γ SST (3-eq, local).
  • Bluff body, massive separation, shedding: → SST-DDES (transient). Budget permitting: SST-IDDES.
  • Aeroacoustics: → IDDES or WMLES + FW-H acoustic analogy.
  • Rotating machinery, turbomachinery: → k-ω SST + MRF (steady) or sliding mesh (unsteady). For curvature: SA-RC.
  • Fundamental research / benchmark data: → WALE or Sigma LES. DNS for Re < 10⁵.
Key OpenFOAM Solver-Model Compatibility
SolverRegimeCompatible Models
simpleFoamIncompressible RANS steadySA, k-ε, k-ω SST, RSM
pimpleFoamIncompressible RANS/DDES transientSST, DDES, SAS, IDDES
pisoFoamIncompressible LESSmag, WALE, dyn-k, Sigma, Vreman
rhoPimpleFoamCompressible transientk-ω SST, SA, DDES
rhoSimpleFoamCompressible steadyk-ω SST, k-ε, SA
rhoCentralFoamSupersonic/hypersonick-ω SST + Wilcox compressibility
buoyantSimpleFoamBuoyant RANS steadyk-ε, k-ω SST (add buoyancy source)