Skip to main content

Velocity profile near smooth boundary

Short answer

Prandtl's mixing length theory gives τ_t = ρ(κy)²(du/dy)². With constant τ₀ near the wall, du/dy = u*/(κy). Integrating and matching to the laminar sublayer at y = 11.6 ν/u* (where u = 11.6 u*) yields u/u* = (1/κ) ln(u* y/ν) + 5.5. With κ=0.4, 1/κ=2.5, converting ln to log₁₀ gives 5.75 log₁₀(u* y/ν) + 5.5.

Explain Prandtl mixing length theory. Starting from the expression for turbulent shear stress derive the velocity distribution for turbulent flow near hydrodynamically smooth boundaries in the form uu∗=5.75log⁡10(u∗yν)+5.5\frac{u}{u_*} = 5.75\log_{10}\left(\frac{u_* y}{\nu}\right) + 5.5.

Prandtl's mixing length theory

In turbulent flow, lumps of fluid are moved across layers by eddies. Prandtl assumed a lump keeps its momentum over a distance ll (mixing length) before mixing. The velocity fluctuations are then u′≈v′≈l dudyu' \approx v' \approx l\,\dfrac{du}{dy}, and the turbulent shear stress is

τt=−ρu′v′‾=ρ l2(dudy)2\tau_t = -\rho\overline{u'v'} = \rho\,l^2\left(\frac{du}{dy}\right)^2

Near a wall the mixing length increases with distance from it: l=κyl = \kappa y with κ≈0.4\kappa \approx 0.4 (von Karman's constant).

Velocity distribution near a smooth boundary

Let u∗=τ0/ρu_* = \sqrt{\tau_0/\rho}. Close to the wall the shear stress is practically constant, τ=τ0\tau = \tau_0, and the viscous stress is negligible outside the laminar sub-layer:

τ0=ρκ2y2(dudy)2  ⇒  dudy=u∗κy\tau_0 = \rho\kappa^2y^2\left(\frac{du}{dy}\right)^2 \;\Rightarrow\; \frac{du}{dy} = \frac{u_*}{\kappa y}

Integrating:

u=u∗κln⁡y+C⋯(1)u = \frac{u_*}{\kappa}\ln y + C \quad\cdots(1)

Constant from the laminar sub-layer. In the laminar sub-layer of thickness δ′\delta', τ0=μu/y\tau_0 = \mu u/y, so u=τ0yμ=u∗2yνu = \dfrac{\tau_0y}{\mu} = \dfrac{u_*^2y}{\nu}. At its edge y=δ′y = \delta', experiments give δ′=11.6νu∗\delta' = \dfrac{11.6\nu}{u_*}, so

uδ′=u∗2ν⋅11.6νu∗=11.6 u∗u_{\delta'} = \frac{u_*^2}{\nu}\cdot\frac{11.6\nu}{u_*} = 11.6\,u_*

Applying (1) from y=δ′y = \delta' to yy:

u−11.6u∗=u∗κln⁡yδ′=u∗κln⁡u∗y11.6νu - 11.6u_* = \frac{u_*}{\kappa}\ln\frac{y}{\delta'} = \frac{u_*}{\kappa}\ln\frac{u_*y}{11.6\nu}

With κ=0.4\kappa = 0.4:

uu∗=11.6+2.5ln⁡u∗yν−2.5ln⁡11.6=2.5ln⁡u∗yν+11.6−6.13\frac{u}{u_*} = 11.6 + 2.5\ln\frac{u_*y}{\nu} - 2.5\ln 11.6 = 2.5\ln\frac{u_*y}{\nu} + 11.6 - 6.13

Convert ln⁡\ln to log⁡10\log_{10}: 2.5ln⁡x=2.5(2.303)log⁡10x=5.75log⁡10x2.5\ln x = 2.5(2.303)\log_{10}x = 5.75\log_{10}x:

uu∗=5.75log⁡10(u∗yν)+5.5\boxed{\frac{u}{u_*} = 5.75\log_{10}\left(\frac{u_*y}{\nu}\right) + 5.5}

(the exact constant 5.47 is rounded to 5.5, a value that agrees with Nikuradse's measurements).

✓ Answer checked by Claude Sonnet 5.5

All Hydraulics chapter 1 questions →