Y⁺ Calculator

UNITS:

Flow Parameters

Density ρ kg/m³

Dynamic Viscosity μ Pa·s

Freestream Velocity U∞ m/s

Geometry

Characteristic Length L m

Target y⁺ –

Growth Ratio r –

❯_Enter parameters and press Calculate

Flow Parameters

Density ρ kg/m³

Dynamic Viscosity μ Pa·s

Freestream Velocity U∞ m/s

Geometry

Characteristic Length L m

Target y⁺ –

Growth Ratio r –

❯_Enter parameters and press Calculate

Grid Convergence Index — ASME V&V 20

ℹ Enter results from three systematically refined meshes (coarse → medium → fine). The monitored variable φ should be a key flow quantity (Cf, pressure, velocity, etc.).

Refinement Ratio r₁₂ h₂/h₁

Refinement Ratio r₂₃ h₂/h₃

ϕ₁ — Coarse Mesh Value –

ϕ₂ — Medium Mesh Value –

ϕ₃ — Fine Mesh Value –

Monitored Variable Name optional

Theoretical Order p₀ –

❯_Enter three mesh values and press Calculate

Flow & Mesh Inputs

ℹ Enter your flow velocity, minimum cell size, and solver type. The tool outputs required Δt, flow-through time, CFL, and acoustic CFL for compressible solvers.

Flow Velocity U m/s

Min. Cell Size Δx m

Domain Length L m

Solver Type

Flow Regime

Mach Number Ma for compressible

Target CFL –

Iterations / Timestep for CPU estimate

Total Cells (millions) for CPU estimate

❯_Enter parameters and press Calculate

Geometry & Domain Inputs

ℹ Validates domain sizing against AIAA/NASA external aero guidelines. Computes blockage ratio, recommends upstream/downstream distances, and estimates domain cell count.

Simulation Type

Body Reference Length L m

Body Height / Span H m

Body Width / Depth W m

Domain Length X m

Domain Height Y m

Domain Width Z m

Inlet from body leading edge m

Target bulk cell size m

❯_Enter geometry and press Calculate

Flow & Surface Inputs

ℹ Estimates laminar-to-turbulent transition location using Michel (1951), Abu-Ghannam & Shaw (1980), and Menter γ-Reθ (2009) criteria. Used for airfoil, wing, body, and internal duct flows.

Geometry

Density ρ kg/m³

Dynamic Viscosity μ Pa·s

Freestream Velocity U∞ m/s

Length / Chord L m

Freestream Turbulence Intensity Tu %

Pressure Gradient

Pressure Gradient Parameter Λ for non-ZPG

❯_Enter parameters and press Calculate

Flow & Mesh Parameters

ℹ CFL = U·Δt/Δx. For explicit solvers CFL ≤ 1 is mandatory. For implicit (SIMPLE/Coupled) CFL can be larger but affects convergence rate.

Freestream / Bulk Velocity U m/s

Minimum Cell Size Δx mm

Target CFL Number –

Solver Type

Flow Regime

Domain Length L m

Speed of Sound a m/s

Iterations per Timestep –

Total Cells (millions) –

❯_Enter parameters and press Calculate

Body Dimensions

ℹ Blockage ratio = A_body / A_domain × 100%. External aero: < 3%. Wind-tunnel simulation: < 5%. Bluff bodies: < 1%.

Flow Configuration

Chord / Length C m

Span / Width W m

Max Thickness / Height t m

Domain Dimensions

Domain Length (streamwise) m

Domain Height (cross-stream) m

Domain Width (spanwise) m

Inlet Distance (body LE to inlet) m

❯_Enter body and domain dimensions and press Calculate

Flow Parameters

ℹ Estimates laminar-to-turbulent transition location using Michel (1952), Abu-Ghannam & Shaw (1980), and Granville (1953) criteria. Valid for attached 2D boundary layers with low-to-moderate freestream turbulence.

Density ρ kg/m³

Dynamic Viscosity μ Pa·s

Freestream Velocity U∞ m/s

Surface Length L m

Freestream Turbulence Intensity Tu %

Pressure Gradient Character

Geometry / Application

❯_Enter parameters and press Calculate

CFL Courant-Friedrichs-Lewy Condition

CFL Number — Courant, Friedrichs & Lewy (1928)

CFL= U · Δt / Δx (convective)

CFL_a= (U + a) · Δt / Δx (acoustic, compressible)

Δt= CFL · Δx / U (required timestep)

Ref: Courant R., Friedrichs K. & Lewy H. (1928), Math. Ann. 100, 32–74. Stability condition: explicit solvers require CFL ≤ 1. Implicit solvers allow CFL > 1 but convergence per step degrades. Acoustic CFL must also be ≤ 1 for compressible explicit solvers.

Flow-Through Time & Simulation Duration

T_ft= L / U [s] (one convective pass through domain)

N_steps= T_total / Δt (total timesteps)

T_total= 5 · T_ft (minimum for turbulent statistics)

For LES and time-averaged RANS, accumulate statistics over at least 5 flow-through times after initial transients have cleared (typically 2–3 T_ft for washout). Ref: Pope, S.B. (2000), Turbulent Flows, Cambridge.

CFL Limits by Solver Class

ExplicitCFL ≤ 1.0 (hard stability limit)

SIMPLE/SIMPLECCFL typically 5–50

Coupled PB implicitCFL 50–200

Density-based implicitCFL 5–50

Local time stepping (pseudo-transient) for steady-state: start at CFL = 1, increase to 100+ once residuals stabilise. OpenFOAM: set maxCo in controlDict; Fluent: set Courant Number in solver controls.

DOMAIN Blockage & Domain Sizing

Solid Blockage Ratio — Maskell (1965)

BR= A_body / A_domain × 100 [%]

BR (2D)= t / H_domain × 100 (cross-section ratio)

BR (3D)= (W · t) / (H · D_domain) × 100

Limits: External aero: BR < 3% (NASA/AIAA standard). Wind-tunnel simulation: BR < 5%. Bluff bodies: BR < 1% recommended. Above 5%, Maskell solid blockage correction factor C_B = 1 + K·BR applies to Cd measurements. Ref: Maskell E.C. (1965), RAE TR 65-237; Barlow, Rae & Pope (1999) Low Speed Wind Tunnel Testing, Wiley.

Domain Extent Guidelines (AIAA/NASA)

Upstream≥ 5C (airfoil); ≥ 10C (bluff body)

Downstream≥ 15C (airfoil); ≥ 20C (bluff body, wake recovery)

Lateral≥ 5C each side (external aero)

C = characteristic length (chord, diameter, or body length). Domain must be large enough for outlet/side boundary conditions not to influence the solution near the body. Pressure outlet: reflection-free if placed ≥ 15C downstream. Ref: AIAA CFD Validation Workshop guidelines; NASA CFL3D User Manual.

Blockage Correction — Maskell (1965)

U_corr= U_meas · (1 + C_B · BR)

Cd_corr= Cd_meas / (1 + 2·BR/100)²

Applied when comparing CFD (or wind tunnel) results at BR > 3% to free-air data. ESDU 98007 provides more comprehensive corrections. C_B ≈ 0.96 for bluff bodies; 0.52 for streamlined bodies.

TRANSITION Laminar-Turbulent Transition Criteria

Michel Criterion (1952) — ZPG, Low Tu

Re_θ_tr≥ 1.174 · (1 + 22400/Re_x)^0.46 · Re_x^0.46

Re_θ(x)= 0.6641 · √Re_x (Blasius laminar)

Transition occurs at the x-location where the Blasius Re_θ curve crosses the Michel criterion curve. Valid for ZPG, Tu < 0.5%, smooth wall, attached flow. Ref: Michel R. (1952), ONERA NT 1/16.

Abu-Ghannam & Shaw (1980) — Tu-modified

Re_θ_onset= 163 + exp(6.91 + 12.986λ − Tu·F/A)

λThwaites' pressure gradient parameter (ZPG: λ = 0)

Accounts for freestream turbulence intensity Tu (%) and pressure gradient. FPG (λ > 0) stabilises BL (raises Re_θ_onset); APG (λ < 0) promotes early transition. Ref: Abu-Ghannam B.J. & Shaw R. (1980), J. Mech. Eng. Sci. 22(5), 213–228.

Granville / Mayle Bypass Correlation (1991)

Re_x_tr= (52 / Tu%)² (Tu in percent)

Bypass transition for Tu > 0.3%. Tollmien-Schlichting waves are bypassed; Klebanoff streaks trigger turbulence directly. Valid for 0.3% ≤ Tu ≤ 5%. Ref: Mayle R.E. (1991), ASME J. Turbomachinery 113(4), 509–537.

Combined Cf (Michel transition correction)

Cf= 0.074/Re^0.2 − 1742/Re (ZPG, tr at Re_x = 5×10⁵)

Accounts for laminar leading-edge portion then turbulent remainder. Overestimates Cf if used with assumed tr at 5×10⁵ when actual tr differs. Use the general form: Cf = Cf_turb − (Cf_turb_tr − Cf_lam_tr)·Re_x_tr/Re_L. Ref: Schlichting & Gersten, BL Theory 8th ed.

e^N Method — Arnal (1994)

N_tr= 9 (free flight, Tu < 0.1%)

N_tr= 7 (low-turbulence wind tunnel)

N_tr= −8.43 − 2.4·ln(Tu/100) (Mack's correlation)

The e^N method integrates the spatial growth rate of Tollmien-Schlichting waves; transition when amplification reaches e^N. N = 9 is the standard for free-flight aero (Lockheed, ISRO, HAL design practice). N decreases with Tu. Ref: Arnal D. (1994), AGARD-R-793.

γ-Reθ SST Model — Langtry & Menter (2009)

γIntermittency (γ = 0: laminar; γ = 1: fully turbulent)

Reθ_tTransition momentum-thickness Re — transport equation

Two additional transport equations for γ and Reθ_t coupled to k-ω SST. Correlates Tu and pressure gradient to transition onset. Most widely used transitional RANS model (ANSYS Fluent: Transition SST; OpenFOAM: gammaReThetatSST). Ref: Langtry R.B. & Menter F.R. (2009), AIAA J. 47(12).

CORE Universal Definitions & Derived Quantities

Reynolds Number — Osborne Reynolds (1883)

Re= ρ · U∞ · L / μ = U∞ · L / ν

Re_eff= ρ · U∞ · L_eff / μ

Ratio of inertial to viscous forces. Governs laminar/turbulent transition. Re < 2300 (pipe): laminar; Re > 4000 (pipe): turbulent. Re < 5×10⁵ (flat plate): laminar BL.

Kinematic Viscosity

ν= μ / ρ [m²/s]

Momentum diffusivity. Sets the thickness of the viscous sublayer. Standard air (20°C): ν ≈ 1.516×10⁻&sup5; m²/s.

Wall Shear Stress & Friction Velocity

τ_w= Cf · ½ · ρ · U∞² [Pa]

u_τ= √(τ_w / ρ) [m/s]

l_ν= ν / u_τ [m]

u_τ is the friction velocity — the fundamental velocity scale for near-wall turbulence. l_ν is the viscous length scale; it sets the physical thickness of the viscous sublayer (~5·l_ν). All non-dimensional wall coordinates are normalised by these two scales.

Non-dimensional Wall Distance y⁺

y⁺= ρ · u_τ · y / μ = y / l_ν

Δy= y⁺_target · μ / (ρ · u_τ) = y⁺_target · l_ν

y is the physical distance from the wall to the first cell centroid. The rearranged form gives the required first-layer height Δy for a specified y⁺_target. Note: Δy is the cell-centre distance, NOT the face distance.

Inflation Layer Count N & Total Stack Height T

N= ⌈ log[(δ·(r−1)/Δy + 1)] / log(r) ⌉

T= Δy · (r^N − 1) / (r − 1)

N is the minimum number of inflation layers for the stack to reach the BL thickness δ at growth ratio r. T should be ≥ δ. Derived from geometric series sum. Recommended r = 1.1–1.2 for accuracy; r > 1.3 introduces wall-normal numerical diffusion.

LAW OF THE WALL Near-Wall Velocity Profile

Viscous Sublayer — Prandtl (1925)

u⁺= y⁺ (valid 0 ≤ y⁺ ≤ 5)

Linear velocity profile; viscous stresses dominate; turbulent stresses negligible. Solver must resolve this region for accurate τ_w with low-Re turbulence models (k-ω SST, Launder-Sharma k-ε). First cell must sit inside this layer (y⁺ ≤ 1).

Buffer Layer

5 < y⁺ < 30Transition zone — neither law valid

Both viscous and turbulent stresses are significant. No analytical expression applies. Wall functions are inaccurate here; direct resolution is difficult without sufficient mesh density. Always avoid placing the first cell in this region.

Log-Law Region — von Kármán & Prandtl (1930s)

u⁺= (1/κ) · ln(y⁺) + B (valid 30 ≤ y⁺ ≤ 300)

κ= 0.41 (von Kármán constant)

B= 5.0 (smooth wall constant)

Logarithmic overlap region where inertial and viscous effects balance. Standard wall functions (k-ε family) exploit this region. B = 5.0 for hydraulically smooth walls; B decreases with roughness. κ = 0.41 is the experimentally determined von Kármán constant (some sources use 0.4 or 0.42).

Spalding's Composite Law — Spalding (1961)

y⁺= u⁺ + e^(−κB) · [e^(κu⁺) − 1 − κu⁺ − (κu⁺)²/2 − (κu⁺)³/6]

Single implicit equation valid across the entire inner layer (sublayer + buffer + log-law). Solved numerically for u⁺ given y⁺. Used in enhanced wall treatments and RANS models. Ref: Spalding, D.B. (1961), A single formula for the 'law of the wall', J. Appl. Mech. 28(3).

Outer Layer (Wake Region) — Coles (1956)

y⁺ > 300Outer layer; Coles wake function applies

Turbulent fluctuations dominate; viscosity negligible. Log-law deviates; Coles wake term Π/κ·W(y/δ) corrects this. Beyond the scope of standard wall functions — place first cell here only for very coarse engineering estimates.

LAMINAR Skin Friction & BL Thickness Correlations

flat_plate · airfoil · curved_plate

Blasius Flat-Plate Solution — H. Blasius (1908)

L_eff= L (plate / chord / arc length)

Cf= 1.328 / √Re_eff (average skin friction)

δ‑99= 5.0 · L_eff / √Re_eff

δ*= 1.7208 · L_eff / √Re_eff (displacement)

θ= 0.6641 · L_eff / √Re_eff (momentum)

H= δ*/θ = 2.591 (shape factor, canonical Blasius value)

Validity: Re < 5×10⁵, zero pressure gradient (ZPG), incompressible, smooth wall, BL starts at leading edge. Exact self-similar solution to the BL equations assuming Falkner-Skan β=0. Ref: Blasius, H. (1908), Grenzschichten in Flüssigkeiten mit kleiner Reibung, Z. Math. Phys. 56, 1–37.

pipe · duct

Hagen-Poiseuille Law — G.H.L. Hagen (1839) & J.L.M. Poiseuille (1840)

L_eff= D_h = 4A/P (hydraulic diameter)

f_Darcy= 64 / Re_Dh (Darcy-Weisbach friction factor)

Cf= f_Darcy / 4 (Fanning friction factor)

δ= D_h / 2 (BL grows to fill centreline)

Δp/L= f_Darcy · ρU² / (2D_h) (pressure gradient)

Validity: Re < 2300, fully developed (x > 0.06·Re·D_h), circular cross-section. Exact analytical solution; parabolic velocity profile u/U_c = 1−(r/R)². For non-circular ducts, D_h = 4A/P is an approximation; accuracy degrades at aspect ratio >4:1. Ref: Schlichting & Gersten, Boundary Layer Theory, 8th ed., Springer.

cylinder · sphere

Blasius Applied to Curved Arc Length

L_eff= π · D / 2 (stagnation-point arc to separation)

Re_eff= ρ · U∞ · L_eff / μ

Cf= 1.328 / √Re_eff

δ‑99= 5.0 · L_eff / √Re_eff

Engineering approximation only. Treats the front hemisphere arc as a flat plate of equivalent length. Ignores streamline curvature, favourable/adverse pressure gradients, and the actual separation angle (sub-critical cylinder: ~80°; super-critical: ~140°). Ref: White, F.M., Fluid Mechanics, 8th ed., McGraw-Hill.

TURBULENT Skin Friction & BL Thickness Correlations

flat_plate · airfoil · curved_plate

Prandtl-Schlichting Formula — L. Prandtl & H. Schlichting (1932/1950)

L_eff= L

Cf= 0.370 / (log₁₀ Re_eff)^2.584

δ= 0.37 · L_eff / Re_eff^0.2 (Prandtl 1/5-power law)

δ*= 0.048 · L_eff / Re_eff^0.2

θ= 0.037 · L_eff / Re_eff^0.2

H= δ*/θ ≈ 1.3 (attached turbulent BL)

Validity: 5×10⁵ < Re < 10⁹, ZPG, smooth wall, fully turbulent from leading edge, incompressible (Ma < 0.3). Derived from the 1/7 power-law velocity profile u/U∞ = (y/δ)^(1/7). Ref: Schlichting, H. & Gersten, K. (2000), Boundary Layer Theory, 8th ed., Springer. Also: cfd-online.com/Wiki/Skin_friction_coefficient.

pipe · duct — Re < 10⁵

Blasius Pipe Friction Formula — H. Blasius (1913)

f_Darcy= 0.3164 · Re_Dh^(−0.25)

Cf= f_Darcy / 4 (Fanning friction factor)

δ= D_h / 2

Validity: 4×10³ < Re < 10⁵, smooth pipe, fully developed turbulent flow. Empirical correlation with 1/4-power law. For Re > 10⁵, use Petukhov-Filonenko. Ref: Blasius, H. (1913), Das Ähnlichkeitsgesetz bei Reibungsvorgängen in Flüssigkeiten, VDI Forschungsheft 131.

pipe · duct — Re ≥ 10⁵

Petukhov-Filonenko Formula — Filonenko (1954) / Petukhov (1970)

f_Darcy= (0.790 · ln Re_Dh − 1.64)^(−2)

Cf= f_Darcy / 4

δ= D_h / 2

Validity: 10⁵ < Re < 5×10⁶, smooth pipe, fully developed turbulent flow. More accurate than Blasius at high Re. Also known as the Filonenko correlation. For rough pipes use Colebrook-White: 1/√f = −2·log(k_s/(3.7D) + 2.51/(Re·√f)). Ref: Petukhov, B.S. (1970), Advances in Heat Transfer, Vol. 6, Academic Press.

cylinder · sphere

Prandtl-Schlichting on Effective Arc Length

L_eff= π · D / 2 (stagnation to equator)

Re_eff= ρ · U∞ · L_eff / μ

Cf= 0.370 / (log₁₀ Re_eff)^2.584

δ= 0.37 · L_eff / Re_eff^0.2

Engineering approximation. Flat-plate turbulent correlation applied to the arc length from stagnation point to separation. Ignores curvature, pressure gradient effects, and the actual separation angle. Use only as first estimate for mesh sizing. Ref: White, F.M., Fluid Mechanics, 8th ed.

BL INTEGRALS Integral Thicknesses, Shape Factor & Derived Re

Displacement Thickness δ* (von Kármán, 1921)

δ*= ∫₀∞ (1 − u/U∞) dy [general definition]

δ* (turb)= 0.048 · L_eff / Re_eff^0.2 (1/7 power-law)

δ* (lam)= 1.7208 · L_eff / √Re_eff (Blasius exact)

Equivalent distance by which the external streamlines are displaced outward due to the retarded flow in the BL. Key parameter in airfoil effective camber and wind-tunnel blockage corrections.

Momentum Thickness θ (von Kármán, 1921)

θ= ∫₀∞ (u/U∞)(1 − u/U∞) dy [general]

θ (turb)= 0.037 · L_eff / Re_eff^0.2

θ (lam)= 0.6641 · L_eff / √Re_eff

Equivalent thickness of fluid with zero momentum. Directly related to drag via the von Kármán integral relation: dC_D/dx = 2(dθ/dx). Used in Michel's transition criterion: transition occurs when Re_θ ≥ 1.174(1 + 22400/Re_x)^0.46 · Re_x^0.46.

Shape Factor H

H= δ* / θ

H = 1.29Turbulent, attached (1/7 power-law)

H = 2.59Laminar, attached (Blasius)

H ≥ 2.5Turbulent separation imminent

H ≥ 3.5Laminar separation imminent

Indicator of BL health. H increases as the BL approaches separation. Rotta (1950) extended the integral method using H as a parameter. Separation occurs at H ≈ 2.5–3.5 for turbulent; H ≈ 3.5–4 for laminar.

Derived Reynolds Numbers

Re_θ= ρ · U∞ · θ / μ (momentum-thickness Re)

Re_τ= u_τ · δ / ν (friction Re)

Re_θ is the key parameter in transition prediction (e.g. e^N method, Abu-Ghannam & Shaw criterion). Re_τ governs the turbulence length-scale hierarchy — the ratio of outer to inner scales. LES becomes feasible at Re_τ ≤ ~2000; DNS requires Re_τ ≤ ~500–1000 for practical domains. For pipe flow: Re_τ ≈ 0.09·Re^0.875 (McKeon et al. 2004).

pipe · duct (turbulent — 1/7 power-law)

Fully Developed Pipe BL Integrals

δ*= D_h / 16 (= R/8)

θ= 7·D_h / 144 (= 7R/72)

H≈ 1.286 (turbulent pipe, 1/7 law)

MESHING STRATEGY y⁺ Target by Turbulence Model

Resolve Sublayer — y⁺ ≤ 1

k-ω SSTy⁺ ≤ 1 required. Menter (1994). Best for APG, separation, airfoil stall, turbomachinery.

k-ω (Wilcox)y⁺ ≤ 1. Wilcox (1988, 2006). Good free-stream sensitivity; integrates to wall.

LES (wall-resolved)y⁺ ≤ 1, Δx⁺ ≤ 100, Δz⁺ ≤ 20. Dynamic Smagorinsky or WALE model.

DNSy⁺ ≤ 1 (stricter). All length scales resolved; no model. Re_τ ≤ ~1000 practical.

Low-Re k-εy⁺ ≤ 1. Launder-Sharma (1974), Lam-Bremhorst (1981). Require viscous damping terms.

Heat transfery⁺ ≤ 1 regardless of turbulence model (Prandtl number effects near wall).

Standard Wall Functions — y⁺ = 30–300

Std k-εy⁺ = 30–100. Launder & Spalding (1974). Log-law assumed at first cell.

Realisable k-εy⁺ = 30–100. Shih et al. (1995). Better for rotating flows and strong streamline curvature.

RNG k-εy⁺ = 30–100. Yakhot & Orszag (1986). Renormalization group; good for transitional effects.

Scalable wall fn. (Fluent)Clips y⁺ to ≥ 11.225 automatically. Avoids buffer-layer placement.

Standard wall fn. (OpenFOAM)nutWallFunction, kqRWallFunction, epsilonWallFunction on wall patches.

LES (wall-modelled)y⁺ = 50–200. Piomelli & Balaras (2002). WMLES reduces cost at high Re.

Spalart-Allmaras — Flexible

SA modelSpalart & Allmaras (1992). y⁺ ≤ 1 (resolved) or y⁺ ≥ 30 (wall-fn). Avoid buffer layer.

1-equation model; robust for external aerodynamics (aircraft, missiles). Less accurate for strongly separated flows or free shear layers. Standard in aerospace industry (NASA CFL3D, DLR TAU).

Buffer Layer — Always Avoid

5 < y⁺ < 30Neither wall functions nor direct resolution is accurate here.

Placing the first cell in the buffer layer causes systematic errors in Cf, heat flux, and separation prediction regardless of turbulence model. If using wall functions, ensure y⁺_min ≥ 30 everywhere on walls; if resolving, ensure y⁺_max ≤ 1 everywhere.

GLOSSARY All Variables & Physical Significance

Input Variables

ρDensity [kg/m³] — fluid mass per unit volume. For air at 15°C, 1 atm: 1.225 kg/m³. For water at 20°C: 998 kg/m³.

μDynamic viscosity [Pa·s] — resistance to shear deformation. Air at 15°C: 1.789×10⁻&sup5; Pa·s. Temperature dependent: Sutherland's law for gases.

νKinematic viscosity [m²/s] = μ/ρ — momentum diffusivity. Sets thickness of viscous sublayer. Air: 1.46×10⁻&sup5; m²/s.

U∞Freestream / bulk velocity [m/s] — far-field speed or mean pipe velocity. Used as the reference velocity for Re and Cf.

LCharacteristic length [m] — plate length, chord, D_h, or diameter. Sets Re and BL development length.

L_effEffective length [m] — = L for flat surfaces; = πD/2 for cylinder/sphere (arc from stagnation to separation).

rInflation growth ratio [–] — each successive layer height = r × previous. Typical value 1.1–1.2. r > 1.3 increases numerical diffusion.

Wall & Friction Quantities

CfSkin friction coefficient [–] = τ_w / (½ρU∞²). Non-dimensional wall shear. Integrated along the surface gives total friction drag coefficient C_D,f = ∫Cf dx/L.

τ_wWall shear stress [Pa] — tangential traction at the wall surface. Drives near-wall production of turbulent kinetic energy. k = 0 at the wall; τ_w = μ·(du/dy)|_wall for laminar or via wall functions.

u_τFriction velocity [m/s] = √(τ_w/ρ) — fundamental inner-layer velocity scale. All turbulent wall quantities are normalised by u_τ and ν. Typical value ≈ 1–5% of U∞ for attached BL.

l_νViscous length scale [m] = ν/u_τ — thickness scale for the viscous sublayer. The sublayer spans 0–5·l_ν from the wall. For air at Re=10⁶: l_ν ≈ 5–50 µm.

Mesh Output Variables

y⁺Non-dimensional wall distance [–] = y/l_ν. The single most important mesh quality metric for wall-bounded flows. Determines which wall treatment the solver can use.

ΔyFirst cell height [m] — physical distance from wall face to cell centroid. This is what you enter in ANSYS Meshing (First Layer Thickness) or blockMesh/snappyHexMesh (firstLayerThickness).

NNumber of inflation layers [–] — minimum layers needed for stack to grow from Δy to δ at ratio r. Add 10–20% margin for actual simulations.

TTotal inflation stack height [m] — should be ≥ δ to capture the full boundary layer. T = Δy·(r^N − 1)/(r−1) from geometric series.

CFL Courant-Friedrichs-Lewy Number & Timestep

CFL Condition — Courant, Friedrichs & Lewy (1928)

CFL= U · Δt / Δx (convective; incompressible)

CFL_a= (U + a) · Δt / Δx (acoustic; compressible)

Δt= CFL · Δx / (U + a) (required timestep)

Original ref: Courant, R., Friedrichs, K. & Lewy, H. (1928), Math. Ann. 100, 32–74. The CFL number measures how many cells a wave crosses per timestep. CFL ≤ 1 is the stability limit for explicit schemes; implicit schemes are unconditionally stable but lose accuracy at high CFL. Ref: Blazek (2015) CFD: Principles and Applications, §6.3.

Solver-Specific Limits & Guidance

ExplicitCFL ≤ 1.0 strictly (forward Euler, Runge-Kutta). Violation = divergence.

PISO/PIMPLECFL ≤ 1 for time-accuracy (URANS, LES). CFL ≤ 5 if time accuracy not critical.

SIMPLEPseudo-transient steady; no strict CFL limit. Under-relaxation governs convergence.

CoupledLocal CFL = 5–100 for startup; increase to 200+ near convergence (Fluent Coupled solver).

Density-implicitCFL = 1–50. Higher CFL reduces time-accuracy. Ref: ANSYS Fluent TG §26.4.

For LES: CFL ≤ 0.5–1.0 everywhere in the domain (not just average). Use max CFL, not mean, as the limiting metric. Ref: Ferziger, Peric & Street (2020) Computational Methods for Fluid Dynamics, §7.3.

Flow-Through Time & Simulation Duration

T_flow= L / U [s] — time for fluid to traverse domain once

N_steps= T_flow / Δt [–] — timesteps per flow-through

T_sim≥ 5–10 · T_flow for statistically converged LES/URANS

For periodic or vortex-shedding flows, simulate at least 5 shedding cycles after initial transient. Strouhal number St = f·D/U gives shedding frequency f. Typical: St ≈ 0.2 for circular cylinder (Kármán vortex street). Ref: Bearman (1984) Annu. Rev. Fluid Mech. 16.

DOMAIN Domain Sizing & Blockage Ratio

Blockage Ratio — Barlow, Rae & Pope (1999)

BR= A_body / A_domain × 100 [%]

A_body= H_body × W_body (frontal projected area)

A_domain= H_domain × W_domain (cross-sectional area)

Limit< 3% streamlined; < 3% bluff body; < 5% internal

High blockage artificially accelerates flow around the body, increasing measured drag and lift. Maskell (1965) correction: C_D,corrected = C_D,measured / (1 + BR·k), k ≈ 1 for bluff bodies. Ref: Barlow, Rae & Pope (1999) Low-Speed Wind Tunnel Testing, §9.4; Maskell (1965) RAE Report 1665.

Domain Extent Guidelines — NASA/TM-2006-214293

StreamlinedUpstream: 5–10L; Downstream: 10–20L; Lateral: 5L each side

Bluff bodyUpstream: 5–10L; Downstream: 15–25L; Lateral: 5–10L each side

BuildingUpstream: 5H; Downstream: 15H; Lateral: 5H; Vertical: 5H

Internal pipeInlet: ≥ 10D_h (developed inlet) or ≥ 2D_h + profiled inlet BC

Larger domains reduce boundary condition influence but increase cell count. The outlet should be in a region of zero or near-zero velocity gradient. Pressure outlet (zero gauge) is the standard BC; avoid inlet condition (backflow) at the outlet boundary. Ref: Vassberg & Jameson (2010) AIAA-2010-4438; COST Action 732 (2007).

Cell Count Estimate

N_bulk= V_domain / (Δx_bulk)³ (uniform hex)

N_total≈ 3–5 × N_bulk (including BL inflation & surface refinement)

Actual count depends heavily on local refinement boxes, surface curvature adaptation, and inflation layers. For external aero: typical 5–50M cells for RANS; 100M–10B for wall-resolved LES. Use AMR (adaptive mesh refinement) to concentrate cells in wake and separated regions.

TRANSITION Laminar-Turbulent Transition Criteria

Michel Criterion — Michel (1951)

Re_θ_tr≥ 1.174 · (1 + 22400/Re_x)^0.46 · Re_x^0.46

Re_θ= ρ U θ / μ = 0.6641 √Re_x (Blasius)

Empirical correlation; compares local momentum-thickness Re with a critical value that depends on Re_x. Transition occurs where the left side first exceeds the right side as x increases. Accurate to ±15% for ZPG flows; less accurate with pressure gradients. Ref: Michel, R. (1951) ONERA Report 1/1578A; Schlichting & Gersten (2000) §16.2.

Abu-Ghannam & Shaw (1980) — Turbulence-Intensity-Dependent

Re_θ_onset= 163 + exp[F(Λ) − Tu]

F(Λ)= 6.91 + 12.75Λ + 63.64Λ² (Λ ≥ 0, FPG)

F(Λ)= 6.91 + 2.48Λ − 12.27Λ² (Λ < 0, APG)

Re_θ_end= 540 + exp[7.57 − 2Tu/3] (fully turbulent)

Λ= Thwaites pressure-gradient parameter = θ²/ν · dU/dx

Tu is freestream turbulence intensity in percent. Λ ≈ −0.09 at laminar separation; Λ ≈ +0.1 at strong acceleration. Accounts for bypass transition at high Tu. Ref: Abu-Ghannam, B.J. & Shaw, R. (1980) J. Mech. Eng. Sci. 22(5), 213–228.

Menter γ-Reθ Model — Langtry & Menter (2009)

Re_θ_t= 1173.51 − 589.428Tu + 0.2196/Tu² (Tu ≤ 1.3%)

Re_θ_t= 331.5 · (Tu − 0.5658)^(−0.671) (Tu > 1.3%)

γIntermittency factor: 0 (laminar) → 1 (fully turbulent)

Transport equations for intermittency γ and transition onset momentum-thickness Reθ_t. Coupled with k-ω SST. Handles natural, bypass, separation-induced, and wake-induced transition. Ref: Langtry & Menter (2009) AIAA J. 47(12), 2894–2906. Two transport equations added to SST; requires y⁺ ≤ 1.

Tollmien-Schlichting Instability — Schlichting (1933)

Re_δ*_crit≈ 520 (neutral stability, ZPG flat plate)

Re_x_crit≈ 9.1×10⁴ (TS wave onset)

TS waves are the first instability mechanism in low-Tu environments. Amplification begins at Re_x ≈ 10⁵ and grows exponentially (e^N method, N = 9 for flight; N = 6 for wind-tunnel Tu). Full transition at Re_x ≈ 3–5×10⁵ for ZPG. Ref: Schlichting & Gersten (2000) §15.3; Drela & Giles (1987) AIAA J. (ISES e^N method).

ALL CASES General Assumptions & Fundamental Limits

Flow Regime Assumptions

•Steady, time-averaged flow — all correlations are derived for time-averaged (RANS) or steady laminar conditions. Instantaneous turbulent fluctuations are not represented. Transient simulations require unsteady RANS (URANS) or LES and may need different mesh requirements at each timestep.

•Incompressible flow (Ma < 0.3) — density ρ and dynamic viscosity μ are treated as constants throughout the domain. For gases, this is valid when density variation < 5%, corresponding to Ma < 0.3. Above this, compressibility corrections are required.

•Newtonian fluid — linear constitutive relation τ = μ·(du/dy). Not applicable to non-Newtonian fluids: polymer solutions, blood, drilling mud, food products, or slurries. These require generalised Newtonian or viscoelastic models.

•Hydraulically smooth wall — all Cf correlations assume zero surface roughness (k_s⁺ < 5, where k_s⁺ = k_s·u_τ/ν). Surface roughness increases Cf substantially and shifts the log-law intercept B. For rough surfaces, the Colebrook-White equation applies to pipes; Nikuradse roughness corrections apply to flat plates.

•Single-phase flow — no multiphase effects, cavitation, free surfaces, or dissolved gases. Near-wall physics changes significantly in bubbly flow, film boiling, or cavitating flows.

Mesh Sizing Assumptions

•Δy is cell-centroid distance from the wall — not the cell face. ANSYS Fluent uses centroid distance internally. In OpenFOAM, blockMesh firstLayerThickness refers to cell height (= 2×centroid distance for uniform cells). Verify convention for your solver before applying Δy directly.

•Uniform geometric growth at constant ratio r. The tool assumes all N layers grow at exactly the same ratio from Δy. Variable-ratio or hybrid inflation stacks (common in snappyHexMesh) are not modelled.

•2D flat boundary layer assumption — the y⁺ estimate is based on average Cf over the whole surface. In 3D geometries, y⁺ will vary significantly: lower at the leading edge/inlet, higher downstream. Check y⁺ contour in post-processing after the first simulation and adjust Δy accordingly.

•First-solve adjustment — this tool provides a pre-simulation estimate. The calculated Δy assumes the Cf from the empirical correlation. After the first solver run, use the actual wall shear stress field to recompute y⁺ = ρ·u_τ·Δy/μ and refine the mesh if needed.

LAMINAR Blasius & Hagen-Poiseuille Limitations

flat_plate · airfoil · curved_plate

Blasius (1908) — Critical Limitations

✗Re ≥ 5×10⁵ — transition to turbulence begins. The Blasius solution is invalid above this threshold; turbulent Cf will be higher. Transition Re depends on freestream turbulence intensity: high Tu lowers transition Re to ~10⁵; low-turbulence environments may sustain laminar flow to Re ≈ 3×10⁶.

✗Zero pressure gradient (ZPG) only — Blasius is the exact solution for a flat plate with no external pressure gradient (d p/dx = 0). Favourable pressure gradients (accelerating flow) thin the BL and reduce Cf; adverse pressure gradients (decelerating flow) thicken it, increase Cf, and can cause laminar separation. For non-ZPG geometries, use Falkner-Skan solutions or Thwaites' method.

✗BL starts at leading edge (x = 0) — assumes a sharp leading edge with no upstream BL. Pre-existing BL from upstream geometry (e.g. wind-tunnel contraction, nozzle wall) is not accounted for. Use displacement thickness δ* to account for virtual origin shift.

✗2D, attached flow — Blasius is strictly 2D. 3D effects (crossflow, swept wings, corner vortices) are not captured. For swept wings, the Falkner-Skan-Cooke solution applies.

⚠Airfoil application — using chord length as L_eff gives an average Cf for the entire chord. The actual BL thickness varies strongly along the chord due to pressure distribution. Suction side has APG near trailing edge; pressure side may have FPG near mid-chord. Use panel methods (XFOIL) for more accurate airfoil BL estimates.

⚠Heat transfer coupling — Blasius assumes isothermal flow. Significant wall heating/cooling alters fluid viscosity and density near the wall, modifying the BL profile. Use the Eckert reference temperature method for coupled flows.

pipe · duct

Hagen-Poiseuille — Critical Limitations

✗Re ≥ 2300 — transitional flow begins. Fully turbulent pipe flow typically establishes at Re ≈ 4000. Between 2300–4000 flow is intermittently turbulent (slug flow). Hagen-Poiseuille is strictly invalid above Re = 2300.

✗Entry length not covered — Hagen-Poiseuille is the fully-developed solution valid for x > L_h = 0.06·Re·D (laminar hydrodynamic entry length). Near the pipe inlet, the BL is developing and local Cf is much higher. For x < L_h, use the Graetz entry-length solution.

✗Circular cross-section — exact solution valid only for circular pipes. Rectangular, elliptical, or annular ducts require separate solutions. The hydraulic diameter D_h = 4A/P approximation degrades at aspect ratio >4:1; for very flat channels use the exact rectangular duct solution (Shah & London 1978).

⚠Straight pipe only — bends, elbows, and T-junctions create secondary flows (Dean vortices in bends) that significantly increase friction factor. Use the Dean number correlation: f/f_straight = 1 + 0.033(log De)^4 for De > 11.6 (De = Re·√(D/R_c)).

cylinder · sphere

Curved Body Limitations

✗Arc length approximation only — L_eff = πD/2 is a crude estimate of the BL development length. The actual stagnation-to-separation arc depends on Re (separation for laminar cylinder: ~80° from stagnation; turbulent: ~140°). This is an order-of-magnitude mesh sizing estimate only.

✗Streamline curvature ignored — the flat-plate Blasius solution does not account for centrifugal effects of curved streamlines. Concave surfaces destabilise the BL (Görtler instability); convex surfaces stabilise it. Both alter Cf and δ significantly.

TURBULENT Prandtl-Schlichting & Pipe Correlation Limitations

flat_plate · airfoil · curved_plate

Prandtl-Schlichting (1932/1950) — Critical Limitations

✗Re outside 5×10⁵–10⁹ — below the lower bound, flow may still be laminar or transitional; above 10⁹ the correlation is an extrapolation. Most engineering applications fall within this range; aircraft wings at cruise typically 10⁷–10⁹.

✗Zero pressure gradient (ZPG) only — the most critical limitation. Real airfoils have suction peaks near the leading edge (FPG) and adverse pressure gradients (APG) near the trailing edge. APG reduces Cf and can trigger separation. Use CFD with k-ω SST or SA model, or integral methods (Head's method, Drela's XFOIL) for pressure-gradient flows.

✗Fully turbulent from leading edge — the correlation assumes the BL is turbulent everywhere. In reality a laminar portion exists from x = 0 to the transition point x_tr. Using fully turbulent Cf overestimates friction drag by 20–50% for Re ≈ 10⁶–10⁷. Use the combined transition correlation: Cf = 0.074/Re^0.2 − 1742/Re (for transition at Re_x ≈ 5×10⁵).

✗Incompressible only (Ma < 0.3) — above Ma = 0.3, compressibility alters the density field, temperature near the wall (adiabatic wall temperature), and effective viscosity. Use the van Driest II compressible transformation or Spalding-Chi correlation for Ma > 0.3. See Compressibility card for full Mach thresholds.

⚠1/7 power-law for δ — the BL thickness formula δ = 0.37L/Re^0.2 comes from the 1/7 power-law profile. It overestimates δ at low Re and underestimates at high Re. For flows with strong pressure gradients, δ can differ from this estimate by a factor of 2–3.

⚠Smooth wall only — for rough surfaces (sand grains, machined surfaces, riblets), use the fully rough or transitionally rough forms of the Moody chart. The equivalent sand-grain roughness k_s must satisfy k_s⁺ = k_s·u_τ/ν > 70 for fully rough flow, where Cf is independent of Re and increases with k_s/L.

pipe · duct

Blasius (1913) & Petukhov-Filonenko Limitations

✗Blasius pipe: 4×10³ < Re < 10⁵ — outside this range the correlation is inaccurate. For Re < 4000 use Hagen-Poiseuille (laminar). For Re > 10⁵ use Petukhov-Filonenko. Both are smooth-pipe correlations.

✗Petukhov-Filonenko: smooth pipe only — for rough pipes, use Colebrook-White: 1/√f = −2·log(k_s/(3.7D) + 2.51/(Re·√f)), or the explicit Swamee-Jain approximation. The Moody diagram covers smooth and rough regimes.

✗Fully developed flow required — turbulent hydrodynamic entry length L_h ≈ 4.4·Re^(1/6)·D. For Re = 10⁶: L_h ≈ 40D. Within this entry region, local Cf is elevated above the fully-developed value. Entry effects are critical in short heat exchangers, orifice meters, and manifolds.

⚠D_h for non-circular ducts — hydraulic diameter D_h = 4A/P works well for aspect ratios up to ~4:1. For very flat rectangular channels (aspect > 8:1), the flow is dominated by the wide flat walls; use the dedicated correlation for parallel plates: f = 0.3164/Re^0.25 with D_h = 2H (channel half-height). For annuli: D_h = D_outer − D_inner.

⚠Straight duct only — curved ducts, bends, and S-bends generate secondary flow (Dean vortices) that increase pressure drop and modify the wall shear distribution. Local Cf peaks at the outer wall of bends.

COMPRESSIBILITY Mach Number Effects & Corrections

Mach Number Regimes

✓Ma < 0.3 — Incompressible: Density variation <5%. All correlations in this tool valid. For gases: standard Prandtl-Schlichting, Blasius, Hagen-Poiseuille apply directly.

⚠0.3 ≤ Ma < 0.8 — Subsonic Compressible: Density variation 5–50%. Use Eckert's reference temperature T* = T∞(1 + 0.032Ma² + 0.58(T_w/T∞−1)) to compute effective properties, then apply incompressible correlations at T*. Alternatively, use the van Driest II transformation for Cf.

✗0.8 ≤ Ma < 1.2 — Transonic: Shockwaves form on surfaces. Shock-induced boundary layer separation (SBLI) fundamentally alters the near-wall flow. Local Cf spikes at the shock foot. The incompressible flat-plate correlations are completely invalid. Require transonic CFD (density-based solver) with appropriate compressible turbulence models.

✗1.2 ≤ Ma < 5 — Supersonic: Oblique shocks, expansion fans, and strong aerodynamic heating dominate. Adiabatic wall temperature T_aw = T∞(1 + r·γ−1/2·Ma²) where r ≈ Pr^(1/3) ≈ 0.88 (turbulent recovery factor). Correlations are invalid; use CFD with real-gas or calorically perfect gas models.

✗Ma ≥ 5 — Hypersonic: Aerodynamic heating causes high-temperature gas effects (vibrational excitation, dissociation, ionisation). Sutherland viscosity law breaks down. Requires thermochemical non-equilibrium CFD. Completely outside the scope of this tool.

Compressible Cf Corrections (for reference)

→Van Driest II (1956): Cf_comp / Cf_incomp = f(Ma, T_w/T∞, Pr). Most accurate for adiabatic and non-adiabatic walls at Ma = 0.3–5. Ref: Van Driest, E.R. (1956), J. Aero. Sci. 23(11).

→Eckert Reference Temperature Method (1955): Evaluate all incompressible correlations at T* instead of T∞. Simple and reasonably accurate for Ma < 3. Ref: Eckert, E.R.G. (1955), J. Aero. Sci. 22.

→Spalding-Chi Transformation (1964): Transforms compressible problem to equivalent incompressible using velocity profile transformation. Used in many engineering textbooks.

Mach Number Estimation in This Tool

⚠Ma = U/a where a = √(γRT), T derived from ideal gas: T = p/(ρR) at standard pressure p = 101 325 Pa, R = 287 J/(kg·K), γ = 1.4. For non-standard altitudes or elevated-pressure applications (e.g. wind tunnels at 3 atm), input the actual measured density to obtain an accurate Mach number estimate.

SOLVER Wall Treatment & y⁺ Guidance (ANSYS Fluent & OpenFOAM)

ANSYS Fluent — Wall Treatment Options

→Standard Wall Functions (Launder & Spalding 1974): y⁺ = 30–300. Log-law applied at first cell. Activate in: Cell Zone Conditions → Fluid → Turbulence. Most robust for high-Re flows with attached BL.

→Scalable Wall Functions: Clips y⁺ to min 11.225 to prevent placement in buffer/sublayer. Safer for cases where y⁺ varies widely. Recommended over standard WF for general use.

→Enhanced Wall Treatment (EWT): Blends wall functions with low-Re formulation. y⁺ ≤ 1 gives best accuracy; also works with y⁺ up to ~300 (two-layer model). Required for heat transfer, natural convection, and APG flows. Activates with k-ε models.

→k-ω SST (Menter 1994): Automatically integrates to wall if y⁺ ≤ 1. No separate wall function activation needed. Default recommended model for most external aerodynamics and turbomachinery.

OpenFOAM — Wall Boundary Conditions

→Wall functions (y⁺ = 30–300): On wall patches set: nut → nutkWallFunction, k → kqRWallFunction, epsilon → epsilonWallFunction (for k-ε); or k → kLowReWallFunction, omega → omegaWallFunction (for k-ω SST). U → noSlip always.

→Low-Re (y⁺ ≤ 1): For k-ω SST set: nut → nutLowReWallFunction (or nutk), omega → omegaWallFunction, k → kLowReWallFunction. OpenFOAM's k-ω SST integrates to the wall automatically when y⁺ ≤ 1.

→nutk vs nutkWallFunction: nutk uses k to compute nut at the wall (more stable for k-ω models); nutkWallFunction uses log-law directly (standard for k-ε with wall functions).

Buffer Layer — Always Avoid (5 < y⁺ < 30)

✗Placing the first cell in the buffer layer (5 < y⁺ < 30) is the most common meshing mistake in CFD. In this region: (1) the log-law assumed by wall functions is not yet established, (2) direct resolution of the viscous sublayer is incomplete, (3) neither treatment is accurate. The resulting errors in Cf, heat flux, skin friction drag, separation point prediction, and turbulent kinetic energy can exceed 20–50%. Always check the y⁺ contour on all walls after solving and ensure y⁺ < 1 (resolved) or y⁺ > 30 (wall functions) everywhere.

ROUGH WALL Surface Roughness Effects

Roughness Regimes (Nikuradse 1933)

✓Hydraulically smooth: k_s⁺ = k_s·u_τ/ν < 5. Roughness elements are submerged in the viscous sublayer. Cf is independent of roughness; smooth-wall correlations apply.

⚠Transitionally rough: 5 ≤ k_s⁺ ≤ 70. Partial interaction of roughness elements with the turbulent buffer and log layers. Cf increases with both Re and k_s. Log-law shift: ΔB = (1/κ)·ln(k_s⁺) − 3.3 (Nikuradse).

✗Fully rough: k_s⁺ > 70. Cf becomes independent of Re (viscosity is irrelevant). Cf = [1.89 + 1.62·log(L/k_s)]^(−2.5) (Schlichting). All smooth-wall correlations in this tool are invalid. This applies to paint-roughened hulls, compressor blades with deposits, or machined surfaces at high Re.

Equivalent Sand-Grain Roughness k_s for Common Surfaces

•Polished metal / glass: k_s ≈ 0.0015–0.005 µm

•Hydraulically smooth drawn tubing: k_s ≈ 0.0015–0.007 mm

•Commercial steel pipe: k_s ≈ 0.046 mm

•Concrete / cast iron: k_s ≈ 0.3–3 mm

•Aircraft paint surface: k_s ≈ 0.01–0.1 mm

Check k_s⁺ before assuming smooth-wall correlation validity. Even “smooth” engineering surfaces can be hydraulically rough at high Re.

ACCURACY BL Integral & Inflation Stack Caveats

Integral Thickness Accuracy

⚠δ*, θ, H from power-law profiles — the 1/7 power-law (turbulent) and Blasius similarity (laminar) give engineering approximations. Under favourable pressure gradients, the profile is fuller (lower H); under adverse pressure gradients, it is fuller initially then becomes inflected near separation. Errors in δ* and θ of 10–30% are common in pressure-gradient flows.

✗Shape factor H > 2.5 (turbulent) or H > 3.5 (laminar) — indicates separation is imminent. Once H exceeds these values, integral methods and Cf correlations all break down. The BL has effectively separated and must be treated with full N-S CFD including turbulence model capable of capturing separation (k-ω SST recommended over k-ε).

⚠Re_θ and transition — the momentum-thickness Reynolds number Re_θ is used in Michel's transition criterion and the e^N method. The tool computes Re_θ from the correlation-derived θ, not from an actual velocity profile. Use as first estimate only.

Inflation Stack & Mesh Quality

⚠Add 10–20% margin to N — the computed N targets exactly δ from the correlation. The actual BL in the simulation may be thicker or thinner due to pressure gradients, 3D effects, inlet turbulence, or surface curvature. A margin ensures the inflation stack fully covers the BL.

⚠Growth ratio r = 1.1–1.2 recommended — r > 1.3 causes aspect ratio discontinuities at the inflation-to-bulk mesh interface, introducing numerical diffusion in the wall-normal direction. For high-accuracy LES or DNS, r ≤ 1.1. For engineering RANS, r = 1.15–1.2 is standard.

⚠Aspect ratio at wall — with very small Δy (y⁺ ≈ 1 at high Re), cells become highly stretched in the wall-normal direction (aspect ratio ~10³–10⁴). This is normal and expected for wall-resolved RANS. ANSYS Meshing reports this as high aspect ratio but it is acceptable for boundary layers — solvers are designed for it. Maximum skewness should remain < 0.85.

•Check y⁺ post-solve — always generate a y⁺ contour plot on all wall surfaces after the first simulation. If y⁺ > 1 anywhere for a resolved-wall model, Δy is too large; if y⁺ < 30 for a wall-function model, Δy is too small. Use post-processing results to guide mesh refinement.

CFL Timestep & Solver Limitations

Convective CFL Assumption

✗Explicit solvers: CFL > 1 causes instability — the CFL = 1 limit is a necessary but not sufficient stability condition. For multi-dimensional explicit schemes, the effective stability limit is CFL ≤ 1/√d where d is the number of spatial dimensions (CFL ≤ 0.577 for 3D). Violating this causes exponential divergence within a few timesteps.

⚠Implicit does not mean unconditionally stable — implicit pressure-based solvers (SIMPLE) are conditionally stable; very high CFL degrades the linear system conditioning and inner-iteration convergence. CFL > 100 in SIMPLE typically gives poor mass imbalance.

⚠Minimum Δx governs timestep — Δx is the smallest cell face-to-face distance in the mesh, not the average. Near walls and in regions of high resolution, Δx can be orders of magnitude smaller than the bulk mesh, forcing very small Δt globally in explicit solvers. This is why local time stepping (LTS) is used for steady-state explicit solutions.

•Acoustic CFL for compressible flows — in compressible flows at low Mach number (Ma < 0.3) the acoustic CFL dominates: CFL_a = (U + a)Δt/Δx where a is the speed of sound. This forces very small Δt even when convective CFL is small. Use pressure-based compressible solvers (Fluent: pressure-based + compressibility; OpenFOAM: sonicFoam, rhoPimpleFoam) to avoid acoustic restriction.

Flow-Through Time Assumptions

⚠5 T_ft is a minimum for turbulent statistics — rule of thumb only. Highly separated flows, recirculation zones, and bluff body wakes may require 20–50 T_ft for converged mean and RMS quantities. Always check running averages of key quantities (Cd, Cl, Nu).

⚠CPU estimate is order-of-magnitude only — actual runtime depends on hardware, parallel scaling efficiency, linear solver type (AMG vs ILU), and convergence of inner iterations. The tool assumes 100% parallel efficiency, which is unrealistic for > O(1000) cores.

DOMAIN Blockage & Domain Sizing Limitations

Blockage Ratio Assumptions

✗3% limit is for incompressible external aero only — the 3% blockage limit applies to low-speed, attached, incompressible external flows (airfoils, streamlined bodies). For bluff bodies (cylinders, buildings), < 1% is recommended because the separated wake is highly sensitive to confinement. For internal flows (ducts, passages), blockage has a different meaning and the limit does not apply.

✗Projected area approximation — the tool uses the frontal projected area W×t as the blocked area. For non-convex bodies, highly swept wings, or bodies with internal flow paths, the effective blockage differs. Use the maximum cross-sectional area perpendicular to the flow direction.

⚠Domain extent rules are for subsonic incompressible flows — for compressible flows (> Ma 0.3), characteristic waves reflect from far-field boundaries and can corrupt the solution near the body. Use non-reflective (far-field) boundary conditions (Riemann invariant-based) rather than fixed-pressure outlets. The required domain may be larger.

⚠Downstream extent depends on Strouhal number — for bluff bodies with periodic vortex shedding (St ≈ 0.2 for a cylinder), the wake extends 10–20D. For a body at high angle-of-attack with large separated wake, 30C downstream may be required. Always check that pressure has recovered to free-stream at the outlet boundary.

Boundary Condition Assumptions

•Pressure outlet assumes recovered pressure — placing a pressure outlet too close to the body means the boundary condition forces a pressure gradient that does not exist in free-air. This corrupts the wake and can cause backflow. Always confirm no backflow at outlet in post-processing.

•Symmetry BC assumes zero normal gradient — a symmetry plane imposes zero normal velocity and zero normal gradient of all scalars. This is valid only if the flow is genuinely symmetric. Using symmetry for asymmetric flows (post-stall airfoils, yawed bodies) introduces significant errors in Cd and Cl.

TRANSITION Transition Estimator Limitations

Michel Criterion Limitations

✗ZPG only — Michel (1952) is calibrated for zero pressure gradient flat plate flows. For airfoils with strong suction peaks (FPG near leading edge) or adverse pressure gradients (APG near trailing edge), the criterion can mispredict transition location by 20–50% of chord. Use the Abu-Ghannam-Shaw criterion or γ-Reθ SST for pressure-gradient flows.

✗Tu < 0.5% only — Michel is calibrated for natural transition (Tollmien-Schlichting pathway). For Tu > 0.5%, bypass transition mechanisms dominate and Michel overestimates x_tr (predicts transition too late).

✗3D effects not included — crossflow instability on swept wings, Görtler vortices on concave surfaces, and corner flow transitions are not captured by any of the 2D criteria in this tool. For swept wings (> 25° sweep), crossflow transition can dominate at cruise; use linear stability theory or PSE methods (HAL/NAL practice for LCA wing design).

Abu-Ghannam & Shaw Limitations

⚠Simplified Thwaites λ implementation — the full AGS criterion requires computing the Thwaites pressure gradient parameter λ = (θ²/ν)(dU/dx) from the actual velocity gradient along the body. This tool uses a user-specified λ as a constant over the surface, which is an approximation. For accurate results, extract d(U)/dx from a preliminary inviscid or RANS solution.

⚠Tu is assumed uniform — freestream Tu is assumed constant along the surface. In practice Tu decays with downstream distance from the turbulence-generating source. The AGS correlation is most accurate when Tu is measured at the location of interest.

General Transition Assumptions

•Smooth, clean wall assumed — surface roughness (paint, fasteners, steps) promotes early transition. A roughness Reynolds number k_s⁺ = k_s·u_τ/ν > 5 can trigger turbulence upstream of any predicted natural transition location. For rough surfaces, assume fully turbulent.

•Attached flow only — all criteria assume an attached laminar BL. In separated flow (laminar separation bubble), transition occurs inside the bubble via Kelvin-Helmholtz instability and the re-attached flow is turbulent. This is not captured by these criteria. Laminar separation bubbles are common on low-Re airfoils (ISRO SLV fins, NAL Saras wing).

•2D flat surface assumed — all Blasius-based computations of Re_θ(x) assume a flat plate. For curved surfaces, the actual Re_θ(x) differs due to surface curvature effects on the velocity gradient. This introduces errors of up to 15% for highly cambered airfoils.