A Criterion of Heat Transfer Deterioration for Supercritical Organic Fluids Flowing Upward and Its Heat Transfer Correlation

: The main objective of this study was to develop the supercritical heat transfer correlation applicable for organic ﬂuids when ﬂowing upward in smooth tubes based on the available experimental data. The organic ﬂuids contain R-22, R-134a, R-245fa and Ethanol and the associated heat transfer characteristics were compared with non-organic ﬂuids like water and carbon-dioxide (CO 2 ). It was found that the limit heat ﬂux may result in heat transfer deterioration (HTD) of organic ﬂuid and the corresponding values are much smaller than water or CO 2 . A new criterion to predict the HTD was developed and this criterion yields the best predictive ability against database. It was found that HTD occurs can be well described by the acceleration parameter evaluated at the wall condition rather than at bulk condition. For estimation of the supercritical heat transfer coe ﬃ cient (HTC) for organic ﬂuid, the present study proposes a new correlation with a physically based correction factor, which gives satisfactory predictions against the HTC of supercritical organic ﬂuid. The new correlation can o ﬀ er the smallest average deviation of 0.007 and standard deviation of 0.181 among the existing correlations.


Introduction
The cycling process of an Organic Rankine Cycle (ORC) is similar to Rankine cycle except the working fluid is organic. ORC can harvest low grade energy such as geothermal energy, solar energy, etc. [1]. However, the thermal efficiency for typical ORC is comparatively low. By changing the subcritical cycle to the trans-critical cycle, the efficiency of ORC can be improved because a better temperature can be reached in the evaporator, and higher exergy is attainable [2][3][4]. However, the estimation to heat transfer performance in the supercritical heating process is much more difficult than the subcritical heating process due to the immense variations of physical properties in the critical/pseudo-critical regions [3]. Figure 1 shows the enormous variations of properties for R-22 [5] at a certain temperature, called the pseudo-critical point. At a pseudo-critical point at a certain specified pressure above the critical point, the specific heat capacity may achieve to the maximum, and the density, viscosity, and thermal conductivity decrease remarkably. Yet, the variations become even more pronounced when the pressure is close to the critical point.  [5], where ρ is density, Cp is specific heat capacity, μ is viscosity, and k is thermal conductivity, respectively.
The gigantic change of physical properties also results in huge change in heat transfer coefficients (HTC). Figure 2 shows the comparison of measured HTC against Dittus-Boelter correlation [5] for R-22. As depicted in the Figure, the correlation is applicable only at a lower heat flux (10 kW/m 2 ) or at a temperature that is some distance away from the pseudo-critical point. A further rise in heat flux to 20 kW/m 2 causes a significant over-prediction of the HTC, and it becomes more conspicuous with the heat flux. This phenomenon of reduction of HTC is called heat transfer deterioration (HTD) [6,7] and the lowest heat flux that leads to HTD is called limit heat flux (LHF) [5,8]. The phenomenon occurs due to the rapid change of density, which induces significant effects regarding buoyancy and acceleration [9][10][11]. In this regard, it is imperative to have an accurate correlation to predict the occurrence of HTD for supercritical fluid.  [5], where ρ is density, Cp is specific heat capacity, µ is viscosity, and k is thermal conductivity, respectively.
The gigantic change of physical properties also results in huge change in heat transfer coefficients (HTC). Figure 2 shows the comparison of measured HTC against Dittus-Boelter correlation [5] for R-22. As depicted in the Figure, the correlation is applicable only at a lower heat flux (10 kW/m 2 ) or at a temperature that is some distance away from the pseudo-critical point. A further rise in heat flux to 20 kW/m 2 causes a significant over-prediction of the HTC, and it becomes more conspicuous with the heat flux. This phenomenon of reduction of HTC is called heat transfer deterioration (HTD) [6,7] and the lowest heat flux that leads to HTD is called limit heat flux (LHF) [5,8]. The phenomenon occurs due to the rapid change of density, which induces significant effects regarding buoyancy and acceleration [9][10][11]. In this regard, it is imperative to have an accurate correlation to predict the occurrence of HTD for supercritical fluid.  [5], where ρ is density, Cp is specific heat capacity, μ is viscosity, and k is thermal conductivity, respectively.
The gigantic change of physical properties also results in huge change in heat transfer coefficients (HTC). Figure 2 shows the comparison of measured HTC against Dittus-Boelter correlation [5] for R-22. As depicted in the Figure, the correlation is applicable only at a lower heat flux (10 kW/m 2 ) or at a temperature that is some distance away from the pseudo-critical point. A further rise in heat flux to 20 kW/m 2 causes a significant over-prediction of the HTC, and it becomes more conspicuous with the heat flux. This phenomenon of reduction of HTC is called heat transfer deterioration (HTD) [6,7] and the lowest heat flux that leads to HTD is called limit heat flux (LHF) [5,8]. The phenomenon occurs due to the rapid change of density, which induces significant effects regarding buoyancy and acceleration [9][10][11]. In this regard, it is imperative to have an accurate correlation to predict the occurrence of HTD for supercritical fluid. There have been many available correlations regarding HTD for inorganic fluids such as water and CO 2 and can be found in typical review articles (e.g., Jackson [11], Pioro et al. [12], and Yoo et al. [13]). However, there is few research concerning the heat transfer performance of organic fluids. The investigation for organic fluids in smooth tubes are listed on Table 1 [5,[14][15][16][17]. Most investigations applied the correlations that were developed from water or carbon dioxide to predict the performance of organic fluids. Yamashita et al. [5] commented that the Watts and Chou [18] correlation developed from water gives the best agreement against the experimental data for R-22. Zhang et al. [16] suggested that the correlation proposed by Jackson [19] shows the best agreement against the experimental data for R-134a. He et al. [17] recommended two correlations that were originally developed from water proposed by Yamagata et al. [8] and Jackson [19] upon the R-245fa experimental data. On the other hand, the existing correlations developed for water or carbon dioxide all fail to predict the experimental data regarding the onset of HTD [16,17]. Several studies have developed new empirical correlations for organic fluids. For example, Kang and Chang [14] proposed a modified correlation based on the Jackson and Fewster correlation [20], and reported a mean deviation of 20% in predicting the Nusselt number for 94.6% of their R-134a experimental data. Zhang et al. [16] proposed a new simple modified correlation Dittus-Boelter correlation [21] by adding an acceleration parameter. The correlation is in good agreement with their R-134a test data. Notice that the foregoing correlations are only applicable for their test data only and no general heat transfer correlation for organic fluids yet. Table 1. Investigations for performance of the heat transfer coefficient (HTC) and heat transfer deterioration (HTD) for organic fluids in smooth tube at supercritical region [5,[14][15][16][17].

Reference
Organic Fluid  [16] R-134a Yes Jackson [19] Yes No He et al. [17] R-245fa Yes Yamagata et al. [8], Jackson [19] No No There were also some review papers regarding the heat transfer performance and the occurrence of HTD for supercritical water and carbon dioxide in the smooth tube [7,22,23]. Figure 3 compares the LHF against the mass flux of water, carbon dioxide, and the organic fluids investigated from the existing literature [5,[14][15][16][17]23]. As shown in Figure 3, water contains the highest LHF among these fluids with the same mass flux, followed by the carbon dioxide, and lastly by the organic fluid. Comparing the criteria with the data of LHF, Yamagata's criterion [8], which is developed by water, shows good agreement only with water. Kim's criterion [24], which is developed based on carbon dioxide, is also only applicable for carbon dioxide. As can be clearly seen, both criteria over-predict the LHF of organic fluids, especially Yamagata's criterion.
Moreover, the quantitative comparison with experimental data of LHF and the existing criteria in terms of polynomial form against the mass flux is listed in Table 2. The average deviation (AD) and standard deviation (SD) are defined [10] below: The criteria developed by water and carbon dioxide were discussed in the review by Huang et al. [23]. The result shows that the criteria can provide good predictions only for its database, extrapolations to other working fluids are normally futile. The criteria developed for water (e.g., Yin et al. [25], Yamagata et al. [8], Styrikovich et al. [26], and Mokry et al. [27]) could provide an absolute average deviation of below 0.6 and a standard deviation of less than 0.3. The criteria developed for carbon dioxide such as Kim et al. [24] could provide an absolute average deviation less than 0.3 with a standard deviation being lower than 0.6. However, it also reveals that these criteria over-predict the data of organic fluids, especially for those criteria developed for water. The average deviations are at least 0.5. Based on the aforementioned reviews, it was concluded that no appropriate criteria of LHF is applicable yet. Moreover, the quantitative comparison with experimental data of LHF and the existing criteria in terms of polynomial form against the mass flux is listed in Table 2. The average deviation (AD) and standard deviation (SD) are defined [10] below: The criteria developed by water and carbon dioxide were discussed in the review by Huang et al. [23]. The result shows that the criteria can provide good predictions only for its database, extrapolations to other working fluids are normally futile. The criteria developed for water (e.g., Yin et al. [25], Yamagata et al. [8], Styrikovich et al. [26], and Mokry et al. [27]) could provide an absolute average deviation of below 0.6 and a standard deviation of less than 0.3. The criteria developed for carbon dioxide such as Kim et al. [24] could provide an absolute average deviation less than 0.3 with a standard deviation being lower than 0.6. However, it also reveals that these criteria over-predict the data of organic fluids, especially for those criteria developed for water. The average deviations are at least 0.5. Based on the aforementioned reviews, it was concluded that no appropriate criteria of LHF is applicable yet.   Table 2. The criteria and predicted the average deviation (AD) and standard deviation (SD) of the LHF of water and carbon dioxide [23]. In view of the foregoing review, despite the fact that there were many studies concerning the predictions about the limit heat flux and the occurrence of the HTD in the literature, there were simply no general criteria and correlations for organic fluids. Hence, the objective of this study was to tailor and bridge the gap in developing the associated criteria and correlations for engineering elaborations. By collecting and analyzing exiting experimental data in the literature, the main objective of this study was to propose a new criterion to describe the heat transfer deterioration of supercritical organic fluids. Yet, a corresponding heat transfer correlation based on the experimental data was also developed that will be shown subsequently to offer a superior predictive ability to the existing correlations.

Data Analysis
The related studies [5,[14][15][16][17] in association with organic fluids are tabulated in Table 3 and all the data points are included as the database for developing the criteria and correlations. Yamashita et al. [5] carried out the experiment for supercritical R-22 flowing upward in the vertical tube. The hydraulic diameter is 4.4 mm with the heated length of 2 m. Their range of operating condition included the pressure of 5.5 MPa, the mass flux of 400-2000 kg/m 2 ·s, and wall heat flux of 10-170 kW/m 2 . Kang and Chang [14] carried out the experiment for supercritical R-134a flowing upward in the vertical tube. The diameter of the test tube is 9.4 mm with a heated length of 2 m. The operating conditions were with the inlet pressure of 4.1, 4.3, and 4.5 MPa, respectively, while the wall heat flux is from 10 to 160 kW/m 2 ·s and the mass flux ranges from 600 to 2000 kg/m 2 ·s. Jiang et al. [15] conducted an experiment of supercritical R-22 and ethanol. The diameter of the tube was 1.004 mm. Their wall heat flux ranged from 110 to 1800 kW/m 2 , while the fluid inlet Reynold number ranged from 3500 to 24000, and the pressure spanned from 5.5 MPa to 10 MPa. Zhang et al. [16] performed the experiment for the supercritical R-134a flowing upward in the vertical tube. The hydraulic diameter was 7.6 mm with the heated length of 2. Prior studies addressed some important factors on the thermofluids characteristics such as fluid flow condition, flow direction, geometry, etc. Yet, it was found that the buoyancy and acceleration effects impose severe effect on the heat transfer performance especially for the upward flowing conditions. Hence, this study stresses only for smooth tubes subject to upward flow conditions. To sum up, the fluid investigated includes R-22, R-134a, R-245fa, and ethanol. The total number of data used to develop the correlation is 4260. The diameter ranges from 1 to 10 mm, and the corresponding mass flux and heat flux ranges from 400-4000 kg/m 2 ·s, and 10-1800 kW/m 2 , respectively. Table 3. The experimental data from the literatures [5,[14][15][16][17].

Data Reduction
Since part of the database only discloses the bulk enthalpy and wall temperature, the corresponding bulk temperature and corresponding physical properties were evaluated with the help of REFPROP ver. 8 [28] from the prescribed pressures and the bulk enthalpies. To further investigate the heat transfer performance, the heat transfer coefficient (HTC) was used and analyzed in the following: where q represents the heat flux and T w and T b represent wall and bulk temperature, respectively. To facilitate detailed comparisons amid the proposed and some existing correlations, the HTC is further termed into dimensionless Nusselt number: where D represents the diameter, and k is thermal conductivity. The collected experimental data is further analyzed through Matlab 2018b. Quantitative comparisons with experimental data of Nusselt number and the existing correlation are termed as the average deviation AD and the standard deviation SD [10] in the following: Figure 4 shows the Nusselt number versus the enthalpy at different conditions from Yamashita et al. [5]. As addressed in prior studies, the effect of heat flux occurs only when buoyancy and acceleration effects is in control and this influence is accentuated nearby pseudo-critical regions. For a mass flowrate of 400 kg/m 2 ·s, the Nusselt number increases with enthalpy and then level-off appreciably when it reaches the vicinity of pseudo-critical enthalpy. The peak of Nusselt number becomes less profound when raising the heat flux. However, when heat flux achieves to a threshold value such as 30 kW/m 2 (for G = 400 kg/m 2 ·s), the Nusselt number suddenly drops before the bulk temperature reaches the pseudo-critical temperature. In essence, HTD occurs accordingly.  At a mass flow rate of 1000 kg/m 2 •s, the trend of HTD was similar to low mass flow rate of 400 kg/m 2 •s. However, the Nusselt number of a mass flow rate of 1000 kg/m 2 •s was higher than the one of 400 kg/m 2 •s at lower enthalpy region for its much higher Reynolds number. Interestingly, at the heat flux at 30 kW/ m 2 , HTD did not happen at a mass flow rate of 1000 kg/ m 2 •s as that of G = 400 kg/m 2 •s. Instead, HTD occurred at a higher heat flux such as 90 kW/m 2 . It reveals the lowest heat flux leading the HTD, which is called limit heat flux (LHF) [5,8], which depends on mass flux.

The Criterion of the Lowest Heat Flux for Heat Transfer Deterioration
As aforementioned, there were no appropriate criteria of LHF using polynomial forms for organic fluid. However, Cheng et al. [10] proposed a different polynomial form as the criterion of LHF (Equation (7)). The authors claimed that the property of fluid is included in the criterion: At a mass flow rate of 1000 kg/m 2 ·s, the trend of HTD was similar to low mass flow rate of 400 kg/m 2 ·s. However, the Nusselt number of a mass flow rate of 1000 kg/m 2 ·s was higher than the one of 400 kg/m 2 ·s at lower enthalpy region for its much higher Reynolds number. Interestingly, at the heat flux at 30 kW/ m 2 , HTD did not happen at a mass flow rate of 1000 kg/ m 2 ·s as that of G = 400 kg/m 2 ·s. Instead, HTD occurred at a higher heat flux such as 90 kW/m 2 . It reveals the lowest heat flux leading the HTD, which is called limit heat flux (LHF) [5,8], which depends on mass flux.
As aforementioned, there were no appropriate criteria of LHF using polynomial forms for organic fluid. However, Cheng et al. [10] proposed a different polynomial form as the criterion of LHF (Equation (7)). The authors claimed that the property of fluid is included in the criterion: where Cp PC and β PC is specific heat capacity and thermal expansion coefficient at pseudo-critical condition, respectively. Figure 5 shows a good agreement between the experimental LHF of water and carbon dioxide and the criterion of Cheng et al. [10]. Moreover, the criterion of Cheng et al. [10] contained the effect of pressure in comparison with Yamagata's criterion [8]. Although the criterion by Cheng et al. [10] cannot predict the HTD of the organic fluids well, it revealed that the ratio of β to Cp may affect the occurrence of HTD. The ratios are listed in Table 4, and it was found that the ratios for the organic fluids except ethanol ranging from 13 × 10 −6 to 16 × 10 −6 kg/J. The values were much higher than 1.79 × 10 −6 kg/J of water and 8.80 × 10 −6 kg/J of carbon dioxide at a near critical pressure. The higher ratios lead to the occurrence of HTD at a quite low heat flux due to the higher level of thermal expansion or the higher variety of density for given a heat flux or mass flux. The magnitude of the ratios of β to Cp can explain the LHF of organic fluids that is lower than water and carbon dioxide at a given mass flux.
The correlation is valid for ratios ranging from 13 × 10 -6 to 16 × 10 -6 kg/J and the fluid of R-22, R-134a, and R-245fa. The quantitative comparison amid the experimental data and the criteria proposed by Cheng et al. [10] and the proposed criterion is shown as Table 5. It also reveals a good agreement between the experimental LHF of water and carbon dioxide and the criterion of Cheng et al [10]. However, their correlation showed a deviation of 0.650 for organic fluids, which is much higher than the standard deviation of 0.336 upon the proposed new criterion. Moreover, the absolute average deviation of 0.232 for organic fluids is lower than the deviation provided by the existing criteria form Table 2. In conclusion, the proposed new criterion of LHF for organic fluids shows the best predictive ability.

The onset of Laminarization and HTD
The occurrence of HTD is related to acceleration [9,11]. When the effect of acceleration is increased, the turbulence might be suppressed and the heat transfer performance decreases. The dimensionless form of acceleration factor is proposed by McEligot et al. [9] as Kv (acceleration factor): McEligot et al. [9] suggested that the flow changes from turbulent to laminar when Kv is higher than 3 × 10 −6 . However, Jiang et al. [15] indicated that Kv for R-22 is far less than the threshold value proposed by McEilgot et al. [9]. On the other hand, Cheng et al. [10] proposed a simpler parameter π A to describe the acceleration effect, which is defined as: Obviously, Kv is more complex than π A . Hence, using π A is easier to implement from an engineering perspective and will be adopted in the following developments. Figures 6-8 depict the corresponding Nusselt number and acceleration factor π A at both bulk and wall condition versus enthalpy for R-22, R-134a, and R-245fa. It was found that the peak of π A is increased with a rise of heat flux. When compared with Nusselt number and π A at bulk condition, the Nusselt number will decrease first, followed by an increase against the enthalpy when π A passes through a certain value against the enthalpy. This is the basic cause for deterioration of heat transfer performance. Manipulation and combination of the criterion of LHF by Equation (8) and π A yields the corresponding threshold value: The value shown as dash line in Figure 6 indicates that HTD does occur in the condition when the peak of π A is beyond the threshold. However, it cannot predict the location where HTD occurs in association with π A at bulk condition. Hence, π A at the wall condition should be also taken in account.
By comparing Nusselt number and π A at wall condition, it is found that HTD occurs in the vicinity even after π A at wall condition approaches the threshold value. This correlation can predict the onset of HTD due to laminarization [9,11].       Table 6 lists some correlations that are employed for further comparison against the database from the literature. Table 6. Some correlations used for comparisons against the experimental data subject to supercritical conditions.    Table 7 shows the comparisons between the experimental data and the existing correlations. The average deviation and standard deviation are defined in Eqs. (5)(6) [10], which are used to compare the experimental data. For the data of R-22, Kang and Chang's correlation [14] provides the smallest average deviation of 0.068 and 0.014 against the data of Yamashita et al. [5] and Jiang et al. [15], respectively. The correlations of Watts and Chou [18] and Yamagata et al. [8] provide the smallest standard deviation of 0.169 and 0.462 against the database of Yamashita et al. [5] and Jiang et al. [15], respectively.

Dittus-Boelter [21]
For the data of R-134a, Kang and Chang's correlation [14] also provides the smallest average deviation of 0.117 and −0.040 against the data of both Kang and Chang [14] and Zhang et al. [16], respectively. The correlation of Yamagata et al. [8] provides the smallest standard deviation 0.150 and 0.099 against the database from Kang and Chang [14] and Zhang et al. [16], respectively.
For the data of R-245fa, Yamagata's correlation [8] provides the smallest average deviation of 0.061, and Watts and Chou's correlation [18] provide the smallest standard deviation.
For the data of ethanol from Jiang et al. [15], Zhang's correlation [16] offers the smallest average value, and Yamagata's correlation [8] provides the smallest standard deviation. For all data, Watts and Chou's correlation [18] can provide the smallest average deviation (−0.017) and standard deviation (0.291) among these correlations.
Moreover, the performance of empirical correlation for heat transfer to supercritical organic fluid is listed Table 8. It reveals that Kang and Chang correlation [14] provides the predictive ability against the R-22 and R-134a data. Watts and Chou's correlation [18] provides the best agreement against the R-245fa data. Zhang's correlation [16] offers the best agreement against the ethanol data. Obviously, the best agreement is obtained if the correlation has the smallest either average or standard deviation. It implies that there are no existing correlations that are able to provide satisfactory agreement against all organic data. In fact, Watts and Chou's correlation [18] shows the best overall predictive ability in which 73% of 4260 data falls within the ±30% span.

The Modified Correlation and Corresponding Correction Factor
From the above discussions, Watts and Chou's correlation [18] provides the best overall predictive ability. Hence, efforts are made in modifying the correlation by introducing a rational dimensionless parameter, Gr * , into the correlation. Figure 9 shows the correction factor (CF) proposed by Watts and Chou [18], π A at wall condition and Gr * against enthalpy for R-22 [5]. Note that CF is defined as: It was found that the trend of Gr * is more similar to Nusselt number as compared to π A at the wall condition. Gr * is more appropriate as the correction factor for heat transfer. Moreover, the threshold proposed Watt and Chou [18] is valid for these organic fluids at lower mass flux. However, it is not applicable at a higher mass flux. Figure 9 shows Gr * does not reach the threshold value even when HTD happens. The threshold may be revised to the baseline relative to the boundary condition. It is observed that the baseline of Gr * increase with the heat flux and decrease with the mass flux, implying different resistance that limits HTD happening due to buoyancy force at different conditions. The change of the correction factor, however, almost correspond to the change of Gr * whatever the heat and mass flux are. The baseline can be estimated using the local or bulk condition of the working fluids: where Gr inlet can be written as: where ρ is the density calculated by integration-averaging with temperature in between the wall and bulk conditions. In the subcritical region, the change of density with the temperature is quite small and ρ can be assumed as the arithmetic average density of the wall and bulk. The ratio of the change of density may be written below: The temperature difference between wall and bulk condition can be written as: In subcritical region, the buoyancy effect is much smaller than the inertia effect. The forced convection dominates. Hence, the Dittus-Boelter correlation can be applied to estimate the heat transfer coefficient: According to Equations (13)(14)(15)(16)(17), Gr * base can be written as: (18) Figure 10 shows the effect of these dimensionless parameters such as π A,b , π A,w , Gr * and Gr * /Gr * base . The analysis includes six databases that are considered in this study. It shows the deviations of database pertaining to Gr * /Gr * base is much better than other dimensionless parameters. As a consequence, Gr * /Gr * base is selected as the correction factor for developing the new correlation. Eventually, the derived new correlation is as follows: The correlation takes the basic form from Watts and Chou [18]. Yet Gr * is replaced by Gr * /Gr * base as the correction factor and the fitting constant is slightly higher changing from 0.021 to 0.0219. Figure 11 and Table 9 shows the comparison between the predictive Nusselt number and the experimental data all literature. From the table, the proposed correlation offers an average deviation of 0.007 and a standard deviation of 0.181 for all data (4260 in total). With the proposed correlation, Table 10 also indicates that 90% of the R-22 data, 94% of the R-134a data, and 95% of the R-245fa are within the range of ±30%. In essence, the proposed correlation yields the best predictive ability for the organic fluids except ethanol.
The correlation takes the basic form from Watts and Chou [18]. Yet * Gr is replaced by * * / base Gr Gr as the correction factor and the fitting constant is slightly higher changing from 0.021 to 0.0219.  Figure 11 and Table 9 shows the comparison between the predictive Nusselt number and the experimental data all literature. From the table, the proposed correlation offers an average deviation of 0.007 and a standard deviation of 0.181 for all data (4260 in total). With the proposed correlation, Table 10 also indicates that 90% of the R-22 data, 94% of the R-134a data, and 95% of the R-245fa are within the range of ±30%. In essence, the proposed correlation yields the best predictive ability for the organic fluids except ethanol. R-22 [15] R-134a [14] R-134a [16] R-245fa [17] Ethanol [15] (b)   Figure 11. Comparison of Nusselt numbers between the proposed correlation and all experimental data. Table 9. Comparisons between the experimental data and the proposed correlation.
For ethanol, the study also provides another new correlation to estimate the heat transfer and show as: Figure 12 shows the new correlation for ethanol that is in good agreement with the experimental data from Jiang et al. [15]. Additionally, 96% of the ethanol data are within ±20% of predictions.

Conclusions
This study investigated the heat performance of organic fluids for smooth tubes with fluid flowing upward. The organic fluids include R-22, R-134a, R-245fa, and ethanol. It was found that the limit heat flux for organic fluid near the critical point is much smaller than that of the water and carbon dioxide due to the difference in the ratio of the thermal expansion coefficient and the specific heat capacity. The study presents a new criterion of limit heat flux (LHF) for organic fluids and the proposed criterion of LHF for organic fluid is superior to the existing criteria. The correlation shows an average deviation of -0.232 and a standard deviation of 0.336 for organic fluid.
In the investigation about the onset of the heat transfer deterioration, acceleration parameter

Conclusions
This study investigated the heat performance of organic fluids for smooth tubes with fluid flowing upward. The organic fluids include R-22, R-134a, R-245fa, and ethanol. It was found that the limit heat flux for organic fluid near the critical point is much smaller than that of the water and carbon dioxide due to the difference in the ratio of the thermal expansion coefficient and the specific heat capacity. The study presents a new criterion of limit heat flux (LHF) for organic fluids and the proposed criterion of LHF for organic fluid is superior to the existing criteria. The correlation shows an average deviation of -0.232 and a standard deviation of 0.336 for organic fluid.
In the investigation about the onset of the heat transfer deterioration, acceleration parameter imposes a significant effect on the onset of heat transfer deterioration (HTD) due to laminarization. The acceleration parameter, π A , imposed appreciable influence on the heat transfer deterioration. For accurate assessment of HTD, the estimation for threshold π A should be at the wall condition rather than at the bulk condition. The threshold π A is strongly related to the limit heat flux and property at the pseudo-critical point.
Based on the evaluations of the existing correlations against existing data, it was found that most correlations are only applicable to their own datasets. Through detailed analysis of the prior data and correlations, the non-dimensional parameter Gr * offers better predictive ability of the heat transfer coefficient for organic fluids. The correlation takes the basic form from Watts and Chou [18]. Yet, Gr * is replaced by Gr * /Gr * base as the correction factor. The new correlation can provide the smallest average deviation of 0.007 and standard deviation of 0.181 against all data among these common correlations for supercritical conditions. The 90% of 4260 data are within the range of ±30% of the presented correlation. It is the best among other correlations. On the other hand, the study also presents another new correlation for ethanol due to the worse agreement of previously presented correlation, and 96% of the 262 data for ethanol are within the range of ±20% of the correlation.

Conflicts of Interest:
The authors declare no conflict of interest.