icmm tgzielinski idealflow.slides
TRANSCRIPT
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
1/85
Introduction Ideal ow theory Basic aerodynamics
Fundamentals of Fluid Dynamics:Ideal Flow Theory
Introductory Course on Multiphysics Modelling
TOMASZ G. Z IELINSKI(after: D.J. A CHESON s Elementary Fluid Dynamics )
Institute of Fundamental Technological Research
Warsaw
Poland
http://www.ippt.gov.pl/~tzielins/
mailto:[email protected]:[email protected]:[email protected]:[email protected]:[email protected]:[email protected]://www.ippt.gov.pl/~tzielins/http://www.ippt.gov.pl/~tzielins/mailto:[email protected] -
7/29/2019 ICMM TGZielinski IdealFlow.slides
2/85
Introduction Ideal ow theory Basic aerodynamics
Outline
1 IntroductionMathematical preliminariesPreliminary ideas and denitionsConvective derivative
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
3/85
Introduction Ideal ow theory Basic aerodynamics
Outline
1 IntroductionMathematical preliminariesPreliminary ideas and denitionsConvective derivative
2 Ideal ow theory
Ideal uidIncompressibility conditionEulers equations of motionBernoulli theoremsVorticity
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
4/85
Introduction Ideal ow theory Basic aerodynamics
Outline
1 IntroductionMathematical preliminariesPreliminary ideas and denitionsConvective derivative
2 Ideal ow theory
Ideal uidIncompressibility conditionEulers equations of motionBernoulli theoremsVorticity
3 Basic aerodynamicsSteady ow past a xed wingFluid circulation round a wingKuttaJoukowski theorem and conditionConcluding remarks
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
5/85
Introduction Ideal ow theory Basic aerodynamics
Outline
1 IntroductionMathematical preliminariesPreliminary ideas and denitionsConvective derivative
2 Ideal ow theory
Ideal uidIncompressibility conditionEulers equations of motionBernoulli theoremsVorticity
3 Basic aerodynamicsSteady ow past a xed wingFluid circulation round a wingKuttaJoukowski theorem and conditionConcluding remarks
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
6/85
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
7/85
Introduction Ideal ow theory Basic aerodynamics
Mathematical preliminaries
Theorem (The divergence theorem)
Let the region V
be bounded by a simple surface S
with unit outward normal n . Then:
S
f n dS =V
f dV ; in particular S
f n dS =V
f dV .
Theorem (Stokes theorem)
Let C be a simple closed curve spanned by a surface S with unit normal n . n
S
C
Then:
C
f d x =
S
f n dS .
Greens theorem in the plane may be viewed as a special case ofStokes theorem (with f = u ( x, y), v( x, y), 0 ):
C
u d x+ v d y =S
v
x
u
yd xd y .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
8/85
Introduction Ideal ow theory Basic aerodynamics
Preliminary ideas and denitionsA usual way of describing a uid ow is by means of the owvelocity dened at any point x = ( x, y, z) and at any time t :
u = u ( x, t) = u ( x, t), v( x, t), w( x, t) .
Here, u , v, w are the velocity components in Cartesian coordinates.
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
9/85
Introduction Ideal ow theory Basic aerodynamics
Preliminary ideas and denitions
u = u ( x, t) = u ( x, t), v( x, t), w( x, t) .
Denition (Steady ow)
A steady ow is one for which
u
t= 0 , that is, u = u ( x) = u ( x), v( x), w( x) .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
10/85
Introduction Ideal ow theory Basic aerodynamics
Preliminary ideas and denitions
u = u ( x, t) = u ( x, t), v( x, t), w( x, t) .
Denition (Steady ow)
A steady ow is one for which
u
t= 0 , that is, u = u ( x) = u ( x), v( x), w( x) .
Denition (Two-dimensional ow)
A two-dimensional ow is of the form
u = u ( x, t), v( x, t), 0 where x = ( x, y) .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
11/85
y y
Preliminary ideas and denitions
u = u ( x, t) = u ( x, t), v( x, t), w( x, t) .
Denition (Steady ow)
A steady ow is one for which
u
t= 0 , that is, u = u ( x) = u ( x), v( x), w( x) .
Denition (Two-dimensional ow)
A two-dimensional ow is of the form
u = u ( x, t), v( x, t), 0 where x = ( x, y) .
Denition (Two-dimensional steady ow)
A two-dimensional steady ow is of the form
u = u ( x), v( x), 0 where x = ( x, y) .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
12/85
Streamlines and particle paths (pathlines)
Denition (Streamline)
A streamline is a curve which, at any particular time t , has the samedirection as u ( x, t) at each point. A streamline x = x( s), y = y( s), z = z( s)( s is a parameter) is obtained by solving at a particular time t :
d xd s
u
=
d yd s
v
=d zd s
w
.
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
13/85
Streamlines and particle paths (pathlines)
Denition (Streamline)
A streamline is a curve which, at any particular time t , has the samedirection as u ( x, t) at each point. A streamline x = x( s), y = y( s), z = z( s)( s is a parameter) is obtained by solving at a particular time t :
d xd s
u
=
d yd s
v
=d zd s
w
.
For a steady ow the streamline pattern is the same at all times,and uid particles travel along them.In an unsteady ow, streamlines and particle paths are usually
quite different.
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
14/85
Streamlines and particle paths (pathlines)
Denition (Streamline)
A streamline is a curve which, at any particular time t , has the samedirection as u ( x, t) at each point. A streamline x = x( s), y = y( s), z = z( s)( s is a parameter) is obtained by solving at a particular time t :
d xd s
u
=
d yd s
v
=d zd s
w
.
Example
Consider a two-dimensional ow described as follows
u ( x, t) = u 0 , v( x, t) = at , w( x, t) = 0,
where u 0 and a are positive constants. Now, notice that:in this ow streamlines are (always) straight lines,
uid particles follow parabolic paths: x(t) = u0
t , y(t) = a t2
2 y = ax
2
2u 20.
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
15/85
Streamlines and particle paths (pathlines)
u ( x, t) = u 0 , v( x, t) = at , w( x, t) = 0,
x
y
t = 0
in this ow streamlines are (always) straight lines,
uid particles follow parabolic paths: x(t) = u 0 t , y(t) =a t 2
2 y =ax 2
2u 20 .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
16/85
Streamlines and particle paths (pathlines)
u ( x, t) = u 0 , v( x, t) = at , w( x, t) = 0,
x
y
t = t
in this ow streamlines are (always) straight lines,
uid particles follow parabolic paths: x(t) = u 0 t , y(t) =a t 2
2 y =ax 2
2u 20 .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
17/85
Streamlines and particle paths (pathlines)
u ( x, t) = u 0 , v( x, t) = at , w( x, t) = 0,
x
y
t = 2 t
in this ow streamlines are (always) straight lines,
uid particles follow parabolic paths: x(t) = u 0 t , y(t) =a t 2
2 y =ax 2
2u 20 .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
18/85
Streamlines and particle paths (pathlines)
u ( x, t) = u 0 , v( x, t) = at , w( x, t) = 0,
x
y
t = 0
t = t
t = 2 t
in this ow streamlines are (always) straight lines,
uid particles follow parabolic paths: x(t) = u 0 t , y(t) =a t 2
2 y =ax 2
2u 20 .
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
19/85
Convective derivative
Let f ( x, t) denote some quantity of interest in the uid motion.
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
20/85
Convective derivative
Let f ( x, t) denote some quantity of interest in the uid motion.
Note that f
tmeans the rate of change of f at a xed
position x in space.
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
21/85
Convective derivative
Let f ( x, t) denote some quantity of interest in the uid motion.
Note that f
tmeans the rate of change of f at a xed
position x in space.The rate of change of f following the uid is
D f D t
=dd t
f x(t), t = f x
d xd t
+ f t
where x(t) = x(t), y(t), z(t) is understood to change with time atthe local ow velocity u , so that
d xd t
=d xd t
,d yd t
,d zd t
= u , v, w .
Therefore, f x
d xd t
= f x
d xd t
+ f y
d yd t
+ f z
d zd t
Introduction Ideal ow theory Basic aerodynamics
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
22/85
Convective derivative
Let f ( x, t) denote some quantity of interest in the uid motion.
Note that f
tmeans the rate of change of f at a xed
position x in space.The rate of change of f following the uid is
D f D t
=dd t
f x(t), t = f x
d xd t
+ f t
where x(t) = x(t), y(t), z(t) is understood to change with time atthe local ow velocity u , so that
d xd t
=d xd t
,d yd t
,d zd t
= u , v, w .
Therefore, f x
d xd t
= f x
u + f y
v+ f z
w = (u ) f .
Introduction Ideal ow theory Basic aerodynamics
d
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
23/85
Convective derivative
Let f ( x, t) denote some quantity of interest in the uid motion.
Note that f
tmeans the rate of change of f at a xed
position x in space.The rate of change of f following the uid is
D f D t
=dd t
f x(t), t = f x
d xd t
+ f t
f x
d xd t
= f x
u + f y
v+ f z
w = (u ) f .
Denition (Convective derivative)
The convective derivative (a.k.a. the material , particle ,Lagrangian , or total-time derivative ) is a derivative taken withrespect to a coordinate system moving with velocity u (i.e., followingthe uid):
D f D t
= f t
+ (u ) f .
Introduction Ideal ow theory Basic aerodynamics
C i d i i
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
24/85
Convective derivative
Denition (Convective derivative)
The convective derivative (a.k.a. the material , particle ,Lagrangian , or total-time derivative ) is a derivative taken withrespect to a coordinate system moving with velocity u (i.e., followingthe uid):
D f D t
= f t
+ (u ) f .
By applying the convective derivative to the velocity components u , v,w in turn it follows that the acceleration of a uid particle is
Du
D t=
u
t+ (u )u .
Introduction Ideal ow theory Basic aerodynamics
C i d i i
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
25/85
Convective derivative
By applying the convective derivative to the velocity components u , v,
w in turn it follows that the acceleration of a uid particle isDuD t
=ut
+ (u )u .
Example
Consider uid in uniform rotation with angular velocity , so that:
u = y , v = x , w = 0.
The ow is steady so ut = 0 , but
DuD t
= y
x+ x
y
y, x, 0 = 2 x, y, 0 .
represents the centrifugal acceleration 2 x2 + y2 towards therotation axis.
Introduction Ideal ow theory Basic aerodynamics
O tli
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
26/85
Outline
1 Introduction
Mathematical preliminariesPreliminary ideas and denitionsConvective derivative
2 Ideal ow theoryIdeal uidIncompressibility conditionEulers equations of motionBernoulli theoremsVorticity
3 Basic aerodynamicsSteady ow past a xed wingFluid circulation round a wingKuttaJoukowski theorem and conditionConcluding remarks
Introduction Ideal ow theory Basic aerodynamics
Id l id
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
27/85
Ideal uid
Denition ( Ideal uid )
Properties of an ideal uid :1 it is incompressible : no uid element can change in volume as
it moves;
Introduction Ideal ow theory Basic aerodynamics
Ideal uid
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
28/85
Ideal uid
Denition ( Ideal uid )
Properties of an ideal uid :1 it is incompressible : no uid element can change in volume as
it moves;2 it has constant density : is the same for all uid elements and
for all time t (a consequence of incompressibility);
Introduction Ideal ow theory Basic aerodynamics
Ideal uid
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
29/85
Ideal uid
Denition ( Ideal uid )
Properties of an ideal uid :1 it is incompressible : no uid element can change in volume as
it moves;2 it has constant density : is the same for all uid elements and
for all time t (a consequence of incompressibility);3 it is inviscid , so that the force exerted across a geometrical
surface element n S within the uid is p n S ,
where p ( x, t) is a scalar function of pressure , independent on thenormal n .
Introduction Ideal ow theory Basic aerodynamics
Ideal uid
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
30/85
Ideal uid
Denition ( Ideal uid )
Properties of an ideal uid :1 it is incompressible : no uid element can change in volume as
it moves;2 it has constant density : is the same for all uid elements and
for all time t (a consequence of incompressibility);3 it is inviscid , so that the force exerted across a geometrical
surface element n S within the uid is p n S ,
where p ( x, t) is a scalar function of pressure , independent on thenormal n .
RemarksAll uids are to some extent compressible and viscous.Air , being highly compressible , can behave like anincompressible uid if the ow speed is much smaller thanthe speed of sound .
Introduction Ideal ow theory Basic aerodynamics
Incompressibility condition
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
31/85
Incompressibility condition
S
V
Sn
Let S be a xed closed surface drawnin the uid, with unit outwardnormal n , enclosing a region V .
Introduction Ideal ow theory Basic aerodynamics
Incompressibility condition
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
32/85
Incompressibility condition
S
V
Sn
Let S be a xed closed surface drawnin the uid, with unit outwardnormal n , enclosing a region V .Fluid enters the enclosed region V atsome places on S , and leaves it atothers.
Introduction Ideal ow theory Basic aerodynamics
Incompressibility condition
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
33/85
Incompressibility condition
S
V
Sn
uLet S be a xed closed surface drawnin the uid, with unit outwardnormal n , enclosing a region V .Fluid enters the enclosed region V atsome places on S , and leaves it atothers.
The velocity component along theoutward normal is u n , so thevolume of uid leaving through asmall surface element S in unit timeis u n S .
Introduction Ideal ow theory Basic aerodynamics
Incompressibility condition
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
34/85
Incompressibility condition
S
V
Sn
uLet S be a xed closed surface drawnin the uid, with unit outwardnormal n , enclosing a region V .Fluid enters the enclosed region V atsome places on S , and leaves it atothers.
The velocity component along theoutward normal is u n , so thevolume of uid leaving through asmall surface element S in unit timeis u n S .The net volume rate at which uid isleaving V equals
S
u n dS
and must be zero for an
incompressible uid.
Introduction Ideal ow theory Basic aerodynamics
Incompressibility condition
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
35/85
Incompressibility condition
S
V
Sn
uNow, using the divergence theorem andthe assumption of the smoothness ofow (i.e., continuous velocity gradient)gives
Incompressibility condition
For all regions within an incompressibleuid
S
u n d S =V
u dV = 0
which (assuming that u is continuous)yields the following important constrainton the velocity eld:
u = 0
(everywhere in the uid).
Introduction Ideal ow theory Basic aerodynamics
Eulers equations of motion
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
36/85
Euler s equations of motionThe net force exerted on an arbitrary uid element is
S
p n dS = V
p dV
(the negative sign arises because n points out of S ). Now, providedthat p is continuous it will be almost constant over a small element V . The net force on such a small element due to the
pressure of the surrounding uid will therefore be p V .
Introduction Ideal ow theory Basic aerodynamics
Eulers equations of motion
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
37/85
u e s equat o s o ot oThe net force exerted on an arbitrary uid element is
S
p n dS = V
p dV
(the negative sign arises because n points out of S ). Now, providedthat p is continuous it will be almost constant over a small element V . The net force on such a small element due to the
pressure of the surrounding uid will therefore be p V .
The principle of linear momentum implies that the following forcesmust be equal for a uid element of volume V :
p + g V the total (external) net force acting on theelement in the presence of gravitational body force per unitmass g Nkg =
ms2 (the gravity acceleration),
V DuDt the inertial force , that is, the product of the elementsmass (which is conserved) and its acceleration.
Introduction Ideal ow theory Basic aerodynamics
Eulers equations of motion
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
38/85
q
The principle of linear momentum implies that the following forcesmust be equal for a uid element of volume V :
p + g V the total (external) net force acting on theelement in the presence of gravitational body force per unitmass g N
kg = m
s2 (the gravity acceleration),
V DuDt the inertial force , that is, the product of the elementsmass (which is conserved) and its acceleration.
This results in the Eulers momentum equation for an ideal uid.Another equation of motion is the incompressibility constraint.
Eulers equations of motion for an ideal uid
DuD t
= 1 p + g , u = 0
Introduction Ideal ow theory Basic aerodynamics
Eulers equations of motion
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
39/85
Boundary and interface-coupling conditions
Impermeable surface. The uid cannot ow through the boundary
which means that the normal velocity is constrained(no-penetration condition):
u n = u n .
Here, u n is the prescribed normal velocity of the impermeableboundary ( u
n= 0 for motionless, rigid boundary).
Introduction Ideal ow theory Basic aerodynamics
Eulers equations of motion
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
40/85
Boundary and interface-coupling conditions
Impermeable surface. The uid cannot ow through the boundary
which means that the normal velocity is constrained(no-penetration condition):
u n = u n .
Here, u n is the prescribed normal velocity of the impermeableboundary ( u
n= 0 for motionless, rigid boundary).
Free surface. Pressure condition (involving surface tension effects):
p = p0 2T mean .
Here: p0 is the ambient pressure,T is the surface tension,mean is the local mean curvature of the free surface.
Introduction Ideal ow theory Basic aerodynamics
Eulers equations of motiond d f l d
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
41/85
Boundary and interface-coupling conditions
Impermeable surface. The uid cannot ow through the boundary
which means that the normal velocity is constrained(no-penetration condition):
u n = u n .
Here, u n is the prescribed normal velocity of the impermeableboundary ( u
n= 0 for motionless, rigid boundary).
Free surface. Pressure condition (involving surface tension effects):
p = p0 2T mean .
Here: p0 is the ambient pressure,T is the surface tension,mean is the local mean curvature of the free surface.
Fluid interface. The continuity of normal velocity and pressure at theinterface between uids 1 and 2:
u (1) u (2) n = 0 , p (1) p(2) = 0 2 T (1) + T (2) mean .
Introduction Ideal ow theory Basic aerodynamics
Bernoulli theorems
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
42/85
The gravitational force , being conservative , can be written as thegradient of a potential :
g = where = gz m2
s2 ,
Introduction Ideal ow theory Basic aerodynamics
Bernoulli theorems
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
43/85
The gravitational force , being conservative , can be written as thegradient of a potential :
g = where = gz m2
s2 ,
and the momentum equationut
+ (u ) u
( u ) u + ( 12 u 2 )
= p
+
can be cast into the following form
ut
+ ( u ) u = H where H = p
+12
u 2 + .
Introduction Ideal ow theory Basic aerodynamics
Bernoulli theorems
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
44/85
u
t+ ( u ) u = H where H =
p+
1
2u 2 + .
For steady ow ( ut = 0 ): ( u ) u = H u
(u ) H = 0 .
Theorem (Bernoulli streamline theorem)
If an ideal uid is in steady ow, then H is constant along a streamline.
Introduction Ideal ow theory Basic aerodynamics
Bernoulli theorems
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
45/85
u
t+ ( u ) u = H where H =
p+
1
2u 2 + .
For steady ow ( ut = 0 ): ( u ) u = H u
(u ) H = 0 .
Theorem (Bernoulli streamline theorem)
If an ideal uid is in steady ow, then H is constant along a streamline.
Denition ( Irrotational ow )
A ow is irrotational if u = 0 .
Introduction Ideal ow theory Basic aerodynamics
Bernoulli theorems
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
46/85
u
t+ ( u ) u = H where H =
p+
1
2u 2 + .
For steady ow ( ut = 0 ): ( u ) u = H u
(u ) H = 0 .
Theorem (Bernoulli streamline theorem)
If an ideal uid is in steady ow, then H is constant along a streamline.
Denition ( Irrotational ow )
A ow is irrotational if u = 0 .
For steady irrotational ow :
H = 0 .
Theorem (Bernoulli theorem for irrotational ow)
If an ideal uid is in steady irrotational ow, then H is constant throughout the whole uid.
Introduction Ideal ow theory Basic aerodynamics
Vorticity: rotational and irrotational ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
47/85
Denition ( Vorticity )
Vorticity is a concept of central importance in uid dynamics; it isdened as
= u 1s .
Introduction Ideal ow theory Basic aerodynamics
Vorticity: rotational and irrotational ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
48/85
Denition ( Vorticity )
Vorticity is a concept of central importance in uid dynamics; it isdened as
= u 1s .
For an irrotational ow : = 0 .Vorticity has nothing directly to do with any global rotation of the
uid.
Introduction Ideal ow theory Basic aerodynamics
Vorticity: rotational and irrotational ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
49/85
Denition ( Vorticity )
Vorticity is a concept of central importance in uid dynamics; it isdened as
= u 1s .
For an irrotational ow : = 0 .Vorticity has nothing directly to do with any global rotation of the
uid.Interpretation of vorticity (in 2D)For two-dimensional ow (when u = u ( x, t), v( x, t), 0 ):
= 0, 0, where =v x
u y
.
Introduction Ideal ow theory Basic aerodynamics
Vorticity: rotational and irrotational ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
50/85
Denition ( Vorticity )
Vorticity is a concept of central importance in uid dynamics; it isdened as
= u 1s .
For an irrotational ow : = 0 .Vorticity has nothing directly to do with any global rotation of the
uid.Interpretation of vorticity (in 2D)For two-dimensional ow (when u = u ( x, t), v( x, t), 0 ):
= 0, 0, where =v x
u y
.
Then, at any point of the ow eld 2 represents the average angularvelocity of two short uid line-elements that happen, at that instant,to be mutually perpendicular.
Introduction Ideal ow theory Basic aerodynamics
Vorticity: rotational and irrotational ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
51/85
Interpretation of vorticity (in 2D)For two-dimensional ow (when u = u ( x, t), v( x, t), 0 ):
= 0, 0, where = v x
u y
.
A B
C
x
y
u x
x
v x
x
u y
y
v y
y
average angular velocity:AB + AC
2=
12
v x
u y
=2
Introduction Ideal ow theory Basic aerodynamics
Vorticity: rotational and irrotational ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
52/85
Interpretation of vorticity (in 2D)For two-dimensional ow (when u = u ( x, t), v( x, t), 0 ):
= 0, 0, where = v x
u y
.
A B
C
x
y
u x
x
v x
x
u y
y
v y
y
average angular velocity:AB + AC
2=
12
v x
u y
=2
Thus, the vorticity acts as a measure of the local rotation , or spin,
of uid elements.
Introduction Ideal ow theory Basic aerodynamics
Line vortex ow and uniformly rotating ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
53/85
Flow may be written in cylindrical polar coordinates (r , , z):
u = u r er + u e + u z e z .
Consider the following steady, two-dimensional ows (with u = u (r ),u r = u z = 0 , and being a constant):Line vortex ow: u = r e
Uniformly rotating ow: u = r e
u = u eu
1r
u r
Introduction Ideal ow theory Basic aerodynamics
Line vortex ow and uniformly rotating ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
54/85
Flow may be written in cylindrical polar coordinates (r , , z):
u = u r er + u e + u z e z .
Consider the following steady, two-dimensional ows (with u = u (r ),u r = u z = 0 , and being a constant):Line vortex ow: u = r e
Uniformly rotating ow: u = r e
u = u eu
1r
u r
Introduction Ideal ow theory Basic aerodynamics
Line vortex ow and uniformly rotating ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
55/85
Flow may be written in cylindrical polar coordinates (r , , z):
u = u r er + u e + u z e z .
Consider the following steady, two-dimensional ows (with u = u (r ),u r = u z = 0 , and being a constant):Line vortex ow: u = r e
Uniformly rotating ow: u = r e
u = u eu
1r
u r
Introduction Ideal ow theory Basic aerodynamics
Line vortex ow and uniformly rotating ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
56/85
Flow may be written in cylindrical polar coordinates (r , , z):
u = u r er + u e + u z e z .
Consider the following steady, two-dimensional ows (with u = u (r ),u r = u z = 0 , and being a constant):Line vortex ow: u = r e
Uniformly rotating ow: u = r e
u = u eu
1r
u r
Introduction Ideal ow theory Basic aerodynamics
Line vortex ow and uniformly rotating ow
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
57/85
Flow may be written in cylindrical polar coordinates (r , , z):
u = u r er + u e + u z e z .
Consider the following steady, two-dimensional ows (with u = u (r ),u r = u z = 0 , and being a constant):Line vortex ow: u = r e although the uid is clearly rotating in a
global sense the ow is in fact irrotational since = 0 (except atr = 0 , where neither u nor is dened);
Uniformly rotating ow: u = r e here, = 0, 0, 2 and avorticity meter is carried around as if embedded in a rigid body.
u = u eu
1r
u r
Introduction Ideal ow theory Basic aerodynamics
Rankine vortex
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
58/85
Denition ( Rankine vortex )
Rankine vortex is a steady, two-dimensional ow described as
u = u e with u = r for r a (uniformly rotating ow), a 2
r for r > a (line vortex ow),
where and a are constants.
Introduction Ideal ow theory Basic aerodynamics
Rankine vortex
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
59/85
Denition ( Rankine vortex )
Rankine vortex is a steady, two-dimensional ow described as
u = u e with u = r for r a (uniformly rotating ow), a 2
r for r > a (line vortex ow),
where and a are constants. Therefore,
= e z with =1r
(ru )r
1r
( r 2 )r = 2 for r a (vortex core ),
1r
( a 2 )r = 0 for r > a (irrotational).
Introduction Ideal ow theory Basic aerodynamics
Rankine vortex
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
60/85
r
u
a
a
u =
r for r a a 2r for r > a
r
2
a
= 2 for r a0 for r > a
The Rankine vortex serves as a simple idealized model for areal vortex.Real vortices are typically characterized by fairly small vortexcores (a is small) in which, by denition, the vorticity isconcentrated .Outside the core the ow is irrotational .
Introduction Ideal ow theory Basic aerodynamics
Vorticity equation
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
61/85
ut
+ u = H
t+ ( u ) = 0
t+ (u ) ( ) u + u
0 u
0= 0 .
Here, the fourth term vanishes because the uid is incompressiblewhile the fth term vanishes because = 0 .
Introduction Ideal ow theory Basic aerodynamics
Vorticity equation
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
62/85
ut
+ u = H
t+ ( u ) = 0
t+ (u ) ( ) u + u
0 u
0= 0 .
Here, the fourth term vanishes because the uid is incompressiblewhile the fth term vanishes because = 0 .
Vorticity equation
t+ (u ) = ( ) u , or
D
D t= ( ) u .
The vorticity equation is extremely valuable: as a matter of fact, itinvolves only u since the pressure has been eliminated and = u .
Introduction Ideal ow theory Basic aerodynamics
Vorticity equation
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
63/85
Vorticity equation
t + (u ) = ( ) u , orD
D t = ( ) u .
The vorticity equation is extremely valuable: as a matter of fact, itinvolves only u since the pressure has been eliminated and = u .
In the two-dimensional ow (u = u ( x, t), v( x, t), 0 , = 0, 0, ( x, t) )of an ideal uid subject to a conservative body force g the vorticity of each individual uid element is conserved :
DD t
= 0.
Introduction Ideal ow theory Basic aerodynamics
Vorticity equation
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
64/85
Vorticity equation
t + (u ) = ( ) u , orD
D t = ( ) u .
The vorticity equation is extremely valuable: as a matter of fact, itinvolves only u since the pressure has been eliminated and = u .
In the two-dimensional ow (u = u ( x, t), v( x, t), 0 , = 0, 0, ( x, t) )of an ideal uid subject to a conservative body force g the vorticity of each individual uid element is conserved :
DD t
= 0.
In the steady, two-dimensional ow of an ideal uid subject to aconservative body force g the vorticity is constant along astreamline :
(u ) = 0.
Introduction Ideal ow theory Basic aerodynamics
Outline
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
65/85
1 IntroductionMathematical preliminariesPreliminary ideas and denitionsConvective derivative
2 Ideal ow theoryIdeal uid
Incompressibility conditionEulers equations of motionBernoulli theoremsVorticity
3 Basic aerodynamicsSteady ow past a xed wingFluid circulation round a wingKuttaJoukowski theorem and conditionConcluding remarks
Introduction Ideal ow theory Basic aerodynamics
Steady ow past a xed wing
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
66/85
Steady ow past a wing at small angle of attack (incidence) istypically irrotational .
Introduction Ideal ow theory Basic aerodynamics
Steady ow past a xed wing
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
67/85
Steady ow past a wing at small angle of attack (incidence) istypically irrotational .
There are no regions of closed streamlines in the ow; all thestreamlines can be traced back to .
Introduction Ideal ow theory Basic aerodynamics
Steady ow past a xed wing
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
68/85
Steady ow past a wing at small angle of attack (incidence) istypically irrotational .
There are no regions of closed streamlines in the ow; all thestreamlines can be traced back to .The vorticity is constant along each streamline , and hence equal oneach one to whatever it is on that particular streamline at .
Introduction Ideal ow theory Basic aerodynamics
Steady ow past a xed wing
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
69/85
Steady ow past a wing at small angle of attack (incidence) istypically irrotational .
There are no regions of closed streamlines in the ow; all thestreamlines can be traced back to .The vorticity is constant along each streamline , and hence equal oneach one to whatever it is on that particular streamline at .As the ow is uniform at , the vorticity is zero on all streamlinesthere; hence, it is zero throughout the ow eld around the wing .
Introduction Ideal ow theory Basic aerodynamics
Steady ow past a xed wing
Typical measured pressure distribution on a wing in steady ow:
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
70/85
yp p g ythe pressures on the uppersurface are substantially lowerthan the free-stream value p ;the pressures on the lowersurface are a little higherthan p ;
in fact, the wing gets most of itslift from a suction effect on itsupper surface .
Introduction Ideal ow theory Basic aerodynamics
Steady ow past a xed wing
Typical measured pressure distribution on a wing in steady ow:
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
71/85
yp p g ythe pressures on the uppersurface are substantially lowerthan the free-stream value p ;the pressures on the lowersurface are a little higherthan p ;
in fact, the wing gets most of itslift from a suction effect on itsupper surface .
Why the pressures above the wing are less than those below?
The ow is irrotational and the Bernoulli theorem states that
H = p +12 u 2 is constant throughout 2D irrotational ows .
Explaining the pressure differences, and hence the lift on thewing, thus reduces to explaining why the ow speeds above thewing are greater than those below .An explanation is in terms of the concept of circulation .
Introduction Ideal ow theory Basic aerodynamics
Fluid circulation round a wing
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
72/85
Denition ( Circulation )
The circulation round some closed curve C lying in the uid regionis dened as =
C
u d x m2
s .
IfS
is the region enclosed by the curveC
then the Stokes theoremgives =
C
u d x =S
( u ) n dS ,
or in the two-dimensional context
=
C
u d x+ v d y =S
v x
u y
d xd y .
Notice that in the surface integrals vorticity terms appear.
Introduction Ideal ow theory Basic aerodynamics
Fluid circulation round a wing
= 0 if the closed curve C is spanned by a surface S which lies
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
73/85
p ywholly in the region of irrotational ow, that is, = 0 for anyclosed curve not enclosing the wing .This cannot be stated for any closed curve that does enclosethe wing . What can be stated is that such circuits have thesame value of .
Introduction Ideal ow theory Basic aerodynamics
Fluid circulation round a wing = 0 if the closed curve C is spanned by a surface S which lies
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
74/85
wholly in the region of irrotational ow, that is, = 0 for anyclosed curve not enclosing the wing .This cannot be stated for any closed curve that does enclosethe wing . What can be stated is that such circuits have thesame value of .
Proof.Consider two arbitrary closed curves, ABCA and DFED, enclosing a wing.
0 = ABCADEFDA= ABCA + AD + DEFD + DA
= ABCA + AD DFED AD ABCA = DFED
Introduction Ideal ow theory Basic aerodynamics
Fluid circulation round a wing
= 0 if the closed curve C is spanned by a surface S which lies
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
75/85
p ywholly in the region of irrotational ow, that is, = 0 for anyclosed curve not enclosing the wing .This cannot be stated for any closed curve that does enclosethe wing . What can be stated is that such circuits have thesame value of .
Therefore:
circulation round a wing is permissible in a steady irrotational ow.
Yet, two questions still arise:1 Why there should be any circulation?2 Why it should be negative, corresponding to larger ow speeds
above the wing than below?
The answers are given by the KuttaJoukowski condition .
Introduction Ideal ow theory Basic aerodynamics
KuttaJoukowski theorem and condition
Consider a steady, irrotational ow of uid round a wing. According to
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
76/85
Consider a steady, irrotational ow of uid round a wing. According toideal ow theory, the drag on the wing (the force per unit length of
wing parallel to the oncoming stream) is zero. What is the lift (theforce per unit length of wing perpendicular to the stream) is stated bythe following theorem.
Theorem (Kutta-Joukowski lift theorem )
Let be the uid density and U the ow speed at innity. Then, the lift of the wing is
F y = U Nm ,
where is the uid circulation around the wing.
Introduction Ideal ow theory Basic aerodynamics
KuttaJoukowski theorem and condition
Theorem (K tta Jo ko ski lift theorem )
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
77/85
Theorem (Kutta-Joukowski lift theorem )
Let be the uid density and U the ow speed at innity. Then, the lift of the wing is
F y = U Nm ,
where is the uid circulation around the wing.
Obviously, of a great importance is here the fact that = 0 . In thecase of a wing with a sharp trailing edge this can be explained asfollows: one good reason for non-zero circulation is that there wouldotherwise be a singularity (innity) in the velocity eld. This is statedby the KuttaJoukowski condition .
= 0 = K < 0
Introduction Ideal ow theory Basic aerodynamics
KuttaJoukowski theorem and condition
Obviously, of a great importance is here the fact that = 0 . In the
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
78/85
y, g pcase of a wing with a sharp trailing edge this can be explained as
follows: one good reason for non-zero circulation
is that there wouldotherwise be a singularity (innity) in the velocity eld.
= 0 = K < 0
KuttaJoukowski condition/hypothesis
The circulation is such that the ow leaves the trailing edge
smoothly , or, equivalently, that the ow speed at the trailingedge is nite .The ow speed is nite at the trailing edge only for one valueof the circulation around the wing: the critical value K. Thisparticular ow will correspond to the steady ow that is actuallyobserved.
Introduction Ideal ow theory Basic aerodynamics
KuttaJoukowski theorem and condition
KuttaJoukowski condition/hypothesis
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
79/85
The circulation is such that the ow leaves the trailing edge
smoothly , or, equivalently, that the ow speed at the trailingedge is nite .The ow speed is nite at the trailing edge only for one valueof the circulation around the wing: the critical value K. Thisparticular ow will correspond to the steady ow that is actuallyobserved.
The critical value K depends on the ow speed at innity U , and onthe size, shape, and orientation of the wing.
Thin, symmetrical wings
For a thin and symmetrical wing of length L, making an angle with the oncoming stream, the critical value of circulation is
K U L sin .
Introduction Ideal ow theory Basic aerodynamics
KuttaJoukowski theorem and condition
Thin symmetrical wings
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
80/85
Thin, symmetrical wings
For a thin and symmetrical wing of length L, making an angle with the oncoming stream, the critical value of circulation is
K U L sin .
Using this formula for the lifttheorem gives the following result
F y U 2 L sin ,
which is in excellent accord withexperiment provided that theangle of attack is small , that isonly a few degrees, depending onthe shape of the wing.
F y
c
Inviscid theory < c 6 12
Experiment
Introduction Ideal ow theory Basic aerodynamics
Concluding remarks
KuttaJoukowski hypothesis provides a rational explanation forh l d d h
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
81/85
the circulation round a wing in steady ight.It says nothing about the dynamical process by which thatcirculation is generated when a wing starts from a state of rest.
Introduction Ideal ow theory Basic aerodynamics
Concluding remarks
KuttaJoukowski hypothesis provides a rational explanation forh i l i d i i d i h
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
82/85
the circulation round a wing in steady ight.It says nothing about the dynamical process by which thatcirculation is generated when a wing starts from a state of rest.
Starting vortex
The circulation is generated by the so-called starting vortex , whichis a concentration of vorticity which forms at the trailing edge of
a wing as it accelerates from rest in a uid. It leaves the wing (whichnow has an equal but opposite bound vortex round it), and rapidlydecays through the action of viscosity.
Introduction Ideal ow theory Basic aerodynamics
Concluding remarks
KuttaJoukowski hypothesis provides a rational explanation forth i l ti d i i t d i ht
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
83/85
the circulation round a wing in steady ight.It says nothing about the dynamical process by which thatcirculation is generated when a wing starts from a state of rest.
Starting vortex
The circulation is generated by the so-called starting vortex , whichis a concentration of vorticity which forms at the trailing edge of
a wing as it accelerates from rest in a uid.
Question : Is a starting vortex theoretically explicable?
Introduction Ideal ow theory Basic aerodynamics
Concluding remarks
Question : Is a starting vortex theoretically explicable?
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
84/85
A legitimate conclusion on the basis of ideal ow theory:
If the wing and uid are initially at rest, the vorticity is initially zerofor each uid element.It remains zero since the vorticity is conserved for each uidelement.
Therefore there should be no starting vortex.
Introduction Ideal ow theory Basic aerodynamics
Concluding remarks
Question : Is a starting vortex theoretically explicable?
-
7/29/2019 ICMM TGZielinski IdealFlow.slides
85/85
A legitimate conclusion on the basis of ideal ow theory:
If the wing and uid are initially at rest, the vorticity is initially zerofor each uid element.It remains zero since the vorticity is conserved for each uidelement.
Therefore there should be no starting vortex.
An explanation of the starting vortex:
Ideal ow theory accounts well for the steady ow past a wing.The explanation of how that ow became established involvesviscous effects in a crucial way.But air, in some sense, is hardly viscous at all! Yet, viscouseffects are sufciently subtle that shedding of the vortex , whichis an essentially viscous process, occurs no matter how smallthe viscosity of the uid happened to be.