| Model | Re | Ma/Compress. | APG | Rotation | Tu | Mesh | Notes |
|---|
| Model | Type | Cost | Separation | Heat Xfr | Transition | Acoustics | Rotating | y⁺ Range | Anisotropy | Comp. Ready | Min Cells (Est.) | OF Class |
|---|
| 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. |
Φ_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.
Valid from y⁺=0 to y⁺=300 continuously. Used in OpenFOAM's nutUSpaldingWallFunction.
| Model | y+ < 1 | y+ > 30 |
|---|---|---|
| k-ε | nutLowReWallFn | nutkWallFunction |
| k-ω SST | nutLowReWallFn | nutkWallFunction |
| SA | nutLowReWallFn | nutUSpalding... |
y⁺ = u_τ·Δy / ν. First cell centroid dictates valid assumptions. Buffer layer (y⁺ 5–30) violates both viscous and log-law assumptions.
| Model | Grid (N) | Timesteps | Total Cost |
|---|---|---|---|
| RANS | Re^0 | Re^0 (steady) | O(1) |
| WMLES | Re^1.0 | Re^0.5 | O(Re^1.5) |
| WRLES | Re^1.8 | Re^0.6 | O(Re^2.4) |
| DNS | Re^2.25 | Re^0.75 | O(Re^3.0) |
Uses empirical correlations (e.g. Mayle) linking local Tu to transition onset Re_θt. γ (intermittency) equation multiplies TKE production, switching it on smoothly.
Density-weighted averaging eliminates complex density fluctuation terms in N-S equations.
As Mach increases, compressibility drains energy via shocks/dilatation. Reduces spreading rate of supersonic jets (comp. shear layer anomaly). Corrected via Sarkar/Wilcox modifiers.
The single most authoritative public database for RANS model verification. Contains grid-convergence studies, structured solution files, and experimental comparisons for: Flat Plate (ZPG), Backward Facing Step, RAE2822 Airfoil, NACA 0012, 2D Bump, CRM Wing-Body, Axisymmetric Bump, Juncture Flow, and more.
turbmodels.larc.nasa.gov — free access, updated actively.
European Research Community on Flow, Turbulence And Combustion. 100+ experimental datasets covering: pipe flow, channel flow, backward-facing step, bluff body wakes, free jets, swirling flows, impinging jets, flow over hills (ERCOFTAC SPH11). Each case has structured data files for direct comparison.
cfd.mace.manchester.ac.uk/ercoftac
Structured Verification and Validation test cases for the WIND-US solver; applicable to any CFD code. Includes flat plate, RAE2822, axisymmetric bump, and more.
www.grc.nasa.gov/www/wind/valid
Direct access to DNS datasets (channel flow Re_τ=1000, 5200; isotropic turbulence; MHD; rotating stratified). Petascale data, query-able via API or Python client. Invaluable for SGS model calibration.
turbulence.idies.jhu.edu
Open repository of spectral DNS channel flow data at Re_τ = 180–5200. Fully resolved statistics available: mean profiles, Reynolds stresses, spectra, budgets.
KTH DNS Database
AGARD AR-345 (1998) — "A Selection of Test Cases for the Validation of Large-Eddy Simulations of Turbulent Flows". AGARD AR-291 — Turbulent Shear Flows. Both available via DTIC (US Department of Defense Technical Information Center).
apps.dtic.mil/sti/
| Model | Best Validated Case | Error vs Experiment | Reference Experiment | Known Failure Flow | Peer-Reviewed Source |
|---|---|---|---|---|---|
| k-ω SST | RAE2822 Case 9 (transonic) | Cp within 5%, skin friction ±8% | Cook et al. 1979, ARC R&M 3744 | Massive post-stall separation | Menter (1994) |
| k-ε Standard | Plane mixing layer, pipe flow | Velocity profile ±3% | Rodi (1972) mixing layer | Round jet (15% error in spreading rate) | Launder & Spalding (1974) |
| k-ε Realizable | Round jet, swirling flow | Spreading rate correct | Wygnanski & Fiedler (1969) | APG separation onset | Shih et al. (1995) |
| Spalart-Allmaras | NACA 0012, flat plate BL | Cf within 2% (ZPG) | Gregory & O'Reilly (1970) | Separated free shear, heat transfer | Spalart & Allmaras (1992) |
| v2f | Impinging jet heat transfer | Nu distribution ±5% | Baughn & Shimizu (1989) | Elliptic f solver divergence on coarse mesh | Durbin (1991) |
| γ-Reθ SST | T3A/T3B/T3C flat plates | Transition location ±5–10% | Roach & Brierley (1992), ERCOFTAC | Tu decay between inlet and LE | Langtry & Menter (2009) |
| LRR RSM | Square duct secondary flow | Velocity components ±8% | Gessner & Jones (1965) | Divergence from poor initialisation | Launder et al. (1975) |
| SST-DDES | Tandem cylinder | CL frequency ±3%, PSD shape | Jenkins et al. (2005) NASA TM-2005-213649 | Grey area MSD in thin BL | Spalart et al. (2006) |
| WALE LES | Channel flow Re_τ=395 | u+ within 1%, TKE ±5% | Moser, Kim & Mansour (1999) DNS | HPC required for Re_τ>1000 | Nicoud & Ducros (1999) |
| Smagorinsky LES | Decaying isotropic turbulence | E(k) spectrum shape matches | Comte-Bellot & Corrsin (1971) | Over-dissipative in laminar near-wall regions | Smagorinsky (1963) |
| IDDES | Channel flow Re_τ=395, 2000 | Log-layer mismatch <5% | DNS Moser et al. (1999) | Highly mesh-design dependent | Shur et al. (2008) |
The Grid Convergence Index (GCI) is the AIAA/ASME standard for reporting discretization error:
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.
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.
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.
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.
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.
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.
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.
| Solver | Regime | Compatible Models |
|---|---|---|
| simpleFoam | Incompressible RANS steady | SA, k-ε, k-ω SST, RSM |
| pimpleFoam | Incompressible RANS/DDES transient | SST, DDES, SAS, IDDES |
| pisoFoam | Incompressible LES | Smag, WALE, dyn-k, Sigma, Vreman |
| rhoPimpleFoam | Compressible transient | k-ω SST, SA, DDES |
| rhoSimpleFoam | Compressible steady | k-ω SST, k-ε, SA |
| rhoCentralFoam | Supersonic/hypersonic | k-ω SST + Wilcox compressibility |
| buoyantSimpleFoam | Buoyant RANS steady | k-ε, k-ω SST (add buoyancy source) |