Lecture 3 Part 1 Potential Flow
Lecture 3 Part 1 Potential Flow
Chapter Three
PART ONE
u v
+
=0
x y
u v
+
= (
) + (
)=0
This is true for any smooth
x y x y
y
x
function (x,y)
so that it always satisfies the continuity eq.
using stream
u
function
two unknows
v
} one
unknow
c
conservation of mass will always be satisfied
Why do this?
Physical significance
1.
2.
Change in along
streamline is zero
Let dq represent the volume rate of flow per unit width perpendicular
to the xy plane passing between the two streamlines.
From conservation of mass we know that the inflow, dq, crossing the
arbitrary surface AC must equal the net outflow through surfaces
AB and BC. Thus,
u = 2y
v = 4x
u=
= 2y
y
= y 2 + f1 (x)
v=
= 4x
x
= 2x 2 + f 2 (y)
= 2x2 + y 2 + C
For simplicity, we set C = 0
Example Solution
=0
= 2x 2 + y 2
y2 x 2
=1
/2
tan
lim
lim
z =
2 x y
x =
2 y z
Rotation Vector
1 u w
2 z x
y =
r
1 w v u w v u
= x i + y j + z k =
i +
j + k
2 y z z x x y
2=
w v u w v u
i +
j + k
z
x
y
z
x y
Mathematical Representation
v
i
v
V
x
u
v
j
y
v
v
k
v w v
= i
+
z
y z
w
v u w v v u
j
+ k
z x x y
v
i
v
v
k ( V )
x
u
v
j
y
v
v
k
v u
=
z x y
w
Circulation ( )
r
r r
= V dl = V cos dl
NOTE:
The flow is irrotational if
=0,
=0,
=0
d = udx + ( v +
v
u
) dxdy
y
x
= z . area
v
u
dx ) dy ( u +
dy ) dx vdy
x
y
= (
= . area
Potential Function,
Irrotational approximation:
vorticity is negligibly small
Cartesian
Cylindrical
,
x
,
ur =
r
u=
, w=
z
y
1
, uz =
u =
r
z
v=
Momentum equation
If we can compute from the Laplace equation (which came from
continuity) and velocity from the definition
, why do we
need the NSE? the answer: To compute Pressure.
To begin analysis, apply irrotational approximation to viscous term
of the incompressible NSE
=0
Inviscid
Irrotational ( = 0)
p V2
p V
+
+ gz =
+
+ gz
2
2
Valid for 3D or 2D
u=
v=
z =
v u
=0
x y
= 2
This proves that the Laplace equation holds for the stream
function and the velocity potential
Constant values of :
streamlines
Constant values of :
equipotential lines
and are mutually
orthogonal
is defined by continuity;
2 results from
irrotationality
is defined by irrotationality;
2 results from continuity
u=
v=
d =
dx +
dy = udx + vdy
x
y
And
u
dy
=
v
dx =c
1
dy
=
(dy dx) =c
dx =c
Planar
Axisymmetric
Method of Superposition
1.
Since 2 = 0 is linear, a linear combination of two or more
solutions is also a solution, e.g., if 1 and 2 are solutions,
then (A1), (A+1), (1+2), (A1+B2) are also solutions
2.
3.
u = U=constant , v = 0, w = 0
In Cartesian coordinates
= U x , = U y
Conversion to cylindrical
coordinates can be achieved
using the transformation
= U r cos , = U r sin
Stagnation Flow
= u = Bx,
x
r
V = Bx i By j
= v = By
y
Therefore,
B 2
B 2
2
(
x
y
)
=
r cos 2
And
2
2
= u = Bx,
= v = By
y
x
Therefore
B 2
= Bxy = r sin 2
2
vr
vr =
K
and v = 0
2r
K
=
=
2r r r
=
and v = 0 =
r
r
vr =
= 0 and
r
K
=
r 2r
By integration:
K
ln r
=
2
=0
r
K
K
and
=
=
r 2r
2
By integration:
K
=
If source/sink is moved
to (x,y) = (a,b)
K
K
K
1 $
=
=0
"
$
1
2%
"
& =
&
Free Vortex
Line Vortex
circulation
= V ds =
C
= K
= K ln r
ds = d = 0
K
(rd ) = 2 K
r
ln r
2
Stokess Theorem
The transformation from a line integral to a surface integral in
three-dimensional space is governed by Stokess theorem:
Potential Vortex
The curl of the velocity vector for the potential vortex can be
found using the definition for the curl of V S in cylindrical
coordinates
=
2%
Paths for the calculation of the circulation for a potential vortex: (a)
closed curve C1, which encloses origin; (b) closed curve C2, which does
not enclose the origin.
=
1
=
" 2%
Bernollis equation
p V2
+
+ gz = const
2
at the free surface p=0:
V12 V2 2
=
+z
2g 2g
2
z= 2 2
8 r g