next up previous
Next: Zhukovskii's Hypothesis Up: Incompressible Aerodynamics Previous: Theorem of Kutta and


Cylindrical Airfoils

For the moment, let us work in the complex $ \zeta $ -plane, where $ \zeta=\xi+{\rm i}\,\eta$ . Consider a cylindrical airfoil with a circular cross-section of radius $ a$ , centered on the origin, that is situated in a uniform, high Reynolds number wind of speed $ V$ whose direction subtends a (clockwise) angle $ \alpha $ with the negative $ \xi$ -axis. Let $ {\mit\Gamma}$ be the circulation of air around the airfoil. A slight generalization of the analysis of Section 6.4 reveals that the appropriate complex velocity potential is

$\displaystyle F(\zeta) = V\left(\zeta\,{\rm e}^{\,{\rm i}\,\alpha} + \frac{a^{\...
...}\right) + {\rm i}\,\frac{{\mit\Gamma}}{2\pi}\,\ln\left(\frac{\zeta}{a}\right),$ (9.14)

whereas the associated stream function takes the form

$\displaystyle \psi(r,\theta) = V\left(r-\frac{a^{\,2}}{r}\right)\sin(\theta+\alpha) + \frac{{\mit\Gamma}}{2\pi}\,\ln\left(\frac{r}{a}\right),$ (9.15)

where $ \zeta=r\,{\rm e}^{\,{\rm i}\,\theta}$ . It follows that

$\displaystyle \frac{dF}{d\zeta} = V\,{\rm e}^{\,{\rm i}\,\alpha} + {\rm i}\,\fr...
...}}{2\pi\,\zeta} - \frac{V\,a^{\,2}\,{\rm e}^{\,-{\rm i}\,\alpha}}{\zeta^{\,2}}.$ (9.16)

Comparison with Equation (9.6) (with $ z=\zeta$ ) reveals that

$\displaystyle B = - V\,a^{\,2}\,{\rm e}^{-{\rm i}\,\alpha}.$ (9.17)

Hence, Equations (9.10) and (9.13) yield

$\displaystyle {\bf V}$ $\displaystyle = V\left(-\cos\alpha\,{\bf e}_\xi + \sin\alpha\,{\bf e}_\eta\right),$ (9.18)
$\displaystyle {\bf L}$ $\displaystyle = \rho\,V\,{\mit\Gamma}\left(\sin\alpha\,{\bf e}_\xi+\cos\alpha\,{\bf e}_\eta\right),$ (9.19)
$\displaystyle M$ $\displaystyle = 0,$ (9.20)

where $ {\bf V}$ is the wind vector, $ {\bf L}$ the lift vector, and $ M$ the moment of the lift vector about the origin. We conclude that, for a cylindrical airfoil of circular cross-section, the lift vector is normal to the wind vector, and the line of action of the lift passes through the centroid of the cross-section (because the lift generates zero moment about the origin). (See Figure 9.1.)

Figure 9.1: A cylindrical airfoil of circular cross-section.
\begin{figure}
\epsfysize =3in
\centerline{\epsffile{Chapter09/air.eps}}
\end{figure}

Of course, a cylindrical airfoil of circular cross-section is completely unrealistic, because its back side (i.e., the side opposite to that from which the wind is incident) is not sufficiently streamlined to prevent boundary layer separation. (See Chapter 8.) However, as described in Section 6.7, we can use the conformal transformation

$\displaystyle z = \zeta+ \frac{l^{\,2}}{\zeta}$ (9.21)

to transform a cylinder of circular cross-section in the $ \zeta $ -plane to a cylinder of elliptical cross-section in the $ z$ -plane. (Note that both cross-sections have centroids located at the origin.) Moreover, a cylindrical airfoil of elliptical cross-section that is sufficiently elongated, and whose major axis subtends a sufficiently small angle with the incident wind direction, constitutes a realistic airfoil, because its back side is, for the most part, closely aligned with the external flow.

Figure 9.2: A cylindrical airfoil of elliptical cross-section.
\begin{figure}
\epsfysize =3in
\centerline{\epsffile{Chapter09/ellipse.eps}}
\end{figure}

An elliptical airfoil of width $ c$ and thickness $ \delta<c$ , as shown in Figure 9.2, is obtained when the parameters $ a$ and $ l$ are given the following values:

$\displaystyle a$ $\displaystyle = \frac{1}{4}\,(c+\delta),$ (9.22)
$\displaystyle l$ $\displaystyle =\frac{1}{4}\,(c^{\,2}-\delta^{\,2})^{1/2}.$ (9.23)

In this case, the surface of the airfoil satisfies the parametric equations

$\displaystyle x$ $\displaystyle = \frac{c}{2}\,\cos\theta,$ (9.24)
$\displaystyle y$ $\displaystyle = \frac{\delta}{2}\,\sin\theta.$ (9.25)

In particular, the airfoil's leading and trailing edges correspond to $ \theta=0$ and $ \theta=\pi$ , respectively.

According to Equations (9.16) and (9.21), the complex velocity in the $ z$ -plane is given by

$\displaystyle \frac{dF}{dz} = \frac{dF}{d\zeta}\,\frac{d\zeta}{dz} = \left[V\,{...
...lpha}}{\zeta^{\,2}}\right]\left(\frac{\zeta^{\,2}}{\zeta^{\,2}-l^{\,2}}\right).$ (9.26)

Thus, on the airfoil surface, where $ \zeta=a\,{\rm e}^{\,{\rm i}\,\theta}$ , we obtain

$\displaystyle \frac{dF}{dz} = {\rm i}\left[V\,\sin(\alpha+\theta) + \frac{{\mit...
...ight]\frac{(c+\delta)}{\left(\delta\,\cos\theta+{\rm i}\,c\,\sin\theta\right)}.$ (9.27)

A long way from the airfoil, $ \zeta\simeq z-l^{\,2}/z$ , so that Equation (9.26) reduces to

$\displaystyle \frac{dF}{dz}\simeq V\,{\rm e}^{\,{\rm i}\,\alpha} + {\rm i}\,\fr...
...2}\,{\rm e}^{\,{\rm i}\,\alpha}-a^{\,2}\,{\rm e}^{-{\rm i}\,\alpha})}{z^{\,2}}.$ (9.28)

A comparison with Equation (9.6) reveals that the circulation of air around the airfoil takes the same value, $ {\mit\Gamma}$ , in both the complex $ \zeta $ - and $ z$ -planes. In other words, the conformal transformation (9.21) does not modify the circulation. The transformation also does not modify the external wind speed or direction [because, from (9.16) and (9.28), $ dF/d\zeta= dF/dz=V\,{\rm e}^{\,{\rm i}\,\alpha}$ at very large $ \vert\zeta\vert$ and $ \vert z\vert$ ]. On the other hand, it is clear that the constant $ B$ , which takes the value zero in the complex $ \zeta $ -plane, takes the value

$\displaystyle B = V\,(l^{\,2}\,{\rm e}^{\,{\rm i}\,\alpha} -a^{\,2}\,{\rm e}^{-{\rm i}\,\alpha})$ (9.29)

in the complex $ z$ -plane. Hence, Equations (9.10) and (9.13) reveal that

$\displaystyle {\bf V}$ $\displaystyle = V\,{\bf e}_\parallel,$ (9.30)
$\displaystyle {\bf L}$ $\displaystyle = \rho\,V\,{\mit\Gamma}\,{\bf e}_\perp,$ (9.31)
$\displaystyle M$ $\displaystyle =\frac{\pi}{8}\,\rho\,V^{\,2}\,(c^{\,2}-\delta^{\,2})\,\sin(2\,\alpha),$ (9.32)

where $ {\bf V}$ is the wind vector, $ {\bf L}$ the lift vector, and $ M$ the moment of the lift vector about the origin. Here,

$\displaystyle {\bf e}_\parallel = -\cos\alpha\,{\bf e}_x + \sin\alpha\,{\bf e}_y$ (9.33)

is a unit vector parallel to the incident wind direction, and

$\displaystyle {\bf e}_\perp= {\bf e}_\parallel \times{\bf e}_z= \sin\alpha\,{\bf e}_x+\cos\alpha\,{\bf e}_y$ (9.34)

is a unit vector perpendicular to the wind direction. We conclude that, for a cylindrical airfoil of elliptic cross-section, the lift vector is normal to the wind vector, but the line of action of the lift intersects the major axis of the airfoil a distance

$\displaystyle d = \frac{M}{L\,\cos\alpha} = \frac{1}{4}\,(c-\delta)\,\frac{\sin\alpha}{\widehat{\mit\Gamma}}$ (9.35)

in front of the cross-section's centroid, $ C$ , where $ \widehat{\mit\Gamma}= {\mit\Gamma}/[\pi\,V\,(c+\delta)]$ . (See Figure 9.2.) Incidentally, the point, $ F$ , at which the line of action of the lift intersects the airfoil's major axis is conventionally termed the focus of the airfoil.


next up previous
Next: Zhukovskii's Hypothesis Up: Incompressible Aerodynamics Previous: Theorem of Kutta and
Richard Fitzpatrick 2016-01-22