Skip to main content

First layer height from y+ estimate

Short answer

Re_L = 3×10^6; Cf = 0.00309; wall shear stress = 1.67 Pa; friction velocity = 1.18 m/s. First cell height = 2.55×10^-5 m (0.025 mm). With growth ratio 1.2, about 30 prism layers are needed to cover the boundary layer thickness of 0.028 m.

Air (ρ=1.2\rho = 1.2 kg/m³, μ=1.8×10−5\mu = 1.8\times10^{-5} Pa·s) flows at 30 m/s over a flat plate 1.5 m long. A wall-resolved turbulence model requires y+=1y^+ = 1 at the first cell centre. Estimate (a) the plate Reynolds number, (b) the skin friction coefficient using Cf=0.026 ReL−1/7C_f=0.026\,Re_L^{-1/7}, (c) the wall shear stress and friction velocity, (d) the first cell height, and (e) the number of prism layers with growth ratio 1.2 needed to cover the boundary layer thickness estimated from δ=0.37 L ReL−1/5\delta = 0.37\,L\,Re_L^{-1/5}.

Definition: y+=ρuτyμy^+=\dfrac{\rho u_\tau y}{\mu} with uτ=τw/ρu_\tau=\sqrt{\tau_w/\rho}. The estimate is done before meshing using flat-plate correlations.

(a) Reynolds number

ReL=ρULμ=1.2×30×1.51.8×10−5=3.0×106Re_L=\frac{\rho UL}{\mu}=\frac{1.2\times30\times1.5}{1.8\times10^{-5}}=3.0\times10^{6}

(turbulent over most of the plate).

(b) Skin friction coefficient

Cf=0.026 (3×106)−1/7=0.026×0.1188=3.09×10−3C_f=0.026\,(3\times10^6)^{-1/7}=0.026\times0.1188=3.09\times10^{-3}

(c) Wall shear stress and friction velocity

τw=12CfρU2=0.5×3.09×10−3×1.2×302=1.67 Pa\tau_w=\tfrac12C_f\rho U^2=0.5\times3.09\times10^{-3}\times1.2\times30^2=1.67\ \text{Pa} uτ=τwρ=1.6671.2=1.179 m/su_\tau=\sqrt{\frac{\tau_w}{\rho}}=\sqrt{\frac{1.667}{1.2}}=1.179\ \text{m/s}

(d) First cell height

Kinematic viscosity ν=μ/ρ=1.5×10−5\nu=\mu/\rho=1.5\times10^{-5} m²/s. For y+=1y^+=1 the cell centre is at

y=y+νuτ=1×1.5×10−51.179=1.27×10−5 my=\frac{y^+\nu}{u_\tau}=\frac{1\times1.5\times10^{-5}}{1.179}=1.27\times10^{-5}\ \text{m}

The first cell has its centre at half its height, so the first cell height is

Δy1=2y=2.55×10−5 m≈0.025 mm\Delta y_1=2y=2.55\times10^{-5}\ \text{m}\approx 0.025\ \text{mm}

(e) Number of layers

Boundary layer thickness:

δ=0.37 L ReL−1/5=0.37×1.5×(3×106)−0.2=0.0281 m\delta=0.37\,L\,Re_L^{-1/5}=0.37\times1.5\times(3\times10^6)^{-0.2}=0.0281\ \text{m}

For nn layers with growth ratio g=1.2g=1.2, total thickness is Δy1gn−1g−1≥δ\Delta y_1\dfrac{g^n-1}{g-1}\ge\delta:

gn≥1+δ(g−1)Δy1=1+0.0281×0.22.55×10−5=221.4g^n\ge1+\frac{\delta(g-1)}{\Delta y_1}=1+\frac{0.0281\times0.2}{2.55\times10^{-5}}=221.4 n≥ln⁡221.4ln⁡1.2=29.6 ⇒ n=30 layersn\ge\frac{\ln221.4}{\ln1.2}=29.6\ \Rightarrow\ n=30\ \text{layers}

Check: with 30 layers the total is 2.55×10−5×1.230−10.2=0.03012.55\times10^{-5}\times\dfrac{1.2^{30}-1}{0.2}=0.0301 m, which is above δ\delta.

Answer: ReL=3×106Re_L=3\times10^6; Cf=3.09×10−3C_f=3.09\times10^{-3}; τw=1.67\tau_w=1.67 Pa; uτ=1.18u_\tau=1.18 m/s; first cell height ≈2.5×10−5\approx2.5\times10^{-5} m; about 30 prism layers.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Applied Computational Fluid Dynamics chapter 4 questions →