Potential Flow:: 1.1 General Formulation
Potential Flow:: 1.1 General Formulation
Potential flow:
1
2 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,
∂um ∂um
ǫijk uj ǫklm = (δil δjm − δim δjl )uj
∂xl ∂xl
∂ui ∂ui
= uj − ui
∂xj ∂xj
∂ui 1 ∂(u2i )
= uj − (1.5)
∂xj 2 ∂xj
For an irrotational flow, the left side of the above equation is zero since the
vorticity ∇ × u is zero. Therefore,
∂ui 1 ∂(u2i )
uj = (1.6)
∂xj 2 ∂xj
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 (1.7)
∂t ∂xi 2 ∂xi ρ ∂xi ρ
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
+ ui + + =0 (1.9)
∂xi ∂t 2 ρ ρ
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,
∂φ 1 2 p p0
+ 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 3
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 = −pδij (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 quasi-
static.
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
!
∂2φ
ρ ∂ ∂φ
Z
= dV −φ 2 (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
KE = dV
2 V ∂xj ∂xj
ρ
Z
= φuj nj (1.14)
2 S
4 CHAPTER 1. POTENTIAL FLOW:
where n is the outward unit normal at the surfaces bounding the fluid. There-
fore, the kinetic energy of the fluid can be expressed only in terms of the poten-
tial 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 = u∗i ni on the bounding surfaces of the flow. Consider the integral over
the fluid domain V , Z
I= dV (u∗i − ui )(u∗i − ui ) (1.15)
V
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
I = dV (u∗i − ui )
∂xi
ZV
∂ ∂(u∗i − ui )
Z
= dV ((u∗i − ui )(φ∗ − φ) − dV (φ∗ − φ) (1.16)
V ∂xi V ∂xi
The second integral on the right side of equation 1.16 is zero due to the in-
compressibility condition, while the first integral can be expressed as a surface
integral using the divergence theorem,
Z
I = dS(φ∗ − φ)(u∗i − 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 (u∗i − 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 u∗i , both
of which satisfy the same normal velocity boundary conditions at the bounding
1.2. THREE-DIMENSIONAL POTENTIAL FLOWS 5
ρ
Z
2
KE∗ − KE = dV (u∗2
i − ui )
2 V
ρ
Z
dV (u∗i − ui )2 + 2ui (u∗i − ui )
=
2 V
ρ
Z Z
∗ 2
= dV (ui − ui ) + ρ dV ui (u∗i − ui )
2 V V
ρ ∂φ ∗
Z Z
= dV (u∗i − ui )2 + ρ dV (ui − ui )
2 V V ∂x i
ρ ∂ ∂(u∗i − ui )
Z Z
= dV (u∗i − ui )2 + ρ dV (φ(u∗i − ui )) − φ (1.18)
2 V V ∂xi ∂xi
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 u∗i satisfy the incompressibility condition. Finally, the right
side of equation 1.16 can be simplified using the divergence theorem,
ρ ∂
Z Z
KE∗ − KE = dV (u∗i − ui )2 + ρ dV (φ(u∗i − ui ))
2 V ∂xi
ZV
ρ
Z
= dV (u∗i − ui )2 + ρ dSni φ(u∗i − 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 u∗i 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
KE∗ − KE = dV (u∗i − ui )2
2 V
≥ 0 (1.20)
θ U
R
boundary condition for the fluid velocity ui at the surface of the sphere is,
u i ni = U i ni (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
φ= (1.2)
r3
where A is a constant. The velocity field is then given by,
∂φ
ui =
∂xi
δij 3xi xj
= AUj 3
− (1.3)
r r5
The constant A is determined from the normal velocity boundary condition 1.1.
The normal fluid velocity at the surface is given by,
ui xi
u i ni =
r
2AUj xj
= −
r4
2AUj nj
= − (1.4)
r3
1.2. THREE-DIMENSIONAL POTENTIAL FLOWS 7
1 ∂φ
uθ =
r ∂θ
R3 U sin (θ)
= (1.8)
2r3
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,
ρ φUj xj
Z
KE = − dS
2 S r
sphere
ρ R3 Uk xk Uj xj
Z
= dS
2 S r3 r
sphere
ρR3 Uj Uk xj xk
Z
= dS 4 (1.9)
4 S r
sphere
The surface integral in equation 1.9 can easily be evaluated,
xj xk 4π
Z
dS 4 = (1.10)
S r 3
sphere
8 CHAPTER 1. POTENTIAL FLOW:
Therefore the final expression for the kinetic energy due to the fluid motion is,
2πR3 ρ Uj2
KE = (1.11)
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
2πR3 ρ
Ma = (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 evalu-
ated, !
∂φ u2j
p = p0 − ρ + (1.13)
∂t 2
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 )
φ=− (1.14)
2 |x − x0 |3
Since the sphere is moving with velocity U, the rate of change of potential is,
∂φ ∂φ dUj ∂φ dx0k
= + (1.15)
∂t ∂Uj dt ∂x0k dt
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,
∂φ ∂φ dUj ∂φ
= − Uk
∂t ∂Uj dt ∂xk
∂φ dUj
= − Uk u k (1.16)
∂Uj dt
This is inserted into the equation for the pressure to obtain,
!
u2j ∂φ dUj
p = p0 − ρ − u j Uj (1.17)
2 ∂Uj dt
1.2. THREE-DIMENSIONAL POTENTIAL FLOWS 9
2πR3 ρ
dUi
=
3 dt
dUi
= Mb (1.21)
dt
Here, the outward unit normal to the sphere is given by ni = (xi /r), and we
have used the identity
4πR4 δij
Z
dSxi xj = (1.22)
S 3
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.
evaluated over the volume of the fluid surrounding the fluid. Using the diver-
gence theorem, this volume integral can be written as,
! ! Z !
∂ u2j u2j u2j
Z Z
dV − Uj u j = dSni − Uj u j − dSni − Uj u j
V ∂xi 2 S∞ 2 S 2
object
1.3. TWO-DIMENSIONAL POTENTIAL FLOWS: 11
!
u2j
Z
dSni − Uj u j − Fi (1.3)
S∞ 2
The volume integral on the left side of equation 1.3 can further be simplified as
follows,
!
∂ u2j ∂uj
Z Z
dV − Uj u j = dV (uj − Uj )
V ∂x i 2 V ∂xi
∂ui
Z
= dV (uj − Uj )
∂xj
ZV
∂
= dV ((uj − Uj )ui ) (1.4)
V ∂xj
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
dV ((uj − Uj )ui ) = dSnj (uj − Uj )ui − dSnj (uj − U(1.5)
j )ui
V ∂xj S∞ S
object
In the final step of the above equation, we have used the normal velocity bound-
ary 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,
! !
u2j
Z
Fi = dS ni − Uj uj − nj (uj − Uj )ui (1.6)
S∞ 2
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/r3 ),
whereas the surface area increases proportional to r2 . Therefore, both the inte-
grals 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.
can be expressed as
dF
∆F = ∆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.
2. A transformation, called a conformal mapping, can be effected from the
independent complex variable z to some other variable z ′ in such a way
that the function F (z ′ ) also satisfies the Laplace equation in the z ′ 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,
∂ψ ∂φ
= − (1.6)
∂x ∂y
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 (1.8)
∂x2 ∂y 2
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 = =
∂x ∂y
∂φ ∂ψ
uy = =− (1.9)
∂y ∂x
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
W (z) =
dz
∂φ ∂ψ
= +ı
∂x ∂x
= ux − ıuy (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 = ur cos (θ) − uθ sin (θ)
uy = 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
W (z) =
dz
= mAz m−1
= mArm−1 exp (ı(m − 1)θ) (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 = Arm−1 cos (mθ)
uθ = Arm−1 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
F (z) = log (z) (1.16)
2π
where m is a real constant. The derivative of F (z) with respect to z is,
m
W (z) =
2πz
m
= exp (−ıθ) (1.17)
2πr
Comparing equation 1.17 with 1.12, we find that the components of the velocity
field are,
m
ur =
2πr
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,
ıΓ
F (z) = log (z) (1.19)
2π
where Γ is real. The derivative of F (z) is,
ıΓ
W (z) = −
2πz
ıΓ
= − exp (−ıθ) (1.20)
2πr
Comparing equation 1.20 with equation 1.12, we find the components of the
velocity,
ur = 0
Γ
uθ = (1.21)
2πr
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 = rdθuθ
0
= Γ (1.22)
Ux R 2 ıΓ
F (z) = − − log (z) (1.1)
z 2π
The velocity fields can be inferred from the derivative of this complex potential
with respect to z,
dF
W (z) =
dz
Ux R 2 ıΓ
= 2
−
z 2πz
Ux R 2
ıΓ
= exp (−ıθ) − exp (−ıθ) (1.2)
r2 2πr
16 CHAPTER 1. POTENTIAL FLOW:
The radial and polar components of the velocity can be inferred using equation
1.12,
Ux R 2
ur = cos (θ)
r2
2
Ux R Γ
uθ = 2
sin (θ) + (1.3)
r 2πr
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,
!
u2j
Z
Fi = ρ dSni − Uj u j
S 2
object
Z 2π 2
uθ
= ρ (Rdθ) ni − Ux (ur cos (θ) − uθ sin (θ) (1.4)
0 2 r=R
The ‘drag’ force Fx , which is in the direction of the velocity, can be determined
using nx = cos (θ) in equation 1.4,
2 !
2π
1 Γ
Z
Fx = ρ (Rdθ) cos (θ) Ux sin (θ) + + Ux (ur cos (θ) + uθ sin (θ))
0 2 2πR
2 !
2π
1 Γ
Z
2 2 2
= ρ (Rdθ) cos (θ) Ux sin (θ) + + Ux cos (θ) − sin (θ)
0 2 2πR
= 0 (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,
2 !
2π
1 Γ
Z
2 2 2
Fx = ρ (Rdθ) sin (θ) Ux sin (θ) + + Ux cos (θ) − sin (θ)
0 2 2πR
= ρUx Γ (1.6)
identical to that in section ?? will result in the equation 1.6 for a two dimensional
object, ! !
u2j
Z
Fi = dS ni − Uj uj − nj (uj − Uj )ui (1.1)
S∞ 2
In two dimensions, the surface area (per unit length in the direction perpen-
dicular 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 non-
zero contribution to the force can be caused only by a line vortex with velocity
uθ = (Γ/2πr) where Γ is the circulation. co-ordinates, or ui = ǫijk Γj nk /(2πr)
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
Fi = ρ dSUj (nj ui − ni uj ) (1.2)
S∞
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:
(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.
a2
z = z′ +
z′
(a) Determine the relations between the coordinates (x′ , y ′ ) and (x, y).
(b) Consider a circle of radius a in the z ′ plane. What is the transformed
shape in the z plane? If we consider the flow past the cylinder in the
z ′ plane
a2
F (z ′ ) = z ′ + ′
z
in the z ′ plane, what is the equivalent flow in the z plane?
(c) Consider a circle of radius b > a in the z ′ plane. What is the trans-
formed shape in the z plane? If we consider the flow past the cylinder
in the z ′ plane
b2
F (z ′ ) = z ′ + ′
z
in the z ′ plane, what is the equivalent flow in the z plane?
1.4. FORCE ON A TWO-DIMENSIONAL OBJECT OF ARBITRARY SHAPE:19
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 ξ.