Energy: Cesare Tona, Paolo Antonio Raviolo, Luiz Felipe Pellegrini, Silvio de Oliveira Ju Nior
Energy: Cesare Tona, Paolo Antonio Raviolo, Luiz Felipe Pellegrini, Silvio de Oliveira Ju Nior
Energy: Cesare Tona, Paolo Antonio Raviolo, Luiz Felipe Pellegrini, Silvio de Oliveira Ju Nior
Energy
journal homepage: www.elsevier.com/locate/energy
a r t i c l e i n f o a b s t r a c t
0360-5442/$ – see front matter Ó 2009 Elsevier Ltd. All rights reserved.
doi:10.1016/j.energy.2009.06.052
C. Tona et al. / Energy 35 (2010) 952–959 953
requirements of systems, augmenting the energy efficiency of non- use of the work potential concept [12] to evaluate the performance
propulsive systems [6]. of propulsion systems.
Different studies involving commercial aircraft systems have Therefore, the main goals of this study are:
been developed to apply exergy analysis to evaluate the replace-
ment of the conventional architecture by ‘more-electric’ ones [7,8]. To characterize the power plant operation evaluating the
Since the system energy flows are generated in the power plant, input and output exergy rates and the useful flow cost over
and the useful impact is expressed in terms of fuel consumption a complete mission;
reduction, the analysis of the aircraft engine is of fundamental To define the performance of engine components in each
interest. flight phase.
The focus of this work is the development of a model to evaluate
the performance of the energy conversion processes inside the
engine using exergy and thermoeconomic analysis as assessment
tools. Previous studies have dealt with exergy analysis applied to 2. System description and modeling
scramjet propulsion systems [9] and turbofan with an after burner
[10]; however, these studies considered one or two operational The mission profile used in this study represents a simplified
points (at sea level or in cruise conditions). Roth [11] proposes the typical commercial flight as shown in Fig. 2. The operational points
AIRCRAFT
BFlight BControl Surfaces
AIRFRAME
BDest,Flight BThrust
ENGINE
Fig. 1. Demand and penalties imposed by Air Management Systems to engine and airframe.
954 C. Tona et al. / Energy 35 (2010) 952–959
considered are take-off, climbing, cruise, descent, holding and anti-ice-OFF: air is extracted from the 9th stage of the
landing. compressor.
To develop the engine numerical model, the GSP (Gas Turbine
Simulation Program) software [13] was used in order to calculate The design point for the engine was set on cruise phase: altitude
the fuel consumption, thrust and inlet air, bleed, fan air and outlet 37 000 (ft), Mach 0.77. All other phases were calculated considering
gases thermodynamic characteristics. The power plant is a classical an off-design simulation of the engine.
turbofan with separated flows, a high bypass ratio and two shafts. The numerical results were validated in the design point using
Fig. 3 shows a schematic representation of turbofan analyzed and real reference data and a mathematical model developed for this
implemented in GSP. engine with MatlabÒ [14] software. The average error based on real
Bleed air and extracted power controls were implemented to reference data is 3.3 %, and the maximum error is 8.2 %. Table 1
simulate the different energy demands of the board systems during shows some characteristics for each one of the flight phases.
the whole mission. Bleed controls were located on the fan and
compressor, and hydraulic and electric power extractions were 3. Exergy analysis
located on the turbine stages.
The compressor has two bleed air conditions in order to meet Two cases and models were considered in this study:
the demands of the anti-ice system:
i. Global model: the power plant is considered a unique open
anti-ice ON: air is extracted from the 14th stage of the system; this model allows the evaluation of the exergy
compressor; destruction and loss rates for the different flight conditions.
Fan Bleed
Flow
4
2 3
MIXER
LOW PRESSURE
10
TURBINE
1 11 12
COLECTOR NOZZLE
FAN
5 9
HIGH PRESSURE
COMPRESSOR
TURBINE
7 8
6
Compressor
Bleed Flow
Fuel
Table 1 were used for both references. Fig. 4 presents the control volume
Flight phases characteristics with anti-ice OFF. considered in the global model for a turbofan engine.
Flight phase Height (m) Mach number Time duration (min) Thrust (N) The exergy balance for this control volume is shown in Eq. (1):
Take-off 0 0.20 1 28,268
Climb 6,096 0.57 12 13,652 B1 þ Bfuel ¼ B12;lost þ Wex þ Bcomp bleed þ Bfan bleed þ BT
Climb 10,668 0.66 4 8,634
Cruise 11,277 0.77 34 7,713
þ Bdestroyed (1)
Descent 10,668 0.77 4 2,478
Descent 6,096 0.62 11 2,509
Table 2 presents exergy rates for the different flows during
Holding 4,572 0.40 10 13,193 cruise phase for both reference systems. It can be seen that taking
Landing 0 0.20 0.5 2,095 the reference state on earth gives a negative value for the inlet air
exergy, due to the smaller pressure of this flow compared to the
reference pressure.
ii. Local model: the power plant is studied analyzing its equip- Furthermore, it is important to mention that, when the refer-
ment one by one, using local balances for each component in ence is attached to the engine, the exergy of the thrust is related to
order to evaluate its performance and thermoeconomic cost. the kinetic exergy of the exhaust gases, i.e., the kinetic exergy of the
exhaust gases is equal to the exergy of the thrust plus an exergy
To direct the study, the following assumptions were made: destruction rate associated to the energy conversion process and an
exergy loss rate associated to the exhaust gases. The term B12,lost in
(a) All gases were modelled as ideal gases; Eq. (1) represents this exergy loss rate. When the reference is fixed
(b) The fuel is kerosene (JET-A1), its chemical formulae is C12H23; on the ground, the term B12,lost accounts for the exergy of stream 12
(c) The fuel physical exergy is not considered; that is no further used.
(d) All engine equipment are considered adiabatic;
(e) Heat exchangers (air-fuel, fuel-oil) are not considered. 3.2. Local model
EESÒ [15] software was used to implement the analysis. As for the local model, exergy balances are developed for each
component, defining its efficiency and performance related to the
3.1. Global model power plant global values. Such analysis allows the identification of
the components that contribute most for the exergy destruction
The first step to evaluate the exergy content of an energy flow is inside the engine. Hence, these results might define optimization
the definition of a reference state. Generally this choice is easy, priorities among the system and direct investment to where solu-
because the system under evaluation is fixed upon the earth or its tions will be more cost-effective or exergy-efficient.
altitude is constant. Analyzing a flight vehicle or a jet engine, In this analysis, the reference state for the exergy calculations
operating in high altitude and different environmental conditions was set by the stagnation air temperature and pressure for each
during the flight, makes it more difficult to choose a unique refer- flight phase [16]. Again, as a first approach, standard chemical
ence [16]. Szargut’s environment gives negative values for some exergies given by Szargut et al [17] were used.
exergy flows in these conditions [17]. For this reason, preliminary
analysis were performed using two different reference systems, 3.3. Exergy analysis results
one fixed on the ground and one attached to the power plant. Thus,
when the reference is fixed on the ground, Szargut’s environment 3.3.1. Global balance
was considered, and, when attached to engine, the reference state Table 3 shows the distribution of output exergy flows for all the
was set by the stagnation outside air temperature and pressure for flight phases. The exergy efficiency is defined as the ratio of sum of
each flight phase. In addition, as a first approach to the exergy the thrust and bleed exergies (formed by air bleed and extracted
calculations, standard chemical exergies given by Szargut et al [17] power) to the fuel exergy:
B1 B12
Sea Level
Table 2
Table 4
Exergy rates for the different flows during cruise.
Exergy efficiency in each flight phase.
# Description Exergy flowrate (kW)
Phase Anti-Ice system Exergy efficiency (%)
Ground Engine 1 OFF 10.06
reference reference 2 OFF 20.71
1 Collector inlet 202.8 588.6 3 OFF 24.07
2 Fan inlet 792.0 423.6 4 OFF 26.46
3 Fan outlet to mixer 2114.6 1862.0 5 OFF 20.37
4 Fan bleed 42.3 37.3 6 OFF 13.94
5 Fan outlet to compressor 392.8 340.4 7 OFF 16.34
6 Compressor bleed 43.7 42.5 8 OFF 5.87
7 Compressor outlet 3456.3 3490.0 1 ON 10.55
8 High pressure turbine inlet 8303.6 8734.0 2 ON 21.01
9 Low pressure turbine inlet 4611.1 5055.0 6 ON 14.86
10 Low pressure turbine outlet 2421.0 2875.0 7 ON 17.44
11 Mixer outlet 2984.5 3477.0 8 ON 8.14
12 Nozzle outlet 349.4 3385.0
Fuel 7186.0 7168.8
Thrust 1745.5 1745.5
Fan power 2071.6 2071.6 the exergy loss associated to this flow (since it includes the
Compressor power 3546.6 3546.6 conversion of kinetic energy into thrust). On the other hand, using
Extracted power to hydraulic and electric 16.4 16.4 the reference fixed on the ground gives smaller values for the
systems
thermo–chemical exergy of the exhaust gases, leading to higher
values of exergy destruction. All in all, the relation of total irre-
P versibilities to the exergy input in both cases is very similar, inde-
Bi pendently of the reference chosen.
i To investigate the complete flight, an integral analysis was
hb ¼ P (2)
Bj considered. This study allows relating the destroyed exergy
j
(proportional to thrust and fuel flow) for each phase to its time
The results for the complete mission are presented in Table 4, duration. Thus, it is possible to identify the critical phases during
including anti-ice OFF and ON conditions, and the reference state the mission. These results are presented in Table 5, which shows
fixed on the ground. that cruise is the flight condition of maximum exergy destruction
Table 4 and also Fig. 2 show that: since it is the longest, followed by take-off (maximum thrust
required) and holding phases.
Cruise is the design and longest phase. It presents the
maximum exergy efficiency (26%), demonstrating that the 3.3.2. Local balance
engine design point is the most efficient one; For each flight phase, it was calculated the distribution of irre-
The exergy efficiency, lower than 10% during take-off, versibilities among components. In the nozzle, the exergy balance
increases during the mission reaching its maximum value in considered the variation between the inlet exergy (stream 11 in
cruise. After the cruise, the efficiency decreases reaching its Fig. 3) and the total outlet exergy (stream 12 in Fig. 3).
minimum value during landing (6%);
Exergy losses are a function of the aircraft speed and thrust. B11 ¼ B12 þ Bdestroyed (3)
This leads to an interesting result: a better efficiency is obtained Another exergy balance was set by the difference of stream 12
mainly by reducing exergy losses and not by decreasing the exergy and the thrust exergy to calculate the destroyed exergy
engine destroyed exergy. associated to the thrust generation from the kinetic exergy of
During landing phase with anti-ice system ON the exergetic stream 12 and the lost exergy associated to thermo–chemical and
bleed flow value is near to thrust exergetic value. the residual kinetic exergy (relative velocity of the exhaust gases) of
stream 12.
Considering the reference state attached to the engine, and that
the thermo–chemical and part of the kinetic exergy of the outlet B12 ¼ BT þ Bdestroyed þ B12;lost (4)
gases as an exergy loss, the results are very similar to the ones in
Tables 3 and 4. Hence, the use of different environment references In this way, it is possible to evaluate the exergy destruction of
does not give any misleading results regarding the engine effi- the energy conversion process associated to the thrust generation.
ciency. However, the change in the environment reference leads to Table 6 shows the total destroyed and lost exergy in each
a different distribution of the irreversibilities between destroyed component for the complete flight.
and lost exergy. Whenever the reference is set to the engine, the Components that destroy more exergy are those related to the
thermo–chemical exergy of the exhaust gases is high, increasing more irreversible processes, such as combustor and mixer. Their
Table 3
Exergy rate distribution in anti-ice OFF flight conditions.
Phase Destroyed exergy (%) Lost exergy (%) Thrust exergy (%) Bleed and extracted exergy (%)
1 64 26 9 1
2 68 11 20 1
3 71 5 23 1
4 69 5 25 1
7 71 13 15 1
8 83 11 4 2
C. Tona et al. / Energy 35 (2010) 952–959 957
Table 5 and flowrates of the engine was developed based on the results of
Total destroyed exergy in the various flight phases with anti-ice OFF. the exergy analysis. Two different analyses were performed:
Phase Total destroyed exergy (%)
1 4 Global power plant production cost evaluation, considering
2 25 only the inlet and outlet flows (air, fuel, exhaust gases, thrust
3 5 and bleed exergy rate and extracted power), as shown in Fig. 4;
4 36
Local power plant production cost evaluation, taking into
5 3
6 9 account each power plant component and its relative input and
7 19 output exergy rates and flowrates, as shown in Fig. 3.
8 <1
contribution to the total exergy destruction is between 50% and 66% 4.1. Global power plant production costs evaluation
depending on the flight phase. Also, the destroyed and lost exergy
associated to the generation of thrust is very significant. The sum of Following the same approach used in the exergy evaluation of
the combustor, mixer and thrust generation losses represents up to the power plant, two reference systems were adopted: one fixed
83% of the total irreversibilities, being 80% on average. Besides, upon earth and the other fixed to the engine.
Table 6 shows that the main distribution of irreversibilities does not For both reference systems, the cost balance is:
change considerably during the mission.
More over, the following remarks can be drawn: Cfan bleed þ Ccomp bleed þ Cex þ CT ¼ C1 þ Cfuel þ Cequip (5)
Engine destroyed exergy is inversely proportional to thrust No cost is assigned to stream 12 exergy, since part of it is con-
and represents between 70% and 80% of the total irreversibil- verted into thrust and part is either destroyed or lost in the exhaust,
ities. During phases of high and medium thrust, this value is as discussed previously. Also, Table 7 presents economic data used.
approximately 70%, while during low thrust phases it can reach Making use of the thermoeconomic specific cost (c ¼ C/B) [18],
85%; the equality criterion was used in this analysis, considering the
For the whole flight, the components that destroy more exergy same thermoeconomic specific costs for all the useful products
are the combustion chamber and the mixer; [18]:
The appraisal of optimization studies involving engine compo-
cT ¼ ccomp bleed ¼ cfan bleed ¼ cex ¼ c (6)
nents must consider technological and economics aspects. Thus,
not necessarily the most irreversible components are the better
focus of an optimization process. For instance, considering the Cequip þ cfuel Bfuel
c ¼
‘more-electric’ architectures it seems important to validate the BT þ Bcompr bleed þ Bfan bleed þ Wex
variation on the exergy performance of equipment involved on Cequip c
bleed air and extracted power: fan, compressor and turbines; ¼ þ fuel (7)
BT þ Bcompr bleed þ Bfan bleed þ Wex hGlobal
Fan and compressor destroy more exergy than turbines;
The lost exergy is directly proportional to thrust (or to relative Eq. (7) indicates that the cost of the products is a function of the
speed of exhaust gases). The higher the lost exergy, the higher capital cost and the efficiency of the engine associated with the fuel
the thrust is (or the relative speed of exhaust gases). cost. Thus, the minimization of the product cost is related to the
decrease of the capital cost or the increase of the exergy efficiency
Fig. 5 presents the Grassmann diagram in cruise phase, using the of the system. In the development of Eq. (7) it was considered that
reference attached to the engine to avoid negative values for stream 1 is costless (C1 ¼ c1 ¼ 0).
exergy. This diagram shows that 22.5% of the inlet exergy is used for
thrust and 1.3% as input to systems based on compressor air bleed.
4.2. Local power plant production costs evaluation
4. Thermoeconomic analysis
The local thermoeconomic analysis allows calculating both the
In order to determine the utilities exergy based production useful product costs and the internal exergy rate costs. For this
costs, a complete evaluation of the cost of the different exergy rates evaluation, the local model is the same used in the local exergy
analysis. Again, the equality criterion was used to attribute cost for
more than one stream leaving each control volume.
Table 6 Cost balances applied to components shown in Fig. 3 are as
Relative destroyed and lost exergy in anti-ice OFF conditions. follow:
Destroyed/Lost exergy (%) Flight phase
1 2 3 4 5 6 7 8 Collector:
Inlet 2 3 3 3 3 3 2 <1
Fan 5 5 4 4 6 6 5 3
Compressor 6 6 6 6 7 8 6 10
Fan duct 1 1 1 1 3 3 1 1
c1 B1 þ Ccol ¼ c2 B2 (8)
Combustor 32 33 31 33 34 37 33 44
HP turbine <1 <1 1 1 1 1 1 1 Fan:
LP turbine <1 <1 <1 <1 <1 <1 <1 <1
Turbine duct <1 <1 <1 <1 <1 <1 <1 <1
Mixer 19 20 19 20 20 21 19 22
c2 B2 þ Cfan þ Wfan cfan ¼ c3 B3 þ c4 B4 þ c5 B5 (9)
Nozzle <1 2 1 2 2 2 1 <1
Shaft friction 2 2 2 2 2 2 2 1
Thrust generation 32 27 31 28 21 16 28 15 c3 ¼ c4 ¼ c5 ¼ cfan bleed (10)
958 C. Tona et al. / Energy 35 (2010) 952–959
Table 7 Mixer
Data used for the exergy based costs analysis.
Fuel specific initial cost (US$/GJ) 9.022 c10 B10 þ c3 B3 þ Cmix ¼ c11 B11 (18)
Power plant economic life (years) 25
Annual utilization factor 0.285
Total capital cost (US$) 2,000,000 Nozzle
Compressor:
c11 B11 þ Cnoz ¼ cT BT (19)
c5 B5 þ Ccomp þ Wcomp ccomp ¼ c6 B6 þ c7 B7 (11)
The solution of this set of equations gives the values of each one
of the production costs of the utilities of the power plant.
c6 ¼ c7 ¼ ccomp bleed (12)
4.3. Fuel and power plant components costs
ccomp ¼ c9 (15)
Table 8
Components initial costs.
Table 9 presents the average annual production costs of the power References
plant products for both global evaluation approaches. It is interesting
to verify that the obtained results for electrical, hydraulic and thrust [1] Muñoz JD. Optimization Strategies for the Synthesis/Design of Highly Coupled,
are quite close, since the change in the reference system does not alter Highly Dynamic Energy Systems. Ph.D Dissertation, Faculty of Virginia Poly-
technic Institute and State University, Blacksburg, VA, 2000.
the exergy efficiency of the engine, which is directly related to the cost [2] Moorhouse DJ. Proposed system-level multidisciplinary analysis technique
[18]. On the other hand, exergy flowrates are more sensitive to based on exergy methods. Journal of Aircraft 2003;40(1):11–5.
changes in the reference systems (Table 2). [3] Bejan A, Siems DL. The need for exergy analysis and thermodynamic optimi-
zation in aircraft development. Exergy 2001;1(1):14–24.
Table 10 presents the results obtained with the local evaluation [4] Paulus D, Gaggioli R. Rational objective functions for vehicles. AIAA Paper No.
approach. In this case, as the cost formation process is more accurately 2000-4852, 8th AIAA/USAF/NASA/ISSMO Symposium on multidisciplinary
followed by cost balance equations applied to components of the analysis and optimization, Long Beach, California, September 6–8, 2000, p. 1–11.
[5] Gandolfi R, Pellegrini LF, Silva G, Oliveira Jr. S. Exergy analysis applied to
power plant, the results are more representative than those shown in
a complete flight mission of commercial aircraft. 46th AIAA Aerospace Science
Table 9 for the same products. Besides that, the set of cost balance Meeting and Exhibit 7–10 January 2008, Reno, Nevada.
equations allows one to better understand the energy conversion [6] Ensign TR. Performance and weight impact of electric enviromental control
system and more electric engine on citation CJ2. 45th AIAA Aerospace Science
processes.
Meeting and Exhibit 9–11 January 2007, Reno, Nevada.
Comparing the results of Table 9 to those of Table 10, it is possible [7] Pellegrini LF, Gandolfi R, Silva GA, Oliveira, Jr. S. Exergy analysis as a tool for
to see that the global balance overestimates (more than doubles) the decision making system design. 45th AIAA Aerospace Science Meeting and
cost of bleed air from the fan and the compressor, although the cost of Exhibit 2007, January 8–11, Reno, Nevada. (in CD-ROM).
[8] Gandolfi R, Pellegrini LF, Silva GA, Oliveira Jr. S. Aircraft Air management
thrust is only 3.2% lower than that of local model. This is a direct systems trade-off study using exergy analysis as a design comparison tool.
consequence of the contribution of the thrust exergy to the total 19th International Congress of Mechanical Engineering, November 5–9, 2007,
useful exergy from the engine (95%). Brasilia, DF, Brazil.
[9] Amati V, Bruno C, Simone D, Sciubba E. Development of a novel modular
simulation tool for the exergy analysis of a scramjet engine at cruise condition.
International Journal of Thermodynamics 2006;9(4):1–11.
5. Conclusions [10] Turgut ET, Karakoc TH, Hepbasli A. Exergetic analysis of an aircraft turbofan
engine. International Journal of Energy Research 2007;31(14):1383–97.
[11] Roth BA. A worktransfer perspective of propulsion system performance. 40th
This work presented an application of exergy as a tool to evaluate AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, July 11–14,
the performance of an aircraft power plant. 2004, Fort Lauderdale, Florida, USA.
For the global analysis, two reference systems were proposed: one [12] Roth BA. McDonald R, Mavris D. A method for thermodynamic work potential
analysis of aircraft engines. 38th AIAA/ASME/SAE/ASEE Joint Propulsion
fixed on the ground and one attached to the engine. In the first case, Conference and Exhibit, July 7–10, 2002, Indianapolis, Indiana.
negative values for exergy were obtained, as reported also in [17], [13] GSP Development Team. GSP 10 user manual. Amsterdam, The Netherlands:
while for the second case it does not happen as the reference varies National Aerospace Laboratory NRL; 2004.
[14] MATLAB, The language of technical computing. MathWorks, Inc.; 1994–2007.
according to the stagnation temperature and pressure at each flight
[15] Klein SA, Alvarado FL, EES – Engineering equation solver for microsoft
phase. Although, as a first approach, the chemical exergies used were windows operating systems. F-Chart Software, 2007.
those reported in [17], the use of a variable reference seems to be [16] Göğüs YA, Çamdali U, Kavsaoğlu MS. Exergy balance of a general system with
variation of environmental conditions and some applications. Energy
thermodynamically better.
2002;27:625–46.
In the global analysis, the overall engine performance was evalu- [17] Szargut J, Morris DR, Steward FR. Exergy analysis of thermal, chemical and
ated over a complete flight mission, characterizing the destroyed, lost metallurgical processes. New York: Hemisphere Publishing Corporation; 1988.
and useful exergy trends. The maximum exergy efficiency was found [18] Gaggioli RA, Wepfer WJ. Exergy economics: I. cost accounting applications,
ii. benefit-cost conservation. Energy 1980;5:823–37.
during cruise phase, 26.5%, but this value decreased to 6% during [19] Raymer DP. Aircraft design: a conceptual approach. Washington: American
landing, considering the reference fixed on the ground. With the Institute of Aeronautics and Astronautics; 1992.