Next Article in Journal
On the Mortality of Companies
Next Article in Special Issue
Comparison of the Parameters of the Exergoeconomic Environmental Analysis of Two Combined Cycles of Three Pressure Levels with and without Postcombustion
Previous Article in Journal
Information Theory and Its Application in Machine Condition Monitoring
Previous Article in Special Issue
Thermodynamic Performance of a Brayton Pumped Heat Energy Storage System: Influence of Internal and External Irreversibilities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy, Exergy, Exergoeconomic and Emergy-Based Exergoeconomic (Emergoeconomic) Analyses of a Biomass Combustion Waste Heat Recovery Organic Rankine Cycle

by
Saeed Khojaste Effatpanah
1,
Mohammad Hossein Ahmadi
1,
Seyed Hamid Delbari
2 and
Giulio Lorenzini
3,*
1
Faculty of Mechanical and Mechatronics Engineering, Shahrood University of Technology, Shahrood 3619995161, Iran
2
Department of Renewable Energies and Environmental, Faculty of New Sciences and Technologies, University of Tehran, Tehran 1439957131, Iran
3
Department of Engineering and Architecture, University of Parma, Parco Area Delle Scienze, 181/A, 43124 Parma, Italy
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(2), 209; https://doi.org/10.3390/e24020209
Submission received: 1 December 2021 / Revised: 12 January 2022 / Accepted: 26 January 2022 / Published: 28 January 2022

Abstract

:
In recent decades, there has been an increasing trend toward the technical development of efficient energy system assessment tools owing to the growing energy demand and subsequent greenhouse gas emissions. Accordingly, in this paper, a comprehensive emergy-based exergoeconomic (emergoeconomic) method has been developed to study the biomass combustion waste heat recovery organic Rankine cycle (BCWHR-ORC), taking into account thermodynamics, economics, and sustainability aspects. To this end, the system was formulated in Engineering Equation Solver (EES) software, and then the exergy, exergoeconomic, and emergoeconomic analyses were conducted accordingly. The exergy analysis results revealed that the evaporator unit with 55.05 kilowatts and the turbine with 89.57% had the highest exergy destruction rate and exergy efficiency, respectively. Based on the exergoeconomic analysis, the cost per exergy unit ( c ) , and the cost rate ( C ˙ ) of the output power of the system were calculated to be 24.13 USD/GJ and 14.19 USD/h, respectively. Next, by applying the emergoeconomic approach, the monetary emergy content of the system components and the flows were calculated to evaluate the system’s sustainability. Accordingly, the turbine was found to have the highest monetary emergy rate of capital investment, equal to 5.43 × 10 12   sej / h , and an output power monetary emergy of 4.77 × 10 4   sej / J . Finally, a sensitivity analysis was performed to investigate the system’s overall performance characteristics from an exergoeconomic perspective, regarding the changes in the transformation coefficients (specific monetary emergy).

1. Introduction

Overexploitation of fossil fuels to meet the ever-increasing global energy demand triggers the exhaustion of conventional resources, global warming, air pollution, and ozone depletion [1]. At present, energy conservation and air pollution reduction measures have comprised the main strategies to address the above issues, given the high base population [2]. Waste heat recovery is a promising technology specifically tailored to this end, which has been practiced favorably in recent decades [3]. In this regard, the Organic Rankine Cycle (ORC), known to researchers since the 70s, is considered a potential reliable technology to recover heat efficiently from low to mid-temperature heat sources in sustainable energy utilization systems [4,5], as well as in conventional industrial processes [6,7].
The four main components of the basic ORC are the pump, evaporator, turbine, and condenser [8]. It offers several advantages, including lower operating temperatures and pressure ranges, lower O&M (Operation and Maintenance) costs, safe and autonomous operation, and simplicity, which cultivates interest in the technology [9,10]. Unlike the conventional Rankine cycle, ORC is widely used at micro (below 15 kW) and small (below 100 kW) scales [5,11], as well as in decentralized CHP (Combined Heat and Power) units [12,13]. Typically, the working fluid of an ORC is fluids with low boiling points (for instance, alcohols, ether, refrigerant, etc.), while conventional Rankine cycles operate with steam or water [14]. Utilizing a dry working fluid eliminates the need for superheating and reduces mechanical stresses, thereby lowering the O&M costs and extending the system’s lifetime [15].
ORC power plant coupling with waste heat from industrial processes, biomass incinerators, or geothermal energy has been successfully commercialized. According to the latest report provided in ref. [16], the total capacity of ORC power plants installed at 1754 sites exceeds 2.7 GW, with a biomass share of nearly 11% (301 MW) from 332 sites. It should be noted that most of the biomass combustion waste heat recovery ORC (BCWHR-ORC) units are installed at multi-purpose manufacturing facilities [17].
Numerous low to mid-temperature energy sources such as solar energy [18,19], geothermal energy [20,21], biomass [22,23], and industrial waste heat [24,25] generate clean electricity that reduces environmental pollution while improving energy efficiency. The selection of working fluid [26], operational conditions [27], and integration of ORC with different processes [28] significantly affect the net power output of the system. Additionally, employing optimization methods provides a satisfactory trade-off between the thermal efficiency, rate of the recovered waste heat, and economic factors [29].
A variety of research has been carried out on ORC performance improvement so that they can be widely utilized in waste heat recovery applications over various temperature ranges. In this light, Feng et al. [30] introduced an irreversible regenerative ORC, and optimized the output power and efficiency while keeping the surface area available for heat transfer constant in all heat exchangers. In another study, Braimakis et al. [31] optimized three improved designs of regenerative ORC using various working fluids from an energy analysis perspective. The coupling of a vapor compression cycle with an ORC for heat recovery purposes has been investigated by Zhar et al. [32]; through a multi-objective optimization of their parametric analysis, the proposed system’s return of investment was calculated to be 6.3 years. The implementation of ORC to recover the waste heat capacity in the aluminum industry has been studied by Dokl et al. [33]. The authors concluded that up to 830 kW of electrical power could be generated using the waste heat produced by the Slovenian aluminum manufacturer they investigated. In an experimental study, Wang et al. [34] achieved a turbine isentropic efficiency of 88.6% in a 300 kW low-temperature waste heat recovery ORC. This efficiency corresponded to a waste heat source temperature of 121 °C. Ming et al. [35] extracted electrical power from the waste heat produced from an aluminum melting process using an ORC. Having assessed the thermodynamic performance of the proposed design, they suggested utilizing such systems in aluminum electrolysis facilities as an effective energy conservation measure.
Recently, few researchers have investigated the applicability of mid to high-temperature biomass combustion waste heat recovery with ORC. For instance, Georgousopoulos et al. [36] carried out a thermodynamic and techno-economic assessment of an ORC coupled with an integrated gasification combined cycle under three different scenarios to optimize waste heat recovery. They reported using zeotropic mixtures, as the working fluid yields the highest overall performance in all three scenarios. The Levelized Cost of Electricity of the proposed systems ranged between 35.42–35.67 EUR/MWh.
In another study, Zhang et al. [37] surveyed the thermodynamic performance of ORC coupled with waste heat recovery from the Rankine cycle, the Brayton cycle, and the thermoelectric generator. R123, R245, and R600 organic fluids in five different operational conditions were examined. The results showed that the DORC is the most efficient configuration, and when employed with the R123, the net output, thermal, and exergy efficiencies were 32.63 kW, 26.55%, and 54.36%, respectively. Thermodynamic and exergoeconomic analyses of ORC in combination with biomass-integrated co-firing and biomass-integrated post-firing technologies—using a mixture of natural gas and biomass as a feedstock, and externally fired technology—which solely runs on biomass, have been conducted by Mahramian et al. [38]. Having examined various fluids such as R141b, R123, n-Pentane, HFE7000, and water, they concluded that R141b yields the highest thermal and exergy efficiencies; however, it is the least favorable choice financially.
Wang et al. [39] proposed a novel organic Rankine cycle-based micro-scale cogeneration system operating with a two-stage pressure evaporator. Given the operational characteristics of the system, they examined R141b and R123 and deduced that the former results in lower capital costs and higher performance. Applying a versatile optimization algorithm, they found that the system is capable of producing 1.66 kWe and 37.16 kWth at the thermal efficiencies of 77.8% and 11.28% for CHP and ORC, respectively. The capital rate is estimated to be 0.363 USD/h. Oyekale et al. [40] studied the techno-economic aspects of biomass retrofitting in hybrid concentrated solar power biomass ORC (CSP-Biomass ORC) power plants in both constant and modular states. They used an operating CSP-ORP facility located at Ottana, Itlay, which includes linear Fresnel collectors coupled with two thermal oil storage tanks and a 630 kW ORC unit, as their case study. The results demonstrated that biomass retrofitting at best leads to 5% electrical efficiency improvement and up to 3500 h increase in operational hours. The payback period of the proposed method is 1.4 years with LCOE and NPV of 109 GBP/MWh and GBP 1.83 million, respectively.
Reviewing the literature around ORCs demonstrates that numerous researchers have carried out the optimization and performance improvement of ORCs in various applications. In particular, utilizing ORCs to recover waste heat produced in various processes has been extensively studied from thermodynamic, economic, and environmental viewpoints. However, to the best of the authors’ knowledge, an emergy-based exergoeconomic (emergoeconomic) study of biomass combustion heat recovery combined with an ORC has not yet been carried out. Hence, to fill this knowledge gap, a comprehensive evaluation of a BCWHR-ORC from the energy, exergy, exergoeconomic, and emergoeconomic perspectives have been conducted in this study. The emergoeconomic approach, first introduced by Aghbashlo and Rosen [41], is one of the most robust methods used for analyzing energy systems in recent years. The results not only contain comprehensive information on thermoeconomic characteristics but also provide an appropriate understanding of a system’s stability. The mentioned approach is adopted in this study through the eight steps shown in the Figure 1 flowchart. In the next sections, the following is discussed: Section 2 provides a descriptive account of the proposed system. In Section 3, governing equations and the math behind the modeling are discussed thoroughly, while Section 4 presents the modeling outputs and the results of a sensitivity analysis to account for parameter uncertainties. Finally, the paper is concluded in Section 5.

2. System Description

Figure 2 shows the BCWHR-ORC system studied in this paper. Evaporation takes place within the evaporator, in which the heat is transferred from thermal oil to the organic working fluid at a constant pressure. The thermal oil heats the working fluid from the subcooled liquid phase (point 4) at the pump’s output to the superheated steam required for the turbine operation (point 1). Figure 3 depicts the whole process in the T-S diagram. As the organic fluid passes through the turbine, it expands, generates power by rotating the turbine’s shaft, and builds its pressure up to the working pressure of the condenser (point 2); (note that point 2s in the diagram represents the isentropic expansion process within the turbine). Then, the working fluid is cooled down at the condenser’s pressure by expelling its heat into the cooling water cycle (7→8). Finally, the saturated liquid (point 3) is compressed up to the evaporator’s pressure (point 4), and the process repeats. The thermal oil gains its heat content from the flue gas of a biomass incinerator, enters the evaporator at a high temperature (point 5), and exits at a lower temperature (point 6) after losing its heat to the ORC fluid.

Selection of the Working and Heat Transfer Fluids

Generally, flue gas heat can be recovered either directly or indirectly. In the direct method, the heat is exchanged between the organic working fluid and flue gas through a heat exchanger. However, the indirect method uses a thermal oil loop as a heat transfer medium. The latter has been studied in this paper to avoid working fluid disintegration due to the high flue gas temperature. To select a suitable thermal oil, characteristics such as having high thermal stability, conductivity, and heat capacity should be considered. On the other hand, a low expansion coefficient and viscosity—to minimize the volumetric change and the compression work—low inflammability, and being non-toxic are the other determinants for choosing a thermal oil. Consequently, Therminol VP-1 was used in this system as the thermal oil, the complete specifications of which can be found in reference [42].
Organic fluids with a high critical temperature provide better thermodynamic performance. Lai et al. [43] have suggested linear siloxanes, alkanes, and aromatic compounds suitable for high-temperature ORC systems. In this study, m-xylene, which was previously studied in [44], and water are chosen as the working fluid and the coolant, respectively.

3. Mathematical Modeling

The BCWHR-ORC system mentioned above has been modeled using EES software. The following assumptions were adapted to simplify the developed code:
  • Steady-state operation.
  • Negligible pressure loss in the condenser, the evaporator, and the piping.
  • Negligible heat loss through equipment.
  • Atmospheric pressure and room temperature (298 K) are assumed as the dead state for exergy calculation.
The static parameters introduced in the model are also shown in Table 1.

3.1. First Law of Thermodynamics (Energy Concept)

The mass and energy conservation equations (Equations (1) and (2)), excluding potential and kinetic energy terms, were used to study the system equipment.
m ˙ i = m ˙ e
Q ˙ + m ˙ i h i = W ˙ + m ˙ e h e
Table 2 represents the previous equations, rewritten for each piece of equipment. The thermal efficiency equation of the whole system is given below:
η t h = W ˙ n e t Q ˙ i n
W ˙ n e t = W ˙ T W ˙ P
The total heat transfer coefficient (UA) of the systems’ heat exchangers (condenser and evaporator) is calculated using the following equations:
U A t o t a l = Q ˙ E v a p Δ T L M T D , E v a p + Q ˙ C o n d Δ T L M T D , C o n d
Δ T L M T D , E v a p = ( T 5 T 1 ) ( T 6 T 4 ) l n ( T 5 T 1 ) ( T 6 T 4 )
Δ T L M T D , C o n d = ( T 2 T 8 ) ( T 3 T 7 ) l n ( T 2 T 8 ) ( T 3 T 7 )

3.2. Second Law of Thermodynamics (Exergy Concept)

To assess the performance of a system, particularly from a financial perspective, exergy analysis is a viable method for determining exergy loss and its destruction rates in the system equipment to identify potential performance improvement possibilities. Given that no chemical reaction occurs within the process and that there are negligible contributions from kinetics and potential energies, the exergy flow rate at a certain point in the system is expressed as follows:
E x ˙ t o t a l , i = E x ˙ p h , i = m ˙ [ ( h i h o ) T 0 ( s i s o ) ]
where h o and s o are the enthalpy and entropy values at the dead state, respectively. Ignoring heat loss to the surrounding environment, the rate of exergy destruction for each piece of equipment is calculated using Equation (9):
E x ˙ D , k = E x ˙ F , k E x ˙ P , k
Additionally, the ratio of exergy destruction, which calculates the share of each piece of equipment in the total exergy destruction, is given below:
y D , k = E x ˙ D , k E x ˙ D , t o t a l
where
E x ˙ D , t o t a l = E x ˙ D , E v a p + E x ˙ D , T u r b + E x ˙ D , C o n d + E x ˙ D , P u m p
Finally, the efficiency of the second law (exergy efficiency) for each component and the whole system is calculated using Equations (12) and (13), respectively:
ε k = E x ˙ P , k E x ˙ F , k × 100 = ( 1 E x ˙ D , k E x ˙ F , k ) × 100
ε t o t a l = E x ˙ P , t o t a l E x ˙ i n p u t × 100 = ( W ˙ n e t E x ˙ 5 E x ˙ 4 ) × 100 = ( W ˙ T W ˙ P m ˙ o i l × ( e x 5 e x 4 ) ) × 100
The exergy efficiency and the destruction rate equation of the system components are included in Table 3.

3.3. Exergoeconomic Analysis

Exergoeconomics is an engineering tool that combines exergy and financial analyses to achieve cost-effective system design. Through an exergoeconomic analysis, the efficacy of energy conversion systems is compared, concerning the unit cost of products. Such a comparison is not feasible when conducting separate exergy and financial analyses. Various methods, including the cost of energy theory [45], the average cost approach [46], and the specific exergy cost (SPECO) theory [47], have been proposed for exergoeconomic system analyses. The latter, which consists of the following three steps, has been applied in this study:
  • Determining energy and exergy flows at the component’s boundaries.
  • Determining the fuel and product of each component.
  • Developing cost balance and auxiliary equations for the system components using Equation (14):
    e C ˙ e , k + C ˙ W , k = i C ˙ i , k + C ˙ Q , k + Z ˙ k
    where C ˙ j = c j E j ˙ , so Equation (11) is rewritten as Equation (15):
    e ( c e E x ˙ e ) k + ( c w W ˙ ) k = i ( c i E x ˙ i ) k + ( c q E x ˙ q ) k + Z ˙ k
    where   c i ,  c e ,  c q and c w represent the cost per unit of exergy in USD/GJ.   Z ˙ k   is the total investment cost of the kth component, including the capital cost rate ( Z ˙ k C I ) and the O&M cost ( Z ˙ k O M ). To calculate Z ˙ k , certain parameters such as the capital return factor (CRF), O&M factor ( φ r ), annual operating hours (OH), and capital cost of the kth equipment ( P E C k ) , are to be defined. Finally, Z ˙ k is calculated for all components using Equations (16) and (17) [48]:
    Z ˙ k = Z ˙ k C I + Z ˙ k O M = ( P E C k × C R F × φ r ) O H × 3600
    C R F = i ( i + 1 ) n ( i + 1 ) n 1
    where i and n are the discount ratio and the system’s lifetime, respectively. In this study, the bare module cost method was adopted to calculate the equipment purchase cost of the system [49,50]. Equations (18)–(25) are represent the capital cost calculation of the heat exchangers (evaporator and condenser), the turbine, and the pump, respectively:
    P E C H E = C 0 , H E × [ B 1 , H E + ( B 2 , H E × F M , H E × F P , H E ) ]
    l o g C 0 , H E = [ K 1 , H E + K 2 , H E ( l o g A H E ) + K 3 , H E ( l o g A H E ) 2 ]
    l o g F P , H E = [ C 1 , H E + C 2 , H E ( l o g P H E ) + C 3 , H E ( l o g P H E ) 2 ]
    P E C T u r b i n e = C 0 , T × F M , T
    l o g C 0 , T = [ K 1 , T + K 2 , T ( l o g W ˙ T ) + K 3 , T ( l o g W ˙ T ) 2 ]
    P E C P u m p = C 0 , P × [ B 1 , P + ( B 2 , P × F M , P × F P , P ) ]
    l o g C 0 , P = [ K 1 , P + K 2 , P ( l o g W ˙ P ) + K 3 , P ( l o g W ˙ P ) 2 ]
    l o g F P , P = [ C 1 , P + C 2 , P ( l o g 10 P P ) + C 3 , P ( l o g 10 P P ) 2 ]
In the above equations,   C 0 , F M , F P , A , and P are the initial cost, material factor, pressure factor, surface area, and equipment pressure, respectively. Additionally,   B 1 , B 2 , C 1 , C 2 , C 3 , K 1 , K 2 , and K 3 are material constants—their values are given in Table 4. In this study, carbon steel is assumed to be the equipment fabrication material due to its high durability against water and the organic working fluid. Note that the calculated equipment initial costs via Equations (18)–(25) are based on 2001 data as the reference year. Using Equation (26) and the chemical engineering plant cost index (CEPCI), the costs can be projected to its equivalent values in 2019:
P E C k , 2019 = C E P C I 2019 C E P C I 2001 × P E C k , 2001
After calculating the components’ initial costs using their relevant auxiliary equations derived from the F and P principles of the SPECO method [47], the cost balance equations of all components are solved within the EES simultaneously. The cost balance and auxiliary equations for the system components are provided in Table 5. Additionally, the static parameters employed for the exergoeconomic analysis are listed in Table 6.
To better understand the exergoeconomic performance of individual components, it is necessary to define the key parameters of the analysis. These include average cost per unit exergy of fuel ( c F , k ) and product ( c P , k ), cost flow rate of exergy destruction ( C ˙ D , k ), total cost rate ( C ˙ T o t , k ), relative cost difference ( r k ), and the exergoeconomic factor ( f k ), which are defined respectively in Equations (27)–(31) [48] (losses have not been included):
c F , k = C ˙ F , k E x ˙ F , k   ,     c P , k = C ˙ P , k E x ˙ P , k
C ˙ D , k = C ˙ F , k E x ˙ D , k
C ˙ T o t , k = Z ˙ k + C ˙ D , k
r k = c P , k c F , k c F , k
f k = Z ˙ k Z ˙ k + C ˙ D , k
In Equation (31), r k represents the difference between the specific cost of fuel and the product for a given component due to the existing cost rates of exergy destruction C ˙ D , k and capital investment Z ˙ k . The exergoeconomic factor f k   in Equation (28) quantifies the relative importance of a component’s cost to its exergy destruction cost [48].

3.4. Emergy Concept

Emergy refers to the amount of available energy required, either directly or indirectly, to generate a given output flow or storage of energy or service, expressed as solar equivalent joules (sej). The concept was developed in the 80s by Odum et al. at the University of Florida, based on the thermodynamic principle and general systems theory to evaluate their long-term sustainability [54]. Emergy is a suitable measure for calculating the consumption of natural resources as a function of equivalent solar energy [55]. In this method, a transformity coefficient or Unit Emergy Value (UEV) is defined to convert all process inputs such as energy, materials, workforce, and financial services to a type of available energy, for instance, solar equivalent joules (sej). The transformity values are emergy-based energy; therefore, to conduct an emergoeconomic analysis, these values should be multiplied by the scale factor ( β = 0.93 ) presented in [56] to result in the respective emergy-based exergy values. This factor is the exergetic equivalent of solar energy and represents the maximum extractable work from a solar energy flow, calculated through Equation (32):
β = 1 + 1 3 ( T 0 T S ) 4 4 3 ( T 0 T S )
where T 0   and T S   are the base temperature and solar temperature, respectively.
In this study, the flue gas of the biomass combustor is the driving force of the system, so the determination of the amount of emergy it enters into the system is required. Due to the lack of data on the transformity coefficient of biomass combustion waste heat, the reported Equation (33) for the transformity ratio of coal-fired waste heat by Zhang et al. [57] is employed in this study:
T r a n s f o r m i t y W H = T r b i o m a s s × Δ T η b i o m a s s × T
In Equation (33), T r b i o m a s s is the transformity of the output power of a conventional biomass power plant, expressed as (sej/J), Δ T is the temperature difference of the working fluid in the evaporator and the condenser in terms of K, η b i o m a s s is the ideal (Carnot) efficiency of the power plant, and T is the working fluid temperature in the evaporator (K). The T r b i o m a s s and η b i o m a s s values are determined from [58,59].

3.5. The Emergy-Based Exergoeconomic Analysis (Emergoeconomic)

The emergoeconomic analysis method introduced by Aghbashlo and Rosen [41] is a practical and powerful tool for measuring the long-term sustainability of energy systems. This method has similar rules to exergoeconomic analyses. The SPECO method is calculates the monetary emergy values of all exergy flows. To this end, first, a monetary emergy rate ( M ˙ i ) was defined by the following:
M ˙ i = m i E x ˙ i
where m i denotes the monetary emergy per exergy unit (specific monetary emergy) in terms of sej / exergy   unit , and M ˙ i   is the rate of monetary emergy in terms of sej / time   unit . An emergoeconomic balance equation for a given component is written as:
M ˙ P , k = M ˙ F , k + U ˙ k
where M ˙ P , k and M ˙ F , k , respectively, denote the rate of monetary emergy for product and fuel flows. The emergoeconomic equations of the system components are provided in Table 7. Through Equations (36) and (37), the specific monetary emergy of both the product and the fuel is calculated:
m F , k = M ˙ F , k E x ˙ F , k
m P , k = M ˙ P , k E ˙ x P , k
Accordingly, the monetary emergy rate of exergy destruction is obtained from Equation (38):
M ˙ D , k = m F , k E ˙ D , k
U ˙ k is the investment-related monetary emergy rate for each component given as follows:
U ˙ k = U ˙ k C I + U ˙ k O M = ( E M k × C R F × φ r ) N × 3600
In the above equation, U ˙ k C I and U ˙ k O M denote the rate of monetary emergy of capital investment and the O&M, respectively. E M k represents the emergy content of the k t h component and is obtained using Equations (40)–(42):
E M k = E M c o n s t r u c t , k + E M P u r c h a s e , k
E M c o n s t r u c t , k = M a s s k × E M m a t e r i a l
E M P u r c h a s e , k = P E C k × E M d o l l a r
where E M c o n s t r u c t , k , and E M P u r c h a s e , k , are the emergy content relevant to the construction and purchasing of the k t h component in sej, M a s s k is its respective mass in grams, E M m a t e r i a l is the emergy of the construction material (sej/gr), and E M d o l l a r is the monetary emergy in sej/USD. Note that all emergy values are based on energy, so they have to be multiplied by β = 0.93 to obtain their corresponding exergy-based emergy values.
The total emergy rate of a given component is obtained through the summation of the monetary emergy rates related to both the investment ( U ˙ k )   and the exergy destruction ( M ˙ D , k ) using Equation (43):
M ˙ T o t , k = U ˙ k + M ˙ D , k
The relative monetary emergy between corresponding average values of the product and the fuel is calculated using the following equation:
r m , k = m P , k m F , k m F , k × 100
Additionally, Equation (45) is used to calculate the emergoeconomic factor of each system component:
F m , k = U ˙ k U ˙ k + M ˙ D , k × 100 = U ˙ k M ˙ T o t , k × 100

4. Results and Discussion

4.1. Verification

To validate the developed code, the results derived from modeling the BCWHR-ORC system in this study were compared with those presented in [60], which closely resembles our case study. Two performance indicators, namely, the maximum power output and the maximum thermal efficiency, were compared at five distinct waste heat source temperatures ranging from 325 °C to 365 °C, over three different pinch point temperatures (5/5, 10/10, and 5/10) and varying evaporator pressure from 11 to 36 bar. As shown in Table 8, the results are in good agreement with [60].

4.2. The Energy Analysis Results

In this section, the mass and energy balance equations for the components were developed using the first law of thermodynamics (Table 2), and then the thermodynamic characteristics at various system points were calculated, as shown in Table 9. As mentioned previously, the optimal BCWHR-ORC studied in [46] was chosen as the case study. The system’s net power and the total heat transfer surface area were 160 kW and 16.18 kW/K, respectively.

4.3. The Exergy Analysis Results

To investigate the exergy performance of an energy system, all exergy flows in the system’s components should be determined. Subsequently, the rate of exergy destruction and the exergy efficiency for each system component were obtained by applying the second law of thermodynamics (Table 3). Figure 4 depicts the exergy destruction rate for the equipment. As can be seen, the highest values were attributed to the evaporator, condenser, turbine, and pump, respectively. The evaporator accounted for about 50% of the total exergy destruction rate, owing to the heat transfer at a large temperature difference between the thermal oil and the working fluid. Likewise, the exergy destruction of the condenser was 31.95% of the total value, which is due to the high temperature of the cooling water resulting in a significant exergy flow out of the system. To remedy this and to improve the overall system performance, the cooling water leaving the condenser should be utilized as a high-temperature heat source for a downstream process, as proposed by Aziz et al. [60]. The exergy efficiency of the system equipment is shown in Figure 5, in which the turbine, evaporator, condenser, and pump, at 89.57%, 84.75%, 72.01%, and 70.88%, have the highest efficiency, respectively.

4.4. The Exergoeconomic Analysis Results

To perform the exergoeconomic analysis, first the capital investments required for the equipment purchase were calculated using Equations (18)–(25) and the Table 4 data. Then, to determine cost flows in and out of the components, the corresponding cost balance equations were written, based on Table 5. These equations, along with auxiliary equations compliant with F and P principles and the biomass combustion equation, were solved to calculate the cost of flue gas streamed into the process. The values of the exergy flow, cost rate per exergy unit, and cost rate of the BCWHR-ORC system flow are provided in Table 10. Additionally, Table 11 summarizes the values of the component parameters obtained in the exergoeconomic analysis. Meanwhile, Figure 6 compares the rates of investment cost and the exergy destruction for each piece of equipment. According to Table 10, the cost rate of the output power and the cost rate of the output power per exergy unit were 14.19 USD/h and 24.13 USD/GJ, respectively. Additionally, the cost rate and cost rate per exergy unit calculated for the system were 52.45 USD/h, and 91.47 USD/GJ, respectively.
Based on Table 11, the turbine with 5.238 USD/h had the highest investment cost rate, followed by the evaporator and the pump. C ˙ D , representing the cost rate of exergy destruction, was the highest for both the evaporator and the condenser, with the former leading. This denotes that both components had the highest destruction rate from an exergoeconomic perspective, thus demanding further study to apply measures to reduce their respective cost rates. In contrast, the turbine and pump, with the exergy destruction cost rate of 0.9335 USD/h and 0.0855 USD/h, respectively, indicate a lower potential for improvements in their performance characteristics to reduce their exergy destruction cost rates. To ease the comparison, Figure 6 depicts the Z ˙ k   and C ˙ D , k   of each component. From an exergoeconomic point of view, the higher the total cost rate C ˙ T o t a l for a given component in a system, the more that component should be considered for the evaluation (both in terms of capital cost and the exergy destruction) of performance improvement opportunities. As a result, the turbine, with the highest C ˙ T o t a l , is of interest for performance improvement in this study. A high C ˙ T o t a l for a component could be either due to the high capital cost or the high exergy destruction cost rate; for the turbine, the former is the case. Following the turbine, the condenser and the evaporator, with 3.4862 USD/h and 3.3592 USD/h, respectively, had the highest C ˙ T o t a l in the system, indicating their suitable exergoeconomic potentials for performance improvement.
A high exergoeconomic factor (f) for a component emphasizes the significance of its investment cost rate against its exergy destruction cost rate. As shown in Table 11, the turbine, at 84.87%, had the highest f factor value. This implies that the turbine’s capital cost composed the majority of its cost rates, and from an exergoeconomic perspective, it has to be reduced for performance improvement. The pump, the condenser, and the evaporator are prioritized next for the evaluation of capital investment cost reduction. The low f factor of the evaporator indicates the prevalence of the exergy destruction cost rate in this piece of equipment. Additionally, a high exergy destruction cost rate ( C ˙ D ) in a component results in a correspondingly high relative cost difference ( r ). As it is evident from Table 10, the pump had the highest r value, indicating the significance of its cost rate of exergy destruction over its capital investment cost. The majority of exergy destruction in the pump is attributed to its low efficiency. Replacing the existing pump with an efficient one reduces both the rate of exergy destruction and the r   value. Placed the next in priority for reducing the exergy destruction rate, are the condenser, and the turbine, respectively.

4.5. The Results of the Emergoeconomic Analysis

The emergoeconomic evaluation of the BCWHR-ORC system was carried out according to the methodology described in Section 3.4. Based on Equation (33), the transformity coefficient of the biomass waste heat was 2.77 × 10 4   sej / J . Having determined this coefficient, the emergoeconomic cost balance and the auxiliary equations for the system’s components (Table 7) were solved simultaneously to calculate the values of the monetary emergy per exergy unit and the monetary emergy rate of the equipment. According to Table 12, the monetary emergy per exergy unit of the turbine power was 4.77 × 10 13 sej / GJ . Since the turbine power is the final output of the system, it is expected to have the highest value among the other system flows.
Table 13 shows the monetary energy values of all the system components. To obtain these values, their respective investment costs should be multiplied by 9.95 × 10 11   sej / USD , which is the ratio of the emergy to the dollar [41]. Consequently, with the highest investment cost, the turbine also had the highest monetary emergy value related to equipment purchasing, equal to 3.12 × 10 17   sej .
Additionally, similar to the monetary emergy of equipment’s capital investment, the monetary emergy related to their construction also needs to be calculated. Table 14 provides the energy and exergy-based monetary emergy values of the required materials for the power plant construction. The turbine was the heaviest component, hence requiring more steel for its fabrication. Following the turbine were the condenser, the evaporator, and the pump, also demanding large values of monetary emergy for construction.
Finally, Table 15 summarizes the output of the emergoeconomic analysis of the system, presented graphically in Figure 7. As expected, owing to the high monetary emergy rate of capital cost ( U ˙ k ) and the turbine’s weight, it had the highest monetary emergy rate of investment, with a value of 5.43 × 10 12   sej / h . Ranking from the highest to the lowest total monetary emergy rate ( M ˙ t o t a l )   were the turbine, the evaporator, the condenser, and the pump, with the respective values of 7.79 × 10 12   sej / h , 6.84 × 10 12   sej / h , 6.14 × 10 12   sej / h , and 4.85 × 10 11   sej / h . The high M ˙ T o t a l was attributed to either U ˙ k or M ˙ D , k , implicating the ample potential for performance improvement from an emergoeconomic perspective. In the present system, more than 53% of M ˙ T o t a l was due to the exergy destruction occurring in the equipment. As a result, lowering the irreversibility of the equipment’s monetary emergy, particularly in the evaporator and the condenser, improves the overall performance of the system.
The emergoeconomic factor ( f m ) of equipment varied from 19.68% to 69.71%, with the latter belonging to the turbine. This denotes that U ˙ T u r b prevails over M ˙ D , T u r b , meaning that the main contribution to M ˙ D , T u r b was from the U ˙ T u r b term. Therefore, to lower the f m , T u r b , one should attempt to decrease its capital investment cost. Contrary to the turbine, the evaporator’s exergoeconomic factor was 19.68%, implicating M ˙ D , E v a p dominance over U ˙ E v a p thus requiring the evaporator to be further investigated for lowering the exergy destruction rate. Moreover, increasing the evaporator’s monetary emergy rate of investment due to the utilization of auxiliary equipment to minimize irreversibility is recommended.
The relative monetary emergy difference ( r m ) for the equipment ranged from 22.41% to 117.6%, with the latter attributed to the pump. Following the pump were the condenser, the turbine, and the evaporator, at 54.56%, 38.43%, and 22.41%, respectively. As a result, the pump had the highest potential for a decrease in both the monetary emergy of investment and the monetary emergy of exergy destruction ( M ˙ T o t a l , P u m p ). On the other hand, the evaporator had the lowest r m , indicating the difficulty of improving the equipment’s monetary emergy of the product.
In Figure 8, the Grassmann diagram illustrates all monetary emergy flows entering (left) and leaving (right) the system. The monetary emergy rate of flue gas waste heat ( M ˙ 5 ), as the driving force of the system, was the highest amongst all the other flows, with the value of 4.58 × 10 13   sej / h . The monetary emergy rate of the cooling water ( M ˙ 7 ) was obtained at zero (according to the auxiliary emergoeconomic equations), and the remaining flows were attributed to U ˙ k . Meanwhile, three primary flows left the power plant, between which the turbine’s net power had the highest monetary emergy rate, equal to 2.81 × 10 13   sej / h .

4.6. The Sensitivity Analysis

While applying the emergoeconomic analysis, it may not be possible to accurately obtain the transformity coefficient for some processes owing to their complex and time-consuming calculations. Sometimes, the transformity coefficient of a given process is estimated using a correlation, such as Equation (33), used in this paper for calculating the BCWH transformity coefficient. Consequently, it is necessary to determine the uncertainty of the calculated monetary emergy values of the system due to the estimated transformity coefficient. For this purpose, the ψ parameter is defined to quantify the overall emergoeconomic performance of the system:
ψ = M ˙ T u r b i n e + M ˙ 8 M ˙ 5 + M ˙ 7 + M ˙ P u m p + U ˙ E v a p + U ˙ T u r b + U ˙ C o n d + U ˙ P u m p × 100
Figure 9 depicts the trends of   T r B C W H , T r e m d o l l a r , and T r S t e e l while changing their values ± 50 % with 10 %   increments. As can be seen, ψ was the most sensitive and was inversely correlated to T r B C W H changes, which is the main driving force of the system. The higher the emergy content of BCWH is, the lower the ψ value will be. However, it should be noted that the changes in ψ were less pronounced compared to T r B C W H , as a −50% change in the latter resulted in only a 3% percent decrease in the former. Contrastingly, the T r e m d o l l a r , and T r S t e e l had a direct correlation with ψ . T r S t e e l had hardly any influence on ψ , as it was roughly a straight line drawn parallel to the x-axis. This was because the monetary emergy of the equipment’s construction partook of the least from the total monetary emergy content. T r e m d o l l a r changes also impacted the ψ , since the capital investment constituted the principal share of the equipment’s total emergy content. The ± 50 % change in T r e m d o l l a r , resulted in −1.83% to 1.58% changes in the ψ value.

5. Conclusions

In this paper, a BCWHR-ORC system is presented and evaluated from the thermodynamic, economic, and sustainability aspects by implementing energy, exergy, exergoeconomic, and emergoeconomic analyses. Then, a sensitivity analysis was conducted to account for the effect of the transformity coefficient uncertainties on the overall emergoconomic performance of the system. The key results are concluded as follows:
  • The highest rate of exergy destruction ( E x ˙ d ) occurred in the evaporator and the condenser, with the values of 49.91% and 31.95%, respectively. Meanwhile, the turbine, at 89.57%, had the highest exergy efficiency ( ε ) within the system.
  • According to the exergoeconomic analysis, the cost per exergy unit of the turbine’s power ( c ) was equal to 24.13 USD/GJ, and its cost rate of output power ( C ˙ ) was 14.19 USD/h. Additionally, C ˙ and c for the whole system were 52.45 USD/h and 91.47 USD/GJ, respectively.
  • Amongst the system equipment, the turbine had the highest total cost rate ( C ˙ T o t a l ) and exergoeconomic factors ( f ), with values of 6.1715 USD/h and 84.87%, respectively. It has been deduced that the cost rate of investment for the turbine ( Z ˙ T u r b ) is large, and measures should be taken to reduce its capital investment cost.
  • Conducting an emergoeconomic analysis, the monetary emergy per exergy unit ( m ) and the the monetary emergy rate of the output power ( M ˙ ) were 8.24 × 10 4   sej / J , and 4.87 × 10 13   sej / h , respectively.
  • The highest monetary emergy rate of capital investment ( U ˙ ) belonged to the turbine, the condenser, the evaporator, and the condenser, with the corresponding values of 1.77 × 10 12 , 5.43 × 10 12 , 1.35 × 10 12 , and 3.16 × 10 11 . In addition, ranked from the highest to the lowest, the monetary emergy of exergy destruction ( M ˙ D ) for the evaporator, the condenser, the turbine, and the pumps were 2.64 × 10 12 , 4.90 × 10 12 , 6.18 × 10 12 , and 1.85 × 10 11 , respectively.
  • The overall emergoeconomic factor of the BCWHR-ORC was 44.41%, implying that 55.59% of the total monetary emergy was due to exergy destruction (irreversibility) within the system. The reduction of exergy destruction results in system sustainability and performance improvement from an emergoeconomic perspective.
  • The turbine and the pump had the largest emergoeconomic factor ( f m , T u r b = 67.28 ) and relative monetary emergy difference ( r m , P u m p = 110.8 ), respectively, meaning that U ˙ T u r b was more pronounced than the M ˙ D , T u r b ; however, for the pump, it was the other way around ( M ˙ D , P u m p dominates U ˙ P u m p ).
  • Using the sensitivity analysis, it was found that the overall emergoeconomic performance of the system ( ψ ) was the most sensitive to transformity coefficients of the biomass combustion waste heat ( T r B C W H ) and the emergy per dollar ( T r e m d o l l a r ), respectively.

Author Contributions

Investigation, S.K.E., M.H.A., S.H.D. and G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tiwari, D.; Sherwani, A.F.; Atheaya, D.; Kumar, A.; Kumar, N. Thermodynamic analysis of Organic Rankine cycle driven by reversed absorber hybrid photovoltaic thermal compound parabolic concentrator system. Renew. Energy 2019, 147, 2118–2127. [Google Scholar] [CrossRef]
  2. Little, A.B.; Garimella, S. Comparative assessment of alternative cycles for waste heat recovery and upgrade. Energy 2011, 36, 4492–4504. [Google Scholar] [CrossRef]
  3. Chen, Y.; Han, W.; Jin, H. Thermodynamic performance optimization of the absorption-generation process in an absorption refrigeration cycle. Energy Convers. Manag. 2016, 126, 290–301. [Google Scholar] [CrossRef]
  4. Dong, L.; Liu, H.; Riffat, S. Development of small-scale and micro-scale biomass-fueled CHP systems—A literature review. Appl. Therm. Eng. 2009, 29, 2119–2126. [Google Scholar] [CrossRef] [Green Version]
  5. Franco, A.; Villani, M. Optimal design of binary cycle power plants for water-dominated, medium-temperature geothermal fields. Geothermics 2009, 38, 379–391. [Google Scholar] [CrossRef] [Green Version]
  6. Delgado-Torres, A.M.; García-Rodríguez, L. Preliminary design of seawater and brackish water reverse osmosis desalination systems driven by low-temperature solar organic Rankine cycles (ORC). Energy Convers. Manag. 2010, 51, 2913–2920. [Google Scholar] [CrossRef]
  7. Tian, H.; Shu, G.; Wei, H.; Liang, X.; Liu, L. Fluids and parameters optimization for the organic Rankine cycles (ORCs) used in exhaust heat recovery of Internal Combustion Engine (ICE). Energy 2012, 47, 125–136. [Google Scholar] [CrossRef]
  8. Jiménez-Arreola, M.; Wieland, C.; Romagnoli, A. Direct vs indirect evaporation in Organic Rankine Cycle (ORC) systems: A comparison of the dynamic behavior for waste heat recovery of engine exhaust. Appl. Energy 2019, 242, 439–452. [Google Scholar] [CrossRef]
  9. Zhang, C.; Lin, J.; Tan, Y. A theoretical study on a novel combined organic Rankine cycle and ejector heat pump. Energy 2019, 176, 81–90. [Google Scholar] [CrossRef]
  10. Hung, T.; Shai, T.; Wang, S. A review of organic rankine cycles (ORCs) for the recovery of low-grade waste heat. Energy 1997, 22, 661–667. [Google Scholar] [CrossRef]
  11. Preißinger, M.; Heberle, F.; Brüggemann, D. Thermodynamic analysis of double-stage biomass fired Organic Rankine Cycle for micro-cogeneration. Int. J. Energy Res. 2012, 36, 944–952. [Google Scholar] [CrossRef]
  12. Tchanche, B.F.; Papadakis, G.; Lambrinos, G.; Frangoudakis, A. Fluid selection for a low-temperature solar organic Rankine cycle. Appl. Therm. Eng. 2009, 29, 2468–2476. [Google Scholar] [CrossRef] [Green Version]
  13. Franco, A. Power production from a moderate temperature geothermal resource with regenerative Organic Rankine Cycles. Energy Sustain. Dev. 2011, 15, 411–419. [Google Scholar] [CrossRef]
  14. Narayanan, S.; Cai, C.Y.; Assaraf, Y.G.; Guo, H.Q.; Cui, Q.; Wei, L.; Huang, J.J.; Ashby, C.R., Jr.; Chen, Z.S. Targeting the ubiquitin-proteasome pathway to overcome anti-cancer drug resistance. Drug Resist. Updat. 2020, 48, 100663. [Google Scholar] [CrossRef] [PubMed]
  15. Quoilin, S.; Broek, M.V.D.; Declaye, S.; Dewallef, P.; Lemort, V. Techno-economic survey of Organic Rankine Cycle (ORC) systems. Renew. Sustain. Energy Rev. 2013, 22, 168–186. [Google Scholar] [CrossRef] [Green Version]
  16. Tartiere, T.; Astolfi, M. A World Overview of the Organic Rankine Cycle Market. Energy Procedia 2017, 129, 2–9. [Google Scholar] [CrossRef]
  17. Algieri, A.; Morrone, P. Comparative energetic analysis of high-temperature subcritical and transcritical Organic Rankine Cycle (ORC). A biomass application in the Sibari district. Appl. Therm. Eng. 2012, 36, 236–244. [Google Scholar] [CrossRef]
  18. Marinheiro, M.M.; Coraca, G.M.; Gomez, L.C.; Ribatski, G. Detailed transient assessment of a small-scale concentrated solar power plant based on the organic Rankine cycle. Appl. Therm. Eng. 2022, 204, 117959. [Google Scholar] [CrossRef]
  19. Sinasac, Z.; Jianu, O.A. Parametric study on the exergetic and cyclic performance of a solar-powered organic Rankine cycle coupled with a thermal energy storage and complete flashing cycle. Sustain. Energy Technol. Assess. 2021, 45, 101172. [Google Scholar] [CrossRef]
  20. Astolfi, M.; Romano, M.C.; Bombarda, P.; Macchi, E. Binary ORC (Organic Rankine Cycles) power plants for the exploitation of medium–low temperature geothermal sources—Part B: Techno-economic optimization. Energy 2014, 66, 435–446. [Google Scholar] [CrossRef]
  21. Shokati, N.; Ranjbar, F.; Yari, M. Exergoeconomic analysis and optimization of basic, dual-pressure and dual-fluid ORCs and Kalina geothermal power plants: A comparative study. Renew. Energy 2015, 83, 527–542. [Google Scholar] [CrossRef]
  22. Braimakis, K.; Thimo, A.; Karellas, S. Technoeconomic Analysis and Comparison of a Solar-Based Biomass ORC-VCC System and a PV Heat Pump for Domestic Trigeneration. J. Energy Eng. 2017, 143, 04016048. [Google Scholar] [CrossRef]
  23. Karellas, S.; Braimakis, K. Energy–exergy analysis and economic investigation of a cogeneration and trigeneration ORC–VCC hybrid system utilizing biomass fuel and solar power. Energy Convers. Manag. 2016, 107, 103–113. [Google Scholar] [CrossRef]
  24. Wang, H.; Xu, J.; Yang, X.; Miao, Z.; Yu, C. Organic Rankine cycle saves energy and reduces gas emissions for cement production. Energy 2015, 86, 59–73. [Google Scholar] [CrossRef]
  25. Pereira, E.; Júnior, B.; Diniz, M.; Arrieta, P.; Raul, F.; Arrieta, P. Assessment of a Kalina cycle for waste heat recovery in the cement industry. Appl. Therm. Eng. 2019, 147, 421–437. [Google Scholar]
  26. Meinel, D.; Wieland, C.; Spliethoff, H. Effect and comparison of different working fluids on a two-stage organic rankine cycle (ORC) concept. Appl. Therm. Eng. 2014, 63, 246–253. [Google Scholar] [CrossRef]
  27. He, C.; Liu, C.; Gao, H.; Xu, X.; Xie, H. Parameters optimization and performance analysis of organic rankine cycle for industrial waste heat recovery. J. Eng. Thermophys. 2012, 33, 2042–2046. [Google Scholar]
  28. Yu, H.; Feng, X.; Wang, Y.; Biegler, L.T.; Eason, J. A systematic method to customize an efficient organic Rankine cycle (ORC) to recover waste heat in refineries. Appl. Energy 2016, 179, 302–315. [Google Scholar] [CrossRef]
  29. Yu, H.; Eason, J.; Biegler, L.T.; Feng, X. Simultaneous heat integration and techno-economic optimization of Organic Rankine Cycle (ORC) for multiple waste heat stream recovery. Energy 2017, 119, 322–333. [Google Scholar] [CrossRef]
  30. Feng, H.; Chen, W.; Chen, L.; Tang, W. Power and efficiency optimizations of an irreversible regenerative organic Rankine cycle. Energy Convers. Manag. 2020, 220, 113079. [Google Scholar] [CrossRef]
  31. Braimakis, K.; Karellas, S. Energetic optimization of regenerative Organic Rankine Cycle (ORC) configurations. Energy Convers. Manag. 2018, 159, 353–370. [Google Scholar] [CrossRef]
  32. Zhar, R.; Allouhi, A.; Ghodbane, M.; Jamil, A.; Lahrech, K. Parametric analysis and multi-objective optimization of a combined Organic Rankine Cycle and Vapor Compression Cycle. Sustain. Energy Technol. Assess. 2021, 47, 101401. [Google Scholar] [CrossRef]
  33. Dokl, M.; Gomilsek, R.; Cucek, L.; Abikoye, B.; Kravanja, Z. Maximizing the power output and net present value of organic Rankine cycle: Application to aluminium industry. Energy 2022, 239, 122620. [Google Scholar] [CrossRef]
  34. Wang, R.; Kuang, G.; Zhu, L.; Wang, S.; Zhao, J. Experimental Investigation of a 300 kW Organic Rankine Cycle Unit with Radial Turbine for Low-Grade Waste Heat Recovery. Entropy 2019, 21, 619. [Google Scholar] [CrossRef] [Green Version]
  35. Ming, Y.; Zhou, N. Thermodynamic Performance Analysis of a Waste Heat Power Generation System (WHPGS) Applied to the Sidewalls of Aluminum Reduction Cells. Entropy 2020, 22, 1279. [Google Scholar] [CrossRef]
  36. Georgousopoulos, S.; Braimakis, K.; Grimekis, D.; Sotirios Karellas, S. Thermodynamic and techno-economic assessment of pure and zeotropic fluid ORCs for waste heat recovery in a biomass IGCC plant. Appl. Therm. Eng. 2021, 183, 116202. [Google Scholar] [CrossRef]
  37. Zhang, C.; Shu, G.; Tian, H.; Wei, H.; Liang, X. Comparative study of alternative ORC-based com-bined power systems to exploit high temperature waste heat. Energy Convers. Manag. 2015, 89, 541–554. [Google Scholar] [CrossRef]
  38. Moharamian, A.; Soltani, S.; Rosen, M.A.; Mahmoudi, S.M.S.; Morosuk, T. A comparative thermoeco-nomic evaluation of three biomass and biomass-natural gas fired combined cycles using organic Rankine cycles. J. Clean. Prod. 2017, 161, 524–544. [Google Scholar] [CrossRef]
  39. Wang, Q.; Wu, W.; He, Z. Thermodynamic analysis and optimization of a novel organic Rankine cycle-based micro-scale cogeneration system using biomass fuel. Energy Convers. Manag. 2019, 198, 111803. [Google Scholar] [CrossRef]
  40. Oyekale, J.; Heberle, F.; Petrollese, M.; Brüggemann, D.; Cau, G. Biomass retrofit for existing solar organic Rankine cycle power plants: Conceptual hybridization strategy and techno-economic assessment. Energy Convers. Manag. 2019, 196, 831–845. [Google Scholar] [CrossRef]
  41. Aghbashlo, M.; Rosen, M.A. Consolidating exergoeconomic and exergoenvironmental analyses using the emergy concept for better understanding energy conversion systems. J. Clean. Prod. 2017, 172, 696–708. [Google Scholar] [CrossRef]
  42. Therminol® VP-1. Available online: http://www.therminol.com/products/Therminol-VP1 (accessed on 25 January 2021).
  43. Lai, N.A.; Wendland, M.; Fischer, J. Working fluids for high-temperature organic Rankine cycles. Energy 2011, 36, 199–211. [Google Scholar] [CrossRef]
  44. Aziz, F.; Mudasar, R.; Kim, M.-H. Exergetic and heat load optimization of high temperature organic Rankine cycle. Energy Convers. Manag. 2018, 171, 48–58. [Google Scholar] [CrossRef]
  45. Lozano, M.; Valero, A. Theory of the exergetic cost. Energy 1993, 18, 939–960. [Google Scholar] [CrossRef]
  46. Lazzaretto, A.; Tsatsaronis, G. A General Process-Based Methodology for Exergy Costing. In Proceedings of the 1996 International Mechanical Engineering Congress and Exhibition, Atlanta, GA, USA, 17–22 November 1996; pp. 413–423. [Google Scholar]
  47. Lazzaretto, A.; Tsatsaronis, G. SPECO: A systematic and general methodology for calculating efficiencies and costs in thermal systems. Energy 2006, 31, 1257–1289. [Google Scholar] [CrossRef]
  48. Bejan, A.; Tsatsaronis, G.; Moran, M. Thermal Design and Optimization; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 1996. [Google Scholar]
  49. Turton, R.; Bailie, R.C.; Whiting, W.B.; Shaeiwitz, J.A. Analysis, Synthesis and Design of Chemical Processes; Pearson Education: London, UK, 2008. [Google Scholar]
  50. Karimi, S.H.; Mansouri, S. A comparative profitability study of geothermal electricity production in developed and de-veloping countries: Exergoeconomic analysis and optimization of different ORC configurations. Renew. Energy 2018, 115, 600–619. [Google Scholar] [CrossRef]
  51. Liu, Z.; He, T. Exergoeconomic analysis and optimization of a Gas Turbine-Modular Helium Reactor with new organic Rankine cycle for efficient design and operation. Energy Convers. Manag. 2020, 204, 112311. [Google Scholar] [CrossRef]
  52. Available online: https://www.chemengonline.com/2019-chemical-engineering-plant-cost-index-annual-average/ (accessed on 25 January 2021).
  53. Mosaffa, A.; Mokarram, N.H.; Farshi, L.G. Thermoeconomic analysis of a new combination of ammonia/water power generation cycle with GT-MHR cycle and LNG cryogenic exergy. Appl. Therm. Eng. 2017, 124, 1343–1353. [Google Scholar] [CrossRef]
  54. Reza, B.; Sadiq, R.; Hewage, K. Emergy-based life cycle assessment (Em-LCA) for sustainability appraisal of infrastructure systems: A case study on paved roads. Clean Technol. Environ. Policy 2013, 16, 251–266. [Google Scholar] [CrossRef]
  55. Rugani, B.; Benetto, E. Improvements to Emergy Evaluations by Using Life Cycle Assessment. Environ. Sci. Technol. 2012, 46, 4701–4712. [Google Scholar] [CrossRef] [PubMed]
  56. Bastianoni, S.; Facchini, A.; Susani, L.; Tiezzi, E. Emergy as a function of exergy. Energy 2007, 32, 1158–1162. [Google Scholar] [CrossRef]
  57. Zhang, H.; Guan, X.; Ding, Y.; Liu, C. Emergy analysis of Organic Rankine Cycle (ORC) for waste heat power generation. J. Clean. Prod. 2018, 183, 1207–1215. [Google Scholar] [CrossRef]
  58. Odum, H.T. Environmental Accounting, Emergy and Environmental Decision Making; John Wiley: Hoboken, NJ, USA, 1996; 370p. [Google Scholar]
  59. Sha, S.; Hurme, M. Emergy evaluation of combined heat and power plant processes. Appl. Therm. Eng. 2011, 43, 67–74. [Google Scholar] [CrossRef]
  60. Aziz, F.; Salim, M.S.; Kim, M.-H. Performance analysis of high temperature cascade organic Rankine cycle coupled with water heating system. Energy 2019, 170, 954–966. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the implemented emergoeconomic analysis.
Figure 1. Flowchart of the implemented emergoeconomic analysis.
Entropy 24 00209 g001
Figure 2. Schematic of the BCWHR-ORC system.
Figure 2. Schematic of the BCWHR-ORC system.
Entropy 24 00209 g002
Figure 3. T-S (Temperature-Entropy) diagram of the BCWHR-ORC system.
Figure 3. T-S (Temperature-Entropy) diagram of the BCWHR-ORC system.
Entropy 24 00209 g003
Figure 4. The rate of exergy destruction attributed to each component.
Figure 4. The rate of exergy destruction attributed to each component.
Entropy 24 00209 g004
Figure 5. The exergy efficiency of the system equipment.
Figure 5. The exergy efficiency of the system equipment.
Entropy 24 00209 g005
Figure 6. The bar diagram of the investment ( Z ˙ k ) cost rate and the exergy destruction cost rate ( C ˙ D , k ) cost rates of the system equipment.
Figure 6. The bar diagram of the investment ( Z ˙ k ) cost rate and the exergy destruction cost rate ( C ˙ D , k ) cost rates of the system equipment.
Entropy 24 00209 g006
Figure 7. The bar diagram of the monetary emergy rate of capital investment ( U ˙ k ) and the monetary emergy rate of exergy destruction ( M ˙ D , k ) for the system equipment.
Figure 7. The bar diagram of the monetary emergy rate of capital investment ( U ˙ k ) and the monetary emergy rate of exergy destruction ( M ˙ D , k ) for the system equipment.
Entropy 24 00209 g007
Figure 8. The Grassmann diagram for the total monetary emergy flows of the system.
Figure 8. The Grassmann diagram for the total monetary emergy flows of the system.
Entropy 24 00209 g008
Figure 9. The trend of overall emergoeconomic performance of the system ( ψ ), with respect to the transformity coefficients.
Figure 9. The trend of overall emergoeconomic performance of the system ( ψ ), with respect to the transformity coefficients.
Entropy 24 00209 g009
Table 1. Design parameters used in the model [44].
Table 1. Design parameters used in the model [44].
Input ParameterUnitInput Value
Thermal power inputkW1000
Thermal oil inlet temperature 370
Cooling water inlet temperature 70
Flue gases inlet temperature 450
Flue gases exit temperature 180
Condensation Temperature 85
Turbine isentropic efficiency % 85
Pump isentropic efficiency % 65
Heating process efficiency % 85
Evaporator efficiency % 96
Condenser efficiency % 98
Table 2. Mass and energy balance equations of the system equipment.
Table 2. Mass and energy balance equations of the system equipment.
ComponentBalance Equation
Evaporator Q ˙ E v a p = m ˙ w f ( h 1 h 4 ) = m ˙ o i l ( h 5 h 6 ) = m ˙ o i l C p s ( T 5 T 6 )
Turbine W ˙ T = m ˙ w f ( h 1 h 2 ) = m ˙ w f ( h 1 h 2 s ) η T
Condenser Q ˙ C o n d = m ˙ w f ( h 2 h 3 ) = m ˙ c w ( h 8 h 7 ) = m ˙ c w C p c w ( T 8 T 7 )
Pump W ˙ P = m ˙ w f ( h 4 h 3 ) = m ˙ w f ( h 4 s h 3 ) / η P
Table 3. The exergy destruction rate and efficiency equations of the system components.
Table 3. The exergy destruction rate and efficiency equations of the system components.
ComponentExergy BalanceExergy Efficiency
Evaporator E x ˙ D , E v a p = ( E x ˙ 4 E x ˙ 1 ) + ( E x ˙ 5 E x ˙ 6 ) ( E x ˙ 1 E x ˙ 4 ) / ( E x ˙ 5 E x ˙ 6 )
Turbine E x ˙ D , T u r b = ( E x ˙ 1 E x ˙ 2 ) W ˙ T W ˙ T / ( E x ˙ 1 E x ˙ 2 )
Condenser E x ˙ D , C o n d = ( E x ˙ 2 E x ˙ 3 ) + ( E x ˙ 7 E x ˙ 8 ) ( E x ˙ 8 E x ˙ 7 ) / ( E x ˙ 2 E x ˙ 3 )
Pump E x ˙ D , P u m p = ( E x ˙ 3 E x ˙ 4 ) + W ˙ P ( E x ˙ 4 E x ˙ 3 ) / W ˙ P
Table 4. Constant Values of the equipment’s capital cost equations.
Table 4. Constant Values of the equipment’s capital cost equations.
ConstantsEquipment
Heat ExchangersTurbinePump
B 1 1.6300-1.8900
B 2 1.6600-1.3500
C 1 0.0388-−0.3935
C 2 −0.1127-0.3957
C 3 0.0818-−0.0023
K 1 4.32472.24763.3892
K 2 −0.30301.49650.0536
K 3 0.1634−0.16180.1538
F M 1.00003.50001.5000
Table 5. Cost balance and auxiliary equations of the system equipment.
Table 5. Cost balance and auxiliary equations of the system equipment.
ComponentCost Balance EquationAuxiliary Equations
Evaporator c 4 E x ˙ 4 + m 5 E x ˙ 5 + Z ˙ E v a p = m 1 E x ˙ 1 + m 6 E x ˙ 6 c 5 = c 6
Turbine c 1 E x ˙ 1 + Z ˙ T u r b = c 2 E x ˙ 2 + c W T E x ˙ W T c 1 = c 2
Condenser c 2 E ˙ x 2 + c 7 E x ˙ 7 + Z ˙ C o n d = c 3 E x ˙ 3 + c 8 E x ˙ 8 c 2 = c 3   ;   c 7 = 0
Pump c 3 E x ˙ 3 + c W P E x ˙ W P + Z ˙ P u m p = c 4 E ˙ x 4 c W P = c W T
Table 6. Input parameters of the exergoeconomic analysis.
Table 6. Input parameters of the exergoeconomic analysis.
ParameterUnitValueRef.
Interest rate (i)%10[51]
Plant operational hours (OH)Hours/year7446-
Plant total life time (N)years20[51]
Maintenance factor ( φ ) %6[51]
Chemical engineering plant cost index 2001 (CEPCL2001)-397[52]
Chemical engineering plant cost index 2019 (CEPCL2019)-607.5[52]
Overall heat transfer coefficient of the evaporator (U_Evap) kW / m 2 °C0.6[53]
Overall heat transfer coefficient of the condenser (U_Cond) kW / m 2 °C0.5[53]
Table 7. The auxiliary and emergy balance equations of the system components.
Table 7. The auxiliary and emergy balance equations of the system components.
ComponentEmergy-Based Cost Balance EquationsAuxiliary Equation
Evaporator m 4 E ˙ 4 + m 5 E ˙ 5 + U ˙ E v a p = m 1 E ˙ 1 + m 6 E ˙ 6 m 5 = m 6
Turbine m 1 E ˙ 1 + U ˙ T u r b = m 2 E ˙ 2 + m W T E ˙ W T m 1 = m 2
Condenser m 2 E ˙ 2 + m 7 E ˙ 7 + U ˙ C o n d = m 3 E ˙ 3 + m 8 E ˙ 8 m 7 = 0 ,     m 2 = m 3
Pump m 3 E ˙ 3 + m W P E ˙ W P + U ˙ P u m p = m 4 E ˙ 4 m W P = m W T
Table 8. Comparison of energy analysis results obtained from modeling in the present study (yellow highlighted values) with the values reported in reference [60].
Table 8. Comparison of energy analysis results obtained from modeling in the present study (yellow highlighted values) with the values reported in reference [60].
Heat Source Temperature [°C]ΔTpp = 5/5ΔTpp = 10/10ΔTpp = 5/10
Maximum Work Output Wmax [kW]Maximum Cycle Efficiency ηmax [%]Maximum Work Output Wmax [kW]Maximum Cycle Efficiency ηmax [%]Maximum Work Output Wmax [kW]Maximum Cycle Efficiency ηmax [%]
325127.7 at 11 bar15.6 at 11 bar140.6 at 17 bar17.2 at 17 bar146.0 at 21 bar17.8 at 21 bar
128.7 at 11 bar15.7 at 11 bar141.3 at 17 bar17.3 at 17 bar146.5 at 21 bar17.9 at 21 bar
335127.7 at 11 bar15.6 at 11 bar144.8 at 20 bar17.7 at 20 bar150.0 at 25 bar18.3 at 25 bar
128.7 at 11 bar15.7 at 11 bar145.4 at 20 bar17.8 at 20 bar150.4 at 25 bar18.4 at 25 bar
345127.7 at 11 bar15.6 at 11 bar150.7 at 26 bar18.5 at 26 bar152.2 at 28 bar18.6 at 28 bar
128.7 at 11 bar15.7 at 11 bar151.2 at 26 bar18.5 at 26 bar152.6 at 28 bar18.7 at 28 bar
355------153.5 at 30 bar18.8 at 30 bar154.5 at 32 bar18.9 at 32 bar
------153.8 at 30 bar18.8 at 30 bar154.9 at 32 bar19.0 at 32 bar
365------155.6 at 34 bar19.1 at 34 bar156.4 at 36 bar19.2 at 36 bar
------155.5 at 34 bar19.0 at 34 bar156.5 at 36 bar19.2 at 36 bar
Table 9. Flow type and thermodynamic characteristics at different system points.
Table 9. Flow type and thermodynamic characteristics at different system points.
State No.FluidTemperature (K)Pressure (Bar)Mass Flowrate (kg/s)Specific Enthalpy (kj/kg)Specific Entropy (kj/kg K)
1m-xylene553.6151.199777.81.714
2m-xylene458.80.17641.199641.51.767
3m-xylene358.20.17641.199108.50.3313
4m-xylene359151.199111.30.3341
5Therminol VP-1643.27.3321.830704.31.522
6Therminol VP-1448.67.3321.830267.60.7188
7Water343.2117.63188.20.588
8Water351.8117.63224.50.6923
Table 10. The exergy rate, cost rate per exergy unit, and the cost rate of system flows.
Table 10. The exergy rate, cost rate per exergy unit, and the cost rate of system flows.
State No.Exergy (kW)c (USD/GJ) C ˙ (USD/h)
1320.013.6415.707
2137.613.646.7536
311.6513.640.5717
414.0523.351.1808
5458.510.1416.736
697.5010.143.5590
7227.800
8318.56.9237.9380
Power to Pump3.38424.136.7536
Turbine Power163.424.1315.707
Table 11. The exergoeconomic outputs of the case study.
Table 11. The exergoeconomic outputs of the case study.
Components Z ˙ k (USD/h) C ˙ D , k (USD/h) C ˙ T o t , k (USD/h) c F , k (USD/GJ) c P , k (USD/GJ) r k (%) f k (%)
Evaporator1.34932.00993.359210.1413.1930.0840.17
Turbine5.23800.93356.171513.6424.1376.9684.87
Condenser1.75571.73053.486213.6424.3178.3150.36
Pump0.31500.08550.400524.1370.49192.178.65
Table 12. The values of the monetary emergy per exergy unit, and the monetary emergy rate of the system flows.
Table 12. The values of the monetary emergy per exergy unit, and the monetary emergy rate of the system flows.
State No.m (sej/GJ) M ˙ (sej/GJ)
1 3.45 × 10 13 3.97 × 10 13
2 3.45 × 10 13 1.71 × 10 13
3 3.45 × 10 13 1.45 × 10 12
4 4.63 × 10 13 2.34 × 10 12
5 2.77 × 10 13 4.58 × 10 13
6 2.77 × 10 13 9.73 × 10 12
7 0 0
8 1.52 × 10 13 1.74 × 10 13
Power to Pump 4.77 × 10 13 5.81 × 10 11
Turbine Power 4.77 × 10 13 2.81 × 10 13
Table 13. The values of the monetary emergy of the equipment’s capital investment.
Table 13. The values of the monetary emergy of the equipment’s capital investment.
ComponentCapital Investment Cost (USD)Energy-Based Transformity (sej/USD)Exergy-Based Transformity (sej/USD)Energy-Based Emergy (sej)Exergy-Based Emergy (sej)
Evaporator80,699 1.07 × 10 12 9.95 × 10 11 8.64 × 10 16 8.03 × 10 16
Turbine313,289 1.07 × 10 12 9.95 × 10 11 3.35 × 10 17 3.12 × 10 17
Condenser104,994 1.07 × 10 12 9.95 × 10 11 1.12 × 10 17 1.05 × 10 17
Pump18,839 1.07 × 10 12 9.95 × 10 11 2.02 × 10 16 1.87 × 10 16
Table 14. The values of the monetary emergy of the equipment’s construction.
Table 14. The values of the monetary emergy of the equipment’s construction.
ComponentConstruction MaterialValueUnitEnergy-Based Transformity (sej/g)Exergy-Based Transformity (sej/g)Energy-Based Emergy (sej)Exergy-Based Emergy (sej)
EvaporatorSteel 7.43 × 10 4  ag 2.77 × 10 9  a 2.58 × 10 9 2.06 × 10 14 1.92 × 10 14
TurbineSteel 5.10 × 10 6  ag 2.77 × 10 9  a 2.58 × 10 9 1.41 × 10 16 1.31 × 10 16
CondenserSteel 4.25 × 10 5  ag 2.77 × 10 9  a 2.58 × 10 9 1.18 × 10 15 1.09 × 10 15
PumpSteel 4.74 × 10 4  ag 2.77 × 10 9  a 2.58 × 10 9 1.31 × 10 14 1.22 × 10 14
a Obtained from [57].
Table 15. The output of the emergoeconomic analysis of the case study.
Table 15. The output of the emergoeconomic analysis of the case study.
Components U ˙ k (sej/h) M ˙ D , k (sej/h) M ˙ T o t a l , k (sej/h) m F , k (sej/GJ) m P , k (sej/GJ) r m , k (%) f m , k (%)
Evaporator 1.35 × 10 12 5.49 × 10 12 6.84 × 10 12 2.77 × 10 13 3.39 × 10 13 22.4119.68
Turbine 5.43 × 10 12 2.36 × 10 12 7.79 × 10 12 3.45 × 10 13 4.77 × 10 13 38.4369.71
Condenser 1.77 × 10 12 4.37 × 10 12 6.14 × 10 12 3.45 × 10 13 5.33 × 10 13 54.5628.75
Pump 3.16 × 10 11 1.69 × 10 11 4.85 × 10 11 4.77 × 10 13 1.04 × 10 14 117.665.10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Effatpanah, S.K.; Ahmadi, M.H.; Delbari, S.H.; Lorenzini, G. Energy, Exergy, Exergoeconomic and Emergy-Based Exergoeconomic (Emergoeconomic) Analyses of a Biomass Combustion Waste Heat Recovery Organic Rankine Cycle. Entropy 2022, 24, 209. https://doi.org/10.3390/e24020209

AMA Style

Effatpanah SK, Ahmadi MH, Delbari SH, Lorenzini G. Energy, Exergy, Exergoeconomic and Emergy-Based Exergoeconomic (Emergoeconomic) Analyses of a Biomass Combustion Waste Heat Recovery Organic Rankine Cycle. Entropy. 2022; 24(2):209. https://doi.org/10.3390/e24020209

Chicago/Turabian Style

Effatpanah, Saeed Khojaste, Mohammad Hossein Ahmadi, Seyed Hamid Delbari, and Giulio Lorenzini. 2022. "Energy, Exergy, Exergoeconomic and Emergy-Based Exergoeconomic (Emergoeconomic) Analyses of a Biomass Combustion Waste Heat Recovery Organic Rankine Cycle" Entropy 24, no. 2: 209. https://doi.org/10.3390/e24020209

APA Style

Effatpanah, S. K., Ahmadi, M. H., Delbari, S. H., & Lorenzini, G. (2022). Energy, Exergy, Exergoeconomic and Emergy-Based Exergoeconomic (Emergoeconomic) Analyses of a Biomass Combustion Waste Heat Recovery Organic Rankine Cycle. Entropy, 24(2), 209. https://doi.org/10.3390/e24020209

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