Exergy Based Ecological Performance Analysis of a Waste Heat-Powered Marine Refrigeration Cycle
PDF
Cite
Share
Request
Research Article
VOLUME: 227 ISSUE: 1
P: 34 - 42
May 2025

Exergy Based Ecological Performance Analysis of a Waste Heat-Powered Marine Refrigeration Cycle

J Nav Archit Mar Technol 2025;227(1):34-42
1. Yıldız Technical University Faculty of Naval Architecture and Maritime Deparment of Naval Architecture and Marine Engineering, İstanbul, Türkiye
No information available.
No information available
Received Date: 22.11.2024
Accepted Date: 18.03.2025
Online Date: 30.05.2025
Publish Date: 30.05.2025
PDF
Cite
Share
Request

Abstract

The escalating costs of fuel, global fuel supply disruptions, and stringent international emissions regulations necessitate more efficient ship operations. A key component of environmentally sustainable, economical, and efficient marine transportation is the recovery and utilization of waste heat. This study investigates the potential of converting waste heat from ship exhaust gas into useful energy through an organic Rankine cycle and its subsequent integration into a refrigeration system. In this study, the thermodynamic performance of the refrigerants R717, R152a, R290, and R134a was investigated. The findings demonstrate that R717 exhibits superior performance in terms of utilization factor and second law efficiency (ηII), while R290 emerges as the optimal choice when minimizing total entropy generation (sgen_TOT) is the primary objective. If the priorities are focused on maximizing useful work production, minimizing energy losses, and reducing environmental impacts, the advantages of the R717 refrigerant, which has the highest total ecological coefficient of performance (ECOP) and exergetic performance criteria (EPC) values, are clearly evident. The ideal operating temperatures for evaporator temperature (TEvap) and turbine temperatures (TTurbine) concerning EPC and ECOP, using R717 as the refrigerant, were investigated in the continuation of the study. During this investigation, changes in TEvap and TTurbine were observed to have affected other components. When the TEvap is kept constant, and the TTurbine amount is increased, an improvement was observed in EPCTOT, ECOPTOT, ηII, and utility factor values, while an increase in sgen_TOT and exergy destruction (ExDTOT) values. When the TTurbine is kept constant, a decrease is observed in the EPCTOT value. At the same time, an improvement is seen in ExDTOT, ηII, and ECOPTOT values, while an increase is noted in the sgen_TOT value.

Keywords:
ECOP, EPC, heat-powered refrigeration, temperature management, waste heat recovery

1. Introduction

Contemporary engineering systems must prioritize cost-effectiveness, efficiency, and environmental sustainability. To reach these objectives, various strategies for enhancing efficiency and minimizing environmental impact are implemented throughout land and marine vehicles’ design and operational phases. In the maritime domain, common approaches include resistance reduction techniques, the adoption of alternative fuels, operational optimizations, design parameters such as energy efficiency design index and energy efficiency existing index, and the recovery and utilization of waste heat [1]. Ships often generate substantial waste heat from systems like main engine cooling, lubrication, and flue gas [2]. This heat can be harnessed for various applications, including direct heating of tanks or fuel, steam generation [3], power generation, or cooling through different thermodynamic cycles [4]. However, ships also face significant cooling demands for spaces, food stores, and transported cargo [5]. To meet these cooling needs, engineers explore alternative approaches beyond traditional vapor compression refrigeration (VCR) systems. These include absorption or organic Rankine cycle (ORC) driven VCR systems utilizing waste heat, as well as hybrid methods combining both approaches.

The existing literature encompasses a diverse range of investigations into ORC-driven combined cooling systems or dual-cycle configurations for marine applications and low-temperature heat sources. Comparative performance assessments of various heat-powered refrigeration cycles, including absorption cycles, multi-stage arrangements, and combined heat-power systems, have been reported [6]. Thermodynamic analyses of ORC-VCR systems utilizing low-temperature source heat have been conducted for different working fluids and design parameters, leading to insights into optimal fluid selection and equipment [7-9]. Optimization studies have focused on refining design and operational parameters of ORC-VCR systems based on varying working fluids and heat source temperatures [10, 11]. Exergy-based performance analyses have explored the impact of different working fluids, heat sources, and system layouts on ORC-VCR system efficiency [12-14]. Experimental investigations have been conducted to evaluate the performance of ORC-VCR systems recovering waste heat from internal combustion engines [15].

Exergetic performance analysis and optimization of air refrigeration cycles based on ecological coefficient of performance (ECOP) were studied by [16]. Analyses of exergy-based performance outputs, including exergetic performance criteria (EPC) analyses, of the multipurpose refrigeration system for different design conditions were performed by [17]. A recent research study was explored the optimization of simple Brayton refrigeration models using an exergy-based approach to improve ecological function by [18]. Exergy based thermodynamic analysis and optimizations of low temperature cascade refrigeration systems applications were presented by [19, 20]. The performance analysis of the heat-powered refrigeration system using marine waste heat energy on the second law of thermodynamics was carried out depends on the different design parameters by [21]. This study will evaluate the performance of a refrigeration system driven by an ORC utilizing waste heat from a ship’s diesel engine. The analysis will focus on key performance metrics, including utility factor (Uf), second law efficiency (ηII), total entropy generation (sgen_TOT), exergy destruction ratio (y), ECOP, and EPC. These metrics will be assessed under various operating conditions, considering different working fluids and different evaporator temperatures (TEvap), condenser, and turbine temperatures (TTurbine).

2. Thermodynamic Model

The waste heat-powered ORC-VCR system can be conceptualized as two interconnected subsystems that utilize a single working fluid, which can be applied in such as industrial waste heat recovery systems, combined cooling and power systems, data center cooling, or marine systems. The ORC is a thermodynamic cycle that converts heat energy into mechanical work, similar to a steam Rankine cycle but using an organic working fluid (e.g., R134a, R717, 290, etc) instead of water. The ORC component comprises a condenser, pump, waste heat boiler, and steam turbine. The VCR cycle is a standard refrigeration process that absorbs heat from a low-temperature or cooled source and rejects it at a higher temperature. The VCR subsystem includes an evaporator, compressor, condenser, and expansion valve. Figure 1 illustrates the flow diagram of the ORC-VCR system. The ORC effectively converts waste heat from the main engine exhaust gases into a net work, driving the pump and compressor. In cases where the main engine is not operating or cannot provide exhaust gas with sufficient thermal content, the boiler is activated and thus continuity in the cooling system is ensured. The VCR then leverages this energy to refrigerate the storage compartments.

The fundamental design parameters and associated assumptions are outlined below. The heat addition to the ORC, QH, is calculated based on exhaust gas properties, including mass flow rate and inlet temperature (Texhaust,in). The specific heat capacity at constant temperature (Cp) of the exhaust gases was determined using Equation 1, considering the reference inlet and outlet temperatures.

​​Cp = 956 + ​(​​0.3389   ​T​exhaust,in​​​)​​–​(​​2.476   ​10​​ −5​  ​T​exhaust, in​  2 ​​)​​​​    (1)

Heat transfer (QCon, QL, and QH) values are calculated from the basic balance equations. The Uf, defined by Equation 2, represents the ratio of the refrigeration or cooling load of the VCR system (QL) to the heat input to the ORC (QH).

The ηII calculated with Equation 3 as the ratio of Uf and Carnot Uf (Uf_Carnot).

sgen_TOT, as defined by Equation 6, is the sum of the entropy generations from the environment, Equation 4, (Δsenv) and the system, Equation 5, (Δssys). The temperature difference between relevant environments and system points is denoted by ΔT. The heat transfer from the condenser unit (QCon) is the combined heat transfer from the power loop (PL) and the refrigeration loop (RL). TCon and Tavr represent the temperature of condenser for power loop and the average value between condenser inlet and outlet temperature of RL, respectively.

The exergy density (ρex) is a beneficial tool to compare the size and exergetic capacity of the system point which was defined the ratio of specific exergy to the specific volume of the system point in Equation 8 by [22].

The exergy destruction (ExD) of components is calculated with Equation 9 which depends on the flow rate and physical exergy of inlet and outlet conditions of components. The total ExDTOT is the sum of all components of ExD as in Equation 10.

The y can be defined for each component as the ratio of the ExD of components to the ExDTOT of the system in Equation 11 [23].

ECOP is defined as ratio of power output to the loss rate of availability in Equation 12 [24] and EPC defines the ratio of total exergy output to the loss rate of availability with Equation 13 [25].

Some properties such as global warming potential, refrigerant concentration limit of working fluids (R717, R290, R134a, and R152a) according to the ASHRAE [26] were given in Table 1 Assumptions made for the analysis of the system were also shared in Table 2.

The thermodynamic analysis of the marine refrigeration system (MRS) was conducted for the following refrigerants: R717, R290, R134a, and R152a. The analysis relied on the following simplifying assumptions:

• All system components were assumed to operate under steady-state conditions.

• Chemical, kinetic, and potential energy, as well as their corresponding exergy terms, were disregarded.

• Pressure losses within the system’s pipelines were considered negligible.

• Heat transfer to or from the compressor and expansion valve was assumed to be insignificant.

• The expansion of refrigerants in the expansion valves was assumed to be isenthalpic.

3. Results and Discussion

To evaluate the performance of the system under different refrigerant conditions, a parametric study was conducted using the established model (Equations 1-13). Four commonly employed refrigerants (R717, R290, R134a, and R152a) were selected. Their corresponding performance metrics, including the Uf, the ηII, the ECOPTOT, the EPCTOT, the sgen_TOT, and the ExDTOT, were calculated for different turbine and TEvap. The Python 3.11.4 programming language, in conjunction with the CoolProp 8.3 library [27], was employed for the numerical simulations.

The results of this analysis, the information presented in Figure 2 compares the performance outputs of a waste heat-driven marine refrigeration cycle using four different refrigerants. It was observed that the optimal refrigerant choice is contingent upon the specific performance objective. R717 generally exhibits higher Uf, ηII, ECOPTOT, EPCTOT, and ExDTOT and lower sgen_TOT compared to R290, R152a, and R134a. This suggests that R717 is more efficient in converting waste heat into refrigeration, indicating that ammonia might be a slightly more efficient refrigerant in this specific application, despite its disadvantages for human health.

Temperature (T), pressure (P), specific enthalpy (h), specific entropy (s), mass flow rate (m), specific volume (v), specific exergy (ex), and exergy density (ρEx) values of system points for the R717 are presented in Table 3 as thermophysical properties. Our priorities are to maximize practical work production, minimize energy losses, and reduce environmental impacts. In this case, the advantages of R717 refrigerant, which has the highest total ECOP and EPC values, are apparent.

To evaluate the system’s overall performance, it would be helpful to calculate the performance outputs that would provide a comprehensive understanding of the system’s losses and efficiencies. Table 4 represents performance outputs such as ExD, EPC, y, and ηII for the heat-powered MRS components (steam turbine, compressor, condenser, evaporator, expansion valve, pump, and boiler). The table shows that the evaporator has the highest ExDTOT rate and the highest y value, with an almost 41% ratio, and the lowest ExD and y values belong to the pump unit. It can also be obtained from the table that the heat-related components (condenser, evaporator, and boiler) generate more ExDTOT. Regarding the EPC value, the highest values are for pump and turbine units; the lowest are for heat-related units. By implementing design modifications and operational adjustments to the evaporator, condenser, and boiler units, the system’s ExDTOT can be reduced, leading to a notable increase in overall efficiency.

Figure 3 shows the effect of the TTurbine on the system components’ EPC. As the TTurbine increases, the overall EPC generally decreases. Nevertheless, the individual EPC contributions of different components demonstrate different trends. Significant decreases in EPCTurbine, EPCBoiler, and EPCCondenser between 340 and 355 K temperature and slight decreases in EPCPump are also observed. EPCEvap, EPCCompressor, and EPCValve remain relatively constant, as the evaporator’s temperature does not directly affect the compressor’s and valve’s performance.

Figure 4 shows the effect of TEvap on the EPC of the system components. As the TEvap increases, the overall EPC generally increases. However, the individual EPC contributions of different components show different trends. EPCTurbine, EPCBoiler, and EPCPump remain relatively constant, as the turbine’s, boiler’s, and pump’s performance are not directly affected by the TEvap. EPCEvap and EPCCondenser increase slightly, EPCCompressor and EPCValve increases significantly with TEvap.

Figure 5 shows the relationship between the Uf and the total EPCTOT or EPCTOT where the different lines represent various combinations of evaporator and TTurbine. As the TEvap increases, the EPCTOT generally decrease for a given Uf, and the Uf increases for a given constant TTurbine. However, there is a trade-off between TEvap and cooling capacity. Higher TEvap can reduce the cooling capacity of the system. Increasing the TTurbine generally increases EPCTOT and Uf for constant TEvap. Since improving EPCTOT is the priority, the scenario in which It reaches its maximum improvement should be considered. Therefore, when the TEvap was set to 243.15 K and the TTurbine to 368.15 K, an improvement of nearly 100% in EPCTOT and almost 115% in Uf according to the values in TTurbine is 338.15 K.

Figure 6 depicts the effects of the evaporator and TTurbine on the relationship between ECOPTOT or ECOPTOT and EPCTOT. As the TEvap increases, the ECOPTOT generally increases for a given EPCTOT or at a constant TTurbine. However, increasing in TEvap at constant TTurbine decreases the EPCTOT. TEvap and cooling capacity are inversely correlated, such that an increase in one leads to a decrease in the other. Increasing the TTurbine generally leads to a increase in ECOPTOT for a constant TEvap. The ideal scenario is when the TEvap is 258.15 K, and the TTurbine is 368.15 K to maintain optimal improvements in both ECOPTOT and EPCTOT.

Figure 7 describes the relationship between EPCTOT and ηII of the system using R717. As the TEvap increases, ηII generally increases while EPCTOT decreases for constant TTurbine. Increasing the TTurbine generally increases both EPCTOT and ηII for a given TEvap. Since maximizing EPCTOT is the primary objective, the scenario yielding its highest improvement should be considered. Accordingly, when the TEvap was adjusted to 243.15 K and the TTurbine was 368.15 K, an enhancement of almost 170% in ηII was calculated by comparing it to the lowest scenario.

Figure 8 depicts the relationship between EPCTOT and ExDTOT. The different lines represent various combinations of evaporator and TTurbine. As the TEvap increases, the ExDTOT generally increases for a given EPCTOT, and EPCTOT decreases with higher TEvap at a constant TTurbine value. Increasing the TTurbine generally increases the ExDTOT and EPCTOT. By carefully selecting these temperatures, optimizing the system’s performance in terms of both ExDTOT and EPCTOT is possible. While improving EPCTOT, to keep ExDTOT at its lowest level, two scenarios emerge: The most optimal balance between both values occurs when the TEvap is 243.15 K and the TTurbine is 338.15 K. However, if the increase in ExDTOT can be overlooked, the ideal temperature range for maximizing EPCTOT is when the TEvap is 243.15 K and the TTurbine is 368.15 K.

Figure 9 shows the relationship between EPCTOT and sgen_TOT for the system according to the different TEvap (straight lines) and TTurbine (dotted lines). As the TEvap increases, the sgen_TOT generally increases for a given EPCTOT. However, EPCTOT also decreases with higher TEvap. Higher TEvap can lead to higher irreversibilities in the system, such as heat transfer across finite temperature differences and pressure drops in the components. Increasing the TTurbine generally leads to an increase in sgen_TOT for a given EPCTOT. Besides, EPCTOT also increases with higher TTurbine. The optimal combination of evaporator and TTurbine depends on the specific requirements of the heat-powered refrigeration system, which must balance between minimizing sgen_TOT and maximizing EPCTOT. For example, lower TTurbine and TEvap might be suitable if the priority is to minimize entropy generation. If the priority is to maximize EPCTOT, a lower v and a higher TTurbine might be better. It has been observed that the ideal scenario between EPCTOT and sgen_TOT occurs when the TEvap is 243.15 K, and the TTurbine is 368.15 K.

4. Conclusion

This study investigated the feasibility of converting waste heat from ship exhaust gas into helpful energy through an ORC and subsequently utilizing this converted energy in a VCR cycle. The effects of the varying evaporator and condenser temperatures and different working fluids on the Uf, ηII, and sgen_TOT production were analyzed. The results indicated that ammonia (R717) is the most suitable fluid for this system. Future research endeavors will focus on conducting detailed comparative performance analyses and exploring exergy-based environmental and economic assessments. Analyses have shown that the high ECOP and EPC values of R717 (ammonia) indicate their significant potential to enhance system efficiency. These results suggest that R717 could play an effective role in energy conversion processes, particularly optimizing energy efficiency. Furthermore, the high efficiency of R717 could lead to lower emissions and more sustainable energy solutions from an environmental perspective, while economically, it could reduce operating costs by lowering energy consumption. Future studies will further explore the broader application areas of R717 and its impact on environmental and economic assessments.

The optimal combination of e TEvap and TTurbine depends on the specific requirements of the refrigeration system. For example, if the priority is to maximize cooling capacity, a lower TEvap and a higher TTurbine might be suitable. If the priority is to minimize energy consumption, a higher TEvap and a lower TTurbine might be better. Figures 5-9 show that the EPCTOT of the heat-powered system using R717 is influenced by both the evaporator and TTurbine. To optimize the total EPC, it is advantageous to operate at a lower TEvap and a higher TTurbine. The observations indicate that the optimal thermodynamic conditions for the EPCTOT, occur at an TEvap of 243.15 K and a TTurbine of 368.15 K which have almost 100% increment according to the lowest EPCTOT value. By carefully selecting these temperatures, optimizing the system’s performance for a given application is possible.

Acknowledgement

The short version of this study was presented at the Fourth International Naval Architecture and Maritime Symposium (INTNAM-2023, YTÜ, İstanbul, 11-13.10.2023), the present paper is the improved version in terms of analysis and content. We would like to thank those who contributed to the symposium organisers and the journal publication process.

Authorship Contributions

Concept: O. Yılmaztürk, Design: O. Yılmaztürk, Data Collection or Processing: O. Yılmaztürk, and A. S. Karakurt, Analysis or Interpretation: A. S. Karakurt, Literature Review: O. Yılmaztürk, and A. S. Karakurt, Writing: O. Yılmaztürk, and A. S. Karakurt.
Conflict of Interest: No conflict of interest was declared by the authors.
Financial Disclosure: The authors declared that this study received no financial support.

References

1
L. A. Díaz-Secades, R. González, N. Rivera, E. Montañés, and J. R. Quevedo, “Waste heat recovery system for marine engines optimized through a preference learning rank function embedded into a Bayesian optimizer,”Ocean Engineering, vol. 281, 114747, Aug 2023.
2
A. S. Karakurt, B. Başaran, and E. Z. Çubukçu, “Performance analyses of marine thermal oil system,”Gemi ve Deniz Teknolojisi, no. 222, pp. 196-208, Jan 2023.
3
B. Khan, and M.-H. Kim, “Energy and exergy analyses of a novel combined heat and power system operated by a recuperative organic rankine cycle integrated with a water heating system,”Energies, vol. 15, no. 18, 6658, Jan 2022.
4
T. Koroglu, I. Cavus, and M. U. Bayer, “Marine engine efficiency improvement with supercritical-CO2 Rankine waste heat recovery,”International Journal of Exergy, vol. 40, no. 3, pp. 347-364, Jan 2023.
5
V. Başhan, and A. Parlak, “Economic analysis of a ship refrigeration system in case of variable sea water temperature conditions,”Journal of ETA Maritime Science, vol. 3, no. 2, pp. 67-74, 2015.
6
R. Tozer, and R. W. James, “Heat powered refrigeration cycles,”Applied Thermal Engineering, vol. 18, no. 9, pp. 731-743, Sep 1998.
7
S. Khatoon, N. M. A. Almefreji, and M.-H. Kim, “Thermodynamic study of a combined power and refrigeration system for low-grade heat energy source,”Energies, vol. 14, no. 2, pp. 410, Jan 2021.
8
A. Elbir, F. A. Kodaloglu, I. Ucgul, and M. E. Sahin, “Thermodynamic analysis of refrigerants used in ORC-VCC combined power systems for low temperature heat sources,”Thermal Science, vol. 26, no. 4 Part A, pp. 2855-2863, 2022.
9
R. Zhar, A. Allouhi, A. Jamil, and K. Lahrech, “Performance analysis of combined power and refrigeration: ORC-VCC system,” inWITS 2020, S. Bennani, Y. Lakhrissi, G. Khaissidi, A. Mansouri, and Y. Khamlichi, Eds., in Lecture Notes in Electrical Engineering. Singapore: Springer, 2022, p. 595-605.
10
Y. Liang, Z. Yu, and W. Li, “A waste heat-driven cooling system based on combined organic rankine and vapour compression refrigeration cycles,”Applied Sciences, vol. 9, no. 20, 4242, Jan 2019.
11
X. Xia, Z. Liu, Z. Wang, T. Sun, and H. Zhang, “Multi-layer performance optimization based on operation parameter-working fluid-heat source for the ORC-VCR system,”Energy, vol. 272, 127103, Jun 2023.
12
O. Bounefour, A. Ouadha, and Y. Addad, “An exergy analysis of various layouts of ORC-VCC systems for usage in waste heat recovery onboard ships,”Marine Systems & Ocean Technology, vol. 15, no. 1, pp. 26-44, Mar 2020.
13
C. Aktemur, and G. Hacipasaoglu, “Assessment of an integrated organic rankine cycle (ORC)-vapor compression refrigeration (VCR) system using the energy, conventional exergy, and advanced exergy analysis,”Heat Transfer Research, vol. 52, no. 15, 2021.
14
M.-H. Kim, “Energy and exergy analysis of solar organic rankine cycle coupled with vapor compression refrigeration cycle,”Energies, vol. 15, no. 15, Jan 2022.
15
Y. Liang, A. Mckeown, Z. Yu, and S. F. K. Alshammari, “Experimental study on a heat driven refrigeration system based on combined organic Rankine and vapour compression cycles,”Energy Conversion and Management, vol. 234, 113953, Apr 2021.
16
Y. Ust, “Performance analysis and optimization of irreversible air refrigeration cycles based on ecological coefficient of performance criterion,”Applied Thermal Engineering, vol. 29, no. 1, pp. 47-55, Jan 2009.
17
Y. Ust, A. S. Karakurt, and U. Gunes, “Performance analysis of multipurpose refrigeration system (MRS) on fishing vessel,”Polish Maritime Research, vol. 23, no. 2, pp. 48-56, 2016.
18
L. Chen, Y. Ge, S. Shi, H. Feng, and P. Liu, “Exergetic efficiency and exergy-based ecological function performance optimizations for two irreversible simple Brayton refrigeration cycle models,”Results in Engineering, vol. 22, 101967, Jun 2024.
19
Md. M. Arefin, D. Mondal, and Md. A. Islam, “Optimizing cascade refrigeration systems with low GWP refrigerants for low-temperature applications: A thermodynamic analysis,”Energy Conversion and Management: X, vol. 24,, Oct 2024.
20
Y. Ust, and A. S. Karakurt, “Analysis of a cascade refrigeration system (CRS) by using different refrigerant couples based on the exergetic performance coefficient (EPC) criterion,”Arabian Journal for Science and Engineering, vol. 39, no. 11, pp. 8147-8156, Nov 2014.
21
O. Yilmazturk, and A. S. Karakurt, “Analyses of an ideal combined heat-powered marine refrigeration system,” presented at the 4th International Naval Architecture and Maritime Symposium (Int-Nam 2023), İstanbul, Oct 2023. Accessed: Sep 13, 2024. [Online]. Available: https://avesis.yildiz.edu.tr/yayin/4e11e10f-2bff-4fea-b347-43e870e4a316/analyses-of-an-ideal-combined-heat-powered-marine-refrigeration-system
22
A. S. Karakurt, and B. Sahin, “A new method for the size and performance analyses and optimization of thermal systems: The exergy density,”Sigma Journal of Engineering and Natural Sciences, vol. 37, no. 2, pp. 573-583, Jul 2019.
23
A. Bejan, “Entropy generation minimization: The new thermodynamics of finite‐size devices and finite‐time processes,”Journal of Applied Physics, vol. 79, no. 3, pp. 1191-1218, Feb 1996.
24
Y. Ust, B. Sahin, and A. Kodal, “Ecological coefficient of performance (ECOP) optimization for generalized irreversible Carnot heat engines,”Journal of the Energy Institute, vol. 78, no. 3, pp. 145-151, Aug 2005.
25
Y. Ust, B. Sahin, and T. Yilmaz, “Optimization of a regenerative gas-turbine cogeneration system based on a new exergetic performance criterion: Exergetic performance coefficient,”Proceedings of the Institution of Mechanical Engineers Part A-Journal of Power and Energy, vol. 221, no. 4, pp. 447-456, Jun 2007.
26
ASHRAE,Designation and Safety Classification of Refrigerants, Standart 34, Georgia, US., 2019. [Online]. Available: https://www.ashrae.org/technical-resources/bookstore/standard-34-designation-and-safety-classification-of-refrigerants
27
I. H. Bell, J. Wronski, S. Quoilin, and V. Lemort, “Pure and pseudo-pure fluid thermophysical property evaluation and the open-source thermophysical property library CoolProp,”Industrial & Engineering Chemistry Research, vol. 53, no. 6, pp. 2498-2508, Feb 2014.