Wall Shear Stress
When the velocity profile in the boundary layer is known, the wall shear stress may be evaluated directly from the velocity gradient at the wall.
(29.1a)
The local skin friction coefficient is defined as the ratio of wall shear stress to the dynamic pressure based on the free-stream velocity.
- Substituting from (29.1a) we get the expression for the skin friction coefficient as shown below.
(Skin Friction Coefficient) (29.1B)
- In 1951, Liepmann and Dhawan measured the shearing stress on a flat plate directly. Their results showed a striking confirmation of Eq. (29.1).
- The total frictional force per unit width for a flat plate of length L (measured in the streamwise direction) is obtained by integrating the wall shear stress along the plate.
or
or
and the average skin friction coefficient over the plate length L is
(29.3)
where, Re = U∞ L/ν
For a flat plate of length L in the streamwise direction and width w perpendicular to the flow, the drag D due to skin friction is
(29.4)
Boundary-Layer Thickness
- Because the velocity approaches the free stream value asymptotically, it is customary to define the boundary-layer thickness δ as the distance from the wall where the local velocity reaches 99% of the free stream value, i.e. where u/U∞ ≈ 0.99.
- From Table 28.1, u/U∞ reaches 0.99 at η = 5.0, and thus one may write the practical relation between similarity coordinate η and boundary-layer thickness δ.
η = 5.0 corresponds to u/U∞ = 0.99.
Note that the 0.99 criterion for δ is arbitrary; a more physically meaningful measure is given by the displacement thickness.
- The displacement thickness δ* is defined as the distance by which the external potential flow is displaced outward because of the reduction of velocity inside the boundary layer. Graphically and mathematically it represents the area (per unit span) by which the mass flow deficit causes the effective cross-section available to the outer flow to change.
Fig. 29.1 (a) Displacement thickness (b) Momentum thickness
Displacement thickness (δ*) : It is defined as the distance by which the external potential flow is displaced outwards due to the decrease in velocity in the boundary layer.
- Substituting the similarity expressions for u/U∞ and n from Eqs. (28.21a) and (28.19) into the defining integral for δ* yields the relation shown below.
(29.7)
Following the same reasoning, a momentum thickness is defined to quantify the momentum deficit in the boundary layer relative to the potential flow.
Momentum thickness (δ) : It is defined as the thickness which, when multiplied by the free stream momentum per unit volume, gives the momentum loss due to the boundary layer.
(29.8)
With substitution of the similarity solution expressions one obtains a numerical value for δ for a flat plate:
(29.9)
The relations between δ, δ* and δ are illustrated in Fig. 29.1.
Momentum-Integral Equations for the Boundary Layer
- For engineering applications we often prefer approximate integral methods which give sufficiently accurate answers without solving Prandtl's boundary-layer differential equations exactly.
- Kármán and Pohlhausen developed a convenient method that satisfies the boundary conditions and integral conservation laws rather than the full differential equations. The method leads to the momentum integral equation.
- Consider steady, two-dimensional, incompressible flow and integrate the dimensional momentum equation in the streamwise direction from the wall (y = 0) to y = δ (the free-stream interface of the boundary layer). The integrated form is:
(29.10)
- The second term on the left of (29.10) may be expanded as follows.
(29.11)
- Substituting Eq. (29.11) into Eq. (29.10) gives the next integral relation.
(29.12)
- Substituting the relation between the local outer (inviscid) velocity and the free stream velocity U∞ in Eq. (29.12) yields
which can be reduced to the standard momentum integral form shown next.
- Since the integrands vanish outside the boundary layer, the upper integration limit may be extended to infinity (that is, δ → ∞) when using similarity forms.
(29.13)
Substituting the definitions of δ* and δ [Eqs. (29.6) and (29.7)] into Eq. (29.13) leads to:
(29.14)
Equation (29.14) is known as the
momentum integral equation for a two-dimensional incompressible laminar boundary layer. The same integral balance remains valid for turbulent boundary layers also, although the closure (i.e., the relation between δ*, δ
and δ or velocity profile) differs. The wall shear stress tw will, however, be different for laminar and turbulent flows.
- The term representing dU∞/dx (space-wise acceleration of the free stream) signifies the presence of a pressure gradient in the flow direction. A non-zero value indicates that the outer inviscid flow is accelerating or decelerating along x.
- For external flows the sign and magnitude of dU∞/dx depend on the body shape. In entrance regions of internal flows (pipes, channels) a finite dU∞/dx can exist.
- For the classical flat plate in a uniform stream U∞ is constant so dU∞/dx = 0 and the momentum integral equation reduces to the simpler form:
(29.15)
Separation of the Boundary Layer
- Under certain conditions the flow adjacent to the wall reverses direction; this phenomenon is called boundary-layer separation.
- Separation occurs because the boundary layer loses too much momentum near the wall while attempting to move downstream against an increasing pressure (that is, when the pressure gradient in the flow direction is adverse, ∂P/∂x > 0).
- Figure 29.2 illustrates flow past a circular cylinder in an unbounded medium and the development of the wake behind the cylinder.
Consider flow around a cylinder: up to θ = 90° the streamlines behave like those in a nozzle (favourable pressure gradient); beyond θ = 90° the flow area diverges and the outer flow experiences an adverse pressure gradient.
- Up to θ = 90°, the pressure gradient is negative (favourable) and the pressure force and the inertial (streamwise acceleration) force act in the same direction.
- Beyond θ = 90°, the pressure gradient becomes positive (adverse) and the pressure force opposes the inertial force in the inviscid zone.
- In the viscous region close to the wall, viscous forces oppose the combined pressure and inertial forces. For θ ≤ 90° the viscous forces are overcome by conversion of pressure energy into kinetic energy and the boundary layer remains attached.
- For θ > 90° the adverse pressure increases; viscous and pressure forces together may overcome the forward inertial force and the near-wall fluid particles may decelerate, stop and eventually reverse direction. Where reversal first appears is the point of separation and a broad, often unsteady wake develops downstream.
Fig. 29.2 Flow separation and formation of wake behind a circular cylinder
- Until θ = 90° the pressure force and the streamwise acceleration act together (pressure gradient favourable).
- Beyond θ = 90° the pressure gradient is adverse; pressure and inertial forces oppose each other in the inviscid region.
When viscous effects are included near the wall:
- Up to θ = 90° the viscous force opposes the combined pressure and inertial forces, but the fluid overcomes viscous resistance due to the continuous conversion of pressure head into kinetic energy.
- Beyond θ = 90° the viscous zone experiences both adverse pressure and viscous opposition, which may lead to deceleration and eventual flow reversal near the wall.
- Depending on the magnitude of the adverse pressure gradient, separation typically occurs around θ ≈ 90° (in experiments for a cylinder it is observed near θ ≈ 81°). The separated near-wall layers are carried into the wake where vortices and eddies form.
- The point of separation can be identified mathematically as the streamwise location where the wall shear stress becomes zero, i.e. tw = 0, and beyond which reverse flow occurs.
(29.16)
At the separation point the shear stress at the wall tw = 0. Downstream of this point the adverse pressure causes backflow in the near-wall region.
- Separation may also be understood by examining the second derivative of the streamwise velocity u with respect to y at the wall. From the dimensional momentum equation at the wall (where u = v = 0) one obtains the expression below.
(29.17)
- Under a favourable pressure gradient (dP/dx < 0)="" we="" have="" the="" following="" consequences="" for="" the="" curvature="" of="" the="" velocity="">
- (From Eq. (29.17)).
- As y increases toward the free stream the velocity u approaches U∞ asymptotically, so ∂u/∂y decreases with y.
- Thus ∂²u/∂y² remains negative near the edge of the boundary layer and the curvature of the velocity profile is negative (concave downward) as shown in Fig. 29.3a.
- Under an adverse pressure gradient the curvature near the wall must become positive (since ∂P/∂x > 0 implies positive ∂²u/∂y² at the wall), while near the edge of the boundary layer the curvature is still negative. Therefore a point of inflection exists somewhere within the boundary layer where ∂²u/∂y² = 0. This point of inflection is a precursor to separation, but separation itself occurs when the shear at the wall vanishes, i.e. ∂u/∂y|wall = 0.

- At the wall, the curvature ∂²u/∂y² must be positive under adverse pressure gradient (see ).
- Near the boundary layer edge the curvature returns to negative values (see and ), so there must be a point with ∂²u/∂y² = 0 - the point of inflection (see = 0).
- Separation requires ∂u/∂y|wall = 0 (see = 0 at the wall).
- The wall curvature sign is consistent with separation occurring only under adverse pressure gradient; at the edge of the boundary layer we have = 0. It therefore follows that if separation occurs there must exist an inflection point in the velocity profile.
(a) Favourable pressure gradient,
(b) Adverse pressure gradient,
- Returning to the cylinder wake: the inviscid potential flow pressure distribution (shown by the solid line in earlier figures) predicts symmetric pressures. However, experimentally the boundary layer detaches near θ ≈ 81° (rather than 90°), and the separated wake maintains a pressure close to the separation-point pressure.
- In the wake the eddies and rotational flow prevent efficient conversion of rotational kinetic energy into pressure head. The actual pressure distribution thus departs from the inviscid prediction (indicated by a dotted line in figure). The wake pressure is lower than the forward stagnation pressure and this pressure difference produces a net drag on the body.
- The drag arising from pressure difference between forward and wake regions is called form drag. The drag due to shear stress acting on the body surface is called skin friction drag. The total aerodynamic or hydrodynamic drag is generally the sum of form drag and skin friction drag.