A Comparison of Parameter Variation With Numerical Optimization
A Comparison of Parameter Variation With Numerical Optimization
1111111111111111111111
BREAK
Joachim Kurzke
Engine Performance Department
Mill Munchen GmbH
MOnchen, Germany
ABSTRACT routines. This can be done for example with the Windows-based gas turbine
In gas turbine performance simulations often the question arises: What performance program GasTurb developed by the author.
is the best thermodynamic cycle design point? This is an optimization task Certainly it is not sufficient, to get a single cycle as the best solution
which can be attacked in two ways: One can do a series of parameter from the computer program. There arc two main reasons for that:
variations and pick from the resulting graphs the best solution or one can Firstly, the numerical optimization algorithm will find the optimum of
employ numerical optimization algorithms that produce a single cycle which the mathematical model only as opposed to the "true" optimum. An exotic
fulfills all constraints. cycle as result of an optimization run is mostly a hint to a deficiency of the
The conventional parameter study builds strongly on the engineering model. In such a case most probably a design constraint has been
judgement and gives useful information over a range of parameter overlooked when defining the problem.
selections. However, when values for more than a few variables have to be Secondly, it is always of interest to know about the neighborhood of the
determined while several constraints are existing, then numerical optimum solution. From a parameter study limited to the region of interest
optimization routines can help to find the mathematical optimum faster and it becomes obvious which design variables and constraints have the biggest
more accurately. Sometimes even an outstanding solution is found which impact on the result.
was overlooked while doing a preliminary parameter study. One of the advantages of numerical optimization is, that the region
For any simulation task a sophisticated graphical user interface is of where parameter studies should be performed is narrowed down
great benefit This is especially true for automated numerical optimizations. significantly.
It is quite helpful to see on the screen of a PC how the variables are
changing and which constraints are limiting the design. A quick and clear 2. PARAMETER STUDIES
graphical representation of trade studies is also of great advantage. The The results of parameter studies are normally presented as graphics. In
paper describes how numerical optimization and parameter studies are a single graph one can show for given ranges of two parameters the results
implemented in a Windows - based PC program. for several dependant quantities. In a cycle study for a single spool turbojet
Man example, the cycle selection of a derivative turbofan engine with engine, for example, with compressor pressure ratio and burner exit
a given core shows the merits of numerical optimization. The parameter temperature as design parameters, one can plot the specific fuel
variation is best suited for presenting the sensitivity of the result in the consumption over specific thrust. In the resulting carpet one can additionally
neighborhood of the optimum cycle design point. show lines for other calculated parameters as for example the turbine
pressure ratio and the turbine exit temperature, see figure I.
1. INTRODUCTION Let us assume, for example, that the design aim is a low cost turbojet
The traditional way to select the thermodynamic cycle of a new gas with a single stage turbine and an uncooled turbine exit casing. The feasible
turbine employs extensive parameter variations. For a complex engine with region of design parameter combinations can be marked easily in the carpet
many design variables this is a time-consuming task. One looks for the since it is limited by the following constraints:
optimum solution in a certain respect.
Instead of screening a wide range for the design variables with turbine pressure ratio < 3.5 single stage turbine
systematic parameter variations it is also possible to do an automatic search turbine exit temperature < 1200K uncooled turbine exit casing
for the optimum engine design with the help of numerical optimization
Presented at the International Gas Turbine & Aeroengine Congress & Exhibition
Stockholm, Sweden — June 2—June 6, 1998
This paper has been accepted for publication in the Transactions of the ASME
Discussion of it will be accepted at ASME Headquarters until September 30, 1998
Downloaded From: http://proceedings.asmedigitalcollection.asme.org/ on 03/04/2016 Terms of Use: http://www.asme.org/about-asme/terms-of-use
From figure 1 one can read that with these constraints the cycle with maximum of the figure of merit) faster and more accurately. As shown
the highest specific thrust has a burner exit temperature of 1600K and a above with the turbojet example, in a parameter study with only two
compressor pressure ratio of around 16.5. variables it is easy to find an optimum solution. If there are three variables
the situation is not so clear. With more than three variables the picture may
get obscure. In complex studies the true optimum may never be found with
Compr.
32.5 I M° 7.5 Pressure the conventional parameter study.
Ratio There are many numerical optimization algorithms known from
cotale 1603
e
10 literature. They can be divided basically into two major groups: Methods
— 30
es 03CtISCO 12.5 that use gradient information and others. In the program GasTurb there is
JC
is one method from each group implemented. A short explanation how these
-"
s- 275
-
1Y-lt‘uco
algorithms work is given in the following chapters.
0 No)
25
frstolellfr
Turbine 3.1 Gradient Method
22.5 Pressure moo nl wre Them is a good example for the optimization task: A mountaineer shall
3 Ratio climb the highest peak in a certain region. He has no map and the weather
7
or. 20 005`ne °de is foggy. His only tool is an altimeter. What is he going to do? He will
to
certainly check his surroundings first and then go in the direction of the
17.5 steepest ascent In the end he will come to the top of a mountain. This is a
12 14 16 18 20 place where each step leads downwards.
Net Thrust NM
The steepest ascent may, however, lead toward a border (which is either
Figure 1: Result of a parameter study the lower or upper limit of a design variable) of the region. Then our
mountaineer will walk along the border until he reaches the place where
each step leads downwards or out of the allowed region.
However, most gas turbine cycles are much more complex than the Is that the end of the story? Not necessarily. There might be several
turbojet example discussed above. A two-spool turbofan has the following summits within the region. Our mountaineer may have found the highest
five cycle design variables: peak by chance, but he cannot be sure of that. He has to check other parts
of the region. In mathematical terms there might be "local" optimums
low spool pressure ratio besides the "global" optimum.
high spool pressure ratio Up to now we have not spoken of constraints. They are like fences. A
burner exit temperature part of the region is forbidden to our mountaineer. His task is made more
bypass ratio difficult because on his way to the summit he may have to walk downwards
fan pressure ratio for a while to avoid a forbidden region. The fences (the constraints) often
exclude the summit (where each step leads downwards) as an acceptable
There will also be more constraints than with the turbojet example: solution. They create local optima that would not exist without fences.
Constraints make the task of optimization difficult.
low pressure turbine inlet temperature < limit Let us turn to the mathematical algorithm now. The mountaineer who
uncooted low pressure turbine first makes test steps in several directions uses the "gradient strategy" as a
compressor exit temperature < limit search method. With the test steps he is looking for the partial derivatives
material limit for compressor disk avav, For each optimization variable he must do one test step before he
fan pressure ratio < limit can start his way in the "right" direction.
single stage fan After the first step uphill the local gradients will be different. The test
high pressure turbine pressure ratio < limit steps could now be repeated to find the new direction. Test steps take time,
single stage turbine however, and it is therefore better to go on in the same direction as long as
the altitude increases. Reaching a fence (violating a constraint) could be
In a more detailed study there will be even more design variables as for another reason for stopping the climb. Only then will new gradients be
example the stage numbers for the high and the low pressure turbine. It is sought. The new direction will eventually take you along a fence.
obvious that with a parameter study it will be very difficult and time The gradient search algorithm implemented in GasTurb was derived
consuming to find the optimum values for the design variables. from ref. Pl. The principle is the following (see figure 2). We begin at the
point marked "Start I", looking for the direction of the steepest gradient
3. NUMERICALOPTIMIZATION ("Direction 1"). Following this direction we walk to the highest point. Then
By the way how is the optimum defined in a mathematical sense? In a we change the direction by 90° (orthogonal). This can be done without
parameter study that question must not be answered a priori. In a numerical evaluating the local gradient. We again go for the highest point here. To
optimization, however, a figure of merit must be clearly defined before the define the third direction we use the experience from the first two directions.
calculation can commence. The figure of merit might be the specific fuel We connect the point "Start 1" with the optimum point found along
consumption of a turbofan at cruise which is to be minimized. For a fighter "Direction 2". We follow this direction again as long as altitude increases.
engine it might be that the specific thrust shall be maximized. One can also This procedure can be applied repeatedly until the search steps or the
think of a weighted combination of these parameters. changes in the "figure of merit" become very small. There is also a
When values for more than a few variables have to be determined while maximum limit for the number of optimization steps. In the example of the
several constraints are existing, then numerical optimization routines can figure the optimum is found along search direction 6 (not marked in the
help to find the mathematical optimum (i.e., the minimum respectively figure, perpendicular to direction 5).
2
Downloaded From: http://proceedings.asmedigitalcollection.asme.org/ on 03/04/2016 Terms of Use: http://www.asme.org/about-asme/terms-of-use
The dashed line in the figure shows how optimization would go on, if 4. CYCLE SELECTION FOR A DERIVATIVE TURBOFAN
only local gradient information is used. With this simple strategy, the search A very common design task is to adapt an existing engine for a new
direction would change very often. application. It is quite obvious that in this case there are more constraints
than during the design of a brand new engine. In this chapter at first the
basic engine will be described and then the design variables, the constraints
and the figure of merit for the numerical optimization of a derivative engine.
Stitt 2
Thrust 3.64 kN
Up to now we have only dealt with optimization without constraints. In Burner Exit Temperature 1350 K
the figure there is a shaded zone which suggests a forbidden region. If we
use the strategy just described the search for the optimum will end at the
Overall Pressure Ratio 17.82
point "A" along "Direction 5". We cannot find the global optimum if we
Core Pressure Ratio 12
begin at "Start 1". If we begin at "Start 2", however, we will be at the top of
the hill very quickly. ISA Corrected Mass Flow 60 kg/s
with
new value for optimization variable
value of V; producing the best figure of merit
search region for variable V,
range reduction coefficient (positive integer)
distribution coefficient (positive odd integer)
random number between zero and one
Figure 5: Graphical user interface of GasTurb Bypass 4.5 4.65 5.06 5.23
Ratio
Lo
5
theoretical considerations one can derive, that this ratio should be equal to cc
the product of fan and low pressure turbine efficiency when an unmixed 5
flow turbofan is to be optimized for SEC. Note that the numerical °
al t Thalicksi,
optimization algorithm has automatically found a cycle for which the jet
WI
C.44 8141
1kkS\
%%.141 4‘
11%1114911111111
; 8.75
velocity ratio is near to its theoretically best value. TOC.
4
00 *frit: 4407
470,0411,441
5sF
1925
19
25
Ara
10
3.5 3.75 4 4.*et 4.5 t p475 5 5.25 5.5
u• 17.5'
ie 15
wx 19.6
12.5 38 re1.4 4728a4141
3 1 02
10
3.8 4 4.2 4.4Net Ttlitt 4.8 5 5.2 5.4 t 19.4 5-
all/I
111111■119 ISO
Figure 6: Sensitivity for Booster and Outer Fan Pressure Ratio 0
5 19.2
annm.
WlNMTryj 45 g
parameter on the vertical axis was mainly selected because it spreads the Trust 4.5 IN
carpet nicely. The corrected flow at the low pressure turbine exit is not a 18.84 35
4.4 4.5 4.55
very important engine design parameter. Net4.45
Thrust NM
Actually the limiting line for the fan diameter (0.75m) is identical to the Figure 7: Sensitivity for HPC Map Operating Point
line for the minimum high pressure turbine flow capacity. To the left of the
line for booster pressure ratio 1.79 the lower limit for the high pressure
turbine flow capacity is violated and to the right of this line the fan diameter As an example, the cycle selection of a derivative turbofan engine with
a given core shows the merits of numerical optimization. The parameter
is too big. In this graph only one point is a valid solution: the square which
variation is best suited for presenting the sensitivity of the result in the
fulfills the thrust requirement. One can see, that the optimum solution is
neighborhood of the optimum cycle design point. Sometimes this leads to
pretty much boxed in by the design constraints.
a redefinition of the figure of merit or the constraints imposed on the
In figure 7 two more design variables were systematically varied around
solution. In rare cases even an outstanding solution is found which was
the optimum solution. Again the square marks the only point which fulfills
overlooked while doing a preliminary parameter study.
all constraints.
Figure 8 deals with the operating point of the high pressure compressor.
High values for the auxiliary coordinate beta go with a low surge margin.
The square marks the point with lowest SFC in the region of feasible
designs.