Flow Around a Circular Cylinder

Flow around a circular cylinder can be approached from the previous example by bringing the source and the sink closer. Then we are considering a uniform flow in combination with a doublet. The stream function and the velocity potential for this flow are given by,

Figure 4.28: Schematic for Flow past a Circular Cylinder

$\displaystyle \psi $ $\displaystyle = U_\infty r \sin \theta - {{K \sin \theta} \over r}$ (4.103)
$\displaystyle \phi $ $\displaystyle = U_\infty r \cos \theta + {{K \cos \theta} \over r}$ (4.104)

Streamlines for this flow are plotted in Fig. 4.29.

 


Figure 4.29
: Flow past a Circular Cylinder

Figure 4.30: Stagnation Points for Flow about a Circular Cylinder

The velocity components are given by,

$\displaystyle v_r $ $\displaystyle = {1 \over r} {{\partial \psi} \over \partial \theta} = \cos \theta \left(U_\infty - {K \over r^2}\right)$ (4.105)
$\displaystyle v_\theta $ $\displaystyle = -{{\partial \psi} \over { \partial r}} = -\sin \theta
 \left( U_\infty + {K \over r^2} \right)$ (4.106)

It is seen that the radial velocity is zero when

If we recognise this particular streamline as the surface of the circular cylinder then the radius of the cylinder a is given by,


The equations for the streamline, velocity potential and the velocity components are replaced by,

$\displaystyle \psi $ $\displaystyle = U_\infty r \left( 1 -  {a^2 \over r^2} \right)\sin \theta$ (4.109)
$\displaystyle \phi $ $\displaystyle = U_\infty r \left( 1 +  {a^2 \over r^2} \right)\cos \theta$ (4.110)
$\displaystyle v_r $ $\displaystyle = U_\infty \left( 1 -  {a^2 \over r^2} \right)\cos \theta$ (4.111)
$\displaystyle v_\theta $ $\displaystyle = -U_\infty \left( 1 +  {a^2 \over r^2} \right)\sin
 \theta$ (4.112)

The velocity components on the surface of the cylinder are obtained by putting r = a in the above expressions. Accordingly,


 

$ \sin \theta$ has a zero at 0 and 1800 and a maximum of 1 at = 900 and 2700. The former set denotes the stagnation points of the flow and the later one denotes the points of maximum surface velocity (of magnitude $ 2U_\infty$). Thus the velocity decreases from a value of $ 2U_\infty$ at $ \theta$ equals 900 to $ U_\infty$ as one moves away in a normal direction s shown in Fig 4.30.

The surface pressure distribution is calculated from Bernoulli equation. If we denote the free stream speed and pressure as $ U_\infty$ and $ p_\infty$ we have

Substituting for $ v_{\theta s} = -2 U_\infty \sin \theta$, we have

We can also express pressure in terms of pressure coefficient, Cp,

$\displaystyle C_p =1 - {\left( {v_s \over U_\infty}\right)^2}$ (4.116)
leading to
   
$\displaystyle C_p = 1 - 4 \sin^2 \theta$ (4.117)

Fig. 4.31 shows Cp plotted as a function of $ \theta$. A symmetry about y -axis is apparent. When compared to the experimentally observed Cp
distribution we see that there is some agreement in the region between$ \theta$= 00 and $ \theta$= 900 . But any agreement is lost in the other regions. The reasons for this are obvious. Viscous forces dominate the flow in the region to the right of the centreline giving rise to separation. The pressure tends to plateau out in a separated region, the level depending on whether it is a laminar separation or a turbulent one.

Figure 4.31A: Cp distribution for flow past a circular cylinder.

Figure 4.31B: Cp distribution for flow past a circular cylinder plotted around the cylinder.

Symmetry in the theoretical Cp distribution about both y-axis and x-axis shows that drag and lift forces about the cylinder are each zero. This
may also be proved by integrating pressure around the cylinder, thus,

Drag, $\displaystyle Drag, D = -\int_0 ^{2 \pi} p_s  \cos \theta  a  d \theta $ (4.118)
Lift, $\displaystyle Lift, L = -\int_0 ^{2 \pi} p_s  \sin \theta a d \theta $ (4.119)

By substituting for the surface pressure, ps from Eqn.4.115 we find,

$\displaystyle D$ $\displaystyle = - \int_0 ^ {2 \pi} p_\infty a \cos \theta d \theta -
 \frac{1}{...
...0 ^ {2 \pi} \left( \cos \theta
 - 4 \sin^2 \theta \cos \theta \right ) d \theta$ (4.120)
  $\displaystyle =  -p_\infty a \left [\sin \theta \right ]_0 ^{2
 \pi}-\frac{1}{...
...1}{2} \rho U_\infty
 ^2 a \left [ \frac{4}{3} \sin^3 \theta \right ]_0 ^{2 \pi}$ (4.121)
  = -0 -0 + 0 (4.122)
$\displaystyle L$ $\displaystyle = - \int_0 ^ {2 \pi} p_\infty a \sin \theta d \theta -
 \frac{1}{...
...y ^2 a \int_0 ^ {2 \pi} \left( \sin \theta
 - 4 \sin^3 \theta \right ) d \theta$ (4.123)
  $\displaystyle =  -p_\infty a \left [ \cos \theta \right ]_0 ^{2
 \pi}-\frac{1}...
...fty ^2 a \left [ \frac{4}{3} \cos^3
 \theta - 4 \cos \theta \right ]_0 ^{2 \pi}$ (4.124)
  = -0 -0 + 0 (4.125)

What we have just calculated is in contrast to the experimental results which do predict a significant drag for the flow about a circular cylinder. This seems to have caused in what is called D'Alembert's Paradox in honour of Jean le Rond D'Alembert (1717-1783). Now it is no more a paradox. As we discussed above we calculate a zero drag because we have not taken viscosity into account.

(c) Aerospace, Mechanical & Mechatronic Engg. 2005
University of Sydney