Rayleigh problem explained

In fluid dynamics, Rayleigh problem also known as Stokes first problem is a problem of determining the flow created by a sudden movement of an infinitely long plate from rest, named after Lord Rayleigh and Sir George Stokes. This is considered as one of the simplest unsteady problems that have an exact solution for the Navier-Stokes equations. The impulse movement of semi-infinite plate was studied by Keith Stewartson.[1]

Flow description

Consider an infinitely long plate which is suddenly made to move with constant velocity

U

in the

x

direction, which is located at

y=0

in an infinite domain of fluid, which is at rest initially everywhere. The incompressible Navier-Stokes equations reduce to[2] [3]
\partialu
\partialt

=\nu

\partial2u
\partialy2

where

\nu

is the kinematic viscosity. The initial and the no-slip condition on the wall are

u(y,0)=0,u(0,t>0)=U,u(infty,t>0)=0,

the last condition is due to the fact that the motion at

y=0

is not felt at infinity. The flow is only due to the motion of the plate, there is no imposed pressure gradient.

Self-Similar solution

The problem on the whole is similar to the one dimensional heat conduction problem. Hence a self-similar variable can be introduced[4]

η=

y
\sqrt{\nut
}, \quad f(\eta) = \frac

Substituting this the partial differential equation, reduces it to ordinary differential equation

f''+

1
2

ηf'=0

with boundary conditions

f(0)=1,f(infty)=0

The solution to the above problem can be written in terms of complementary error function

u=Uerfc\left(

y
\sqrt{4\nut
}\right)

The force per unit area exerted on the plate is

F=\mu\left(

\partialu
\partialy

\right)y=0=-\rho\sqrt{

\nuU2
\pit
}

Arbitrary wall motion

Instead of using a step boundary condition for the wall movement, the velocity of the wall can be prescribed as an arbitrary function of time, i.e.,

U=f(t)

. Then the solution is given by[5]

u(y,t)=

t
\int
0
f(\tau)
2\sqrt{\pi\nu
}\frace^d\tau.

Rayleigh's problem in cylindrical geometry

Rotating cylinder

Consider an infinitely long cylinder of radius

a

starts rotating suddenly at time

t=0

with an angular velocity

\Omega

. Then the velocity in the

\theta

direction is given by

v\theta=

a\Omega
2\pii
iinfty
\int
-iinfty
K1(r\sqrt{s/\nu
)}{K
st
1(a\sqrt{s/\nu})}e
ds
s

where

K1

is the modified Bessel function of the second kind. As

t → infty

, the solution approaches that of a rigid vortex. The force per unit area exerted on the cylinder is

F=\mu\left(

\partialv\theta-
\partialr
v\theta
r

\right)r=a=

\rhoa2\Omega
t
-a2
2\nut
e
I
0\left(a2
2\nut

\right)-2\mu\Omega

where

I0

is the modified Bessel function of the first kind.

Sliding cylinder

Exact solution is also available when the cylinder starts to slide in the axial direction with constant velocity

U

. If we consider the cylinder axis to be in

x

direction, then the solution is given by

u=

U
2\pii
iinfty
\int
-iinfty
K0(r\sqrt{s/\nu
)}{K
st
0(a\sqrt{s/\nu})}e
ds
s

.

See also

References

  1. Stewartson, K. T. (1951). On the impulsive motion of a flat plate in a viscous fluid. The Quarterly Journal of Mechanics and Applied Mathematics, 4(2), 182-198.
  2. Batchelor, George Keith. An introduction to fluid dynamics. Cambridge university press, 2000.
  3. Lagerstrom, Paco Axel. Laminar flow theory. Princeton University Press, 1996.
  4. Acheson, David J. (1990) Elementary Fluid Dynamics, Oxford University Press
  5. Dryden, Hugh L., Francis D. Murnaghan, and Harry Bateman. Hydrodynamics. New York: Dover publications, 1956.