 Stokes boundary layer

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 boundarylayer 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.
Contents
Stokes boundary layer for laminar flow near a wall
The oscillating flow is assumed to be unidirectional and parallel to the plane wall. The only nonzero velocity component is called u (SI measure in meter/second, or m/s) and is in the xdirection 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]}
with:
 ρ the mass density of the fluid (kg/m^{3}), taken to be a constant,
 p the fluid pressure (SI: Pa),
 ν the kinematic viscosity of the fluid (m^{2}/s), also taken constant.
and
 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 nonzero component of the vorticity vector is the one in the direction perpendicular to x and z, called ω (in s^{1}) and equal to:^{[3]}
Taking the zderivative of the above equation, ω has to satisfy^{[3]}
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
with
 U_{0} the velocity amplitude of the plate motion (in m/s), and
 Ω the angular frequency of the motion (in rad/s).
The plate, located at z = 0, forces the viscous fluid adjacent to have the same velocity u_{1}( z, t ) resulting in the noslip condition:
Far away from the plate, for z → ∞, the velocity u_{1} 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]}
Such an equation is called a onedimensional heat equation or diffusion equation.
Details on the derivation of the solution This equation can be solved for harmonic motion using complex numbers and separation of variables:^{[7]} with i^{2} = 1 and ℜ{•} denoting the real part of the quantity between brackets. Then from the flow equation, F(z) is required to satisfy:
with the boundary conditions:
The solution for F(z) becomes:
As a result, the solution for the flow velocity is^{[7]}
Here, κ is a kind of wavenumber in the zdirection, associated with a length^{[7]}
which is called the Stokes boundarylayer thickness. At a distance δ from the plate, the velocity amplitude has been reduced to e^{2π} ≈ 0.002 times its value U_{0} at the plate surface. Further, as can be seen from the phase changes Ω t  κ z in the solution u_{1}, the velocity oscillations propagate as a damped wave away from the wall, with wavelength δ and phase speed Ω / κ.
The vorticity ω_{1} is equal to
and, as u_{1}, 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 farfield 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_{∞}:
which satisfies the flow equations for the Stokes boundary layer, provided it is driven by a pressure gradient
Subtracting the solution u_{1}( z, t ) from u_{∞}( z, t ) gives the desired solution for an oscillating flow near a rigid wall at rest:^{[3]}
which is zero at the wall z = 0, corresponding with the noslip condition for a wall at rest. Further the velocity u_{2} oscillates with amplitude U_{0} 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
Notes
 ^ Wang, C. Y. (1991). "Exact solutions of the steadystate NavierStokes equations". Annual Review of Fluid Mechanics 23: 159–177. Bibcode 1991AnRFM..23..159W. doi:10.1146/annurev.fl.23.010191.001111.
 ^ Landau & Lifshitz (1987), pp. 83–85.
 ^ ^{a} ^{b} ^{c} ^{d} Phillips (1977), p. 46.
 ^ ^{a} ^{b} Batchelor (1967), p. 179.
 ^ 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.
 ^ Batchelor (1967), p. 190.
 ^ ^{a} ^{b} ^{c} Batchelor (1967), p. 192.
References
 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.
Categories: Fluid dynamics
 Boundary layers
Wikimedia Foundation. 2010.