Potential Flow
Potential Flow
Potential Flow
1.1 Introduction:
Mathematical analyses of problems in generally possible only if
certain simplifying assumptions are made. And complete solution
of problems can seldom be obtained without resource to
experiment.
The most common simplifying assumption is that the fluid is
ideal / perfect, thus eliminating the complicating viscous effects.
The other may be assuming that the flow doesn’t change with
time.
1
Two methods of study of the fluid flow
A. In the Lagrangian method
a single particle is followed over the flow
field, the co-ordinate system following the
particle.
B. In the Eularian method,
the description of flow is on fixed coordinate
system based and the description of the
velocity etc. are with reference to location and
time i.e., V = V (x, y, z, t) and not with
reference to a particular particle.
2
3
Types of fluid motion
4
I. Laminar and Turbulent flow
6
III. Uniform and Non-uniform flow
• Uniform flow – the flow said to be uniform
when there is no variation in the magnitude
and direction of the velocity vector from one
point to another along the path flow.
i.e ds
dV
0
dP
ds
0
7
IV. Incompressible and Compressible
flow
8
V. Irrotational and Rotational flow
9
One, two and three dimensional flow
(1D, 2D and 3D flow)
• One- Dimensional flow (1D flow) – in one dimensional flow, all the flow parameters
like velocity, pressure etc may be expressed as function of time and one space
coordinate only V, P ……. etc…. f (x, t)
• Two – Dimensional flow (2D flow) - in these case the flow parameters are the
function of time and two space coordinates (x and y) only. i.e f(x, y, t)
• Three – Dimensional flow (3D-flow) - in 3D- flow the parameters are function of
time and three coordinates system. i.e. f(x, y, z, t)
13
Definition of Terms
• Streak line – a dye or smoke is frequently injected into
a fluid in order to trace its subsequent motion. The
resulting dye or smoke trails are called streak line.
In the steady flow of a
liquid, a colored dye reveals
the streamlines.
14
Stream tube, Path lines and Streak
lines
The flow lines passing through the edge of
an imaginary area such as A is called a
Stream tube
16
BASIC SCIENTIFIC LAWS USED IN THE
ANALYSIS OF FLUID FLOW
• Laws of Mechanics
– Written for systems
– System = arbitrary quantity
of mass of fixed identity
– Fixed quantity of mass, m
Q: Have you ever used your thumb to control the water flowing
from the end of a hose?
A: When the end of a hose is partially closed off, thus reducing its
cross-sectional area, the fluid velocity increases.
This kind of fluid behavior is described by the equation of
continuity.
18
Example 1.1
and
19
Equation of Continuity
20
Example 1.2
and
21
Equation of Continuity
• Conservation of
mass is applied to
the control-volume
element δx δy δz
(Fig. 1.2) with
center at (x, y, z),
• where the velocity
components in the
x, y, z directions
are u, v, w,
respectively, and ρ
is the density.
22
Equation of Continuity
• Consider first the flux through the pair of faces normal to the x direction.
On the right-hand lace the flux outward is
• Since both ρ and u are assumed to vary continuously throughout the fluid.
• ρu δy δz is the mass flux through the center face normal to the x axis.
• The second term is the rate of increase of mass flux, with respect to x multiplied by
the distance δx/2 to the right-hand face.
23
On the left-hand lace
• The other two directions yield similar expressions; the net mass outflow (MNET)
is
24
the rate of change of the control volume is
the net mass outflow is equal to the rate of change of the control volume.
Then
and if their sum does equal zero, we call the equation satisfying continuity.
26
The Bernoulli Equation
for constant density the Bernoulli equation yields
When it is divided by g,
27
1.2 Circulation, Vorticity, Definition of
potential flow
• Considering a closed path in a flow field as
shown in Fig.a,
• CIRCULATION is defined as the line integral of
velocity about this closed path. The symbol
used is or K=
28
circulation
29
Definition of Vorticity
• Vorticity is defined as circulation per unit
area. i.e.,
• Vorticity = circulation per unit area, here area
is dx dy, so
30
Example 1.4
• Consider the two-dimensional velocity distribution u = –By, v = +Bx, where
B is a constant, evaluate the circulation Γ around the rectangular closed
curve defined by (x, y) = (1, 1), (3, 1), (3, 2), and (1, 2).
31
Vorticity
32
N.B For irrotational flow, vorticity and circulation are both zero.
33
34
35
36
Occurrence of irrotational and rotational regions for flow
through a pipe.
37
1.3 Basic equations for a 2-D potential flow
The Stream Function
Stream function is a mathematical expression that describes a flow
field.
The stream function is defined as
For 2 D Flow in Re c tan gular coordinate For 2 D flow in Polar Coordinates
1
u ;v ur ;V
y x r r
Does it satisfy the equation of continuity? YES
u v 2 2
0
x y xy xy
If the flow is irrotational (no vorticity) then it satisfies Laplace
u v 2 2
Z 2 2 0
y x y x
and is therefore additive 1 2 3 ...... 38
Properties of Stream Function
A) =constant is a streamline
B) d between two streamlines is proportional to the
Volumetric Flow
Let us consider a line given by , a constant
as shown in Fig. We have
i.e.,
giving
The scalar function of space and time such that its derivative with respect to any
direction gives the fluid velocity in that direction is called velocity potential .
In inviscid, incompressible, irrotational the continuity equation breaks down to:-
u v 2 2
.V . 2
2 2 0 it satisfies Laplace Equation
x y x y
and is therefore additive 1 2 3 ......
These are second order Laplace equations for which there are many solutions and
analytical techniques known to solve them 41
Revise
42
Laplace Equation and is encountered in many branches
of physics and engineering.
A flow governed by this equation is called a Potential
Flow.
43
Example 1.5
.
44
Example 1.6
•Given that
, determine the stream function and potential function for the flow described.
45
Comparison
• Velocity potential φ
• Stream function
• Flow field velocities by
• By differentiating normal to the
other?
Eq. 1 showing that equipotential lines and streamlines are orthogonal to each
other. This enables one to calculate the stream function when the velocity
potential is given and vice versa.
47
Stream lines and Equi-potential lines for flow through a bend
48
1.4 Velocity potential and stream functions of elementary
potential-flow models:
There are only a few types of flow for which stream and
are described
Examples of elementary potential-flow models:-
Uniform flow,
source/sink flow,
vortex; 49
Elementary Planar Irrotational Flows
I. Uniform Flow with constant velocity U
The flow can be described by the condition, u = constant
and v =0
50
Generalised linear flow(Uniform Flow)
u 0 y v0 x
u 0 U cos
v0 U sin
y
= U(y cos – x sin ) U
v0
u0
N.B A uniform flow is a physically
possible incompressible flow (i.e it
satisfies .V=0) and that is x
52
On Integrating
Sources
y
m is the source strength in m2s-1 r
m x
(r , )
2
m y
( x, y ) tan 1
2 x
m x m y
u ; v
y 2 x y
2 2
x 2 x 2 y 2
m The radial velocity decreases inversely with distance from
ur
2r the source, but the flow rate, m, remains constant
m x m x m y m y
u u r cos ; v u r sin 54
2 r 2
2r r 2 r 2
2r r
Source at the origin
3
1
=0
0 y
-1
-2
-3
-4
-5 -4 -3 -2 -1 0 1 2 3 4 5
x 55
Sources not at the origin
For a source in a generalised position (x0,y0):
m 1 y y0
tan
2 x x0
m x x0
u
2 x x0 2 y y 0 2
m y y0
v
2 x x0 2 y y 0 2
56
Sink Flow
Sink is the opposite of source and the radial streamlines are
directed inwards to a common point, origin, where the fluid is
absorbed at a constant rate.
57
Sinks
Basically a source with a negative m y
r
m is the source strength in m2s-1
m x
(r , )
2
m 1 y
( x, y ) tan
2 x
m x m y
u ; v
y 2 x y
2 2
x 2 x 2 y 2
m
ur
2r
m x m y
u ur cos ; v ur sin
2 r 2
2 r 2
58
III) Vortex
•A vortex has the pathlines being circles centered on
the origin, and fluid particles move along these circles.
•The strength of a vortex is measured by the circulation.
On Integrating
59
Whirlpool
60
Elementary Planar Irrotational Flows
Line Vortex
• If vortex is moved to
(x,y) = (a,b)
Example 1.7
• A source with strength 0.2 m3/s·m and a vortex with strength 1 m2/s
are located at the origin. Determine the equations for velocity potential
and stream function. What are the velocity components at x = 1 m, y =
0.5 m?
Solution:
• The velocity potential for the source is
62
• The velocity potential for the vortex is
The streamline pattern for a doublet is shown in the right figure above.
64
Doublet - Theory
m m
(r , ) d
2 2
y
md
(r , )
2
d
r
sin d sin r-dr
By the sine rule:
ds r dr +d
ds sin source ds sink x
As ds 0, sin d d and r dr r d
r
md mds sin sin
Then (r , ) K
2 2 r r
A Doublet is formed when the source and sink approach each other,i.e a0 &
m.Such that mds/2π=constant
The strength of the doublet is K=mds/2π.=constant
65
Doublet - Theory y
r
mds x
K is the doublet strength
2
In cartesian coordinates:
y 2 x2
u u ( x, y ) K
y x 2
y
2 2
xy
v v ( x , y ) 2 K
x x y
2 2 2
66
Doublet - example
K= 100m3s-1
13
12
11
10
9 14-16
8 12-14
7 10-12
6
8-10
5
4 6-8
3 4-6
2
2-4
1
0
0-2
-1 -2-0
-2 -4--2
-3
-6--4
-4
-5
-8--6
-6 -10--8
-7 -12--10
-8
-14--12
-9
-10 -16--14
-11
-12
10
12
14
-8
-6
-4
-2
8
-20
-18
-16
-14
-12
-10
67
1.5 Superposition of elementary potential-flow models
70
I) Source + Uniform Flow = Flow Past a
Half Body
uniform uy m
2
source
For the combined flow ,
m m
Uy Ur sin
2 2
71
stagnation point for flow past Half Body
At a point x=−b the velocities due to the two flows cancel each
other, and this is identified as the stagnation point.
Uniform flow velocity U=radial outward flow due to source Ur
(r=b)
m m m
U U b m 2 bU
2 r 2 b 2 U
r
72
Flow Past a Half Body
Stagnation streamline(Dividing stream line)
The stagnation streamline defines the shape of the aerodynamic
‘body’.
The value of the stream function on this streamline is the same as
the value at the stagnation point(y=0,=) since Y=r sin
m m. m
Stagnation streamline, s u 0 y u 0 .0 bU
2 2 2
For y=0 on the x-axis. There are two solutions, one for the upper and
one for the lower streamlines.
s=m/2
m
s
2
s=-m/2 73
Flow Past a Half Body
The stagnation streamline (the dividing streamline that passes
through the stagnation point is given by s)
74
Flow Past a Half Body
V
The velocity components for this flow are given by
U
If the pressure in the free stream is it follows from Bernoulli Equation that
75
Source and horizontal flow
=3 3
=2 2
=1 1
=0
=0 0 y
=-1 -1
=-2 -2
=-3 -3
=-4 -4
-5 -4 -3 -2 -1 0 1 2 3 4 5
x 76
(Example 1.8)Rankine half-body
Streamlines may not be
crossed by the fluid, so the
streamline can represent a
body in the flow given by:
y
y 10 tan 1
x 10
Determine the distance to the
stagnation point (b)?
77
Solution
m 1 y y0 uniform u0 y
source tan
2 x x0
y y0
For a source and horizontal flow: U 0 y m tan 1
2 x x0
m x x0
u u u0
y 2 x x0 2 y y0 2
At a stagnation point u=0
m x x0
u0
2 x x0 2 y y0 2
If the stagnation point is on the x-axis, y=y0 m
x x0
2u0
For the example over where:- m
u0 1m / s , 0 10m and 10m 2 / s
y 2
y 10tan 1
x 10 10 x 20m b
x 10 78
Example 1.9 Y
81
Rankine Bodies
13
12
=10 11
10
=8 9
8
=6 7
6
10-12
8-10
=4 5
6-8
4
=2 3
2
4-6
2-4
1
=0 0 0-2
-1 -2-0
=-2 -2
-3
-4--2
-4
-6--4
=-4 -5 -8--6
-6 -10--8
=-6 -7
-12--10
-8
=-8 -9
=-10 -10
-11
-12
-20
-18
-16
-14
-12
-10
10
12
14
-8
-6
-4
-2
82
Rankine Bodies
linear flow + source+ sink
m 1 y y0 m 1 y y1
u0 y tan
tan
2 x x0 2 x x1
u
y
u u0
m x x0
m x x1
2 x x0 y y0 2 x x1 2 y y1 2
2 2
83
Rankine Bodies - stagnation points
y
s s
-a a x
To find the stagnation points, put y0=0, y1=0 and x0=-a, x1=a
m x a x a
0 u0
2 x a x a 2
2
84
Rankine Bodies - stagnation points
m 1 1
0 u0
2 x a x a
m x a x a
0 u0
2 x 2 a 2 x 2 a 2
m 2a
0 u0
2 x 2 a 2
ma
x2 a2
u0
85
Rankine Bodies – example 1.11
• A Rankine body is formed with a source and
sink, 5 m apart, each of strength 10 m2/s, and
a uniform flow of 5 m/s as shown in Fig. Find
the length L and thickness t of the body.
86
Solution
V=u0= 5ms-1, m=10m2s-1, a=2.5m ma
x a
2 2
u0
10 2.5
x 2 2.52
5
x 2.8m
The body is therefore L=5.6m long.
m 1 y y0 m 1 y y1
u0 y tan
tan
2 x x0 2 x x1
10 1 t /2 10 1 t /2
0 5 (t / 2 ) tan tan
2 0 ( 2 .5 ) 2 0 2 .5
10 1 t
0 5 (t / 2 ) tan
5 87
1.5.2 Lifting flow over a cylinder and
an airfoil, Kutta-Joukowski theorem
As the source and the sink combine to become a doublet, the
Rankine Oval becomes a circular cylinder.
As the flow past a circular cylinder is of fundamental interest, let
us examine the flow in some detail.
A. Doublet + Uniform Flow = Flow Past a Circular Cylinder
88
Doublet + Uniform Flow = Flow Past a Circular
Cylinder
(U=V, Free stream velocity)
89
Examples of Irrotational Flows Formed by
Superposition(Flow over a cylinder)
Laminar
separation
Irrotational
flow
B) Free Vortex + Doublet + Uniform Flow =
(Flow Past a Rotating Circular Cylinder )
Lifting flow over a cylinder
91
Lifting flow over a cylinder
Consequently the velocity components will be,
The tangential velocity on the surface of the cylinder (when r=a) is given by,
The streamline pattern for this flow depend upon the location of the stagnation points given by,
92
.
Stagnation Points for a lifting circular cylinder
The stagnation points we saw in Fig. above are for the case when the
circulation imposed on the cylinder was such that (see fig b)
With zero circulation the stagnation points lie at = o,2 (see fig.a)
If circulation is further increased the stagnation point will no longer be found on
the cylinder surface, but will appear in the flow as shown in Fig.d
93
(Flow Past a Rotating Circular Cylinder ) Lifting flow over a cylinder (cont’d)
94
Kutta-Joukowski theorem
:- It states that lift per unit span on a two-dimensional body is directly proportional to the
circulation around the body.
:-It is simply an alternative way of expressing the consequences of the surface pressure
distribution; it is a mathematical expression that is consistent with the special tools we
have developed for the analysis of inviscid, incompressible flow.
95
Summary
KEY CONCEPTS
• Elementary Flows (Building Blocks, why such a name?)
1. Uniform Flow
2. Source / Sink Flow
3. Doublet Flow
4. Vortex Flow
• What is the purpose? → Simulate real shapes in a simple manner
– Combine (1) + (2) → flow over half-body or oval
– Combine (1) + (3) → flow over a cylinder
– Combine with (4) → flow over a lifting cylinder
– Kutta-Joukowski Theorem
– Combinations of sources, vortex, uniform flow.
– Why can we combine so easily (simply add)?
• Know how to set up and for all cases and combined flows (no
time to solve)
• Know how to get velocity components u and v
• How would you model some basic shapes using these tools?
96
SUMMARY OF STREAM AND POTENTIAL FUNCTIONS
TABLE
R r
2
V r sin 1 2 ln
r 2 R
L V Kutta-Joukowski Theorem
98
Summary
Irrotational (potential) Flow Approximation
2D Flows
99
The Finite Wing Theory
100
HOW DOES AN AIRFOIL GENERATE LIFT?
• Lift due to imbalance of pressure distribution over top and bottom
surfaces of airfoil (or wing)
– If pressure on top is lower than pressure on bottom surface, lift is generated
– Why is pressure lower on top surface?
Lift = PA
Introduction
Streamtube A is squashed
most in nose region
(ahead of maximum thickness)
A
B
102
Introduction
103
AIRFOILS VERSUS WINGS
How do we design?
What is limit of behavior?
Introduction
HOW DOES WING GENERATE LIFT?
When the wing produces positive lift, the pressure of the
lower surface of wing is bigger than the pressure of the upper
surface.
Under the function of pressure difference between the upper
surface and lower surface of wing, the airflow of lower surface
flows around wing tip to the upper surface and forms the tip
eddy.
Look forward from the back of wing, the left tip eddy rotates
clockwise and the right tip eddy rotates anti-clockwise.
105
a. The creation of induced drag
the tip eddy (vortex):
Y
107
Application of Vortex flow pattern
The flow over an airfoil is two-dimensional. The flow over the finite wing is three-dimensional
w w
downwash speed
108
Lift on an airplane wing…. Flow lines are more
crowded above the wing, which means the velocity of
the air moves faster on top. Increased velocity means
lower pressure on top, producing an upward force
________
we call “Lift”. If “Lift” > weight
_________, the plane “takes off”,
accelerating up.
If angle of wing
becomes too
great, plane can
“stall”.
109
But, is this REALLY the correct “story” to explain
how airplanes can fly??
Believe it or not, there’s some controversy over
the issue…
http://www.amasci.com/wing/airfoil.html
Flow around an airfoil: the dots move
with the flow. Note that the velocities
are much higher at the upper surface
than at the lower surface. The black
dots are on timelines, which split into
two — an upper and lower part — at
the leading edge. The part of a
timeline below the airfoil does not
catch up with the one above. Colors of
the dots indicate streamlines
From Wikipedia.com 110
Aerodynamic force includes: lift, drag and side force .
Aerodynamic force definition : The acting force of air on an object when the object
moves relatively with the air.
111
Section One: The Characteristics
and main geometric parameters of Wing
Airfoil
112
I. The Shape of Wing
(I)the commonly used Airfoil
113
(II)The commonly used plane shape of wing
114
Section Two: The parameters of wing
115
Camber line: The
connection line of a series
of incenters of the airfoil
116
2.Relative camber
• Maximum camber :Maximum distance between the
mean camber line of the airfoil and the wing chord
fmax
bb
Relative camber :the ratio of the maximum camber
(fmax) and the wing chord (b)
117
3.Wing area:
The projected area of wing
on plane xoz is called wing
area which is represented
by S。
118
4.Wingspan:
The distance between
the left and right wing
end (wingtip) is called
wingspan (L) 。
PT-6 aircraft,
L =10.22m。
119
5.Aspect ratio:
The ratio of wingspan to
the average wing chord
(baverage) is called aspect
ratio (λ)
L / b L2
PT-6 aircraftav, λ=6;/S
120
6.Wing taper ratio:
the ratio of wing root
chord length to wing tip
chord length , is called
wing taper ratio (η )。
rectangular wing:η=1
Delta wing η=∞。
121
The influencing factors of the Lift
There are three factors which influencing
the Lift:
Y=Cy(1/2)ρV2S
(1)wing area S
Lift is in direct proportion along with the
wing area S.
122
(2) Dynamic pressure
Dynamic pressure includes two factors: air
density and airflow speed.
( 3) Lift coefficient Cy
123
The changing regularity of the lift
coefficient Cy along with the angle of
attack α
the relationship between α and Cy can be drawn
into a diagram which is called lift coefficient
curve.
Cymax
α αc
124
0 r