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

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×1012 sej/h, and an output power monetary emergy of 4.77×104 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).


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 ergoeconomic 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 kW e and 37.16 kW th 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.  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. Step 1. Thermodynamic modeling and energy analysis

System Description
Step 2. Exergy analysis and specifying the exergy flows Step 3. Equipment cost estimation and specifying cost per exergy unit Step 4. Identifying the required emergy related to the construction and investment cost of system components Step 5. Collecting and calculating required transformity values Step 6. Solving exergoeconomic balance equations and determine exergoeconomic parameters Step 7. Solving emergoeconomic balance equations and determine emergoeconomic parameters

Final
Step. Performing sensitivity analysis Emergy-based exergoeconomic (Emergoeconomic) analysis approach implemented.  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.  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. Step 1. Thermodynamic modeling and energy analysis

System Description
Step 2. Exergy analysis and specifying the exergy flows Step 3. Equipment cost estimation and specifying cost per exergy unit Step 4. Identifying the required emergy related to the construction and investment cost of system components Step 5. Collecting and calculating required transformity values Step 6. Solving exergoeconomic balance equations and determine exergoeconomic parameters Step 7. Solving emergoeconomic balance equations and determine emergoeconomic parameters

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.

Mathematical Modeling
The BCWHR-ORC system mentioned above has been modeled using EES software. The following assumptions were adapted to simplify the developed code: 1. Steady-state operation. 2. Negligible pressure loss in the condenser, the evaporator, and the piping. 3. Negligible heat loss through equipment. 4. 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.

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.

Mathematical Modeling
The BCWHR-ORC system mentioned above has been modeled using EES software. The following assumptions were adapted to simplify the developed code:

2.
Negligible pressure loss in the condenser, the evaporator, and the piping.

3.
Negligible heat loss through equipment.

4.
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. The mass and energy conservation equations (Equations (1) and (2)), excluding potential and kinetic energy terms, were used to study the system equipment. Table 2 represents the previous equations, rewritten for each piece of equipment. The thermal efficiency equation of the whole system is given below: .
The total heat transfer coefficient (UA) of the systems' heat exchangers (condenser and evaporator) is calculated using the following equations: Table 2. Mass and energy balance equations of the system equipment.

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: .
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): .
Additionally, the ratio of exergy destruction, which calculates the share of each piece of equipment in the total exergy destruction, is given below: where .
Finally, the efficiency of the second law (exergy efficiency) for each component and the whole system is calculated using Equations (12) and (13), respectively: The exergy efficiency and the destruction rate equation of the system components are included in Table 3. Table 3. The exergy destruction rate and efficiency equations of the system components.

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:

1.
Determining energy and exergy flows at the component's boundaries.

2.
Determining the fuel and product of each component. 3. Developing cost balance and auxiliary equations for the system components using Equation (14): (11) is rewritten as Equation (15): where c i , c e , c q and c w represent the cost per unit of exergy in USD/GJ. Z k , certain parameters such as the capital return factor (CRF), O&M factor (ϕ r ), annual operating hours (OH), and capital cost of the k th equipment (PEC k ), are to be defined. Finally, . Z k is calculated for all components using Equations (16) and (17) [48]: 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: logC 0,HE = K 1,HE + K 2,HE (logA HE ) + K 3,HE (logA HE ) 2 (19) logF P,HE = C 1,HE + C 2,HE (logP HE ) + C 3,HE (logP HE ) 2 (20) logF P,P = C 1,P + C 2,P (log 10 P P ) + C 3,P (log 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: 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. Table 4. Constant Values of the equipment's capital cost equations.

Constants Equipment
Heat Exchangers Turbine Pump 3.5000 1.5000 Table 5. Cost balance and auxiliary equations of the system equipment.

Component Cost Balance Equation Auxiliary Equations
Evaporator 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 Tot,k ), relative cost difference (r k ), and the exergoeconomic factor ( f k ), which are defined respectively in Equations (27) . .
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].

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): 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: In Equation (33), Tr biomass 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, η biomass is the ideal (Carnot) efficiency of the power plant, and T is the working fluid temperature in the evaporator (K). The Tr biomass and η biomass values are determined from [58,59].

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: 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: Accordingly, the monetary emergy rate of exergy destruction is obtained from Equation (38): .
. U k is the investment-related monetary emergy rate for each component given as follows: .
EM construct,k = Mass k × EM material (41) where EM construct,k , and EM Purchase,k , are the emergy content relevant to the construction and purchasing of the k th component in sej, Mass k is its respective mass in grams, EM material is the emergy of the construction material (sej/gr), and EM dollar 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. Table 7. The auxiliary and emergy balance equations of the system components.

Emergy-Based Cost Balance Equations Auxiliary Equation
Evaporator The total emergy rate of a given component is obtained through the summation of the monetary emergy rates related to both the investment ( The relative monetary emergy between corresponding average values of the product and the fuel is calculated using the following equation: Additionally, Equation (45) is used to calculate the emergoeconomic factor of each system component:

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].

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.

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.

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.

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/hr and 24.13 USD/GJ, respectively. Additionally, the cost rate and cost rate per exergy unit calculated for the system were 52.45 USD/hr, 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. ̇ , 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 ̇ and ̇ , of each component. From an exergoeconomic point of view, the higher the total cost rate ̇ 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 ̇ , is of interest for performance improvement in this study. A high ̇ 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 ̇ 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

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  C Total 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 Total , is of interest for performance improvement in this study. A high . C Total 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 Total 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.

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 total ) 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 Total 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 Total 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.    U Evap 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 Total,Pump ). 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 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 timeconsuming 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: Figure 9 depicts the trends of , , and while changing their values ±50% with 10% increments. As can be seen, was the most sensitive and was inversely correlated to changes, which is the main driving force of the system.

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 timeconsuming 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:

U Pump
× 100 (46) Figure 9 depicts the trends of Tr BCW H , Tr emdollar , and Tr Steel while changing their values ±50% with 10% increments. As can be seen, ψ was the most sensitive and was inversely correlated to Tr BCW 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 Tr BCW H , as a −50% change in the latter resulted in only a 3% percent decrease in the former. Contrastingly, the Tr emdollar , and Tr Steel had a direct correlation with ψ. Tr Steel 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. Tr emdollar changes also impacted the ψ, since the capital investment constituted the principal share of the equipment's total emergy content. The ±50% change in Tr emdollar , resulted in −1.83% to 1.58% changes in the ψ value.
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 Tr , as a −50% change in the la er resulted in only a 3% percent decrease in the former. Contrastingly, the , and had a direct correlation with . 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.
changes also impacted the , since the capital investment constituted the principal share of the equipment's total emergy content. The ±50% change in , resulted in −1.83% to 1.58% changes in the value. Figure 9. The trend of overall emergoeconomic performance of the system ( ), with respect to the transformity coefficients.

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: 1. The highest rate of exergy destruction ( ̇ ) 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. 2. According to the exergoeconomic analysis, the cost per exergy unit of the turbine's power ( ) was equal to 24.13 USD/GJ, and its cost rate of output power ( ) was 14.19 USD/h. Additionally, ̇ and for the whole system were 52.45 USD/h and 91.47 USD/GJ, respectively. 3. Amongst the system equipment, the turbine had the highest total cost rate ( ̇ ) and exergoeconomic factors ( ), with values of 6.1715 USD/h and 84.87%, respectively. It has been deduced that the cost rate of investment for the turbine ( ̇ ) is large, and measures should be taken to reduce its capital investment cost. 4. Conducting an emergoeconomic analysis, the monetary emergy per exergy unit (m) and the the monetary emergy rate of the output power ( ̇) were 8.24 × 10 sej/J, and 4.87 × 10 sej/h, respectively. Figure 9. The trend of overall emergoeconomic performance of the system (ψ), with respect to the transformity coefficients.

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: 1.
The highest rate of exergy destruction ( . Ex 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.

2.
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,