Scribd
Upload
45 views

NPTEL Web Course On Complex Analysis: A. Swaminathan

This document is a module on complex analysis from the NPTEL web course taught by A. Swaminathan and V.K. Katiyar. It discusses two-dimensional fluid flow and stream functions. Specifically, it defines the velocity and rotation of a fluid, derives that the velocity potential must satisfy Laplace's equation, and introduces the stream function as the harmonic conjugate of the velocity potential. Streamlines represent curves of constant stream function value.

Uploaded by

Aditya Koutharapu
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
45 views

NPTEL Web Course On Complex Analysis: A. Swaminathan

This document is a module on complex analysis from the NPTEL web course taught by A. Swaminathan and V.K. Katiyar. It discusses two-dimensional fluid flow and stream functions. Specifically, it defines the velocity and rotation of a fluid, derives that the velocity potential must satisfy Laplace's equation, and introduces the stream function as the harmonic conjugate of the velocity potential. Streamlines represent curves of constant stream function value.

Uploaded by

Aditya Koutharapu
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 21

NPTEL web course

on
Complex Analysis

A. Swaminathan
I.I.T. Roorkee, India

and

V.K. Katiyar
I.I.T. Roorkee, India

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 1 / 21


Complex Analysis

Module: 9: Applications of Conformal Mapping


Lecture: 3: Two dimensional flow and Stream function

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 2 / 21


Applications of Conformal Mapping

Two-dimensional fluid flow

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 3 / 21


Two-Dimensional Fluid Flow

In this section we see the application of conformal mapping in


Two-Dimensional Fluid Flow.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 4 / 21


Two-Dimensional Fluid Flow

Let us assume that the motion of the fluid is the same in all planes
parallel to the xy plane and the velocity is parallel to that plane
and independent of time. So it is sufficient to consider the motion
of a sheet of fluid in xy the plane.
Let V = p + iq be the velocity of a particle at any point (x, y ),
where p(x, y ), q(x, y ) are real valued continues functions and its
first order partial derivatives are also continuous at points which
are interior to a region of flow without any sources or sinks of the
fluid.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 5 / 21


Two-Dimensional Fluid Flow

Let C be a contour. Then the circulation of the fluid along C is


defined as Z
VT (x, y )dσ
C
where σ is the arc length of the tangential component VT (x, y ) of
the velocity vector along C.
Now, the above integral can be written as
Z Z
VT (x, y )dσ = p(x, y )dx + q(x, y )dy
C C

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 6 / 21


Two-Dimensional Fluid Flow

If C is a positively oriented simple closed contour which lies inside


a simply connected domain of flow containing no sources or sinks,
then using Green’s theorem, we have
Z Z Z
p(x, y )dx + q(x, y )dy = [qx (x, y ) − py (x, y )]dA
C R

where R is the closed region consisting of points interior to and on


C.
Hence,
Z Z Z
VT (x, y )dσ = [qx (x, y ) − py (x, y )]dA
C R

for such a contour.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 7 / 21


Two-Dimensional Fluid Flow

Let us consider a circle counterclockwise at x0 , y0 and denote it by


C. Now if we divide the circumference 2πr , we will get mean
speed along C. To get the corresponding angular speed of the
fluid about the centre of the circle, we have to divide the mean
speed by r : Z Z
1 1
2
[qx (x, y ) − py (x, y )]dA
πr R 2
So, mean value of the function
1
ω(x, y ) = [qx (x, y ) − py (x, y )]dA
2
over the circular region R and can be expressed by the above
integral.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 8 / 21


Two-Dimensional Fluid Flow

Consequently, ω(x, y ) is called the rotation of the fluid which


represents the limiting angular speed of a circular element of the
fluid as the circle shrinks to its centre (x,y), the point at which ω is
evaluated.
The flow is said to be irrational in a simply connected domain if
ω(x, y ) = 0 in each point of the domain. Here we have considered
only irrational flows of an incompressible fluid which is free from
viscosity. It can be shown that the fluid pressure p(x, y ) satisfies
the following special cases of Bernouli’s equation

p 1 2
+ |V | = c
ρ 2

where c is constant.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 9 / 21


Two-Dimensional Fluid Flow

It can be observed that the pressure is greatest where the speed


|V | is least.
Now we consider an irrational flow in a simply connected domain
D. From the value of ω(x, y )

py = qx
R
in D. Hence, C p(s, t)ds + q(s, t)dt along a contour C lying
entirely in D and if we join two points x0 , y0 and x, y in D, it will be
independent of the path.
So, for fixed x − 0, y0 ,
Z (x,y )
φ(x, y ) = p(s, t)ds + q(s, t)dt
(x0 ,y0 )

is well defined on D.
A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 10 / 21
Two-Dimensional Fluid Flow

Taking partial derivatives on both sides of φ(x, y ), we have

φx (x, y ) = p(x, y ), φy (x, y ) = q(x, y )

According to φx (x, y ), velocity vector V = p + iq is the gradient of


φ and the directional derivative of φ in any direction represents the
component of the velocity flow in that direction. φ(x, y ) is called
velocity potential.
From φ(x, y ), it is clear that when the reference point (x0 , y0 ) is
changed, φ(x, y ) also changes. The curves φ(x, y ) = c1 are called
equipotentials.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 11 / 21


Two-Dimensional Fluid Flow

y
V

1 x

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 12 / 21


Two-Dimensional Fluid Flow

In the case of the heat flow, the condition that the incompressible
fluid enter or leave an element of volume only by floating through
the boundary of that element requires that φ(x, y ) must satisfy
Laplace equation

φxx (x, y ) + φyy (x, y ) = 0

in a domain when the fluid is free from sources or sinks.


From φ(x, y ) and the continuity of p, q and its first order partial
derivatives it is clear that first and second order partial derivatives
of φ are continuous in such a domain. consequently, the velocity
potential φ is Harmonic function in that domain.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 13 / 21


Applications of Conformal Mapping

Stream function

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 14 / 21


Stream Function

In this section, we see the application called stream functions.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 15 / 21


Stream Function

Let φ be the velocity potential and V = p(x, y ) + iq(x, y ) be the


velocity vector for a simply connected domain in which the flow is
irrational.
Then
V = φx (x, y ) + iφy (x, y ) = gradφ(x, y ).
It is normal to an equipotential passing through the point (x, y ) if
the velocity vector is not the zero vector.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 16 / 21


Stream Function

Let ψ(x, y ) be the harmonic conjugate of φ(x, y ) and the velocity


vector is tangent to a curve ψ(x, y ) = c2 , which are called stream
lines of the flow and ψ is called the stream function. A particular
example can be given as a boundary across which fluid cannot
flow is a streamline.
The complex potential of a flow is given by

F (z) = φ(x, y ) + iψ(x, y )


⇒ F 0 (z) = φx (x, y ) + iψx (x, y )

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 17 / 21


Stream Function
Using Cauchy-Riemann equations, we have

F 0 (z) = φx (x, y ) + iφy (x, y )

From V = φx (x, y ) + iφy (x, y ) = gradφ(x, y ),

V = F 0 (z)

The speed or magnitude is given by

|V | = |F 0 (z)|.

Since φ is harmonic in a simply connected domain D. A harmonic


conjugate of φ is given by
Z (x,y )
ψ(x, y ) = −φt (s, t)ds + φs (s, t)dt
(x0 ,y0 )

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 18 / 21


Stream Function

(x, y )
V

O x

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 19 / 21


Stream Function

So,
Z
ψ(x, y ) = −q(s, t)ds + p(s, t)dt
C

where C is any contour in D from (x0 , y0 ) to (x, y ).


We know that the right-hand side of ψ(x, y ) represents the integral
with respect to the arc length σ along C of normal component
VN (x, y ) of the vector whose x and y components are p(x, y ) and
q(x, y ) respectively.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 20 / 21


Stream Function

Hence, ψ(x, y ) can be written as


Z
ψ(x, y ) = VN (s, t)dψ.
C

Physically ψ(x, y ) represents the rate of flow by volume across a


surface of unit height which stands perpendicular to the xy plane
on the curve.

A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 21 / 21

You might also like