[go: up one dir, main page]

0% found this document useful (0 votes)
121 views5 pages

MIT2 25F13 Solution6.4 PDF

Uploaded by

funk singh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
121 views5 pages

MIT2 25F13 Solution6.4 PDF

Uploaded by

funk singh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

MIT Department of Mechanical Engineering

2.25 Advanced Fluid Mechanics

Kundu & Cohen 6.4


This problem is from “Fluid Mechanics” by P. K. Kundu and I. M. Cohen

(a) Take a plane source of strength m at point (−a, 0), a plane sink of equal strength at (a, 0), and
superpose a uniform stream U directed along the x-axis.
(b) Show that there are two stagnation points located on the x-axis at points

m 1/2
±a +1 .
πaU

(c) Show that the streamline passing through the stagnation points is given by ψ = 0. Verify that the line
ψ = 0 represents a closed oval-shaped body, whose maximum width h is given by the solution of the
equation
πU h
h = a cot .
m
The body generated by the superposition of a uniform stream and a source-sink pair is called a Rankine
body. It becomes a circular cylinder as the source–sink pair approach each other.

2.25 Advanced Fluid Mechanics 1 Copyright @


c 2010, MIT
Potential Flow Kundu & Cohen 6.4

Solution:

(a)

W (z) = Wuniform flow + Wsource + Wsink

where W = φ + iψ, φ is the potential function, and ψ the stream function.

Recap from Lecture: W satisfies the Laplace equation which is linear. Therefore, one can superimpose
its solutions as above.

Wuniform flow = U∞ (x + iy) y arbitrary


  point (x, y)
m m mθ
Wsource = ln reiθ = ln r + i r
2π 2π 2π

  θ
m m mθ
y0
Wsink = − ln reiθ = − ln r + i
2π 2π 2π location of
source/sink
Substitute expressions for r and θ in terms of x

and y (see figure):


x0 x
  _
m m mθ
Wsource = ln reiθ = ln r + i r= (x − x0 )2 + (y − y0 )2
2π 2π 2π  
  y − y0
m m mθ θ = arctan
Wsink = − ln reiθ = − ln r + i x − x0
2π 2π 2π
   
m (x + a)2 + y 2 m y m y
⇒ Wtotal = U∞ x + ln 2 2
+i U∞ y + arctan − arctan
4π (x − a) + y 2π x+a 2π x−a
φ ψ

(b) Obtain the velocity field (vx , vy ) by invoking v = ∇φ

2 ��
∂φ m � −�
(x� a)� + y2 2(x + a) (x + a)2 + y 2
vx = = U∞ + · � � − · 2(x − a)
∂x 4π (x + a)2 + y 2 � −�
(x� a)�
2
+ y2 [(x − a)2 + y 2 ] �
2

m x+a x−a
⇒ vx = U∞ + −
2π (x + a)2 + y 2 (x − a)2 + y 2
2 ��
∂φ m � −�
(x� a)� + y2 2y (x + a)2 + y 2
vy = = · � � � − · 2y
∂y 2
4π (x + a) + y 2 (x�
� −� 2
a) + y 2
[(x − a)2 + y 2 ] �
2

my 1 1
⇒ vy = −
2π (x + a)2 + y 2 (x − a)2 + y 2

∂ψ
Alternatively, one can find v by using: vx = ∂y
, vy = − ∂ψ
∂x
.

Find the stagnation point(s) by finding (x, y) such that vx = vy = 0.

vy = 0 at y = 0, ∀x

2.25 Advanced Fluid Mechanics 2 Copyright @


c 2010, MIT
Potential Flow Kundu & Cohen 6.4

Plug in y = 0 into vx and find x that lets vx = 0:


 
m x+a x−a
vx (x, y = 0) = U∞ + − =0
2π (x + a)2 (x − a)2
 
m 1 1
U∞ + − =0
2π x + a x − a
ma
OR (after some algebra. . . ) : x2 − a2 − =0
πU∞

Using the quadratic formula1 ,


m
x = ±a 1 +
aπU∞

(c) Going back to ψ:

θ1 θ2
(x, y)

 

 

m y m y
ψ = U∞ y + arctan − arctan
2π x+a 2π x − a
θ2
  θ1
m 2ay

= U∞ y − arctan a a
2π x2 + y 2 − a2



−(θ1 −θ2 )

A “Rankine oval” is defined by a curve ψ = 0, or


 
m 2ay
U∞ y − arctan =0 (1)
2π x2 + y 2 − a2

Maximum half-width, h, is obtained when x = 0:



m 2ah h
U∞ h = arctan
2π h2 − a2
 2  −a a
2πU∞ h h − a2
= arccot
m 2ah

 a 2   
2πU∞ h
⇒ h 1− = 2a cot (2)
h m

Problem Solution by Sungyon Lee, Fall 2005



1 ax2 −b± b2 −4ac
+ bx + c = 0, x = 2a

2.25 Advanced Fluid Mechanics 3 Copyright @


c 2010, MIT
Potential Flow Kundu & Cohen 6.4

Figure 1: MATLAB® plot of streamlines for a Rankine oval.

2.25 Advanced Fluid Mechanics 4 Copyright @


c 2010, MIT
MIT OpenCourseWare
http://ocw.mit.edu

2.25 Advanced Fluid Mechanics


Fall 2013

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

You might also like