Stokes boundary layer

Stokes boundary layer in a viscous fluid due to the harmonic oscillation of a plane rigid plate. Velocity (blue line) and particle excursion (red dots) as a function of the distance to the wall.

In fluid dynamics, the Stokes boundary layer, or oscillatory boundary layer, refers to the boundary layer close to a solid wall in oscillatory flow of a viscous fluid. Or, it refers to the similar case of an oscillating plate in a viscous fluid at rest, with the oscillation direction(s) parallel to the plate. For the case of laminar flow at low Reynolds numbers over a smooth solid wall, George Gabriel Stokes – after whom this boundary layer is called – derived an analytic solution, one of the few exact solutions for the Navier–Stokes equations.[1] [2] In turbulent flow, this is still named a Stokes boundary layer, but now one has to rely on experiments, numerical simulations or approximate methods in order to obtain useful information on the flow.

The thickness of the oscillatory boundary layer is called the Stokes boundary-layer thickness.

An important observation from Stokes' solution for the oscillating Stokes flow is, that vorticity oscillations are confined to a thin boundary layer and damp exponentially when moving away from the wall.[3] This observation is also valid for the case of a turbulent boundary layer. Outside the Stokes boundary layer – which is often the bulk of the fluid volume – the vorticity oscillations may be neglected. To good approximation, the flow velocity oscillations are irrotational outside the boundary layer, and potential flow theory can be applied to the oscillatory part of the motion. This significantly simplifies the solution of these flow problems, and is often applied in the irrotational flow regions of sound waves and water waves.


Stokes boundary layer for laminar flow near a wall

The oscillating flow is assumed to be uni-directional and parallel to the plane wall. The only non-zero velocity component is called u (SI measure in meter/second, or m/s) and is in the x-direction parallel to the oscillation direction. Moreover, since the flow is taken to be incompressible, the velocity component u is only a function of time t (in seconds) and distance from the wall z (in meter). The Reynolds number is taken small enough for the flow to be laminar. Then the Navier–Stokes equations, without additional forcing, reduce to:[4]

\frac{\partial u}{\partial t} = -\frac{1}{\rho} \frac{\partial p}{\partial x} + \nu \frac{\partial^2 u}{\partial z^2},


  • ρ the mass density of the fluid (kg/m3), taken to be a constant,
  • p the fluid pressure (SI: Pa),
  • ν the kinematic viscosity of the fluid (m2/s), also taken constant.


  • u the velocity of the fluid along the plate (m/s)
  • x the position along the plate (m)
  • z the distance from the plate (m)
  • t the time (s)

Because the velocity u is not a function of position x along the plate, the pressure gradient ∂p/∂x is also independent of x (but the pressure p varies linearly with x). Moreover, the Navier–Stokes equation for the velocity component perpendicular to the wall reduces to ∂p/∂z = 0, so the pressure p and pressure gradient ∂p/∂x are also independent of the distance to the wall z. In conclusion, the pressure forcing ∂p/∂x can only be a function of time t.[4]

The only non-zero component of the vorticity vector is the one in the direction perpendicular to x and z, called ω (in s-1) and equal to:[3]

\omega = \frac{\partial u}{\partial z}.

Taking the z-derivative of the above equation, ω has to satisfy[3]

\frac{\partial \omega}{\partial t} = \nu \frac{\partial^2 \omega}{\partial z^2}.

As usual for the vorticity dynamics, the pressure drops out of the vorticity equation.[5]

Oscillation of a plane rigid plate

Harmonic motion, parallel to a plane rigid plate, will result in the fluid near the plate being dragged with the plate, due to the viscous shear stresses. Suppose the motion of the plate is

u_0(t) = U_0\, \cos\left( \Omega\, t \right),\,


The plate, located at z = 0, forces the viscous fluid adjacent to have the same velocity u1zt ) resulting in the no-slip condition:

u_1(0,t) = u_0(t) = U_0\, \cos\left( \Omega\, t \right) \quad \text{ at }\; z = 0.

Far away from the plate, for z → ∞, the velocity u1 approaches zero. Consequently, the pressure gradient ∂p/∂x is zero at infinity and, since it is only a function of time t and not of z, has to be zero everywhere:[6]

\frac{\partial u_1}{\partial t} = \nu \frac{\partial^2 u_1}{\partial z^2}.

Such an equation is called a one-dimensional heat equation or diffusion equation.

As a result, the solution for the flow velocity is[7]

  u_1(z,t) = U_0\, \text{e}^{-\kappa\, z}\, \cos\left( \Omega\, t\, -\, \kappa\, z\right)
  \quad \text{ with }\; \kappa\, =\, \sqrt{\frac{\Omega}{2\nu}}.

Here, κ is a kind of wavenumber in the z-direction, associated with a length[7]

\delta = \frac{2\pi}{\kappa} = 2\pi\, \sqrt{\frac{2\nu}{\Omega}}

which is called the Stokes boundary-layer thickness. At a distance δ from the plate, the velocity amplitude has been reduced to e-2π ≈ 0.002 times its value U0 at the plate surface. Further, as can be seen from the phase changes Ω t - κ z in the solution u1, the velocity oscillations propagate as a damped wave away from the wall, with wavelength δ and phase speed Ω / κ.

The vorticity ω1 is equal to

    = \frac{\partial u_1}{\partial z} 
    = -\kappa\, U_0\, \text{e}^{-\kappa\, z}\, 
         \cos\left( \Omega\, t\, -\, \kappa\, z \right)\, 
         \sin\left( \Omega\, t\, -\, \kappa\, z \right)\,

and, as u1, dampens exponentially in amplitude when moving away from the plate surface.

Flow due to an oscillating pressure gradient near a plane rigid plate

The case for an oscillating far-field flow, with the plate held at rest, can easily be constructed from the previous solution for an oscillating plate by using linear superposition of solutions. Consider a uniform velocity oscillation u:

u_\infty(z,t) = U_0\, \cos\left( \Omega\, t \right), \,

which satisfies the flow equations for the Stokes boundary layer, provided it is driven by a pressure gradient

\frac{\partial p_2}{\partial x} = \rho\, \Omega\, U_0\, \sin\left( \Omega\, t \right).

Subtracting the solution u1zt ) from uzt ) gives the desired solution for an oscillating flow near a rigid wall at rest:[3]

  u_2(z,t) = U_0\, \Bigl[\, \cos\left(  \Omega\, t \right)\, -\, \text{e}^{-\kappa\, z}\, \cos\left( \Omega\, t\, -\, \kappa\, z \right)\, \Bigr],

which is zero at the wall z = 0, corresponding with the no-slip condition for a wall at rest. Further the velocity u2 oscillates with amplitude U0 far away from the wall, z → ∞. This situation is often encountered in sound waves near a solid wall, or for the fluid motion near the sea bed in water waves.

The vorticity, for the oscillating flow near a wall at rest, is equal to the vorticity in case of an oscillating plate but of opposite sign: ω2 = - ω1.

See also


  1. ^ Wang, C. Y. (1991). "Exact solutions of the steady-state Navier-Stokes equations". Annual Review of Fluid Mechanics 23: 159–177. Bibcode 1991AnRFM..23..159W. doi:10.1146/annurev.fl.23.010191.001111. 
  2. ^ Landau & Lifshitz (1987), pp. 83–85.
  3. ^ a b c d Phillips (1977), p. 46.
  4. ^ a b Batchelor (1967), p. 179.
  5. ^ Since the vorticity equation is obtained by taking the curl of the Navier–Stokes equations, and the curl of the pressure gradient equals zero, see vector calculus identities.
  6. ^ Batchelor (1967), p. 190.
  7. ^ a b c Batchelor (1967), p. 192.


  • Batchelor, G.K. (1967). An Introduction to Fluid Dynamics. Cambridge University Press. ISBN 0521663962. 
  • Lamb, H. (1994). Hydrodynamics (6th ed.). Cambridge University Press. ISBN 9780521458689.  Originally published in 1879, the 6th extended edition appeared first in 1932.
  • Landau, L.D.; Lifshitz, E.M. (1987). Fluid Mechanics. Course of theoretical physics. 6 (2nd ed.). Pergamon Press. ISBN 0 08 339932 8 
  • Phillips, O.M. (1977). The dynamics of the upper ocean (2nd ed.). Cambridge University Press. ISBN 0 521 29801 6. 

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Boundary layer — For the anatomical structure, see Boundary layer of uterus. Boundary layer visualization, showing transition from laminar to turbulent condition In physics and fluid mechanics, a boundary layer is that layer of fluid in the immediate vicinity of… …   Wikipedia

  • Planetary boundary layer — Depiction of where the planetary boundary layer lies on a sunny day The planetary boundary layer (PBL), also known as the atmospheric boundary layer (ABL), is the lowest part of the atmosphere and its behavior is directly influenced by its… …   Wikipedia

  • Stokes drift — For a pure wave motion in fluid dynamics, the Stokes drift velocity is the average velocity when following a specific fluid parcel as it travels with the fluid flow. For instance, a particle floating at the free surface of water waves,… …   Wikipedia

  • Navier–Stokes equations — Continuum mechanics …   Wikipedia

  • George Gabriel Stokes — Infobox Scientist box width = 300px name = George Stokes image size = 220px caption = Sir George Gabriel Stokes, 1st Baronet (1819–1903) birth date = birth date|1819|8|13|df=y birth place = Skreen, County Sligo, Ireland death date = death date… …   Wikipedia

  • Divergent boundary — In plate tectonics, a divergent boundary or divergent plate boundary (also known as a constructive boundary or an extensional boundary) is a linear feature that exists between two tectonic plates that are moving away from each other. Divergent… …   Wikipedia

  • Convergent boundary — In plate tectonics, a convergent boundary, also known as a destructive plate boundary (because of subduction), is an actively deforming region where two (or more) tectonic plates or fragments of lithosphere move toward one another and collide. As …   Wikipedia

  • Airy wave theory — In fluid dynamics, Airy wave theory (often referred to as linear wave theory) gives a linearised description of the propagation of gravity waves on the surface of a homogeneous fluid layer. The theory assumes that the fluid layer has a uniform… …   Wikipedia

  • fluid mechanics — an applied science dealing with the basic principles of gaseous and liquid matter. Cf. fluid dynamics. [1940 45] * * * Study of the effects of forces and energy on liquids and gases. One branch of the field, hydrostatics, deals with fluids at… …   Universalium

  • D'Alembert's paradox — In fluid dynamics, d Alembert s paradox (or the hydrodynamic paradox) is a contradiction reached in 1752 by French mathematician Jean le Rond d Alembert.[1] D Alembert proved that – for incompressible and inviscid potential flow – the drag force… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.