Ladd 03
Ladd 03
Ladd 03
MATERIAL
w dl
L
Kelvins Theorem:
Assume that:
the volume force field f potential,
the fluid is inviscid and barotropic
the flow is stationary.
Then: the circulation of the velocity field along any closed material line L(t) is constant
in time, i.e.
d
dt
(t ) dtd
(t , x) d l 0
L (t )
P 1 p
Thus, the acceleration (which consists of the convective part only) can be expressed as
a ( ) ( P )
In order to evaluate the time derivative of the circulation along the material line, it is
convenient to use Lagrange approach. Thus, the circulation can be expresses as
(t )
(t , x) d l V (t , ) J (t , )d l0
L (t )
L (t )
0
where
Lagrangian coordinates.
d
d
(
t
,
x
)
d
l
V (t , ) J (t , )d l0 a (t , ) J (t , )d l0
dt L (t )
dt L (t )
L (t )
0
L (t )
V (t , ) V (t , ) d l0
L (t )
a (t , x ) d l
L (t )
( 21 V V )(t , ) d l0
0
int. of the grad . along the closed loop
( P ) d l 0
L (t )
0
int. of the grad . along
the closed loop
VO
OR
RT
TIIC
CIIT
TY
Y
As we already know, the vorticity is defined as the rotation of the velocity:
rot
Definitions:
A vortex line is the line of the vorticity vector field. At each point of such line, the
vorticity vector is tangent to this line.
The vortex tube is the subset of the flow domain bounded by the surface made of the
vortex lines passing through all point of a given closed contour (the contour L on the
picture below)
vorticity field
vortex line
STTRREENNG
GT
TH
HO
OF
FT
TH
HE
EV
VO
OR
RT
TE
EX
XT
TU
UB
BE
E
S2
S n d l dx
Sside
n1
S1
0 dx n ds n ds
S1
S2
n ds 0
side
Note that the last integral vanishes because the surface Sside is made of the vortex lines and
thus at each point of Sside the normal versor n is perpendicular to the vorticity vector.
Note also that the orientations of the normal versors at S1 and S2 are opposite (in order to
apply the GGO Theorem, the normal versor must point outwards at all components of the
boundary ).
Reversing the orientation of n at S2, we conclude that
S n ds S n ds
1
R
RD
D
HEELLM
MH
HO
OL
LT
TZ
Z (3 ) TH
HE
EO
OR
RE
EM
M
Assume that:
the flow is inviscid and barotropic,
the volume force field is potential.
Then: the vortex lines consist of the same fluid elements, i.e. the lines of the vorticity field
are material lines.
Proof:
We need the transformation rule for the vectors tangent to a material line.
Let at initial time t = 0 the material line be described parametrically as l0 :a a ( s) .
At some later time instant t > 0, the shape of the material line follows from the flow mapping
Ft : R 3 a
x R 3 , i.e., l : x x( s) Ft [a( s)].
The corresponding transformation of the tangent vector can be evaluated as follows
(s) dds x(s) dds Ft (a( s)) [ x ](a( s)) dds a( s) [ x ](a( s)) 0 ( s)
Jacobi matrix
In the above, the following vector identity, written for p and q , is used
( p q) (q ) p ( p ) q ( q) p ( p) q
Next, one can calculate the Lagrangian derivative of the vector field
1 DtD 1 DtD
D 1
Dt
/ as follows
1 a ( ) ( )
1 1 a ( 1 )
From the equation of motion and assumed flow properties that the acceleration field is
potential and thus
Then, the equation for the vector field
a 0
/ reduces to
( 1 )
D 1
Dt
i xij c j
, or equivalently,
[ x ] c .
Jacobi
matrix
dt
R ( 1 ) [ x ] c [ x ][ x ] c
c 0
c(t , ) c(0, ) c0
Note that for the initial time t = 0 the transformation between Lagrange and Euler
descriptions reduces to identity.
[ x ]
Therefore
c0 10 0
and since
t 0
c(t ) c0 we get
1 [ x ] 1
0 0 .
The last equality has the form of the transformation rule for the vectors tangent to
material lines. Since the vector 0 / 0 is tangent to the vortex line passing through the
point at t = 0, it follows that the vector / is tangent to image of this line at some later
time t. But / is also tangent to the vortex line passing through the point x, which means
that the vortex lines must be material.
Since the vortex lines are material, so are the vortex tubes. If we define a closed, material
contour lying on the vortex tubes surface (and wrapped around it), then such a contour
remains on this surface for any time. It follows from the Kelvin Theorem that the circulation
along such contour remains constant. Consequently, the strength of any vortex tube also
remains constant in time. It is important conclusion showing that the vortex motion of the
inviscid, barotropic fluid exposed to a potential force field cannot be created or destroyed.
EQ
QU
UA
AT
TIIO
ON
NO
OF
FT
TH
HE
EV
VO
OR
RT
TIIC
CIIT
TY
YT
TR
RA
AN
NSSP
PO
OR
RT
T
In fluid mechanics the vorticity plays a very important role, in particular in understanding of
the phenomenon of turbulence. In this section we derive the differential equation governing
spatial/temporal evolution of this field.
Recall that the equation of motion of an inviscid fluid can be written in the following form
t ( 21 2 ) 1 p f
Thus, the application of the rotation operator yields
t ( ) ( 1 p) f
The pressure term can be transformed as follows
( 1 p) ( 1 ) p 1 p 12 p
Note: the above term vanishes identically when the fluid is barotropic since the gradients of
pressure and density are in such case parallel.
t ( ) ( ) 12 p f
or, using the full derivative
D
Dt
( )
vortex stretching
term
p f
baroclinic
term
nonpotential
volume force
term
For the barotropic (in particular incompressible) motion of inviscid fluid, the vorticity
equation reduces to
D ( )
Dt
D 0
Dt
We conclude that in any 2D flow the vorticity is conserved along trajectories of fluid
elements.
If the fluid is viscous, the vorticity equation contains the diffusion term. We will derive this
equation assuming that the fluid is incompressible. Again, we begin with the Navier-Stokes
equation in the Lamb-Gromeko form
t ( 21 2 ) 1 p f
If the rotation operator is applied, we get the equation
t ( ) ( ) f
which reduces to
t ( ) ( )
rot v rot ( grad divv rot rotv ) rot rot grad div rot rot
showing that the vector Laplace and rotation operators commute.
The vorticity equation can be also written equivalently as
D
Dt ( )
The viscous term describes the diffusion of vorticity due to fluid viscosity. This effect
smears the vorticity over the whole flow domain. Thus, in the viscous case the vortex lines are
not material lines anymore.
There exists a relation between the streamfunction and vorticity. Since the flow is 2D, the
vorticity field is perpendicular to the flows plane and can be expressed as
(12 21 ) e3 e3
Then, the streamfunction satisfies the Poisson equation
11 22 (12 21 )
Two dimensional motion of an incompressible viscous fluid can be described in terms of the
purely kinematical quantities: velocity, vorticity and streamfunction. The pressure field is
eliminated and the continuity equation div 0 is automatically satisfied. The complete
description consists of the following equations:
t 1 1 2 2
1 2 , 2 1
12 21