Next Article in Journal
Entropy Production in the Theory of Heat Conduction in Solids
Previous Article in Journal
iDoRNA: An Interacting Domain-based Tool for Designing RNA-RNA Interaction Systems
Previous Article in Special Issue
The Impact of Entropy Production and Emission Mitigation on Economic Growth
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exergy and Thermoeconomic Analysis for an Underground Train Station Air-Conditioning Cooling System

Department of Energy and Refrigerating Air-conditioning Engineering, National Taipei University of Technology, 1, Sec. 3, Zhong Xiao E. Road, Taipei 10608, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 25 November 2015 / Revised: 22 February 2016 / Accepted: 23 February 2016 / Published: 7 March 2016
(This article belongs to the Special Issue Entropy and the Economy)

Abstract

:
The necessity of air-conditioning causes the enormous energy use of underground train stations. Exergy and thermoeconomic analysis is applied to the annual operation of the air-conditioning system of a large underground train station in Taiwan. The current installation and the monitored data are taken to be the base case, which is then compared to three different optimized designs. The total revenue requirement levelized cost rate and the total exergy destruction rate are used to evaluate the merits. The results show that the cost optimization objective would obtain a lower total revenue requirement levelized cost rate, but at the expense of a higher total exergy destruction rate. Optimization of thermodynamic efficiency, however, leads to a lower total exergy destruction rate, but would increase the total revenue requirement levelized cost rate significantly. It has been shown that multi-objective optimization would result in a small marginal increase in total revenue requirement levelized cost rate, but achieve a significantly lower total exergy destruction rate. Results in terms of the normalized total revenue requirement levelized cost rate and the normalized total exergy destruction rate are also presented. It has been shown by second law analysis when applied to underground train stations that lower annual energy use and lower CO2 emissions can be achieved.

1. Introduction

The necessity of air-conditioning and also long operation hours is a major cause of the enormous amount of energy use for an underground train station. In this study, exergy and thermoeconomic analysis is applied to an underground train station that was constructed in Taiwan in recent years. The current status of the station is taken as the base case. Exergy indicates both the quality and the quantity in the energy conversion. Exergy destruction in reverse is the indicator of energy loss in terms of quality and quantity. Exergy is in fact indicative of thermodynamic second law efficiency. For the practical concerns of cost and economy, the total revenue requirement levelized cost rate can be used to evaluate the merit [1].
Some research literature can be found in the field of exergy and thermoeconomic analysis. However, only some is on the optimization of air-conditioning systems. Moreover, in some of the literature, only one of the exergetic or thermoeconomic criteria is considered. The exergetic criteria were discussed in many articles. Kodala et al. [2] studied the irreversibility of the thermoeconomic performance of refrigerators and heat pumps. Kodala et al. [3] investigated the thermoeconomic optimization of a two-stage combined heat pump system. Yumrutas et al. [4] performed exergy analysis on vapor compression refrigeration systems. Rezaie et al. [5] used exergy analysis for residential buildings’ assessment. Lu et al. [6] studied the patterns of domestic energy consumption using exergy analysis. Marletta et al. [7] applied exergy analysis to assess different constructions of air-conditioning systems.
In the field of thermal systems, Yilmaz et al. [8] applied second law analysis to evaluate the performance of heat exchangers. Bejan [9] analyzed the optimized geometry of engineering flow in terms of thermodynamics. Bisio and Rubatto [10] studied different aspects of irreversible processes in closed and steady-state open systems. Chen et al. [11] studied the optimal performance of an irreversible absorption refrigerator. Wall and Gong [12] presented the concept of exergy in sustainable development. Abhilash et al. [13] applied exergy efficiency to evaluate the inlet air cooling system for a centrifugal compressor.
There are several studies on thermoeconomic criteria. Sciubba [14] presented the concept of extended exergy accounting for the design of thermal systems. Li et al. [15] studied the economic benefits of automated fault detection and diagnosis for rooftop air-conditioners. Sivasakthivel et al. [16] discussed the optimized operating parameters for ground source heat pump systems. Liu and Niu [17] developed an optimal design analysis for heat recovery devices in buildings. Piacentino and Talama [18] proposed an improved thermoeconomic approach to the increase the reliability of the refrigeration system. The above studies show that multi-objective optimization can be applied to deal with different and competing objectives.
Multi-objective optimization of energy systems has attracted some attention in the related fields of research in recent years. Multi-objective optimization was used by Sayyadi and Nejatolahi [19] for a cooling tower-assisted refrigeration system. The total revenue requirement method and the total exergy destruction of the system were applied in their study. Sayyadi and Nejatolahi [20] also carried out multi-objective optimization for a cooling tower-assisted ground source heat pump. Three optimized objectives were included in their study. Sun et al. [21] investigated the relationship between the economic cost and the reliability of city electric supply.
More published literature on multi-objective optimization can be found recently. Rosen [22] discussed the economic benefits of applying the concept of exergy. In the area of thermal systems, published literature on multi-objective optimization can be found for CO2/NH3 refrigeration, integrated molten carbonate fuel cells, the organic Rankine cycle, combined gas turbines and solar energy [23,24,25,26,27].
Performance models of chillers and pumps for air-conditioning systems can be found in the research literature [28,29,30,31]. For dynamic energy saving control, Ma et al. [32] presented a real-time supervisory control strategy for building cooling water systems.
In the area of increasing the energy efficiency of air-conditioning, Du et al. [33] present a control-perfect index based on exergy analysis. They use the index to evaluate the control of an air-conditioning system. A novel cascade air-conditioning system of a 100-kW capacity of a sample building was proposed by Hojjat Mohammadi and Ameri [34]. Zagoet et al. [35] analyzed different combinations of heating systems and evaluated the primary energy efficiency for residential buildings in northern Italy.
In summary, most of the research on exergy and thermoeconomic analysis for air-conditioning systems has been of simple systems and at fixed state conditions. However, for an underground station, the air-conditioning operation is dynamic, also with the variation in ambient temperature. The energy use would be enormous when the underground train station is situated in a subtropical region, like Taiwan. Design analysis to balance the construction cost and energy efficiency is of upmost importance. This study considers the operation and maintenance cost due to the selection of systems that were overlooked in most of the published literature. In this study, the objective functions are the total revenue requirement levelized cost rate of the system C ˙ T R R L and the total exergy destruction rate I ˙ t o t . The monitored operation of the present system was taken to be the base case, against different cases of optimization. The results presented are based on an actual underground train station.

2. Air-Conditioning System and Methods

2.1. The Underground Train Station and Its Air-Conditioning System

The underground train station is located in Taipei, Taiwan, which began operation in 2011. The station building has dimensions of 330 m in length, 24 m in width and 21 m underground. The floor plan of the station includes the U-2F (2 floors underground) platform area, the U-1F (1 floor underground) concourse area, the G + 1F (ground level) lobby and retail and the G + 2F (second floor) retail store area. The G + 2F retail store area has a separate air-conditioning system and is metered separately; therefore, it was not included in this study. The floor plan of the station and the different zones served by a general and a 24-h system are shown in Figure 1. Zones enclosed by blue lines were served by the general system, and those enclosed by green lines were served by the 24-h system.
The 24-h air-conditioning system serves specific areas, including the central control room, the telecom room, the UPS room, the main power station, the power relay room and also the ticket area. The rest of the station is served by the general air-conditioning system, such as the G + 1F lobby and retail area, the U-1F concourse area and the U-2F platform area, as shown in Figure 2.

2.2. Energy Exergy and Thermoeconomic Analysis

The hourly monitored data used for the base case (BC) were collected from the building automation system (BAS), supplemented by information from the equipment suppliers, the Central Weather Bureau of Taiwan and also the local power company for the year 2012.

2.2.1. Energy Analysis

The cooling capacity Q ˙ (kW) of a chiller can be calculated as:
Q ˙ = m ˙ w e c p ( T e i T e o )
In Equation (1), m ˙ w e is the chilled water flowrate. T e i   and T e o are the chilled water temperatures flowing in and out through the evaporator, respectively. The specific heat capacity of water C P = 4.186 kJ·kg−1·K−1. Figure 2 shows that there is more than one chiller for each system. The average coefficient of performance C O P a v g of the chillers for either the general or the 24-h systems can be calculated using Equation (2):
C O P avg = 1 j 1 k [ m ˙ w e , k j c p ( T e i , k j T e o , k j ) ] / 1 j 1 k W ˙ a c , k j
In Equation (2), k is the number of chillers, and j is each of the operating hour in a year. W ˙ a c , k j is the compressor of the acting chiller in the hour. In times of lower cooling demand, only one chiller would be operating.
The operating power (kW) of a pump, either for a chilled water pump or a condensing water pump, can be calculated by using Equation (3):
W ˙ p u = ν ˙ Δ p η
where ν ˙ is the volume flow rate m3·s−1 and Δ p is the pressure jump (Pa) though the pump. The total power consumption of an air-conditioning system can be calculated as in Equation (4):
W ˙ tot = W ˙ p u , c h + W ˙ p u , c d + W ˙ A H U , F C + W ˙ c h + W ˙ C T
In Equation (4), W ˙ p u , c h , W ˙ p u , c d , W ˙ A H U , F C , W ˙ c h and W ˙ C T are respectively the power load of chilled water pumps, condensing water pumps, air handling units (AHU) or fan coils (FC), chillers and cooling towers. In this study, a system performance factor (SPF) for the entire system can be expressed as Equation (5):
S P F = Q ˙ e W ˙ t o t
Q ˙ e is the cooling capacity of the air-conditioning system.

2.2.2. Exergy Analysis

Kummel [36] relates the second law to economics. Exergy can be taken as the flow availability of a system that performs the maximum reversible work. Regarding the scope of this study, the kinetic and the potential work can be neglected. Hence, exergy can be expressed as H − H0 − T0(S − S0), where H is the enthalpy and S is the entropy of the respective states. T0(S − S0) is the entropy generation term or exergy destruction against the dead state. Sayyaadi and Neiatolahi [19] stated that a more convenient form of the exergy of a flowing system consists of a temperature component and a pressure component Ex = ExT + ExP. In this study the pressure component is neglected. The exergy of an air-conditioning system can be defined to be the maximum useful work attainable from a heat transfer process. Lu and Wu [6] mentioned that in differential form, exergy can be expressed as in Equation (6):
d E x = d E x Δ T = d E ˙ x = ( T 0 T T ) d Q ˙ e = ( T 0 T T ) m ˙ C p d T
Integrating from the given state ( T e i , T e o ) to the dead state ( T 0 ), which is the ambient condition in this study, a general equation for calculating exergy rate E ˙ x is obtained as in Equation (7):
E ˙ x = Q ˙ e ( T 0 T e i T e o ln T e i T e o 1 )
In Equation (7), Q ˙ e is the refrigeration effect (kW), T0 is the outdoor temperature (the reference environment temperature), Tei is inlet chiller water temperature and Teo is the outlet chiller water temperature. The thermodynamic second law efficiency (exergy efficiency) can be expressed as Equation (8):
η Π = E ˙ x W ˙ t o t
In this study, the total exergy destruction rate is the sum of the general and the 24-h systems, as expressed in Equation (9):
I ˙ t o t = I ˙ t o t , 24 h + I ˙ t o t , g e n
where I ˙ t o t , 24 h = E ˙ x 24 h + W ˙ c h , 24 h + W ˙ p u , 24 h + W ˙ A H U & F C , 24 h + W ˙ C T , 24 h and I ˙ t o t , g e n = E ˙ x g e n + W ˙ c h , g e n + W ˙ p u , g e n + W ˙ A H U & F C , g e n + W ˙ C T , g e n . I ˙ t o t , I ˙ t o t , 24 h and I ˙ t o t , g e n are the exergy destruction rate for the total, the 24-h system and the general air-conditioning system, respectively.

2.2.3. Thermoeconomic Analysis

The thermoeconomic analysis has to consider the construction cost, amortization, maintenance and electricity consumption. In this study, the total revenue requirement levelization (TRRL) method was applied. In this study, the TRRL method takes into account the estimated total carrying charge, along with assumptions for economic, financial and market input parameters calculated on a yearly basis. Moreover, the non-uniform annual monetary values of the carrying charges, maintenance cost and electricity cost are levelized and converted to an equivalent series of constant payments (annuities) [37].
TRRL can be computed as in Equation (10):
T R R L = C R F 1 n T R R j ( 1 + i e f f ) j
In Equation (10), CRF is the capital recovery factor. T R R j is the total revenue requirement in the j-th year of system operation, and i e f f is the interest rate (%). In applying Equation (10), each monetary transaction is assumed to occur at year end. The capital recovery factor C R F is calculated by Equation (11):
C R F = i e f f ( 1 + i e f f ) n ( 1 + i e f f ) n 1
In Equation (11), n is the economic lifetime in years. The series of payments of the annual electrical cost E C j can be calculated using Equation (12):
E C j = E C 0 ( 1 + r E C ) j
In Equation (12), r E C is the constant escalation. The levelized value E C L can then be determined multiplying the first year electricity cost E C 0 (NT) by the constant escalation levelization factor (CELF) as Equation (13):
E C L = E C 0 C E L F = E C 0 k E C ( 1 k E C n ) ( 1 k E C ) C R F
In Equation (13), k EC = ( 1 + r EC ) / ( 1 + i eff ) . The terms r E C and C R F are respectively the annual escalation rate of electricity cost and the capital recovery factor. E C 0 can be computed by Equation (14):
E C 0 = C e l e c t ( E p u , c h + E p u , c d + E A H U , F C + E c h + E C T )
In Equation (14), Celect is the electricity price per kWh in Taiwan taken as 3.1 NT/kWh. The total annual electric consumption (kWh) includes that of the chiller water pump (Epu, ch), condensing water pumps (Epu, cd), air side equipment (EAHU, FC), such as air handling units (AHU) and fan coils (FC), water chillers (Ech) and cooling towers (ECT).
The levelized annual operating and maintenance costs OMCL are given in Equation (15):
O M C L = O M C 0 C E L F = O M C 0 k O M C ( 1 k O M C n ) ( 1 k O M C ) C R F
where k O M C = 1 + r O M C 1 + i e f f , and the term r O M C is a constant representing the nominal escalation rate for the operating and maintenance costs.
The economic operating lifetime of the system is taken to be 15 years. Therefore, the magnitude of the economic constant, such as i e f f , r O M C and r E C are taken to be 0.02, 0.03 and 0.03, respectively, in the analysis.
For the air-conditioning system, the annual total revenue requirement levelization (TRRL) is equal to the sum of the carrying charge levelization (CCL), electricity costs levelization (ECL) and operating and maintenance cost levelization (OML), as shown in Equation (16):
T R R L = C C L + E C L + O M C L
The levelized cost rate of the total revenue requirement can be calculated by Equation (17):
C ˙ TRRL = C ˙ C C L 24 h + C ˙ E C L 24 h + C ˙ C C L g e n + C ˙ E C L g e n + C ˙ O M L = C C L 24 h m 24 h + E C L 24 h m 24 h + C C L g e n m g e n + E C L g e n m g e n + O M L m O M
In Equation (17), C ˙ TRRL is the hourly cost of the total revenue requirement levelization. The hourly rates of CCL and ECL are to be calculated separately for the general and the 24-h systems. The annual operating hours of the 24-h and the general systems, m 24 h and m g e n , are 8784 h and 4870 h, respectively, according to the actual operation of the train station.
In this study, the reduction of carbon dioxide (CO2) emissions due to energy saving was estimated. The emission of carbon dioxide per kWh of electricity used given by Taipower [38] for the year of 2012 was 0.522 kg/kWh.

2.2.4. Empirically-Based Performance Models of Equipment

The energy performance models of chillers refer to Lee et al. [29] and Jiang and Reddy [30]. The regression models of the chillers and pumps are fitted using data provided by equipment suppliers. These models are used in thermoeconomic and exergy analysis in this study. Prediction of chiller COP refers to the modified Gordon-Ng universal model proposed by Jiang and Reddy [30]. The model considers the refrigerant flow rate that may change the internal entropy production in the compressor. The independent variables of chiller COP are then the cooling load ratio Q ˙ e / Q ˙ e ( m a x ) , evaporator outlet temperature T e o and condenser inlet water temperature T c i . The functional form of the model is shown in Equation (18):
T e o T c i ( 1 + 1 C O P ) 1 = ( β 1 + β 2 Q ˙ e Q ˙ e ( max ) ) T e o Q ˙ e + β 3 T c i T e o T c i Q ˙ e + β 4 Q ˙ e T c i ( 1 + 1 C O P )
In Equation (18), the regression constants β1–β4 are determined by regressing the performance data provided by equipment suppliers.
Similarly, the power consumption of chilled water pumps and condensing water pumps W ˙ p u can be expressed in regression form, as shown in Equation (19):
W ˙ p u = α 0 + α 1 V ˙ + α 2 V ˙ 2 + α 3 V ˙ 3
The independent variable is only the volumetric flow rate V ˙ , whether of the chilled water or the condensing water. The regression constants α1–α3 are determined by regressing the performance data provided by equipment suppliers.

2.2.5. Multi-Objective Optimization

The two objective functions studied are the economic and thermodynamic objective functions shown in Equation (20) and Equation (21):
C ˙ TRRL = C ˙ T R R L C ˙ T R R L , min C ˙ T R R L , max C ˙ T R R L , min
I ˙ tot = I ˙ t o t I ˙ t o t , min I ˙ t o t , max I ˙ t o t , min
The multi-objective optimization refers to the method of the Pareto optimal frontier [39]. In Equations (20) and (21), C ˙ TRRL , min and C ˙ TRRL , max are the minimum and maximum values of C ˙ TRRL in the Pareto frontier, similarly for I ˙ t o t , min and I ˙ t o t , max . C ˙ TRRL = C ˙ TRRL , min and I ˙ tot = I ˙ t o t , min correspond to the limits of economic and thermodynamic optimization. When Equations (20) and (21) are separately used in single objective optimization, C ˙ TRRL and I ˙ tot would lies between 0 and 1.
In this study, the decision variables of selecting chillers are:
(1)
Using screw or centrifugal compressors, constant or variable speed drive.
(2)
The operating refrigeration capacity.
(3)
The coefficient of performance (COP).
(4)
Inlet/outlet chilled water and condensing water temperatures.
(5)
The cost of chillers.
Similarly, the decision variables of selecting pumps are:
(1)
Pumps with constant speed drive or variable speed drive.
(2)
The efficiency of pumps.
(3)
The pump pressure and water volumetric flow rate (m−3·s−1).
(4)
The cost of pumps.
Three optimization cases are compared to the base case, with the operating parameters described in Table 1 and Table 2. The three cases are namely cost consideration (CC), thermodynamic efficiency (TE) and multiple objectives of efficiency and cost (MO). The operation constraints of the 24-h and the general air-conditioning systems are listed in Table 1 and Table 2. Constraints of these decision variables are partly the limitations emanating from the technical data of equipment suppliers. For the general system, the total chilled water flow rate is 4032 kg/min (equivalent to liters per minute (lpm) for water) and is 1210 kg/min for the 24-h system. It is noted that the variable frequency drive (VFD) is applied to the chillers and pumps in Case TE and Case MO. VFD systems general cost more, but it would better match the cooling demand and the overall operation of the system.

3. Results

3.1. Electricity Consumption and Electricity Cost of the Underground Train Station

The electricity use of an underground train station is vastly different from a commercial building. The train station studied would begin operation at 4:30 A.M. with the general air-conditioning system started at about 6:30 A.M. The hourly monitored power use of the station (the base case (BC)) is presented in Figure 3. Figure 3 shows that for summer months, the peak power occurs in the afternoon of the day. The monthly electricity cost of the underground station is presented in Figure 4, which clearly indicates that the peak power use of air-conditioning happens in the summer months. It is then obvious that air-conditioning constituted a significant part of the total electricity cost of the underground train station and is also the major cause of the peak power load.

3.2. Refrigerating Energy and Exergy of Underground Train Station

The hourly refrigerating capacity of the air-conditioning system is calculated using Equation (1) for each operating hour. The daily cooling load is shown in Figure 5. The 24-h air-conditioning system appears to have an almost constant load through the year. The general system has high peak demand in summer months. The daily refrigeration exergy is determined from Equation (7). The computed results are shown in Figure 6. The chilled water temperature Te in Equation (6) is from the monitored data of the station. The ambient temperature T0 was obtained from the Central Weather Bureau of Taiwan [40].
Comparing the results in Figure 5 and Figure 6, it is noted that exergy is about one tenth of the refrigeration effect indicating that the ideal COP of a chiller would be larger than 10, significantly higher than commercial products. The chiller’s performance used in this study is presented in Table 1 and Table 2.
The monthly average COP of the chillers computed is shown in Figure 7. Chillers of the general air-conditioning system are shown to have higher COP. Larger capacity centrifugal chillers normally are more energy efficient. However, the average COP is significantly less than that obtained from the above exergy analysis. It is also noted that Case TE has the highest average COP followed by Case MO, with Case CC having the lowest average COP.
The annual energy consumption of pumps calculated from the performance models is shown in Table 3 and Table 4. The performance models for each of the optimization cases are based on the pump energy efficiency shown in Table 1 and Table 2. Case TE has the lowest pump energy consumption followed by Case MO. Case CC has the highest energy consumption as lower energy efficiency pumps are used.
Table 5 shows the annual energy use of the air-conditioning systems. Pump energy is the second largest item of energy use. This is the characteristics of an underground train station as long pipelines are needed to deliver chilled water and condensing water. The air side equipment, such as AHU and FC, are not optimized in this study, nor are the cooling towers. It is also noted in Table 5 that chillers make up about half of the total energy use.
The monthly average SPF and the second law efficiency of the air-conditioning system are shown in Figure 8. The results are computed using the monitored data of the underground train station in 2012. It can be seen in Figure 8 that SPF is higher than 2.5 for most of the months. The results indicate one unit of energy input would have more than 2.5 units of refrigeration effect. SPF is calculated using Equation (5). The second law efficiency of the system was evaluated using Equation (8).
The second law efficiency (efficiency of exergy) is evaluated against the maximum useful work. Therefore, the second law efficiency has a maximum value of 1.0. It can be seen in Figure 8 that the second law efficiency is lower than 0.3. In related works found in the research literature, the second law efficiency was evaluated for each component of the air-conditioning system and mostly calculated for ideal cycle operation. In this study, the actual operation of the entire system is evaluated for either SPF or the second law efficiency.

3.3. Multi-Objective Optimization

In this study, the multi-objective functions are the normalized total revenue requirement levelization C ˙ TRRL and the normalized exergy destruction I ˙ tot . The method of the Pareto optimal frontier [39] is applied to obtain the normalized C ˙ TRRL and I ˙ tot . The normalized C ˙ TRRL and I ˙ tot of the four cases are shown in Table 6.
Figure 9 presents the normalized form of the Pareto optimal frontier in the multi-objective optimization (MO). In Figure 9, multi-objective optimization is at the lowest point of the root mean squares of C ˙ TRRL and I ˙ tot . For Case MO optimal C ˙ TRRL is calculated to be about 0.27 and I ˙ tot at about 0.26. Figure 9 also shows that better cost benefits can be obtained when C ˙ TRRL is lower than 0.26 and better thermodynamics when C ˙ TRRL is higher than 0.26. It is noted in Table 6 that for Case MO, C ˙ TRRL equals 0.26, as compared to that of Case BC at 0.21. However, I ˙ tot of Case MO is significantly lower than that of Case BC.

3.4. Energy Exergy and Thermoeconomic Analysis

The results presented in Table 7 are reevaluated for an economic lifetime of 15 years. For the Case CC, in order to lower the initial construction cost, lower energy efficiency equipment is selected. Consequently, Case CC has the highest electricity cost, the lowest exergy efficiency, the lowest system performance factor (SPF) and the highest CO2 emission. At the same time, Case CC has the lowest carrying charge rate and the lowest total revenue requirement levelized rate. In comparison, when thermodynamic efficiency is the sole concern (Case TE), the highest system performance and exergy efficiency are obtained. However, Case TE has the highest total revenue requirement levelized rate. For the case of multi-objective optimization (Case MO), a balance of cost and efficiency is the concern. The results show that the construction cost of MO is slightly higher than the base case (BC), but would result in a system performance factor (SPF) and second law efficiency (exergy efficiency) next to Case TE. It is also noted that Case BC, compared to Case TE and Case MO, has higher electricity consumption, lower SPF and higher annual CO2 emission. It is known that Case BC was designed for energy saving. However, when exergy and thermodynamic analysis is applied for further analysis, BC is not an optimal design.
The comparison of the exergy destruction I ˙ t o t for the components of the air-conditioning system is shown in Figure 10. The chiller is the largest item of exergy destruction, followed by pumps and air side equipment. As the air side components are not included in the analysis, the exergy destruction of AHU and FC is the same for the four cases. The main reason that air side equipment is not included in the evaluation is due to the uncertainty in the monitored data. When the entire system is concerned, Case CC has exergy destruction significantly higher than all other cases. When cost and energy saving are both considered, Case MO is superior to Case BC in terms of exergy destruction.
The comparison of the levelized cost rate is shown in Figure 11. Case CC has the lowest carrying charge cost rate due to the lower construction cost, but has significantly higher electricity cost over the other cases. When the total revenue requirement levelized cost rate is compared, Case CC is less than Case BC by 33 NT/h. Case MO has a total revenue requirement levelized cost rate slightly higher than Case BC, but has a lower electricity cost rate. In comparison, Case TE has the highest total revenue requirement levelized cost rate.
Figure 12 shows the percentage comparison of construction cost, total exergy destruction rate, total revenue requirement levelized cost rate and CO2 emission. Case BC is taken to be 100% in all of the items. The construction cost of Case CC is 95.7%, but is 99.24% in terms of C ˙ T R R L . The result shows that Case CC may not be a good selection. Case TE is the best selection in terms of reducing CO2 emission and low exergy destruction. However, Case TE has C ˙ T R R L 2.76% higher than that of Case BC. The comparison shows that Case TE is a good selection when reducing CO2 emission is a prime concern. Case MO has C ˙ T R R L almost as an equal percentage as for Case BC, but has lower exergy destruction and lower CO2 emission. Case MO appears to be a superior selection when economics and second law efficiency are both a concern.

4. Discussion

The above results show that applying exergy and thermoeconomic analysis to the air-conditioning system design would be of practical importance. It is also noted that when the total revenue requirement levelized cost rate is evaluated, a lower construction cost may not result in good economics. Optimization of multiple objectives would provide a balanced consideration of both energy efficiency and investment return. Case MO has a total revenue requirement levelized cost rate comparable to Case BC, but performs better in reducing CO2 emission and low exergy destruction.
The results obtained also show that both the system performance factor (SPF) and second law efficiency can be applied to the entire air-conditioning system. The results of exergy efficiency analysis as given in Table 7 show that system exergy efficiency is lower than 0.2 in all four cases. The results indicate high potential for further improvement in equipment efficiency. SPF is an effective indicator of the system performance. For Case TE, SPF equals about 3.0; this means one unit of electrical energy input would result in about three units of refrigeration effect in energy terms.

5. Conclusions

It has been shown that thermoeconomic and exergy analysis can be used to improve the energy performance and cost savings of the air-conditioning system of an underground train station. It is also shown that lower initial construction cost may have only a marginal benefit when the total revenue requirement levelized rate is evaluated. Optimization of thermodynamic efficiency would bring the benefits of a lower exergy destruction rate and lower energy use. The results also show that multi-objective optimization has the benefit of a lower total revenue requirement levelized rate, but performs well in reducing CO2 emission and low exergy destruction. The annual system exergy efficiency obtained is lower than 0.2, which points to high potential in equipment efficiency improvement. A system performance factor is introduced to evaluate the efficiency of the entire air-conditioning system. In summary, this study shows that thermoeconomic and exergy analysis can be applied in an annual analysis, so as to achieve optimized hourly operation of an air-conditioning system in general and for an underground train station in particular. Moreover, it has been shown that a higher system performance factor (SPF) can be achieved along with lower annual CO2 emission.

Acknowledgments

This research is supported by the Ministry of Science and Technology of Taiwan with research Project MOST 103-2221-E-027-110, and this is gratefully acknowledged hereby. Comments from the reviewers and the editorial office help to improve the paper and are also acknowledged hereby.

Author Contributions

Ke Yang Liao performed the analysis. Yew Khoy Chuah contributed in the presentation of the results and partly the writing of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

CCLCarrying charge levelization (NT/y)
CELFConstant escalation levelization factor
COPCoefficient of performance
CRFCapital recovery factor
C ˙ T R R L Total revenue requirement levelized cost rate (NT/h)
C ˙ TRRL Normalized form of C ˙ T R R L
ECElectrical cost (NT)
ExExergy (kWh)
HEnthalpy (kWh)
I ˙ t o t Total exergy destruction rate (kW)
I ˙ tot Normalized form of I ˙ t o t
KECThe constant of electricity cost (Equation (13))
NTNew Taiwan dollars, 1 USD = 33 NT
OMCOperating and maintenance cost (NT/y)
Q ˙ Cooling capacity (kW)
QHeat transfer (kWh)
SEntropy (kWh/K)
SPFSystem performance factor
UInternal energy (kWh)
TCRTotal capital recovery
TRRLTotal revenue requirement levelization
W ˙ Power (kW)
Greek symbols
η Efficiency
η Π The thermodynamic second law efficiency (the efficiency of exergy)
ν ˙ The volume flowrate (m3·s−1)
Subscripts and superscripts
acActing
AHUAir handling unit
cdCondenser
CTCooling tower
eEvaporator
electElectrical
FCFan coil
iInlet
jThe operating hour in a year of chillers (h)
kThe number of chillers
mAnnual operating time (h)
nEconomic life time (y)
oOutlet
puPump
totTotal

References

  1. Fu, Q.S. Thermodynamic Analysis Method in Energy System; Xi’an Jiaotong University Press: Xi’an, China, 2005; Chapter 15. [Google Scholar]
  2. Kodal, A.; Sahin, B.; Yilmaz, T. Effects of internal irreversibility and heat leakage on the finite time thermoeconomic performance of refrigerators and heat pumps. Energy Convers. Manag. 2000, 41, 607–619. [Google Scholar] [CrossRef]
  3. Kodal, A.; Sahin, B.; Oktem, A.S. Performance analysis of two stage combined heat pump system based on thermoeconomic optimization criterion. Energy Convers. Manag. 2000, 41, 1989–1998. [Google Scholar] [CrossRef]
  4. Yumrutas, R.; Kunduz, M.; Kanoglu, M. Exergy analysis of vapor compression refrigeration systems. Exergy Int. J. 2002, 2, 266–272. [Google Scholar] [CrossRef]
  5. Rezaie, B.; Dincer, I.; Esmailzadeh, E. Exergy options for residential buildings assessment. Energy Convers. Manag. 2013, 65, 637–647. [Google Scholar] [CrossRef]
  6. Lu, S.; Wu, J.Y. Optiomal selection among different domestic energy consumption patterns based on energy and exergy analysis. Energy Convers. Manag. 2010, 51, 1398–1406. [Google Scholar] [CrossRef]
  7. Marletta, L. Air condition systems from a 2nd law perspective. Entropy 2010, 12, 859–877. [Google Scholar] [CrossRef]
  8. Yilmaz, M.; Sara, O.N.; Karsli, S. Performance evaluation criteria for heat exchangers based on second law analysis. Exergy Int. J. 2001, 1, 278–294. [Google Scholar] [CrossRef]
  9. Bejan, A. Thermodynamic optimization of geometry in engineering flow systems. Exergy Int. J. 2001, 1, 269–277. [Google Scholar] [CrossRef]
  10. Bisio, G.; Rubatto, G. Work and production aspects of irreversible processes in closed and steady-state open systems. Exergy Int. J. 2002, 2, 192–205. [Google Scholar] [CrossRef]
  11. Chen, L.G.; Li, Y.; Sun, F.R.; Wu, C. Optimal performance of an irreversible absorption refrigerator. Exergy Int. J. 2002, 2, 167–172. [Google Scholar] [CrossRef]
  12. Wall, G.; Gong, M. On exergy and sustainable development-Part 1: Conditions and concepts. Exergy Int. J. 2001, 1, 128–145. [Google Scholar] [CrossRef]
  13. Abhilash, S.; Pradeep, A.; Gyuwan, K.; Dong, K.H. Thermodynamic analysis of inlet air cooling system for a centrifugal compressor. In Proceedings of the 7th International Exergy, Energy and Environment Symposium, Valenciennes, France, 27–30 April 2015. C12.
  14. Sciubba, E. Beyond thermoeconomics? The concept of Extended Exergy Accounting and its application to the analysis and design of thermal systems. Exergy Int. J. 2001, 1, 68–84. [Google Scholar] [CrossRef]
  15. Li, H.R.; Braun, J.E. Economic Evaluation of Benefits Associated with Automated Fault Detectionand Diagnosis in Rooftop Air Conditioners. ASHRAE Trans. 2007, 113, 200–210. [Google Scholar]
  16. Sivasakthivel, T.; Murugesan, K.; Thomas, H.R. Optimization of operating parameters of ground source heat pump system for space heating and cooling by Taguchi method and utility concept. Appl. Energy 2014, 116, 76–85. [Google Scholar] [CrossRef]
  17. Liu, X.P.; Niu, J.L. An optimal design analysis method for heat recovery devices in building Applications. Appl. Energy 2014, 129, 364–372. [Google Scholar] [CrossRef]
  18. Piacentino, A.; Talamo, M. Innovative thermoeconomic diagnosis of multiple faults in air conditioning units: Methodological improvements and increased reliability of results. Int. J. Refrig. 2013, 36, 2343–2365. [Google Scholar] [CrossRef]
  19. Sayyaadi, H.; Neiatolahi, M. Multi-objective optimization of a cooling tower assisted vapor compression refrigeration. Int. J. Refrig. 2013, 34, 243–256. [Google Scholar] [CrossRef]
  20. Sayyaadi, H.; Neiatolahi, M. Multi-objective optimization of a cooling tower assisted ground source heat pump. Geothermics 2011, 40, 221–232. [Google Scholar]
  21. Sun, T.Q.; Wang, X.H.; Ma, X.G. Relationship between the economic cost and the reliability of the electric power supply system in city: A case in Shanghai of China. Appl. Energy 2009, 86, 2262–2267. [Google Scholar] [CrossRef]
  22. Rosen, M.A. Exergy and economics: Is exergy profitable? Exergy Int. J. 2002, 2, 218–220. [Google Scholar] [CrossRef]
  23. Aminyavari, M.; Najafi, B.; Shirazi, A.; Rinaldi, F. Exergetic, economic and environmental (3E) analyses, and multi-objective optimization of a CO2/NH3 refrigeration system. Appl. Energy 2014, 65, 42–50. [Google Scholar] [CrossRef]
  24. Mamaghani, A.H.; Najafi, B.; Shirazi, A.; Rinaldi, F. 4E analysis and multi-objective optimization of an integrated MCFC (molten carbonate fuel cell) and ORC (organic Rankine cycle) system. Energy 2015, 82, 650–663. [Google Scholar] [CrossRef]
  25. Mansureh, K.; Rahim, K.S.; Keyvan, B. Muti-objective optimization of power and heat in a combined gas turbine and organic Rankine cycle. In Proceedings of the 7th International Exergy, Energy and Environment Symposium, Valenciennes, France, 27–30 April 2015. C15.
  26. Keyvan, B.; Rahim, K.S. Thermoeconomic muti-objective optimization of an ammonia-water power/cooling cycle coupled with a HCCI engine. In Proceedings of the 7th International Exergy, Energy and Environment Symposium, Valenciennes, France, 27–30 April 2015. C30.
  27. Hanifi, K.; Javaherdeh, K.; Yari, M. Exergy and exergoeconomic analysis and optimization of the cogeneration cycle under solar radiation dynamic model using genetic algorithm. In Proceedings of the 7th International Exergy, Energy and Environment Symposium, Valenciennes, France, 27–30 April 2015. C41.
  28. Lee, K.P.; Cheng, T.A. A simulation–optimization approach for energy efficiency of chilled water system. Energy Build. 2012, 54, 290–296. [Google Scholar] [CrossRef]
  29. Lee, T.S.; Liao, K.Y.; Lu, W.C. Evaluation of the suitability of empirically-based models for predicting energy performance of centrifugal water chillers with variable chilled water flow. Appl. Energy 2012, 93, 583–595. [Google Scholar] [CrossRef]
  30. Jiang, W.; Reddy, T.A. Reevaluation of the Gordon-Ng performance models for water-cooled chillers. ASHRAE Trans. 2003, 109, 272–287. [Google Scholar]
  31. Tirmizi, S.A.; Gandhidasan, P.; Zubair, S.M. Performance analysis of a chilled water system with various pumping schemes. Appl. Energy 2012, 100, 238–248. [Google Scholar] [CrossRef]
  32. Ma, Z.J.; Wang, S.W.; Xu, X.H.; Xiao, F. A supervisory control strategy for building cooling water systems for practical and real time applications. Energy Convers. Manag. 2008, 49, 2324–2336. [Google Scholar] [CrossRef]
  33. Du, Z.; Jin, X.; Fan, B. Evaluation of operation and control in HVAC (heating, ventilation and air conditioning) system using exergy analysis method. Energy 2015, 89, 372–381. [Google Scholar] [CrossRef]
  34. Mohammadi, S.M.H.; Ameri, M. Energy and exergy comparison of a cascade air conditioning system using different cooling strategies. Int. J. Refrig. 2013, 41, 14–26. [Google Scholar] [CrossRef]
  35. Zago, M.; Casalegono, A.; Marchesi, R.; Rinaldi, F. Efficiency analysis of independent and centralized heating systems for residential buildings in northern Italy. Energies 2011, 4, 2115–2131. [Google Scholar] [CrossRef]
  36. Kummel, R. The Second Law of Economics: Energy, Entropy, and the Origins of Wealth; Springer: New York, NY, USA, 2011; pp. 100–102. [Google Scholar]
  37. Bejan, A.; Tsatsaronis, G.; Morran, M. Thermal Design and Optimization; Wiley: New York, NY, USA, 1996. [Google Scholar]
  38. Taiwan Power Company, Jan. to Dec. Taiwan Railway Administration Songshan Station Electricity Price List (1 January 2012–31 December 2012). Available online: http://www.taipower.com.tw/ (accessed on 29 March 2014).
  39. Konak, A.; Coit, D.W.; Smith, A.E. Multi-objective optimization using genetic algorithms: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 992–1007. [Google Scholar] [CrossRef]
  40. Central Weather Bureau. Taipei Outdoor Temperature (1 January 2012–31 December 2012). Available online: http://www.cwb.gov.tw/V7/observe/real/46692.htm (accessed on 24 April 2014).
Figure 1. Station floor plan and zones served by a general and a 24-h AC system. (a) The ground floor plan; (b) the underground U-1 floor plan; (c) the underground U-2 floor plan.
Figure 1. Station floor plan and zones served by a general and a 24-h AC system. (a) The ground floor plan; (b) the underground U-1 floor plan; (c) the underground U-2 floor plan.
Entropy 18 00086 g001
Figure 2. The schematics of the 24-h and the general air-conditioning systems. (a) The schematics of the 24-h air-conditioning system; (b) the schematics of the general air-conditioning system. AHU, air handling unit; FC, fan coils.
Figure 2. The schematics of the 24-h and the general air-conditioning systems. (a) The schematics of the 24-h air-conditioning system; (b) the schematics of the general air-conditioning system. AHU, air handling unit; FC, fan coils.
Entropy 18 00086 g002
Figure 3. Hourly electricity consumption of the underground train station in 2012.
Figure 3. Hourly electricity consumption of the underground train station in 2012.
Entropy 18 00086 g003
Figure 4. Monthly electricity cost of the underground train station in 2012.
Figure 4. Monthly electricity cost of the underground train station in 2012.
Entropy 18 00086 g004
Figure 5. Daily cooling load of the underground train station under study.
Figure 5. Daily cooling load of the underground train station under study.
Entropy 18 00086 g005
Figure 6. Daily refrigerating exergy of the underground train station under study in 2012.
Figure 6. Daily refrigerating exergy of the underground train station under study in 2012.
Entropy 18 00086 g006
Figure 7. The monthly average COP of the chillers. (a) The 24-h air-conditioning system; (b) the general air-conditioning system.
Figure 7. The monthly average COP of the chillers. (a) The 24-h air-conditioning system; (b) the general air-conditioning system.
Entropy 18 00086 g007
Figure 8. The monthly system performance factor and second law efficiency.
Figure 8. The monthly system performance factor and second law efficiency.
Entropy 18 00086 g008
Figure 9. The normalized Pareto optimal frontier.
Figure 9. The normalized Pareto optimal frontier.
Entropy 18 00086 g009
Figure 10. Exergy destruction of the air-conditioning system.
Figure 10. Exergy destruction of the air-conditioning system.
Entropy 18 00086 g010
Figure 11. Levelized costs rate for the four cases.
Figure 11. Levelized costs rate for the four cases.
Entropy 18 00086 g011
Figure 12. Percentage comparison of the four cases.
Figure 12. Percentage comparison of the four cases.
Entropy 18 00086 g012
Table 1. Operating range of the 24-h air-conditioning system. BC, base case; CC, cost consideration; TE, thermodynamic efficiency; MO, multiple objectives.
Table 1. Operating range of the 24-h air-conditioning system. BC, base case; CC, cost consideration; TE, thermodynamic efficiency; MO, multiple objectives.
ItemBCCCTEMO
1.  422-kW screw chillerChiller controlCFD 1CFDVFD 2VFD
COP4.16–5.123.66–4.304.51–5.354.36–5.32
COPavg4.523.885.134.91
2.  Zone chiller pump water head (kPa): 440Pump controlCFDCFDVFDVFD
Flowrate (lpm)12101210726–1210726–1210
Efficiency η0.750.70.72–0.750.71–0.73
3.  Condenser pump water head (kPa): 350Pump controlCFDCFDVFDVFD
Flowrate (lpm)15121512907–1512907–1512
Efficiency η0.750.70.72–0.750.71–0.73
Chiller operating at Tci (°C): 26–30; Tco (°C): 29–37; Tei (°C): 9.8–12.5; Teo (°C): 6.5–7.5
1 CFD = constant frequency drive; 2 VFD = variable frequency drive.
Table 2. Operating range of the general air-conditioning system.
Table 2. Operating range of the general air-conditioning system.
ItemBCCCTEMO
1.  2110-kW centrifugal chillersChiller controlCFD 1CFDVFD 2CFD
COP4.6–6.54.52–5.705.12–7.274.80–6.67
COPavg5.845.246.395.89
2.  Primary chiller pumps water head (kPa): 190Pump controlCFDCFDVFDVFD
Flowrate (lpm)403240322419–40322419–4032
Efficiency η0.770.720.73–0.770.72–0.75
3.  Zone 1 chiller pumps water head (kPa): 370Pump controlVFDCFDVFDVFD
Flowrate (lpm)2612–435443542612–43542612–4354
Efficiency η0.810.760.77–0.810.75–0.80
4.  Zone 2 chiller pump water head (kPa): 270Pump controlVFDCFDVFDVFD
Flowrate(lpm)1019–169316931019–16931019–1693
Efficiency η0.72–0.750.70.72–0.750.71–0.74
5.  Condenser pumps water head (kPa): 385Pump controlCFDCFDVFDVFD
Flowrate (lpm)756075604356–75604356–7560
Efficiency η0.840.790.81–0.840.80–0.83
Chiller operating at Tci (°C): 26–30; Tco (°C): 29–37; Tei (°C): 9.1–12.5; Teo (°C): 6.5–15.7
1 CFD = constant frequency drive, 2 VFD = variable frequency drive.
Table 3. Energy consumption of the pumps for the 24-h system.
Table 3. Energy consumption of the pumps for the 24-h system.
ItemBCCCTEMO
Condenser water pump (kWh)86,74893,03463,74671,654
Chilled water pump (kWh)108,255116,33366,91680,259
Total (kWh)195,003209,367130,662151,912
Table 4. Energy consumption of pumps for the general system.
Table 4. Energy consumption of pumps for the general system.
ItemBCCCTEMO
Condenser pump (kWh)417,128443,867368,903389,143
Primary chilled water pump (kWh)118,344126,563105,474111,168
Zone 1 Pump (kWh)216,224248,865216,224225,037
Zone 2 Pump (kWh)63,94177,67663,94167,649
Total (kWh)815,637896,971754,542792,998
Table 5. Annual energy use of the major components of the air-conditioning system.
Table 5. Annual energy use of the major components of the air-conditioning system.
Electric Energy Use Items (kWh)BCCCTEMO
24-h system chillers521,201607,417480,118527,682
General system chillers1,601,7991,785,0791,464,4791,588,376
24-h pumps195,003209,367130,662151,912
General system pumps815,637896,971754,542792,998
AHU and FC949,551949,551949,551949,551
Cooling towers126,607126,607126,607126,607
Total system4,209,7984,574,9913,905,9594,137,126
Table 6. C ˙ TRRL and I ˙ t o t and their optimized normalized values.
Table 6. C ˙ TRRL and I ˙ t o t and their optimized normalized values.
Item C ˙ TRRL I ˙ tot C ˙ TRRL I ˙ tot
Base case (BC)43716460.210.44
Cost consideration (CC)433871201
Thermodynamic efficiency (TE)449259410
Multiple objectives, cost and efficiency (MO)43786260.260.27
Table 7. Energy exergy and thermoeconomic analysis of the air-conditioning system.
Table 7. Energy exergy and thermoeconomic analysis of the air-conditioning system.
ItemBCCCTEMO
Air-conditioning system construction cost (NT)136,001,250130,156,950144,767,700138,527,235
Electricity levelized cost rate (NT/h)1797195716701733
Maintenance levelized cost rate (NT/h)562455680517
Carrying charge levelized cost rate (NT/h)2012192621422058
Total revenue requirement levelized cost rate (NT/h)4371433844924378
Total exergy destruction rate (kW)647712594626
Total refrigeration effect (kWh/y)11,719,26511,719,26511,719,26511,719,265
Total refrigeration exergy (kWh/y)722,226717,814722,925722,627
Total electricity consumption(kWh/y)4,209,7984,574,9913,905,9594,094,428
CO2 emission (kg/y)2,197,5142,388,1452,038,9102,137,291
System performance factor (SPF)2.782.563.002.83
Exergy efficiency (thermal second-law efficiency)0.1720.1570.1850.175

Share and Cite

MDPI and ACS Style

Liao, K.Y.; Chuah, Y.K. Exergy and Thermoeconomic Analysis for an Underground Train Station Air-Conditioning Cooling System. Entropy 2016, 18, 86. https://0-doi-org.brum.beds.ac.uk/10.3390/e18030086

AMA Style

Liao KY, Chuah YK. Exergy and Thermoeconomic Analysis for an Underground Train Station Air-Conditioning Cooling System. Entropy. 2016; 18(3):86. https://0-doi-org.brum.beds.ac.uk/10.3390/e18030086

Chicago/Turabian Style

Liao, Ke Yang, and Yew Khoy Chuah. 2016. "Exergy and Thermoeconomic Analysis for an Underground Train Station Air-Conditioning Cooling System" Entropy 18, no. 3: 86. https://0-doi-org.brum.beds.ac.uk/10.3390/e18030086

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop