Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Ladd 03

Download as pdf or txt
Download as pdf or txt
You are on page 1of 17

KINEMATICS OF FLUID ADDITIONAL

MATERIAL

CIRCULATION, VORTICITY AND STREAMFUNCTIONS


CIIRRCCUULLAATTIIO
ON
N
Definition: Circulation of the vector field w along the (closed) contour L is defined as

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 )

Proof of the Kelvin Theorem:


Since the flow is barotropic and the volume force field is potential, we can write

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

J (t , ) x denotes the Jacobi matrix of the transformation between Eulerian and

Lagrangian coordinates.

Then, the time derivative of the circulation is evaluated as follows

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

where the relation t J (t , ) V (t , ) has been used.

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

It is defined as the flux of vorticity through a cross-section of the


n2 tube. Using the Stokes Theorem we can write:

S n d l dx
Sside
n1

We see that the strength of the vortex tube is equal to the


circulation of the velocity along a closed contour wrapped
around the tube.

S1

The above definition does not depend on the choice of a


particular contour. Indeed, since the vorticity field is divergence-free, the flux of the
vorticity is fixed along the vortex tube. To see this, consider the tube segment located
between two cross-section S1 and S2 .
From the GGO theorem we have

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

Lets now write the acceleration in the Lamb-Gromeko form:


D ( 1 2 )
a Dt
t
2

The rotation of a can be expressed as


D ( ) ( )
a t ( ) ( ) Dt

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

Define the vector field c such that

i xij c j

, or equivalently,

[ x ] c .
Jacobi
matrix

In the above, the symbol denotes the Lagrange variables.


The left-hand side of the above equation can be transformed as follows
D 1 d [ x ]c [ x ] d c [ ]c [ x ] d c [ ][ x ]c
L Dt
dt
x
dt

dt

The right-hand side can be written as

R ( 1 ) [ x ] c [ x ][ x ] c

Since L = R, we conclude that


d
dt

c 0

Thus, c is constant along trajectories of the fluid elements.

c(t , ) c(0, ) c0

Using the Lagrange description, we can write

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.

The equation of the vorticity transport can be written in the form

t ( ) ( ) 12 p f
or, using the full derivative

D
Dt

( )
vortex stretching
term

p f
baroclinic
term

nonpotential
volume force
term

The change of the vorticity appears due to the following factors:


Local deformation of the pattern of vortex lines (or vortex tubes) known as the vortex
stretching effect. This mechanism is believed to be crucial for generating spatial/temporal
complexity of turbulent flows. The vortex stretching term vanishes identically for 2D
flows.
Presence of baroclinic effects. If the flow is not barotropic then the gradients of pressure
and density field are nonparallel. It can be shown that in such situation a torque is
developed which perpetuates rotation of fluid elements (generates vorticity).
Presence of nonpotential volume forces. This factor is important e.g. for electricityconducting fluids.

For the barotropic (in particular incompressible) motion of inviscid fluid, the vorticity
equation reduces to
D ( )
Dt

In the 2D case it reduces further to

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 ( ) ( )

when the field of the volume forces f is potential.

In the above, the following operator identity has been used

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:

Equation of the vorticity transport (2D)


Equation for the streamfunction
Relation between the streamfunction and velocity
Definition of vorticity (2D)

t 1 1 2 2

1 2 , 2 1
12 21

accompanied by appropriately formulated boundary and initial conditions.

You might also like