Boundary layer

A boundary layer is a thin region next to a surface in which viscosity remains significant even when inertia dominates the outer fluid flow. Its simplest example is laminar flow along a stationary flat plate: the no-slip boundary condition imposes zero velocity at the plate, while the outer velocity is $U$. The resulting thin velocity profile contains a large normal velocity gradient, so small dynamic viscosity can still produce appreciable wall shear stress. The Blasius boundary layer makes this competition quantitative and explains how a plate aligned with the incoming fluid flow acquires drag.

The same balance appears in a singular perturbation of a differential equation: large derivatives in a thin region can keep a term with a small coefficient significant. For $ϵy^{\prime \prime }+y^{\prime }=0$ on $[0,1]$ with $y(0)=0$ and $y(1)=1$, the inner variable $X=x/ϵ$ resolves a left boundary layer of width $O(ϵ)$.

In an attached two-dimensional boundary layer with zero pressure gradient, use the stream function $\Psi$ as transverse coordinate, assuming ${\Psi }_Y=u>0$ away from the wall. If $u=U(x,\Psi )$ and ${\Psi }_x=-V$, the chain rule gives $u_x=U_x-V{U}_{\Psi }$, $u_Y=U{U}_{\Psi }$ and $u_{YY}=U(U{U}_{\Psi }{)}_{\Psi }$. Substitution into $u{u}_x+V{u}_Y=u_{YY}$ proves the displayed degenerate nonlinear diffusion equation. Its flat-plate similarity form $U=g(\Psi /\sqrt{x})$ obeys $(g{g}^{\prime }{)}^{\prime }+\eta g^{\prime }/2=0$, $g(0)=0$, $g(\infty )=1$.

The wall coordinate is singular: finite physical shear corresponds to $g\sim C\sqrt{\eta }$, rather than a finite transformed derivative. A monotone profile exists by putting $q=g{g}^{\prime }$, $q(0)=a>0$. Then $q^{\prime }/q=-\eta /(2g)$ and $(g^2{)}^{\prime }=2q$. Since $q\le a$, $g\le \sqrt{2a\eta }$; this forces $q$ to decay at least as $\exp [-c{\eta }^{3/2}]$. Thus $g$ has a finite positive limit, which can be normalized to one by $g(\eta )\mapsto Ag(\eta /\sqrt{A})$.

A correction of amplitude $O(\delta )$ varying on width $O(\delta )$ can be negligible in a function's leading outer approximation while producing an order-one change in its derivative. For example $-\delta e^{-x/\delta }$ has derivative $e^{-x/\delta }$, allowing an endpoint derivative condition to be met without changing the leading outer boundary value. Such weak layers matter in matched asymptotic expansions involving derivative boundary conditions, even when only another endpoint has an order-one layer in the function itself.

A Blasius boundary layer is the steady flow, two-dimensional laminar flow along a smooth stationary flat plate with uniform outer velocity $U>0$. Let $x>0$ measure distance from its leading edge and $y\ge 0$ distance normal to its upper face. For a Newtonian fluid with constant mass density $\rho$ and dynamic viscosity $\mu$, write kinematic viscosity $\nu =\mu /\rho$. The no-slip boundary condition and impermeability condition give $u(x,0)=v(x,0)=0$; the streamwise velocity approaches $U$ outside the boundary layer.

For example, choose $U=10\mathrm{m}{\mathrm{s}}^{-1}$ and $\nu =1.5\times 10^{-5}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. At $x=0.30\mathrm{m}$, the local Reynolds number is ${\mathrm{Re}}_x=Ux/\nu =2.0\times 10^5$ and the boundary layer thickness is about $3.29\mathrm{mm}$. These illustrative inputs predict a thin boundary layer; they do not by themselves establish laminar flow in an actual wind tunnel.

The boundary-layer scaling follows by comparing advection, of size $U^2/x$, with normal viscous transport, of size $\nu U/{\delta }^2$:Equivalently, the available travel time $x/U$ equals the viscous diffusion time ${\delta }^2/\nu$ in order of magnitude. Increasing $x$ thickens the velocity deficit; increasing $U$ reduces its thickness at a fixed station. This boundary-layer scaling requires ${\mathrm{Re}}_x\gg 1$. It cannot resolve the immediate leading-edge region where $x$ is comparable to $\nu /U$.^[ref]

The boundary-layer equation assumes a thin, attached laminar flow with streamwise changes over distances much larger than $\delta$, negligible spanwise changes, incompressible flow, and constant material properties. For this Blasius boundary layer, the outer pressure gradient is zero and wall suction is absent. Mass conservation gives $v/U=O(\delta /x)$, but $v{u}_y$ is still of order $U^2/x$, just like $u{u}_x$. Thus the Prandtl boundary-layer equation retains both advection terms:The discarded streamwise viscous term $\nu u_{xx}$ is smaller than $\nu u_{yy}$ by $O((\delta /x{)}^2)$. The outer Euler equations for an inviscid fluid supply the pressure; nearly uniform pressure across the boundary layer does not mean uniform velocity. Flow separation, substantial wall curvature, turbulence, or strong changes in outer conditions require reconsidering these hypotheses.

Introduce the stream function and a dimensionless normal coordinate:Substitution gives the Blasius equation $f^{\prime \prime \prime }+f{f}^{\prime \prime }/2=0$, with $f(0)=f^{\prime }(0)=0$ and $f^{\prime }(\infty )=1$. Its numerical boundary value problem gives $f^{\prime \prime }(0)\simeq 0.332057$. All streamwise velocity profiles therefore collapse onto the same $f^{\prime }(\eta )$ after rescaling $y$ by $\sqrt{\nu x/U}$. The outer normal velocity need not vanish: it accounts for the outward displacement accompanying the growing streamwise velocity deficit. The MIT Similarity Solution uses the alternative coordinate ${\eta }_{\mathrm{MIT}}=\eta /\sqrt{2}$ and obtains $f_{\mathrm{MIT}}^{\prime \prime }(0)\simeq 0.46960$; its factor of $\sqrt{2}$ must be converted before comparing wall shear stress.

The historical Motion of fluids with very little viscosity by L. Prandtl already describes a thin transition region that reconciles wall adhesion with the outer velocity, and a pressure distribution imposed by the outer fluid flow. Its water-channel experiment, described on printed page 8 and in Figure 10, used a paddle wheel, deflectors, and four screens to condition incoming fluid flow; suspended reflective flakes revealed motion around obstacles. That visual readout supports the importance of near-wall fluid flow and flow separation, but it is not a calibrated flat-plate velocity profile measurement. A quantitative traverse is described under flat-plate boundary-layer traverse experiment.

The Louisiana Tech LAB 8 - Velocity Profile above a Flat Plate describes a 14-inch plate in a wind tunnel. A Pitot-static tube of 0.10-inch outside diameter traverses normally about 12 inches downstream; its screw scale has 0.001-inch increments. A water-column manometer reads the difference between stagnation pressure and static pressure. The procedure establishes a wall position, zeros the tunnel indicator, sets fan speed, and records at least ten positions inside the boundary layer, together with room temperature and barometric pressure. The instrument list includes a flow-stream thermometer implemented with a temperature-sensitive resistor. The manual requests comparison of outer velocity with the tunnel indicator and of the recovered velocity profile with laminar or turbulent predictions; it provides no observed dataset. These are apparatus specifications and instructions, not a reported confirmation of laminar flow.

Under incompressible flow conditions, an ideal aligned Pitot tube gives $q(y)=p_0(y)-p(y)=\rho u(y{)}^2/2$. Convert the manometer head $h$ to $q=({\rho }_{\mathrm{water}}-\rho )gh$, then infer $u=\sqrt{2q/\rho }$. Here $g$ is gravitational acceleration; mass density $\rho$ can be estimated from temperature and absolute pressure using the ideal gas law. A normalized velocity profile rises from the unresolved near-wall region to a plateau, whose value defines $U$; downstream profiles should broaden and collapse when plotted against $\eta$ if the Blasius boundary layer assumptions hold. Integrating a recovered velocity profile gives displacement thickness and momentum thickness.

Additional controls for interpreting this experiment follow from the measurement model: check differential-pressure calibration with equal port pressures and known liquid heads; leak-check tubing; align the Pitot tube with the incoming velocity; convert wall contact to the inlet-centre coordinate rather than treating outside contact as $y=0$; repeat traverses and outer velocity readings to detect drift. Keep the plate smooth and aligned, and check outer static pressure at multiple downstream stations before assuming zero pressure gradient. Repeating at another station or outer velocity tests the predicted boundary-layer scaling more directly than fitting one curve. These controls are proposed checks, not actions documented as completed in the manual.

Finite probe size, wall interference, flow fluctuations, and weak near-wall dynamic pressure limit inference. Even the manual's probe-based Reynolds number criterion, ${\mathrm{Re}}_d=ud/\nu >1000$ for its ideal Bernoulli equation conversion, must be checked with local $u$, not just $U$. In the illustrative model at $x=0.3048\mathrm{m}$, ${\delta }_{99}\simeq 3.32\mathrm{mm}$; the quoted outside diameter is $2.54\mathrm{mm}$, a substantial fraction of that thickness. A real inlet's sensing footprint and correction cannot be deduced from outside diameter alone. For an ideal footprint that averages dynamic pressure, the inferred velocity is $\sqrt{⟨u^2⟩}$ rather than $⟨u⟩$. Differentiating such sparse readings to obtain wall shear stress magnifies measurement error; the momentum thickness route is less sensitive to this differentiation but still needs the unresolved near-wall contribution. The no-slip boundary condition is a model constraint here, not a direct reading from a probe centred on the wall.

The conventional boundary layer thickness ${\delta }_{99}$ is the normal distance at which the streamwise velocity first reaches $0.99U$ in a monotone velocity profile. A Blasius boundary layer approaches $U$ smoothly, so its apparent edge depends on the chosen percentage; there is no discontinuity. Solving $f^{\prime }({\eta }_{99})=0.99$ in the Blasius equation gives ${\eta }_{99}\simeq 4.910$, henceThe factor $4.910$ belongs to this definition, whereas the boundary-layer scaling specifies only order of magnitude. Near this edge the velocity changes slowly, making a threshold position sensitive to measurement error; displacement thickness and momentum thickness summarize the whole velocity profile.^[ref]

For constant-mass density incompressible flow, the displacement thickness measures the lost volume flux per unit span divided by the outer velocity:An outer inviscid flow encountering the same volume flux deficit behaves, at leading order, as though the wall were displaced outward by ${\delta }^{\ast }$. In a Blasius boundary layer, numerical integration gives ${\delta }^{\ast }\simeq 1.720788\sqrt{\nu x/U}$. This equivalent displacement is smaller than the boundary layer thickness and is not the position of a physical surface.^[ref]

The momentum thickness weights the streamwise velocity deficit by the local velocity carrying it:Thus $\rho U^2\theta$ measures a momentum flux deficit per unit span relative to transporting the same local mass flow rate at $U$. It is not $\rho \int (U^2-u^2)dy$, which compares a different reference flux. A Blasius boundary layer has $\theta \simeq 0.664115\sqrt{\nu x/U}$; the ratio of displacement thickness to momentum thickness is about $2.591$.^[ref]

For the steady flow over an impermeable flat plate at constant $U$, integrating the boundary-layer equation and using mass conservation givesThis connects growth of the momentum thickness to local wall shear stress. Consequently the drag per unit span on one face up to $x=L$ is $D^{\prime }=\rho U^2\theta (L)$ in the ideal Blasius boundary layer with $\theta (0)=0$. Measured velocity profiles can therefore estimate integrated drag without differentiating noisy near-wall readings. The identity assumes no wall suction and zero outer pressure gradient.^[ref]

For a stationary flat plate in a Newtonian fluid, the streamwise shear stress exerted by the fluid on its upper face isThe sign here is positive when the fluid pulls the plate downstream; the plate exerts the opposite tangential traction on the fluid. In a Blasius boundary layer, ${\tau }_w\simeq 0.332057\rho U^2/\sqrt{{\mathrm{Re}}_x}$. A growing boundary layer thickness therefore accompanies decreasing local wall shear stress, because its near-wall velocity gradient becomes gentler. The formal divergence as $x\to 0$ is integrable, but that leading-edge limit lies outside the thin-layer hypothesis.^[ref]

The local skin friction coefficient normalizes wall shear stress by the outer dynamic pressure:For the ideal Blasius boundary layer on one face of length $L$, averaging the local skin friction coefficient gives ${\overline{c}}_f\simeq 1.328230/\sqrt{{\mathrm{Re}}_L}$. The factor of two follows from integrating $x^{-1/2}$, not from counting both plate faces. Drag on two equally exposed faces doubles when the reference area remains that of one face.

A steady two-dimensional laminar jet driven by force $F$ per span has conserved momentum flux $\rho \int u^{2}dy=F$. The boundary-layer pressure gradient is zero. Balancing inertia and viscosity with this flux gives width proportional to $x^{2/3}$ and centre velocity proportional to $x^{-1/3}$. The streamfunction $(F\nu x/\rho {)}^{1/3}f(\eta )$, $\eta =y[F/(\rho {\nu }^2{x}^2){]}^{1/3}$, obeys $3f^{\prime \prime \prime }+(f^{\prime }{)}^2+f{f}^{\prime \prime }=0$. The displayed profile has $\int (f^{\prime }{)}^{2}d\eta =2A^3/9$, so momentum equality fixes $A$. Merely requiring independence of $x$ does not fix the amplitude. Entrainment makes volume flux $2A(F\nu x/\rho {)}^{1/3}$ grow downstream.

A steady attached laminar boundary layer with slip speed $u_t$ and downstream distance $s$ has thickness $\delta \sim \sqrt{\nu s/u_t}$. The wall shear is ${\tau }_w\sim \mu u_t/\delta$, giving the displayed scaling, where $\mu =\rho \nu$ is dynamic viscosity. Multiplying by a nearly constant wetted perimeter gives local drag per unit length proportional to $u_t^{3/2}{s}^{-1/2}$, and integration gives drag proportional to $u_t^{3/2}\sqrt{L}$. A force density weighted by $s^{1/2}$ instead produces $L^{3/2}$; it can be a separately prescribed phenomenological model but is not this local laminar shear law. Strong curvature, unsteadiness, separation or turbulent transition can invalidate the estimate.

An axisymmetric sink inside a cone has outer speed $U=-A/x^2$. With local distance $x$ along the wall and normal distance $y$, leading continuity is $u_x+u/x+v_y=0$ because the circumference is proportional to $x$. Balancing convection with normal viscous diffusion gives the displayed thickness. With streamfunction $\Psi =-\sqrt{A\nu }{x}^{1/2}F(y/\delta )$, $xu={\Psi }_y$ and $xv=-{\Psi }_x$, the similarity equation is $F^{\prime \prime \prime }-F{F}^{\prime \prime }/2+2(1-F^{\prime 2})=0$, with $F(0)=F^{\prime }(0)=0$ and $F^{\prime }(\infty )=1$. The displacement-flux correction scales as $\sqrt{A\nu }{R}^{1/2}$, and the corresponding outer velocity correction as $\sqrt{A\nu }{R}^{-3/2}$.

In a symmetric thin radial free jet, $F=\int r{u}_r^{2}dz$ is conserved by the conservative momentum equation. Balancing advection and viscosity gives thickness $\delta =Cr$, with $C=({\nu }^2/F{)}^{1/3}$. The normalized similarity streamfunction uses $A=(F^2/\nu {)}^{1/3}$ and $f=2k\tanh (k\eta )$, $k=(3/16{)}^{1/3}$, so $\int f^{\prime 2}d\eta =1$. The displayed velocity follows, and the thin-layer requirement is $F\gg {\nu }^2$.

A steady planar viscous jet whose conserved momentum flux and convection-diffusion balance give the displayed scaling. Choosing $U{\delta }^2=\nu x$ and $U$ as centre-line speed gives streamfunction profile $f^{\prime \prime \prime }+(f{f}^{\prime \prime }+(f^{\prime }{)}^2)/3=0$, with $f(0)=f^{\prime \prime }(0)=0$, $f^{\prime }(0)=1$ and $f^{\prime }(\pm \infty )=0$. Integration yields $f=\sqrt{6}\tanh (\eta /\sqrt{6})$ and velocity profile $f^{\prime }={\mathrm{sech}}^2(\eta /\sqrt{6})$. Its momentum integral is $4\sqrt{6}/3$, fixing the dimensional prefactors. Transverse entrainment does not require $f(\infty )=0$.

The leading momentum and mass equations in a thin laminar shear region with flow predominantly along a surface or jet axis. Longitudinal viscous derivatives are smaller than transverse ones, and pressure is supplied by the outer flow. The transverse velocity remains necessary for mass conservation and entrainment even when it is small compared with the axial velocity.

For a fluid occupying $z<0$, match an outer stagnation point flow $u=x,w=-z$ to a flat impermeable interface with dimensionless shear stress $\beta x$. With boundary layer thickness $\delta ={\mathrm{Re}}^{-1/2}$, inner variable $\zeta =z/\delta$ and $\gamma =\delta \beta$, the leading boundary-layer equation admits $u=x{f}^{\prime }(\zeta )$ and $w=-\delta f(\zeta )$. Incompressibility is automatic, and tangential momentum conservation gives the displayed ordinary differential equation, with $f(0)=0$, $f^{\prime \prime }(0)=\gamma$ and $f^{\prime }\to 1$ at $\zeta \to -\infty$. The sign of this far-field endpoint is fixed by the fluid domain. This is a prescribed-traction interface; it does not impose a no-slip boundary condition.

For the stress-driven stagnation-point boundary layer on $\zeta <0$, a regular perturbation $f=\zeta +\gamma f_1+O({\gamma }^2)$ gives $g=f_1^{\prime }$ satisfying $g^{\prime \prime }+\zeta g^{\prime }-2g=0$, $g^{\prime }(0)=1$ and $g(-\infty )=0$. Reduction of order with the known solution $1+{\zeta }^2$ givesIndeed, writing $g=(1+{\zeta }^2)v$ gives $v^{\prime }=C{e}^{-{\zeta }^2/2}/(1+{\zeta }^2{)}^2$; decay at $-\infty$ fixes the lower endpoint and $g^{\prime }(0)=1$ fixes $C=1$. The integral at zero is $\sqrt{2\pi }/4=\sqrt{\pi /8}$, as follows by differentiating the antiderivativeHere $\mathrm{erf}$ is the error function. Thus positive imposed outward shear stress increases the outward surface velocity by the displayed first-order factor. Replacing decay at $-\infty$ by decay at $+\infty$ while leaving the equation and wall derivative unchanged reverses this coefficient and describes a different problem.

A boundary layer in which turbulent velocity fluctuations transport momentum and scalars. A near-wall constant-stress region supports the law of the wall; buoyancy modifies it through Monin-Obukhov similarity theory. A viscous or roughness sublayer and an outer-layer scale bound the range of either approximation.

The near-wall region where molecular momentum diffusion dominates turbulent transport. Using the shear velocity $u_{\ast }=\sqrt{{\tau }_d/\rho }$ gives the viscous wall length $\nu /u_{\ast }$. The precise extent depends on the chosen crossover criterion and wall condition. In this region constant stress and no slip give $U\simeq u_{\ast }^{2}z/\nu$; the logarithmic law of the wall applies outside it.

A local similarity description for approximately stationary homogeneous surface-layer turbulence with nearly constant momentum and buoyancy fluxes. Besides height and friction velocity, the Obukhov length enters as a stability parameter. The hypothesis fixes the dimensional shear prefactor and function argument, while empirical or dynamical closures determine the functions themselves.

If the similarity functions of a turbulent boundary layer are differentiable at zero and have finite positive neutral values, expanding $f=a+dz/\ell +O((z/\ell {)}^2)$ gives a shear proportional to $1/z-d/(a\ell )$. Integration yields a logarithmic term plus a linear correction. Dimensional reasoning alone does not fix the coefficients or guarantee this regular neutral limit. The gradient Richardson number is of order $z/\ell$, so $|z/\ell |\ll 1$ is a shear-dominated, nearly neutral regime on either side of stability.

In local continuously turbulent scaling at $z/{\ell }_O\gg 1$, dependence on the distance to the wall disappears at leading order. The velocity profile is linear and the buoyancy frequency is constant. Local energy balance gives $ϵ=(c_m-1)|B_0|$ and ${\mathrm{Ri}}_f=1/c_m$, with positive dissipation requiring $c_m>1$. The coefficients depend on a closure, and this regime should not be imposed on fully collapsed turbulence.

The signed wall-turbulence length formed from stress velocity and the vertical turbulent buoyancy flux. With upward buoyancy transport positive, stable cooling has $B_0<0$ and ${\ell }_O>0$; heating from below gives ${\ell }_O<0$. The magnitude marks the transition between shear-dominated and buoyancy-dominated scaling.

In a neutral turbulent constant-stress region, the stress velocity and wall distance give mean shear $U_z=u_{\ast }/(\kappa z)$. Its logarithmic integral uses the roughness length. The law applies outside the inner sublayer and below the outer boundary-layer scale, not all the way to a smooth wall.

The dimensionless coefficient in the logarithmic wall law and its mixing-length form $\ell =\kappa z$. Its numerical value is determined from the turbulence closure and measurements, rather than from dimensional analysis alone.

The extrapolation height at which a logarithmic wall profile vanishes. For a fully rough wall it is set mainly by the geometry; for a smooth wall it scales with $\nu /u_{\ast }$. It is not the thickness of the molecular viscous sublayer or the literal location of a no-slip surface.

A thermal boundary layer is a thin region in which temperature changes rapidly and thermal conduction remains important. Its thickness can be estimated by balancing thermal diffusion against advection or by requiring a local Rayleigh number to be near its critical value.

In a renewal model, a diffusive thermal boundary layer grows until a convective instability overturns it, after which growth restarts. A layer reaching thickness $\delta \sim \sqrt{\kappa \tau }$ has mean heat flux of order $k\Delta T/\delta$. A diffusive $t^{-1/2}$ flux averaged over one cycle differs only by an order-one coefficient. Distinct thermal and material layers need not share a renewal time.

Balancing $U^2/x$ with $\nu U/{\delta }^2$ gives $\delta \sim \sqrt{\nu x/U}$.

For steady incompressible flow with imposed outer pressure gradient,The normal momentum equation makes pressure approximately uniform across the layer.

For a sheet moving with $u(x,0)=\alpha x$ through otherwise stationary fluid, the layer thickness is $\delta =\sqrt{\nu /\alpha }$. Withthe boundary-layer equations reduce toThe exact solution is $f=1-e^{-\eta }$.

Boundary-layer separation occurs where near-wall flow detaches from a surface, commonly after an adverse pressure gradient drives the wall shear to zero and then reverses it.