1. Introduction
Heat exchangers have applications in many industries (i.e., chemical, energy, food, pharmaceuticals, aviation, transportation, refrigeration, HVAC, etc.). Optimized heat exchanger designs reduce the carbon footprint, improve energy efficiency, reduce production costs, extend equipment lifespan, improve product quality, lower operating costs, produce more energy, reduce the footprint of the unit, etc. Therefore, studies of newly developed enhanced heat transfer tubes are important in order to be able to design high efficiency heat transfer systems.
Heat transfer performance of a heat exchanger tube is influenced by several parameters (i.e., geometric shape of the tube, surface enhancement structure of the tube, type of refrigerant, mass flow rate, saturation temperature, vapor quality, etc.). Previous condensation studies have been typically performed on smooth tubes; some studies have evaluated enhanced tubes. Longo et al. [
1] performed a comparative study on the condensation heat transfer characteristics of micro fin tubes using refrigerants R410A and R32 tubes with an inner diameter of 4.2 mm. This study investigated the effects of saturation temperature, vapor quality, and mass flow rate on the heat transfer coefficient (HTC) of micro fin tubes. Results indicate that micro fin tubes exhibited the best enhancement (when compared to smooth tubes) for mass flow rates ranging from 150 to 300 kg/m
2·s. Additionally, they determined that when using R410A the condensation HTC and pressure drop were greater than when using R32. Ma et al. [
2] conducted an experimental study that evaluated the condensation heat transfer performance of smooth tubes, micro-fin tubes, and enhanced three-dimensional (3D) tubes (with an inner diameter of 9.52 mm) using R410A and R32 as the working fluids. Under varying mass flow rates and vapor qualities, the results showed that the condensation heat transfer coefficient (when using R410A) in the micro-fin was 2.0~2.2 times the HTC found in smooth tubes, and for the 3D enhanced tubes, it was 1.4~1.5 times greater, while when using R32 it was 1.5~2.0 times (micro fin) and 1.5~1.6 times (3D tube) the HTC found in smooth tubes. By analyzing the flow pattern, they determined that the enhanced tubes promoted the occurrence of intermittent flow and annular flow and this produces an improvement to the condensation heat transfer coefficient (when compared with a smooth tube).
The emergence of hydrophobic surfaces provides a new direction and research area to consider when evaluating the performance of enhanced heat transfer tubes. Alwazzan et al. [
3] studied the condensation heat transfer characteristics of a copper tube composed of alternating sections of hydrophobic and hydrophilic areas (studying tubes with different area ratios). As a result of their experimental evaluation and visual observations, they concluded that the use of alternating surface structures can significantly enhance the condensation heat transfer performance. When surfaces (composed of different width ratios) are compared with hydrophobic (with a uniform stripe width) and hydrophilic surfaces (with uniform stripe width), the heat transfer coefficients (for the surfaces composed of different width ratios) were increased by 480% and 180%, respectively. Additionally, it appears that there is an optimal width ratio that maximizes the enhancement effect. Wan et al. [
4] experimentally studied the heat transfer performance of a composite surface enhanced tube with a hydrophobic inner cavity wall and a hydrophilic outer fin surface. They found that the condensation heat transfer coefficient and boiling heat transfer coefficient of the tube with a hydrophobic inner cavity wall and hydrophobic outer fin surface were 3 times and 1.8 times the HTC of a smooth tube, respectively. They go on to conclude that hydrophilic surface enhanced tubes should be preferred for boiling heat transfer, and hydrophobic surface enhanced tubes should be preferred for condensation heat transfer. Additional work in this area is needed.
Scholars have previously studied correlation formulas that are used to determine condensation heat transfer coefficients in tubes; some studies have proposed empirical formulas that can be used in the design and optimization of heat exchangers. Yun et al. [
5] proposed a new correlation (based on previous research) for heat transfer coefficients and friction factors. They introduced correction factors (C1 and C2 for use with R134a and R410A) that are used in a comprehensive correlation (within a 25% error range) and had prediction accuracies of 88.1% and 73.2%, respectively. Hosseini et al. [
6] collected in-tube condensation heat transfer data (6521 data points) from 40 published works and evaluated eight existing condensation correlations. They developed a new general correlation that can accurately estimate the condensation heat transfer coefficient in both microchannels and conventional channels. Li et al. [
7] studied the condensation heat transfer characteristics of R410A in smooth tubes, herringbone micro fin tubes, dimpled tubes, and herringbone-dimpled composite tubes. They found that the condensation heat transfer coefficient and pressure drop increased with the refrigerant mass flow rate and average vapor quality. An improved Nusselt equation was developed for smooth tubes. Additionally, they introduced the mass flux factor (EUR) and the enhanced tube structure factor (Y) that could be used in order to predict the condensation HTC in the enhanced tubes. Shah [
8] proposed modified correlation equations for condensation heat transfer in microchannels and conventional channels for use with horizontal or downward flows. Shah’s equations have a wide range of applications, covering 51 types of fluids (including water, low-temperature media, refrigerants, chemicals, etc.) for tubes with diameters ranging from 0.08 to 49.0 mm, with a mass flux range from 1.1 to 1400 kg/m
2·s for use with various channel cross-sections (including circular, rectangular, triangular, etc.). Additional conditions considered include (i) single and multiple channels, (ii) annular flow, and (iii) horizontal or vertical downward flow conditions. These modified equations were validated using a database consisting of 8298 data points from 130 sources; the equations had a high prediction accuracy that had an average absolute error of 17.9%.
There are few studies on enhanced tube condensation heat transfer performance correlations that describe the condensation HTC for enhanced tubes. In order to advance the knowledge base for enhanced heat transfer tubes, a modification (for use in enhanced tubes) of an existing smooth tube condensation correlation was undertaken in the present study; the modified correlation will allow enhanced tubes to be used when trying to optimize the design of heat exchangers. This will improve heat transfer efficiency and reduce costs (production and operational). The condensation heat transfer performance of smooth (ST), hydrophobic (HYD), herringbone micro fin (HB), and hydrophobic/herringbone (HYD/HB) tubes was evaluated for various mass flow rates. Several previously reported empirical correlation equations have been assessed using the data from the current study. Finally, one correlation was selected for modification and use. The prediction accuracy of the modified equation was found to be high; this was verified after finding a small deviation between the predicted values from the modified correlation and experimental values.
3. Data Processing and Discussion
3.1. Experimental Determination of the Heat Transfer Coefficient
Data are collected using the data acquisition system; this experimental data can be used to calculate the heat transfer coefficient for each heat transfer tube under different operating conditions. Since the tube in the tube heat exchanger is well insulated, the heat loss of the test section can be neglected. The total heat transfer in the test section can be calculated from the heat transfer of the water in the shell of the heat exchanger,
.
where
, is the heat transfer of the water in the shell, W;
, the specific heat capacity corresponding to the average from the inlet and outlet temperature of the water, J/(kg · K);
, the mass flow rate of the water, kg/s;
, the outlet temperature of the water, K;
, the outlet temperature of the water, K.
The total heat transfer of the refrigerant,
is composed of two parts: the sensible heat,
, and the latent heat of vaporization,
.
Based on the thermal balance of the test section, the inlet quality,
, of the refrigerant can be calculated using Equation (5).
, is the saturation temperature of the refrigerant, K; , the latent heat of vaporization of the refrigerant, J/kg;, the mass flow rate of the refrigerant, kg/s; and is the inlet temperature of the refrigerant, K.
The test section refrigerant outlet vapor quality
is calculated using Equation (6).
The logarithmic mean temperature difference LMTD of the test section is calculated using:
with
representing the inlet and outlet temperatures of the refrigerant, K;
representing the inlet and outlet temperatures of water, K.
All of the test tubes are new tubes, so the effect of fouling thermal resistance can be ignored, and the condensation heat transfer coefficient inside the tubes can be calculated using Equation (8).
represents the inner and outer surface areas of the tube, m2; represents the inner and outer diameters of the tube, mm; k represents the thermal conductivity of the tube wall, W/(m · K); represents the inside heat transfer coefficient on the refrigerant side the tube and the heat transfer coefficient on the water side outside the tube, W/(m2·K).
The heat transfer coefficient of water on the shell side is calculated using the Gnielinski [
10] correlation equation:
The fanning friction factor is calculated using the Petukhov [
11] correlation equation:
3.2. Condensation Heat Transfer Coefficient as a Function of Refrigerant Mass Flux
Figure 2 compares the variation of the condensation heat transfer coefficient with the mass flow rate for the smooth and enhanced tubes.
As can be seen in
Figure 2, the condensation heat transfer coefficients of R410A and R32 in the smooth tube increase linearly with the mass flow rate; however, in the three enhanced tubes (HB, HYD, and HYD/HB) the HTC first decreases and then increases. This may be due to the fact that when compared to smooth tubes, enhanced heat transfer is produced in the three enhanced tubes as a result of turbulence and an increase in the heat transfer surface area. Compared to smooth tubes, herringbone tubes have some special protrusions on the inner wall. These structural protrusions disrupt the fluid boundary layer, enhancing fluid disturbance and mixing. Additionally, the extended surface area generated by these protrusions is beneficial for enhanced heat transfer. The enhanced tubes with hydrophobic surfaces create a thin liquid film on the wall, which has lower surface tension. Therefore, when condensing fluid flows through the tubes, the presence of the liquid film makes it easier for the fluid to transition into turbulent flow within the tube. Turbulent flow disrupts the boundary layer and removes the thicker thermal resistance layer, thereby increasing the heat transfer coefficient. At low mass flow rates (below 60 kg·m
−2s
−1), the refrigerant flow velocity is low, and the condensation process inside the tube is mainly the result of surface heat transfer. Turbulent heat transfer may only have a small contribution to the overall heat transfer process, so the heat transfer coefficient is high. As the mass flow rate increases, the refrigerant flow velocity increases; the heat transfer process gradually transitions from surface heat transfer to turbulent heat transfer. Turbulence suppresses the heat transfer process, resulting in a decrease in the heat transfer coefficient. As the flow rate continues to increase, the turbulence intensity reaches its peak and then gradually decreases. The degree of mixing between the refrigerant and the heat transfer surface gradually increases and, therefore, the heat transfer coefficient also gradually increases. This indicates that the three enhanced tubes are more suitable for condensation heat transfer under low mass flow rate (below 60 kg·m
−2s
−1) and high mass flow rate conditions (above 120 kg·m
−2s
−1). Additionally, the heat transfer coefficient is slightly lower for medium mass flow rates (60~120 kg·m
−2s
−1) and condensation heat transfer designs in this range should be avoided.
The condensation heat transfer coefficient of refrigerant R32 is larger than that of R410A in all the tubes (under identical conditions). This can be explained from the physical properties of the two refrigerants; R32 has a smaller molecular weight and has a lower density and higher thermal conductivity. This means that in the condensation process inside the tube (when using R32) can release heat quicker and transfer it to the surface of the tube. Additionally, R410A has a higher density and viscosity; this can lower its flow rate and flow state inside the tube. This will reduce the condensation heat transfer coefficient inside the tube.
When two refrigerants are condensed under different conditions, the heat transfer coefficients of the three enhanced tubes are all enhanced to varying degrees. The tube with the largest enhancement ratio (enhanced tube HTC/smooth tube HTC) is the HB tube (1.55), followed by the HYD tube (1.4) and the HYD/HB tube (combines the hydrophobic surface and the herringbone microfin) has the lowest enhancement ratio (1.37). Both hydrophobic surfaces and herringbone microfins can enhance the condensation heat transfer inside the tube because they can increase the apparent contact angle of the surface and also reduce the contact area of the liquid condensate on the surface. This reduces the liquid film resistance and heat resistance of the heat transfer surface. Specifically, the hydrophobic surface can reduce surface tension and cause the condensate droplets to spontaneously detach from the surface. This reduces the thickness of the surface liquid film and increases the contact angle between the droplets and the surface. This improves the heat transfer coefficient. Herringbone microfins increase the heat transfer surface area and reduce the thickness of the liquid film at the surface. This reduces the heat resistance of the heat transfer surface.
However, the combination of the hydrophobic surface and the herringbone microfin results in a decrease in the condensation heat transfer coefficient. This may be because the combination increases the evaporation of the surface liquid film, causing a consumption of heat that lowers the condensation heat transfer coefficient. At the same time, the condensate droplets on the hydrophobic surface are blocked by the microfins and the droplets cannot completely detach from the surface; therefore, the diffusion and evaporation of the liquid film between the microfins consume more heat. This combination also increases the surface heat resistance and reduces the heat transfer coefficient.
3.3. Empirical Correlation for Smooth Tube Condensation
Table 2 lists several previously reported empirical correlations that model the condensation heat transfer coefficient for smooth tubes.
Figure 3 compares the condensation heat transfer coefficient (based upon experimental data) in smooth tubes and the predicted results using five empirical models. The data points were collected under experimental conditions obtained using refrigerants R410A and R32 for saturated temperatures of 35 °C and 45 °C and a vapor quality ranging from 0.8 to 0.2.
As shown in
Figure 3, the majority of the data points fall in an approximate error range of ±40%; therefore, the empirical models provide a reasonable prediction of results. Data points with large errors appear for lower experimental conditions. In this range the Cavallini et al. [
13] and Shah [
12] relations significantly underestimated the heat transfer coefficient. At those conditions, the deviation of some data points is more than −40%. Predicted values, when using the correlation of Bohdal et al. [
15] for R32, showed a good linear relationship (closer to the experimental values) for the experimental values at conditions of both 35 °C and 45 °C. Predicted values when using R410A were lower than the experimental values. Here, the overall predicted values are relatively accurate (typical deviation ranging from −20% to −40%). Values predicted using the Cavallini and Zecchin [
14] relation showed a large variation trend with increasing HTC. Predicted values were lower than experimental values in the low HTC region (maximum deviation of −40%). Predicted values were higher in the high heat transfer coefficient region (maximum deviation being approximately +20%).
Figure 3.
Comparison of experimental condensation heat transfer coefficient values in ST tubes with values obtained using various empirical models [
12,
13,
14,
15,
16].
Figure 3.
Comparison of experimental condensation heat transfer coefficient values in ST tubes with values obtained using various empirical models [
12,
13,
14,
15,
16].
3.4. Development of New HTC Correlations for Enhanced Tubes
Cavallini et al. [
13] developed a prediction model using 5648 data points from eighteen references that covered a wide range of operating conditions. They also considered the influence of different flow patterns on the convective heat transfer coefficient inside the tube. Since the prediction accuracy for convective condensation heat transfer inside the tube was high, this relation has been widely recognized.
Conditions (refrigerants used (R410A and R32), saturated temperatures, vapor qualities, and tube diameters) considered in the current experiment were within the range of the Cavallini et al. model. These facts made the model developed by [
13] the correlation that could be modified for use with the enhanced tubes from this study (HYD, HB, and HYD/HB tubes).
There were differences between the predicted values (using the unmodified Cavallini model) and experimental HTC values for the conditions considered here. In order to account for the enhanced tubes used in this study, an enhancement factor E (as shown in Equations (17) and (18)) is introduced to correct the condensation heat transfer coefficient calculated using the Cavallini et al. [
13] correlation. The enhancement factor E includes parameters to characterize different tube features (i.e., contact angle α and the helix angle β), as well as additional coefficients (B1 and B2) that were introduced in order to make the Cavallini et al. [
13] formula more accurate when using enhanced tubes. The specific parameters of the enhanced tubes used in this study (HYD, HB, and HYD/HB tubes) are listed in
Table 1.
where as
,
,
, , and
, respectively, represent the fin height, fin pitch, contact angle, helix angle, inner diameter, and length of the tubes. The detailed parameters are listed in
Table 1 and include the following: (i) when applied to HYD tubes, the values of B1 and B2 for use in Equation (17) are 1.76 and 0, respectively, (ii) for HB tubes, the values of B1 and B2 for use in Equation (17) are 0 and 0.93, respectively, and (iii) when applied to HYD/HB tubes, the values of B1 and B2 for use in Equation (17) are 0.602 and 0.1, respectively.
A comparison of differences between predicted and experimental HTC values of the enhanced (HYD, HB, and HYD/HB) tubes is shown in
Figure 4,
Figure 5 and
Figure 6. Good results are obtained from the modified Cavallini et al. [
13] formula. Deviations are within (i) −10%~+12% for the HYD tube (see
Figure 4); (ii) −15%~+17% for the HB tube (see
Figure 5); and (iii) −15%~+14% for the HYD/HB tube (see
Figure 6). Predicted deviations for all the enhanced tubes are within ±17%. This indicates good prediction accuracy. This modified Cavallini et al. [
13] formula has the best prediction accuracy when predicting the condensation of R32 in the HB and HYD/HB tubes. For these tubes, the predicted values are close to the experimental values. When using the modified condensation equation for use with R410A, the predicted values (for the three enhanced tubes in the low heat transfer coefficient region) are relatively close to the experimental values. As the HTC increases, the deviation between the predicted and experimental value increases; however, the error does not exceed +17% and this is still in the acceptable range.
The evaluation of correlation equations typically uses statistical parameters such as the mean absolute deviation (MAD) as a measure of how well the correlation performs (see Equation (19)). This parameter measures the prediction accuracy of the correlation equation. It is the measurement of the average absolute difference between actual observed values and the predicted values. An additional parameter, mean relative deviation (MRD), can also be used in order to determine whether the predictive correlation is biased high or low (see Equation (20)). This correlation equation evaluation index measures the average relative deviation between predicted values and the actual data. The use of these two indicators can effectively evaluate the predictive ability of the correlation and will determine whether it is applicable to specific operating conditions and applications.
The calculated MAD and MRD parameters for the three enhanced tubes are shown in
Table 3. For all of the tubes, the MAD of the experimental data points is less than 6% and the MRD is less than 2.4%. This indicates that the modified Cavallini et al. [
13] correlation has high prediction accuracy. Values of MAD for the three enhanced tubes are less than 6%, with the maximum MAD value (5.76%) associated with the HB tube and the minimum value (4.18%) coming from the HYD/HB tube. The MRDs of the three enhanced tubes are all greater than zero, indicating that the predicted values of the modified correlation will be slightly higher than the experimental values. The maximum MRD value comes from the HYD tube (2.33%) and the minimum value comes from the HB tube (0.12%). When considering the MAD parameter, the modified equation produces the best results for the HYD/HB tube; however, if you consider the MRD parameter, the modified formula produces the best results for the HB and HYD/HB tubes.