Generalized Kutta-Joukowski Theorem
Generalized Kutta-Joukowski Theorem
Generalized Kutta-Joukowski Theorem
KEYWORDS
Incompressible ow;
Induced drag;
Induced lift;
Multi-airfoils;
Vortex
Abstract For purpose of easy identication of the role of free vortices on the lift and drag and for
purpose of fast or engineering evaluation of forces for each individual body, we will extend in this
paper the KuttaJoukowski (KJ) theorem to the case of inviscid ow with multiple free vortices and
multiple airfoils. The major simplication used in this paper is that each airfoil is represented by a
lumped vortex, which may hold true when the distances between vortices and bodies are large
enough. It is found that the KuttaJoukowski theorem still holds provided that the local freestream
velocity and the circulation of the bound vortex are modied by the induced velocity due to the outside vortices and airfoils. We will demonstrate how to use the present result to identify the role of
vortices on the forces according to their position, strength and rotation direction. Moreover, we will
apply the present results to a two-cylinder example of Crowdy and the Wagner example to demonstrate how to perform fast force approximation for multi-body and multi-vortex problems. The
lumped vortex assumption has the advantage of giving such kinds of approximate results which
are very easy to use. The lack of accuracy for such a fast evaluation will be compensated by a rigorous extension, with the lumped vortex assumption removed and with vortex production included,
in a forthcoming paper.
2014 Production and hosting by Elsevier Ltd. on behalf of CSAA & BUAA.
Open access under CC BY-NC-ND license.
1. Introduction
For a two-dimensional incompressible ow around a single airfoil with a sharp trailing edge at incidence, it is well known that
the KuttaJoukowski (KJ) hypothesis holds at least for steady
* Corresponding author. Tel.: +86 10 62784116; fax: +86 10
62772480.
E-mail address: ziniuwu@mail.tsinghua.edu.cn (Z. Wu).
Peer review under responsibility of Editorial Committee of CJA.
unseparated ow. In this hypothesis, viscosity is explicitly ignored but implicitly incorporated in the Kutta condition (see
Refs. 1,2). The Kutta condition imposes a circulation or a bound
vortex attached to the airfoil and creates a starting vortex, of an
opposite sign, which moves in the downstream direction.
The lift thus predicted by the KuttaJoukowski theorem
within the framework of inviscid ow theory is quite accurate
even for real viscous ow, provided the ow is steady and
unseparated.
There are a number of applications where we encounter
multiple vortices and multiple airfoils. Streitlien and Triantafyllou3 considered a single Joukowski airfoil surrounded
with point vortices convecting freely and derived a force formula. In this formula the force is related to the relative velocity
1000-9361 2014 Production and hosting by Elsevier Ltd. on behalf of CSAA & BUAA. Open access under CC BY-NC-ND license.
http://dx.doi.org/10.1016/j.cja.2013.07.022
Generalized KuttaJoukowski theorem for multi-vortex and multi-airfoil ow (a lumped vortex model)
of the point vortices. Smith et al.4,5 studied multi-blade ows
with interaction. Panel methods have also been developed
for multi-element airfoils (see Ref. 6). Crowdy7 found the complex potential to calculate the lift for a nite stack of cylinders
with imposed circulation. Force decomposition for multi-body
ows has also been considered using an integral approach (see
Refs. 8,9). None of this work considered the force formula for
each airfoil in the form of the KuttaJoukowski theorem. Wu
et al.10 extended the Lagally theorem to multi-body ows with
free point vortices, but not considering the existence of any
bound vortex.
Through nding the complex potential and using the
Blasius theorem, Katz and Plotkin6 (see chapter 6.9, generalized KuttaJoukowski theorem) developed a generalized
KuttaJoukowski theorem for an airfoil in interaction with
another airfoil represented by a lumped vortex of opposite
circulation.
For problems with multiple airfoils and multiple free vortices, it is always possible to derive an exact theory or use
numerical computation to obtain the forces on each body.
However, for easy identication of the role of each free vortex
and for fast evaluation of force approximation as required by
engineers, we require the force formulas to be explicit, capturing the main physics and easy to use, almost as simple as the
classical KuttaJoukowski theorem. For this purpose, we
make an extension of the generalized KuttaJoukowski theorem to the case of multiple airfoils with multiple free vortices,
using the lumped vortex assumption as by Katz and Plotkin.6
In this assumption each airfoil is represented by a lumped vortex at the position of the center of the airfoil. The accuracy of
this assumption will be assessed in conclusion. The extension
to multiple airfoils without the lumped vortex assumption will
be considered in a forthcoming paper. This paper will be organized as follows. In Section 2, we derive force formula in terms
of induced velocity for a single airfoil interacting with multiple
free vortices. In Section 3 we extend the force formula to the
case of multiple body in such a way that the force formula
holds individually for each airfoil. Supplemented by a two cylinder example and the Wagner problem, which are presented
in Appendix A and B respectively, we will assess the accuracy
and usefulness of the present results in Section 4. In Section 4,
we will also discuss how to use the present results to identify
the role of outside vortices on the lift and drag.
2. Single airfoil with multiple free vortices
Consider an incompressible two-dimensional ow around an
airfoil with a velocity eld v = (u, v) at constant density q, in
an unbounded domain Rf. The freestream velocity V1 is
assumed horizontal. The circulation of the bound vortex is
dened as Cb =oA(udx + vdy) for the closed curve oA along
the airfoil, with a counterclockwise path, so that a clockwise
circulation has a negative sign. The classic KuttaJoukowski
theorem expresses the lift (L) and the drag (D) per unit span
as:
L qV1 Cb
1
KJ theory
D0
Katz and Plotkin6 showed that for an airfoil in interaction
with another airfoil, represented by a lumped vortex of
opposite circulation, the force for the airfoil is given by:
35
L qV1 us Cb
Generalized KJa
D qvs Cb
where (us, vs) is the velocity at the (center of the) current airfoil
induced by an outside vortex.
Now, we consider M free vortices, each with a position
(xj, yj) and circulation Cj for j = 1, 2, , M around an airfoil
represented by a bound vortex of circulation Cb at the center of
the airfoil. A point vortex j, whether it is a free or bound vortex, induces a velocity at (x, y) given by:
8
y yj
Cj
>
v x; y
>
>
< x
2p x xj 2 y yj 2
3
>
Cj
x xj
>
>
: vy x; y
2p x xj 2 y yj 2
It is well-known that a single point vortex in relative motion, at the velocity (dxj/dt, dyj/dt), to the uid offers a force
to the body given by11:
8
dxj
>
>
L
q
V
Cj
>
j
1
<
dt
Saffman
4
>
>
dyj
>
: Dj q
Cj
dt
In the case of a bound vortex (bounded to a body) surrounded by M free vortices, the forces due to each vortex
sum up to give:
8
M
X
>
dxj
>
>
Cj
V1
> L qV1 Cb q
>
<
dt
j1
5
M
>
X
dyj
>
>
>
D
q
C
j
>
:
dt
j1
The force formulas Eq. (5) or formulas in similar forms
have been frequently used to study the interaction of a cylinder
with M free vortices (see Ref. 12 for more details and
references).
Now we relate the velocity (dxj/dt, dyj/dt) for each vortex to
the induction by other vortices to nd an explicit and new
force formula through Eq. (5).
With r2ij xj xi 2 yj yi 2 and r2jb xj xb 2
yj yb 2 , the velocity of the free vortices due to free
convection induced by the bound vortex and other free
vortices is:
8
M
X
>
Ci yj yi
>
> dxj V1 Cb yj yb
>
>
2
< dt
2prjb
2pr2ij
i1;ij
6
M
>
X
dyj Cb xj xb
>
Ci xj xi
>
>
>
: dt
2pr2jb
2pr2ij
i1;ij
Inserting the above equations into the force formula
Eq. (5), we obtain:
8
M
M X
M
X
X
Cb Cj yj yb
Cj Ci yj yi
>
>
>
q
> L qV1 Cb q
2
>
2pr
2pr2ij
<
jb
j1
j1 i1;ij
>
M
M X
M
X
X
>
Cb Cj xj xb
Cj Ci xj xi
>
>
q
>
: D q
2
2pr
2pr2
j1
jb
j1 i1;ij
ij
The last terms on the right hand sides in both forces are due
to mutual interactions between the free vortices. For each pair
36
C. Bai, Z. Wu
and vs) along this contour remains the same whether the vortices outside of this contour are bound (xed) or free.
For the problem considered here, the pressure is related to
the velocity through the unsteady Bernoulli equation
1
p /t u2 v2 C. The only difference between the
2
pressure along the contour oB for an outside FV and an
outside BV, both supposed to have a circulation C which induces a velocity potential / = Ch/(2p), comes from the term
/t. We have /t 0 for an outside FV and /t = 0 for an outside
BV. However, for a closed curve oB not enclosing this BV or
H
FV, we always have @B!0 /n dl 0 (n is the unit outward
normal and l is a vector along the tangent direction) since
the potential /, which is due to an outside vortex, has no singularity inside oB. Therefore, there is no difference in the integrated force due to this term. Moreover, this term is linear so
that the conclusions still hold when there are many other outside vortices (BV or FV).
Hence the integrated forces, along the contour oB, are the
same if some of the free vortices outside of oB are replaced
by bound ones. This means that the force formula Eq. (7) still
holds for each airfoil individually in the presence of other airfoils represented by bound vortices. Now we can write the
above force formulas for a multi-airfoil case. For the ith
airfoil,
(
Generalized KJc
Db;i qvs;i i
3. Multiple airfoils with multiple free vortices
As displayed in Fig. 1, we consider N airfoils (assumed nonrotating and xed in this paper) interacting with M free vortices of given circulation Cm with m = 1, 2, , M. The nth airfoil at location xnb ; ynb has an unperturbed circulation Cb,n,
i.e., a circulation in an oncoming uniform ow without other
airfoils, supposed to be known from traditional theories.
For a single airfoil with a bound vortex (BV) interacting
with a number of free vortices (FVs), we have proved, without
working with pressure, the force formula Eq. (7). Hence using
pressure integration along the bound vortex (namely, a contour oB of vanishing size just enclosing the bound vortex
and excluding the free vortices) also leads to the same formula.
The instantaneous ow pattern (thus the induced velocities us
where us,i = ub,i + uf,i and vs,i = vb,i + vf,i are the total
induced velocity components at the ith airfoil due to other
airfoils and free vortices. Here:
ub;i
N
X
yn yi
n;
2pr2in
n1;ni
vb;i
N
X
xn xi
n
2pr2in
n1;ni
and
uf;i
M
X
ym yb;i
m1
2pr2im
Cm ;
vf;i
M
X
xm xb;i
m1
2pr2im
Cm
10
Fig. 1
11
Generalized KuttaJoukowski theorem for multi-vortex and multi-airfoil ow (a lumped vortex model)
Under lumped vortex assumption we may assume us,i, vs,i to
be small compared to V1, hence we obtain a system of N linearized equations for solving i . In this case:
i pcAi V1 us;i ai vs;i =V1 /Bi
12
For vertically aligned airfoils without considering free vortices, we have vs,i = 0 and as,i = 0, so that i Cb;i 1 us;i =V1 ,
and then we have from Eqs. (9) and (12):
i
N
X
n1;ni
Cb;i
n Cb;i
2pV1 yn yi
13
Similarly, for horizontally aligned airfoils without considering free vortices, we have us,i = 0 and i Cb;i pcAi ms;i , and
then:
i
N
X
n1;ni
cAi
n Cb;i
2xn xi
14
37
q
1
L
C
b;2 Cb;1 5
>
< b;1
2ph
2ph
>
1
25
>
>
: Lb;2 q 1
Cb;1 Cb;2 5
2ph
2ph
A1
5
1
1
>
<
2ph
2ph 2ph
1
>
>
5
5
5
>
>
: 2 5 1
1
2ph
2ph 2ph
A2
38
1
F
b;1
>
<
2ph
2ph
2ph 2ph
A3
1
>
>
5
5
5
5
>
>
5 1
1
: Fb;2 1
2ph
2ph
2ph 2ph
The lift forces as a function of h are displayed in Fig. A1.
When h > 3, the present result (see Curve b in Fig. A1)
determined by Eq. (A1) compares well with that of Crowdy
(see Curve c in Fig. A1), in the case of xed circulations. When
h < 3, for which the cylinders can no longer be represented by
lumped vortices, the present theory is no longer valid. When
the circulations are allowed to depend on the distance according to Eq. (A3), the lift forces (see Curve a in Fig. A1) are higher than those given by the xed circulation model.
Appendix B. Wagner problem
Consider the well-known problem of Wagner13, which is the
initial lift for a at plate of length cA, set instantaneously into
motion with a constant velocity V1 and at a small angle of attack a. For this problem, an accurate and thus mathematically
difcult solution is required to consider a vortex sheet downstream to the trailing edge. Here we consider the lumped vortex approximation, by assuming a bound vortex of circulation
Cb(t) at x = xb, the center of the plate, and a starting vortex of
circulation Ca(t) = Cb(t) downstream to the airfoil at x = xs.
The velocity induced at x = xb by the starting vortex is
Cs t
us ; vs 0;
. Hence according to Eq. (12),
2pxs xb
Cs t
Cb pcA V1 a
2pV1 xs xb
Cb t
pcA V1 a
2pV1 xs xb
If at the initial time t = 0, we assume the starting vortex is just
at the trailing edge (xs(0) = cA)and the bound vortex is at the
1
center of the airfoil xb cA , and then with dxs/dt = V1
2
and thus xs = cA + V1t, the above equation can be solved to
give:
1
1
cA
1
Cb t 1
C1 1
C1
2xs xb
1 2V1 t=cA
C. Bai, Z. Wu
where C1 = pcAV1a is the circulation of the bound vortex
when the starting vortex has moved innitely downstream to
the airfoil. Then, with Eq. (7), the lift and (induced) drag are
found to be:
8
1
>
1
>
>
L
L1
>
<
1 2V1 t=cA
1
>
>
q
1
>
>
1
C21
:D
pcA 2pV1 t
1 2V1 t=cA
where L1 = qV1C1a is the nal lift after a steady state is
reached. Thus at t = 0, we have:
8
1
>
< L L1
2
>
: D p cA V2 a2
1
4
Thus even with the lumped vortex approximation we can
recover the well known results of Wagner (see Refs. 11,14) that
for an impulsively starting ow of a at plate, the initial lift is
exactly one-half of the nal lift (the Wagner effect) and the inip
tial drag is D cA V21 a2 . Hence the lumped vortex approach,
4
though approximate but quite simpler to use than the Wagner
approach based on the Laplace transform, gives reasonably
good results even for the initial behavior. We remark that
the initial lift has also been studied for airfoils with thickness
effect (see Refs. 14,15).
Appendix C. Extension of the results
By using a special momentum approach and with the help of
interchange between singularity velocity and induced ow
velocity, we have recently derived in a physical way the explicit
force formulas for two-dimensional inviscid ow involving
multiple bound and free vortices, multiple airfoils and vortex
production16, without use of lumped vortex assumption so that
the force formulas are exact for arbitrary distance between
bodies. As in the present work, these new force formulas hold
individually for each airfoil thus allowing for force decomposition and the contributions to forces from singularities (such
as bound and image vortices, sources and doublets) and bodies
out of an airfoil are related to their induced velocities at the
location of singularities inside this airfoil. The force contribution due to vortex production will related to the vortex production rate and the distance between each pair of vortices in
production. For the Crowdy example, the comparison will
be good even the distance between the two cylinders is small.
We remark that force formulas for multibody problems and
in the form of induced velocities have been also studied by
Wu et al.10 They followed Lagally theorem to obtain the forces
and their results are restricted to problems without bound vortices and without vortex production. The results of Bai et al.16,
which are obtained through a different approach, are more
general since bound vortex and vortex production are
considered.
References
1. Batchelor FRS. An introduction to uid dynamics. Cambridge: Cambridge University Press; 1967.
Generalized KuttaJoukowski theorem for multi-vortex and multi-airfoil ow (a lumped vortex model)
2. Grighton DG. The Kutta condition in unsteady ow. Annu Rev
Fluid Mech 1985;17(1):41145.
3. Streitlien K, Triantafyllou MS. Force and moment on a Joukowski prole in the presence of point vortices. AIAA J
1995;33(4):60310.
4. Smith FT, Timoshin SN. Planar ows past thin multi-blade
congurations. J Fluid Mech 1996;324:35577.
5. Bowles RGA, Smith FT. Lifting multi-blade ows with interaction. J Fluid Mech 2000;415:20326.
6. Katz J, Plotkin A. Low speed aerodynamics. Cambridge: Cambridge University Press; 2001.
7. Crowdy D. Calculating the lift on a nite stack of cylindrical
aerofoils. Proc R Soc A 2006;462(2069):1387407.
8. Chang CC, Yang SH, Chu CC. A many-body force decomposition
with applications to ow about bluff bodies. J Fluid Mech
2008;600:95104.
9. Hsieh CT, Kung CF, Chang CC. Unsteady aerodynamics of
dragony using a simple wing-wing model from the perspective of
a force decomposition. J Fluid Mech 2010;663:23352.
10. Wu CT, Yang FL, Young DL. Generalized two-dimensional
Lagally theorem with free vortices and its application to uidbody interaction problems. J Fluid Mech 2012;698:7392.
11. Saffman PG. Vortex dynamics. New York: Cambridge University
Press; 1992.
39
Bai Chenyuan is a Ph.D. student of Fluid Mechanics and Aerodynamics in the School of Aerospace Engineering at Tsinghua University. Her main interest is vortex dynamics for aerodynamic.
Wu Ziniu is a professor of Fluid Mechanics in the School of Aerospace
Engineering at Tsinghua University. His main interests are vortex
dynamics, shock dynamics. He with his students obtained a number of
analytical results for vortex ow and shock ow, with 10 papers
published in Journal of Fluid Mechanics during the years 20062013.