Chapter 1

Potential flow:
1.1

General formulation

Inviscid and irrotational flows in the limit of high Reynolds number are referred
to as potential or ideal flows. The term inviscid refers to flows where viscous
forces are small compared to inertial forces, so that the fluid viscosity can be
neglected in comparison to fluid inertia. Potential or ideal flows are a class of
inviscid flows in which the vorticity , which is the curl of the velocity vector,
is zero, i. e.
uk
=0
(1.1)
i = ijk
xj
Since the curl of the velocity is zero, the velocity can be expressed as the gradient
of a potential ,

ui =
(1.2)
xi
and hence the name potential flow.
Using equation 1.2 for the velocity field, the Navier-Stokes mass and momentum equations can be written in terms of the velocity potential. The mass
conservation equation, expressed in terms of the velocity potential, is,
ui
2
=
=0
xi
x2i

(1.3)

Thus, the mass conservation equation simply states that the velocity potential
satisfies the Laplace equation, 2 = 0. Therefore, for potential flow we can
use all the techiques developed for solving the Laplace equation earlier. The
momentum conservation equation an inviscid flow is,
ui
1 p
fi
ui
+ uj
=
+
t
xj
xi

(1.4)

where fi is the body force per unit volume acting on the fluid. The second term
on the right side of the above equation can be simplified for an irrotational flow
1

CHAPTER 1. POTENTIAL FLOW:

in which the vorticity is zero. Consider the cross product of the velocity and
vorticity, u = u u, which can be written in indicial notation as,
ijk uj klm

um
xl

= (il jm im jl )uj
= uj

ui
ui
ui
xj
xj

= uj

ui
1 (u2i )

xj
2 xj

um
xl

(1.5)

For an irrotational flow, the left side of the above equation is zero since the
vorticity u is zero. Therefore,
uj

1 (u2i )
ui
=
xj
2 xj

(1.6)

Using the above substitution for the second term on the left side of equation 1.4,
and also expressing the first term on the left in terms of the velocity potential
using equation 1.2, we obtain,
1 (u2i ) 1 p
fi

+
+

=0
t xi
2 xi
xi

(1.7)

If the body force fi is conservative, it can be expressed as the gradient of a

potential V as follows,
V
fi =
(1.8)
xi
When this is inserted into equation 1.7, and simplified, we obtain,



1 2 p V
=0
+ ui + +
xi t
2

(1.9)

Since all components of the gradient of the term in brackets is zero in the above
equation, the term in brackets has to be equal to a constant,


1 2 p V
p0
+ ui + +
(1.10)
=
t
2

where the constant p0 is the reference pressure at a location where the velocity
and potential V are zero. In a gravitational field, the potential is equal to gz,
where z is the height above the zero-potential reference. Equation 1.10 can then
be written as,


p0
1 2 p
+ ui + + gz =
(1.11)
t
2

Equation 1.11 is referred to as the Bernoulli equation for the potential flow.
An important point to note is that in potential flows, it is possible to satisfy
only the normal velocity and stress conditions at a surface, and it is not possible

1.1. GENERAL FORMULATION

to satisfy the tangential velocity and stress conditions. This is because we have
neglected the viscous terms in the momentum conservation equation 1.4, which
contain the second spatial derivatives of the velocity field. Consequently, we
have reduced the equations from a second order to a first order equation in
the spatial coordinates. And therefore, it is possible to satisfy only the normal
velocity and stress boundary conditions at the surface. If we neglect viscous
effects, the only contribution to the stress is the isotropic contribution due to
the pressure,
ij = pij
(1.12)
It should be noted that the Laplace equation 1.3 for the velocity potential
is a linear equation, and if normal velocity boundary conditions are prescribed
at the bounding surface, the potential is a linear function of the velocity. If the
velocity of all surfaces is changed by a constant factor, the potential and the
velocity at all points in the fluid also change by the same factor. In addition,
the solutions for the velocity potential are quasi-static for imposed normal
velocities of the bounding surfaces, because the potential at a given instant
in time depend only on the velocity of the bounding surfaces at that instant.
However, the Bernoulli equation 1.10 for the pressure is a non-linear equation,
since it has a contribution proportional to the square of the velocity, and is also
explicitly dependent on time, since it contains a term proportional to the time
derivative of the potential. Therefore, the solution for the velocity potential
under imposed normal stresses at the bounding surfaces is not linear or quasistatic.
There are some general results that can be derived for potential flows. The
first is that the kinetic energy in a potential flow can be expressed as a function
of the potential and normal velocity at the bounding surfaces. The kinetic
energy of the fluid in a volume V with bounding surfaces S is defined as,
Z

KE =
dV u2j
2 V


Z

=
dV
2 V
xj xj
!


Z
2

2
dV
=
(1.13)
2 V
xj xj
xj
where integration by parts has been used in the final step. Since the potential
satisfies the Laplace equation 1.3, the second term in the integrand on the right
side of equation 1.13 is zero, and the expression for the kinetic energy reduces
to,


Z

dV
KE =
2 V
xj xj
Z

=
uj nj
(1.14)
2 S

CHAPTER 1. POTENTIAL FLOW:

where n is the outward unit normal at the surfaces bounding the fluid. Therefore, the kinetic energy of the fluid can be expressed only in terms of the potential and normal velocity at the bounding surfaces. It follows that if the normal
velocity at all bounding surfaces is zero, then the kinetic energy is zero, which
implies that the fluid velocity is zero throughout the domain.
It is possible to prove a uniqueness theorem for potential flows, which
states that the potential flow solution is unique if the normal velocity boundary
conditions are specified on all boundaries. To prove this, we first assume that
the potential flow solution is not unique, and that there are two potential flow
solutions, u and u , which satisfy the same normal velocity boundary conditions,
ui ni = ui ni on the bounding surfaces of the flow. Consider the integral over
the fluid domain V ,
Z
I=

dV (ui ui )(ui ui )

(1.15)

Since the integrand is always positive, the integral also has to be positive. For a
potential flow, this integral can be rewritten by expressing the velocity in terms
of the velocity potential, and using integration by parts as follows.
Z
( )
dV (ui ui )
I =
xi
ZV
Z

(ui ui )
dV
=
dV ( )
((ui ui )( )
(1.16)
xi
xi
V
V

The second integral on the right side of equation 1.16 is zero due to the incompressibility condition, while the first integral can be expressed as a surface
integral using the divergence theorem,
Z
I =
dS( )(ui ui )ni
S

= 0

(1.17)

Here, we have used the condition that the two velocity fields satisfy the same
normal velocity conditions on the boundaries, so that (ui ui )ni is zero on all
boundaries.
We had earlier assumed that u and u are different velocity fields, in which
case the integral I has to be positive. However, we have proved in equation 1.17
that I is zero if the normal velocity fields at all surfaces are equal. Therefore,
this implies, that it is not possible to have two different potential flow solutions
which satisfy the same normal velocity boundary conditions at all the bounding
surfaces, and the solutions for the potential flow equations is unique if the normal
velocity conditions are specified at all boundaries.
Another general result is the minimum energy theorem, which states that
the kinetic energy of a potential flow is less than or equal to the kinetic energy
of any other flow that satisfies the same normal velocity boundary conditions
at the bounding surfaces. Let us consider two velocity profiles, ui and ui , both
of which satisfy the same normal velocity boundary conditions at the bounding

1.2. THREE-DIMENSIONAL POTENTIAL FLOWS

surfaces (ui ni = ui ni ). The velocity field ui is a solution of the potential

flow equations, whereas the velocity field ui is not necessarily a solution of
the potential flow condition, though it satisfies the mass conservation condition
(ui /xi ) = 0. The difference in the kinetic energies of the two flows is,
Z

2
dV (u2
KE KE =
i ui )
2 V
Z

=
dV (ui ui )2 + 2ui (ui ui )
2 V
Z
Z

2
dV ui (ui ui )
dV (ui ui ) +
=
2 V
V
Z
Z

=
(ui ui )
dV (ui ui )2 +
dV
2 V
x
i
V


Z
Z

(ui ui )
dV
=
dV (ui ui )2 +
((ui ui ))
(1.18)
2 V
xi
xi
V
where integration by parts has been used in the final step. The second term
on the second integral on the right side of equation 1.16 is zero because the
velocities ui and ui satisfy the incompressibility condition. Finally, the right
side of equation 1.16 can be simplified using the divergence theorem,


Z
Z

dV
dV (ui ui )2 +
((ui ui ))
KE KE =
2 V
xi
Z
ZV

dV (ui ui )2 + dSni (ui ui ))

(1.19)
=
2 V
S
Note the surface integral on the right side of equation 1.19 is zero, since it was
assumed that the velocities ui and ui satisfy the same normal velocity boundary
conditions at the surface. Therefore, the kinetic energy of the velocity field u i
that does not necessarily satisfy the potential flow equations is always equal to
or greater than that of velocity field ui which does satisfy the potential flow
equations,
Z

dV (ui ui )2
KE KE =
2 V
0
(1.20)

1.2
1.2.1

Three-dimensional potential flows

Motion of a sphere in an infinite fluid

The simplest three-dimensional potential flow is the motion of a sphere in a fluid

that it at rest at a large distance from the sphere. The configuration consists
of a sphere of radius R moving with a velocity U in a fluid that is at rest at
a large distance from the sphere, as shown in figure 1.1. The normal velocity

CHAPTER 1. POTENTIAL FLOW:

R
Figure 1.1: Sphere of radius R moving with a velocity U in a fluid which is at
rest at large distance from the sphere.
boundary condition for the fluid velocity ui at the surface of the sphere is,
ui n i = U i n i

(1.1)

The equation for the velocity potential, 1.3, has to be solved subject to the
boundary conditions 1.1 in order to determine the velocity field. Since equation
1.3 is linear, the potential is a linear function of the velocity of the particle. In
addition, the potential is a solution of the Laplace equation which decays to zero
at a large distance from the sphere, so it is a linear combination of the spherical
harmonic solutions. It is possible to construct only one scalar function which is
linear in the velocity and in one of the spherical harmonics,
=

AUj xj
r3

(1.2)

where A is a constant. The velocity field is then given by,

ui

xi

= AUj

ij
3xi xj

3
r
r5

(1.3)

The constant A is determined from the normal velocity boundary condition 1.1.
The normal fluid velocity at the surface is given by,
ui n i

u i xi
r
2AUj xj
=
r4
2AUj nj
=
r3

(1.4)

1.2. THREE-DIMENSIONAL POTENTIAL FLOWS

Equation the normal velocity 1.4 to Uj nj at the surface r = R, we obtain the

solution A = (R3 /2). Therefore, the final solution for the velocity potential
which satisfies the boundary conditions at the surface is,
=

R 3 U j xj
2r3

(1.5)

It is convenient to express the potential in a spherical co-ordinate system with

co-ordinates (r, , ), in which the azimuthal angle is measured from the direction of the velocity U. In this co-ordinate system, the potential is independent
of the meridional angle, and is given by,
=

R3 U cos ()
2r2

(1.6)

The components of the velocity in the r and directions are,

ur

r
R3 U cos ()
=
r3
=

=
=

1
r
R3 U sin ()
2r3

(1.7)

(1.8)

The total kinetic energy due to the fluid flow can be calculated using equation
1.14. Note that the unit normal n is the outward unit normal at the bounding
surfaces of the fluid. The integral in equation 1.14 has no contribution due to the
surface at a large distance from the sphere, since the velocity decreases to zero.
At the surface of the fluid, the outward unit normal to the fluid is directed
into the sphere, and therefore ni = (xi /r) at the surface of the sphere. In
addition, the normal velocity of the fluid at the surface uj nj is equal to the
normal velcoity of the sphere Uj nj , and therefore the kinetic energy is given by,
Z
Uj xj

dS
KE =
2 S
r
sphere
Z
R 3 U k xk U j xj

dS
=
2 S
r3
r
sphere
Z
R3 Uj Uk
xj xk
=
dS 4
(1.9)
4
r
S
sphere
The surface integral in equation 1.9 can easily be evaluated,
Z
4
xj xk
dS 4 =
r
3
S
sphere

(1.10)

CHAPTER 1. POTENTIAL FLOW:

Therefore the final expression for the kinetic energy due to the fluid motion is,


2R3 Uj2
(1.11)
KE =
3
2
If the kinetic energy of due to the moving fluid is expressed as (Ma Uj2 /2), where
Ma is the added mass which has to be moved in addition to the mass of the
moving fluid, then the added mass is given by
2R3
(1.12)
3
It can easily be seen that the added mass is equal to half the mass of fluid that
is displaced by the sphere.
The pressure exerted on the fluid at the surface of the sphere can be evaluated,
!
u2j
(1.13)
p = p0
+
t
2
Ma =

It is necessary to exercise some care in calculating the time derivative of the

potential in the above equation. Though the velocity of the sphere is a constant,
the potential in a fixed reference frame is not a constant. This is because the
radius vector from the center of the sphere to a fixed location in space changes as
the sphere moves, and consequently the potential at this location also changes.
If x0 is the location of the center of the particle, and x is the observation point,
the equation for the potential 1.2 can be written as,
=

R3 Uj (xj x0j )
2 |x x0 |3

(1.14)

Since the sphere is moving with velocity U, the rate of change of potential is,

dUj
dx0k
=
+
t
Uj dt
x0k dt

(1.15)

The derivative of the potential with respect to x0k can be written as the negative
of the derivative with respect to xk in the above equation, because the potential
is only a function of the distance (x x0 ). Therefore, the change in x x0
due to a displacement of x in the position x0 is identical to that due to a
displacement of x in the position x. Therefore, (/x0k ) = (/xk ),
and equation 1.16 reduces to,

=
=

Uj

Uj

dUj

Uk
dt
xk
dUj
U k uk
dt

This is inserted into the equation for the pressure to obtain,

!
u2j
dUj
u j Uj
p = p0
2
Uj dt

(1.16)

(1.17)

1.2. THREE-DIMENSIONAL POTENTIAL FLOWS

In the above equation, we have returned to a co-ordinate system in which the

origin is located at the center of the particle. The radial component of the fluid
velocity, ur , at the surface is zero, and the only contribution to the kinetic energy
is due to the tangential velocity at the surface, u (equation 1.8). Therefore,
the pressure at the surface r = R is given by,
 2

dU
u
p = p0
+ U sin ()u R cos ()
2
dt
!
2
R
dU
5 sin ()
2
(1.18)
2 cos ()
= p0 U
16
U
dt
In equation 1.19, we have substituted U = U sin () and u = (U sin ()/2)
at r = R.
The force exerted on the sphere can be calculated by integrating the pressure
over the surface of the sphere,
Z
dSij nj
Fi =
S
sphere
Z
dSpni
(1.19)
=
S
sphere
By symmetry, the force perpendicular to the velocity is zero. The force along
the direction of the velocity is,
Z
Fz =
dSpnz
(1.20)
S
sphere
where Fz is the force in the direction of the velocity vector, and nz is the
component of the unit normal to the sphere in the direction of the velocity
vector. For the steady flow of a particle with (dU/dt) = 0, it can easily be
shown that the force Fi is identically zero due to symmetry considerations. The
pressure p is a symmetric function of about the mid-plane cutting the sphere
perpendicular to the direction of flow at = (/2). However, the unit normal
nz = cos () is anti-symmetric about the plane = (/2), and therefore the
integral of the pressure and nz over the surface of the sphere is identically zero.
This result, that the steady motion of a sphere through a fluid does not exert a
force on the fluid, is referred to as dAlemberts paradox.
If the particle is accelerating, there is a net force exerted on the particle by
the fluid. This is calculated by integrating the contribution to the pressure due
to acceleration over the surface of the sphere.

 3

Z
R xj dUj 
Fi =
dSni
2r3 dt 
S
sphere
r=R
Z
dUj
=
dSxi xj
2R dt S
sphere

10

CHAPTER 1. POTENTIAL FLOW:

2R3
=
3
dUi
= Mb
dt


dUi
dt
(1.21)

Here, the outward unit normal to the sphere is given by ni = (xi /r), and we
have used the identity
Z
4R4 ij
dSxi xj =
(1.22)
3
S
sphere
to obtain the third step in equation 1.22. Not surprisingly, the force requried to
accelerate the particle through the fluid is the product of the added mass and
the acceleration.

1.2.2

General three-dimensional potential flows:

Some general results can be obtained for the form of the velocity potential for an
object of arbitrary shape translating through a fluid in potential flow. The solution for the velocity potential is, in general, more complicated than the solution
1.6 for the motion of a sphere. However, the leading contribution to the potential is still the dipole contribution which decays proportional to (1/r 2 ) in three
dimensions. There is no source term proportional to (1/r) due to the incompressibility condition and the constant volume of the object. However, higher
order terms which decay faster than (1/r 2 ) could be present in the solution.
Equation 1.18 for the pressure is valid for an object of arbitrary shape, since
no assumption was made regarding the specific form of the velocity potential
while deriving the pressure.
It is possible to show that the force required for the steady motion of an
object of arbitrary shape through the fluid,
Z
Fi =
dSpni
S
object
!
Z
u2j
dSni
=
(1.1)
U j uj
2
S
object
is equal to zero. Consider the volume integral
!
Z
u2j

dV
U j uj
xi
2
V

(1.2)

evaluated over the volume of the fluid surrounding the fluid. Using the divergence theorem, this volume integral can be written as,
!
! Z
!
Z
Z
u2j
u2j
u2j

U j uj
U j uj
U j uj
dV
=
dSni
dSni
xi
2
2
2
V
S
S
object

11

1.3. TWO-DIMENSIONAL POTENTIAL FLOWS:

Z

dSni
S

u2j
U j uj
2

Fi

The volume integral on the left side of equation 1.3 can further be simplified as
follows,
!
Z
Z
u2j

uj
dV
dV (uj Uj )
=
U j uj
x
2
xi
i
V
V
Z
ui
=
dV (uj Uj )
xj
ZV

((uj Uj )ui )
(1.4)
=
dV
xj
V
In the second step of the above equation, we have used the fact that the stress
tensor is symmetric for an irrotational flow, so that (ui /xj ) = (uj /xi ). In
the final step, we have used the incompressibility condition (uj /xj ) = 0. The
volume integral in equation 1.4 can further be converted into a surface integral
by first using the incompressibility condition
Z
Z
Z

dSnj (uj U(1.5)

dSnj (uj Uj )ui
dV
((uj Uj )ui ) =
j )ui
xj
S
V
S
object
In the final step of the above equation, we have used the normal velocity boundary condition, uj nj = Uj nj on the surface of the object, so that the integral
over the surface of the object is zero. Equating the right sides of equation 1.3
and 1.6, the force on the object can be expressed as,
!
!
Z
u2j
dS ni
Fi =
(1.6)
Uj uj nj (uj Uj )ui
2
S
The right side of the above equation contains terms that are linear or quadratic
in the velocity, and the integral is evaluated at a large distance from the object
in the limit r . In this limit, the velocity decreases proportional to (1/r 3 ),
whereas the surface area increases proportional to r 2 . Therefore, both the integrals on the right side are zero in the limit r . This shows that the force
on an object of arbitrary shape is zero in a potential flow.

1.3

Two-dimensional potential flows:

The analysis of potential flows in two dimensions is considerably simplified by

the use of complex functions. The simplification is due to two important properties of complex functions.
1. If a complex function F (z) of a complex variable z is analytic, i. e. the
change in the function DeltaF due to a small change in the variable z

(1.3)

12

CHAPTER 1. POTENTIAL FLOW:

can be expressed as

dF
z,
(1.1)
dz
then the real and imaginary parts of the complex function satisfy the
Laplace equation in the complex plane. Since the velocity potential is also
a solution of the Laplace equation, we can associate the real parts of all
complex functions to the velocity potential. Thus, there is a potential flow
field associated with all complex functions subject to the normal velocity
boundary condition at the surfaces.
F =

2. A transformation, called a conformal mapping, can be effected from the

independent complex variable z to some other variable z 0 in such a way
that the function F (z 0 ) also satisfies the Laplace equation in the z 0 plane.
In this way, it is possible to map the flow in a complicated domain onto
the flow in a much simpler domain, and solve the Laplace equation in the
simpler domain.
Both the above properties, and their application to two-dimensional complex
flows, will be explained in this section.
First, we will show that if the if a complex function is analytic, then the real
and imaginary parts satisfy the Laplace equation. Consider a complex function

F (z) which is a function of the independent variable z = x + y, where = 1.

The function can be written as the sum of its real and imaginary parts,
F (z) = (x, y) + (x, y)

(1.2)

The variation in the function F when the independent variable z is displaced

by z = x + y is,





F (z + z) F (z) =
x +
y
+
+
x
x
y
y





=
x +
y (1.3)
+
+
x
x
y
y
The function is analytic if the variation F is proportional to z in the limit
z 0, or if F can be written as,
F =

dF
z
dz

(1.4)

Therefore, for an analytic function, the real and imaginary parts of the coefficients of x and y in equation 1.3 are equal, and the functions and have
the relations,

(1.5)

(1.6)

1.3. TWO-DIMENSIONAL POTENTIAL FLOWS:

13

The above conditions are referred to as the Cauchy Reimann conditions for the
analyticity of a complex function. If we sum the partial derivative of equation
1.5 with respect to x and the partial derivative of equation 1.6 with respect
to y, we find that the function (x, y) satisfies the Laplace equation in two
dimensions,
2 2
+ 2 =0
(1.7)
x2
y
Similarly, if we take the difference of the partial derivative of equation 1.5 with
respect to y and equation 1.6 with respect to x, we find that the function also
satisfies the Laplace equation in two dimensions,
2 2
+
=0
x2
y 2

(1.8)

Therefore, if a complex function F (z) is analytic, then its real and imaginary
parts satisfy the Laplace equation.
Since the real part (x, y) of any complex function F (z) is a solution of
the Laplace equation 1.7, it is also a valid solution for the velocity potential in
a two dimensional domain provided it satisfies the normal velocity boundary
conditions at the prescribed boundaries. Thus, every complex function has a
two dimensional velocity field associated with it that satisfies the potential flow
equations. The velocities in the x and y direction are then given by,
ux

uy

=
x
y

=
y
x

(1.9)

where the Cauchy-Reimann conditions have been used to express derivatives

of the potential in terms of derivatives of the imaginary part . It is clear,
from the relations 1.9 between ux , uy and , that the imaginary part is the
stream function for the flow. Thus, the real and imaginary parts of the complex
potential are the velocity potential and the stream function respectively. The
function F (z) is often referred to as the complex potential. The derivative of
F (z) with respect to z can be written as,
dF
dz

=
+
x
x
= ux uy

W (z) =

(1.10)

It is also convenient to express the function W (z) in terms of the components

(ur , u ) in a polar co-ordinate system. The relation between the velocities in a
Cartesian and polar co-ordinate systems are,
ux
uy

= ur cos () u sin ()
= ur sin () + u cos ()

(1.11)

14

CHAPTER 1. POTENTIAL FLOW:

Using these it can easily be verified that the complex potential is,
W (z) = (ur u ) exp ()

(1.12)

Since we have shown that there are potential flow velocity fields associated
with any analytic complex function, we can examine the potential flow fields
due to some fundamental forms of the function F (z). We first consider the form,
F (z) = Az m

(1.13)

where A is assumed to be real without loss of generality. The derivative of F (z)

with respect to z is,
dF
dz
= mAz m1
= mArm1 exp ((m 1))

W (z) =

(1.14)

Comparing the form of W (z) in equation 1.14 with that in equation 1.12, we
find that the velocities ur and u in a polar co-ordinate system are,
ur

= Arm1 cos (m)

= Arm1 sin (m)

(1.15)

Equations 1.15 represent a flow that, for positive A, is directed radially outward
for = 0, radially inward for = (/m), and along the co-ordinate for
= (/2m), as shown in figure ??. Since the normal velocity is zero along the
lines = 0 and = (/m), the velocity field 1.15 represents the flow in a corner
of subtended angle (/m). Several special cases can be considered for specific
values of m. For m = 1, the subtended angle is , and the flow is a steady
linear flow with velocity A in the x direction, as shown in figure ??. The flow
for m = 2 is the stagnation point flow in a corner of angle (/2), as shown in
figure ??. The flow for m > 2 is the flow in a corner with a subtended acute
angle, while the flow for 1 < m < 2 is the flow in a corner with a subtended
obtuse angle, as shown in figure ??.
Next, we consider the form
m
log (z)
(1.16)
F (z) =
2
where m is a real constant. The derivative of F (z) with respect to z is,
m
W (z) =
2z
m
=
exp ()
(1.17)
2r
Comparing equation 1.17 with 1.12, we find that the components of the velocity
field are,
m
ur =
2r
u = 0
(1.18)

1.3. TWO-DIMENSIONAL POTENTIAL FLOWS:

15

This is the flow from a point source of fluid in two dimensions, from which the
total volume of fluid generated (per unit length in the third dimension) Q at
any radius R is equal to m, and is independent of radius.
A logarithmic function with an imaginary coefficient represents circulation
about the origin. Consider the form,

log (z)
2

F (z) =

(1.19)

where is real. The derivative of F (z) is,

2z

=
exp ()
2r

W (z) =

(1.20)

Comparing equation 1.20 with equation 1.12, we find the components of the
velocity,
ur

= 0

2r

(1.21)

Thus, the flow due to the potential 1.21 is a circulating flow around the origin.
For this flow, the circulation is the integral of the tangential velocity along a
closed curve around the origin, is equal to .
Z 2
Circulation =
rdu
0

1.3.1

(1.22)

Flow around a cylinder

The motion around a cylinder translating with a velocity Ux in the x direction

in a fluid which is at rest at a large distance from the cylinder is expressed by
a complex potential of the form,
F (z) =

Ux R 2

log (z)
z
2

(1.1)

The velocity fields can be inferred from the derivative of this complex potential
with respect to z,
dF
dz
Ux R 2

2
z
2z



Ux R 2
exp ()
=
exp ()
r2
2r

W (z) =

(1.2)

16

CHAPTER 1. POTENTIAL FLOW:

The radial and polar components of the velocity can be inferred using equation
1.12,
ur

Ux R 2
cos ()
r2
2
Ux R

sin () +
2
r
2r

(1.3)

The above velocity field satisfies the condition ui ni = Ui ni at the surface of

a cylinder with radius R, since the normal velocity ui ni , which is the radial
velocity ur in a polar co-ordinate system, is equal to Ui ni = Ux nx = Ux cos (),
where the component nx of the unit normal to the surface is equal to cos ().
The force on the cylinder can be determined from the equation 1.17 for the
force on an object moving with velocity Ui through the fluid,
!
Z
u2j
dSni
Fi =
U j uj
2
S
object
 2

Z 2

u
Ux (ur cos () u sin () 
=
(Rd) ni
(1.4)
2
0
r=R
The drag force Fx , which is in the direction of the velocity, can be determined
using nx = cos () in equation 1.4,

Fx

=
=

Z
Z

(Rd) cos ()

1
2

Ux sin () +
2R

2

(Rd) cos ()

1
2

Ux sin () +
2R

2

2
0
2
0

+ Ux (ur cos () + u sin ())

+

Ux2

= 0

cos () sin ()

!


(1.5)

As in three dimensions, the drag force due to the potential flow in two dimensions
is also equal to zero. The lift force perpendicular to the direction of gravity
can be determined using ny = sin () in equation 1.4,
Fx

(Rd) sin ()
0

= Ux

1.4

1
2

Ux sin () +
2R

2

Ux2

cos () sin ()

(1.6)

Force on a two-dimensional object of arbitrary shape:

The net force exerted on a two-dimensional object in potential flow can be

calculated in a manner similar to that for a three-dimensional flow. A procedure

1.4. FORCE ON A TWO-DIMENSIONAL OBJECT OF ARBITRARY SHAPE:17

identical to that in section ?? will result in the equation 1.6 for a two dimensional
object,
!
!
Z
u2j
(1.1)
Uj uj nj (uj Uj )ui
Fi =
dS ni
2
S
In two dimensions, the surface area (per unit length in the direction perpendicular to the plane of flow) of a surface with radius r increases proportional
to r in two dimensions. Therefore, the force is non-zero only if the integrand
decreases proportional to r in the limit r . In two dimensions, the most
slowly decaying velocity field is due to a line source of fluid equation 1.18, and
due to a line vortex equation 1.21, both of which decay proportional to (1/r)
in the limit r . In the absence of a source of fluid within the object, a nonzero contribution to the force can be caused only by a line vortex with velocity
u = (/2r) where is the circulation. co-ordinates, or ui = ijk j nk /(2r)
in Cartesian co-ordinates, where is the circulation, and i = iz . In equation
1.1, there is no contribution due to the terms proportional to u2j and ui uj , since
these decay proportional to (1/r) in the limit r . Therefore, it is sufficient
to consider the components of the integrand of equation 1.1 linear in the velocity
ui ,
Z
dSUj (nj ui ni uj )
(1.2)
Fi =
S

If we choose the direction of the mean velocity to be along the x direction

without loss of generality, it is clear that the drag force along the x direction
is zero from equation 1.2. However, the lift force along the y direction is not
zero, and can be calculated by writing the velocity u = (/2r) (equation
1.2) in terms of the Cartesian components, ux = ( sin ()/(2r)) and uy =
( cos ()/(2r)),
Z
Fy =
dSUx (nx uy ny ux )
S
Z
=
dSUx /(2r)
S

= U

(1.3)

Problems:
1. Find the added mass per unit length of an infinite cylinder in potential
flow, using a procedure identical to that for a three-dimensional object
derived here.
2. Determine the total kinetic energy due to a moving sphere under potential
flow conditions. From this, calculate the added mass. Also, find the
total kinetic energy of the fluid when a sphere moves in the limit of zero
Reynolds number. How do the two compare?

18

CHAPTER 1. POTENTIAL FLOW:

3. Determine the energy dissipation due to the motion of a sphere at high
Reynolds number, assuming that the fluid velocity field is given by the
potential flow solution. Using this, find the drag force on the sphere. How
does it compare to the drag force on a sphere at zero Reynolds number?
The rate of dissipation of energy is given by,
Z
D = dV (u:u)
4. Consider a bubble with internal pressure pb expanding in a fluid in which
the pressure a large distance from the bubble is p0 . The radius of the
bubble R(t) is a function of time as it expands.
(a) Solve the potential flow equations to determine the fluid velocity field
due to the expanding bubble.
(b) Determine the pressure at the surface of the bubble.
(c) Form a pressure balance condition, find the equation for the evolution
of the bubble radius.
5. Consider two line sources of strength m and m separated by a distance
d along the x axis, in a fluid flowing with a constant velocity U in the
x direction. What is the equation for the shape of the object which is
equivalent to these two sources? Show that in the limit d 0 and (md)
finite, the object assumes the shape of a cylinder, and find its radius.
6. Consider the transform from the z to the z 0 plane given by
z = z0 +

a2
z0

(a) Determine the relations between the coordinates (x0 , y 0 ) and (x, y).
(b) Consider a circle of radius a in the z 0 plane. What is the transformed
shape in the z plane? If we consider the flow past the cylinder in the
z 0 plane
a2
F (z 0 ) = z 0 + 0
z
in the z 0 plane, what is the equivalent flow in the z plane?
(c) Consider a circle of radius b > a in the z 0 plane. What is the transformed shape in the z plane? If we consider the flow past the cylinder
in the z 0 plane
b2
F (z 0 ) = z 0 + 0
z
in the z 0 plane, what is the equivalent flow in the z plane?

1.4. FORCE ON A TWO-DIMENSIONAL OBJECT OF ARBITRARY SHAPE:19

7. A three dimensional irregular body is moving with a velocity U1 in the x1
direction near a wall which is perpendicular to the x2 axis as shown in the
figure. The Reynolds number is large so that the potential flow equations
are applicable. The wall is impermeable so that the fluid flow at the wall
is tangential to it. Find the force on the body in the x1 and x2 directions
as a function of the velocity of the object U1 and the fluid velocity at the
wall.
(Hint: Consider a control volume bounded by the surface of the object
S, the surface of the wall Sw and the surface at infinity S . Use methods
similar to the derivation of the dAlembert paradox for a general body in
irrotational flow.)
8. Consider the dynamics of waves on the surface of a liquid of wavelength
, frequency and amplitude 0 . Find the conditions (at high Reynolds
number) under which the uj j ui term in the momentum conservation
equation can be neglected compared to the t ui term. Under these conditions, find the frequency of the surface waves on the surface of a liquid
of infinite depth as a function of the wavelength of the waves. The height
of the surface fluctuations is given by the equation:
= 0 exp (ikx + it)
where k = (2/) is the wave number, and is the frequency of the waves.
In addition, the normal velocity at the surface z = 0 is given by the time
rate of change of the displacement .