ries 02

Flight Mechanics 03

Optimum Cruise Performance of

Subsonic Transport Aircraft

E. Torenbeek

Optimum Cruise Performance of

Subsonic Transport Aircraft

Bibl iotheeL TU Delft

IIIII 1111

C 3021896

~I

Series 02: Flight Mechanics

03

Optimum Cruise Performance of

Subsonic Transport Aircraft
E. Torenbeek

Delft University Press / 1998

2392
330
2

Published and distributed by:

Delft University Press
Mekelweg 4
2628 CD Delft
The Netherlands
Telephone +31 (0)152783254
Fax +31 (0)152781661
e-mail: [email protected]

by order of:
Faculty of Aerospace Engineering
Delft University of Technology
Kluyverweg 1
P.O. Box 5058
2600 GB Delft
The Netherlands
Telephone +31 (0)152781455
Fax + 31 (0) 15 278 18 22
e-mail: [email protected]
website: http://www.lr.tudelft.nl/

Cover: Aerospace Design Studio, 66 .5 x 45.5 cm, by:

Fer Hakkaart, Dullenbakkersteeg 3, 2312 HP Leiden, The Netherlands
Tel. + 31 (0)71 51267 25

90-407-1579-3
Copyright

1998 by Faculty of Aerospace Engineering

All rights reserved.

No part of the material protected by this copyright notice may be
reproduced or utilized in any form or by any means, electronic or
mechanical, including photocopying, recording or by any information storage
and retrieval system, without written permission from the publisher: Delft
University Press.
Printed in The Netherlands

Contents
1 Cruise Performance for Preliminary Design
1.1 Types of Performance Problems
1.2 Optimum Cruise Performance
1.3 Scope of this Report ..

2 The Ratio of Lift to Drag

2.1 Partial and Unconstrained Optima
2.2 Constraints on Altitude or Thrust .

5
9

1
2
3

13

3 The Parameter MLjD

3.1 Maximum ML/D at Subcritical Speeds
3.2 ML/D with Compressibility Effects

15

4 The Specific Range

4.1 Generalized Engine Performance .
4.2 Maximum range parameter. . . .
4.3 Interpretation of the Derivative 1]M
4.4 Refinements in the Optima.

21

5 Analysis of Cruising Flight

5.1 The Generalized Range Equation
5.2 Range at Subcritical Speeds
5.3 Range at High Speeds ..

34

16
18
22
23
28
31
34
35
38

Prediction of the FUel Load

41

6.1 Mission Fuel . . . . . .

6.2 Reserve and Total Fuel

44
46

7 Summary of Results

49

A Statistical Derivation of the Range Parameter

53

Bibliography

58

NOMENCLATURE
a

CD
CD,Q

CD
CDL
CDM
CL
CL

Cp
CT
D
F

II, h,h, 14
g
H
h
he
K
kR
kres

L
M
m

N
n

p
Pbr

""

R
R
RH

Rh
S

sn
T

t
V
Vj

W
y

speed of sound
drag coefficient
zero-lift drag coefficient
minimum drag coefficient
logarithmic derivative of CD w.r.t .CL
logarithmic derivative of CD w.r.t . M
lift coefficient
lift coefficient for minimum drag coefficient
Power Specific Fuel Consumption
Trust Specific Fuel Consumption
drag
fuel weight flow per unit time
generalized engine performance functions
acceleration due to gravity
calorific value of jet engine fuel
altitude
energy height
induced drag factor
factor accounting for flight schedule
reserve fuel fraction
lift
Mach number
mass flow per unit time
engine RPM
exponent of M in approximation for TSFC
Range Parameter
engine shaft power
atmospheric pressure
dynamic pressure
range, distance flown
equivalent all-out range in cruising flight
range-equivalence of fuel calorific value, RH
harmonie or nominal range
reference wing area
Specific Range
net T hrust
time
true airspeed
mean exhaust jet velocity
weight (no index: All Up Weight)
CL/CL,md

(mis)
(- )
(-)
(-)
(- )
(- )

(- )
(- )
((Nl s)IW)
((NI s)1 N)
(N)
(Nis)
(- )
(mls 2)
(Jlkg)
(m)
(m)
(- )
(- )
(- )

(N)
(-)
(kg/s)
(- )

= Hl 9

(- )
(- )
(W)
(Nlm 2)
(N/m2)
(m)
(m)
(m)
(m)
(m 2)
(miN)
(N)
(s)
(m/s)
(mis)
(N)
(- )

,=

TJprop
TJM
TJT
()

P
<I>

ratio of specific heats of air,

1.40
increment
relative ambient pressure
fuel weight fraction
overall powerplant efficiency
propeller efficiency
logarithmic derivative of TJ w.r.t. Mach number
logarithmic derivative of TJ w.r .t. T/
relative ambient temperature
atmospheric density
normalized slope of payload range diagram

Indices
eb
cl
etg

er
dd
div

eq
F

f
h , es
h,Z
h, M

hoZd

z
land
m

max

md
opt
p
pr

prop
res

sZ
th
to

combustion
climb
contingency
crUIse
drag divergence
diversion
equivalent
fuel
end of cruising
horizontal flight , constant engine setting
horizont al cruise, constant lift coefficient
horizont al flight, constant Mach number
holding
initial cruising
landing
mISSlOn
maXImum
minimum drag for constant M
optimum
payload
propulsion
propeller
reserve fuel
Sea Level, ISA
thermal
take-off

(- )
(- )
(- )
(- )
(- )
(- )
(- )
(- )
(- )
(kg/m 3 )
(-)

Abbreviations
ATC
AUW
BSFC
DOC
EAS
EGT
ESDU
GA
HSC
ISA
LRC
LW
MTOW
MZFW
OEW
RPM
SL
SR
SST
TC
TIT
TSFC
TOW
UL
ZFW

Air Traffie Control

All Up Weight
Brake Specific Fuel Consumption
Direct Operating Cost
Equivalent Air Speed
Exhaust Gas Temperature
Engineering Seiences Data Unit
General Aviation
High Speed Cruise
International Standard Atmosphere
Long Range Cruise
Landing Weight
Maximum Take Off Weight
Maximum Zero Fuel Weight
Operating Empty Weight
Revolutions Per Minute
Sea Level
Specific Range
Super-Sonie Transport
Transport Category
Turbine Inlet Temperature
Thrust Specific Fuel Consumption
Take Off Weight
Useful Load
Zero Fuel Weight

Chapter 1
Cruise Performance for Preliminary
Design
1.1

Types of Performance Problems

During the design process of Transport Category (TC) and General Aviation (GA) aircraft
the designer(s) wiU inevitably be confronted with at least one of the following questions:
What are, for a given aircraftjengine combination, the best cruise altitude and
speed from the point of view of minimum fuel consumed per unit time elapsed or
per di stance travelled ?
For a specified flight mission, how to estimate the (minimum) amount of fuel and
time required? In addition to cruise fuel and time this problem involves an estimation of the fuel quantities required for the ot her phases of the flight, in particular
climb and descent, and reserve fuel estimation.
How to find the best operational flight profile of a transport aircraft , resulting in
minimum fuel consumed, elapsed time or Direct Operating Costs (DOC) with or
without constraints on range or time travelled .
From the flight mechanical point of view all three questions are closely related to the
generalized problem of trajectory optimization. The first problem is the most basic and
is usually treated as an instantanuous performance! problem, assuming equilibrium of
forces in horizontal flight: Lift = Weight and Thrust = Drag. The second problem is
more complicated since it involves integration of the instantanuous performance into path
performance. Moreover , the Take Off Weight (TOW) is an input to the computation of the
fuelload, but it is also an outcome of this calculation in case the Zero Fuel Weight (ZFW)
is specified. Practical numeri cal procedures exist to solve this problem-for example the
Newton-Raphson approach-but for preliminary design it would be nice to avail of a
1 Also

referred to as "point performance" .

closed-form solution. The cruise sector is usually treated as a (quasi) steady fiight, with
either constant or very gradual variations of the altitude and/or the speed.
The third problem is by far the most complicated from the mathematical and numerical points of view, since the fiight mechanics involved have to be treated as a dynamic
problem and the optimization requires the calculus of variations or optimal control theory.
Applications are mostly found in the operational use, with the objective of achieving fuel
reductions or devising optimum Air Traffic Control (ATC) procedures. Since fuel savings
obtained appear to he very modest-at least for medium and long range fiights- it is
not common practice to use this compute-intensive technique in the preliminary design
stage. The optimal trajectory analysis forms a subject in itself, which is not treated in
the present report; for a review the reader is referred to Visser [29] .

1.2

Optimum Cruise Performance

This Report treats the computation of range and fuel required for high speed commercial
aircraft. Cruise performance of propeller aircraft is a classical problem which is adequately
covered in the literature, hence it will not be treated explicitly, although the generalized
method presented applies to propeller aircraft as weIl. Moreover the optimization of
endurance or fiight time is excluded from the analysis; the emphasis is on maximum
range or minimum fuel. Flight conditions for minimum Direct Operating Costs (DOC)
will not be considered primarily because the results are quite specific to the operational
environment. Two more subjects are excluded from the analysis:
the effects of wind on cruise performance are disregarded; for a concise treatment
see Hale (1979) [10];
Super-Sonic Transports (SST's) are not considered explicitly, although the approach
developed is applicable to this aircraft category, provided one is aware of the existance of dual optima at subsonic and supersonic speeds.
This leaves us with the cruise performance problem of high subsonic jet aircraft which fiy
most effectively at high altitudes, unless ATC requirements or non-standard conditionsfor example engine failure--impose a limit on speed or altitude.
The question is justified as to whether the problem of cruise performance needs to be
considered any more after so many years of successful operation of commercial aircraft.
The subject is also treated extensively in several modern texthooks; see for example
Shevell [24], Mair and Birdsall [12] and Ruijgrok [21]. However, there remain several
misunderstandings and questions to be answered about range performance as illustrated
by the following brief review of the literature.
The first analyses of range for turbojet powered aircraft were published soon aft er
WWII, in the UK by Page [18] and Edwards [8], and in the US by Jonas [11], Ashkenas
[4], and Perkins and Hage [20]. The range equations were based either on a cruise/climb
technique with constant speed and lift coefficient-resulting in the Brguet equation for jet

propulsion-or on horizontal flight with constant lift coefficient, resulting in the "squareroot" equation [ll]. Although the cruise/climb technique results in longer ranges than
horizont al cruising fiight , it is usually not a favoured operational procedure due to ATCimposed constraints. But the alternative cruise technique leads to decreasing airspeed and
thrust settings as fuel is consumed during the fiight-which is not practical either-and
inferior range performance. Alternative procedures were developed later [19] for horizont al
cruising flight at constant speed, Mach number or engine rating. They are summarized
in various ESDU Data Items, [2] and [3].
There is agreement in the literature as to the fact that-for specified initial altitude
and speed-a cruise/climb fiight program yields the longe st range. There is, however, no
consensus about the definition of the best initialor mean altitude and speed. In particular
the ratio of optimum fiight speed to the Minimum Drag Speed has been debated by
several authors, even in recent publications. For example, Hale (1976)[9] points out that
different values have been derived for this ratio-namely 21/4 and 31/ 4, dependent on cruise
technique-and concludes from analysis that the lat ter value is the correct one. Ironically,
his derivation presupposes an implicit altitude constraint which inevitably leads to the
following question: which altitude is the best if the constraint is deleted? Bert [6] derived
a "new range equation" , in fact a modification of the "arctan range equation" for flight at
constant altitude and speed that can be found already in Edwards' work [8]. He computed
in an example an optimum speed equal to 1.864 times the speed for minimum drag-an
unusual result, caused by a continuously decreasing Thrust Specific Fuel Consumption
(TSFC) with increasing speed. In a recent article Miller [13] concludes that the best
altitude is one where the drag is a global minimum, "which disproves the theory that
the cruise lift/drag ratio is ,;3/2 times the maximum lift/drag ratio" , corresponding to a
fiight speed equal to 31/ 4 times the Minimum Drag Speed.

1.3

Scope of this Report

Most of the conflicting results from the literature must be ascribed to the fact that below
the drag rise there is no unconstrained optimum fiight condition for maximum range, as
will be shown by the present analysis. Early publications-for example by Page [18]point out by means of examples that optimum cruise Mach numbers might occur in the
drag rise. To the best knowledge of the author the first derivation of an analytical criterion
for optimum flight Mach number has been given by Backhaus [5]. Miele [17] wrote a
thourough treatment of optimum cruise performance at transonic speeds, but practical
applications of his rather complex theory have not been found in the literature. The
generalized criteria for maximum Specific Range derived by Torenbeek and Wittenberg
[27]-which are valid at high subsonic speeds and for arbitrary propulsion systems-have
been referred to as "fairly sophisticated" by Martinez-Val et al [15]. Since they provide
the basis for a comprehensive theoretical fundament, these criteria will be readdressed
in the present Report and further augmented to derive range performance and stage fuel
required.

In the first chapters the effects of compressibility on the liftjdrag ratio and on the
specific range will be treated. Conditions will be derived for the best fl.ight Mach number
and cruise altitude resulting in maximum range of aircraft with various types of gas turbine
based propulsion systems: turboprops, turbojets and turbofans. Analytical conditions will
be derived for the unconstrained optimum lift coefficient and Mach number, as well as for
the case of constraints on the altitude or the engine thrust. The second part deals with the
problem of making a good estimate of the fuel required during the cruising fl.ight and of
the total fuel required to fl.y a specified mission, including reserves. The complication that
fl.ight takes place at transonic speeds makes it necessary to reconsider the fl.ight techniques
for optimum range performance and minimum stage fuel required. Finally a useful method
will be selected on the basis of accuracy, simplicity and suitability to preliminary aircraft
design optimization studies, which require knowledge of the sensitivity of the fuelload to
variations in the aircraft design characteristics.

Chapter 2
The Ratio of Lift to Drag
The lift/drag ratio-sometimes referred to as the aerodynamic fineness ratio or aerodynamic efficiency-is an important performance parameter, since it determines the drag
and thrust required in horizont al flight ,
W

T= D

= L/D

(2.1)

For given Specific Fuel Consumption of the powerplant, the fuel flow rate is therefore
inversely proportional to the lift/drag ratio, which is determining to a large extent:
the endurance and range for a given amount of fuel,
the amount of fuel required to fly during a specified time period or to cover a
specified distance.
The aircraft drag coefficient of high speed transport aircraft,

D
qS

CD = is a function of the lift coefficient

L
CL = qS

(2.2)

(2.3)

and the flight Mach number

M= V/a

(2.4)

where the dynamic pressure is defined as

q = 1/2pV 2 = 1/2,pM2

(2.5)

Most aircraft with low subsonic speeds feature a single drag polar for the en route configuration, independent of the Mach number:

(2.6)

Although variations of Reynolds number and center of gravity location theoretically have
an effect on the drag polar, this is usually ignored in preliminary design and drag polars
are refered to mean conditions of loading and typical flight profiles. Some aircraft feature a
slight variation of the drag coefficient with Mach number at low subsonic speeds, mainly
due to variation in Reynolds number. Compressibility effects come into play at high
subsonic speeds, and the drag coefficient has to be expressed in terms of the lift coefficient
and Mach number

(2.7)
A representative set of drag polars is depicted on Figure 2.1a; an alternative representation
of the same data is presented in Figure 2.1 b. For the example aircraft there is a gradual
drag reduction with Mach number at low speeds, which may be caused by an increasing
Reynolds number andjor decreasing skin friction associated with the heating up of the
boundary layer. For low lift coefficients the drag begins to increase first in the form of a
"drag creep" for Mach numbers in excess of about 0.65, while for lift coefficients in excess
of 0.50 there is also a Mach number effect at low speeds. This drag increase is probably
caused by compressibility effects at locations where pressure peaks occur, e.g. at the nose
of the airfoils. For Mach numbers above :::::J 0.80 the drag rise is much faster: the so-called
"drag divergence". The usual definition of the drag divergence Mach number Mdd is based
on the slope of the drag curve for constant lift coefficient,

8CD
8M

(2.8)

= 0.10

This condition defines a Mach number which depends on the lift coefficient and therefore
the drag initiation is represented in the figure by a curve. By using these drag polars,
values of the liftjdrag ratio have been derived-see Figure 2.2a-and contours of constant
Lj D ratios have been plotted on Figure 2.2b. These contour plots have the advantage
that constant L j D curves do not intersect since each combination of the independent
variables CL and M define only one drag coefficient. This drag representation is therefore
completely unambiguous.
The thrust required in horizontal steady flight follows directly from Equation 2.1

= D = CLjCD

(2.9)

The condition for maximum liftjdrag for given Mach number-resulting in minimumdrag
(md) and hence minimum thrust-is defined by the lift coefficient CL,md .
For low subsonic aircraft designed to cruise at speeds up to about M ~ 0.6, the lift
coefficient CL,md has a unique value defining the Minimum Drag Speed, which is obtained
from Equation 2.3, with L = W:
2WjS

pCL,md

(2.10)

0.10
M

0.09
[0

0.08

0.01

0.06

0.05

~
0.92

I
0.45

20

a. CL/CD versus CL

15

0.45

0.65

Hach numbor

\
\

\
\

(DL =1
0.4

0.5

0.6

0.7

0.8

0.9

1.0

(L

1.0

0.9
(L

b. L/ Din the CL versus M plane

12

13
14

0.8

15
16

0.7
17
0.6

18

LID
0.5

----

_ (DL =1

(L.md

-----------------

0.45
0.40

18
0.4

0.3

0.2

0.1

0.35

17

0.30

16
15
14
13
12

05

0.6

0.7

0.9

0.8

Hach numbor

Figure 2.2: Different representations of the lift to drag ratio

or the Mach number 1 :

2WjS

(2.11)

ipCL,md

In this case the Minimum Drag Speed is thus a constant Equivalent Air Speed (EAS),
provided the All Up Weight (AUW) is specified.

For high subsonic aircraft CL,md is a function of the Mach number; the example in
Figure 2.2b indicates that it is nearly constant for Mach numbers up to about 0.65. Above
this speed the effects of compressibility on drag re sult in rapidly decreasing values of the
Lj D ratio. The contours of constant Lj D values are almost straight at low speeds but
become closed at high Mach numbers, and a dip occurs in CL,md at high Mach numbers.

2.1

Partial and Unconstrained Optima

The mathematical condition of maximum Lj D is-for the case of general drag polars
according to Equation 2.7-obtained from logarithmic differentiation:
dlog(CL/CD )

= dlogCL -

dlogCD

=0

(2.12)

with
dlog CD

o log CD
o log CL dlog CL

o log CD

+ o log M

dlog M

(2.13)

The following logarithmic derivatives are now introduced,

CD ~ o log CD
L
o log CL

CLOCD
CD oCL

~ o log CD

---

CD
M

ologM

MOCD

CD oM

(constant M)

(2.14)

(constant CL)

(2.15)

From their definitions it is clear that logarithmic derivatives can be interpreted as a

percentage change in a dependent variabIe due to a given percentage change of the independent variabIe. Log-derivatives are therefore non-dimensional; their numeri cal value
is often between zero and plus or minus two or three, which is convenient for numeri cal
treatment.
After substitution of the log-derivatives Equation 2.1 3 reads as follows,
(2.16)
1 In the case of a constraint on the thrust the Mach number for minimum drag is obtained from the
equilibrium Thrust=Drag.

Equation 2.12 may now be written alternatively by substitution of this result:

(2.17)
The partial optima for CL or Mare obtained by setting dCL=O or dM=O-alternatively
dlog CL=O or dlog M = O-resulting in
constant M: CD L = 1 -+

(2.18)

D
constant CL : CDM = 0 -+ oC = 0
oM

(2.19)

The first of these partial optima can be allocated in the drag polar(s)-see Figure 2.3aby means of a tangent from the origin, the second is relevant only for high speed aircraft
where CDM = 0 defines a relationship between CL and M as illustrated by Figure 2.2b.
However, in most cases the drag at subcritical speeds is defined as a unique polar-hence
CDM = O-up to a certain Mach number. Above this speed CDM > 0 and the liftjdrag
ratio is generally degraded.
The condition for an unconstrained maximum value of Lj Dis that both partial optima
occur at the same time. Hence the combination of CDL = 1 and CD M = 0 defines a unique
condition for (Lj D)max only in the case that a maximum or a minimum exists for a given
lift coefficient. Figure 2.2b this condition is identified as Point A, with Lj D = 18.40. H,
on the ot her hand, for low subsonic speeds there is only one drag polar, the aircraft has
a constant maximum liftjdrag ratio defined by CL = CL.md .

For the case of parabolic drag polars, generally with coefficients dependent on the
Mach number,
CD = CD.o(M)

+ K(M)CI

with K(M) ~f dCDjd(Cl)

(2.20)

the conditions for minimum drag for given Mach number are

CL md

= JCD.o(M)jK(M)

(2.21 )

CD.md

= 2CD.O(M)

(2.22)

(LjD)

_
md -

1
2y'K(M)CD.O(M)

(2.23)

For this case the logarithmic derivative CDL is

CDL

10

2K(M)Ci
CD.o(M) + K(M)Ci

= 2 lnduced

Drag
Total Drag

(2.24)

n
M,

a . Curve I: GDL

=1

b. Curve 11: GDM

=1

M2 M3 M4 Ms M6
Mach number

Figure 2.3: Graphical construction of conditions GDL

= 1 and

GDM

=1

11

2.0

1.0

LI D
(L / Dl md

0.8

0.6

0.4

0.2

0.5

1.0

Figure 2.4: Generalized drag properties for parabolic polars

and for CDL = 1 the well known result is obtained that the induced drag and the zero-lift
drag are equal. The relationship between CDL and the lift coeflicient is:
(2.25)
and the lift to drag ratio-in relation to its maximum value-amounts to

LID
(LID)md

(2.26)

These interrelationships are depicted in a generalized form in Figure 2.4, which applies
to low subsonic as well as high subsonic speeds, provided that drag polars can be approximated by Equation 2.20 at any Mach number (including the drag rise). However,
for a cambered wing the term drag coeflicient at zero lift is not very meaningful. For the
purpose of performance computation at subcritical speeds Equation 2.20-representing a
symmetrie drag polar- is usually adequate. For high speeds a bet ter approximation is
for example

CD

= CD(M) + K * (M)(CL - Ci:?

(2.27)

with Ci: denoting the lift coefficient for minimum drag coefficient CD' and K* denotes a
modified (higher) induced drag factor than K used in Equation 2.20.

12

2.2

Constraints on Altitude or Thrust

Altitude constraint: in a practical situation where for an aircraft (given AUW) flying
at a specified altitude the condition of minimum drag is sought, the results of the previous
section are not valid if the resulting Mach number is in the drag rise. Since for this
condition L = W dictates that
CL M 2

WjS
-1-

= constant

(2.28)

the following constraint on CL and M applies

dlog CL

+ 2d log M = 0

(2.29)

Substitution into Equation 2.17 yields

(2.30)
For low-subsonic speeds CDM is usually negligible-except for high lift coefficients-and
this result is identical to Equation 2.18. As soon as the drag rise is entered, CD M increases
and hence C DL > 1 and CL > CL,md. This is illustrated in Figure 2.2b where curves of
constant CLM 2 have been drawn and their points of tangency to the Lj D contours have
been constructed, defining the maximum according to Equation 2.30. The example shows
that for M > 0.65 Equation 2.30 begins to diverge from Equation 2.18; this di vergen ce
starts at point A, corresponding to CL M 2 = 0.215. In the present example the altitude for
which Equation 2.30 has to be applied is defined by p < 2Wj(0.215S), or p < 6.64WjS.
For higher cruise altitudes-hence lower ambient pressures-than point A the con dition for maximum LI D at given altitude is found by intersecting the appropriate curve
CL M 2 = constant with Equation 2.30, as depicted in Figure 2.2b. The flight Mach number for maximum Lj D corresponding to this altitude is lower than the Mach nu mb er
for CL,md defined by Equation 2.11. For example, in the same figure a selected condition CL M 2 = 0.35 intersects Equation 2.30 in Point B-the combination CL = 0.57,
M = 0.78-where Lj D :::::: 17.0. This lift coefficient is considerably higher than CL,md for
the same Mach number and hence the conditiuons for minimum drag are not correctly
calculated when Equations 2.10 and 2.11 are used.

Thrust constraint: this case is important when a flight condition is to be determined

which results in minimum drag for given engine rating, since it will result in minimum
thrust and hourly fuel consumption. An engine thrust rating limit is approximately
equivalent to a constraint on the corrected thrust 2 ,
D
J = 1"
= 1I 2Ps/M2 CDS = constant

T
2 For

(2.31 )

details see Section 4.1.

13

This constraint can be interpreted as

CD M 2

= constant

-+

d log CD

+ 2d log M = 0

(2.32)

Expanding dlog CD in terms of the CL - and M-derivatives in accordance with Equation 2.16 yields the following condition for minimum drag:
(2.33)
After substitution into Equation 2.17 the result appears-perhaps not surprisingly-to
be identical to Equation 2.30:

. For a selected value of CD M 2 = 0.02 point C in Figure 2.2b corresponds to this condition.
If the aircraft flies in the drag ri se this condition again does not result in the Minimum
Drag condition defined by Equations 2.10 and 2.11.
In conclusion we find that for a given Mach number the condition for maximum L/ D
is CDL = 1, whereas for an altitude or thrust constraint one finds CD L = 1 + 1/2CDM" At
subcritical speeds these two conditions are compatible only for oCn/oM =
At high
subsonic speeds they result in different combinations of lift coefficient and Mach number,
but always with a lower L/ D than at subcritical Mach numbers. However, the conclusion
that optimum flight should take place at Point A in Figure 2.2b is not justified, because for
jet and turbofan powered aircraft the engine efficiency is also variabIe with Mach number.
For long range flights the condition for minimum drag is not the most interesting case, as
wil! be shown in the following chapter.

o.

14

Chapter 3
The Parameter ML/D
The Specific Range 1 sn is the instantanuous value of distance covered per unit quantity
of fuel consumed for a given aircraft weight, speed and altitude,

sn~f b.R = b.Rlb.t = V

(3.1)
b.WF
b.WFlb.t
F
wh ere F denotes the fuel weight consumed per unit time. Integration of the sn over a
given fiight trajectory-for a change of aircraft weight equal to the fuel consumed-gives
the range. Introducing the Thrust Specific Fuel Consumption TSFC for jet or turbofan
engines:

CT

~ f.-

(3.2)

T
and noting that in steady cruising flight it can be assumed that Lift is equal to Weight
and Thrust is equal to Drag, the SR can be rewritten as follows:

SR

asl

M LI D

(3.3)

W CTI,[
where asl and () denote the speed of sound at Sea LeveliSA and the relative ambient
temperature, respectively. The term cTI.,f is usually referred to as the "Corrected
Specific Fuel Consumption". For a given engine rating and Mach number this parameter
is constant in the stratosphere, and its variation below the tropopause is generally small for
altitudes above 9150 m (30000 ft). Although the variation of TSFC with Mach number
is by no means negligible-especially for high bypass ratio turbofans-the parameter
M LI Dis often used for defining the Mach number for maximum SRThe resulting fiight
condition appears to be fairly accurate for idealized jet engine propulsion, provided the
Mach number range under consideration is not too large 2 In this section the condition
for maximum M LI D will be analyzed in terms of logarithmic derivates:
d log M

+ d log CL -

d log CD

=0

(3.4)

1 Since the effects of wind are not considered in this report the Specific Range is identical to the Specific
Air Range.
2The Specific Range treated in Chapter 4 is a more general criterion for optimum range performance.

15

M~=4

1.2

10

12

1.0

14

0.8

0.6

0.2

o L -____

L -____L -____L -__~L__ _~L_~~L__ _~

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

Mach number

Figure 3.1: The factor M LI D at subcritical speeds

Expanding d log GD in terms of the lift coefficient and Mach number using Equation 2.16
one finds

(3.5)
where the logarithmic derivatives GDL and GDM are defined by Equations 2.14 and 2.15.

3.1

Maximum ML/D at Subcritical Speeds

Although compressibility drag will be present at high subsonic Mach numbers, the case
GDM = 0 wiU be treated as an instructive example. Figure 3.1 shows contours of constant
M LID in the lift coefficient versus Mach number plane for the case of a single parabolic
drag polar. For given lift coefficient the lift/drag ratio is independent ofthe Mach number,
and M LID increases monotonically with the Mach number, leading to the following
observations:
For each Mach number the minimum drag has the same value---defined by GDL
or GL = GL,md-and M(LI D)md increases proportional to M.

16

= 1,

1.0

I'\~=b

n
1

0.9

CL

0.5

0.2

0.1 ~_---'-_ _::-':-_---''--_~_ _...l..-_---,J,=-_-,-_-=n:....,..t._--L.---L-..L-_.l..-_--.J

0.4
0.5
0.6
0.7
0.8 Mach number 0.9

Figure 3.2: The factor M LI D with compressibility effects

17

 For a given lift coefficient there is no partial optimum for the Mach number and the
parameter M L / D does not have an unconstrained maximum.
This means that-if the engines deliver unlimited thrust and TSFC is not affected by the
Mach number-the aircraft will continue to improve its sn with increasing Mach number,
and for constant lift coefficient the altitude will also increase. However, in practice there
will always be at least one of the following constraints.

Altitude constraint: CLM 2 is constant and Equation 2.29 will apply. Substitution
into Equation 3.5 yields

CDL = 1/2

(for CDM = 0)

(3.6)

This result is identical to the weIl known maximum of CL/Cb- For a parabolic drag polar
Equation 2.25 yields accordingly: CL = CL,md/v3, as indicated by Line III in Figure 3.1.

Thrust constraint: the condition C D M 2 = constant applies and Equation 2.32 will
apply. Substitution into Equation 3.5 yields

CDL = 2/3

(3.7)

(for CDM = 0)

which is identical to the weIl known maximum of cl/Cb. For a parabolic drag polar the
re sult is: CL = CL ,md/V2, indicated by Line IV in Figure 3.1.

In conclusion we find that the case CDM = 0 provides useful results only for constrained optimization. For parabolic drag polars the result is summarized in the following
tabie.

I Curve in Fig. 3.1 I

I case
M=constant
W/=constant
T / =constant

1
1/2
2/3

CL,md
CL,md/.J3
CL,md/.J2

I
III
IV

These conditions are clearly incompatible, confirming the absence of an unconstrained

optimum for this hypothetical case, which is treated frequently in textbooks on flight
performance.

3.2

ML/D with Compressibility Effects

Partial and unconstrained optima: the situation depicted in Figure 3.1 remains
valid up to the Mach number where compressibility effects on the drag are becoming

18

manifest, i.e. in the case illustrated by Figure 2.1 up to M ~ 0.65. Figure 3.2 illustrates that contours of constant M L/ D are drastically different in the drag rise: they
become closed curves and there appears an unconstrained optimum for the lift coefficient
(GL = 0.465) and the Mach number (M = 0.795)-approximately 0.015 below the drag
divergence Mach number-where M L/ D reaches a value of about 14.0: Point A. For this
flight condition L/ D = 17.61, or 95.7% of (L/ D)max = 18.40.
The partia1 optimum for given Mach number (Curve I) is obvious1y the same as that
for maximum L/ D at given Mach number, but due to compressibility effects there is
now a1so a partia1 optimum for given lift coefficient: CD M = 1 (Curve II). lts governing
equation can be derived from Equation 3.5. Curves land 11 intersect at Point A, the
unconstrained maximum of M L / D.

Altitude constraint: for constant CL M 2 the optimum is found by substitution of

Equation 2.29 into Equation 3.5, resulting in
(3.8)
identified in Figure 3.2 as Curve IIL For low Mach numbers (CDM = 0) this equation is
identica1 to the resu1t found in Subsection 3.1 but for M > 0.65-in the present examp1e-there is a rapid increase in CDL and in CL, resulting in a bent-up branch of Curve 111,
intersecting Curves land 11 in point A. It should also be noted that for Mach numbers
in excess of Point A the optimum for specified altitude occurs at CL > CL,md, as opposed
to the case of subcritical speeds.

Thrust constraint: the optimum can be found by substitution of Equations 2.32 and
2.16 into Equation 3.5,
(3.9)
This Curve IV-see Figure 3.2- again intersects the ot her partial optima in point A.
This can be verified by combining Equations 3.8 and 3.9 resulting in CDL = CDM = 1.

Graphical solution: although the location of the unconstrained optimum (point A) in

the lift coefficient versus Mach number plane has been illustrated by means of iso-M L/ D
contours, its determination does not require these contours. A graphical solution using
the drag polars-see Figure 2.3-provides a direct solution as follows.
The condition GDL = 1 is identical to aGD/aCL = CD/CL for constant Mach
number. The solution is found graphically by means of a tangent construction, see
Figure 2.3a.
Similarlya tangent construction is applied to aCD/aM
see Figure 2.3b.

= CD/M for

constant CL,

Both conditions are transferred to the CL versus M plane (see Figure 3.2): Curves I
and 11, intersecting in Point A.

19

In conclusion an unconstrained fiight condition for maximum M LI D exists only when

compressibility effects are present , and as a consequence it is located in the drag rise.
For the example aircraft the criterion CD I M = CD I M defines a slope of approximately 0.033, as opposed to the usual definition for the drag di vergen ce Mach number:
CDlM = 0.10. In the next section it will be shown that a similar criterion applies to
the maximum
of high speed jet aircraft.

sn

20

Chapter 4
The Specific Range
The parameter M LID treated in the previous section is useful in the case that variations
in the TSFC are negligible, but this is usually not the case when accurate results are
required or Mach number variations are large. In that case it is appropriate to introduce
the variation of the overall powerplant efficiency,
def

.,., =

Net Propulsive Power

Heat Content of Fuel Flow

TV

TV

= mFH =

FH/g

(4.1)

where mF denotes the mass flow rate of the fuel with calorific value H. The relation
between the overall powerplant efficiency and the TSFC for jet and turbofan engines is
as follows,
M
.,.,=--H/gCTlVO
asl

(4.2)

Substitution into Equations 3.1 and 3.3 yields

SR= V
F

= RH.,.,L/D
W

(4.3)

where RH ~f H / g amounts to approximately 4400 km (2376 nm) for typical gas turbine
engine fueP. The last equation shows that for given aircraft AUW the maximum value of
the Range Parameter,

(4.4)
results in the maximum Specific Range. Since the overall powerplant efficiency has basically the same definition for propeller aircraft and for turbojet or turbofan aircraft, the
Range Parameter is a more significant parameter than M LI D since it covers the complete
range of transport aircraft powerplants. In particular for making a comparison between
lIn that case Equation 4.1 yields 1/ = O.2788Mj(CTj.../).

21

subsonic and supersonic cruising aircraft it is not accept ab Ie to use M LI D instead of the
Range Parameter.
The objective of the cruise point optimization is to find a maximum of the product of
the powerplant efficiency TJ and the aerodynamic fineness LID. Since both quantities are
mutually unrelated functions of the airspeed and altitude, the combined optimum does
not generally coincide with any of the maxima for TJ or LID, provided these exist.

4.1

Generalized Engine Performance

The engine thrust and fuel consumption-and hence the overall powerplant efficiencyare basically functions of the altitude 2 , the airspeed and the engine rating. Textbooks on
gasturbine engine performance show that the generalized performances of turbojet and
turbofan engines can be described in generalized equations as follows,
Corrected Thrust:
Corrected Fuel Flow:

J = fl(M, NIVO)

Ffi

Jy8

= h(M, NI.,fB)

(4.5)
(4.6)

where J denotes the relative ambient pressure and Njv' is the corrected High Pressure
Rotor RPM. The latter can be eliminated to yield

J~ = h(M, TIJ)

(4.7)

and by substition into Equation 4.1 we obtain

-M!!:.!!...
TjJ
-f(MTjJ)
RH Fj(Jv') - 4
,

TJ -

(4.8)

Hence, for variabie engine rating it is found that the overall powerplant efficiency is determined by the flight Mach number and the Corrected Thrust. It can also be shown [25J
that Equation 4.8 holds for any gasturbine-based propulsion system with propeller power
andjor jet thrust output, provided effects due to Reynolds number variation and installation effects are ignored. Figure 4.1 is a representative example of the overall powerplant
efficiency of a present day turbofan engine, showing the effects of Mach number variation
and the engine rating parameter T j (JTto ). The latter compares the Corrected Thrust
to the maximum take off Thrust for the Sea Level ISA static condition, Tto . For typical
turbofan engine cruise conditions this rating parameter varies between approximately 0.8
and 1.1. The maximum cruise rating has also been indicated in the diagram for altitudes
near 35000 ft (10668 m)3. Figure 4.1 shows that the Mach number has a much larger
effect on the powerplant efficiency variation than the engine rating. Hence the case that
the powerplant efficiency is a function of the Mach number only will be considered first.
2For standard atmospheric conditions the altitude defines a combination of ambient pressure and
temperature.
31n fact this curve is a limitation by the engine RPM.

22

0.40
(

0.35

max. cruise rating

@ alt.~ 30 000 ft

'1

1.0

0.30

0.25

0.20
Figure 4.1: Engine performance of a high bypass turbofan engine

4.2

Maximum range parameter

In order to find conditions of maximum SR it appears appropriate to introduce the

following log-derivative, for the case that the efficiency is a function of the Mach number
only:
der

TJM =

d log TJ
M dTJ
=-dlogM
TJ dM

(4.9)

This parameter generally varies between zero for prop uI sion with efficiency independent
of the flight speed, and one for the (theoretical) case of constant TSFC4 . The condition
for maximum SR may now be derived from logarithmic differentiation of Equation 4.4:
d log P

= d log TJ + d log CL -

d log CD

=0

(4.10)

which can be expanded in terms of partial derivatives:

(4.11)
It is noted that for TJM = 1-and hen ce constant TSFC-this equation is the equivalence
of Equation 3.5, the condition for maximum M LI D.
4See Section 4.3 for further explanation .

23

log (0
log IJ

log IJ
shifted
verticaily

log M

Figure 4.2: Graphical construct ion of the condition GDM

= TJM

Partial and unconstrained optima: for CL and M the partial optima are obtained
from Equation 4.11,
constant Mach number:
constant lift coefDcient:
The first of these partial optima is identical to the case of maximum M L / D, the
second ean be rewritten alternatively,
(4.12)
Since both the Mach number and the drag coefDcient are positive and generally CDM is
either zero or positive, there can be two solutions:
1. TJM

>0

-t aCD/aM

> 0, defining an optimum condition in the drag risej

2. TJM = 0 -t aCD/aM = 0, with an optimum defined by the condition CDM = 0 (see

Figure 2.2b), or with an undertermined optimum Mach number in case there is one
(low speed) drag polar.

24

~~~~~------- - ~~~-

The interpretation of this result

is as follows:

If the overall efficiency is independent of the Maeh number, the range parameter
has its maximum value--for any subcritical Mach number-when CL = CL,md' The
maximum specific range has the same value for different altitudes, provided the
airplane tEes at the Minimum Drag Speed and the efficiency is also independent of
the altitude.
If the overall efficiency reaches a maximum value for a subcritical Mach number,
say M~-opt, the maximum specific range occurs at this flight speed, and the pressure
altitude is defined by

2WIS

(4.13)

If the overall efficiency increases continuously with Mach number, the condition
CDM = TJM must be combined with the condition CDL = 1, in order to find the
unconstrained optimum.

The flight condition defined by CDM = TJM can be found numerically or graphically.
An example of a graphical approach is shown on Figure 4.2 where a log-log scale has been
used in accordance with the use of logarithmic derivatives. The points where the slop es
of log CD are equal to the slopes of log TJ define the solution of Equation 4.12.
An example of iso-P curves in the CL versus M plane is shown on Figure 4.3 with
partial optima for CL and M, indicated as Curves I and II, respectively. Although the
general shape of this diagram is similar to Figure 3.2 for M LID there are some noticeable
differences.
Since TJM < 1, the unconstrained optimum (point A) is at a lower Mach number
than point A for maximum M LI D-approximately 0.045 below the drag divergence
Mach number-while the optimum itself is not so sharp.
For subcritical Mach numbers the constrained optima for constant W I and TI Curves III and IV respectively-occur at higher values of the lift coefficient compared to Figure 3.2. This will be explained in the following paragraphs.
Altitude constraint: the maximum value of Pis obtained by substituting CL M 2 =constant (or dlogCL + 2dlogM = 0) into Equation 4.11, which yields
(4.14)
In Figure 4.3 this equation is represented by Curve lIL lts intersection with the approppriate value of CLM 2 defines the constrained optimum, Point B. In the case of a drag
potar without eompressibitity eiJeets (CD M = 0) Equation 4.14 has the following solution,
(4.15)

25

1.0

n
0.9

3.0

CL

0.At<'

"0
0<,

0.8

l\~:3S

0.7

0.6

0.5

0.4

DJ
0.2

n
0.1
0.4

O.S

0.6

0.7

0.8

0.9
Mach number

Figure 4.3: Variation of the Range Parameter

26

Usually the drag polar is a parabola at low speeds, hence from Equation 2.25

= M md

1/4
(

2+T]M
2 - T]M )

( 4.16)

with M md defined by Equation 2.11. The corresponding lift/drag ratio according to Equation 2.26 is

L/ D

= 1/2.)4 - TJ~ (L/ D)md

( 4.17)

This result covers various classicaloptima: CL = CL,md/v'3 for constant corrected TSFC
(T]M = 1), and CL = CL,md for constant overall propulsive efficiency (T]M = 0). For
modern turbofan engines TJM is typically between 0.4 and 0.65 . Equations 4.14 and 4.16
are therefore useful generalisations of the classical flight-mechanical criteria, applicable to
arbitrary propulsion systems, though applicable to low speeds only.
Thrust constraint: the maximum of P is obtained by combining Equation 4.11 and
Equation 2.32, yielding the following result:

_ 2+ CD M
2+TJM

(4.18)

DL -

An example of a constrained optimum is Point C in Figure 4.3. In the case of zero

compressibility Equation 4.18 has the following solution:
(4.19)
At subcritical speeds the drag polar is usually parabolic. In that case Equation 2.25 applies
and the result is the following condition for the optimum lift coefficient,

1 + T]M
1+T)M/2

( 4.20)

where in this case the Minimum Drag Mach number is determined by the drag coefficient
for the minimum drag condition instead of the lift coefficient 6 :

2T/S

(4.21 )

pCD,md

The lift/drag ratio according to Equation 2.26 is

L/D= 1+T]M/2
~(L/D) md
This result is again in accordance with classical criteria: CL
TSFC, and CL = CL,md for constant powerplant efficiency.

(4.22)

CL,md/V'i for constant

5The actual value is not entirely independant of the Mach number; see Section 4.3.
6For further explanation reference is made to [27] .

27

For turboprop engines the engine rating is based on the engine output shaft power
Pbr instead of the thrust. In that case a more logical combination of independent variables
is the Mach number and the corrected power,
(4.23)
The results for the various partial optima may now be derived analoguous to the previous
case. Only one result will be mentioned here: the optimum for constant Corrected Power
Pbr/(VO) = CDhps/a s/M3S, hence dlog CD + 3dlog M = O. For this case Equation 4.11
yields

_ 3+ CDM

DL -

(4.24)

3 +"'M

For aparabolie drag polar and zero compressibility this is equivalent to:
(4.25)
aresult that is directly comparable to Equation 4.20. This case may be assumed to apply
to turboprop or prop-fan powered aircraft where "'M has a small positive or negative value
near the subcritical cruise Mach number.
Multivalued optimum solutions. It is possible that the conditions for partial optima
do not result in a unique solution. For example, supersonic cruising aircraft-not treated
in the present report-may have maxima for LID and VI F both at subsonic and supersonic speeds. Subsonic aircraft mayalso have anomalies in the drag polars, for example
when a shock-free condition with low drag occurs in a small region of Mach numbers.
This behaviour will be clearly manifest in the drag polars, making the analist aware of
the situation. In such a case both local maxima for ."L/ D must be studied in detail and
compared to find the highest value [17] . For supersonic cruising aircraft good subsonic
range performance is also of paramount importance, and both the subsonic and supersonic
optima are of practical significance.

4.3

Interpretation of the Derivative

1]M

The logarithmic derivative TJM introduced in the previous section effectively denotes the
corresponding percentage change in efficiency divided by the corresponding percentage
change in Mach number . It can be related to the corrected TSFC as follows,
_ 1_
M d(CT/VO)
TJM CT/VO dM

(4.26)

This parameter has been used extensively in the analysis to define the optimum flight
Mach number and deserves further explanation.

28

The overall powerplant efficiency may be expressed as the product of a combustion

efficiency 'TJcb, a thermal efficiency 'TJth, and a propulsive efficiency 'TJpr,
'TJ

= 'TJcb

(4.27)

'TJth 'TJpr

Logarithmic differentiation yields

TJM

d log TJcb
dlog M

d log TJth

d log TJpr

(4.28)

+ dlog M + --=-~
dlog M

The fi.rst of these terms is generally nearly zero, while the second will have a small positive
value due to the ram effect at high speeds. The third contribution is the most important
one; its value can be associated with the well known expres sion for the Froude efficiency

2V
'TJpr

(4.29)

= V + Vj = 1 + Vj/V

with Vj denoting the mean exhaust jet velo city. The variation of Vj with the flight Mach
number is not readily determined since off-design engine performance analysis is a rather
complex problem. However, in general it can be observed that for powerplant installations
with a high jet velocity (e.g. idealized jets) this variation is small, while for low exhaust
velocities-in the case of very high bypass engines and propeller propulsion-the speed
increment Vj - V varies little with airspeed. The following hypothetical cases can therefore
be distinguished:
constant

Vj :

constant

Vj -

dlog'TJpr
dlogM

dlog'TJpr
dlogV

= 1- _V_ = 1- 'TJ r/ 2
V +Vj

(4.30)

and

V :

dlog 'TJpr
dlogM

dlog TJpr
2V
d I V = 1 - -V-- = 1 og
+Vj

(4.31 )

'TJpr

Thus it is found that 'TJM = 1 - factor x TJpr, where the factor depends primarily on the
specific thrust. Since for all cases TJpr = 0 for M = 0, TJM approaches 1.0 for M O. For
subsonic cruise Mach numbers the propulsive efficiency of turbofan engines is generally
between 0.4 and 0.8 and hence TJM will vary between 0.2 and 0.8. For propeller powered
aircraft the propulsive efficiency varies only slightly with airspeed and 'TJM is approximately
zero 7 , unless the propellers operate at supersonic tipspeeds. In that case 'TJM will become
negative.
As an example Figure 4.1 shows that for a high-bypass turbofan the constant-T / Jcurves have a nearly constant slope for Mach nu mb ers between 0.6 and 0.9. As the
diagram is a double logarithmic plot, this indicates that an exponential relatioDship for
the corrected TSFC caD be used as an approximatioD,

CT /.JO = constant x Mn

(4.32)

29

0. 8

---_.. -- --'
--..
--- _ -~--

Upper limit of data for Hp =20,000 ft ~ ~

0.6

..... 4'

~-

(.

-:-

0.4

.... " ~"

",,"0 (
I..00
0<0
0.2
(.

1/....

0
./

."

."

./

~_

--~

.... " "

.... .. "

....

Data obtained from engine brochure

for maximum cruise rating 0.6!: M!: 0.9

~ Lower limit

of data
for Hp =r089 ft

-~

o
I

Indicates data points for Hp =30.000 ft

10

12

Bypass ratio

Figure 4.4: The engine efficiency power law exponent

Substitution into Equation 4.26 yields
TJM = 1 - n

(0::; n ::; 1)

(4.33)

Statistical values for the exponent nare given in [2] and have been reproduced in Figure 4.4. In accordance with the derivation given above there appears to be a correlation
between n and the bypass ratio-which is related to the specific thrust. An alternative
relationship between the corrected TSFC and the Mach number,

CT/VO = Cl

+ C2 M

(4 .34)

can be used in Equation 4.26 and it follows that

1

(4.35)

resulting in a value of TJM which is not independent of the Mach number. Although
Equation 4.35 provides an accurate approximation over a large range of Mach numbers,
Equation 4.33 is usually adequate for cruise performance analysis to derive TJM for a
given engine. Figure 4.5 depicts the relation between the overall powerplant efficiency,
the corrected TSFC, the log-derivative and the flight Mach number. It illustrates the
(in)accuracy of the approximation of a constant TJM for high subsonic speeds.
7Note that since TJpr op < 1 Equation 4.31 gives a positive valuefor
value of the simple analysis.

30

TJM;

this demonstrates the limited

4.4

Refinements in the Optima

In the previous section only the primary variabie for the thrust (the Mach number) has
been taken into account. For a more accurate analysis, accounting for the effects of
engine rating on efficiency and of Mach number on the maximum cruise thrust, a further
extension of the analysis may be required.
Partial and unconstrained optima: the effect of thrust setting is taken into account
by means of the following logarithmic derivative:
der

TJT

Olog 17
olog(T10)

TI 0

017

= ---:;J o(T I 0)

(4.36)

(constant Mach number)

and since 17 now depends on two variables we have to modify the definition of TJM accordingly,
TJM

der

a log 17

alog M

M aTJ
17 oM

= - - (constant Tlo)

(4.37)

The condition for maximum SR according to Equation 4.11 must now be modified into8:
dlogP

(I-CDL+TJTCDL)dlogCL

+( TJM

- CDM

+ 2TJT + TJTCDM)d log M = 0

(4.38)

The partial optima with respect to CL and Mare obtained by setting the bracketed terms
equal to zero:

-+

(4.39)

constant CL -+

(4.40)

constant M

The unconstrained optimum can be found when these partial optima are combined. This
will generally require a numerieal iterative solution.
Altitude constraint: the condition CLM 2 constant, or dlog CL
substituted into Equation 4.38 yielding
CDL

= 1 + 1/2

(CD M

-~)
1 - TJT

+ 2dlog M

0,

IS

(4.41)

The solution for the optimum flight condition may now be obtained iteratively by first
assuming TJT = 0, resulting in a first order approximation for the lift coefficient and the
Mach number, and henee Tlo. The engine data-see Figure 4.1-will then provide an
actual value for TJT and an improved approximation ean be obtained.
BThe derivation can be found in [27] .

31

0.5

1.0

(T/{
IJ

[Nih)
N 0.8

0.4

IJ M

0.3

0.6

0.2

0.2

0.1

max. cruise rating

alt.~ 30,000 ft

oIL---...l.----'-----'-:----:-':---~O

0.2

0.4

0.6

0.8

tO

Mach number

Figure 4.5: Efficiency and TSFC of a high bypass turbofan engine

Constrained optimum for given engine rating: although the case of constant TIJ
gives areasonabie approximation for a thrust rating limit, a more accurate analysis may
be required for the following reasons:
An engine rating may refer to an RPM limit, a maximum Turbine Inlet Temperature
(TIT) , or an Exhaust Gas Temperature (EGT) limit. Not only will the result be
a value of Nlv'O which varies with the altitude in the troposphere, but for a given
altitude a constant RPtvi results in a variation of the corrected thrustj see Figure 4.1
as an example.
For accurate performance analysis it will be necessary to account for installation
effects due to inlet total pressure loss, bleed air and power extraction, scrubbing
drag and Reynolds number effects. All these effects will not only re duce the overall
powerplant efficiency, but they mayalso introduce an altitude effect on the efficiency
for given TIJ.
In this case it will be necessary to first calculate the Mach number versus altitudehence lift coefficient-limit by means of a numeri cal or graphical procedure, using the
two equilibrium equations: Thrust=Drag, and Lift=Weight. Since this operationallimit
depends on the aircraft AUW the result is no longer a unique rating limit in the CL versus
M diagram.

32

Attitude

45

1000 FT

m rl
I

40

35

30
/
,/
,/

25

,/
,/

/
/

80%(~)
F max

,/

,/

20

Iw =constant I

m
0.6

0.7

0.8

0.9
Mach number

Figure 4.6: The Specific Range affected by operating conditions

Figure 4.6 shows an example of range performance with the SR plotted versus altitude
and Mach number, for a given AUW. For low AUW the maximum SR appears to be
unconstrained (Point A), but for high weights the cruise thrust limit may force the fiight
condition to a Mach number and altitude below the 100% SR condition.

Non standard conditions: the engine thrust rating and efficiency are dependent on
the Mach number, the pressure and the non standard temperature. Similar to the previous
case a numeri cal analysis is required to identify the best cruise condition for each aircraft
weight and the maximum SR will vary during the cruising fiight. Depending on the
ThrustjWeight ratio of the aircraft, a temperature above standard may force the aircraft
to a considerably reduced altitude or Mach number compared to the 100% maximum SR
condition, resulting in a range performance degradation. This applies in particular to
cruising fiight after failure of an engine.

33

Chapter 5
Analysis of Cruising Flight
Different from the problem treated in the previous chapters, dealing with (optimum)
instantaneous conditions, the cruise range is an integrated performance. Since there are
two control variables at each moment of the flight-we have used the lift coefficient and
the Mach number-optimal control laws can be proposed for these variahles in order to
maximize the range for a given amount of fuel. For the flight crew these controllaws have
different characteristics and it will be shown that several cruise techniques have serious
disadvantages.
Since equilibrium of forces will he assumed during the flight , the result is not necessarily an absolute maximum range. For example, cyclic fiight with dolphin type climb
and descent profiles may, at least theoretically, increase the range; see [16] and [22]. However, in practice civil transport aircraft fiy the cruise sector as a quasi steady fiight , with
small variations in the control variables. We are therefore interested in approximate solutions for the integrated (optimum) range performance, since they will be adequate for
preliminary design purposes. It will be shown that the cruise range can be calculated by
means of a simple analytical equation. Conversely, for a specified range the amount of
cruise fuel required can be derived in closed form, and from these results the block fuel
and the reserve fuel will be obtained readily.

5.1

The Generalized Range Equation

Classical derivations of range are hased on either constant propulsive efficiency for propeller aircraft, or constant TSFC for idealized jet propulsion. The pertinent analytical
equations can be found in most textbooks on performance. They are usually quite accurate, provided the conditions do not vary considerably during the fiight. For example, the
following well known expression for the range-ascribed to Louis Brguet (1880-1955)was derived for propeller aircraft cruising at constant altitude and constant angle of attack,
(5 .1)

34

where "i" and "f" denote the initial and final conditions of the flight . Later, it was found
that a very similar equation applies to the range of jet aircraft in a cruise/climb with
constant speed and constant angle of attack,
(5.2)

In the following text a general expression will be derived for the range of aircraft with
various types of gas turbine based powerplants. The range in quasi steady cruising flight
is obtained from integration of the specific range,

R=

WiV
-dW
w, F

For constant

(5.3)

sn during the flight , the range would simply be:

R= F(Wi -Wj )= FWF

(5.4)

with WF denoting the fuel weight. However, due to the decreasing AUW during the flight,
the fuel flow is actually decreasing and V / F is increasing. This can be taken into account
by using Equation 4.4 for the Specific Range,

R= RH {Wi pdW
Jw, W

(5.5)

with P denoting the Range Parameter. One special solution is found for the case of
constant lift coefficient and Mach number-hence constant P,
(5.6)
This equation is a generalization of Equations 5.1 and 5.2, sin ce it applies to propeller as
weil as jet propulsion.
Due to practical operational considerations it is usually impossible to control the
flight so that P is constant. Several alternative techniques have therefore been proposed,
resulting in somewhat short er ranges than Equation 5.6. Although these derivations are
covered in many textbooks-see also [2]-for the case of zero compressibility drag, the
significant effect of compressibility will be demonstrated as weil.

5.2

Range at Subcritical Speeds

This low speed case is obviously most important for propeller aircraft, since jet aircraft
have their best performance in the drag rise. But even jet aircraft may be forced to fly at
reduced speed and/or altitude, for example when ATC requires this or after engine failure
[15]. Especially in the latter case good off-design performance is important in order to
reduce the required reserve fuel penalty.

35

In the following derivations certain combinations of terms are often used and will be
designated by the following notations ,

(~WFjWi = 1- WfjWi

(5.7)

y ~f CL)CL,md

(5.8)

For subsonic speeds below the drag rise the drag coefficient is usuallya unique function
of the lift coefficient, hen ce CL,md is a single number.
Cruise/climb (cic): flying with constant CL and M results in a constant CD-hence
constant liftjdrag ratio-and constant overall efficiency and corrected TSFC. Thus Pis
constant and equal to its initial value, whence
WidW
Rele = RHPi

W,

= RHPi 10g(1 -

(5.9)

()-1

This flight profile requires the cruise altitude to be steadily increased as fuel is consumed
in such a way that W / and T / remain constant. The relative ambient pressure must
therefore decrease proportional to the AUW and in the (isothermal) stratosphere the
aircraft flies with constant engine setting. The extra thrust required to increase altitude
can be taken into account readily, but for a fair comparison with alternative schedules-at
constant altitude-it is preferred to avoid this, assuming that after the cruising flight the
aircraft can glide back to the initial altitude, thus regaining the small range loss. As
mentioned before the cruise/climb technique will not be acceptable in many situations
because of the requirements of ATC.
Horizontal cruise, constant lift coefficient (h,l): this requires the speed of the
aircraft to be steadily reduced as fuel is consumed, in such a way that V is proportional
to VW. Again L j D is constant and for jet aircraft the Range Parameter varies only due
to the decrease in Mach number, and hence the overall efficiency. The engines have to be
throttled back so that T varies proportional to W in order to keep T jW constant. If we
assume that for that for relatively small speed variations the log-derivative TJM remains
constant, the following result is found [28]:
Rh,l

= 2RHPi {1 -

(1- ()~M/2} /TJM

(5 .10)

This equation covers idealized jet propulsion with TJM = 1, propeller powered aircraft
with TJM = 0, as well as turbofans with intermediate values of TJM; see Section 4.3. For jet
propulsion Equation 5.10 complies with the classical "square root range equation" which
can be written as
R

= 2RHPi{1 - .JI=(}

For TJM .!. 0 the numerical values become equal to Equation 5.9
This flight procedure has several major disadvantages:
lThis can be proven by means of de I'Hpital's rule.

36

(5.ll )
1

 For jet propulsion there is a considerable loss in range relative to the cruise/climb
procedure due to the deterioration of total efficiency with reducing speed.
The flight time is increased.
Steady decrease in speed and continuous variation of the engine setting during the
cruise are unlikely to be acceptable to airlines as a normal operational procedure.
Horizontal cruise, constant Mach number (h,M): this is a practical procedure
from the operational point of view, because for given altitude a constant Mach number
results in constant air speed. During the flight the drag is decreasing due to the weight
reduction, and the engines have to be throttled back. This is easily accomplished with
modern flight control systems. Integration of the SR results in the following equation
[19] for a parabolic drag polar,

Rh,M

= 2RHTJi(L/ D)md{ arctan y -

arctan y(1 - ()}

(5.12)

The difference between the two angles in this equation can be evolved by means of goniometric equations into a single angle. After introduction of the Range Parameter for the
initial cruise condition one finds:

= RHPi (y-1 + y) arctan {y-1 + ;(1 _

()}

(5.13)

Since for optimum cruising usually CL,; < CL,md, hen ce y < 1, the angle defined by the
arctan is small, and the arctan can be set equal to its argument. This yields a very good
approximation of Equation 5.13,
(5.14)
For y ::::::; 1 the range is only slightly shorter than for the cruise/climb, but the penalty
due to the altitude constraint increases when CL,i < CL,md' If the available engine
thrust limits the initial cruise altitude appreciably, it is the usual practice to execute
intermediate step climbs to increased flight levels, the co-called "stepped cruisejclimb, "
in order to approximate the continuous cruise/climb as closelyas possible.
Horizontal cruise, constant engine setting (h,es): the decreasing AUW results in
an increasing flight speed during the cruise. Formally, the range can be written as follows:

Rh,es = RHPi

J1 (TJ/TJi)d(W/Wi )

(5.15)

1-(

Due to the increasing speed the overall efficiency increases as weIl, but an analytical solution is rather complicated. For the special case of constant TSFC a graphical procedure
has been proposed by Peckham [19]. No attempt has been made here to derive a general
solution, since the following considerations would make such a solution virtually useless.

37

 ATC requirements and practical flight procedures disfavour the variation in flight
speed .
The present section deals with subcritical speeds, but it is very likely that during
the flight the aircraft will enter the drag rise.

Initial conditions: Equations 5.9 and 5.10 show that the maximum range is obtained
for initial flight altitude and Mach number such that the initial SR is maximal. For the
horizontal cruise with constant Mach number, however, Equation 5.13 contains not only
Pi but the ratio y = CL,dCL,md as weIl. In other words, the mean value of the SR should
be maximized, not the initial value.
Differentiation of Equation 5.13, assuming a specified Mach number, yields an optimum
altitude defined by
1

Y=

(5.16)

.11-

For a specified a/titude the optimum Mach number is determined by

2 -7]M

1/2

y = { (2 + 7]M )(1 - () }

(5.17)

These two conditions are incompatible for jet propulsion (7]M > 0), the same finding of
Section 3.1. The aircraft will thus tend to improve its performance by flying faster and
higher, until an engine thrust limit is met, which forms the ultimate constraint. The best
range performance will than be obtained for an initiallift coefficient determined by
1
Y

= -..,;r;=:(l=+=7]=M:;=;:)(=1_=(~)

(5.18)

Contrary to what is usually suggested in the literature this equation shows that for long
range flight the optimum initial lift coefficient can be quite close to the minimum drag
condition. For example, a fuel fraction of 30% for a turbofan powered airplane, with
'f/M = 0.60, yields an optimum y = 0.945. The optimum flight speed is than only a few
percents above the Minimum Drag Speed.

5.3

Range at High Speeds

Figure 5.1 shows an enlarged sector of Figure 4.3 near the unconstrained maximum of the
Range Parameter. The flight condition for maximum SR is Point A, corresponding to
Point A in Figure 4.6. The absolute maximum range will be obtained in a cruise/climb
with constant P for every AUW, resulting in
(5.19)

38

buffet
limit
engine
rating

-H

Figure 5.1: Flight schedules for high speed cruising

Although in the present example the available thrust allows the aircraft to fly continuously at this optimum, this procedure is usually not selected in operational practice. For
example, long range cruising is carried out in a flight condition where the specific range
is 98% or 99% of its maximum value. The small fuel penalty is considered acceptable in
view of the time saved due to the higher cruise speed. For short haul fiights a larger fuel
penalty can be accepted and the aircraft is flown at the maximum cruise rating of the
engines at a lower altitude. These practical flight conditions are referred to as Long Range
Cruise (LRC) and High Speed Cruise (HSC) conditions, respectively. Gusty conditions
may as weIl force the pilot to reduce the cruise lift coefficient below its optimum value
by flying at a reduiced flight level in order to avoid buffeting problems and to improve
comfort.
The initial condition ofthe flight-denoted "i" in Figure 5.1-has been selected so that
several flight schedules are possible with only slightly different ranges. In a stepped climb
the idealized cruise/climb can be approximated, where the aircraft stays in Point "i". A
horizont al cruise with constant lift coefficient, line "a", carries the aircraft almost "over
the top" of the sn and yields a good range, as opposed to the low speed case. In spite of
this, the speed reduction during the flight is usually objectionable for practical operation.
A horizont al cruise with constant Mach number, line "b" , is another schedule "over the
top", which may even be superior to the cruise/ climb staying in Point "i". Cruising
with constant engine setting, line "c", is the only schedule which brings the aircraft to

39

significantly lower Range Factors because it penetrates the drag rise deeply.
Similar to flying at subcritical speeds the constant altitude and Mach number cruise
appears to be preferable. Different from the low speed case, the latter flight schedule
does not necessarily have a significant unfavourable effect on the range. Equation 5.13
will therefore be selected 2 for further use to calculate the fuel load. Again the question
is: at what altitude and Mach number should the cruising flight be started? Although
for given Mach number Pi is maximum for CDL = 1, hence y = 1, the reduction in lift
coefficient during the flight brings the optimum theoretically to a higher value. However,
the high "optimal initial altitude" may not be achievable in practice due to a thrust limit,
a pressure cabin limit or a required buffet margin. For practical performance analysis it
is appropriate to select the intersection of the 98% or 99% maximum VI F condition and
the cruise thrust limit as the initial flight condition.

2This equation was derived for parabolic drag polars. Since during the fiight the Mach number is
constant, the drag polar may be different from the low speed case, but it usually remains a parabola.

40

Chapter 6
Prediction of the Fuel Load
A typical mission profile for an international flight is shown on Figure 6.1. Total fuel is
subdivided into

block Juel required to fly the complete mission, including fuel to taxi out before
take-off and taxi in after landing, and
reserve Juel, in accordance with the pertinent rules for the operation under consideration.
For the determination of the Take Off Weight (TOW) taxi fuel is not included since
the TOW is defined at the runway threshold, before starting the take off. Taxi fuel after
landing is usually taken from the reserve fuel , except for multi mission operation. Mission
Juel includes fuel required for
take off, acceleration and climb to cruise altitude,
cruising flight (with or without steps),
descent, approach and landing,
manoeuvering.
The AUW of a transport aircraft is subdivided into Operating Empty Weight (OEW),
Payload Weight and Fuel Weight. The payload plus fuelload is referred to as the Useful Load (UL). For a given aircraft, the OEW is considered to be independent of the
range. The payload is limited (see Figure 6.2) by the Volumetrie Payload-determined
by the number of seats and volume of the cargo holds-or by the Maximum Zero Fuel
Weight (MZFW), a structurallimit. In the early design stage maximum payload is determined directly by the design specifications. For short ranges the required mission fuel
is determined by the Landing Weight (LW), which in turn depends on the payload.
Figure 6.2 illustrates that the mission fuel versus range is located within an envelope,
initially determined by the case of zero payload and by the MZFW, assuming that the
latter forms the critical payload limit. The AUW increases with range until it reaches

41

MISSION

t:

RESERVES

Figure 6.1: A typical flight mission

42

FUEL LOAO

reserve tue I
RANGE
I
I

I
I

I
I

PAYLOAO
\

21 cu
~ I go

,,
,, I
,

:UI~

~I

MZFW

~I
~I

-=1

...
111

~I

~I

cu

'":1
1:1

~I

~I

111

:;:

EI
I
I
I

1:'
EI
I

,I

:;

I
I
I

I
I

Ie

cruise
techniques

cu
c:

en

cu'
enl

I
' ,

~, I
"+
. '
~ I',

"

"'I

I(

""-

"-

""-,10

RANGE

Rh

Figure 6.2: Payload and fuelload envelopes versus range

43

its limit, assumed here as the MTOW. The maximum range with Maximum Payload is
referred to as the Nominalor Harmonic Range, Point A in Figure 6.2. For increasing
ranges the fuel required is determined by the MTOW, which is independent of the range.
The highest payload is than obtained by subtracting the Fuel Weight plus the OEW from
the MTOW.
A further increase in range and fuel is limited by the fuel tank capacity corresponding
to the maximum useful range: Point B. The maximum range, Point C, is obtained by
decreasing the MTOW through progressive payload reduction, since the fuel fraction can
be increased only by reducing the initial weight. In the absence of a fuel tank capacity
limit the ultimate range 1 , Point D, would be reached, where the total fuelload is equal
to the MTOW minus the OEW.
The payload versus range and fuel versus range diagrams dep end on the cruise technique. In particular there is a considerable difference between Long Range Cruise (LRC)
and High Speed Cruise (HSC) techniques; see also Section 5.3. Furthermore there are
effects of atmospheric conditions, assumptions with respect to wind, climb and descent
flight profiles, reserve fuel policy and various payload configurations of the aircraft, which
make it difficult to compare otherwise similar aircraft on their payload versus range diagrams. Nevertheless, it will be shown in the Appendix that these diagrams can provide
useful information.

6.1

Mission Fuel

Accurate calculation of the contributions to the fuel load required for various flight segments can only be done when sufficiently detailed data are available in the form of drag
polars, engine thrust and fuel flow diagrams, design weights, etc. However, in the conceptual design stage the AUW is to be estimated from elementary information such as
the maximum payload, the range, some preliminary engine data and basic aircraft dimensions. In that case it will often be acceptable to estimate the fuel weight fraction WF/Wto
on the basis of primarily the cruise fuel requirements. Additional allowances can be given
for other flight segments, in particular climb and descent. A major factor to be used as
input for such a method will be the range factor P = TlL/ D at the initial cruise altitude.
A practical procedure to calculate the mission fuel is to determine the amount of
fuel required for a hypothetical cruising fiight over the complete specified (still air) mission
range, starting the cruising flight with TOW. An extra amount offuel ~WF is than added
to allow for excess fuel consumed during take off, climb to cruise altitude and acceleration
to cruise speed, referred to as "lost fuel" [2].
The analysis of range performance has resulted in expressions for the range according
to several flight schedules. Recommended expressions are
lSometimes referred to as the "Ferry Range", flown for special purposes with extra fuel tanks installed
in the fuselage.

44

 Equation 5.14: subcritical speeds, single parabolic drag polar, for a fiight schedule
with constant altitude and constant Mach number.
Equation 5.19: the maximum achievable range for a high speed cruise/climb schedule, a significant result for theoretical work.
The first of these equations is also applicable to high speed cruising. A slightly more refined version is proposed to cover all practical cases by using a "cruise control factor kR "
This factor has been developed to offer the possibility to compare different cruise procedures, using the same basic equation for the cruise Ju el load,

(6.1 )
where Pi denotes the Range Factor TJLI D at the initial cruise condition, and the cruise
control factor is
RIRH
k R = 1 + 6P

for a cruise/climb

RIRH 2y 2

kR = (1 - -6P )---2 for honzontal cruise, constant Mach number

i

l+y

For a step/climb kR will have a value in between these two; a simple assumption kR = 1
will be adequate.
The following advice for the initial conditions is based up on the previous analysis.
For subcritical cruise speeds use Equation 5.18 to find the optimum initiallift coefficient for a horizontal cruise with a thrust constraint. The initial Mach number
is than obtained from the equilibrium condition Thrust=Drag, Equation 4.21. For
reasons of speed stability it should be above the Minimum Drag Mach number.
For high speeds the condition for the unconstrained maximum Range Factor should
be determined first, as weIl as the contour of 99% or 98% of this value. The intersection of this curve with the thrust limit is a good initial cruise condition; see also
Figure 5.1.
For a more superficial approximation one may simply use kR = 1, and the Range Parameter for LRC is assumed equal to 98% of its maximum value.
The lost Ju el may be obtained from an energy balance [28],

6,.WF TJcl R H

Wto

= Wtohcr + 2;v.:r = Wtohe,cr

(6.2)

where TJcl denotes the mean value of the overall efficiency during take-off and climb to
cruise altitude. Although TJcl is dependent on the speed schedule selected for the climb

45

and on the type of propulsion, it may be assumed equal to a certain fraction of the value
during cruising fiight. In terms of the energy height at the top of the climb,

v;,;

(6.3)

hecr= h cr+'
29
one finds [28]
~WF

- - ~ (1.1
W to

he cr
+ 0.5'TlM)-R'

(6.4)

'Tlcr H

The additional Iud to account for manoeuvering [24] amounts to

~WF

-w: ~ 0.0025/'Tlcr

(6.5)
to
Adding these losses to the cruise fuel results in the mission fuel for the mission range

Rm,
Rm/ RH
Pi

+ (1.1 + 0.5'TlM )he, cr / RH + 0.0025

+ 0.5kRRm/ RH

'Tlcr

(6.6)

In accordance with procedures sometimes adopted in airline evaluations of proposed aircraft the fuel consumed during descent, approach and landing is assumed to be equal to
the fuel used during a cruising fiight over the same distance. It is thus assumed that a
hypothetical cruising fiight extension up to the field of destination is made without any
further allowances for the actual fuel used, generally a conservative approach.
A further simplification can be made by writing Equation 6.6 as follows:

WF,m
W to

Req/RH
Pi

(6.7)

The equivalent range Req is obtained from Equation 6.6

Rm

Req = 1 + 0.5k (Rm/ RH )/P + {(1.1 + 0.5'TlM )he,cr + 0.0025RH }(L/ D)i
i
R

(6 .8)

The second term, accounting for the "lost fuel," amounts to approximately 5% for very
long range fiights (~ 12000 km) up to 25% of the first term for relatively short ranges
(~ 1200 km). For eveny short er ranges this term will become dominant, but due to its
low accuracy a more accurate (numerical) procedure will than be required.

6.2

Reserve and Total Fuel

Reserve fuel will usually include fuel required for

a diversion fiight over a specified distance, including a missed approach and overshoot manoeuvre,

46

 a holding flight of specified duration at a specified altitude,

contingency fuel, a certain percentage of the block fuel.
an extended duration of the flight.
A reserve policy will combine several of these contributions, in accordance with government al or company rules and dependent on the type of operation. There are important
differences in reserve fuel policy between domestic and international fiights on the one
hand, and between different nations and operators on the ot her hand. For example, the
Association of European Airlines (AEA) specifies
a 200 nm (370 km) diversion flight for short and medium range aircraft, or a 250 nm
(463 km) divers ion for long range aircraft,
30 minutes holding at 1500 ft (457 m) altitude,
5% of the mission fuel for contingency reserve.
For flights in the USA typical reserves are less , for example 130 nm (241 km) diversion
and 30 minutes holding at 1500 ft (457 m) altitude. For business aircraft the reserves
are in total frequently specified to be equal to a 3/4 hour extension of the cruising flight.
Long range international flights mayalso require an extension of the cruise, for example
one hour at LRC flying.
Reserve fuel estimation for preliminary design should be based on a method which
allows the user to insert the actual reserve policy for the case under consideration. This
results in the selection of several of the following allowances in accordance with the design
specification:
1. A diversion distance equal to D.Rdiv flown at the landing weight, expressed as an
increase of the mission range,

D.R

= (r D.R)div ~;nd

(6.9)

YVto

where the factor rdiv accounts for the fuel penalty due to the overshoot, and the
deterioration in the
due to flying at reduced speed and altitude.

sn

2. An allowance for a holding period of D.thold flown at the landing weight, in terms of
a mission range increment ,

D.R

= (r V D.thold ~;nd

(6.10)

YVto

where the factor rhold accounts mainly for the reduced powerplant efficiency due to
the low altitude and RPM during holding.

47

3. A contingency juelfraction t..WF,etg/WF,m equal to 0.05 or 0.10.

4. An extension time of the cruising flight t..t en resulting in an increment of the design
range t..R = Vcrt..ter.
The total fuel load is obtained from summation of mission fuel (Equation 6.6) and
the various reserve fuel allowances discussed in the previous paragraph. It is appropriate
here to introduce the equivalent all-out range Pi for the case that diversion, holding
and contingency fuel are specified,

n = Req { 1 + t..WF,et
W

g}

~m

l-Viand + (V t..t )cr

+ {( r V t..t )hold + (r t..R)div} -w:

(6.11 )

The holding speed is generally about 50% of the cruising speed, while Wland/Wto =
1- WF,m/Wto. The penalty factors rhold and rdiv may require a rather complex calculation
procedure, but the following rule of the thumb is proposed as an alternative:

rhold

= rdiv ::::: 1.10 + 0.57]M

(6.12)

If the reserves are specified as a continued duration of the cruising flight , the equivalent
all-out range is simply

n=

Req

{I +

er
t..t }
ter

(6.13)

Finally the total fuel fraction is obtained from

WF = n/RH = R/RH

Wto

Pi

(TJL/D)i

(6.14)

The advantage of these analytical equations is that they provide a closed-form, noniterative, and therefore computationally very efficient method for calculating the total
fuelload fraction.

48

- - -- -- - - - - - - - - - - - - - - - - - -

- - - - - - -- - - - - - - - - - - -

Chapter 7
Summary of Results
In spite of the large amount of literature on optimum range performance there is still no
undisputed generalized approach. The present report forms a basis for such a solution.
Analytical criteria are presented for the unconstrained optimum of the Specific Range,
which-in the case of tubofan powered aircraft-occurs at high subsonic speeds. For
subcritical speeds only constrained optima exist. All criteria are given in relation to the
powerplant overall efficiency and its derivative with respect to the Mach number and
engine rating. They apply to any presently existing gas turbine based propulsion system
for subsonic aircraft. Most criteria found in the literature are special cases of the general
theory presented in the present report.
Logarithmic derivates have proved to be very useful for the present analysis to derive
optimum cruise conditions; they are defined as follows:

CDL

= Blog CD/Blog CL

CDM = 8logCD/BlogM

TiM

= Blog Ti/Blog M

The present approach obviates the classical assumptions made in many previous publications on cruise performance: compressibility effects are absent, the drag polar is parabolic,
the overall powerplant efficiency or TSFC is constant.
The range parameter

is a non-dimencional generalized cruise performance index,

= TiL/D

The Specific Range SR appears to be proportional to this parameter, according to Equation 4.4,

v = RH p

49

The combined aerodynamic/propulsive quality of aircraft with widely different AUW

can be compared on the basis of this Range Parameter. For a given aircraft/engine
combination P is a unique function of the lift coefficient and the Mach number only,
and it appears advantageous to represent it in a CL versus M diagram; see for example
Figure 4.3.

For subcritical speeds where a single drag polar is used-e.g. for propeller aircraftan unconstrained optimum is possible only for TJM = 0, hen ce constant overall efficiency or
an efficiency which reaches a maximum value. In the lat ter case both the lift/drag ratio
and the overall efficiency can have a maximum value at the same altitude and speed,
provided the altitude is chosen appropriately, according to Equation 4.13. For turbofan
engines TJM > 0, and there will only be a constrained optimum,

= 1- TJM/2

or, in the case of a parabolic drag polar:

= 2/(2 + TfM)

or, in the case of a parabolic drag polar:

attitude constraint: CVL

thrust constraint: CVL

C - C
1
L- L,md)l+TfM
The latter is the most significant of the two constraints. For TJM
accordance with classical criteria for constant TSFC.

= 1 these results

are in

Unconstrained optima for high speeds: the drag coefficient, the lift/drag ratio and
the range parameter are essentially a function of the lift coefficient and the Mach number.
Optimum conditions have to be defined in terms of both CL and M, as opposed to the
case of a single drag polar. Conditions for the unconstrained optimum are

CVM

= TJM

-+

aCv/aM

= TJMCv/M

defining a flight condition in the drag rise, slightly below the drag divergence Mach number.
Reference is made to Figure 7.1, depicting the following unconstrained optima:
propeller aircraft with constant overall efficiency: TJM
Point Al ,

defining

(L/ D)max at

idealized jet propulsion with constant TSFC: TfM = 1 resulting in aCv/aM = Cv / M,

hence (M L/ D)max, Point A 2 This curve in Figure 7.1 can also be found directly
from the drag polars; see Figure 2.3.

50

0.9

0.8

0.7

0.6

0.5

[0 -

L - 1. Y2 ([0 -I] )
O.4 ---_
M M

98% (I] L/Ol max

--

__~!:.= Y;i:-(~--lt.- ---~=====_//

0.3

-------'-------

0.2 L - _ . - L_ _...L_ _L - _ . - L_ _...L---'---LJL--L......L...lL_..l.-_---'

0.9
0.7
0.8
0.5
0.6
Mach number

Figure 7.1: Overview of the partial and constrained optima

turbofan propulsion: 0 < TJM < I, the maximum for the SR at Point A3. Typical
values for TJM are between 0.4 and 0.8; see Figure 4.4.
Although Point A3 is theoretically the condition for maximum Specific Range it is usual
to fly long distances at a somewhat higher Mach number, thereby accepting one or two
percent increase in fuel consumption: the Long Range Cruise condition, LRC.
Constrained optima for high speed flight: operational considerations~ . g. Air
Traffic Control or a cabin pressure limit-may force the aircraft to fly at a specified
altitude, imposing a constraint on GL M2. In combination with the general requirements
for maximizing P defined by Equation 4.11 this constraint is: GDL = 1 + (GDM - TJM)/2.
In Figure 7.1 this equation must be intersected with the specified value of GL M 2 , resulting
in the optimum Point B.
A similar constraint due to the engine thrust limit on TIJ, or on GD M 2 , results in
the following condition: GDL = (2 + GDM )/(2 + TJM) , which defines the optimum Point G
aft er intersection with the specified value of GD M 2 .
Ot her factors affecting the optimum may be associated with:
the occurrance of multiple optima, see Section 4.2,

51

 refinements in the analysis associated with realistic engine rating limits and installation effects, see Section 4.4,
operational and/or economie considerations which may result in off-optimum fiight
conditions.
For example, Point D (M = 0.8, CL = 0.45) in Figure 7.1 could be a suitable initial value
for long range cruising, since excursions to higher or lower lift coefficients will not degrade
the SR by more than a few percent of its maximum value. If for operational reasons the
cruise Mach number has to be reduced, the SR is in fact slightly improved. Final selection
of the High Speed Cruise condition (HSC) will be based on economie considerations which
are outside the scope of this report.
Range in cruising flight is obtained from integration of the Specific Range. Several
flight schedules have been used for deriving the integrated range in Section 5.2. The
cruise/climb schedule is often used in theoretical work but the constant altitude and
Mach number flight is more practical from the operational point of view. An accurate
approximation for the range is
(7.1 )
where ( denotes the fuel fraction WF/Wi and kR is defined by Equation 6.1 for several
flight schedules. Usually kR is close to one, but lower values apply to constrained flight
conditions.
Calculation of the Fuel Load: using the approach presented in Sections 6.1 and 6.2
the mission and reserve fuel loads can be computed by means of closed form analytical
expressions, avoiding the usual labourious iterative procedure. The essential input required is the initial cruise Range Parameter Pi, which can be selected on the basis of the
described optimization process. If detailed data of the design are lacking, the Appendix
proposes a method for obtaining this factor from existing aircraft data, using an available
payload versus range diagram.

52

Appendix A
Statistical Derivation of the Range
Parameter
Section 4.2 was devoted to the analysis of flight conditions resulting in the maximum value
of the Range Parameter. Unconstrained as weU as constrained maxima have been treated
and the availability of drag polars and some basic engine data wi11 suffice to calculate
the (initial) Range Parameter required for the mission fuel computation. In some cases,
however, it may be useful to estimate Pi from statistical data for the Range Parameter
of existing aircraft derived from published payload versus range diagrams. The most
significant part of this diagram, see Figure A.l, is sector AB, corresponding to take off
with MTOW. At point A, referred to as the Harmonic or Nominal Range Rh, the slope of
this line is equal to minus the increase in cruise fuel required for a small range increment,
which is found from differentiation of Equation 5.14, for y = 1,
d(R IR )
h H

=P

d(WFIWt o )
'(I-O.5WF IWto F

(A.I)

Substitution of the fuel fraction according to Equation 6.l with kR

equation in P , for which an approximate solution is

= 1 yields aquadratic
(A.2)

where the normalized slope of the payload versus range diagram has been introduced,

~ ~f d( RI RH) ~
d(WFIWto )

Wt o ~ for R
RH -.6.Wp

= Rh

(A.3)

The identity .6.WF = -.6.Wp is strictly only valid when the reserve fuel is a constant
fraction of the MTOW; since this is not exactly true a prediction error of plus or minus
five percent may be introduced. The (negative) slope of the range versus the payload
(.6.RI.6. W p ) can be measured in the payload versus range diagram; see Figure A.I. Although Equation A.2 has been derived for flight at constant altitude and Mach number, it
is an equally good approximation when the cruise fuel is derived from the Brguet range
equation.

53

PAYLOAD

~R

~--------------------~--------------------~--RANGE

Rh

Figure A.l: Derivation of 'Pi from a payload versus range diagram

The method derived in this paragraph has been used to obtain a mean value of the
Range Parameter for a number of existing aircraft, for which a payload versus range
diagram was available to the author. The result is found in Tables 1 and 21 , from which
the following observations can be made:
1. Large long range turbofan powered transport aircraft achieve the highest range
factors: presently (1995) about 6.5 to 7.0 for modern, very long range transports
for LRC conditions, and about 10% less for HSC.

2. The Range Factor is aft"ected by the size of the aircraft: large aircraft have higher
values than small ones.
3. Short range aircraft are considerably less efficient in fuel usage due to their low
range factor. This is aresult mainly of their :fl.ight taking place at lower altitudes
and relatively short cruise sectors, contributing to a reduced mean value of the SR
for the overall :fl.ight.
4. In spite of the high propulsive efficiency of propellers the most fuel efficient propeller
aircraft reach their limit at about 4.5 for LRC conditions, and 15-20% below this
value for HSC.
1 Due to differences in reserve fuel policy, aircraft should not be compared on the basis of only these
derived data.

54

5. For several aircraft categories the effect of improvements in the state-of- the-art
can be clearly seen. For example, for the Boeing 707/320 which was designed
during the fifties, the range parameter was about 4.5 for LRC, while the more
recent Boeing 757 and 767 achieve 5.2 and 5.7, respectively. The present day, very
long range transport aircraft establish the practical limit of current state-of-the art:
a range factor of about 6.5 to 7, complying with a lift/drag ratio of about 20 and
an overall powerplant efficiency near to 35%.

55

Aircraft type
Airbus A 300-B2
Airbus A 310-220
Airbus A 320-200
Airbus A 330
Airbus A 340-200
Airbus A 340-300
Boeing 707-320
Boeing 727-200
Boeing 737-300
Boeing 747-300
Boeing 747-400
Boeing 757
Boeing 767-200
Boeing 777
Br.Aerospace 146-100
Br.Aerospace 146-200
BAC 1-11-400
BAC 1-11-500
Canadair Challenger
Fokker F-28/2000
Fokker 100
Lockheed L-1011/100
Lockheed L-1011/500
McD.Douglas DC-8/61
McD.Douglas DC-8/63
McD.Douglas DC-9/50
McD.Douglas DC-9/80
McD.Douglas DC-10/10
McD.Douglas DC-10/30
VFW-614

MTOW
(lb)
302,030
291,000
158,510
467,400
566,600
566,600
312,000
169,000
124,500
833,000
850,000
220,000
335,000
535,000
84,000
93,000
88,500
99 ,650
43,100
65,000
91,500
426,000
496,000
325,000
350,000
121,000
140,000
440,000
572,000
44,000

Rh
(nm)
860
1,245
1,265
3,233
5,475
4,965
4,280
948
952
4,750
5,600
1,100
2,930
3,100
925
1,175
980
880
2,557
440
780
2,500
3,660
2,550
3,400
765
820
2,330
4,167
658

t>.WF/t>.R
(lb/nm)
25.82
21.48
11 .80
23.76
25.54
26.66
19.61
15.13
10.37
48.39
43.97
16.41
19.75
26.47
9.65
11.83
8.84
10.49
3.45
12.71
9.08
28.92
30.03
25.40
23.32
11.53
14.07
30.07
34.98
7.81

(T/L/D)i
(-)

Remarks

4.55
5.16
5.10
6.78
6.65
6.57
4.55
4.29
4.64
4.85
5.28
5.16
5.78
7.12
3.26
2.77
3.78
3.62
4.04
1.96
3.90
5.05
5.21
4.18
4.68
4.09
3.83
5.09
4.83
2.08

HSC: 3.56
A 320-100: 4.60
GE CF-6-81 eng.

LRC

99% max.range
CF-34, LRC M=0.74
HSC
LRC
HSC: 4.96

Table 1: The derived Range Parameter for turbofan powered aircraft

56

Aircraft type
ATR-42/100
ATR-42/200
ATR-72
B.Aerosp. ATP
B.Aerosp. Jetstream
D.H.Canada Dash-7
D.H.Canada Dash-8
Dornier Do-228/200
Embraer Brasilia
Fokker F-27-Mk500
Fokker 50 Basic
Fokker 50 incr.MTOW
Hawker Siddeley 748
Piaggio P-180
SAAB-Fairchild 340
Fairchild Metro III

MTOW
(1b)
32,440
34,280
44,070
49,500
15,322
44,000
34,500
12,500
24,250
45,900
41,865
45,900
46,500
10,150
27,275
16,000

Rh
(nm)
776
910
0
560
600
678
550
580
260
600
100
900
755
800
300
760

-6.WF /6.R
(lb/nm)
3.00
2.97
4.47
4.95
2.11
5.79
3.33
2.33
2.98
5.64
3.98
4.26
4.95
1.03
3.35
1.79

(ryL/D)i
(-)

Remarks

4.22
4.46
4.15
3.97
2.93
2.90
4.12
2.00
3.32
3.16
4.39
4.14
3.42
3.93
3.30
3.43

HSC: 3.58
HSC: 3.82
LRC
LRC
HSC: 2.60
HSC: 3.74

LRC
HSC: 3.47
LRC
320 kts
LRC @ 25000 ft

Table 2: The derived Range Parameter for turboprop aircraft

57

Bibliography
[1] Anonymus, "Non-dimensional Graphical Method for the Analysis of Measurements
of Steady Level Speed, Range and Endurance: Aircraft with Turbo-jet and Turbo-fan
Engines," ESDU Item No. 70022, October 1973, Performance Volume 10.

[2] Anonymus, "Approximate Methods for Estimation of Cruise Range and Endurance:
Aeroplanes with Turbo-jet and Turbo-fan Engines," ESDU Item No. 73019, May
1982, Performance Volume 6.

[3] Anonymus, "Introduction to Estimation of Range and Endurance," ESDU Item No.
73018, October 1980, Performance Volume 6.

[4] Ashkenas, I.L., "Range Performance of Turbojet Airplanes," Journalof the Aeronautical Sciences, Vol. 15, Feb. 1948, pp. 97-101.

[5] Backhaus, G., "Grundbeziehungen fur den Entwurf optimaler Verkehrsflugzeugen,"

Jahrbuch der WGLR, 1958, pp. 201-213.

[6] Bert, C.W., "Prediction of Range and Endurance of Jet Aircraft at Constant Altitude," Journalof Aircraft, Vol. 18, 1981, pp. 890-892.
[7] Dommasch, D.O., Sherby, S.S., and Connolly, T.F., Airplane Aerodynamics, Pitman,
New Vork, 1957, pp. 286-287.
[8] Edwards, A.D., "Performance Estimation of Civil Jet Aircraft," Aircraft Engineering,
Vol. 22, No. 254, April 1950, pp. 94-99.
[9] Hale, F.J., "Best-Range Flight Conditions for Cruise-Climb Flight of a Jet Aircraft,"
NASA Conference Proceedings No 2001, Vol. 4,1976, pp. 1713-1719.
[10] Hale, F.J. and Steiger, A.R., "Effects of Wind on Aircraft Performance," Journalof
Aircraft, Vol. 16, June 1979, pp. 382-387.

[ll] Jonas, J. , "Jet Airplane Range Considerations," Journal ofthe Aeronautical Sciences,
Vol. 14, Feb. 1947, pp. 124-128.
[12] Mair, W.A. and Birdsall, D.L., Aircraft Performance, Cambridge Aerospace Series,
Cambridge University Press 1992, ISBN 0-521-36264-4.

58

[13J Milier, L.E., "Optimal Cruise Performance," Joumal of Aircraft, Vol. 30, No. 3,
May-June 1993, pp. 403-405.
[14] Martinez-Val, R. and Prez, E., "Optimum Cruise Lift Coefficient for Constant Altitude Constraint in Initial Design of Aircraft," Joumal of Aircraft, Vol. 29, No. 4,
1992, pp. 712-714.
[15J Martinez-Val, R. and Prez, E., "Extended Range Operations of Two and Three
Turbofan Engined Airplanes," Joumal of Aircraft, Vol. 30, No. 3, May-June 1993,
pp. 382-386.
[16] Menon, P.K.A., "A Study of Aircraft Cruise," AIAA Paper No. 86-2286 (AIAA
Conference Proceedings CP-8687, 1986).
[17] Mieie, A., Flight Mechanics, Volume 1: Theory of Flight Paths, Addison Wesley
Publishing Inc., Reading, Massachusets, 1962.
[18] Page, R.K., "Performance Calculation for Jet-Propelled Aircraft," Joumal of the
Royal Aeronautical Society, Vol.51, 1947, pp. 440-450.
[19] Peckham, D.H., "Range Performance in Cruising Flight," Royal Aircraft Establishment Technical Report 73164, March 1974.
[20] Perkins, C.D. and Hage, R.E., Airplane Performance, Stability and Control, Wiley,
New York, 1949, pp. 183-194.
[21] Ruij grok , G.J., Elements of Airplane Performance, Delft University Press, 1990,
ISBN 90-6275-608-5.
[22] Sachs, G. and Christodoulou, T., "Reducing Fuel Consumption of Subsonic Aircraft
by Optimal Cyclic Cruise," Joumal of Aircraft, Vol. 24, 1987, pp. 616-622.
[23] Schairer, G.S. and Olason, M.L. , "Some Performance Considerations of a Jet Transport Airplane," First Turbine Powered Air Transportation Meeting, Institute of the
Aeronautical Sciences, Seattle, Washington, August 9-11, 1954.
[24] Shevell, R.S., Fundamentals of Flight, Prentice-Hall, 1983, ISBN 0-13-339093-4.
[25] Torenbeek, E., "Fundamentals of Conceptual Design Optimization of Subsonic Transport Aircraft," Dept. of Aerospace Engineering, Delft University of Technology, Report LR-292, 1980 (Abstract in STAR N81- 13951).
[26] Torenbeek, E., "Some Fundamental Aspects of Transport Aircraft Conceptual Design," AGARD Conference Proceedings, No. 280, Paper 5, 1979.
[27] Torenbeek, E. and Wittenberg, H., "Generalized Maximum Specific Range Performance," Joumal of Aircraft, Vol. 20, No. 7, July 1983, pp. 617-622.

59

[28] Torenbeek, E., "The Initial Calculation of Range and Mission Fuel during Conceptual
Design," Dept. of Aerospace Engineering, Delft University of Technology, Report LR525, 1987.

[29] Visser, H.G., "Terminal Area Traffic Management," ProgTess in Aerospace Sciences,
Vol. 28, pp. 323-368, 1991.

60

Series 01: Aerodynamics

01.

02 .

03.

04.

05.

06.

07.
08.

09.

10.

11 .

12.

13.

F. Motallebi, 'Prediction of Mean Flow Data for Adiabatic 2-0 Compressible

Turbulent Boundary Layers'
1997 / VI + 90 pages / ISBN 90-407-1564-5
P.E. Skre, 'Flow Measurements for an Afterbody in a Vertical Wind
Tunnel'
1997 / XIV + 98 pages / ISBN 90-407-1565-3
B.W. van Oudheusden, 'Investigation of Large-Amplitude 1-DOF Rotational
Galloping'
1998 / IV + 100 pages / ISBN 90-407-1566-1
E.M. Houtman / W .J . Bannink / B.H. Timmerman, 'Experimental and
Computational Study of a Blunt Cylinder-Flare Model in High Supersonic
Flow'
1998 / VIII + 40 pages / ISBN 90-407-1567-X
G.J.D . Zondervan, 'A Review of Propeller Modelling Techniques Based on
Euler Methods'
1998 / IV + 84 pages / ISBN 90-407-1568-8
M.J. Tummers / D.M. Passchier, 'Spectral Analysis of Individual Realization
LDA Data'
1998 / VIII + 36 pages / ISBN 90-407-1569-6
P.J.J. Moeleker, 'Unear Temporal Stability Analysis'
1998 / VI + 74 pages / ISBN 90-407-1570-X
B.W. van Oudheusden, 'Galloping Behaviour of an Aeroelastic Oscillator
with Two Degrees of Freedom '
1998/ IV + 128 pages / ISBN 90-407-1571-8
R. Mayer, 'Orientation on Ouantitative IR-thermografy in Wall-shear Stress
Measurements'
1998 / XII + 108 pages / ISBN 90-407-1572-6
K.J.A. Westin / R.A.W.M. Henkes, 'Prediction of Bypass Transition with
Differential Reynolds Stress Modeis'
1998/ VI + 78 pages / ISBN 90-407-1573-4
J.L.M. Nijholt, 'Design of a Michelson Interferometer for Ouantitative
Refraction Index Profile Measurements'
1998/ 60 pages / ISBN 90-407-1574-2
R.A.W.M. Henkes / J.L. van Ingen, 'Overview of Stability and Transition in
External Aerodynamics'
1998/ IV + 48 pages / ISBN 90-407-1575-0
R.A.W.M. Henkes, 'Overview of Turbulence Models for External Aerodynamics'
1998/ IV + 40 pages / ISBN 90-407-1576-9

Series 02: Flight Mechanics

01.

02.

03.

E. Obert, 'A Method for the Determination of the Effect of Propeller Slipstream on a Statie Longitudinal Stability and Control of Multi-engined
Aircraft'
1997/ IV + 276 pages / ISBN 90-407-1577-7
C. Bill / F. van Dalen / A. Rothwell, 'Aircraft Design and Analysis System
(ADAS)'
1997 / X + 222 pages / ISBN 90-407-1578-5
E. Torenbeek, 'Optimum Cruise Performance of Subsonic Transport Aircraft'
1998 / X + 66 pages / ISBN 90-407-1579-3

Series 03: Control and Simulation

01.

02.

03.

J.C. Gibson, 'The Definition, Understanding and Design of Aircraft Handling

Qualities'
1997 / X + 162 pages / ISBN 90-407-1580-7
E.A. Lomonova, 'A System Look at Electromechanical Actuation tor Primary
Flight Control'
1997 / XIV + 110 pages / ISBN 90-407-1581-5
C.A.A.M. van der Linden, 'DASMAT-Delft University Aircraft Simulation
Model and Analysis Tooi. A Matlab/Simulink Environment for Flight Dynami cs and Control Analysis'
1998 / XII + 220 pages / ISBN 90-407-1582-3

Series 05: Aerospace Structures and

Computional Mechanics
01.

02.

03.

04.

A.J. van Eekelen, 'Review and Selection of Methods for Structural Reliability Analysis'
1997 / XIV + 50 pages / ISBN 90-407-1583-1
M.E. Heerschap, 'User's Manual for the Computer Program Cufus. Quick
Design Procedure for a CUt-out in a FUSelage version 1 .0'
1997 / VIII + 144 pages / ISBN 90-407-1584-X
C. Wohlever, 'A Preliminary Evaluation of the B2000 Nonlinear Shell
Element Q8N.SM'
1998/ IV + 44 pages / ISBN ,90-407-1585-8
L. Gunawan, 'Imperfections Measurements of a Perfect Shell with Specially
Designed Equipment (UNIVIMP)
1998/ VIII + 52 pages / ISBN 90-407-1586-6

Series 07: Aerospace Materials

01.

02.

03.

04.

05.

06.
07.

08.

A. Vaek / J. Schijve, 'Residual Strenght of Cracked 7075 T6 AI-alloy

Sheets under High Loading Rates'
1997 / VI + 70 pages / ISBN 90-407-1587-4
I. Kunes, 'FEM Modelling of Elastoplastic Stress and Strain Field in Centrecracked Plate'
1997 / IV + 32 pages / ISBN 90-407-1588-2
K. Verolme, 'The Initial Buckling Behavior of Flat and Curved Fiber Metal
Laminate Panels'
1998 / VIII + 60 pages / ISBN 90-407-1589-0
P.W.C. Prov6 Kluit, 'A New Method of Impregnating PEl Sheets for the InSitu Foaming of Sandwiches'
1998/ IV + 28 pages / ISBN 90-407-1590-4
A. Vlot I T. Soerjanto / I. Yeri / J.A. Schelling, 'Residual Thermal Stresses
around Bonded Fibre Metal Laminate Repair Patches on an Aircraft Fuselage'
1998 I IV + 24 pages I ISBN 90-407-1591-2
A. Vlot, 'High Strain Rate Tests on Fibre Metal Laminates'
1998 / IV + 44 pages / ISBN 90-407-1592-0
S. Fawaz, 'Application of the Virtual Crack Closure Technique to Calculate
Stress Intensity Factors tor Through Cracks with an Oblique Elliptical Crack
Front'
1998 / VIII + 56 pages I ISBN 90-407-1593-9
J. Schijve, 'Fatigue Specimens tor Sheet and Plate Material'
1998 / VI + 18 pages / ISBN 90-407-1594-7

Series 08: Astrodynamics and Satellite Systems

01.
02.

03.

E. Mooij, 'The Motion ot a Vehicle in a Planetary Atmosphere'

1997 / XVI + 156 pages / ISBN 90-407-1595-5
G.A. Bartels, 'GPS-Antenna Phase Center Measurements Performed in an
Anechoic Chamber'
1997 / X + 70 pages / ISBN 90-407-1596-3
E. Mooij, 'Linear Quadratic Regulator Design tor an Unpowered, Winged Reentry Vehicle'
1998/ X + 154 pages / ISBN 90-407-1597-1

A unified analytical treatment of the cruise performance of subsonic

transport aircraft is derived, valid for arbitrary gas turbine
powerplant installations: turboprop, turbojet and turbofan powered
aircraft. Different from the classical treatment the present report
takes into account compressibility effects on the aerodynamic
characteristics. Analytical criteria are derived for optimum cruise lift
coefficient and Mach number, without and with constraints on the
altitude and/or the engine rating. A simple alternative to the Brquet
range equation is presented which applies to aircraft cruising at
constant altitude and a high Mach number. A practical non-iterative
procedure for calculating the mission and reserve fuel loads in the
conceptual design stage is proposed as a conclusion of this report.

ISBN 90-407-1579-3

9 799040 715791

DELFT AEROS~""~":..
ENGINEERING

& TECHNOLOGY