Numerical Investigations on Supercritical Heat Transfer Characteristics of Environmental Friendly Refrigerants

Adnan Ibrahim¹, Peng Hu^1*, Yiran Jiang¹, Ali Riaz², Farrukh Saleem³, Yifang Dong³, Lei Jia^4*, Panpan Zhao⁴

¹Department of Thermal Science and Energy Engineering, University of Science and Technology of China, Hefei, China.
²Department of Mechanical Engineering, Pakistan Institute of Engineering and Applied Sciences, Islamabad, Pakistan.
³School of Engineering Science, University of Science and Technology of China, Hefei, China.
⁴Hefei General Machinery Research Institute Co., Ltd., Hefei, China.

DOI: 10.4236/ijcce.2021.102002                

Abstract

The heat transfer of supercritical fluids is a vastly growing field, specifically to find suitable alternatives to replace conventional R134a, which can be beneficial for climate change. Most of the experimental and numerical investigations have been conducted to explore supercritical water, carbon dioxide and R134a as heat transfer working fluids. Hydrofluoroolefin (HFO) and refrigerants blends have been considered the most environment-friendly refrigerants to replace Chlorofluorocarbons (CFCs), Hydrochlorofluoro-carbons (HCFCs) and Hydrofluorocarbons (HFCs). Their main advantage of zero Ozone Depletion Potential (ODP) and comparatively lower Global Warming Potential (GWP) have attracted growing amount of attention to mitigate environmental issues. This work adopts the computational method and takes the environmentally friendly refrigerants to investigate the heat transfer characteristics under widely used shear-stress transport (SST) model. A comprehensive comparison was performed at reduced pressure of 1.10 for supercritical fluids R515A, R1234ze(E) and R134a. The peaks of heat transfer coefficient occurred in the vicinity of pseudo critical temperature for all of these considered fluids; however, R134a resulted in higher heat transfer coefficient, Reynolds number and Prandtl number in comparison with R515A and R1234ze(E). The higher heat transfer coefficient of supercritical fluid R134a is owing to its thermophysical properties and the specific heat plays crucial role in the heat transfer of supercritical fluids. Owing to environmental issues, R515A can be a considerable replacement of R134a. R1234ze(E) is also promising alternative to R134a; however, safety issues should thoroughly concern its mild flammable characteristics.

Keywords

Heat Transfer, Supercritical Fluids, Simulations, R515A, R1234ze(E), Environmental Friendly Refrigerants

1. Introduction

Supercritical fluids have wide range of industrial applications owing to the substantial impact of their heat transfer characteristics [1] - [8]. Supercritical fluids, in comparison with conventional fluids, have attracted growing amount of attention because of relatively higher heat transfer rate and lower exergy losses [9]. The thermophysical properties of supercritical fluids considerably vary near the critical (T_c) or pseudo critical temperature (T_pc). The heat transfer coefficient, owing to dramatic variations near pseudo-critical point, depends upon pressure, tube diameter, flow direction, heat flux and type of working fluid [10]. Therefore, this results in complex heat transfer characteristics which may account for heat transfer enhancement or deterioration [11].

Most of the experimental and numerical investigations have been conducted to explore supercritical water and carbon dioxide (sCO₂) [12] - [18]. Dang and Hihara [19] [20] investigated the effects of tube diameter on heat transfer coefficient of sCO₂ under cooling conditions and proposed the modified Gnielinski equation. Zhang and Hu [2] measured the effects of buoyancy and tube diameter for sCO₂. The influence of mass flux, pressure and tube diameter were plotted against heat transfer coefficient and pressure drop. Further to that, dimensionless diameter was incorporated in the development of correlation which can precisely estimate heat transfer in large-diameter tube. Wang and Guan [21] computationally investigated underlying mechanism of buoyancy effects for supercritical carbon dioxide flowing through large horizontal tube. At higher heat flux, the buoyancy is more pronounced and can cause considerable difference in the temperature at top and bottom walls. Zahlan and Groeneveld [22] performed extensive experimental tests for sCO₂ under vertical conditions.

However, supercritical organic fluids have not been thoroughly investigated for in tube heat transfer. Zhao and Jiang [23] examined that fluid temperature, mass flux and pressure can considerably impact the in-tube cooling heat transfer and flow of supercritical fluid R134a. Experimental data predicated (using the least square curve-fitting method) a modified Gnielinski’s correlation which can give heat transfer coefficient within ±15% accuracy. Wang and Tian [24] conducted experimental investigations for supercritical fluid R134a flowing through micro-fin and smooth tube under horizontal position. These measurements under different mass fluxes, heat fluxes and pressures suggested that micro-fin tube resulted in higher heat transfer coefficient than that of smooth tube. Herein,

buoyancy criteria of $\frac{G{r}_b}{{\mathrm{Re}}_b^2}\left(\frac{{\rho }_b}{{\rho }_w}\right)\frac{x}{d}$ were suggested to accurately predict results.

Further to that, micro-fin tube can significantly reduce the buoyancy effects. In more recent work, Wang and Tian [25] suggested that internally ribbed tube resulted in higher heat transfer coefficient than that of smooth tube under similar working conditions.

Kang and Chang [26] performed experiments for steady-state and transient-pressure in upward flow of supercritical fluid R134a. The study suggested that pressure transient rates have slight impact upon heat transfer rate. Cui and Wang [27] experimentally examined supercritical fluid R134a for different flow directions in a vertical tube. The data suggested good heat transfer in downward flow as compared to upward direction. He and Dang [28] [29] experimentally investigated supercritical fluid R245fa in vertical tube under heating condition. The experimental results revealed 70% data can be calculated by Yamagata’s correlation within ±30% accuracy. The experimental data of supercritical fluid R1233zd(E) showed good agreement with Petukhov’s correlation. In comparison with supercritical fluid R245fa, supercritical fluid R1233zd(E) can bring higher heat transfer coefficient. Jiang et al. [30] compared supercritical fluid R-22 and ethanol using smaller tube (1.004 mm) under higher heat flux (110 - 1800 kW·m⁻²). Ethanol was suggested for better flow and heat transfer performance; therefore, it’s reasonable for cooling applications in combustion chambers.

Xiong and Gu [31] performed experiments and numerical simulations to evaluate the intermittent heating effects for supercritical fluid R134a. After analyzing experimental data and simulation models, SST k-ω model was suggested to accurately predict heat transfer enhancement as well as heat transfer deterioration. The decrease in velocity for near-wall region can cause heat transfer deterioration. Liu and Xu [32] compared nine turbulence models with experimental results of sCO₂ passing through helical tube and suggested the Shear Stress Transport model for best prediction to heat transfer characteristics. The comparisons of various turbulent models were performed in previous research works for different supercritical fluids including sCO₂ [32] [33] [34] [35] [36], supercritical water [37] - [42], supercritical methane [43], supercritical nitrogen [44], supercritical fluid R134a [31] [45] [46] and supercritical fluid R1234ze(E) [1]. These findings suggested good agreement between simulations (performed by SST k-ω model) and experimental data. This model can provide most accurate prediction to heat transfer coefficient, wall and bulk temperatures [36]; therefore, the present simulations of supercritical fluid R515A were performed using SST k-ω model.

R515A is non-flammable and azeotrope replacement of R134a [47], and the mixture information is shown in Table 1. It has a lower global warming potential (GWP) of 403 than that of R134a (1300 GWP of R134a). R515A /R1234yf system was suggested to lower emissions and increase energy efficiency as compared to R744 system [48]. In this context, the potential alternatives are being investigated regarding heat transfer behavior, cyclic performance and safety problems. The heat transfer characteristics are almost similar. Moreover, the environmental issues caused by R134a can be mitigated with these potential eco-friendly alternatives.

Table 1. Refrigerants information.

However, most of the studies have explored their flow boiling or condensation heat transfer. There are few investigations at supercritical pressure, which attempted to evaluate the feasibility of HFO. R134a has been widely used for the organic Rankine cycle due to its better thermal performance. But global environmental issues demand to replace of the existing R134a with environment-friendly refrigerants which can serve potential alternatives. In this regard, R515A and R1234ze(E) have been compared with R134a at reduced pressure (P/P_cr = 1.10). This work is a step forward to study and explore the environment-friendly refrigerants. The simulations performed in this study can provide details about heat transfer of supercritical fluids under different mass fluxes, heat flux and tube diameters.

2. Numerical Simulations

2.1. Physical Model

Thermophysical properties of supercritical fluids vary considerably near pseudo-critical point, as shown in Figure 1. The thermophysical properties of supercritical fluids R515A, R1234ze(E) and R134a at different temperatures and pressures were taken from REFPROP 9.1 and were plotted using OriginPro tools. Therefore, it is crucial to investigate the supercritical heat transfer in the vicinity of T_pc under different pressure rates. A 3D physical model is employed in the simulations to consider the effects of buoyancy for supercritical fluids, as shown in Figure 2. Most of the commercial heat exchangers, which are employing organic Rankine cycle, are using horizontal flow direction rather than vertical [1] [2] [25]. Therefore, the present simulations adopted horizontal flow to explore the heat transfer. An adiabatic section (200 mm) is considered to eliminate the entrance effect, and constant heat flux boundary (q) is used for the wall (1000 mm) with different diameters. The relevant refrigerants information regarding critical parameters has been described in Table 1.

2.2. Mathematical Model

The detailed mathematical model is described below [35].

The continuity equation is described as:

Figure 1. Supercritical R515A, R1234ze(E) and R134a at reduced pressure of 1.10 (a) Density (b) Specific heat (c) Thermal conductivity (d) Viscosity (REFPROP 9.1).

Figure 2. 3D physical model.

$\frac{\partial }{\partial x_i}\left(\rho u_i\right)=0$ (1)

The momentum equation is described as:

$\frac{\partial }{\partial x_j}\left(\rho u_i{u}_j\right)=\frac{\partial }{\partial x_j}\left[{\mu }_{eff}\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right)-\frac{2}{3}{\mu }_{eff}\frac{\partial u_k}{\partial x_k}\right]-\frac{\partial p}{\partial x_j}+\rho g_i$ (2)

The energy equation is described as:

$\frac{\partial }{\partial x_i}\left(\rho u_i{c}_{p}T\right)=\frac{\partial }{\partial x_i}\left(\lambda \frac{\partial T}{\partial x_i}\right)+\Phi$ (3)

where ${\mu }_{eff}$ describes effective viscosity coefficient, and $\Phi$ describes energy dissipation.

The turbulent kinetic energy equation is described as [1] [35]:

$\frac{\partial }{\partial t}\left(\rho k\right)+\frac{\partial }{\partial x_i}\left(\rho k{u}_i\right)=\frac{\partial }{\partial x_j}\left({\Gamma }_k\frac{\partial k}{\partial x_j}\right)+G_k-Y_k+S_k$ (4)

The dissipation rate equation is described as:

$\frac{\partial }{\partial t}\left(\rho \omega \right)+\frac{\partial }{\partial x_i}\left(\rho \omega u_i\right)=\frac{\partial }{\partial x_j}\left({\Gamma }_{\omega }\frac{\partial \omega }{\partial x_j}\right)+G_{\omega }-Y_{\omega }+D_{\omega }+S_{\omega }$ (5)

where $G_k$ and $G_{\omega }$ denotes the generation of k and $\omega$ , ${\Gamma }_k$ and ${\Gamma }_{\omega }$ denotes the effective diffusivity of k and $\omega$ , respectively, $Y_k$ and $Y_{\omega }$ denotes the dissipation of k and $\omega$ due to turbulence, $D_{\omega }$ defines the cross-diffusion term, $S_k$ and $S_{\omega }$ are user-defined source terms.

2.3. Boundary Conditions

ANSYS FLUENT was employed for 3D simulation of turbulent flow. The thermophysical properties of supercritical fluids at different temperatures were taken from REFPROP 9.1 and input by piecewise-linear function. SST model was adopted for present simulations owing to relatively accurate results for a range of supercritical fluids. This model has been widely used for predicting reliable results. The detailed working conditions are described in Table 2. The reference values including inlet velocity are computed from inlet for each case using ANSYS FLUENT. The following boundary conditions were adopted: mass flow inlet, outflow boundary, and constant wall heat flux. SIMPLE algorithm is used for pressure and velocity coupling.

The bulk temperature and heat transfer coefficient were calculated as follows:

$T_b={\displaystyle {\int }_0^A\rho uT\text{d}A}/{\displaystyle {\int }_0^A\rho u\text{d}A}$ (6)

$h=\frac{q}{T_b-T_w}$ (7)

where T_b is the bulk temperature, T_w is the wall temperature, u is the local velocity and A is the cross-sectional area of the tube.

2.4. Mesh Independence Verification and Model Validation

ANSYS ICEM is used to generate high-quality hexahedral mesh as shown in Figure 3. Keeping all the working conditions same, h is plotted for different mesh sizes as illustrated in Figure 4. The deviation in h values obtained from different mesh sizes is trivial and further details have been described in Table 3. A reasonable compromise is to use mesh 2 for further simulations which can bring satisfactory accuracy and calculation speed.

The model verification is performed against experimental data presented by Dang and Hihara [19] and Jiang and Hu [1]. The present simulations resulted in a reliable heat transfer performance and better consistency with the experimental results (Figure 5) and can be employed for supercritical fluids.

Table 2. Working conditions considered for comparison analyses.

Table 3. Mesh independence for different cell numbers.

Figure 3. Details of mesh.

Figure 4. Mesh independence for different cell numbers.

3. Results

R1234yf and R1234ze(E) have been suggested as drop-in replacements for R134a in many previous research works. The performance parameters and heat transfer characteristics are almost similar. Moreover, the environmental issues caused by R134a can be mitigated with these potential eco-friendly alternatives [51] - [56]. In this section, heat transfer characteristics at supercritical pressure are investigated using wider range of operating conditions. The comprehensive comparison was performed at reduced pressure of 1.10 for supercritical fluids R515A, R1234ze(E) and R134a, as shown in Figure 6-8.

Figure 5. Model validation by comparing previous experimental data (a) Dang and Hihara [19] and (b) Jiang and Hu [1].

Figure 6. Comparison of heat transfer coefficient for supercritical fluids R515A, R1234ze(E) and R134a.

Figure 7. Comparison of T_w for supercritical fluids R515A, R1234ze(E) and R134a.

Figure 8. Comparison of the dimensionless numbers for supercritical fluids R515A, R1234ze(E) and R134a.

The peaks of h occurred in the vicinity of T_pc for these considered fluids, however, R134a resulted in higher heat transfer coefficient in comparison with R515A and R1234ze(E). Herein, specific heat plays crucial role in the heat transfer of supercritical fluids. R134a has comparatively higher specific heat as shown in Figure 1. The higher heat transfer coefficient of supercritical fluid R134a is owing to its thermophysical properties.

In order to study the influence of flow rate on the supercritical heat transfer process, different flow rates were selected in the 9.44 mm diameter heat exchange tube, and their convective heat transfer coefficients were compared and analyzed. The effect of mass flow on the heat transfer process is shown in Figure 6(c) and Figure 6(d). When the heat flux and tube diameter of the supercritical heat exchange process remains unchanged, it can be seen that the supercritical heat transfer coefficient increases with the increase of the mass flow rate. When the diameter of the heat exchange tube is the same, the larger flow rate causes the Reynolds number to increase, so that the thickness of the boundary layer on the inner surface of the heat exchange tube decreases during fluid flow. As a result, the fluid heat transfer coefficient can be improved.

Tube geometry, concerning different diameter, was considered for further simulations. The heat transfer coefficient may slightly lower with relatively large diameter tube as demonstrated in Figure 6(b) and Figure 6(c). The ratio of boundary layer flow with total stream becomes lower in case of large diameter tube and this can slightly decrease the h. The considered range of heat flux in the present simulations showed slight impact upon h (Figure 6(a) and Figure 6(b)). The rest of the simulations were performed under heat flux of 50 kW/m².

The influential nondimensional parameters (Figure 8) have been illustrated to further comprehend the heat transfer ability. Al the nondimensional parameters showed considerable influence of pseudo critical temperature, specifically in the vicinity of T_pc. It can be seen from Figure 8(a) that these supercritical fluids are affected by the buoyancy in heat exchange tubes for the considered conditions. Richardson numbers are from 0.1 to 10, as a result, the heat transfer process was affected by natural convection and force convection.

R134a showed beneficial values of Reynolds number and Prandtl number under the same operating conditions. Owing to environmental issues, the present simulations suggest that R515A is a considerable replacement of R134a. Table 1 showed that the GWP of R134a is more than three times higher. Both R515A and R134a have safety classification of A1-non Flammable. R1234ze(E) is also promising alternative to R134a, however, safety issues should thoroughly concern the flammability.

4. Conclusions

Hydrofluoroolefin and refrigerants blends have been attracting huge attentions owing to their environment-friendly characteristics. These can replace CFCs, HCFCs and HFCs, because the heat transfer characteristics are almost similar. The present simulations attempted to investigate supercritical fluids R515A, R1234ze(E) and R134a under cooling conditions flowing through horizontal tube. Herein, this work investigated the influence of different heat fluxes, mass fluxes and tube diameters as follows:

· The increase in mass flux can enhance the heat transfer owing to increase in Reynolds number. The influence of G is considerably prominent around T_pc, specifically when the T_b is slightly higher than T_pc. Higher values of G resulted in increased Re with thin boundary layer, consequently, increase in heat transfer and higher h values.

· The 9.44 mm diameter tube showed slightly lowered heat transfer coefficient than that of 4.12 mm. For the considered range of heat flux, heat transfer coefficient changed slightly.

· A comprehensive comparison was performed at reduced pressure of 1.10 for supercritical fluids R515A, R1234ze(E) and R134a. The peaks of heat transfer coefficient occurred in the vicinity of pseudo critical temperature for all of these considered fluids, however, R134a resulted in higher heat transfer coefficient, Reynolds number, Prandtl number and Nusselt number in comparison with R515A and R1234ze(E).

· Owing to environmental issues, the present simulations suggest that R515A is a considerable replacement of R134a. R1234ze(E) is also promising alternative to R134a; however, safety issues should thoroughly concern its mild flammable characteristics.

Further investigations are required to thoroughly explore the heat transfer characteristics of potential alternatives in cooling and heating conditions. The effect of vertical tube and flow directions can also be considered. The simulations are reliable enough; however, experimental investigations under actual conditions can strongly affirm the results.

Acknowledgements

This work is supported by the National Natural Science Foundation of China (Grant No. 51576187).

Nomenclature

NOTES

*Corresponding authors, E-mail: [e-mail] (Peng Hu), [e-mail] (Lei Jia).

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Jiang, Y.-R., Hu, P. and Ibrahim, A. (2020) Experimental and Numerical Investigation on Heat Transfer Characteristics of Supercritical R1234ze (E) Cooled in Horizontal Tubes. International Journal of Heat and Fluid Flow, 85, Article ID: 108650. https://doi.org/10.1016/j.ijheatfluidflow.2020.108650
[2]	Zhang, G.-W., Hu, P., Chen, L.-X. and Liu, M.-H. (2018) Experimental and Simulation Investigation on Heat Transfer Characteristics of in-Tube Supercritical CO2 Cooling Flow. Applied Thermal Engineering, 143, 1101-1113. https://doi.org/10.1016/j.applthermaleng.2018.07.108
[3]	Ibrahim, A., Peng, H., Riaz, A., Basit, M.A., Rashid, U. and Basit, A. (2021) Molten Salts in the Light of Corrosion Mitigation Strategies and Embedded with Nanoparticles to Enhance the Thermophysical Properties for CSP Plants. Solar Energy Materials and Solar Cells, 219, Article ID: 110768. https://doi.org/10.1016/j.solmat.2020.110768
[4]	Tian, R., Wang, D., Zhang, Y., Ma, Y., Li, H. and Shi, L. (2019) Experimental Study of the Heat Transfer Characteristics of Supercritical Pressure R134a in a Horizontal Tube. Experimental Thermal and Fluid Science, 100, 49-61. https://doi.org/10.1016/j.expthermflusci.2018.08.027
[5]	Mishra, S., Nayak, M. and Misra, A. (2020) Thermal Conductivity of Nanofluids—A Comprehensive Review. International Journal of Thermofluid Science and Technology, 7, Article ID: 070301. https://doi.org/10.36963/IJTST.2020070301
[6]	Zhou, Y., Zhang, L., Bu, S., Sun, C., Xu, W., Xiao, Y. and Liu, L. (2020) Study on Heat Transfer Characteristics of the Whole Plate Fin Tube Cooler. International Journal of Thermofluid Science and Technology.
[7]	Islam, T., Parveen, N. and Asad, M. (2020) Hydromagnetic Natural Convection Heat Transfer of Copper-Water Nanofluid within a Right-Angled Triangular Cavity. International Journal of Thermofluid Science and Technology, 7, Article ID: 070304. https://doi.org/10.36963/IJTST.2020070304
[8]	Marques, S.P. and Campo, A. (2019) Finite Strip Method Applied to Steady Heat Conduction and Thermal Radiation in a Planar Slab: Absorbing-Emitting Gray Material and Parallel Diffuse Surfaces. International Journal of Thermofluid Science and Technology, 6, Article ID: 19060102.
[9]	Cui, Y., Wang, H. and Wang, Y. (2019) Experimental and Numerical Studies on Convective Heat Transfer of Supercritical R-134a in a Horizontal Tube. International Journal of Heat and Mass Transfer, 136, 34-45. https://doi.org/10.1016/j.ijheatmasstransfer.2019.02.083
[10]	Wang, D., Tian, R., Zhang, Y., Li, L., Ma, Y., Shi, L. and Li, H. (2019) Heat Transfer Investigation of Supercritical R134a for Trans-Critical Organic Rankine Cycle System. Energy, 169, 542-557. https://doi.org/10.1016/j.energy.2018.12.034
[11]	Ye, K., Zhang, Y., Yang, L., Zhao, Y., Li, N. and Xie, C. (2019) Modeling Convective Heat Transfer of Supercritical Carbon Dioxide Using an Artificial Neural Network. Applied Thermal Engineering, 150, 686-695. https://doi.org/10.1016/j.applthermaleng.2018.11.031
[12]	Ehsan, M.M., Guan, Z. and Klimenko, A. (2018) A Comprehensive Review on Heat Transfer and Pressure Drop Characteristics and Correlations with Supercritical CO2 under Heating and Cooling Applications. Renewable and Sustainable Energy Reviews, 92, 658-675. https://doi.org/10.1016/j.rser.2018.04.106
[13]	Huang, D., Wu, Z., Sunden, B. and Li, W. (2016) A Brief Review on Convection Heat Transfer of Fluids at Supercritical Pressures in Tubes and the Recent Progress. Applied Energy, 162, 494-505. https://doi.org/10.1016/j.apenergy.2015.10.080
[14]	Wang, H., Leung, L.K., Wang, W. and Bi, Q. (2018) A Review on Recent Heat Transfer Studies to Supercritical Pressure Water in Channels. Applied Thermal Engineering, 142, 573-596. https://doi.org/10.1016/j.applthermaleng.2018.07.007
[15]	Rahman, M.M., Dongxu, J., Beni, M.S., Hei, H.C., He, W. and Zhao, J. (2016) Supercritical Water Heat Transfer for Nuclear Reactor Applications: A Review. Annals of Nuclear Energy, 97, 53-65. https://doi.org/10.1016/j.anucene.2016.06.022
[16]	Zhang, S., Xu, X., Liu, C. and Dang, C. (2020) A Review on Application and Heat Transfer Enhancement of Supercritical CO2 in Low-Grade Heat Conversion. Applied Energy, 269, Article ID: 114962. https://doi.org/10.1016/j.apenergy.2020.114962
[17]	Knez, Ž., Markočič, E., Leitgeb, M., Primožič, M., Hrnčič, M.K. and Škerget, M. (2014) Industrial Applications of Supercritical Fluids: A Review. Energy, 77, 235-243. https://doi.org/10.1016/j.energy.2014.07.044
[18]	Huang, D. and Li, W. (2018) A Brief Review on the Buoyancy Criteria for Supercritical Fluids. Applied Thermal Engineering, 131, 977-987. https://doi.org/10.1016/j.applthermaleng.2017.12.042
[19]	Dang, C. ad Hihara, E. (2004) In-Tube Cooling Heat Transfer of Supercritical Carbon Dioxide. Part 1. Experimental Measurement. International Journal of Refrigeration, 27, 736-747. https://doi.org/10.1016/j.ijrefrig.2004.04.018
[20]	Dang, C. and Hihara, E. (2004) In-Tube Cooling Heat Transfer of Supercritical Carbon Dioxide. Part 2. Comparison of Numerical Calculation with Different Turbulence Models. International Journal of Refrigeration, 27, 748-760. https://doi.org/10.1016/j.ijrefrig.2004.04.017
[21]	Wang, J., Guan, Z., Gurgenci, H., Hooman, K., Veeraragavan, A. and Kang, X. (2018) Computational Investigations of Heat Transfer to Supercritical CO2 in a Large Horizontal Tube. Energy Conversion and Management, 157, 536-548. https://doi.org/10.1016/j.enconman.2017.12.046
[22]	Zahlan, H., Groeneveld, D. and Tavoularis, S. (2015) Measurements of Convective Heat Transfer to Vertical Upward Flows of CO2 in Circular Tubes at Near-Critical and Supercritical Pressures. Nuclear Engineering and Design, 289, 92-107. https://doi.org/10.1016/j.nucengdes.2015.04.013
[23]	Zhao, C.-R. and Jiang, P.-X. (2011) Experimental Study of In-Tube Cooling Heat Transfer and Pressure Drop Characteristics of R134a at Supercritical Pressures. Experimental Thermal and Fluid Science, 35, 1293-1303. https://doi.org/10.1016/j.expthermflusci.2011.04.017
[24]	Wang, D., Tian, R., Zhang, Y., Li, L. and Shi, L. (2019) Experimental Comparison of the Heat Transfer of Supercritical R134a in a Micro-Fin Tube and a Smooth Tube. International Journal of Heat and Mass Transfer, 129, 1194-1205. https://doi.org/10.1016/j.ijheatmasstransfer.2018.10.052
[25]	Wang, D., Tian, R., Li, L., Dai, X. and Shi, L. (2020) Heat Transfer of R134a in a Horizontal Internally Ribbed Tube and in a Smooth Tube under Super Critical Pressure. Applied Thermal Engineering, 173, Article ID: 115208. https://doi.org/10.1016/j.applthermaleng.2020.115208
[26]	Kang, K.-H. and Chang, S.-H. (2009) Experimental Study on the Heat Transfer Characteristics during the Pressure Transients under Supercritical Pressures. International Journal of Heat and Mass Transfer, 52, 4946-4955. https://doi.org/10.1016/j.ijheatmasstransfer.2009.06.005
[27]	Cui, Y.-L. and Wang, H.-X. (2018) Experimental Study on Convection Heat Transfer of R134a at Supercritical Pressures in a Vertical Tube for Upward and Downward Flows. Applied Thermal Engineering, 129, 1414-1425. https://doi.org/10.1016/j.applthermaleng.2017.10.120
[28]	He, J., Dang, C. and Hihara, E. (2018) Experimental Investigation of Heat Transfer to Supercritical R245fa Flowing Vertically Upward in a Circular Tube. International Journal of Heat and Mass Transfer, 127, 286-295. https://doi.org/10.1016/j.ijheatmasstransfer.2018.06.126
[29]	He, J., Dang, C. and Hihara, E. (2018) Supercritical Heat Transfer Characteristics of R1233zd(E) in Vertically Upward Flow. International Journal of Heat and Mass Transfer, 127, 497-505. https://doi.org/10.1016/j.ijheatmasstransfer.2018.07.078
[30]	Jiang, P.-X., Zhao, C.-R. and Liu, B. (2012) Flow and Heat Transfer Characteristics of r22 and Ethanol at Supercritical Pressures. The Journal of Supercritical Fluids, 70, 75-89. https://doi.org/10.1016/j.supflu.2012.06.011
[31]	Xiong, Z., Gu, H., Zhang, S. and Chen, J. (2017) Effect of Intermittent Heating on the Heat Transfer Performance of Supercritical R134a Flowing in a Pipe. Experimental Thermal and Fluid Science, 88, 434-443. https://doi.org/10.1016/j.expthermflusci.2017.06.014
[32]	Liu, X., Xu, X., Liu, C., Ye, J., Li, H., Bai, W. and Dang, C. (2017) Numerical Study of the Effect of Buoyancy Force and Centrifugal Force on Heat Transfer Characteristics of Supercritical CO2 in Helically Coiled Tube at Various Inclination Angles. Applied Thermal Engineering, 116, 500-515. https://doi.org/10.1016/j.applthermaleng.2017.01.103
[33]	Zhang, S., Xu, X., Liu, C., Liu, X., Ru, Z. and Dang, C. (2020) Experimental and Numerical Comparison of the Heat Transfer Behaviors and Buoyancy Effects of Supercritical CO2 in Various Heating Tubes. International Journal of Heat and Mass Transfer, 149, Article ID: 119074. https://doi.org/10.1016/j.ijheatmasstransfer.2019.119074
[34]	Wang, K.-Z., Xu, X.-X., Liu, C., Bai, W.-J. and Dang, C.-B. (2017) Experimental and Numerical Investigation on Heat Transfer Characteristics of Supercritical CO2 in the Cooled Helically Coiled Tube. International Journal of Heat and Mass Transfer, 108, 1645-1655. https://doi.org/10.1016/j.ijheatmasstransfer.2017.01.004
[35]	Xiang, M., Guo, J., Huai, X. and Cui, X. (2017) Thermal Analysis of Supercritical Pressure CO2 in Horizontal Tubes under Cooling Condition. The Journal of Supercritical Fluids, 130, 389-398. https://doi.org/10.1016/j.supflu.2017.04.009
[36]	Wang, K., Xu, X., Wu, Y., Liu, C. and Dang, C. (2015) Numerical Investigation on Heat Transfer of Supercritical CO2 in Heated Helically Coiled Tubes. The Journal of Supercritical Fluids, 99, 112-120. https://doi.org/10.1016/j.supflu.2015.02.001
[37]	Palko, D. and Anglart, H. (2008) Theoretical and Numerical Study of Heat Transfer Deterioration in High Performance Light Water Reactor. Science and Technology of Nuclear Installations, 2008, Article ID: 405072. https://doi.org/10.1155/2008/405072
[38]	Liu, L., Xiao, Z., Yan, X., Zeng, X. and Huang, Y. (2013) Heat Transfer Deterioration to Supercritical Water in Circular Tube and Annular Channel. Nuclear Engineering and Design, 255, 97-104. https://doi.org/10.1016/j.nucengdes.2012.09.025
[39]	Wen, Q. and Gu, H. (2011) Numerical Investigation of Acceleration Effect on Heat Transfer Deterioration Phenomenon in Supercritical Water. Progress in Nuclear Energy, 53, 480-486. https://doi.org/10.1016/j.pnucene.2011.02.012
[40]	Podila, K. and Rao, Y. (2015) CFD Analysis of Flow and Heat Transfer in Canadian Supercritical Water Reactor Bundle. Annals of Nuclear Energy, 75, 1-10. https://doi.org/10.1016/j.anucene.2014.07.039
[41]	Zhao, H., Li, X. and Wu, X. (2017) Numerical Investigation of Supercritical Water Turbulent Flow and Heat Transfer Characteristics in Vertical Helical Tubes. The Journal of Supercritical Fluids, 127, 48-61. https://doi.org/10.1016/j.supflu.2017.03.016
[42]	Jaromin, M. and Anglart, H. (2013) A Numerical Study of Heat Transfer to Supercritical Water Flowing Upward in Vertical Tubes under Normal and Deteriorated Conditions. Nuclear Engineering and Design, 264, 61-70. https://doi.org/10.1016/j.nucengdes.2012.10.028
[43]	Han, C.-L., Zhang, Y.-N., Yu, H., Lu, Y.-P. and Jiao, B. (2018) Numerical Analysis on Non-Uniform Flow and Heat Transfer of Supercritical Cryogenic Methane in a Heated Horizontal Circular Tube. The Journal of Supercritical Fluids, 138, 82-91. https://doi.org/10.1016/j.supflu.2018.04.007
[44]	Zhang, P., Huang, Y., Shen, B. and Wang, R.Z. (2011) Flow and Heat Transfer Characteristics of Supercritical Nitrogen in a Vertical Mini-Tube. International Journal of Thermal Sciences, 50, 287-295. https://doi.org/10.1016/j.ijthermalsci.2010.06.014
[45]	Xiao, Y., Pan, J. and Gu, H. (2018) Numerical Investigation of Spacer Effects on Heat Transfer of Supercritical Fluid Flow in an Annular Channel. International Journal of Heat and Mass Transfer, 121, 343-353. https://doi.org/10.1016/j.ijheatmasstransfer.2018.01.030
[46]	Yang, Z., Cheng, X., Zheng, X. and Chen, H. (2019) Numerical Investigation on Heat Transfer of the Supercritical Fluid Upward in Vertical Tube with Constant Wall Temperature. International Journal of Heat and Mass Transfer, 128, 875-884. https://doi.org/10.1016/j.ijheatmasstransfer.2018.09.049
[47]	Kedzierski, M. and Lin, L. (2020) Pool Boiling of R515A, R1234ze(E), and R1233zd(E) on a Reentrant Cavity Surface. International Journal of Heat and Mass Transfer, 161, Article ID: 120252. https://doi.org/10.1016/j.ijheatmasstransfer.2020.120252
[48]	Petersen, M., Pottker, G., Sethi, A. and Yana Motta, S.F. (2018) Refrigerants with Low Environmental Impact for Commercial Refrigeration Systems. International Refrigeration and Air Conditioning Conference. https://docs.lib.purdue.edu/iracc/1882
[49]	Saengsikhiao, P., Taweekun, J., Maliwan, K., Sae-ung, S. and Theppaya, T. (2020) Investigation and Analysis of R463A as an Alternative Refrigerant to R404A with Lower Global Warming Potential. Energies, 13, 1514. https://doi.org/10.3390/en13061514
[50]	Tian, R., Xu, Y., Shi, L., Song, P. and Wei, M. (2020) Mixed Convection Heat Transfer of Supercritical Pressure R1234yf in Horizontal Flow: Comparison Study as Alternative to R134a in Organic Rankine Cycles. Energy, 205, 118061. https://doi.org/10.1016/j.energy.2020.118061
[51]	Longo, G.A., Mancin, S., Righetti, G. and Zilio, C. (2019) R1234yf and R1234ze(E) as Environmentally Friendly Replacements of R134a: Assessing Flow Boiling on an Experimental Basis. International Journal of Refrigeration, 108, 336-346. https://doi.org/10.1016/j.ijrefrig.2019.09.008
[52]	Jemaa, R.B., Mansouri, R., Boukholda, I. and Bellagi, A. (2017) Energy and Exergy Investigation of R1234ze as R134a Replacement in Vapor Compression Chillers. International Journal of Hydrogen Energy, 42, 12877-12887. https://doi.org/10.1016/j.ijhydene.2016.12.010
[53]	Mota-Babiloni, A., Navarro-Esbrí, J., Molés, F., Cervera, á.B., Peris, B. and Verdú, G. (2016) A Review of Refrigerant R1234ze(E) Recent Investigations. Applied Thermal Engineering, 95, 211-222. https://doi.org/10.1016/j.applthermaleng.2015.09.055
[54]	Li, J., Liu, Q., Ge, Z., Duan, Y. and Yang, Z. (2017) Thermodynamic Performance Analyses and Optimization of Subcritical and Transcritical Organic Rankine Cycles Using R1234ze(E) for 100-200°C Heat Sources. Energy Conversion and Management, 149, 140-154. https://doi.org/10.1016/j.enconman.2017.06.060
[55]	Yataganbaba, A., Kilicarslan, A. and Kurtbaş, İ. (2015) Exergy Analysis of R1234yf and R1234ze as R134a Replacements in a Two Evaporator Vapour Compression Refrigeration System. International Journal of Refrigeration, 60, 26-37. https://doi.org/10.1016/j.ijrefrig.2015.08.010
[56]	Janković, Z., Atienza, J.S. and Suárez, J.A.M. (2015) Thermodynamic and Heat Transfer Analyses for R1234yf and R1234ze (E) as Drop-in Replacements for R134a in a Small Power Refrigerating System. Applied Thermal Engineering, 80, 42-54. https://doi.org/10.1016/j.applthermaleng.2015.01.041