1. Introduction
In recognition environmental challenges such as ozone depletion and global warming, there is a growing interest among manufacturers, end-users, and scientific researchers regarding the utilization of natural operating fluids in air conditioning, heat pump, and refrigeration systems [
1].
However, the use of natural refrigerants has become increasingly complex, involving countless experimental research on heat transfers, pressure drops, and fluid flow patterns. Natural refrigerant research has become extremely active, with the goal of establishing more precise design methodologies and more energy-efficient cycles using natural refrigerants.
Natural refrigerants such as ethanol, hydrocarbons such as propane (R290), butane (R600)/isobutene (R600a), and others have been studied over the last decade as a replacement for CFCs, HCFCs, and HFCs in refrigeration, air conditioning, and heat pump systems [
1,
2,
3].The ozone depletion potentials (ODPs) of these natural refrigerants are zero, and most of them have close-to-zero global warming potentials (GWPs) in comparison to CFCs and HCFCs [
1,
2,
3].
Hydrocarbons had earlier been employed as refrigerants. Because of increased concern for the environment, their use has recently been reassessed. Hydrocarbon (HC) refrigerants are already used in small freezers in some European countries. Propane, isobutane, n-butane, perfluorocyclobutane, cyclopropane, propylene, and other common hydrocarbons are under investigation [
4]. In low-temperature applications, R290 is thought to be a suitable alternative for R22 and R502. In household refrigerators, R290 mixtures with R600a or R600 are proposed as R12 substitutions.
For those reasons, hydrocarbons such as Propane (R290) are commonly utilized as a refrigerant in a large cascade refrigeration system in ethylene manufacturing facilities. Germany has recently considered using propane and propane/n-butane mixes as a replacement for CFCs, HCFCs, and HFCs [
5]. Furthermore, iso-butane is now used in the majority of new residential refrigerators in Germany. As a result, the usage of hydrocarbons as refrigerants is increasing due to their effectiveness.
The only issue is that hydrocarbon refrigerants are highly flammable. However, most refrigerators use a hermetic compressor, which seals the entire system. Furthermore, because hydrocarbons have a higher specific volume of liquid when compared to traditional refrigerants, the refrigerant charge is reduced.
The utilization of natural refrigerants has become progressively intricate, necessitating numerous experimental studies on heat transfer of these fluids. Along with the goal of developing more precise design techniques and energy-efficient cycles that use natural refrigerants, research in this area has become highly active. Nevertheless, many aspects of the present correlations require enhancement to attain the aimed level of accuracy for the design of refrigerant evaporators, as observed by Thome (1996) [
6], from a predictive perspective.
Therefore, the aim of the present study is to present a comprehensive review on two-phase flow boiling heat transfer coefficient of ethanol and hydrocarbons (R290, R600, and R600a) refrigerants. Firstly, the experimental work done in literature was discussed. Then, a description of the available flow boiling heat transfer characteristic of the mentioned fluids and the two-phase heat transfer coefficient prediction correlations was assessed. Since each correlation is derived based on its own data, fluids, geometry, and operating conditions. Therefore, there are no specific prediction methods for ethanol and hydrocarbon. Furthermore, studying the operating fluids and using their data to evaluate the performance of the prediction correlations will be definitely an efficient way that might help in better understanding the complexity of natural refrigerants in the field of boiling two-phase flow. Leading to developing better designs and optimization methods for manufacturing and heat transfer equipment.
3. Assessment of Previous Correlations
Many researchers offered a number of correlations to predict the heat transfer coefficient of flow boiling. The applicability of these prediction approaches, on the other hand, is often dependent on a specific database and is derived based on their operating conditions such as saturation temperature, pressure, heat flux, and mass flux ranges. The correlations are also dependent on their own experimental geometry, which includes tube orientation (e.g., horizontal, vertical, inclined), tube diameter, etc. In general, these correlations can be classified into two main categories: (A) the superposition model, which combines nucleate boiling and two-phase forced convection mechanisms, and (B) functions of the most significant dimensionless groups.
Review of Flow Boiling Heat Transfer Coefficient Correlations
Chen Correlation [
30] is considered to be a leading method for calculating heat transfer coefficients in convective boiling Chen employed a superposition model (Category A) summarized in
Table 3, which combines nucleate pool boiling with convective heat transmission. Nevertheless, convection is responsible for the suppression of nucleate boiling. Chen declared that when flow rate increases, the temperature gradient near the wall decreases, reducing the temperature difference between the tube inner wall and bubbles forming outside the wall. Numerous studies of flow boiling heat transfer correlations based on the Chen superposition model have been proposed to date, including the Liu-Winterton [
22], Gungor-Winterton [
19], and Saitoh et al. [
31], Jung et al. [
32], Bennett and Chen [
33] and Choi et al. [
34] Correlations.
Bennett and Chen [
30] employed Chen superposition model for over 1,000 data sets for forced convective boiling of pure water, ethylene glycol, and aqueous ethylene glycol mixes. The majority of these measurements were made in the annular flow regime. Their study's aims were to collect forced convective boiling heat transfer data for a liquid with a relatively high Prandtl number and for aqueous mixtures, and to develop a correlation based on these data sets that would be applicable to other pure com-ponents and binary mixtures. The authors expected that their proposed correlation would be valid for boiling of binary mixtures anticipated to be in the annular and annular dispersed flow regimes. It should also be applicable for boiling liquids with Prandtl values between unity and ten under saturated circumstances. Furthermore, because it reduces to the normal Chen correlation for pure liquids with Prandtl values close to unity, it should be applicable for the same range of flow conditions as the Chen equation.
Lazarek and Black [
38] introduced a new method for constructing correlation by fitting experimental data which consist of 728 data points and R113 was used a working fluid as a function of the most important dimensionless groups (Category B) as shown in
Table 4. The authors used mini and micro-channels with 3.1 mm diameter and range G = 125 - 750 kg/m
2s, q = 14 – 380 kW/m
2 and P = 1.3 - 4.1 bar. The liquid only Reynolds number and the boiling number were utilized by the authors. The influence of mass flux on the proposed correlation was negligible.
In contrast to Chen's superposition concept, Kenning and Cooper [
44] suggested that the heat transfer mechanism in narrow tubes could be more like an enhancement of single-phase forced convection heat transfer. The expression of enhanced convective heat transfer was then proposed. To express the effect of two-phase flow, they used the LockhatMartinelli parameter ‘Xtt’. The correlation of Kenning and Cooper [
44] takes the form of the product of single-phase forced convection and the function of the Lockhart-Martinelli parameter. However, Kenning and Cooper's experimental correlation data were still based on water as working fluid in a macro-scale tube with range of diameter of 9.6 - 14.1mm, x = 0 - 1, p = 0.16 - 0.6 MPa and G = 123 - 630 kg/m2s
Saitoh et al. [
31] developed a correlation for the boiling heat transfer of R-134a in horizontal tubes including the effect of tube diameter were developed here for both the pre- and post-dry out regions. The authors modified Chen correlation model for flow boiling heat transfer was established, which includes the effect of tube diameter. The Weber number was used to describe the effect of tube diameter on flow boiling heat transfer. For a wide range of tube diameters from 0.51 to 10.92 mm, the correlation corresponded reasonably well with experimental results.
Jung et al. [
32] studied the mixture effect on a horizontal flow boiling heat trans-fer. More than 2000 local heat transfer coefficients are obtained with the azeotropic R12/R152a mixture and compared against the previously measured data with the non-azeotropic R22/Rll4 mixture. He modified the superposition model suggested by chen using only phase equilibrium data to consider mixture effects to develop his correlation. Firstly, they replaced Forster and Zuber’s [
48] nucleate pool boiling correlation shown in Equation (1) by Stephan and Abdelsalam correlation [
49] as shown in Equation (2).
Secondly, they discussed that in a given quality in the partial boiling regime, nu-cleate boiling coefficient is a strong function of heat flux. It is also a function of mass flow rate for a given heat flux. Consequently, it becomes a function of quality, heat flux, and mass flow rate. After Since the derived their own suppression factor ‘S’ which is a strong function of heat flux which is differs from other similar factors such as the other suppression factors by Chen [
30] and by Gungor and Winterton [
19] fac-tors which are independent of heat flux.
Choi et al. [
34] conducted experiments on the convective boiling heat transfer in horizontal mini channels with CO2. The test section is made of stainless-steel tubes with inner diameters of 1.5 and 3.0 mm and range of G= 200 – 600 kg/m2s, x= 0 – 1, Tsat = 0, 10, -5, -10 °C. In their developed correlation, the have used chen superposition model, however they developed their own suppression factor and enhancement factor by including the two-phase frictional multiplier that is based on pressure gradient for liquid only flow according to Chisholm [
50]. They discovered that the heat transfer co-efficient at low vapor qualities (up to 0.2) is independent of mass flux and vapor quality but reliant on heat flux, and the boiling regime was designated as nucleate boiling. In the convective boiling regime, the heat transfer coefficient increases with mass flux and vapor quality. At high vapor quality (0.5-1), a drop in heat transfer coefficient happened more rapidly for higher mass fluxes due to the annular regime occurring earlier at a higher mass flow as shown in Equation (3).
ElFaham and Tang [
51] conducted a review study for the flow boiling heat transfer of pure ethanol in various experimental work. A wide range of parameters has been used in their study, Where Tsat = 4.-86.6 °C, G = 0.33 - 290 kg/m2s, q = 2.8 - 104 kW/m2 and x= 0.11 – 0.91. As well as the range of internal diameter is 5-10 mm. Their findings were Liu-Winterton [
22] and Chen [
30] correlations showed a good prediction among 14 correlations in the study. The authors have used the same approach as Liu-Winterton [
22] correlation to derive their proposed correlation. Furthermore, they modified the suppression factor ‘S’ by introducing a suppression factor multiplier (Ms) to control the performance of suppressing the nucleate boiling contribution. They also classified their dataset into 4 regions using boiling number intervals to identify the value of (Ms) in each region [
35]. Also, they used Cooper’s correlation as shown in Equation (4) and Dittus-Boelter correlation shown in Equation (5) to calculate the nucleate boiling and single-phase heat transfer coefficient.
Warrier et al. (2002) [
43] conducted sub-cooled and saturated nucleate boiling experiments in a multi micro-channel section with 5 rectangular channels with a hydraulic diameter of 0.75 mm. FC-84 is the test fluid. Their saturated boiling studies revealed that for a given Boiling number, the boiling heat transfer coefficient falls with increasing vapor quality. They presented an empirical correlation for saturated nucleate boiling and vapor qualities up to 0.55, with the heat transfer coefficient associated with the liquid heat transfer coefficient, boiling number, and vapor quality, with a maximum divergence with experimental data of 28%.
The flow boiling heat transfer characteristics were investigated by Tran [
45] for R12 as working fluid. The authors used in their experiments small circular and rectangular tubes with diameters of d
in= 2.46 mm, d
h = 2.4 mm respectively. Moreover, a range of vapor quality (x = 0 - 0.94), Mass Flux (G = 44 – 832 kg/m
2s) and Heat flux (q = 7.5 – 129 kW/m
2). He emphasized that nucleate boiling contribution predominates in high-wall superheat conditions, while forced convection is the primary heat transfer mechanism in low-wall superheat conditions.
Wattelet et al. [
37] The authors conducted their experiment for R-134a, MP-39, and R-12 in horizontal smooth 7.04 mm internal diameter copper tube, with testing parameters with ranges of G = 25-100 kg/m
2s, q = 2-10 kW/m
2, x=0.1 - 0.9, T
sat = -15 – 5 °C. Experimental heat transfer coefficients were reported from their experimental work. In addition, an empirical correlation developed for annular flow data using an asymptotic form was modified to account for the decrease in heat transfer due to the wavy-stratified flow pattern in the low mass flux cases.
Yoon et al. [
36] measured heat transfer coefficients and pressure drop during evaporation process of carbon dioxide in a horizontal smooth stainless-steel tube with internal diameter 7.53 mm. Experiments were conducted at Tsat = 4 – 20 °C, q = 12 - 20 kW/m
2 and G = 200 - 530 kg/m
2s. They developed a correlation for carbon dioxide during evaporation was developed by considering the critical quality. They proposed that the correlation can predict the critical quality at which the liquid film breaks down at top of tube. Before the critical quality, the Liu and Winterton correlation is used to predict the heat transfer coefficient. After the critical quality, the Dittus–Boelter equation for vapor flow and the Gungor and Winterton correlation for liquid flow is superposed to predict the heat transfer coefficient of carbon dioxide.
Kew and Cornwell [
39] introduced the confinement number in their model to take into consideration for the effect of bubbles limited by the mini-channel wall. They identified three types of flow patterns in mini-scale tubes: isolated bubble regime, restricted bubble regime, and annular-slug flow. In terms of bubble regime, nucleate boiling is the primary heat transfer mechanism, whereas convection is more important in the confined bubble regime and annular-slug flow. However, the authors presented the concept of confinement number for the correlation, but the final version proposed by Kew and Conwell is an improved correlation based on Lazarek and Black [
38] correlation that takes vapor quality into account. They found that nucleate boiling in pre and post dry out zones is related to mass flux, boiling number, and vapor quality. It was shown that forced convection is the major mechanism of flow boiling heat transfer in small channels under their experimental conditions. They also considered the impact of equilibrium quality.
Yu et al. [
42] used his experimental results to propose a correlation that was the same as Tran.[
45] Water, ethylene glycol, and aqueous mixtures of ethylene glycol was the experimental fluids at high temperatures up to 250 °C, pressure (<345 kPa), G= 50 – 200 kg/m
2s. The proposed correlation for two-phase heat transfer coefficient based on water flow boiling in a small 2.98 mm diameter horizontal channel which is also developed based on boiling results with refrigerants under similar conditions.
Sun and Mishima [
41] used a relatively large database of 2050 data points from 20 published studies to propose a correlation based on the Lazarek and Black [
36] correlation which is a function of the most important dimensionless groups as shown in
Table 5.
The experimental results of Tran [
45] reveal that nucleate boiling is the major process and that the effect of mass flux is insignificant. Furthermore, Tran's formula is totally empirical, with R12 as the working fluid in a circular channel with a hydraulic diameter of 2.46 mm and a tiny rectangular channel with a hydraulic diameter of 2.4 mm.
Agostini and Bontemps [
46] have tested flow boiling of R-134a in multiple mini-channels experimentally. The tube is made up of 11 rectangular channels with a hydraulic diameter of 2.01 mm. They discovered a convective to nucleate boiling regime transition for a wall superheat greater than 3 K and heat flux greater than 14 kW. m2, which was consistent with Tran et al. [
45] findings. They also observed dry-out at a critical quality of 0.4, without considering heat or mass flux.
Hamdar et al. [
40] investigated the boiling heat transfer coefficient and pressure drop of R-152a in a 1-mm square mini-channel under test circumstances of G= 200-600 kg m-2 s-1, q= 10-60 kW/m
2, P= 600 Kpa, and T
sat = 25.2 °C. The Nusselt number and the dimensionless groups defined by Tran et al. [
45] are used to correlate the experimental data matching their results. For the data set, the correlation coefficients are determined using a linear multi-regression technique.