Document Type : Research Article


1 Division of Thermal Science & Energy Systems, Department of Mechanical Engineering, Faculty of Technology & Engineering, University of Qom, Qom,

2 Division of Thermal Sciences & Energy Systems, Department of Mechanical Engineering, Faculty ‎of Technology & Engineering, University of Qom, Qom, Iran


In this paper, Multi-Objective Optimization of a 160 MW combined cycle power plant has been performed from Exergetic, Economic, Environmental aspects simultaneously. In this regard, thermodynamic, exergetic and thermoeconomic modeling and simulation for this case have been done. The Multi Objective Evolutionary Algorithms (MOEAs) are used and evaluated by MINLP approach to find optimum design in view of exergetic, exergoeconomic and environmental impacts. In order to facilitate selecting optimum sets, the environmental impact objective function has been defined and expressed in cost terms and added to the economic objective. Furthermore, extended combined pinch and exergy analysis has been performed to demonstrate the system graphically for each case in the base case and the optimum case. So, the performance of a different component of the system can be demonstrated better. In addition, the feasible region for optimization problem has been indicated by extended combined pinch and exergy method. Results show overall exergetic efficiency increases about 7.5% through MOEA method. Finally, the exergetic product cost of electricity reduces to 0.0183 $/MJ consecutively


Main Subjects

1. Introduction

Exergy, thermoeconomic and pinch approaches can be used for evaluation of process plant. Exergy analysis usually estimates the thermodynamic performance. Entropy generation is calculated by the entropy balance, and exergy balance calculates irreversibility. Also, thermoeconomic analysis predicts the unit cost of products such as power, steam and calculates monetary loss (Sanjay, Singh, & Prasad, 2007; Zhang, Wang, Zheng, & Lou, 2006). Furthermore, it supplies a tool for the optimum design and operation optimization of thermal systems. Nowadays, such analysis is vital because accurate prediction of the production costs is necessary for companies to work profitably(Modesto & Nebra, 2006).

Furthermore, combined pinch and exergy analysis guide us to the better understanding of the system by graphical representation as Feng introduced in 1997(Feng & Zhu, 1997). The power of pinch analysis is that the system can be represented by simple diagrams.  So, targets are easily achieved prior to the design phase. The strength of exergy analysis is that it can demonstrate the major inefficiency of the plant. In this regard, promising modifications can be indicated easily. For better improvement of base design, by applying the advantages of these analyses, the whole system can be shown in one diagram which guides to find the efficient modifications quickly.

Heuristic rules are often applied in the design and improvement of energy systems to simplify the problem. Also, combined pinch and exergy analysis as a graphical tool can guide us to an indication of best modification and constraints for optimization(Kotas, 2013). Rules taken from the fields of artificial intelligence (Frank Cziesla, 2000) and computational intelligence fuzzy systems (F. Cziesla & Tsatsaronis, 2002) and evolutionary algorithms can help the designer to achieve an optimum cost-effective case.

Generally, the electricity cost is more sensitive to modifies in the configuration of the process structure than to modified values of the variables. The objective of this research is to optimization with an evolutionary algorithm ("Matlab 7.1 Tutorial "). 

Uhlenbruck and Lucas (Uhlenbruck & Lucas, 2004) applied an exergoeconomic approach (Bejan, Tsatsaronis, Moran, Moran, & Moran, 1996) with an evolution concept to speed up the search for an optimum economic design of simple power plant concepts. By the evolution approach, only process variables are modified. Emmerich et al integrate a knowledge-based system based on an evolutionary algorithm which produces feasible design configurations (M. Emmerich, Gr, tzner, Sch, & tz, 2001; Michael Emmerich, Grötzner, Groß, & Schütz, 2000). So, the requirement for superstructure model can be removed.

Sepelling et al proposed a thermoeconomc optimization for a combined cycle with a solar tower power plant. A dynamic model of a solar combined-cycle power plant has been proposed.   Two objectives namely minimal investment costs and minimal electricity costs will be considered (Spelling, Favrat, Martin, & Augsburger, 2012).

Bracco et al proposed a MILP optimization model for a combined heat and power plant based on economic and environmental objectives. The multi-objective optimization of operation and capital costs and carbon dioxide emissions have been considered (Bracco, Dentici, & Siri, 2013).

Multi-objective exergoeconomic optimization of a CGAM solar-hybrid cogeneration cycle using genetic algorithm has been performed by Soltani et al (Soltani et al., 2014). Solar power tower plant has been considered.

 Mahmoudi et al focused on thermoeconomic modeling, evaluation and optimization of a novel combined supercritical carbon dioxide recompression based on Brayton/Kalina cycle (Mahmoudi, Delkhah Akbari, & Rosen, 2016).

Energy recovery and overall energy efficiency improvement in the gas transmission networks has been studied by Safarian & Mousavai. In this regard, they extended a model with detailed characteristics of compressor and pressure reduction stations. Furthermore, they determined three different scenarios with gas turbine including an organic Rankine cycle, air bottoming cycle, and air bottoming cycle along with steam injection (Safarian & Mousavi, 2015).

Multi-objective thermoeconomic optimization for combined-cycle power plant using particle swarm optimization has been performed by Abdalisousan et al (Abdalisousan, Fani, Farhanieh, & Abbaspour, 2015). Thermoeconomic optimization and parametric investigation and of a combined cycle for recovering the waste heat from nuclear-closed Brayton cycle have been performed by Luo et al (Luo, Gao, Liu, & Xu, 2016).

The energetic, exergetic and exergoeconomic evaluation of using different inlet air cooling systems in warm dry and wet climate stations have been performed by Khoshgoftar Manesh et al (Khoshgoftarmanesh, Vazini Modabber, & Mazhari, 2016).

Khoshgoftar Manesh and Ameryan used Cuckoo Search algorithm to optimal synthesis of a hybrid solar CHP cycle (M. H. Khoshgoftar Manesh & Ameryan, 2016). 

Memon et al investigated on thermo-environmental and economic evaluation of a combined cycle power plants based on regression modeling and optimization. In this regard, a relationship between optimal efficiency and the total cost is defined. In addition, multiple polynomials based on regression models are developed (Memon, Memon, & Qureshi, 2017).

The scope of this work is to find optimal configuration and process condition of a 160 MW gas-fired combined cycle based on thermoeconomic and environmental objectives. In this regard, the optimum structure and main process variables have been achieved by MOEA.  Furthermore, extended combined pinch and exergy approach has been employed to show the performance of cycle before and after optimization. In addition, the extended combined pinch and exergy method help us to determine the limits of the optimization problem to reduce the superstructure of an optimization problem.

In this paper, three configurations for HRSG have been considered. In addition, multi-objective optimization of a 160 MW combined based on the minimization of the cost of electricity using GA as well as maximizing exergy efficiency and minimization of environmental pollutions has been performed.

2. Case Study: A 160 MW Combined Cycle

A 160-MW combined cycle power plant was studied here. As demonstrated in the schematic diagram (Fig.1), this power plant has one gas turbine, one compressor, one HRSG, one deaerator, one steam turbine and a cooling system. The exhaust flue gas at 5460C enters the HRSG. The rated output power of the steam turbine is   49.638 MW (at 100% Load) at the design condition i.e. the ambient temperature of 250C. The stream data of a 160 MW combined cycle power plant is shown in Table 1.

3. Thermodynamic Modeling

Thermodynamic modeling has been performed as indicated in (Bejan et al., 1996; Dincer, Rosen, & Ahmadi, 2017; Koch, Cziesla, & Tsatsaronis, 2007; Tsatsaronis, 1993). The equations of the thermodynamic model can be found in Appendix A.  Furthermore, exergy equations have been illustrated in Appendix B. The assumptions of thermodynamic modeling are as follows (Dincer et al., 2017):

  • The plant operates at steady state.
  • Ideal-gas mixture assumption is applied to the combustion products and air.

So, we will have slight differences in simulation results rather than plant data.


4. Combined Pinch and Exergy Evaluation




Figure 1. PFD of a 160 MW combined cycle power plant



Figure 2. Exergy transformation from CC to ECC and EGCC (Feng & Zhu, 1997)

Table 1. Stream Data at 100% Load for 160MW Combined Cycle 



T  (۫°C)


m (kg/s)


Air in





Air out










GT inlet





GT Out





Stack gas





ST in





ST out





Condenser out





Pump in





CW in





CW out





Condenser  in

















Pinch technology has become a general approach for targeting and design of process plants and power stations (Feng & Zhu, 1997). The composite curves (CC) and the grand composite curve (GCC) are two main graphical curves, which are indicated by temperature versus enthalpy axes (Feng & Zhu, 1997; Mohammad Hasan Khoshgoftar Manesh & Amidpour, 2008; M. H. Khoshgoftar Manesh & Amidpour, 2009; Manesh & Rosen, 2018). The CC and GCC supply the targets set. The CC and GCC have been developed for use in the heat and power systems. Therefore, based on the Carnot factor (η) versus enthalpy, the exergy composite curves (ECC) and the exergy grand composite curve (EGCC) were proposed (Feng & Zhu, 1997; M. H. Khoshgoftar Manesh & Amidpour, 2009; Manesh & Rosen, 2018).

The CCs for the thermal system can be converted into the ECCs and the GCC. The shaded areas determine the exergy destruction related to the heat transfer. By a combination of pinch analysis and exergy analysis, it is possible to estimate the power demand or production for both power systems and refrigeration (Feng & Zhu, 1997). The combined pinch and exergy analysis was extended for targeting of shaft work. Particularly, only processes associated with heat transfer can be demonstrated on the ηc-H diagram but not the processes related to composition and pressure modifications. Therefore, the curve is constructed associated with temperatures:



The turbines and heat transfer systems are major parts of a thermal power plant. The chemical energy in the natural gas supplies the total exergy for the plant, as an exergy source. Part of the exergy from the fuel is lost in the heat transfer system, including the combustion chamber, heat recovery steam generator, and the condenser. The rest of the exergy enters into the turbine and compressor as the exergy input for power production. Some of the exergy destruction is associated with the driver such as turbines, compressors, and pumps which are related to by their machine efficiency. In addition, the exhausted gas lost the evident amount of input exergy. The remaining exergy gives the shaft work which is received by electrical generators, which becomes the final exergy sink (Feng & Zhu, 1997; Mohammad Hasan Khoshgoftar Manesh & Amidpour, 2008).

Based on pinch and exergy analysis, the Energy Level Curves (ELC) draws on the earlier strategies of the thermodynamic method to process integration, namely the concept of CCs and the combined pinch and exergy approach (Manesh & Rosen, 2018). Figure 2 shows the exergy transformation from CC to ECC and EGCC.

In this regard, the graphical representation as energy level (Ώ) defined as (Feng & Zhu, 1997; Mohammad Hasan Khoshgoftar Manesh & Amidpour, 2008; Manesh & Rosen, 2018):



Thus, for work



and for heat




and for a steady-state-flow system



5. Economic Analysis

The cost of the components has been taken into account in the economic model, including amortization, maintenance, and the cost of fuel demand. In order to introduce a function of cost which related to optimization parameters, the cost of components has to be defined as functions of thermodynamic parameters (Sanjay et al., 2007). These relationships can be constructed by statistical correlations between costs and the main thermodynamic. In this paper, cost equations have been considered based on (Dincer et al., 2017).

In this paper, to predict each equipment capital cost the Moran method was applied to calculate the total annualized cost (Bejan et al., 1996; Sanjay et al., 2007). The cost of amortization cost for equipment can be estimated as follows:





The annualized cost is calculated from the present worth of each component by applying the capital recovery factor CRF (i,n) (Sanjay et al., 2007). The capital cost for the kth component of the plant for 8000 operating hours per year is calculated as:



By assumption of 30 years of plant life, the maintenance and operating cost factor   are for each plant component (Sanjay et al., 2007).

6. Exergoeconomic Modeling and Analysis

The results of exergy analysis are used for exergoeconomic modeling and analysis (Bejan et al., 1996; Dincer et al., 2017). For the whole system an exergoeconomic balance as (Bejan et al., 1996; Dincer et al., 2017):




So, for a component receiving a heat transfer and generating power, we would write (Bejan et al., 1996; Dincer et al., 2017):




To solve for the unknown variables, it is necessary to develop a system of equations by applying Eq. (10) to each component, and in some cases we need to apply some additional equations, to fit the number of unknown variables with the number of equations (Bejan et al., 1996; Dincer et al., 2017). The general equations applied to each component are the following, according to Fig.1.  The cost balance equations for these components are as follows:

Air compressor






Gas turbine







Figure 3. Cost structure of a 160 MW gas combined cycle power plant



Figure 4. Basic concept of an evolutionary algorithm (Coello et al., 2002)


Steam turbine





Heat recovery steam generator











Feed water pump






Condenser pump






The cost balance equation for each component has been done. The exergoeconomic model of the case study is shown in Fig.3.

7. Optimization Strategy

7.1. Mathematical Programming

The equations have been integrated into an optimization framework developed to analysis the exergoeconomic optimization of the problem. The product cost of power generation is defined as an objective function. The optimization framework can modify process conditions to minimize exergetic production cost. In addition, three configurations for HRSG have been proposed. In addition, LINGO software as mathematical optimization program has been linked to Matlab environment by Excel. The iteration continues until solution convergence has been achieved.



Figure 5. Comparison between the single and the multi-objective approaches to the energetic and economic design optimization of energy systems (Coello et al., 2002)


Figure 6. Integration of different techniques and optimization procedure



7.2. Evolutionary Algorithms

7.2.1. Introduction

An iterative, stochastic search strategy as an evolutionary algorithm was applied to find the optimum structure and condition (Fig. 4) (Coello, Van Veldhuizen, & Lamont, 2002). For the plant, a black box model and one objective were applied to evaluate individual fitness. Pairs of individuals are selected to create new individuals based on their performance to optimize the objective function. 

Each individual is evaluated to calculate its fitness. The thermodynamic simulation, estimation of equipment costs, and an economic analysis have been evaluated (Deb, 2001).

If an individual violates any, additional penalty terms are added to the fitness value. After a fitness value has been assigned to each individual in the initial population, some of the individuals are selected for the mating pool. Selected individuals have more chances to produce in the next generation. In the next step, recombination and mutation operators are performed by the individuals in the mating pool, generating the off springs. These operators integrate randomly and change slightly the decision parameters of different individuals in the mating pool so that an offspring might get a better fitness than its parents. As shown in Fig.4, the iteration loop is repeated until the maximum number of generations is achieved. Optimization of evolutionary procedures includes evolution strategies and genetic algorithms. A detailed description of evolutionary computation and introduction is presented in (Abdalisousan et al., 2015; Coello et al., 2002; Deb, 2001; Michael Emmerich et al., 2000; M. H. Khoshgoftar Manesh & Amidpour, 2009; Lazzaretto & Toffolo, 2004).

Selection of appropriate variable settings for the evolutionary algorithm is a time-consuming task. An individual with a new process structure, even the optimal one, might need a few generations to improve the values of its process variables and to be competitive with the best individuals in the population. Regardless, the best existing individual produces some good performance off springs. The individual with the new process structure is often inferior to these solutions and is quickly removed from the population. Niching approaches within a population are capable to find and maintaining multiple solutions.

In each generation, 95% of the individuals are replaced by offspring. Since the initial population is generated randomly, all optimization runs were repeated several times with a different initial population. The optimal solutions for all cases differed only slightly. In the search for an optimal solution, an evolutionary algorithm only makes use of the value of the objective function. No gradient information is required (Coello et al., 2002; Deb, 2001).

7.2.2. Multi-Objective Optimization

The MOEA method applied as demonstrated in Fig. 5. In the first step, with each individual having the same probability of being parent selection has been done. In the next step, parents enter the reproduction step, producing μ offspring with a crossover approach. So, the values of the decision variable of the offspring convert a range defined by the decision variable values of the parents. Some decision variable values of the offspring are also randomly mutated with a probability Pmutt. Then, checking of the whole population of μ+μ individuals is performed for possible clones. Then, to encourage the finding of the search space, the clones are crashed and replaced with new randomly produced individuals. In the next step, the values of an objective function of the μ offspring are evaluated. According to the scheme of Goldberg, a Pareto ranking of the  individuals is done in multi-objective optimization (Coello et al., 2002; Deb, 2001; Lazzaretto & Toffolo, 2004).

Lastly, the number of generations go away is evaluated to create maximum number of generations; otherwise, the population for the next generation is begun. The thermo-economic objective function to be minimized expresses the total cost rate as the sum of fuel and investment (equipment/maintenance) cost rates(Deb, 2001; Dincer et al., 2017; M. H. Khoshgoftar Manesh & Amidpour, 2009; Lazzaretto & Toffolo, 2004):

Ctotal = C fuel (ε ( X )) + Cinv ( X )


Where ε is the exergetic efficiency, and X is a vector containing the design variables. Ctotal includes information in the fuel cost rate (C fuel = c fuel .E fuel) through E fuel (fuel exergy flow rate).This dependence does not change the nature of the objective function, which still accounts for the total money to be spent. Thus, the two terms of the function, that in principle are associated with two different objectives, merge into a single objective Ctotal. If Ctotal is minimized for a particular unit cost of fuel, only the extremum of the Pareto front corresponding to the economic minimum is found. On the other hand, in a single-objective approach considering both the thermodynamic and economic objectives (ε and Ctotal ), the overall single-objective function (to be maximized) would be constructed by combining the two functions through appropriate weights, as follows (Lazzaretto & Toffolo, 2004):

F ( X ) = w1 ⋅ ε ( X ) − w2 ⋅ Ctotal ( X , ε ( X )) = w1 ⋅ ε ( X ) − w2 ⋅ Ctotal ( X )


With w1 ≥ 0, w2 ≤ 1 and w1 + w2 = 1.

By modification of w1 and w2, the optimal conditions regarding the Pareto front of the multi-objective method are achieved as shown in Fig 5.

For multi-objective optimization, the Pareto approach is the type of multi-objective evolutionary algorithms (MOEAs). This method has been extended over the past decade (Lazzaretto & Toffolo, 2004); difficult tests on complex problems and on real engineering problems have determined that MOEAs can remove difficulties of classical approaches.

Due to MOEAs apply a population of solutions during explore; a single run will obtain multiple Pareto optimal solutions.

7.3. Environmental Objective

The minimization of environmental pollution regarding the operating of the combined cycle was considered.  So, the other objective function was introduced based on environmental emissions as environmental impact [23, 24]

Gülder proposed that the adiabatic flame temperature in the primary zone of the combustion chamber is derived as follows (Dincer et al., 2017; Gülder, 1986):

Tpz = A σ α exp( β (σ + λ ) 2 ) π x θ y ψ z


where π is a dimensionless pressure p/pref. In this regard, p is related to the combustion pressure p3 and pref is equal to 101.325 kPa. In addition, θ is T/Tref as a dimensionless temperature and T is related to the inlet temperature T3 and Tref is equal to 300 K. Furthermore, ψ is the H/C atomic ratio ψ = 4, the fuel is assumed to be pure methane; Also, σ = φ for φ ≤ 1 which φ related to the fuel to air equivalence ratio and σ = φ − 0.7 for φ > 1 (It is assumed that φ = 0.64); x, y, and z are quadratic functions of σ related to the following equations:

0.3 ≤ φ ≤ 1.0    and   0.92 ≤ θ < 2.0


0.3 ≤ φ ≤ 1.0    and   2.0 ≤ θ ≤ 3.2


1.0 < φ ≤ 1.6    and   0.92 ≤ θ < 2.0


1.0 < φ ≤ 1.6    and   2.0 ≤ θ ≤ 3.2 


As Rizk and Mongia were proposed, the pollutant emissions in grams per kilogram of fuel were calculated as (Rizk & Mongia, 1993):





The minimization of the total product cost of a plant as the main objective function has been considered.

In addition, another objective function is an environmental impact to be minimized. The environmental objective has been formulated in cost terms by multiplying the respective flow rates of pollutions damage cost (Lazzaretto & Toffolo, 2004)(cCO and cNOx are levelized back to 1994 to keep all cost values homogeneous and they are equal to 0.02086 $/kgCO and 6.853 $/kg NOx, respectively and then combined with economic objective. Therefore, in this research, the pollution damage cost has been considered in conjunction with economic objective and the optimization problem, named as environmental optimization, is a single-objective as follows:

CP = CF + ∑ Z k + Cenvtot,k


7.4. Decision Variables

In the optimization problem, the decision variables may be varied. All other variables are dependent variables. Their values are calculated from independent variables using a thermodynamic model. In this paper, decision variables have been selected as follows:

  • Gas turbine inlet;(T)
  • High-pressure steam  (m ; T)
  • Low-pressure steam   (m ; T)
  • Hot reheat steam   (m ; T)
  • HRSG exhaust  (T)
  • ηAC   (efficiency)
  • ηGT   (efficiency)
  • Compressor pressure ratio  (p2/p1)
  • HRSG configuration topology  (one,  two or triple pressure)

The range of the decision variable is as follow:

6 ≤ p2 / p1 ≤ 16


0.6 ≤ η GT≤ 0.92


0.6 ≤ η st ≤ 0.92


700 ≤ T3 ≤ 1000 K


1200 ≤ T4 ≤ 1550 K


The heat exchange between hot and cold streams in the HRSG must satisfy the following feasibility constraint:



8. Computer Program

Integration of different techniques and optimization procedures has been shown in Fig.6. As shown, in the first step thermodynamic simulation that is necessary for another analysis has been performed in Matlab code. In the next step, exergy and thermoeconomic analysis have been applied for the demonstration of irreversibility of equipment and estimation of electricity exergetic cost that is produced in steam and gas turbines. Objective function has been achieved by calculation of exergetic cost of electricityon and environmental pollution. After iteration best solution has been found. Iteration started from base case and it has been finished in optimum case. Three structures for HRSG have been proposed. Three configurations for HRSG and different process condition have been examined by evolutionary algorithms and best configuration and process condition have been selected as an optimal solution. Moreover, combined pinch-exergy representation demonstrates thermal performance of the system in base and optimum case. Also, it can be used for analysis and evaluation of the system.



Table 2. Main Data at Full Load for 160-MW combined cycle (Plant Data and Simulation Result) 

Parameter (Unit)

Plant data

Simulation Result

Ambient Temperature (oC)



Atmospheric pressure (bar)



Relative humidity (%)



Fuel :Natural gas, LHV (kJ/kg)



Combustion efficiency (%)



Gas turbine efficiency (%)



HRSG efficiency (%)



Net heat rate (kJ/kWh)



Gross heat rate



Gross electric efficiency (%)



Steam turbine shaft work (MW)



Gas turbine shaft work (MW)



Net total electric power (MW)



Net plant efficiency (%)



















Air Compressor














Gas Turbine







Steam Turbine





















CW Pump







FW Pump







 Table 3. Exergy and exergoeconomic parametes in each component in Base Case



By using stream data, the extended combined pinch and exergy representations have been constructed.

A computer program for thermodynamic, economic, exergy, thermoeconomic and combined pinch and exergy analyses of the 162-MW combined power plant has been developed in Matlab environment. The data inputs of this program as:

(a) Standard pressure (P0) and Standard temperature (T0);

(b) Compositions of fuel and fuel costs

(c) The composition of air and air relative humidity;

(d) Load conditions

(d) Gross shaft power of steam turbine and gas turbine;

(e) Required power for compressor and pumps;

(f) For streams at the inlet and outlet of each component pressure (MPa), mass flow rate (kg/s) and temperature (°C);

(g) Economic inputs such as capital cost, interest rate, and salvage value factor;

By applying these input data; composition of combustion products and adiabatic flame temperature can be calculated. Also, enthalpy and entropy for each fluid streams at various conditions can be computed. Also, computer code has been written for thermodynamic simulation and evaluation of plant. The enthalpy and entropy of non-interacting gas species were calculated by using appropriated polynomials fitted to the JANAF thermochemical tables ("JANAF thermochemical tables, 1975 supplement," 1975). Furthermore, by the International Association for the Properties of Water and Steam, enthalpy and entropy of water and steam were calculated (Wagner & Kretzschmar, 2007).

Using the values of these properties, the thermodynamic simulation and exergy analysis of plant‘s components have been performed. In addition, exergy balances for the plant boundary and components were created. Then, the unit cost of products was calculated by solving the cost balance equations simultaneously. Furthermore, this program can generate improved combined pinch and exergy representation. To perform the optimization of the 160-MW case study, the evolutionary algorithms of Multi-Objective Evolutionary through NSGA-II are interfaced with a MATLAB m-file model in which the thermodynamic, economic and environmental equations of the case problem are implemented and solved and then returning the values of the objectives for a given set of decision variables.  In addition, the thermoeconomic and environment MATLAB m-file model have been exported to LINGO to find optimal solution by MINLP.

9. Result

Exergy and thermoeconomic analysis have been performed for a 160-MW combined cycle power plant. Matlab code has been developed for the simulation of this power plant. Main data of 160-MW combined cycle in plant data and simulation case are shown in Table 2. The slight differences existing between the two columns show the accuracy of simulation results and due to some assumptions that have been used for thermodynamic modeling and simulation. Exergy and exergoeconomic parameters such as exergy destruction, the cost rate of fuel, cost rate of the product, exergy of fuel, exergy of product and efficiency of each component are shown in Table 3.



Figure 7. Optimum solution through multi-objective function


Figure 8. Sensitivity analysis of multi-objective optimization by modification of specific fuel cost corresponding to Pareto optimal sets


Figure 9. PFD of power plant-after optimisation



Multi-Objective Thermoeconmic optimization through Evolutionary Algorithm (NSGA-II) has been applied for 160- MW combined cycle power plant to minimization of generated electricity cost. In this research, three methods have been used for finding optimum configuration and process condition of combined cycle plants. The combined pinch exergy approach has been employed to demonstrate the limitation of problem. It helps us to reduce the superstructure of problem. Also, different analysis has been applied for better showing thermal behavior and the evaluation of the system in the base and optimum case by best optimization method.

In this case, the setting of the evolutionary algorithm is as follows:

Population size: 500

No. of Generations: 900

Pc (Probability of crossover): 0.7

No. of crossover points: 2

Pm (Probability of mutation): 0.01

Selection process: Tournament

Tournament Size: 2

The Pareto Front for this multi-objective optimization problem, which shows the best trade-off values of two objective functions, has been presented in Fig.7.

As shown in this figure, while the total exergetic efficiency of the plant increased to about 55%, the total cost rate of the products increases very slightly. Increasing of total exergetic efficiency from 55% to 57.3% is corresponding to a moderate increase of cost rate of the product. Eventually, from 57.3% to a higher value of ε tot, the total cost rate increases sharply. It should be noted that the selection of the optimum solution is depending on preferences and criterions of each decision maker. Therefore, another person may select a different point as the optimum solution which better suits his desires. Therefore, in the optimum solution exergetic efficiency is about 56.8% and cost of electricity is 0.01384 $/MJ.

Also, mathematical optimization has been applied by LINGO software. Table 4 shows the comparison of mathematical programming and one objective GA and MOEA methods. As shown in this table, MOEA gives us a better solution rather than other methods.

Sensitivity analysis of multi-objective optimization by modification of specific fuel cost corresponding to Pareto optimal sets has been shown in Fig.8.   As shown in these figures, by increasing specific fuel cost or interest rate, specific production cost has been increased.

The value of characteristic process variables after the calculation has been determined in Table 5.

The best structure of HRSG after optimization has been shown in Fig.9.  In this study, the configuration of HRSG has been modified through the Evolutionary algorithm. The profiles of HRSG temperature in the base case and the optimum case has been illustrated in Fig.10 and Fig.11.

In addition, for better evaluation of component’s performance in the base case and optimum case, extended combined pinch and exergy approach has been applied to demonstrate the thermal behavior of the plant.

The energy level curve helps us to a better understanding of the thermal behavior of the power plant's components in the base and optimal case as shown in Fig.12 and Fig.13 consequently.



Table 4. Results Comparison of Different Methods (MINLP, GA-One Objective and MOEA)




GA One Objective


Net total electric power (MW)




Net plant efficiency (%)




CP steam turbine($/MJ)





Table 5. Values of characteristic process variables after 60 generation





Gas turbine inlet

m ; T

372.51 kg/s ; 597.49 °C

High pressure steam

m ; T

38.01 kg/s ; 319.31 °C

Low pressure steam

m ; T

2.102 kg/s ; 107.34 °C

 Hot reheat steam

m ; T

47,94 kg/s ; 539.18 °C

HRSG exhaust


160.12 °C



84.17 %



92.16 %

Compressor pressure ratio






Figure 10.HRSG temperature profile-base case


Figure 11.HRSG temperature profile-optimal case



Energy level curve that is generated by Matlab code shows the effects of this improvement in different component and thermal interaction in whole plant graphically.

In addition, comparison of exergy destruction in the base and optimum case after thermo economic optimization with evolutionary algorithm are shown in Fig.14. Also, figure 15 compares the cost of exergy destruction in the base and optimum case.




Figure 12.Energy level representation- base case


Figure 13. Energy level representation- optimal case


Figure 14.Comparison of exergy destruction in base and optimal case


Figure 15.Comparison of exergy cost destruction in base and optimal case



10. Conclusion

In this paper, Matlab m-file program has been developed for the thermodynamic simulation of 160-MW combined cycle power plant. An exergoecomic has been applied to this plant to predict the unit costs of power produced from gas and steam turbines.

Also, this program that has been developed which shows that the exergy, thermo economic, combined pinch-exergy analysis and thermo economic optimization of configurations in HRSG and process variables via NSGA-II. This algorithm presented here can be applied to combined cycle power plant systematically.

In addition, extended combined pinch and exergy approach have the strength for promising modifications and feasible region for optimization.  It can help us to identify major improvement. Also, it can be applied for evaluation and comparison of the base and optimum case. As shown in the results, overall exergetic efficiency increases about 7.5% by MOEA approach. In addition, the exergetic product cost of electricity reduces to 0.0183 $/MJ consecutively. As sensitivity analysis has been shown, by increasing specific fuel cost or interest rate, specific production cost has been increased.

The proposed evolutionary algorithm has been shown to be a powerful and effective tool in finding the set of the optimal solutions rather than multi-objective mathematical programming and GA one objective optimization for finding optimum design variables in this plant. In this study, the combined pinch exergy analysis has been performed to determine the optimization potentials and limitations.  Therefore, this evolutionary-based procedure will be very useful for the optimization of complex thermal systems.


A     constant coefficient

e      exergy rate per mass

E     time rate of exergy

T     temperature

P     pressure

m    mass flow rate

W    shaft work

C     cost rate

Z     cost rate of capital investment and O&M

x     quadratic functions

y     quadratic functions

z     quadratic functions

Greek letters

    carnot factor

W        energy level

π     dimensionless pressure

θ      dimensionless temperature

ψ     H/C atomic ratio

φ     the fuel to air equivalence ratio

ε      exergetic efficiency

σ    constant related to fuel to air equivalence ratio

   time constant

  constant coefficient

  constant coefficient

   constant coefficient


CI      capital investment

OM    operating and maintenance cost


p      product

f       fuel

D     destruction

L      loss

pz  the adiabatic flame temperature in the primary zone of the combustion chamber


Appendix A (Thermodynamic model)


The enthalpy and entropy of gas species were calculated by using polynomials fitted to the thermo physical data in the JANAF Tables [29]. In addition, the values of as enthalpy and entropy for water and steam were evaluated by using by the International Association for the Properties of Water and Steam (IAPWS-IF97) [30].

Governing Equations

A set of governing equations can be developed as follows [10]:

Combustion Chamber

Denoting the fuel-air ratio on molar basis on a molar basis as λ, the molar flow rates of the fuel, air, and combustion products are related by [10]:



where the subscripts F, P and a denote, respectively, fuel, combustion products and air. For complete combustion of methane, the chemical equation takes the form [10]:




The molar analysis of the combustion products is fixed once the fuel-air ratio λ has been determined. The fuel air ratio can be obtained from an energy rate balance as follows [10]:





As the heat loss is assumed to be 2% of the lower heating value, we have [10]:





Turbine and Compressor

For this control volume the energy rate balance takes the form [10]:





The term (h1 − h2) of Eq. (45) is evaluated using the compressor isentropic efficiency [10]:



where h2sdenotes the specific enthalpy for an isentropic compression from inlet state 1 to the specified exit pressure p2. The state 2s fixed using s2ss1 = 0.

The value of h2 determined from Eq. (45) is used to calculate the value of T2 by solving an equation derived from Eq. (46). The term (h7h8) of Eq. (47) is evaluated using the turbine isentropic efficiency [10]:




Heat Recovery Steam Generator

For this control volume the energy rate balance takes the form [10]:




Appendix B (Exergy Analysis)


The exergy destruction in the combustor is calculated as [10]:



Ein is sum of fuel exergy and air exergy that input to the combustion chamber. Eout is exergy of combustion that produces in combustor [6].


In this cycle we have one gas and steam turbine. Exergy destruction for turbines is defined as:



 is the shaft work. The exergetic efficiency of the turbines introduced as the ratio of the minimum work input to the actual work input [6, 7] as follows:




Heat Exchanger

Heat recovery steam generator and condensers are essentially heat exchangers designed to perform different tasks. The exergy model is defined as follows:



The exergetic efficiency of a heat exchanger is defined as follows [6]:




Compressor or Pump

The exergy destruction in compressor or pump can be defined as:



The exergetic efficiency of the compressor or pump can be defined as:



3.5. Cycle

The overall exergetic efficiency of the cycle can be defined as:



 Abdalisousan, A., Fani, M., Farhanieh, B., & Abbaspour, M. (2015). Multi-objective thermoeconomic optimisation for combined-cycle power plant using particle swarm optimisation and compared with two approaches: an application. International Journal of Exergy, 16(4), 430-463. doi:10.1504/IJEX.2015.069112
Bejan, A., Tsatsaronis, G., Moran, M., Moran, M. J., & Moran, P. G. M. C. M. (1996). Thermal Design and Optimization: Wiley.
Bracco, S., Dentici, G., & Siri, S. (2013). Economic and environmental optimization model for the design and the operation of a combined heat and power distributed generation system in an urban area. Energy, 55, 1014-1024. doi:
Coello, C. A. C., Van Veldhuizen, D. A., & Lamont, G. B. (2002). Evolutionary Algorithms for Solving Multi-Objective Problems: Kluwer Academic.
Cziesla, F. (2000). Produktkostenminimierung beim Entwurf komplexer Energieumwandlungsanlagen mit Hilfe von wissensbasierten Methoden (PhD), Berlin, Techn. Univ, German Düsseldorf
Cziesla, F., & Tsatsaronis, G. (2002). Iterative exergoeconomic evaluation and improvement of thermal power plants using fuzzy inference systems. Energy Conversion and Management, 43(9), 1537-1548. doi:
Deb, K. (2001). Multi-Objective Optimization Using Evolutionary Algorithms: Wiley.
Dincer, I., Rosen, M. A., & Ahmadi, P. (2017). Optimization of Energy Systems: Wiley.
Emmerich, M., Gr, M., tzner, Sch, M., & tz. (2001). Design of Graph-Based Evolutionary Algorithms: A Case Study for Chemical Process Networks. Evolutionary Computation, 9(3), 329-354. doi:10.1162/106365601750406028
Emmerich, M., Grötzner, M., Groß, B., & Schütz, M. (2000). Mixed-Integer Evolution Strategy for Chemical Plant Optimization with Simulators, London.
Feng, X., & Zhu, X. X. (1997). Combining pinch and exergy analysis for process modifications. Applied Thermal Engineering, 17(3), 249-261. doi:
Gülder, O. L. (1986). Flame Temperature Estimation of Conventional and Future Jet Fuels. Journal of Engineering for Gas Turbines and Power, 108(2), 376-380. doi:10.1115/1.3239914
JANAF thermochemical tables, 1975 supplement. (1975). Journal of Physical and Chemical Reference Data, 4(1), 1-176. doi:10.1063/1.555517
Khoshgoftar Manesh, M. H., & Ameryan, M. (2016). Optimal design of a solar-hybrid cogeneration cycle using Cuckoo Search algorithm. Applied Thermal Engineering, 102, 1300-1313. doi:
Khoshgoftar Manesh, M. H., & Amidpour, M. (2008). New Graphical Methodology for Energy Integration in Nuclear Steam Power Plant. (48140), 139-155. doi:10.1115/ICONE16-48432
Khoshgoftar Manesh, M. H., & Amidpour, M. (2009). Multi-objective thermoeconomic optimization of coupling MSF desalination with PWR nuclear power plant through evolutionary algorithms. Desalination, 249(3), 1332-1344. doi:
Khoshgoftarmanesh, M. H., Vazini Modabber, H., & Mazhari, V. (2016). Techno-Economic Assessment of Different Inlet Air Cooling Systems in Warm Dry & Wet Climate Stations. Gas Processing, 4(2), 19-31. doi:10.22108/gpj.2017.100042.1001
Koch, C., Cziesla, F., & Tsatsaronis, G. (2007). Optimization of combined cycle power plants using evolutionary algorithms. Chemical Engineering and Processing: Process Intensification, 46(11), 1151-1159. doi:
Kotas, T. J. (2013). The Exergy Method of Thermal Plant Analysis: Elsevier Science.
Lazzaretto, A., & Toffolo, A. (2004). Energy, economy and environment as objectives in multi-criterion optimization of thermal systems design. Energy, 29(8), 1139-1157. doi:
Luo, L., Gao, H., Liu, C., & Xu, X. (2016). Parametric Investigation and Thermoeconomic Optimization of a Combined Cycle for Recovering the Waste Heat from Nuclear Closed Brayton Cycle. Science and Technology of Nuclear Installations, 2016, 12. doi:10.1155/2016/6790576
Mahmoudi, S. M. s., Delkhah Akbari, A., & Rosen, M. (2016). Thermoeconomic Analysis and Optimization of a New Combined Supercritical Carbon Dioxide Recompression Brayton/Kalina Cycle (Vol. 8).
Manesh, M. H. K., & Rosen, M. A. (2018). Combined Cycle and Steam Gas-Fired Power Plant Analysis through Exergoeconomic and Extended Combined Pinch and Exergy Methods. Journal of Energy Engineering, 144(2), 04018010. doi:doi:10.1061/(ASCE)EY.1943-7897.0000506 . Matlab 7.1 Tutorial
Memon, A. G., Memon, R. A., & Qureshi, S. (2017). Thermo-environmental and economic analyses of combined cycle power plants with regression modelling and optimization. Applied Thermal Engineering, 125, 489-512. doi:
Modesto, M., & Nebra, S. A. (2006). Analysis of a repowering proposal to the power generation system of a steel mill plant through the exergetic cost method. Energy, 31(15), 3261-3277. doi:
Rizk, N. K., & Mongia, H. C. (1993). Semianalytical Correlations for NOx, CO, and UHC Emissions. Journal of Engineering for Gas Turbines and Power, 115(3), 612-619. doi:10.1115/1.2906750
Safarian, S., & Mousavi, M. (2015). Improvement of Overall Efficiency in the Gas Transmission Networks: Employing Energy Recovery Systems. Gas Processing, 3(2), 1-24. doi:10.22108/gpj.2015.20403
Sanjay, Y., Singh, O., & Prasad, B. N. (2007). Energy and exergy analysis of steam cooled reheat gas–steam combined cycle. Applied Thermal Engineering, 27(17), 2779-2790. doi:
Soltani, R., Mohammadzadeh Keleshtery, P., Vahdati, M., KhoshgoftarManesh, M. H., Rosen, M., & Amidpour, M. (2014). Multi-objective optimization of a solar-hybrid cogeneration cycle: Application to CGAM problem (Vol. 81).
Spelling, J., Favrat, D., Martin, A., & Augsburger, G. (2012). Thermoeconomic optimization of a combined-cycle solar tower power plant. Energy, 41(1), 113-120. doi:
Tsatsaronis, G. (1993). Thermoeconomic analysis and optimization of energy systems. Progress in Energy and Combustion Science, 19(3), 227-257. doi:
Uhlenbruck, S., & Lucas, K. (2004). Exergoeconomically—aided evolution strategy applied to a combined cycle power plant. International Journal of Thermal Sciences, 43(3), 289-296. doi:
Wagner, W., & Kretzschmar, H. J. (2007). International Steam Tables - Properties of Water and Steam based on the Industrial Formulation IAPWS-IF97: Tables, Algorithms, Diagrams, and CD-ROM Electronic Steam Tables - All of the equations of IAPWS-IF97 including a complete set of supplementary backward equations for fast calculations of heat cycles, boilers, and steam turbines: Springer Berlin Heidelberg.
Zhang, C., Wang, Y., Zheng, C., & Lou, X. (2006). Exergy cost analysis of a coal fired power plant based on structural theory of thermoeconomics. Energy Conversion and Management, 47(7), 817-843. doi: