WOLFRAM|DEMONSTRATIONS PROJECT

Steady Flow over a Rotating Disk: von Kármán Swirling Flow

​
velocity profile
3D trajectory
kinematic viscosity
0.44
angular velocity
0.85
The steady flow induced by a infinite disk that rotates in its own plane at
z=0
is a classical problem in fluid mechanics. It is one of the few examples of a viscous flow that involves all three components of velocity and admits an exact solution to the Navier–Stokes equations. The velocity field for the swirling flow is given by:
v(r,z)=u(r,z)
e
r
+v(r,z)
e
θ
+w(r,z)
e
z
. The Navier–Stokes equations reduce to:
1
r
∂
∂r
(ru)+
∂w
∂z
=0
,
u
∂u
∂r
+w
∂u
∂z
-
2
v
r
=
-1
ρ
∂p
∂r
+ν
2
∂
u
∂
2
r
+
1
r
∂u
∂r
+
2
∂
u
∂
2
z
-
u
2
r
,
u
∂v
∂r
+w
∂v
∂z
+
vu
r
=ν
2
∂
v
∂
2
r
+
1
r
∂v
∂r
+
2
∂
v
∂
2
z
-
v
2
r
,
u
∂w
∂r
+w
∂w
∂z
=
-1
ρ
∂p
∂z
+ν
2
∂
w
∂
2
r
+
1
r
∂w
∂r
+
2
∂
w
∂
2
z
.
The boundary conditions at the disk surface are no slip and impenetrability conditions:
u(r,0)=0,v(r,0)=Ω,w(r,0)=0
,
where
Ω
is the rotational speed of the disk. Far away from the surface of the disk, viscous effects are negligible, such that
u(r,z)0,v(r,z)0asz∞
.
However, the viscous pumping action of the disk is balanced by a uniform axial inflow at infinity:
w(r,z)constantasz∞
.
von Kármán was able to show that the above equations admit a self-similar solution defined by:
u(r,z)=ΩrF(η)
,
v(r,z)=ΩrG(η)
,
w(r,z)=
Ων
H(η)
, and
p(r,z)=ρΩP(η)
, where
η=z
Ω
ν
, and the functions
F
,
G
,
H
, and
P
are determined by:
H'+2F=0
,
F''=-
2
G
+
2
F
+F'H
,
G''=2FG+HG'
,
and
P'=2FH-2F'
,
subject to the boundary conditions:
F(0)=H(0)=P(0)=0,G(0)=1
,
and
F(∞)=G(∞)=0
.
It is evident from the structure of the equations that the equation for the pressure field
P
is decoupled from the equations that define the velocity field given by the functions
F
,
G
, and
H
. The boundary value problem for the velocity field can be solved by a shooting method using Mathematica. Also of interest are the fluid particle trajectories in the swirling flow; they can be computed from the velocity field by solving:
dr
dt
=ΩrF(η)
,
dz
dt
=
Ων
H(η)
,
r
dθ
dt
=ΩrG(η)
.
It is noteworthy that the boundary layer thickness in the above swirling flow is constant over the disk surface; this feature of the flow is exploited in electrochemistry to study mass transfer.
This Demonstration shows the velocity components of the swirling flow as a function of the angular velocity of the disk and the kinematic viscosity of the fluid. The thickness of the boundary layer is also displayed. Fluid particle trajectories are also shown for different starting positions in the flow.