Boundary layers, lift, and drag

Start with a plate in a stream

Hold a thin, smooth plate parallel to a stream moving right at 2 m/s. At the stationary surface, our no-slip model gives fluid speed zero; a short distance away, it approaches 2 m/s. How can the same flow have both speeds? Why does the plate feel a drag force even though an ideal slipping-wall model would have none?

We will first follow the changing speed across the layer, compare three positions along the plate, and distinguish wall friction from pressure drag. Then we will turn to a lifting surface. The growth/profile diagram, separation comparison and supplied pressure-load graph form the first-reading path. The final similarity and momentum-integral derivation is an optional return path.

Vectors and vector spaces and Derivatives, Taylor series, and local approximation are useful published preparation. Conservation of mass, momentum, and energy and Euler’s equation and Bernoulli’s principle provide a published return path for force and pressure accounts; Integrals and multivariable calculus helps with the final integrals. Essential terms are defined here. Potential flow and Viscosity, Reynolds number, and Navier–Stokes equations are adjacent lessons; this article does not assume that you have read them.

Keep the wall and outer flow distinct

Choose x along the stationary plate in the direction of the outer flow, and y away from its upper face. Write u for streamwise fluid velocity and v for normal velocity. No slip means that tangential fluid velocity at a contacting wall equals the wall’s velocity. No penetration means their normal velocities match. Here u(x,0) = v(x,0) = 0.

These are macroscopic continuum wall conditions for our model. They do not mean that every molecule remains permanently attached, or that all fluid near a stationary object is motionless. Rarefied gases and some microscopic interfaces need different slip descriptions.

A boundary layer is a region next to the surface where the velocity adjusts strongly between the wall value and the outer velocity U_e(x). The profile changes smoothly: it has no material membrane at its edge. Outside an attached thin layer, viscous terms may be small enough for an inviscid approximation, which omits viscous stress from the momentum balance. Inside it, the large wall-normal gradients keep viscosity important.

An inviscid wall condition prevents penetration but generally permits tangential slip. It cannot simultaneously replace the viscous no-slip profile. An outer flow need not be irrotational merely because it is approximated as inviscid; a potential-flow model makes that extra velocity-field assumption. Nor does a thin layer guarantee a negligible wake: separation can strongly alter the outer flow.

Follow the wall’s influence into the flow

Viscosity transports streamwise momentum between neighbouring regions moving at different speeds. In an orderly laminar flow of a Newtonian fluid, shear stress is proportional to the local velocity gradient: approximately τ_xy = μ ∂u/∂y in this thin layer. Dynamic viscosity μ has units Pa·s. Stress is force per area; 1 Pa = 1 N/m². Near a stationary flat wall, the positive gradient of a forward-flowing profile makes the fluid pull the plate downstream while the plate exerts the opposite force on the fluid.

Let density ρ be mass per volume and define kinematic viscosity ν = μ/ρ, with units m²/s. It sets a momentum-diffusion scale. Write U∞ for the uniform outer-stream speed. A parcel reaching distance x has a travel-time scale x/U∞. During that time, diffusion influences a wall-normal distance of order √(νx/U∞). This gives the growth scale; it is not a claim that a parcel travels at U∞ at every height.

Write δ for a thickness measure of the layer. The layer becomes thicker downstream while occupying a smaller fraction of the distance from the leading edge: δ/x is of order 1/√Re_x, where Re_x = U∞x/ν is the Reynolds number based on that distance. This dimensionless number compares inertial and viscous scales. A thinner layer can still exert appreciable shear because its speed change occurs across a shorter distance.

The normal velocity v is small compared with u in the thin-layer approximation, but it is not generally zero away from the wall. As the streamwise profile grows, mass conservation deflects the outer stream slightly away from the surface. This displacement will have a different thickness measure from the visible velocity-adjustment region.

Declare the laminar flat-plate case

Supply U∞ = 2.000 m/s, ρ = 1.200 kg/m³ and μ = 0.000018 Pa·s, so ν = 0.000015 m²/s. These are assigned constant properties, not measurements identifying a particular air sample. Use a steady two-dimensional incompressible Newtonian model over one stationary, smooth, impermeable flat face at zero incidence. The outer speed is uniform and the streamwise pressure gradient is zero.

Prescribe an attached laminar thin boundary layer. Select a strip of length L = 0.300 m and width b = 0.200 m from the ideal two-dimensional solution. The width is a force multiplier; we have not added finite side edges. At L, Re_L = U∞L/ν = 40000. This number specifies our scale comparison, not proof that a real plate remains laminar.

The Blasius similarity solution gives a common profile shape at different x when height is rescaled by √(νx/U∞). Its numerical solution gives δ_99 ≈ 4.910√(νx/U∞), τ_w ≈ 0.332057ρU∞²/√Re_x, and local skin-friction coefficient c_f,x = τ_w/q∞ ≈ 0.664115/√Re_x. Here q∞ = ρU∞²/2 = 2.400 Pa is the reference dynamic pressure, and τ_w is the downstream shear force per area on the plate.

Velocity thickness δ_99 is the height where u first reaches 99% of U∞ in this monotone profile. The 99% choice is a convention, not a sharp boundary; the velocity continues to approach U∞ above it. The growth/profile diagram uses the calculated profiles, rather than a linear interpolation from wall to outer speed.

The ideal thin-layer approximation fails sufficiently close to the mathematical leading edge, where δ/x is no longer small. Its shear formula also grows without bound as x tends to zero. Do not treat that limit as the stress on a real finite-thickness tip. Our displayed stations are downstream, and the later finite drag integral is the conventional ideal leading-edge approximation, not a resolved tip calculation.

Compare three thickness meanings

At x = 0.300 m, √(νx/U∞) = 1.500 mm. Similarity height η = y/√(νx/U∞) is height divided by this local scale; here η = 1 means y = 1.500 mm. The table samples the actual similarity profile there. Fluid at y = 1.500 mm moves at about 0.660 m/s, not at the 0.400 m/s that a straight-line profile reaching 2 m/s at 7.5 mm would suggest. The profile shape matters for both shear and transported momentum.

Displacement thickness is δ* = ∫₀∞[1 − u/U_e]dy. At a station, U_e δ* is the volume-flow deficit per unit width relative to uniform U_e over the same heights. It is the height by which an ideal uniform stream would be displaced to account for that deficit. It is not the height at which velocity reaches 99%.

Momentum thickness is θ = ∫₀∞(u/U_e)[1 − u/U_e]dy. Multiplying by ρU_e² gives a momentum-deficit measure per unit width. The factor u/U_e accounts for the amount of fluid actually crossing the station; simply squaring the velocity deficit would give a different quantity.

For this solution, δ* ≈ 1.720788√(νx/U∞) and θ ≈ 0.664115√(νx/U∞). At L, δ_99 = 7.364984 mm, δ* = 2.581181 mm and θ = 0.996172 mm. They describe three different relationships in the same profile. The thin-layer ratio δ_99/L is about 0.02455 at L; it does not measure a transition threshold or a real-device error bound.

Integrate shear and name the area

The local wall shear decreases as x⁻¹ᐟ² even though the total drag accumulated from the leading edge increases. Integrate over one face: D_f = b∫₀ᴸτ_w dx ≈ 0.664115ρU∞²bL/√Re_L. For our strip, D_f = 0.000956325 N = 0.956325 mN, downstream on the plate. The wall exerts the opposite force on the fluid.

The integral is finite because ∫₀ᴸx⁻¹ᐟ²dx = 2√L. It uses the ideal leading-edge asymptotic solution over the strip; it does not remove its local breakdown or describe a finite tip. Independent identical flow on two faces would double the friction force. A finite body’s additional edge and pressure forces would need their own model.

For one face, the wetted reference area is S_w = bL = 0.0600 m². Average skin-friction coefficient is C_f = D_f/(q∞S_w) ≈ 1.328229/√Re_L = 0.00664115. At the downstream end, local c_f,L ≈ 0.00332057. The average is twice the downstream local value in this particular x⁻¹ᐟ² model, not in every boundary layer.

More generally define lift and drag coefficients by C_L = L_force/(q∞S_ref) and C_D = D_force/(q∞S_ref). State the reference area: wing planform, frontal projected area and wetted area are not interchangeable. A coefficient is a dimensionless summary of a specified flow and geometry, not a material constant.

Separate pressure and friction drag

With the incoming flow rightward, drag on the stationary body is the force component along that flow; lift is the component perpendicular to it, here chosen upward. They describe directions of the resultant, not two different kinds of fluid. In the equivalent frame of a moving body, drag opposes its motion relative to the undisturbed fluid.

Pressure acts normally on each surface element. Viscous stress supplies another contact-force contribution, commonly dominated by tangential wall shear in a thin boundary layer. Sum their actual directional components over the complete body. With n pointing from the solid into the fluid, pressure traction on the body is −pn. A large pressure at one point does not determine the resultant without surface normals and the rest of the pressure distribution.

Skin-friction drag comes from the streamwise component of viscous surface traction. Pressure drag comes from the streamwise component of the normal pressure load. On our ideal zero-thickness parallel face, pressure has no streamwise component, so the calculation gives friction alone. A blunt body can have a large pressure imbalance between its front and separated rear flow even when its skin-friction contribution is relatively small.

These components are coupled physically. Viscosity changes the profile and may trigger separation; separation changes the surrounding pressure and wake. An inviscid attached-flow pressure solution that predicts no drag does not establish zero drag on a real body. Conversely, adding a flat-plate friction estimate to an unchanged ideal pressure pattern is not generally sufficient once the flow separates.

Follow an adverse pressure gradient

Now leave the constant-speed flat-plate case. Let the attached outer velocity U_e(x) vary along a smooth surface. For a steady incompressible inviscid outer balance with no streamwise body force, U_e dU_e/dx = −(1/ρ)dp_e/dx. A decreasing outer speed therefore accompanies increasing pressure downstream: dp_e/dx > 0 is called an adverse pressure gradient. A falling pressure, dp_e/dx < 0, is favourable for forward motion.

To leading order in a thin attached layer, pressure at each x is nearly constant across its thickness and matches outer pressure p_e(x). It is not zero there. The same adverse pressure gradient acts against the slow near-wall fluid as against the faster outer fluid. Viscous transfer can supply momentum from farther out, but near-wall forward motion may still weaken, develop zero wall shear and reverse.

In a classical steady two-dimensional flow over a stationary smooth wall, incipient separation is associated with τ_w = μ(∂u/∂y)_wall tending to zero as forward flow changes toward near-wall reverse flow. After separation, a recirculating region and detached shear layer can develop. No slip still holds at the wall; fluid does not start passing through the solid or acquire a finite tangential wall speed.

The separation comparison is schematic, not a calculation of a location or pressure-gradient threshold. Zero wall shear alone is not a universal detector: stagnation, unsteadiness and three-dimensional flows need more careful criteria. An inflected profile under an adverse gradient does not by itself prove that separation has occurred.

A separated wake often prevents the pressure recovery of an attached ideal solution, producing a larger front–rear pressure imbalance and pressure drag. The attached thin-layer equations cannot be continued as a complete downstream wake solution; strong displacement, reversed flow and unsteadiness may require a coupled or full viscous-flow treatment.

Compare laminar and turbulent layers

Laminar flow has an orderly resolved velocity field. A turbulent boundary layer has fluctuating motion that transports momentum as well as molecular viscous stress. In its time-averaged momentum equation, correlations of fluctuating velocities represent this extra transport. A turbulent mean profile is usually fuller away from the wall than a comparable laminar profile; no slip still makes mean and instantaneous wall velocity match the stationary wall.

At otherwise comparable conditions, turbulent mixing commonly increases wall friction, but it can also move higher streamwise momentum toward the wall and delay separation under an adverse pressure gradient. Those changes can raise friction drag while lowering pressure drag. Which effect wins depends on the body, pressure distribution and flow conditions; turbulence is not automatically either a drag-reduction method or a drag penalty for total body force.

Transition depends on disturbances, roughness, pressure gradient, geometry and other conditions as well as Reynolds number. We have prescribed laminar flow in the numerical plate example, not assigned a universal transition value. A laminar boundary layer can border a turbulent wake, and a flow can be unsteady without being turbulent.

Streamlining can make pressure recovery more gradual and help a layer remain attached. Suction or deliberate surface disturbances can alter momentum transport and separation, but their energetic cost and friction effects also matter. Fluid instabilities and convection and Turbulence continue these questions; predicting transition or a turbulence closure is outside this lesson’s numerical model.

Connect lift to pressure and flow turning

A wing section, called an airfoil, turns the surrounding flow and develops a pressure distribution. Its chord is the line from leading to trailing edge; angle of attack is the angle between that chord and the undisturbed incident flow. Camber is the curvature of the section’s mean line. A symmetric section or a thin flat plate can generate lift at nonzero incidence: unequal upper and lower path lengths are not required.

When flow speed is small compared with the speed of sound, an appropriate steady incompressible inviscid outer approximation lets Bernoulli’s relation connect local outer speed and pressure along a streamline. If a common Bernoulli constant is justified, faster outer flow has lower pressure at the same height. This relation does not determine the velocity distribution from shape alone, and it must not be applied across a dissipative boundary layer as though nothing changed.

There is no rule that fluid parcels split at the leading edge must meet again at the trailing edge. Equal transit time is an invented constraint, not mass conservation. The complete flow and its boundary conditions determine speeds and paths. Pressure on both upper and lower surfaces contributes to lift; a description involving only air bouncing off the underside also misses essential flow turning.

A surface-force account integrates pressure and viscous traction to obtain the force on the wing. A momentum account obtains the force on the fluid from changes in momentum and takes the opposite wing force. For a steady control volume, include momentum flux and pressure/other boundary forces consistently; for an unsteady volume, include momentum storage too. Downward turning and an upward wing force are compatible parts of the same balance, not competing explanations.

Optional circulation return path: Vorticity and circulation is published preparation for this paragraph. In a steady incompressible two-dimensional attached high-Reynolds-number approximation, circulation helps represent the lifting outer flow. Circulation Γ_ccw = ∮v·dl is the velocity line integral around the section, with units m²/s and positive counterclockwise orientation. For a rightward uniform stream, the ideal Kutta–Joukowski relation is L′_y = −ρU∞Γ_ccw, where L′ is lift per span. A sharp trailing edge adds the Kutta condition: the selected outer solution leaves that edge smoothly with finite velocity. Starting flow, viscosity and the wake help establish the physical circulation; an arbitrary added swirl is not a prediction of wing lift.

These attached-flow approximations do not predict stall by themselves. Separation can substantially change the pressure distribution and reduce lift as incidence increases. The stall condition is not one universal angle. Finite-span wings also create trailing vortices and downwash, with a drag component associated with producing lift, called induced drag. It is part of the same total surface-force/momentum account, not a third independent traction added on top of pressure and viscosity.

Integrate a separately supplied lift load

Use a separate illustrative load, not the 2 m/s friction plate. Represent two faces by one horizontal projected rectangle with chord c = 0.250 m and span b = 0.400 m. Supply pressure difference Δp(x) = p_lower − p_upper = (30 Pa)(1 − x/c), uniform across the span. It decreases linearly from 30 Pa at the leading edge to zero at the trailing edge.

Take p_ref = 100000 Pa and set p_lower = p_ref + Δp/2, p_upper = p_ref − Δp/2. Both pressures are positive. Pressure on the lower face pushes upward and pressure on the upper face pushes downward. The common reference cancels between the two faces; their difference determines the net upward load.

The lift from this supplied thin-surface pressure load is L_force = b∫₀ᶜΔp dx = bcΔp₀/2 = 1.500 N. The pressure-difference graph is a triangle of average height 15 Pa; multiplying by projected area bc = 0.100 m² gives the same force. The largest pressure difference alone times the entire area would overestimate the force by two.

For reference density 1.200 kg/m³ and speed 10.000 m/s, q∞ = 60.000 Pa and C_L = 1.500/[(60.000)(0.100)] = 0.250. The pressure load has zero streamwise component in this horizontal projected-face idealisation. That does not predict zero drag for a real lifting wing.

The pressures are supplied, not measured or solved from an airfoil geometry. Thickness, surface slopes, viscous traction, finite-span flow and induced drag are not included in this load calculation. A real lift prediction would have to determine the flow and pressure first; this example teaches how a known distributed pressure load becomes a resultant.

Optional return: the thin-layer equations

For the declared steady two-dimensional incompressible laminar layer with constant μ and ρ, mass conservation is ∂u/∂x + ∂v/∂y = 0. The leading streamwise momentum equation is u∂u/∂x + v∂u/∂y = −(1/ρ)dp_e/dx + ν∂²u/∂y². The normal balance gives ∂p/∂y ≈ 0 at this order. Apply u = v = 0 at the stationary wall and match u to U_e outside the layer.

Use a streamfunction ψ, defined by u = ∂ψ/∂y and v = −∂ψ/∂x, so continuity is automatic. For uniform U∞, choose η = y√(U∞/(νx)) and ψ = √(νU∞x)f(η). Then u = U∞f′(η) and v = (1/2)√(νU∞/x)[ηf′ − f]. Primes here mean derivatives with respect to η, not x.

Substitution into the zero-pressure-gradient momentum equation gives f‴ + (1/2)ff″ = 0, with f(0) = f′(0) = 0 and f′(∞) = 1. Numerically solving this boundary-value problem gives f″(0) ≈ 0.332057 and f′(4.909990) = 0.99. The shear follows by differentiating u in y: τ_w = μU∞f″(0)√(U∞/(νx)). The displayed thickness and friction constants therefore share one solved profile.

For a more general smooth attached layer, integrating momentum gives dθ/dx + [(2θ + δ*)/U_e]dU_e/dx = τ_w/(ρU_e²). With constant U_e = U∞ this reduces to ρU∞²dθ/dx = τ_w. Integrating from the ideal leading edge gives D_f/b = ρU∞²θ(L), matching the wall-shear calculation. With a pressure gradient, the additional term must be retained; one cannot reuse the constant-speed drag relation unchanged.

These are approximations conditioned on a thin attached layer and the stated continuum model, not replacements for all Navier–Stokes solutions. Strong separation, significant compressibility, rapid unsteadiness, three-dimensional effects or a non-Newtonian stress law change the problem. The neighbouring permanent subjects develop those regimes rather than extending the numerical flat-plate case beyond its assumptions.

Worked examples

Try it yourself

1. In the supplied laminar model compare x = 0.030 m with x = 0.120 m. What happens to δ_99, local τ_w and Re_x? Does total accumulated friction decrease?

   Show the answer and reasoning

   δ_99 doubles from 2.329013 to 4.658025 mm; τ_w halves from about 0.025201 to 0.012601 Pa; Re_x quadruples from 4000 to 16000. Accumulated ideal drag from the leading edge doubles, not decreases.

   Distance increases by four, so √x doubles and x⁻¹ᐟ² halves. Reynolds number is proportional to x. Integrating the positive shear gives D_f(0→x) proportional to √x: a smaller local addition can still make a larger positive total.

2. Set U∞ = 4.000 m/s and retain every other supplied flat-plate input and the assumed laminar regime. Compare Re_L, δ_99(L), τ_w(L) and one-face drag with the original case. Are you predicting transition?

   Show the answer and reasoning

   Re_L = 80000; δ_99(L) ≈ 5.207830 mm; τ_w(L) ≈ 0.022540799 Pa; one-face drag ≈ 2.704896 mN. No transition prediction is made.

   At fixed x and ν, thickness scales as U∞⁻¹ᐟ². Shear and the integrated friction scale as U∞³ᐟ², so doubling speed multiplies them by 2√2, not four. The regime remains a condition of the calculation; disturbances and transition have not been modelled.

3. Double μ to 0.000036 Pa·s at fixed ρ, U∞, L and b, keeping the same assumed laminar flat-plate model. Find Re_L, δ_99(L) and one-face drag. Can a thicker layer exert more friction?

   Show the answer and reasoning

   Re_L = 20000; δ_99(L) ≈ 10.415661 mm; one-face drag ≈ 1.352448 mN. The thicker layer has greater friction in this comparison.

   Kinematic viscosity doubles, so thickness increases by √2 and Reynolds number halves. The wall gradient decreases by 1/√2, but stress includes μ times that gradient. Doubling μ therefore raises stress and drag by √2. Thickness alone does not determine friction.

4. Assume two independent identical wetted faces, each governed by the supplied plate model. Find total friction drag. What coefficient results if the reference area is total wetted area, and what if it is the original one-face planform area?

   Show the answer and reasoning

   Total friction is approximately 1.912650 mN. With total wetted area 0.1200 m² the coefficient is about 0.00664115; with one-face planform area 0.0600 m² it is about 0.01328229.

   Doubling the face count doubles force. Normalizing by an area that also doubles leaves the coefficient unchanged; keeping the smaller reference area doubles it. Both numbers describe the same force if their definitions are explicit. No additional finite-body edge force is included.

5. A schematic steady two-dimensional separation progression has wall gradients +20, 0 and −5 s⁻¹ while the stationary wall has u = 0 in all three profiles. Which state has near-wall reverse flow? Does zero gradient alone prove separation in every flow?

   Show the answer and reasoning

   The negative-gradient profile has upstream velocity just away from the wall; the positive-gradient profile has forward velocity there. The zero-gradient state marks incipient separation only in the declared forward-to-reverse progression, not in every possible flow.

   Near a smooth wall, u(y) begins with its wall slope times y unless that first slope vanishes. Negative slope therefore gives negative near-wall velocity while u(0) remains zero. Separation does not violate no slip. Stagnation, unsteady or three-dimensional situations show why a zero-shear value alone is insufficient as a universal criterion.

6. Double the supplied pressure difference in the separate lift case while keeping chord, span, reference density and 10.000 m/s reference speed fixed. Find lift and C_L. Does this change imply doubled upper-flow speed or equal parcel transit times?

   Show the answer and reasoning

   Lift doubles to 3.000 N and C_L doubles to 0.500. Neither doubled upper-flow speed nor equal transit times follows.

   The pressure-load integral is linear in Δp, and q∞S_ref is unchanged. No velocity field was supplied or solved, so this integration cannot determine upper and lower speeds. A Bernoulli closure would need its own valid flow assumptions; equal transit time would remain an unjustified constraint.

What to carry forward

• No slip fixes fluid velocity at the wall; viscous momentum transport creates a profile that adjusts toward the outer flow without a sharp material edge.
• A thin attached layer can retain strong viscous gradients beside an approximately inviscid outer region. Inviscid and irrotational are different assumptions.
• Velocity, displacement and momentum thickness describe different properties of one profile. The supplied laminar flat-plate solution grows as √x and its local shear decreases as x⁻¹ᐟ².
• Integrate shear over the actual wetted faces and state the area used by a force coefficient. The ideal leading-edge approximation is not a finite-tip solution.
• An adverse pressure gradient can weaken forward near-wall motion and lead to separation, changed wake pressure, pressure drag and reduced lift.
• Transition is conditional, not fixed by one universal Reynolds number. More turbulent friction can coexist with delayed separation and lower total drag on a particular body.
• Lift from surface pressure and lift from momentum change are consistent accounts of one interaction. Neither requires parcels above and below a wing to reunite at the trailing edge.
• A supplied pressure distribution can be integrated into lift, but predicting that distribution, stall or finite-span induced drag needs an appropriate flow model.