← Portfolio

Y⁺ Calculator

UNITS:
Flow Parameters
❯_Enter parameters and press Calculate
Flow Parameters
❯_Enter parameters and press Calculate
Compressible Flow Parameters
ℹ For compressible flows (Ma ≥ 0.3). Uses Eckert reference temperature method with Sutherland's viscosity law. Valid for calorically perfect air (Ma < 5).
ℹ Total (stagnation) pressure p₀ is the correct input for compressible boundary layer calculations. Static pressure is derived via the isentropic relation: p_static = p₀ / (1 + (γ−1)/2 · Ma²)γ/(γ−1).
Gas Properties (Air)
Sutherland: μ_ref = 1.716×10⁻&sup5; Pa·s
T_ref = 273.15 K, S = 110.4 K
❯_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.).
❯_Enter three mesh values 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.
❯_Enter parameters and press Calculate
Body Dimensions
ℹ Blockage ratio = A_body / A_domain × 100%. External aero: < 3%. Wind-tunnel simulation: < 5%. Bluff bodies: < 1%.
Domain Dimensions
❯_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.
❯_Enter parameters and press Calculate
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.
• Reynolds, O. (1883), Phil. Trans. R. Soc. Lond. 174, 935–982. DOI: 10.1098/rstl.1883.0029
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.
• Prandtl, L. (1925), Z. Angew. Math. Mech. 5, 136–139.
• CFD-Online: Skin Friction Coefficient
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.
• CFD-Online: Y+ Wall Distance Estimation
• CFD-Online: Dimensionless Wall 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.
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.
• Blasius, H. (1908), Grenzschichten in Flüssigkeiten mit kleiner Reibung, Z. Math. Phys. 56, 1–37.
• Schlichting & Gersten (2000) BL Theory, §7.3. Springer Link
• CFD-Online: Flat Plate Boundary Layer
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.
• Hagen, G.H.L. (1839), Ann. Phys. Chem. 46, 423–442.
• Poiseuille, J.L.M. (1840), C. R. Acad. Sci. Paris 11, 961–967, 1041–1048.
• Schlichting & Gersten (2000) BL Theory, §5.2. Springer Link
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°).
• White, F.M. (2016) Fluid Mechanics, 8th ed., McGraw-Hill, §7.3.
• Achenbach, E. (1968), J. Fluid Mech. 34(4), 625–639. DOI: 10.1017/S0022112068002120 (cylinder separation angles)
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).
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).
• Prandtl, L. (1927), Ergebnisse der Aerodynamischen Versuchsanstalt zu Göttingen, III, 1–5.
• Schlichting, H. & Gersten, K. (2000) Boundary Layer Theory, 8th ed., Springer, §21.2. Springer Link
• CFD-Online: 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.
• Blasius, H. (1913), Das Ähnlichkeitsgesetz bei Reibungsvorgängen in Flüssigkeiten, VDI Forschungsheft 131.
• CFD-Online: Fanning Friction Factor
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)).
• Filonenko, G.K. (1954), Teploenergetika 1(4), 40–44.
• Petukhov, B.S. (1970), Advances in Heat Transfer, Vol. 6, Academic Press, 503–564. DOI: 10.1016/S0065-2717(08)70153-9
• Colebrook, C.F. (1939), J. Inst. Civ. Eng. 11(4), 133–156. DOI: 10.1680/ijoti.1939.13150
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.
• White, F.M. (2016) Fluid Mechanics, 8th ed., McGraw-Hill, §7.3.
• Schlichting & Gersten (2000) BL Theory, §21.2. Springer Link
COMPRESSIBLE Compressible Boundary Layer & Cf Correlations
Sutherland's Viscosity Law — Sutherland (1893)
μ(T)= μ_ref · (T/T_ref)^(3/2) · (T_ref + S) / (T + S)
μ_ref= 1.716×10⁻&sup5; Pa·s    (air)
T_ref= 273.15 K
S= 110.4 K    (Sutherland constant for air)
Ref: Sutherland, W. (1893), Phil. Mag. Ser. 5, Vol. 36, pp. 507–531. Valid for T = 100–1900 K for air. Above 1900 K, dissociation occurs and Sutherland breaks down.
• CFD-Online: Sutherland's Law
Recovery Temperature & Adiabatic Wall Temperature
T_aw= T∞ · (1 + r · (γ−1)/2 · Ma²)
r (turb)= Pr^(1/3) ≈ 0.896    (turbulent recovery factor)
r (lam)= Pr^(1/2) ≈ 0.845    (laminar recovery factor)
T_0= T∞ · (1 + (γ−1)/2 · Ma²)    (total/stagnation temperature)
T_aw is the wall temperature when no heat transfer occurs (adiabatic). Recovery factor r accounts for viscous dissipation in the BL. For Pr = 0.72 (air): r_turb ≈ 0.896, r_lam ≈ 0.849. Ref: Schlichting & Gersten (2000) §10.4; Anderson (2017) §18.3.
• NASA: Boundary Layers
Eckert Reference Temperature Method — Eckert (1955)
T*= 0.5·(T∞ + T_w) + 0.22·(T_aw − T∞)
ρ*= p∞ / (R·T*)    (constant pressure across BL)
μ*= Sutherland(T*)
Re*= ρ*·U∞·L / μ*
Cf_comp= Cf_incomp(Re*) · (ρ*/ρ∞)    [based on ½ρ∞U²]
Evaluate all fluid properties at the reference temperature T*. The incompressible correlation with Re* gives the local skin friction at reference conditions; multiply by ρ*/ρ∞ to convert to the conventional Cf definition. Accurate to ±5% for Ma < 5, both laminar and turbulent. Ref: Eckert, E.R.G. (1955), J. Aero. Sci. 22, pp. 585–587.
• NASA: Viscosity & Temperature
Van Driest II Transformation — Van Driest (1956)
Cf_comp= Cf_incomp(Re·F_Re) / F_c
F_c= (T_w/T∞)^0.5 · [arcsin(A) + arcsin(B)]^2 / (A²+B²)
A, B= functions of Ma, T_w/T∞, γ, r
Most accurate compressible Cf transformation for adiabatic and non-adiabatic walls at Ma = 0.3–5. Transforms the compressible problem to an equivalent incompressible one using velocity profile transformation. Ref: Van Driest, E.R. (1956), J. Aero. Sci. 23(11), pp. 1007–1011.
• CFD-Online: Skin Friction Coefficient
Compressible y⁺ — Wall Properties
y⁺= ρ_w · u_τ · y / μ_w
ρ_w= p∞ / (R·T_w)    (thin BL: p_w ≈ p∞)
μ_w= Sutherland(T_w)
u_τ= √(τ_w / ρ_w)
Δy= y⁺_target · μ_w / (ρ_w · u_τ)
In compressible flows, y⁺ is defined using wall density and viscosity, not freestream. The wall temperature (adiabatic or specified) determines ρ_w and μ_w via ideal gas and Sutherland's law. For adiabatic walls at high Ma, T_w >> T∞, so ρ_w << ρ∞ and μ_w >> μ∞, leading to larger Δy than incompressible.
• CFD-Online: Y+ Wall Distance Estimation
Compressible BL Thickness
δ (lam)= 5.0·L / √Re* · √(T*/T∞)
δ (turb)= 0.37·L / Re*^0.2 · (T*/T∞)^0.5
The BL thickness increases in compressible flow due to reduced density near the hot wall. The (T*/T∞)^0.5 factor accounts for volumetric expansion. At Ma = 3 (adiabatic), δ is roughly 2–3× the incompressible value. Ref: Anderson (2017) Eq. 18.27; White (2006) §7.5.
Online Verification Links
•CFD-Online: Skin Friction Coefficient — incompressible and compressible Cf correlations
•CFD-Online: Sutherland's Law — viscosity-temperature relation for gases
•CFD-Online: Y+ Wall Distance Estimation — compressible y+ definition
•NASA GRC: Boundary Layers — compressible BL overview
•NASA GRC: Viscosity — temperature-dependent viscosity
•Anderson, J.D. (2017) Fundamentals of Aerodynamics, 6th ed., McGraw-Hill — Ch. 18 (Compressible BL)
•White, F.M. (2006) Viscous Fluid Flow, 3rd ed., McGraw-Hill — Ch. 7 (Compressible BL)
•Schlichting, H. & Gersten, K. (2000) Boundary Layer Theory, 8th ed., Springer — Ch. 10 (Compressible BL)
COMPRESSIBLE Total Pressure, Isentropic Relations & Correct Input for y+
Total (Stagnation) Pressure p₀ — Isentropic Relation
p₀= p∞ · (1 + (γ−1)/2 · Ma²)^(γ/(γ−1))
p∞= p₀ / (1 + (γ−1)/2 · Ma²)^(γ/(γ−1))    ← derived from p₀
Isentropic factor= (1 + (γ−1)/2 · Ma²)^(γ/(γ−1))
Why total pressure is the correct input: In external compressible flows, the measured or known quantity is the total (stagnation) pressure p₀ — this is what a Pitot probe reads in subsonic flow, and what is set by the upstream reservoir in supersonic nozzles. Static pressure p∞ varies with Ma and must be derived. For boundary layer sizing, p∞ (the static pressure across the BL, which is essentially constant normal to it) is needed to compute ρ∞, ρ*, and ρ_w via the ideal gas law. The derivation uses isentropic relations which are exact for calorically perfect gas, adiabatic, inviscid flow.
• Anderson, J.D. (2017) Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, §8.4 — Isentropic flow relations.
• Anderson, J.D. (2003) Modern Compressible Flow, 3rd ed., McGraw-Hill, §3.4.
7-Step Temperature Chain for Compressible y+ Sizing
STEP 1: T∞Freestream static temperature (known/measured)
STEP 2: T₀= T∞ · (1 + (γ−1)/2 · Ma²)    — total/stagnation temperature
STEP 3: r= Pr^0.5 (laminar)   or   Pr^(1/3) (turbulent)    — recovery factor
STEP 4: T_aw= T∞ · (1 + r · (γ−1)/2 · Ma²)    — adiabatic wall temperature
STEP 5: T_w= T_aw (adiabatic)   or   user-specified T_w (isothermal)
STEP 6: T*= 0.5(T∞+T_w) + 0.22(T_aw−T∞)    — Eckert reference temperature
STEP 7: ρ*, μ*ρ* = p_static/(R·T*), μ* = Sutherland(T*)    — reference-state properties
Why T_aw is correct (not T behind shock): For an inviscid outer flow, the wall temperature in an adiabatic case is the adiabatic wall temperature T_aw, which accounts for viscous dissipation in the boundary layer via the recovery factor r. This is NOT the post-shock static temperature. Using the post-shock temperature equalized to T_w is physically incorrect because: (1) shocks occur at the nose, not everywhere along the body; (2) after the shock, the flow re-accelerates; (3) the wall itself does not "know" about the shock — it only sees the local BL physics. The correct T_w for y+ sizing is T_aw (if adiabatic) or the actual wall temperature (if radiatively or conductively cooled/heated).
• Eckert, E.R.G. (1955), J. Aero. Sci. 22(8), 585-587.
• White, F.M. (2006) Viscous Fluid Flow, 3rd ed., §7.6.
• Schlichting & Gersten (2000) BL Theory, 8th ed., §10.4.
ROCKET BODIES Mangler Transformation, Blunt Bodies & Rocket CD Physics
cone · ogive nose · rocket nose · axisymmetric body
Mangler Transformation — W. Mangler (1948)
Re_eff= Re_L · sin(α)^(2/3)    — equivalent flat-plate Re
α= cone half-angle (default 10° in this calculator)
Mangler factor= sin(10°)^(2/3) ≈ 0.353    — BL thinner than flat plate at same Re
Cf_cone (lam)= 1.328 / √Re_eff
Cf_cone (turb)= 0.370 / (log₁₀ Re_eff)^2.584
The Mangler transformation maps the axisymmetric BL on a cone (or ogive nose) to an equivalent 2D flat-plate problem by using the effective Reynolds number Re_eff = Re_L · sin(α)^(2/3). This accounts for the fact that streamlines are diverging on a cone — the BL is thinner than on a flat plate of the same length. L = axial length of the cone nose.
• Mangler, W. (1948), Z. Angew. Math. Mech. 28(4), 97-103.
• White, F.M. (2006) Viscous Fluid Flow, 3rd ed., §7.4 (Axisymmetric BL).
• Schlichting & Gersten (2000) BL Theory, 8th ed., §11 (Axisymmetric BL).
blunt body · re-entry capsule · hemisphere nose
Hiemenz Stagnation-Point BL — Hiemenz (1911)
K= U∞ / R_n    — velocity gradient at stagnation, R_n = nose radius [m⁻¹·s⁻¹]
δ_stag= 2.4 · √(ν/K)    — stagnation BL thickness [m]
τ_stag= μ · U∞ / δ_stag    — wall shear at stagnation [Pa]
At the stagnation point of a blunt body (capsule, hemisphere, ICBM nose), the velocity field is u(x) = K·x near the stagnation point. The boundary layer is characterized by the Hiemenz flow solution. The BL thickness is set by the strain rate K = U∞/R_n, NOT by the plate length. For re-entry capsules, drag is dominated by pressure drag (~70-80% of Cd), so y+ sizing matters primarily for accurate heat-flux prediction (aerothermal engineering).
• Hiemenz, K. (1911), Dingl. Polytechn. J. 326, 321-326.
• White, F.M. (2006) Viscous Fluid Flow, 3rd ed., §3.3 (Stagnation-point flow).
• Fay, J.A. & Riddell, F.R. (1958), J. Aero. Sci. 25(2), 73-85 (aerothermal heating).
sounding rocket · ballistic missile · orbital LV · SRB
Rocket CD Physics — Why y+ Matters
CD_total= CD_wave + CD_pressure + CD_friction + CD_base
CD_friction= 2 · Cf · A_wetted / A_ref    — skin friction drag component
Subsonic rocketCD_friction ≈ 40-60% of CD_total    — dominant at low Ma
Supersonic rocketCD_wave ≈ 40-70% of CD_total    — wave drag dominates
Composite mesh: nosey+ ≤ 1  (Mangler BL, heat flux critical)
Composite mesh: bodyy+ ≤ 1  (skin friction drag)
Composite mesh: boat-taily+ ≤ 1 + 50% extra layers  (APG, separation)
Composite mesh: finsy+ ≤ 1 using root chord, turbulent from LE
Why correct y+ sizing is critical for rocket CD: At subsonic speeds, skin friction drag contributes 40-60% of total drag on a well-streamlined rocket. Incorrect y+ (too large) causes the solver to under-resolve the viscous sublayer, leading to wall function error in Cf and hence error in CD_friction. At supersonic speeds, wave drag is larger, but the boat-tail and base region still have significant pressure-drag sensitivity to boundary layer separation — which is controlled by the near-wall mesh quality. For accurate CD prediction, y+ ≤ 1 with k-ω SST is the industry standard.
• Barrowman, J.S. (1967), NACA-TM-X-57853 (sounding rocket aerodynamics).
• ESDU 87034 (1987): Drag of bodies of revolution in compressible flow.
• Anderson, J.D. (2017) Fundamentals of Aerodynamics, 6th ed., Ch. 17 (Compressible drag).
• Menter, F.R. (1994), AIAA J. 32(8), 1598-1605 (k-ω SST — recommended for rocket CFD).
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.
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).
• Prandtl, L. (1925), Z. Angew. Math. Mech. 5(2), 136–139.
• Schlichting, H. & Gersten, K. (2000) Boundary Layer Theory, 8th ed., Springer, §17.1. Springer Link
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).
• von Kármán, T. (1930), Nachr. Ges. Wiss. Göttingen, Math.-Phys. Kl., 58–76.
• Schlichting & Gersten (2000) BL Theory, §17.2. Springer Link
• CFD-Online: Law of the Wall
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.
• Spalding, D.B. (1961), A single formula for the 'law of the wall', J. Appl. Mech. 28(3), 455–458. DOI: 10.1115/1.3641728
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.
• Coles, D. (1956), J. Fluid Mech. 1(2), 191–226. DOI: 10.1017/S0022112056000135
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)
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.
• Courant, R., Friedrichs, K. & Lewy, H. (1928), Math. Ann. 100(1), 32–74. DOI: 10.1007/BF01448839
• Blazek, J. (2015) Computational Fluid Dynamics: Principles and Applications, 3rd ed., Elsevier, §6.3. Elsevier
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.
• Ferziger, J.H., Perić, M. & Street, R.L. (2020) Computational Methods for Fluid Dynamics, 4th ed., Springer, §7.3. Springer Link
• ANSYS Fluent Theory Guide 2023 R2, §26.3–26.4 (time advancement).
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).
• Bearman, P.W. (1984), Annu. Rev. Fluid Mech. 16, 195–222. DOI: 10.1146/annurev.fl.16.010184.001211
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.
• Barlow, J.B., Rae, W.H. & Pope, A. (1999) Low-Speed Wind Tunnel Testing, 3rd ed., Wiley, §9.4.
• Maskell, E.C. (1965), A theory of the blockage effects on bluff bodies in a closed wind tunnel, RAE Report & Memoranda No. 3400.
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.
• Vassberg, J.C. & Jameson, A. (2010), In pursuit of grid convergence for two-dimensional Euler solutions, AIAA-2010-4438. DOI: 10.2514/6.2010-4438
• NASA/TM-2006-214293 (external aero best practice guidelines).
• COST Action 732 (2007), Best Practice Guideline for the CFD simulation of flows in the urban environment.
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.
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.
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.
ROCKET / COMPRESSIBLE Rocket Geometries & Total Pressure Assumptions
Total Pressure Input — Why It Is Correct
✓Total pressure p₀ is the correct input for compressible BL calculations. Static pressure p∞ is derived via the isentropic relation p∞ = p₀ / (1 + (γ−1)/2 · Ma²)^(γ/(γ−1)). Pitot probes read stagnation pressure directly; wind tunnel reservoir pressure is the total pressure. These are the physically measured quantities. Ref: Anderson (2017) §8.4.
✗Do not use static pressure as input when the true measurement is total pressure. If you input static pressure directly, you are assuming the flow is at rest everywhere, which is wrong for any compressible flow. The correct procedure is to input p₀ and let the calculator derive p_static for density computation.
⚠In subsonic incompressible flow (Ma < 0.3), the distinction between static and total pressure is small (isentropic factor ≈ 1), so both are approximately equal. The correction becomes important above Ma ≥ 0.3 — the Compressible tab is designed for this regime and correctly handles the distinction.
Wall Temperature for y+ Sizing — Adiabatic vs Post-Shock
✓Adiabatic wall: use T_aw (recovery temperature), not the temperature behind the shock. T_aw = T∞ · (1 + r · (γ−1)/2 · Ma²). The recovery factor r accounts for viscous dissipation in the BL. This is what the wall "sees" when there is no active cooling or heating. Ref: Eckert (1955), White (2006) §7.6.
✗Do not equalize the wall to the post-shock temperature. The post-shock static temperature is a property of the flow at the shock location (e.g., at the nose), not the thermal condition of the wall. After the shock, the flow re-accelerates and cools. The wall at a distance from the nose does not experience the post-shock temperature — it experiences T_aw, which is determined by the local Mach number, recovery factor, and freestream temperature, not by the shock. Equalizing T_w to T_post-shock overcools the wall in supersonic calculations and produces wrong ρ_w and μ_w.
⚠Isothermal wall: If the wall is actively cooled (ablative cooling, film cooling, radiation equilibrium), specify T_w directly. For radiation equilibrium: T_w = (q_rad/(sigma*epsilon))^0.25. For ablative: T_w is the ablation temperature of the material (e.g. ~3600 K for carbon-carbon). These must be user-specified; the tool does not compute them.
Composite Rocket Body — Mesh Sizing Strategy
✓Nose cone (sharp cone / ogive): Use the Mangler transformation. BL is thinner by sin(α)^(2/3) vs flat plate. Run the calculator in Compressible tab with geometry = "Cone" or relevant rocket nose type. L = axial nose length.
✓Cylindrical body: Use the Compressible tab with flat plate baseline. BL thickens along the body. For y+ sizing, use the design condition (maximum dynamic pressure, typically Ma ≈ 0.8–1.5 at max-Q for sounding rockets). L = total body length.
✗Boat-tail / afterbody: Adverse pressure gradient causes BL deceleration and potential separation. The flat-plate Cf OVERESTIMATES the actual skin friction. Use 50% more prism layers than the calculated N for safety margin. Separation is likely near the boat-tail shoulder.
✓Fins: BL is TURBULENT from the leading edge of the fin (body wake is already turbulent). Use root chord as L, Turbulent tab (or Compressible for Ma > 0.3). The Laminar tab underestimates Δy for fins in flight by 2–5×.
⚠Base region: This tool does not compute base drag or the base BL. The base (aft face) has a separated recirculation zone. Near-base y+ should still satisfy y+ ≤ 1 for the base wall to capture reattachment shock if present. Use the calculated Δy from the body as a conservative estimate for the base.
Mangler Transformation Assumptions
✓Valid for: Slender cones with half-angle ≤ 15°, attached laminar or turbulent BL, no crossflow or swirl. Extends the 2D flat-plate Blasius / Prandtl-Schlichting correlations to axisymmetric bodies by computing an effective flat-plate Re. Ref: Mangler (1948).
✗Not valid for: Half-angles > 20° (flow may separate at the tip), non-conical ogive shapes (use numerical BL integration), supersonic shock-detachment conditions, crossflow components (spinning rockets, angle-of-attack > 5°). At angle of attack, crossflow-induced BL transition occurs much earlier.
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.
COMPRESSIBLE Compressible Flow Assumptions & Limitations
Perfect Gas & Sutherland Validity
•Calorically perfect gas assumed — γ is constant (1.4 for air). Above ~800 K, vibrational modes activate and γ decreases. Above ~2500 K, dissociation begins (O&sub2; → 2O). For temperatures above 800 K, use thermally perfect gas tables or CEA.
•Sutherland's law valid for 100–1900 K — below 100 K, the power-law model is more appropriate. Above 1900 K, chemical reactions invalidate the single-species viscosity model. For hypersonic flows (Ma > 5), use Blottner or Gupta-Yos curve fits.
•Air only — Sutherland constants (μ_ref, T_ref, S) are for dry air. For other gases (N&sub2;, CO&sub2;, He), different constants are required.
Reference Temperature Method Accuracy
•Eckert reference temperature: ±5% for Ma < 5 — validated against direct numerical simulations and experimental data for flat plates. For Ma > 5, real-gas effects dominate and T* loses accuracy.
•Constant pressure across BL assumed — ρ* = p∞/(RT*). This is valid for thin BLs where dp/dy ≈ 0. For strong shock-BL interaction (e.g. reflected shocks), the pressure varies significantly across the BL.
•Flat plate or slender body geometry — no pressure gradient effects. For bodies with strong curvature, shock impingement, or adverse pressure gradients, the reference temperature method should be supplemented with CFD pre-analysis.
Compressible y⁺ Definition
•Wall properties used for y⁺ — y⁺ = ρ_w · u_τ · y / μ_w. This is the correct compressible definition. Using freestream ρ∞ and μ∞ (as in the incompressible tabs) would produce incorrect first-layer heights.
•Adiabatic wall default — if no wall temperature is specified, T_w = T_aw (recovery temperature). For cooled walls (e.g. turbine blades), specify T_w explicitly. For actively heated walls, T_w > T_aw is possible but unusual.
•Δy is larger than incompressible — at Ma = 2 with adiabatic wall, the compressible Δy is typically 2–4× the incompressible value because μ_w >> μ∞ and ρ_w << ρ∞. This means compressible meshes can use coarser first-layer cells while maintaining the same y⁺.
SOURCES References & Literature
Boundary Layer & Skin Friction
•Blasius, H. (1908), Grenzschichten in Flüssigkeiten mit kleiner Reibung, Zeitschrift für Mathematik und Physik 56, 1–37.
•Schlichting, H. & Gersten, K. (2000), Boundary-Layer Theory, 8th ed., Springer.
•Prandtl, L. & Schlichting, H. (1932), Das Widerstandsgesetz rauher Platten, Werft, Reederei, Hafen 13.
•White, F.M. (2006), Viscous Fluid Flow, 3rd ed., McGraw-Hill.
•Hagen, G. (1839) / Poiseuille, J.L.M. (1840) — laminar pipe flow, Hagen–Poiseuille law.
•Petukhov, B.S. (1970), Heat transfer and friction in turbulent pipe flow, Adv. Heat Transfer 6, 503–564.
•Mangler, W. (1948), Zusammenhang zwischen ebenen und rotationssymmetrischen Grenzschichten, Z. Angew. Math. Mech. 28(4), 97–103.
•Hiemenz, K. (1911), Dinglers Polytechnisches Journal 326, 321–326 — stagnation-point boundary layer.
Transition & Turbulence Modelling
•Michel, R. (1951), Étude de la transition sur les profils d'aile, ONERA Report 1/1578A.
•Abu-Ghannam, B.J. & Shaw, R. (1980), Natural transition of boundary layers, J. Mech. Eng. Sci. 22(5), 213–228.
•Menter, F.R., Langtry, R.B. et al. (2006), A Correlation-Based Transition Model Using Local Variables — Part I, J. Turbomach. 128(3), 413–422. doi:10.1115/1.2184352
•Langtry, R.B. & Menter, F.R. (2009), Correlation-Based Transition Modeling for Unstructured Parallelized CFD Codes, AIAA J. 47(12), 2894–2906.
•Pope, S.B. (2000), Turbulent Flows, Cambridge University Press.
•ANSYS Inc., ANSYS Fluent Theory Guide, Release 23 R2.
Compressibility & Gas Properties
•Eckert, E.R.G. (1955), Engineering Relations for Friction and Heat Transfer to Surfaces in High Velocity Flow, J. Aeronautical Sciences 22(8), 585–587.
•Sutherland, W. (1893), The viscosity of gases and molecular force, Phil. Mag. Series 5, 36(223), 507–531.
•Van Driest, E.R. (1956), The Problem of Aerodynamic Heating, Aeronautical Engineering Review 15(10), 26–41.
•Coleman, G.N. & Mansour, N.N. (1991), CTR Annual Research Briefs, NASA Ames/Stanford.
Grid Convergence (GCI)
•Celik, I.B., Ghia, U., Roache, P.J., Freitas, C.J., Coleman, H. & Raad, P.E. (2008), Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications, ASME J. Fluids Eng. 130(7), 078001. doi:10.1115/1.2960953