Next Article in Journal
Optimized Power Dispatch for Smart Building and Electric Vehicles with V2V, V2B and V2G Operations
Next Article in Special Issue
Experimental Analysis of an Innovative Electrical Battery Thermal Management System
Previous Article in Journal
Research and Analysis of Promotional Policies for Small Hydropower Generation in Taiwan
Previous Article in Special Issue
Revisiting the Corresponding-States-Based Correlation for Pool Boiling Critical Heat Flux
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Modified Correlative Model for Condensation Heat Transfer in Horizontal Enhanced Tubes with R32 and R410A Refrigerants

1
College of Electromechanical Engineering, Qingdao University of Science and Technology, Qingdao 266061, China
2
Ningbo Jintian Copper Tube Co., Ltd., Ningbo 315000, China
3
Department of Energy Engineering, Zhejiang University, 38 Zheda Road, Hangzhou 310027, China
4
Department of Engineering Technology 1300 Elmwood Avenue, State University of New York College at Buffalo, Buffalo, NY 14222, USA
*
Authors to whom correspondence should be addressed.
Energies 2023, 16(13), 4883; https://doi.org/10.3390/en16134883
Submission received: 10 May 2023 / Revised: 7 June 2023 / Accepted: 20 June 2023 / Published: 22 June 2023
(This article belongs to the Special Issue Enhanced Two-Phase Heat Transfer)

Abstract

:
An experimental study was performed that compared tube side condensation heat transfer characteristics of enhanced tubes (hydrophobic surface tubes (HYD), herringbone micro fin tube (HB), and a composite hydrophobic/herringbone (micro fin) tube (HYD/HB)) to the performance of a smooth tube (ST). The condensation heat transfer coefficient (HTC) was calculated from data that were recorded for smooth and enhanced tubes that had an outer diameter (OD) of 12.7 mm. Data were collected (as a function of mass flow rate) using a couple of refrigerants (R410A and R32), for saturated temperatures of 35 °C and 45 °C, with vapor qualities that ranged from 0.8 to 0.2. Several previously reported smooth tube HTC models were used to calculate values that could be compared to experimentally obtained HTC values. The correlation model that demonstrated the best accuracy (for the conditions considered) was then modified for use with the enhanced tubes from this study. Results from the modified correlation show differences with experimental values that ranged from −10% to +17%; the new modified correlation demonstrates high prediction accuracy. An accurate correlation allows the evaluation of enhanced heat transfer tubes for use in high-efficiency heat exchanger systems. The development of this new model is significant in the study of enhanced heat transfer.

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/m2·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/m2·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.

2. Experimental Apparatus and Test Conditions

2.1. Experimental System

An experimental process was carried out on a high-precision test rig that can use various refrigerants (i.e., R410A, R32, R22, etc.) and can perform single-phase, condensation, and evaporation heat transfer performance tests on various types of heat exchange tubes. The test rig is able to test circular tubes, rectangular channel tubes, internally enhanced tubes, and three-dimensional tubes with outer tube diameters (OD) ranging from 6 mm to 19 mm. Various experimental conditions can be controlled: (i) inlet and outlet vapor qualities (0 to 1); (ii) refrigerant mass flow rate (10 kg/h to 300 kg/h); (iii) boiling conditions ranging from 5 °C to 15 °C; and (iv) condensation conditions ranging from 35 °C to 45 °C. The schematic diagram of the test setup is shown in Figure 1 and comprises a main refrigerant loop with three auxiliary loops.
The refrigerant flows out of the storage tank and enters the preheating section after passing through a digital gear pump and a Coriolis mass flow meter. After being preheated, the refrigerant flows into the test section (the heat transfer tube being evaluated is located here). The refrigerant flows out of the test section and enters the supercooling section. Finally, the refrigerant enters the storage tank to complete the cycle.
The testing section is the most critical part of the experimental system. It consists of a single tube counter-flow heat exchanger (the tube being evaluated is the inner tube of the heat exchanger). The length of the tube heat exchanger can be changed according to the experimental requirements. The inlet and outlet conditions of the heat exchanger are measured with temperature, pressure, and differential pressure sensors. The preheating section consists of the water bath heating system, and functions to control the required test conditions of the refrigerant and the gear pump deliver the fluid to the entrance of the experimental section. The subcooling section cools the refrigerant flowing out of the test section and reduces its temperature and pressure; this brings the refrigerant to a subcooled state and sends it to the storage tank for the next cycle. Temperatures in the test section, preheating section, and subcooling section are all controlled by PID temperature controllers with an accuracy of up to 0.01 °C. The pressure sensor (working temperature ranging from −40 °C to 100 °C) has a range of 0~6 MPa with an accuracy of 0.075% full-scale (FS). The differential pressure sensor is used to measure the pressure difference of the working medium between the inlet and outlet of the test section; it has a range of 50 kPa with an accuracy of 0.075% FS and a working temperature ranging from −40 to 100 °C. Additional experimental details and parameter uncertainties are given in [9].

2.2. Parameters and Operating Conditions

Tubes evaluated include smooth (ST), hydrophobic surface tubes (HYD), a herringbone micro fin tube (HB), and a composite hydrophobic herringbone micro fin tube (HYD/HB). All are 304 L stainless steel with a 12.7 mm OD with a test length of 1.7 m. Two types of herringbone micro fin enhanced tubes were considered, HB and HYD/HB. They have a herringbone pitch of 0.636 mm, a herringbone fin height of 0.052 mm, and a helix angle of 18°. Experimental results show that HB tubes significantly improve condensation heat transfer characteristics. Increasing the fin height increases the inner wall surface area and enhances heat transfer. The fin pitch affects the fluid flow pattern, and the helix angle influences flow behavior and velocity distribution, all contributing to the heat transfer performance. These parameters all affect the heat transfer performance of the tubes and make them suitable parameters to include in the correlation model. Additional tube parameters are given in Table 1. Condensation experiments of the enhanced tubes were carried out using R410A and R32, with saturation temperatures of 35 °C and 45 °C, respectively.

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, Q w , e x p .
Q w , e x p = c p l , w , e x p m w , e x p T w , e x p , o u t T w , e x p , i n
where Q w , e x p , is the heat transfer of the water in the shell, W; c p l , w , e x p , the specific heat capacity corresponding to the average from the inlet and outlet temperature of the water, J/(kg · K); m w , e x p , the mass flow rate of the water, kg/s; T w , e x p , o u t , the outlet temperature of the water, K; T w , e x p , i n , the outlet temperature of the water, K.
The total heat transfer of the refrigerant, Q p h , is composed of two parts: the sensible heat, Q l a t , and the latent heat of vaporization, Q s e n s .
Q p h = Q l a t + Q s e n s = c p l , w , p h m w , p h T w , p h , i n T w , p h , o u t
Q s e n s = c p l , r e f , p h m r e f T s a t T r e f , p h , i n
Q l a t = m r e f h l v x i n
Based on the thermal balance of the test section, the inlet quality, x i n , of the refrigerant can be calculated using Equation (5).
x i n = Q w , e x p c p l , r e f , p h m r e f T s a t T r e f , p h , i n m r e f h l v
where   T s a t , is the saturation temperature of the refrigerant, K; h l v , the latent heat of vaporization of the refrigerant, J/kg;   m r e f , the mass flow rate of the refrigerant, kg/s; and T r e f , p h , i n is the inlet temperature of the refrigerant, K.
The test section refrigerant outlet vapor quality x o u t is calculated using Equation (6).
x o u t = x i n + Q w , e x p m r e f h l v
The logarithmic mean temperature difference LMTD of the test section is calculated using:
LMT = T w , e x p , i n T r e f , e x p , o u t T w , e x p , o u t T r e f , e x p , i n ln T w , e x p , i n T r e f , e x p , o u t T w , e x p , o u t T r e f , e x p , i n
with T r e f , e x p , i n ,   T r e f , e x p , o u t representing the inlet and outlet temperatures of the refrigerant, K; T w , e x p , i n ,   T w , e x p , o u t 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).
h i = 1 A i LMTD Q w , e x p 1 h o A o d o ln d o / d i 2 k A o
where   A i ,   A o represents the inner and outer surface areas of the tube, m2; d o ,   d i represents the inner and outer diameters of the tube, mm; k represents the thermal conductivity of the tube wall, W/(m · K); h i represents the inside heat transfer coefficient on the refrigerant side the tube and h o 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:
h o = f w / 2 R e w 1000 P r w 1 + 12.7 f w / 2 1 / 2 P r w 2 / 3 1 μ b u l k μ w 0.14 k w d h
The fanning friction factor is calculated using the Petukhov [11] correlation equation:
f w = 1.58 R e w 3.28 2

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].
Energies 16 04883 g003

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).
J G T = 7.5 4.3 X tt 1.111 + 1 3 + C T 3 1 3 C T = { 1.6 Hydrocarbon 2.6 Other   refrigerants
J G = x G g d h ρ v ρ l ρ v 1 2
J G > J G T :   h a = h lo [ 1 + 1.128 x 0.8170 ( ρ l ρ v ) 0.3685 ( μ l μ v ) 0.2363 × ( 1 μ v μ l ) 2.44 Pr l 0.1 ]
J G J G T :   h d = h a J c T J G 0.8 h start J c J G T + h start ,
h stra = 0.725 1 + 0.74 1 x x 0.3321 1 × λ l 3 ρ l ρ l ρ v gh lg μ L D T 0.25 + 1 x 0.087 E h lo
h lo = 0.023 Re lo 0.8 P r L 0.4 λ l D
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.
E = B 1 e c o s α 2 + B 2 f h f p e c o s α d i L e c o s β 1 2
B 1 = 1.76 ,           H Y D           0 ,     H B 0.602 ,       H Y D / H B , B 2 =         0 ,     H Y D     0.93 ,     H B       0.1 ,           H Y D / H B
where as f h , f p , α , β , and d i , L , 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.
M A D = 1 N i = 1 N h i p r e d h i e x p h i e x p
M R D = 1 N i = 1 N h i p r e d h i e x p h i e x p
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.

4. Conclusions

The heat transfer coefficient was measured experimentally for smooth (ST), hydrophobic (HYD), herringbone micro fin (HB), and hydrophobic/herringbone (HYD/HB) tubes for various mass flow rates. Several empirical correlation equations were evaluated using the experimental data points and the Cavallini et al. [13] correlation was modified for use with the enhanced tubes. Predicted values were compared with the experimental values. It was found that the deviation between the predicted values and the experimental values was small. This indicates that the prediction accuracy was high.
  • Condensation heat transfer coefficients (when using R410A and R32) were measured experimentally for smooth and enhanced tubes. For the experimental range considered, the condensation heat transfer coefficient of the smooth tube increased linearly with the increasing flow rate; as the flow rates increased the heat transfer coefficient increased (with a larger slope). It was found that the condensation heat transfer coefficients (for enhanced tubes), when using either refrigerant, initially decreased and then increased with the increasing refrigerant flow rate. For flow rates ranging from 60 to 120 kg·m−2s−1, the heat transfer coefficients were low. Use of these enhanced tubes in the mid-range of flow should be avoided.
  • For the same refrigerant flow rate, the condensation heat transfer coefficient when using R32 was significantly higher than when using R410A. For the same conditions, the condensation heat transfer coefficient of R32 was approximately 1.03 to 1.48 times the HTC when using R410A.
  • Five empirical correlation models were compared for their prediction accuracy in smooth tubes. With the exception of a few data points, the errors of most of the models were within ±40%. Among them, the correlation models of Cavallini and Bohal had better prediction accuracy (differences within −40% to +20%). The correlation model of Cavallini et al. [13] was modified for use with HYD, HB, and HYD/HB tubes. The fitted errors were within −10% to +17%, indicating a high prediction accuracy for the modified correlation.
  • The Cavallini et al. [13] correlation was modified and the modified correlation had a mean absolute deviation (MAD) of <6% and the mean relative deviation (MRD) was <2.4%. The maximum MAD value was 5.76% (for the HB tube) and the minimum MAD value was 4.18% (for the HYD/HB tube). The MRD values of all three enhanced tubes were greater than zero, indicating that the predicted values of the new correlation model were slightly higher than the experimental values. The maximum MRD value was 2.33% (for the HYD tube) and the minimum value was 0.12% (for the HB tube). Overall, the modified Cavallini et al. [13] correlation accurately predicts the condensation heat transfer coefficients of the three enhanced tubes (HYD, HB, and HYD/HB). These are important new results that will allow these tubes to be used in the development of high-performance heat transfer systems.

Author Contributions

Conceptualization: J.Z.; Methodology: G.Z., D.D., L.Z., Y.X. and J.D.; Software: J.Z. and J.D.; Validation, J.D.; Formal analysis: G.Z., D.D., Y.X. and J.D.; Investigation: W.L. and J.Z.; Resources: W.L. and D.J.K. Y.X.; Data curation: G.Z., L.Z. and J.D.; Writing—original draft: G.Z., W.L. and J.D.; Writing—review & editing: D.J.K.; Supervision: W.L. and D.J.K.; Project administration: W.L. and D.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Aarea, m2
ddiameter of tube, m
Eenhancement factor
ffrictional factor
fhfin height
fpfin pitch
Gmass velocity, kg m −2 s −1
hheat transfer coefficient, W m −2 K −1
HBherringbone tube
HYDhydrophobic surface tube
HYD/HBhydrophobic herringbone micro fin composite tube
J G modified Froude number
J G T transition nondimensional gas velocity
LMTDlogarithmic mean temperature difference, K
mmass flux, kg s −1
MADmean absolute deviation
MRDmean relative deviation
Nnumber of data points
Ppressure, Pa
PrPrandtl number
Qheat flow rate, W
Rereynolds number, μd/v,
STsmooth tube
T/ttemperature, K/◦C
xvapor quality
XttMartinelli parameter
Greek symbol
αhelix angle, °
βcontact angle, °
μdynamic viscosity, Pa·s
ρdensity, kg m−3
Subscript
ininlet
bulkbulk temperature
expexperimental
lliquid
loliquid phase with total flow
lvlatent heat, J/kg
outouter
predpredict
refrefrigerant
satsaturation
strafully stratified flow regime
tptwo-phase
vvapor
wwater

References

  1. Longo, G.A.; Mancin, S.; Righetti, G.; Zilio, C. Comparative analysis of microfin vs. smooth tubes in R32 and R410A condensation. Int. J. Refrig. 2021, 128, 218–231. [Google Scholar] [CrossRef]
  2. Ma, L.; Liu, X.; Gao, Y.; Li, W.; Wu, Z.; Luo, X.; Tao, Z.; Kabelac, S. R410A and R32 condensation heat transfer and flow patterns inside horizontal micro-fin and 3-D enhanced tubes. Int. Commun. Heat Mass Transf. 2023, 142, 106638. [Google Scholar] [CrossRef]
  3. Alwazzan, M.; Egab, K.; Peng, B.; Khan, J.; Li, C. Condensation on hybrid-patterned copper tubes (I): Characterization of condensation heat transfer. Int. J. Heat Mass Transf. 2017, 112, 991–1004. [Google Scholar] [CrossRef]
  4. Wan, Z.; Hu, X.; Wang, X.; He, Z. Experimental Study on the Boiling/Condensation Heat Transfer Performance of a Finned Tube with a Hydrophilic/Hydrophobic Surface. Appl. Therm. Eng. 2023, 229, 120494. [Google Scholar] [CrossRef]
  5. Ko, Y.M.; Jung, J.H.; Sohn, S.; Song, C.H.; Kang, Y.T. Condensation heat transfer performance and integrated correlations of low GWP refrigerants in plate heat exchangers. Int. J. Heat Mass Transf. 2021, 177, 121519. [Google Scholar] [CrossRef]
  6. Hosseini, S.; Moradkhani, M.; Valizadeh, M.; Zendehboudi, A.; Olazar, M. A general heat transfer correlation for flow condensation in single port mini and macro channels using genetic programming. Int. J. Refrig. 2020, 119, 376–389. [Google Scholar] [CrossRef]
  7. Li, W.; Chen, H.; Wang, J.; Gao, Y.; Wang, X.; He, Y.; Tang, W.; Kukulka, D.J. Condensation heat transfer in annuli outside horizontal stainless steel enhanced tubes. Int. J. Therm. Sci. 2022, 177, 107479. [Google Scholar] [CrossRef]
  8. Shah, M.M. Improved correlation for heat transfer during condensation in mini and macro channels. Int. J. Heat Mass Transf. 2022, 194, 123069. [Google Scholar] [CrossRef]
  9. Li, W.; Wang, J.; Guo, Y.; Gu, Z.; Wang, X.; Sun, Z.; Tang, W.; Kukulka, D.J. Two-phase heat transfer of R410A in annuli outside enhanced tubes with micro-fin and dimple. Int. J. Heat Mass Transf. 2021, 175, 121370. [Google Scholar] [CrossRef]
  10. Gnielinski, V. New Equations for Heat and Mass Transfer in Turbulent Pipe and Channel Flow. Int. Chem. Eng. 1976, 16, 8–16. [Google Scholar]
  11. Petukhov, B.S. Heat Transfer and Friction in Turbulent Pipe Flow with Variable Physical Properties. Adv. Heat Transf. 1970, 6, 503–564. [Google Scholar]
  12. Shah, M.M. A general correlation for heat transfer during film condensation inside pipes. Int. J. Heat Mass Transf. 1979, 22, 547–556. [Google Scholar] [CrossRef]
  13. Cavallini, A.; Col, D.D.; Doretti, L.; Matkovic, M.; Rossetto, L.; Zilio, C.; Censi, G. Condensation in Horizontal Smooth Tubes: A New Heat Transfer Model for Heat Exchanger Design. Heat Transf. Eng. 2006, 27, 31–38. [Google Scholar] [CrossRef]
  14. Cavallini, A.; Zecchin, R. A dimensionless correlation for heat transfer in forced convection condensation. In Proceedings of the 5th International Heat Transfer Conference, Tokyo, Japan, 3–7 September 1974; Volume 3, pp. 309–313. [Google Scholar]
  15. Bohdal, T.; Charun, H.; Sikora, M. Comparative investigations of the condensation of R134a and R404A refrigerants in pipe minichannels. Int. J. Heat Mass Tran. 2011, 54, 1963–1974. [Google Scholar] [CrossRef]
  16. Vu, P.Q.; Kwang-Il, C.; Jong-Taek, O.; Honggi, C.; Taehum, K.; Jungho, K.; Jaeyoung, C. An experimental investigation of condensation heat transfer coefficient using R-410A in horizontal circular tubes. Energy Proc. 2015, 75, 3113–3118. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic diagram of the experimental equipment [9].
Figure 1. Schematic diagram of the experimental equipment [9].
Energies 16 04883 g001
Figure 2. Comparison of condensation heat transfer coefficient as a function of mass flux using R32 and R410A for smooth and enhanced tubes (at a saturated temperature of 45 °C).
Figure 2. Comparison of condensation heat transfer coefficient as a function of mass flux using R32 and R410A for smooth and enhanced tubes (at a saturated temperature of 45 °C).
Energies 16 04883 g002
Figure 4. Comparison of differences between predicted (using the modified Cavallini et al. [13] formula) and experimental values when using the HYD enhanced tube, using R32 and R410A, with saturation temperatures of 35 °C and 45 °C.
Figure 4. Comparison of differences between predicted (using the modified Cavallini et al. [13] formula) and experimental values when using the HYD enhanced tube, using R32 and R410A, with saturation temperatures of 35 °C and 45 °C.
Energies 16 04883 g004
Figure 5. Comparison of differences between predicted (using the modified Cavallini et al. [13] formula) and experimental values when using the HB enhanced tube, using R32 and R410A, with saturation temperatures of 35 °C and 45 °C.
Figure 5. Comparison of differences between predicted (using the modified Cavallini et al. [13] formula) and experimental values when using the HB enhanced tube, using R32 and R410A, with saturation temperatures of 35 °C and 45 °C.
Energies 16 04883 g005
Figure 6. Comparison of differences between predicted (using the modified Cavallini et al. [13] formula) and experimental values when using the HYD/HB enhanced tube, using R32 and R410A, with saturation temperatures of 35 °C and 45 °C tubes.
Figure 6. Comparison of differences between predicted (using the modified Cavallini et al. [13] formula) and experimental values when using the HYD/HB enhanced tube, using R32 and R410A, with saturation temperatures of 35 °C and 45 °C tubes.
Energies 16 04883 g006
Table 1. Geometric parameters of the tubes evaluated.
Table 1. Geometric parameters of the tubes evaluated.
ParameterST TubeHYD TubeHB TubeHB/HYD Tube
MaterialStainless steel Stainless steel Stainless steelStainless steel
Inner diameter (mm)11.511.511.511.5
Outer diameter (mm)12.712.712.712.7
Thickness (mm)0.610.610.610.61
Length (m)1.71.71.71.7
Fin height (mm)--0.052 0.052
Fin pitch (mm)--0.636 0.636
Helix angle (°)--1818
Addendum angle (°)--161.4161.4
Contact angle (°)-107.3548.2353.37
Table 2. Empirical correlation for condensation heat transfer coefficient of smooth tubes.
Table 2. Empirical correlation for condensation heat transfer coefficient of smooth tubes.
AuthorsEquationsExperimental Conditions
Shah et al. [12] (1979) h t p = h L 1 x 0.8 + 3.8 x 0.76 1 x 0.04 P r 0.38
h L = 0.023 Re L 0.8 Pr l 0.4 k l D
Re l = G D μ l ,   Pr l = c p l μ l k l ,   P r = P P c
Fluid: water, R11, R22, alcohol, etc.
G = 10.8~210.6 kg m−2s−1
Tsat = 21~31 °C
do = 7~40 mm
Cavallini et al. [13] (2006) J G T = 7.5 4.3 X tt 1.111 + 1 3 + C T 3 1 3 C T = { 1.6 Hydrocarbon 2.6 Other   refrigerants
J G = x G g d h ρ v ρ l ρ v 1 2
J G > J G T :   h a = h lo [ 1 + 1.128 x 0.8170 ( ρ l ρ v ) 0.3685 ( μ l μ v ) 0.2363 × 1 μ v μ l ) 2.44 Pr l 0.1
J G J G T :   h d = h a J c T J G 0.8 h start J c J G T + h start ,   h lo = 0.023 Re lo 0.8 P r L 0.4 λ l D
h stra = 0.725 1 + 0.74 1 x x 0.3321 1 × λ l 3 ρ l ρ l ρ v gh lg μ L D T 0.25 + 1 x 0.087 h lo
Fluid: R410A, R32, R22, etc.
G = 18~2240kg m−2s−1
Tsat-Twall = 0.6~28.7 °C
Tsat = −15~302 °C
do = 3.1~17mm
Cavallini and Zecchin [14] (1974) h l p = 0.05 R e e q o   s Pr i o . 33 k l d i
Re eq = Re v μ v μ l v ρ l ρ v 0.5 + Re l
Re l = G 1 x d i μ v
Fluid: R22, etc.
Re ≥ 1.2 × 103
0.1 ≤ x ≤ 0.9
Bohdal et al. [15] (2011) h = 25.084 R e l 0.258 P r l 0.495 p r 0.288 ( x 1 x ) 0.266 k l D Fluid: R134a, R404A
G = 64~902 kg m−2s−1
do = 0.31~3.30 mm
Tsat = 20~40 °C
Vu et al. [16] (2015) h = 0.023 Re l 0 0.8 P r l 0.4 λ l d i 1 x 0.8 + 3.8 P r 0.52 x 1 x 0.201 Fluid: R410A
G = 200~320 kg m−2s−1
do = 6.61~9.2 mm
Tsat = 48 °C
Table 3. MAD and MRD values of the enhanced tubes (HYD, HB, HYD/HB) considered in this study.
Table 3. MAD and MRD values of the enhanced tubes (HYD, HB, HYD/HB) considered in this study.
HYDHBHYD/HB
MAD4.50%5.76%4.18%
MRD2.33%0.12%1.49%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, G.; Du, D.; Zhang, L.; Xiang, Y.; Li, W.; Zhang, J.; Du, J.; Kukulka, D.J. A Modified Correlative Model for Condensation Heat Transfer in Horizontal Enhanced Tubes with R32 and R410A Refrigerants. Energies 2023, 16, 4883. https://doi.org/10.3390/en16134883

AMA Style

Zhang G, Du D, Zhang L, Xiang Y, Li W, Zhang J, Du J, Kukulka DJ. A Modified Correlative Model for Condensation Heat Transfer in Horizontal Enhanced Tubes with R32 and R410A Refrigerants. Energies. 2023; 16(13):4883. https://doi.org/10.3390/en16134883

Chicago/Turabian Style

Zhang, Gangan, Dehui Du, Le Zhang, Yanlong Xiang, Wei Li, Jiapei Zhang, Jincai Du, and David J. Kukulka. 2023. "A Modified Correlative Model for Condensation Heat Transfer in Horizontal Enhanced Tubes with R32 and R410A Refrigerants" Energies 16, no. 13: 4883. https://doi.org/10.3390/en16134883

APA Style

Zhang, G., Du, D., Zhang, L., Xiang, Y., Li, W., Zhang, J., Du, J., & Kukulka, D. J. (2023). A Modified Correlative Model for Condensation Heat Transfer in Horizontal Enhanced Tubes with R32 and R410A Refrigerants. Energies, 16(13), 4883. https://doi.org/10.3390/en16134883

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop