4.2 Boundary layer equations

Boundary layer concept

Boundary layer equations describe fluid behavior in the thin region near a solid surface where viscous effects dominate. They're simplified forms of the Navier-Stokes equations, and they're essential for predicting drag, heat transfer, and flow separation on aerodynamic surfaces.

The boundary layer is a thin region near a solid surface where the velocity transitions from zero at the surface (due to the no-slip condition) to the freestream value. Ludwig Prandtl introduced this concept in 1904, and it fundamentally changed fluid mechanics by showing that you could split the flow into two zones: a thin viscous region near the wall, and an essentially inviscid region everywhere else.

Viscous effects near surfaces

Near solid surfaces, fluid particles adhere to the surface due to the no-slip condition, creating a velocity gradient normal to the surface. Viscous forces dominate within the boundary layer, while outside it the flow can be treated as inviscid.

These viscous effects produce two important consequences:

• Skin friction drag from shear stress at the wall
• Heat transfer between the fluid and the surface, driven by the temperature gradient in the boundary layer

Boundary layer thickness

The boundary layer thickness $\delta$ is defined as the distance from the surface where the velocity reaches 99% of the freestream velocity $U_\infty$. This is somewhat arbitrary (why 99% and not 98%?), but it's the standard convention.

The boundary layer starts with zero thickness at the leading edge and grows along the surface in the flow direction. Three main factors influence how quickly it grows:

• Reynolds number: Higher $Re$ means a thinner boundary layer relative to the body length
• Surface roughness: Rougher surfaces promote faster boundary layer growth
• Pressure gradient: Adverse pressure gradients cause the boundary layer to thicken more rapidly

Displacement thickness

The displacement thickness $\delta^*$ quantifies how much the boundary layer pushes the external inviscid flow outward. Think of it this way: because the boundary layer slows the fluid near the wall, there's a mass flow deficit compared to what you'd get with uniform freestream velocity. The displacement thickness is the effective distance the surface would need to move outward to account for that deficit.

$\delta^* = \int_0^{\delta} \left(1 - \frac{u}{U_{\infty}}\right) dy$

Momentum thickness

The momentum thickness $\theta$ measures the momentum deficit in the boundary layer relative to inviscid flow. Where displacement thickness captures the mass deficit, momentum thickness captures the momentum deficit, making it directly useful for drag calculations.

$\theta = \int_0^{\delta} \frac{u}{U_{\infty}} \left(1 - \frac{u}{U_{\infty}}\right) dy$

The ratio $H = \delta^*/\theta$ is called the shape factor and is a useful indicator of whether the boundary layer is close to separating (higher $H$ values signal greater separation risk).

Boundary layer equations

The governing equations for boundary layer flows come from conservation of mass and momentum. They're simplified versions of the full Navier-Stokes equations, tailored to exploit the thinness of the boundary layer.

Continuity equation

The continuity equation enforces conservation of mass. For two-dimensional, incompressible boundary layer flow:

$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$

This says that if the streamwise velocity $u$ changes in the $x$-direction, there must be a compensating change in the wall-normal velocity $v$ in the $y$-direction so that mass is conserved.

Navier-Stokes equations

The full Navier-Stokes equations govern viscous fluid flow and come from applying Newton's second law to a fluid element. For two-dimensional, incompressible flow:

• x-momentum: $\rho \left(\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y}\right) = -\frac{\partial p}{\partial x} + \mu \left(\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2}\right)$
• y-momentum: $\rho \left(\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y}\right) = -\frac{\partial p}{\partial y} + \mu \left(\frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2}\right)$

These are exact but very difficult to solve in general. The boundary layer approximations make them tractable.

Boundary layer approximations

The key insight is that the boundary layer is thin. This geometric fact lets you drop several terms from the Navier-Stokes equations. The standard approximations are:

1. The flow is steady and two-dimensional
2. The boundary layer thickness $\delta$ is much smaller than the characteristic length $L$ of the surface ($\delta \ll L$)
3. The wall-normal velocity is much smaller than the streamwise velocity ($v \ll u$)
4. The pressure gradient across the boundary layer is negligible: $\frac{\partial p}{\partial y} \approx 0$, so pressure is "imposed" by the outer inviscid flow
5. Streamwise diffusion is negligible compared to wall-normal diffusion: $\frac{\partial^2 u}{\partial x^2} \ll \frac{\partial^2 u}{\partial y^2}$

These approximations hold for high Reynolds number flows with thin boundary layers.

Simplified boundary layer equations

Applying the approximations above reduces the Navier-Stokes equations to the Prandtl boundary layer equations:

• x-momentum: $u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = -\frac{1}{\rho} \frac{dP}{dx} + \nu \frac{\partial^2 u}{\partial y^2}$
• Continuity: $\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$

Notice what changed compared to the full Navier-Stokes equations: the time derivative is gone (steady flow), the $\frac{\partial^2 u}{\partial x^2}$ term is gone (streamwise diffusion negligible), and the y-momentum equation has collapsed to $\frac{\partial p}{\partial y} = 0$. The pressure gradient $\frac{dP}{dx}$ is now an ordinary derivative because pressure doesn't vary across the boundary layer; it's determined entirely by the outer inviscid solution.

Boundary conditions:

• At the wall ($y = 0$): $u = 0$, $v = 0$ (no-slip, no penetration)
• At the boundary layer edge ($y \to \delta$): $u \to U_\infty(x)$

Laminar boundary layers

Laminar boundary layers feature smooth, orderly flow with no mixing between fluid layers. They occur at lower Reynolds numbers and produce less skin friction drag than turbulent boundary layers, but they're also more susceptible to separation.

Blasius solution

The Blasius solution is the classic exact solution for a laminar boundary layer over a flat plate with zero pressure gradient ($dP/dx = 0$).

The approach uses a similarity transformation. Blasius introduced the similarity variable:

$\eta = y \sqrt{\frac{U_\infty}{\nu x}}$

and a dimensionless stream function $f(\eta)$ such that the velocity profile at every streamwise station collapses onto a single curve. This transforms the two partial differential equations (momentum + continuity) into a single third-order nonlinear ODE:

$f''' + \frac{1}{2} f f'' = 0$

This ODE is solved numerically (it has no closed-form solution). Key results from the Blasius solution:

• Boundary layer thickness: $\delta \approx \frac{5.0 x}{\sqrt{Re_x}}$
• Displacement thickness: $\delta^* \approx \frac{1.72 x}{\sqrt{Re_x}}$
• Momentum thickness: $\theta \approx \frac{0.664 x}{\sqrt{Re_x}}$
• Local skin friction coefficient: $C_f = \frac{0.664}{\sqrt{Re_x}}$

Falkner-Skan solution

The Falkner-Skan solution generalizes Blasius to include a pressure gradient. It assumes the freestream velocity varies as a power law: $U_\infty(x) \propto x^m$, which introduces a pressure gradient parameter $\beta = \frac{2m}{m+1}$.

• $\beta > 0$: favorable pressure gradient (accelerating flow)
• $\beta = 0$: zero pressure gradient (recovers the Blasius solution)
• $\beta < 0$: adverse pressure gradient (decelerating flow)

The Falkner-Skan equation is solved numerically for different values of $\beta$. As $\beta$ becomes more negative, the velocity profile near the wall becomes less full, and at a critical value the wall shear stress reaches zero, indicating the onset of separation.

Thermal boundary layer

When there's a temperature difference between the surface and the fluid, a thermal boundary layer develops alongside the velocity boundary layer. It's the region where the temperature transitions from the surface value $T_w$ to the freestream temperature $T_\infty$.

The thermal boundary layer thickness $\delta_T$ is defined as the distance where the temperature reaches 99% of the freestream temperature. The ratio of the velocity to thermal boundary layer thickness depends on the Prandtl number $Pr = \nu / \alpha$:

• $Pr > 1$ (most liquids): $\delta_T < \delta$ (thermal boundary layer is thinner)
• $Pr < 1$ (liquid metals): $\delta_T > \delta$ (thermal boundary layer is thicker)
• $Pr \approx 1$ (most gases): $\delta_T \approx \delta$

Heat transfer in laminar flow

Heat transfer in laminar boundary layers occurs through conduction and convection. The local heat transfer is characterized by the Nusselt number $Nu_x$, which is the ratio of convective to conductive heat transfer.

For flow over a flat plate with constant surface temperature:

$Nu_x = \frac{h_x x}{k} = 0.332\, Re_x^{1/2}\, Pr^{1/3}$

where $h_x$ is the local heat transfer coefficient and $k$ is the fluid thermal conductivity. This correlation is valid for $Pr \gtrsim 0.6$. The average Nusselt number over a plate of length $L$ is obtained by integrating, giving $\overline{Nu}_L = 0.664\, Re_L^{1/2}\, Pr^{1/3}$.

Turbulent boundary layers

Turbulent boundary layers feature chaotic, fluctuating flow with intense mixing between fluid layers. They occur at higher Reynolds numbers and are far more common in real engineering applications. Compared to laminar boundary layers, turbulent ones have higher skin friction drag and higher heat transfer rates, both due to the enhanced mixing.

Transition from laminar to turbulent

The transition from laminar to turbulent flow occurs when the local Reynolds number $Re_x = \frac{U_\infty x}{\nu}$ exceeds a critical value, typically around $5 \times 10^5$ for a flat plate in a low-disturbance environment.

The transition process involves several stages:

1. Small disturbances (Tollmien-Schlichting waves) grow in the laminar boundary layer
2. These disturbances become nonlinear and form three-dimensional structures
3. Turbulent spots appear and grow as they convect downstream
4. The spots merge into a fully turbulent boundary layer

Factors that promote earlier transition include surface roughness, higher freestream turbulence intensity, and adverse pressure gradients. Favorable pressure gradients tend to delay transition.

Turbulent velocity profile

The turbulent velocity profile is much fuller near the wall than the laminar profile, meaning more momentum is carried close to the surface. It consists of three distinct regions:

• Viscous sublayer ($y^+ < 5$): Velocity varies linearly with distance from the wall ($u^+ = y^+$). Viscous stresses dominate, and turbulent fluctuations are damped.
• Buffer layer ($5 < y^+ < 30$): A transition zone where both viscous and turbulent stresses are significant. Neither the linear nor the log law applies cleanly here.
• Log-law region ($30 < y^+ < 500$): Turbulent mixing dominates, and the velocity follows a logarithmic profile.

Here $y^+ = \frac{y u_*}{\nu}$ and $u^+ = \frac{u}{u_*}$ are the wall-normal distance and velocity in wall units, non-dimensionalized by the friction velocity $u_* = \sqrt{\tau_w / \rho}$.

Logarithmic law of the wall

In the log-law region, the velocity profile follows:

$\frac{u}{u_*} = \frac{1}{\kappa} \ln \left(\frac{y u_*}{\nu}\right) + B$

where $\kappa \approx 0.41$ is the von Kármán constant and $B \approx 5.0$ for a smooth wall. For rough walls, $B$ decreases depending on the roughness height relative to the viscous sublayer thickness.

The log law is valid for roughly $30 < y^+ < 500$. Beyond $y^+ \approx 500$, the velocity profile transitions into the wake region, where it deviates from the log law and approaches the freestream velocity.

Turbulent heat transfer

Turbulent mixing enhances heat transfer significantly compared to laminar flow. The Stanton number $St$ characterizes the heat transfer performance:

$St_x = \frac{h_x}{\rho c_p U_{\infty}} = 0.0296\, Re_x^{-1/5}\, Pr^{-2/3}$

This correlation applies to flow over a flat plate with constant surface temperature. The Reynolds analogy connects skin friction and heat transfer in turbulent flow: higher wall shear stress generally correlates with higher heat transfer rates.

Boundary layer separation

Boundary layer separation occurs when fluid near the wall can no longer push against an increasing pressure (adverse pressure gradient) and detaches from the surface. This creates a recirculating flow region that dramatically increases drag and can alter the entire flow field around a body.

Adverse pressure gradient

An adverse pressure gradient exists when pressure increases in the flow direction ($dP/dx > 0$). The rising pressure acts like a headwind on the slow-moving fluid near the wall, decelerating it further.

If the adverse pressure gradient is strong enough, the near-wall fluid decelerates to zero velocity and can even reverse direction. At that point, the boundary layer separates from the surface.

Separation point

The separation point is where the boundary layer detaches. It's defined by the condition:

$\tau_w = \mu \left.\frac{\partial u}{\partial y}\right|_{y=0} = 0$

At this location, the velocity gradient at the wall is zero. Downstream of this point, the flow near the wall reverses direction. The separation point location depends on the pressure distribution, Reynolds number, and whether the boundary layer is laminar or turbulent. Turbulent boundary layers resist separation better than laminar ones because turbulent mixing brings high-momentum fluid closer to the wall.

Separated flow regions

Downstream of the separation point, a recirculation zone forms with reversed flow near the surface. This region is highly unsteady, featuring vortex shedding and intense turbulent mixing. The size of the separated region depends on the Reynolds number, the strength of the adverse pressure gradient, and the body geometry.

In some cases, the separated flow can reattach downstream, forming a separation bubble. This is common on airfoils at moderate angles of attack.

Effects of separation on drag

Separation dramatically increases pressure drag (also called form drag). The recirculation zone creates a low-pressure region behind the body, so there's a large pressure imbalance between the front (high pressure) and the rear (low pressure).

For streamlined bodies like airfoils, separation can cause a sudden and severe increase in drag along with loss of lift (stall). For bluff bodies like cylinders, separation is the dominant source of drag. Controlling or delaying separation is one of the primary goals of aerodynamic design.

Boundary layer control

Boundary layer control techniques manipulate the boundary layer to delay separation, reduce drag, or enhance heat transfer. These methods fall into two categories: active methods (requiring external energy input) and passive methods (using geometric features with no energy input).

Suction vs blowing

Suction removes low-momentum fluid from near the wall through small holes or porous surfaces. This thins the boundary layer and makes it more resistant to separation. It's effective but requires a pumping system.

Blowing (or tangential injection) adds high-momentum fluid near the surface, energizing the boundary layer. This is commonly used on aircraft flaps and engine nacelles.

The effectiveness of both methods depends on the location along the surface, the flow rate, and the distribution of suction/blowing slots. Suction is generally more effective at delaying transition, while blowing is better at preventing separation in regions with strong adverse pressure gradients.

Vortex generators

Vortex generators are small vanes or tabs mounted on a surface that create streamwise vortices. These vortices mix high-momentum fluid from the outer boundary layer down toward the wall, energizing the near-wall region and delaying separation.

They're a passive method, so they require no energy input, but they do add a small amount of parasitic drag. You'll find them on aircraft wings (especially near the trailing edge), wind turbine blades, and inside diffusers. Their placement is critical: too far upstream and the vortices dissipate before reaching the separation-prone region; too far downstream and separation has already occurred.

Riblets and surface roughness

Riblets are tiny streamwise grooves (typically on the order of tens of microns) that reduce turbulent skin friction drag by 5-8%. They work by restricting the spanwise motion of near-wall streamwise vortices, which reduces turbulent momentum transfer to the wall.

Surface roughness serves a different purpose. Roughness elements (sand grains, dimples, trip strips) promote earlier transition from laminar to turbulent flow. While this increases skin friction, the resulting turbulent boundary layer is more resistant to separation. This is the principle behind dimples on a golf ball: the turbulent boundary layer stays attached longer, reducing the wake size and overall drag.

Laminar flow control techniques

Laminar flow control aims to maintain laminar flow over as much of the surface as possible, since laminar skin friction is much lower than turbulent. Techniques include:

• Surface shaping: Designing the pressure distribution to maintain a favorable pressure gradient over a large portion of the surface (natural laminar flow, or NLF)
• Suction: Removing boundary layer fluid to keep the flow stable and laminar
• Active wave cancellation: Introducing disturbances that destructively interfere with growing instabilities

These techniques are most valuable in high-speed applications (commercial aircraft wings, for example) where skin friction is a major fraction of total drag. The challenge is that laminar flow is fragile: surface contamination (insects, ice) or manufacturing imperfections can trigger premature transition.

Computational methods

Computational methods solve the boundary layer equations (or the full Navier-Stokes equations) numerically. They range from fast approximate methods to extremely expensive but highly accurate simulations. The right choice depends on the flow regime, geometry complexity, and the accuracy you need.

Integral boundary layer methods

Integral methods don't solve for the full velocity profile. Instead, they assume a parametric profile shape (polynomial, power law, etc.) and solve for integrated quantities like $\delta^*$ and $\theta$ using the von Kármán integral momentum equation:

$\frac{d\theta}{dx} + (2 + H)\frac{\theta}{U_\infty}\frac{dU_\infty}{dx} = \frac{C_f}{2}$

where $H = \delta^*/\theta$ is the shape factor. These methods are fast and give reasonable estimates of boundary layer development, but they can't capture detailed profile shapes or complex separation behavior.

Finite difference schemes

Finite difference methods discretize the boundary layer equations on a grid. The derivatives are replaced with algebraic approximations (central differences, upwind schemes, etc.), producing a system of equations solved at each grid point.

For boundary layer equations, a common approach is the Crank-Nicolson scheme, which marches the solution downstream station by station. This is efficient because the boundary layer equations are parabolic in the streamwise direction. Body-fitted grids allow these methods to handle curved surfaces.

Turbulence modeling for boundary layers

Directly resolving all turbulent scales is impractical for most applications, so turbulence models approximate the effect of turbulent fluctuations on the mean flow. Common approaches for boundary layers:

• Algebraic models (e.g., Baldwin-Lomax): Simple, fast, but limited to attached flows
• One-equation models (e.g., Spalart-Allmaras): Solve one transport equation for a turbulence variable; widely used in aerospace
• Two-equation models (e.g., $k$-$\varepsilon$, $k$-$\omega$ SST): Solve transport equations for turbulence kinetic energy and a dissipation-related quantity; the $k$-$\omega$ SST model is particularly good for boundary layers with adverse pressure gradients
• Reynolds Stress Models (RSM): Solve transport equations for each component of the Reynolds stress tensor; more accurate for complex flows but computationally expensive
• Large Eddy Simulation (LES): Resolves large-scale turbulent structures and models only the smallest scales; much more expensive than RANS but captures unsteady separation and transition better

Direct numerical simulation (DNS)

DNS solves the full Navier-Stokes equations without any turbulence modeling, resolving every scale of motion from the largest energy-containing eddies down to the smallest Kolmogorov scales. The grid must be fine enough to capture these smallest scales, which means the computational cost scales roughly as $Re^3$.

This makes DNS feasible only for low Reynolds numbers (typically $Re_\theta < 5000$ or so with current computing power) and simple geometries. DNS is primarily a research tool used to study turbulence physics and to generate benchmark data for validating turbulence models.

Experimental techniques

Experimental techniques measure and visualize boundary layer flows in real conditions. They provide validation data for computational methods and reveal flow physics that simulations might miss. The choice of technique depends on what you need to measure (velocity, temperature, wall shear stress) and the required spatial and temporal resolution.

Hot-wire anemometry

Hot-wire anemometry measures instantaneous velocity by exploiting the cooling effect of the flow on a heated wire. A thin wire (typically tungsten or platinum, about 5 μm in diameter) is heated electrically and placed in the flow. As the flow velocity increases, more heat is convected away from the wire, changing its temperature and electrical resistance.

Two operating modes are common:

• Constant Temperature Anemometry (CTA): A feedback circuit adjusts the current to keep the wire at a fixed temperature. The required current is related to the flow velocity.
• Constant Current Anemometry (CCA): The current is held fixed, and the wire temperature (and thus resistance) varies with velocity.

Hot-wire anemometry offers excellent temporal resolution (up to hundreds of kHz), making it ideal for measuring turbulent fluctuations in boundary layers. Its main limitations are that it's intrusive (the probe must be placed in the flow) and it can be fragile in high-speed or particle-laden flows.