Property values and thermal, and mechanical exergy flows and entropy production rates at various state points in the gas-turbine power plant at 100% load condition.
Exergy costing to estimate the unit cost of products from various power plants and refrigeration system is discussed based on modified-productive structure analysis (MOPSA) method. MOPSA method provides explicit equations from which quick estimation of the unit cost of products produced in various power plants is possible. The unit cost of electricity generated by the gas-turbine power plant is proportional to the fuel cost and inversely proportional to the exergetic efficiency of the plant and is affected by the ratio of the monetary flow rate of non-fuel items to the monetary flow rate of fuel. On the other hand, the unit cost of electricity from the organic Rankine cycle power plant with heat source as fuel is proportional to the unit cost of heat and the ratio of the monetary flow rate of non-fuel items to the generated electric power, independently. For refrigeration system, the unit cost of heat is proportional to the consumed electricity and inversely proportional to the coefficient of performance of the system, and is affected by the ratio of the monetary flow rate of non-fuel items to the monetary flow rate of consumed electricity.
- unit exergy cost
- power plant
- refrigeration system
Exergy analysis is an effective tool to accurately predict the thermodynamic performance of any energy system and the efficiency of the system components and to quantify the entropy generation of the components [1, 2, 3]. By this way, the location of irreversibilities in the system is determined. Furthermore, thermoeconomic analysis provides an opportunity to estimate the unit cost of products such as electricity and/or steam from thermal systems [4, 5] and quantifies monetary loss due to irreversibility for the components in the system . Also, thermoeconomic analysis provides a tool for optimum design and operation of complex thermal systems such as cogeneration power plant  and efficient integration of new and renewable energy systems . Recently, performance evaluation of various plants such as sugar plant , drying plant , and geothermal plant  has been done using exergy and thermoeconomic analyses. In this chapter, a procedure to obtain the unit cost of products from the power plants and refrigeration system is presented by using modified-productive structure analysis (MOPSA) method. The power plants considered in this chapter are gas-turbine power plant and organic Rankine cycle power plant. These systems generate electricity as a product by consuming the heat resultant from combustion of fuel and by obtaining heat from any hot stream as fuel, respectively. In addition, MOPSA method is applied to an air-cooled air conditioning system, which removes heat like a product while the consumed electricity is considered as fuel. Explicit equations to estimate the unit cost of electricity generated by the gas-turbine power plant and organic Rankine cycle plant, and the unit cost of heat for the refrigeration system are obtained and the results are presented.
2. A thermoeconomic method: modified productive structure analysis (MOPSA)
2.1. Exergy-balance and cost balance equations
A general exergy-balance equation that can be applied to any component of thermal systems may be formulated by utilizing the first and second law of thermodynamics . Including the exergy losses due to heat transfer through the non-adiabatic components, and with decomposing the material stream into thermal and mechanical exergy streams, the general exergy-balance equation may be written as 
The fourth term in Eq. (1) is called the neg-entropy which represents the negative value of the rate of lost work due to entropy generation, which can be obtained from the second law of thermodynamics. The term in Eq. (1) denotes the rate of exergy flow of fuel, and in the fourth term denotes heat transfer interaction between a component and the environment, which can be obtained from the first law of thermodynamics.
However, the quantity for each component, which is usually not measured, may be obtained from the corresponding exergy-balance equation with the known values of the entropy flow rate at inlet and outlet.
Exergy, which is the ability to produce work, can be defined as the differences between the states of a stream or matter at any given particular temperature and pressure and the state of the same stream at a reference state. The exergy stream per unit mass is calculated by the following equation:
where T is temperature, P is pressure, and the subscript ref denotes reference values. The exergy stream per unit mass can be divided into its thermal (T) and mechanical (P) components as follows :
Assigning a unit exergy cost to every exergy stream, the cost-balance equation corresponding to the exergy-balance equation for any component in a thermal system  may be written as
The term includes all financial charges associated with owning and operating the kth component in the thermal system. We call the thermoeconomic analysis based on Eqs. (1) and (7) as modified-productive structure analysis (MOPSA) method because the cost-balance equation in Eq. (7) yields the productive structure of the thermal systems, as suggested and developed by Lozano and Valero  and Torres et al. . MOPSA has been proved as very useful and powerful method in the exergy and thermoeconomic analysis of large and complex thermal systems such as a geothermal district heating system for buildings  and a high-temperature gas-cooled reactor coupled to a steam methane reforming plant . Furthermore, the MOPSA can provide the interaction between the components in the power plant through the entropy flows  and a reliable diagnosis tool to find faulty components in power plants .
2.2. Levelized cost of system components
All costs due to owning and operating a plant depend on the type of financing, the required capital, the expected life of components, and the operating hours of the system. The annualized (levelized) cost method of Moran  was used to estimate the capital cost of components in this study. The amortization cost for a particular plant component is given by
The present worth of the component is converted to annualized cost by using the capital recovery factor CRF(i, n):
The capital cost rate of the kth component of the thermal system can be obtained by dividing the levelized cost by annual operating hours δ.
The maintenance cost is taken into consideration through the factor . It is noted that the operating hours of thermal systems is largely dependent on the energy demand patterns of end users .
3. Gas-turbine power plant
A schematic of a 300 MW gas-turbine power plant considered in this chapter is shown in Figure 1. The system includes five components: air compressor (1), combustor (2), gas turbines (3), fuel preheater (5), and fuel injector (6). A typical mass flow rate of fuel to the combustor at full load condition is 8.75 kg/s and the air–fuel mass ratio is about 50.0. Thermal and mechanical exergy flow rates and entropy flow rate at various state points shown in Figure 1 are presented in Table 1. These flow rates were calculated based on the values of measured properties such as pressure, temperature, and mass flow rate at various state points.
3.1. Exergy-balance equation for gas-turbine power plant
The following exergy-balance equations can be obtained by applying the general exergy-balance equation given in Eq. (1) to each component in the gas-turbine power plant.
The net flow rates of the various exergies crossing the boundary of each component in the gas-turbine power plant at 100% load condition are shown in Table 2. Positives values of exergies indicate the exergy flow rate of “products,” while negative values represent the exergy flow rate of “resources” or “fuel.” The irreversibility rate due to entropy production in each component acts as a product in the exergy-balance equation. The sum of exergy flow rates of products and resources equals to zero for each component and the overall system; this zero sum indicates that perfect exergy balances are satisfied.
|Component||Net exergy flow rates (MW)||Irreversibility rate (MW)|
3.2. Cost-balance equation for gas-turbine power system
When the cost-balance equation is applied to a component, a new unit cost must be assigned to the component’s principle product, whose unit cost is expressed as Gothic letter. After a unit cost is assigned to the principal product of each component, the cost-balance equations corresponding to the exergy-balance equations are as follows:
Applying the general cost-balance equation to the system components, five cost-balance equations are derived. However, these equations present eight unknown unit exergy costs, which are CT, CS, CW, C1P, C2T, CP, C5T, and C6P. To calculate the value of these unknown unit exergy costs, three more cost-balance equations are required. These additional equations can be obtained from the thermal and mechanical junctions and boundary of the plant.
Thermal exergy junction
Mechanical exergy junction
In Table 3, initial investments, the annuities including the maintenance cost, and the corresponding monetary flow rates for each component are given. The cost flow rates corresponding to a component’s exergy flow rates at 100% load condition are given in Table 4. The same sign convention for the cost flow rates related to products and resources was used as the case of exergy balances shown in Table 2. The lost cost due to the entropy production in a component is consumed cost. The fact that the sum of the cost flow rates of each component in the plant becomes zero, as verified in Table 4, shows that all the cost balances for the components are satisfied.
|Component||Initial investment cost|
|Monetary flow rate|
From the above equation, the unit cost of electricity for the gas-turbine power system is given as 
The production cost depends on fuel cost and the exergetic efficiency of the system, and is affected by the ratio of the monetary flow rate of non-fuel items to the monetary flow rate of fuel. With the unit cost of fuel, Co = 5.0 $/GJ, an exergetic efficiency of the gas-turbine power plant, 0.341, and a value of the ratio of the monetary flow rate of non-fuel items to the monetary flow rate of fuel, 0.22, the unit cost of electricity estimated from Eq. (26) is approximately 17.97 $/GJ. However, one should solve Eqs. (17)–(24) simultaneously to obtain the unit cost of electricity and the lost cost flow rate occurred in each component.
4. Organic Rankine cycle power plant using heat as fuel
A schematic of the 20-kW ocean thermal energy conversion (OTEC) plant  operated by organic Rankine cycle, which is considered to apply MOPSA method, is illustrated in Figure 2. Five main components exist in the system: the evaporator (1), turbine (2), condenser (3), receiver tank (4) and pump (5). The refrigerant stream is heated by a heat source in the evaporator, and then the refrigerant stream is divided into two streams. A portion of this stream is passed through the throttling valve and reaches the receiver tank, while the remaining part of the refrigerant stream leaving from evaporator is sent to turbine. A portion of the stream flowing to turbine is throttled and bypassed to turbine outlet. The “pipes” are introduced into the analysis as a component to consider the heat and pressure losses in the pipes and the exergy removal during the throttling processes. Refrigerant of R32 is used as a working fluid in the organic Rankine cycle. At the full load condition, the mass flow rate of the refrigerant is 3.62 kg/s. The warm sea water having mass flow rate of 86.99 kg/s is used as a heat source for the plant, while the cold sea water having mass flow rate of 44.85 kg/s is used as a heat sink for the plant. The reference temperature and pressure for the refrigerant R32 are −40°C and 177.60 kPa, respectively. For water, the reference point was taken as 0.01°C, the triple point of water.
4.1. Exergy-balance equations for the organic Rankine cycle power plant
The α term given in Eq. (32) is the ratio of the bypass streams from state 103 to 108. The value of the α term can be calculated by applying the mass and energy conservation equations to the receiver tank. The stream bypassed from state 103 to 105 may be neglected. An example of exergy calculation for the organic Rankine cycle plant using a stream of warm water at 28°C as a heat source to the evaporator  is shown in Table 5. As mentioned in the previous section, a positive value of exergy flow rate represents “product,” while a negative value of exergy flow rate indicates “fuel.” The last two columns clearly indicate that the electricity comes from expenditure of heat input.
|Component||Refrigerant||Water stream||Irreversibility rate||Heat transfer rate||Work input/output rate|
4.2. Cost-balance equations for the organic Rankine cycle power plant
By assigning a unit cost to every thermal exergy of the refrigerant stream (C1T, C2T, C3T, and CT), mechanical exergy for the refrigerant stream (CP), cold water (C3), neg-entropy (Cs), and electricity (CW), the cost-balance equations corresponding to the exergy-balance equations which are Eqs. (27)–(33) are given as follows. When the cost-balance equation is applied to a specific component, one may assign a unit cost to its main product, which is represented by a Gothic letter.
Seven cost-balance equations for the five components of the plant, pipes, and the boundary were derived with eight unknown unit exergy costs of C1T, C2T, C3T, CT, CP, C3, CS, and CW. We can obtain an additional cost-balance equation for the junction of thermal exergy of the refrigerant stream.
With Eq. (41), we have all the necessary cost-balance equations to calculate the unit cost of all exergies (C1T, C2T, C3T, CT, and C3, neg-entropy (Cs) and a product (electricity, CW) by input (given) of thermal energy (C2) to the evaporator. The overall cost-balance equation for the Rankine power plant can be obtained by summing Eqs. (34)–(41), which is given by
where is the net cost flow rate due to the heat transfer to/from the organic Rankine cycle plant. The term in Eq. (42) may represent the cost flow rate related to the construction of the plant . Rewriting Eq. (42), we have the unit cost of electricity from the Rankine cycle power plant.
where is the net electricity obtained from the organic Rankine cycle plant and abs denotes the absolute value of the quantity in parentheses.
Figure 3 shows that the unit cost of electricity from the organic Rankine cycle plant and the net cost flow rate due to the heat transfer rate to the plant vary depending on the unit cost of warm water in the evaporator, C2, appeared in Eq. (34). As the unit cost of warm water increases, the net cost flow rate due to heat transfer to the plant decreases while the unit cost of electricity increases. The cross point between the line for the unit cost of electricity and the line for the total cost flow rate due to heat transfer determines unit cost of electricity. The unit cost of electricity and the net cost flow rate due to heat transfer for a case whose detailed calculation results shown in Table 6 are $0.205 and −$0.941/kWh, respectively. The value of the unit cost C2 appeared in the cost balance equation, Eq. (34), is approximately $0.117/kWh for this particular case, which may be considered as a fictional one.
Detailed calculation results reveal that the unit cost of electricity from an organic Rankine cycle plant can be obtained from the following equation:
The calculated value of using Eq. (45) is approximately $0.055/kWh, which is quite different from the C2 = $0.177/kWh, a value determined from Figure 3. The value of which is the ratio of the absolute value of net cost flow rate of heat to the produced electricity was found to be a real unit cost of hot water stream . Equation (44) tells us that the unit cost of electricity from the organic Rankine cycle is determined by the sum of the unit cost of heat and the ratio of the monetary flow rate of non-fuel items to the produced electric power.
5. 200 kW air-cooled air conditioning unit
Even though the performance evaluation of a household refrigerator using thermoeconomics was performed , estimation of the unit cost of heat supplied to the room by air conditioning unit was never tried. In this section, the unit cost of heat for a 120-kW air-cooled air conditioning unit is obtained, which is helpful for the cost comparison between air conditioning unit operated by electricity and absorption refrigeration system running by heat .
5.1. Exergy-balance equations for the air conditioning units
The exergy-balance equations obtained using Eq. (1) for each component in an air-cooled air conditioning units shown in Figure 4 are as follows. The heat transfer interactions with environment for the compressor, TXV, and suction line are neglected.
The deposition of heat into the environment and the heat transferred to room are hardly considered to be dissipated to the environment. For such heat delivery system, it may be reasonable that the delivered heat rather than its exergy is contained in the exergy-balance equation. With help of Eqs. (51) and (52), the exergy-balance equation for the condenser and evaporator become
The terms, in Eq. (47’) and in Eq. (49’) represent the irreversibility corresponding to the terms and , respectively. The pair terms in the second bracket in Eqs. (47’) and (49’) are equal in magnitude but opposite in sign to vanish completely because the terms in the first bracket in those equations vanish. This assumption is legitimate since the entropy generation due to the heat transfer between flow streams  in the condenser and evaporator was calculated to be negligibly small.
The simulated data for the difference in the thermal and mechanical exergy flow rates at each component under normal operation for a 120-kW air-cooled air conditioning system  is displayed in Table 7. The cooling capacity of the system () is considered as the heat gained by the refrigerant in the evaporator. However, the heat output to the environment through the condenser was taken to satisfy the exergy-balance equation for the condenser as well as the overall system. The irreversibility rate due to the entropy generation at each component was calculated using the exergy-balance equations for each component given from Eq. (46) to (50). The values in the parentheses in the third and fifth columns represent the first and second quantity inside the second bracket in Eqs. (47’) and (49’), respectively. Note that minus and plus sign indicate the resource or fuel and product exergies, respectively, as usual. The sign of the irreversibility rate is minus at the evaporator, while it is plus at other units which play as boundary.
5.2. Cost-balance equations for the air-cooled air conditioning units
By assigning a unit cost to every thermal and mechanical exergy stream of the refrigerant (CT, CP), lost work (CS), heat (CH), and work (CW), the cost-balance equations corresponding to the exergy-balance equations, i.e., Eqs. (46), (47’), (48), (49’), and (50), are as follows. In this particular thermal system, a unit to a principal product for each component is not applied because the working fluid that flows through all the components makes a thermodynamic cycle.
We now have five cost-balance equations to calculate two unit costs of exergies (CT and CP), neg-entropy (CS), and a product, heat (CH) by input of electricity (CW). So, it is better to combine the cost-balance equation for the evaporator and suction line, which can be written as
Table 8 lists the initial investment, the annuities including the maintenance cost, and the corresponding monetary flow rates for each component of the air-cooled air conditioning system. Currently, the installation cost of an air-cooled air conditioning system with a 120-kW cooling capacity is approximately $17,000 in Korea. The levelized cost of the air conditioning units was calculated to be 0.3122$/h with an expected life of 20 years, an interest rate of 5% and salvage value of $850. The operating hours of the air conditioning system, which is crucial in determining the levelized cost, were taken as 4500 h. The maintenance cost was taken as 5% of the annual levelized cost of the system.
|Component||Initial investment ($)||Annualized cost ($/year)||Monetary flow rate ($/h)|
|Evaporator + Suction line||6000||472.1||0.1102|
The cost flow rates of various exergies and irreversibility rate at each component in the air conditioning system at the normal operation are shown in Table 9. The sign convention for the cost flow rates is that minus and plus signs indicate the resource and product cost flow rates, respectively. Erroneously, reverse sign convention was used in their study on the thermoeconomic analysis of ground-source heat pump systems . The lost cost flow rate due to the entropy generation appears as consumed cost in the evaporator; on the other hand, it appears as production cost in other components. The unit cost of heat delivered to the room or the unit cost of the cooling capacity is estimated to be 0.0344$/kWh by solving the four cost-balance equations given from Eqs. (53) to (59) with unit cost of electricity of 0.120 $/kWh. The unit cost of thermal and mechanical exergies and the irreversibility are CT = 0.1948, CP = 0.1636, and CS = 0.0187 $/kWh at the normal operation. It is noted that the unit cost of heat CH can be obtained from Eq. (60) directly with known values of CW, COP (β) and the ratio of the monetary flow rate of non-fuel items to the monetary flow rate of input (electricity). Table 9 confirms that cost-balance balance is satisfied for all components and the overall system.
|Evaporator+ Suction line||−1.85456||−0.12928||4.16420||−2.07016||−0.11020|
where β is the COP of the air conditioning units. Equation (60) provides the unit cost of cooling capacity as 0.0344 $/kWh with a unit cost of electricity of 0.120 $/kWh, β of 3.77, and a value of 0.081 for the ratio of the monetary flow rate of non-fuel items to the monetary flow rate of consumed electricity.
Explicit equations to obtain the unit cost of products from gas-turbine power plant and organic Rankin cycle plant operating by heat source as fuel and the unit cost heat for refrigeration system using the modified-productive structure analysis (MOPSA) method were obtained. MOPSA method provides two basic equations for exergy-costing method: one is a general exergy-balance equation and the other is cost-balance equation, which can be applicable to any components in power plant or refrigeration system. Exergy-balance equations can be obtained for each component and junction. The cost-balance equation corresponding to the exergy-balance equation can be obtained by assigning a unit cost to the principal product of each component. The overall exergy-costing equation to estimate the unit cost of product from the power plant and refrigeration system is obtained by summing up all the cost-balance equations for each component, junctions, and boundary of the system. However, one should solve the cost-balance equations for the components, junctions, and system boundary simultaneously to obtain the lost cost flow rate due to the entropy generation in each component. It should be noted that the lost work rate due to the entropy generation plays as “product” in the exergy-balance of the component, while the lost cost flow rate plays as “consumed resources” in the cost-balance equation. This concept is very important in the research area of thermoeconomic diagnosis [18, 26, 27, 28].
|C||unit cost of exergy ($/kJ)|
|Ci||initial investment cost ($)|
|CH||unit cost of heat ($/kWh)|
|Co||unit cost of fuel ($/kWh)|
|CS||unit cost of lost work due to the entropy generation ($/kWh)|
|CW||unit cost of electricity ($/kWh)|
|C ̇||monetary flow rate ($/h)|
|COP||coefficient of performance|
|CRF||capital recovery factor|
|ex||exergy per mass|
|E ̇ x||exergy flow rate (kW)|
|h||enthalpy per mass|
|H ̇||enthalpy flow rate (kW)|
|I ̇||irreversibility rate (kW)|
|m ̇||mass flow rate|
|PWF(i,n)||present worth factor|
|Q ̇ cv||heat transfer rate (kW)|
|S ̇||entropy flow rate (kW/K)|
|Sn||salvage value (KRW)|
|To||ambient temperature (°C)|
|W ̇ cv||work production rate (kW)|
|Z ̇ k||capital cost flow rate of unit k ($/h)|
|β||coefficient of performance|
|ϕ k||maintenance factor of unit k|
|W||work or electricity|
|W||work or electricity|