Robust and General Model to Forecast the Heat Transfer Coefﬁcient for Flow Condensation in Multi Port Mini/Micro-Channels

: A general correlation for predicting the two-phase heat transfer coefﬁcient (HTC) during condensation inside multi-port mini/micro-channels was presented. The model was obtained by correlating the two-phase multiplier, ϕ tp with affecting parameters using the genetic programming (GP) method. An extensive database containing 3503 experimental data samples was gathered from 21 different sources, including a broad range of operating parameters. The newly obtained correlation ﬁts the broad range of measured data analyzed with an average absolute relative deviation (AARD) of 16.87% and estimates 84.73% of analyzed data points with a relative error of less than 30%. Evaluation of previous correlations was also conducted using the same database. They showed the AARD values ranging from 36.94% to 191.19%. However, the GP model provides more accurate results, AARD lower than 17%, by considering the surface tension effects. Finally, the effect of various operating parameters on the HTC was studied using the proposed correlation. of vapor qualities, mass velocities, saturation temperatures, channel diameters, and working ﬂuids.


Introduction
Condensers play important roles in different systems, such as air conditioning [1,2], heat pumping [3], refrigeration [4,5], nuclear [6,7], chemical [8,9], etc. Nowadays, condensers with mini/micro-channels are used instead of conventional channels in many industries [10][11][12][13][14][15]. The main reasons are reducing the size of the condensers, refrigeration charge, required space, and overall cost of the process [16]. Therefore, a lower condenser's weight and size are required for a given heat duty in this type of channel. The multi-port mini/micro-channels can raise the contact region between the condenser and the working fluids and leads to a higher evaporating rate as an internal fin [17]. Furthermore, the influence of surface tension becomes important in the smaller channels. That is, the surface tension can decrease the liquid film thickness during condensation in these channels. Furthermore, heat exchange between parallel channels can lead to a difference in HTC in multi-port channels from conventional ones.
In many previous studies, the two-phase Nusselt number, Nu tp , was obtained by multiplying a two-phase correction factor by the single-phase HTC, Nu tp = Nu sp ϕ tp (2) Here, ϕ tp is the two-phase correction factor. The values of the two-phase correction factor, ϕ tp , depend on several parameters, and different studies have been conducted to Processes 2022, 10, 243 3 of 26 estimate this factor. Shah [45] presented the following simple expression for estimating the two-phase correction factor during film condensation in pipes, Based on this correlation, only the vapor quality and the reduced pressure can affect the ϕ tp . Cavallini et al. [46] proposed a general correlation for various fluids. They indicated that the vapor quality, liquid Pr, viscosity, and density of the liquid/vapor are the important parameters for estimating ϕ tp . In addition, they reported that the HTC is related to the wall temperature difference, when J g > J T g . The transition line for dimensionless vapor velocity, J T g was defined as following equation, where C T value is 1.6 for hydrocarbons and 2.6 for other refrigerants.
Dobson and Chato [41] suggested the use of the Lockhart-Martinelli Parameter, X tt for predicting the two-phase correction factor in annular flow, There are several general correlations for estimating the two-phase HTC during condensation that was developed on the basis of large datasets, which are applicable for both single and multi-port channels [47][48][49][50][51]. Kim and Mudawar [52] proposed a correlation approach on the basis of 4045 data from 28 sources for single-multi port mini-channels with 0.424 to 6.22 mm hydraulic diameters. The analyzed database included 1964 data points for single-port channels and 2081 data points for multi-port channels. They suggested two separated correlations for slug/bubbly and annular flows. The transition line between these regimes was defined as follows, We * > 7X tt 0.2 (6) where modified Weber number, We * , can be calculated as Soliman [53] f or Re l > 1250 (8) In the general correlation developed by Dorao and Fernandino [49], only two parameters were used to estimate the Nu tp . The correlation was suggested for horizontal mini/macro-channels with 0.067 to 20 mm diameters. This model presented two different regimes, i.e., high and low mass fluxes (or vapor qualities) for HTC. They indicated that the pipe diameter effect becomes significant when the two-phase Reynolds number, Re tp is greater than transient value, Re T tp , Re T tp = 10 8 D 1.5 (9) Among the data points analyzed by Dorao and Fernandino [49], only the data from Derby et al. [54] and Wang et al. [55] were for multi-port channels, and the model is more suitable for single-port channels. This point was confirmed by Hosseini et al. [56], Processes 2022, 10, 243 4 of 26 which compared the model of Dorao and Fernandino [49] with a large data for single-port channels and observed a good agreement with the experimental values.
Shah [57] presented general correlations based on 5100 data points for different flow regimes and channel orientations. The analyzed data contained single and multi-port channels with HDs from 0.08 mm to 49 mm. This model provided an improvement for estimating the HTC of hydrocarbons. In addition, the author suggested predictive methods for low Froude and Reynolds numbers. It should be mentioned that the machine learning algorithms have been extensively used in various applications such as biomedical [58], clinical [59,60], solar power prediction [61], food production [62], etc., for sustainable solutions. Zhou et al. [63] used the machine learning algorithms for estimating the HTC in mini/micro-channels using 4882 data points. Among the different methods tried by them, the artificial neural network (ANN) with 15 hidden layers and 22 adjusted parameters showed the best agreement with measured data with a total AARD of 6.80%. Recently, Hosseini et al. [56] developed a general model using the genetic programming (GP) approach for predicting the condensation HTC in single channels, with 5728 data points for a broad range of operating conditions. The model predicted the analyzed data points with a good AARD value of 17.82%. In addition, they used a separate database including 712 data points from 6 independent sources for validating the developed model. The results showed that the proposed model provides the best accuracy for single-port channels.

Contributions of the Present Study
It is seen that the previous researchers used the same methods for estimating the HTC during condensation in single and multi-port channels. Therefore, these correlations estimate the same values for these cases. However, as discussed previously, the HTC in parallel channels differs from that for single-port ones because of the contact area effects, surface tension forces, and heat exchange between parallel channels. In the previous work [56], a general correlation was developed for the HTC inside single port mini/micro and macro-channels using an extensive database. The basic objective of this study is to extend the former attainments by developing a general model for the HTC in flow condensation. Therefore, an extensive database containing 3503 data points was gathered from 21 independent studies, covering a broad range of operating conditions and working fluids. In addition, the intelligent and powerful method of genetic programming (GP) was employed for improving a mathematical equation for HTC during the condensation. Furthermore, the performance of previous correlations in estimating the HTC in multiport channels was also evaluated using the same database. Finally, the predictive trends of the new model were evaluated, and the impact of various parameters on the HTC was discussed.

Genetic Programming
Recently, using intelligent methods for modeling, the thermal and hydrodynamic characterizations of different systems have received much attention [63][64][65][66][67]. In the present paper, GP was used for obtaining a general nonlinear correlation for predicting the HTC during the condensation of fluids inside multi-port mini/micro-channels. The flowchart of the present GP model is shown in Figure 1.
Genetic Programming (GP) is based on genetic algorithms thus it comprises a population. An equation in the population, a member of the population, is presented as a tree structure, thus that its branches contain arithmetic functions (i.e., sin, cos, +, *, etc.), and also its leaves are terminal points (i.e., variables). The initial population has a set of equations that each represents an individual. Combining the individuals of the previous population generates the next population. As such, the new equations are produced by applying the crossover and mutation operations to the equations from the last generation in a way that combines the individuals of the previous generation and creates the next individuals of the new generation. A crossover process operates when two individuals are selected based on their evaluations determined by the objective function, and a mutation process operates on a small percentage of the individuals that result from a crossover. Once the individuals of the next generation meet the problem criteria, the GP process is halted. Genetic Programming (GP) is based on genetic algorithms thus it comprises a population. An equation in the population, a member of the population, is presented as a tree structure, thus that its branches contain arithmetic functions (i.e., sin, cos, +, *, etc.), and also its leaves are terminal points (i.e., variables). The initial population has a set of equations that each represents an individual. Combining the individuals of the previous population generates the next population. As such, the new equations are produced by applying the crossover and mutation operations to the equations from the last generation in a way that combines the individuals of the previous generation and creates the next individuals of the new generation. A crossover process operates when two individuals are selected based on their evaluations determined by the objective function, and a mutation process operates on a small percentage of the individuals that result from a crossover. Once the individuals of the next generation meet the problem criteria, the GP process is halted.
This method presents an explicit mathematical model with high accuracy for predicting the experimental data. Our previous studies [68][69][70] used this method to develop general models with high accuracy for thermal and hydrodynamic characterizations of different systems. This method presents an explicit mathematical model with high accuracy for predicting the experimental data. Our previous studies [68][69][70] used this method to develop general models with high accuracy for thermal and hydrodynamic characterizations of different systems.

Experimental Data Samples
An extensive database containing 3503 HTC data samples during condensation inside multi-port mini/micro-channels was collected from 21 independent studies, which covered a wide range of working fluids, HDs, mass fluxes, saturation temperatures, and thermophysical properties. The operating parameters for each analyzed source are summarized in Table 1. All thermophysical characteristics of the fluids are computed with the REFPROP 9.1 software from NIST [71] at the saturation conditions of the working fluids. The analyzed data were gathered from experimental studies of average Nu tp numbers for refrigerants condensation in multi-port mini/micro-channels. Most of these data are available as figures from the database sources listed in Table 1, which have been extracted using the GetData Graph Digitizer. In addition to the range of the non-dimensional parameters used, Table 2 sets out the refrigerant types, channel geometry and dimension, and the operating conditions for the analyzed data. The range of dimensional and dimensionless factors is presented in Table 2. The restrictions of applicability of the new correlation for approximating the HTC during flow condensation inside multi-port channels are shown in this table.
In this paper, the term ϕ tp given by Equation (2) is obtained using the method of GP. The affecting parameters on ϕ tp in multi-port channels, such as vapor quality, saturation pressure, mass flux, surface tension, channel diameter, and thermophysical properties of the working fluids are considered for developing the new model. These parameters are employed to improve a general correlation for approximating the HTC using the GP. The two-phase multiplier is proposed using the following dimensionless expression,

Error Analysis
AARD is used for evaluating the performance of different correlations in this study. Furthermore, the arithmetic average deviation (AAD) is used to assess the over-estimation or under-estimation of the models. These error estimates can be calculated as, Nu tp, calc − Nu tp, exp Nu tp, exp × 100 (12)

Results and Discussion
The former studies showed that the relationships between the two-phase multiplier and the affecting factors were nonlinear. In the present study, a powerful GP method was introduced for developing a nonlinear equation for estimating the two-phase multiplier. Among the many correlations generated by the GP, the best correlation is given by, where We see that Equation (13) describes a nonlinear correlation, consisting of all affecting dimensionless factors.
To show the performance of the established model for estimating the HTC, the predicted Nusselt number values were compared with the measured values in Figure 2. This result verifies that the GP correlation predictions given by Equation (13) were in the desired range, and most of the data were bounded with an error of less than ±30%. The calculations showed the suggested correlation estimates of the 3503 data points with an AARD value of 16.87%. In addition, predicted data with an error of less than 20% and 30% for the GP correlation were 70.37% and 84.73%, respectively. These results showed excellent accuracy of the correlation obtained by the GP for approximating the condensation HTC inside multi-port channels. desired range, and most of the data were bounded with an error of less than ±30%. The calculations showed the suggested correlation estimates of the 3503 data points with an AARD value of 16.87%. In addition, predicted data with an error of less than 20% and 30% for the GP correlation were 70.37% and 84.73%, respectively. These results showed excellent accuracy of the correlation obtained by the GP for approximating the condensation HTC inside multi-port channels.  (13) and the measured data. The solid points represent the outcomes of Equation (13). The solid line represents the best fitting between experimental and estimated data, and the dash lines represent the ±30% error bounds.
As noted before, there are various correlation approaches for predicting the condensation HTC. In most of the previous correlations, the authors did not consider the difference between single-port and multi-port channels in the prediction of the HTC. A list of several well-known correlations for estimating the HTC inside single and multi-port channels is presented in Table 3. Table 3. Correlations for predicting the HTC during condensation.

Reference(s) Correlation Remarks
Kim and Mudawar [52] f or annular f low (We * < 7X tt ) : C and X are calculated by Kim and Mudawar [89] method.
General correlations for condensation in mini/micro-channels.

Condensations for horizontal plain tubes
Crosser [92] Nu tp = 0.0265Pr l Bohdal et al. [ For G ≤ 200 kg m −2 s −1 : For evaluating the capability and generality of the new correlation developed by the GP method, a comparative investigation is presented in this section. Table 4 shows the AARD and AAD values of previous correlations (Table 3) for approximating the condensation HTC as well as the GP correlation for different analyzed data sources. Based on the results shown in Table 4, the GP correlation (Equation (13)) shows the best results for predicting the condensation HTC in multi-port channels among the available correlation with AARD and AAD values of 16.87% and 2.33%, respectively, for all analyzed data. In addition, the AARD values of the new model for estimating the condensation HTC from most of the analyzed data sources were lower than 20%. These results for estimating the condensation HTC confirm the accuracy and generality of the correlation determined by the GP approach. Among the previous models, the best results were obtained by Kim and Mudawar [52] with a total AARD of 36.94%, which was two times greater than that obtained by the GP correlation. It should be noted that most of the data (2081 of 4045 data points) analyzed by Kim and Mudawar [52] were for multi-port channels. The Dorao and Fernandino [49], Shah (2016) [90], Shah (2019) [57] and Crosser [92] correlations showed almost the same results with total AARD values of 41.54%, 43.39%, 43.45%, and 44.07%, respectively. In addition, the Hosseini et al. [56] model, which was improved based on the experimental data for single port channels, also had fairly big deviations in estimating the HTC inside multi-port channels with an AARD of 44.68%. This observation verified that the HTC in multi-port mini/micro-channels differs from that in single port channels. The Bohdal et al. [93] and Akers et al. [91] had the highest deviations with the experimental data among the previous correlations, with the AARD values of, respectively, 164.97% and 191.19%. Therefore, these correlations did not provide reasonable predictions for multi-port mini/micro-channels. Although the previous models exhibit nice accordance with some data sources, none of the models can be used as a general correlation for predicting the condensation HTC in multi-port mini/micro-channels. The new correlation and all of the previous models result in large deviations from the experimental data reported by Belchi et al. [79] and Pham et al. [85] with AARD values of 52.30% and 42.09%, respectively.
To show the accuracy of the other models, the results estimated by Kim and Mudawar [52], Dorao and Fernandino [49], Shah (2019) [57], Shah (2016) [90], Crosser [92], and Hosseini et al. [56] correlations were examined in contrast to the measured values in Figure 3. This figure confirms the above assertion that the previous models have relatively large deviations from the measured values of the condensation HTC in multi-port mini/micro-channels.
The accuracy evaluation of previous models and our model for estimating the HTC inside multi-port micro and mini-channels based on Kandlikar [19] classification was given in Table 5. This table verifies that GP correlation (Equation (13)) can predict the HTC in both micro and mini-channels reasonably well and exhibits the best agreements with the experimental data among all available models with AARD of 16.60% and 16.91% for micromini channels, respectively. In addition, our approach predicts almost 85% of data for both micro-mini channels with an error of less than 30%. The Kim and Mudawar correlation [52] predicts 63.44% of the data for mini-channels with an error of lower than 30%. However, its overall accuracy for mini-channels is not satisfactory, with an AARD of 36.94%. In addition, its deviations for micro-channels are much more significant, and its AARD for micro-channels is 59.10%.
The models developed by Dorao and Fernandino [49], Shah (2016) [90], Shah (2019) [57], and Crosser [92] do not provide acceptable results for multi-port micro and mini-channels, and their AARD values for all cases are more than 35%. Furthermore, similar to Kim and Mudawar [52] model, their deviations for micro-channels are significantly higher than those for mini-channels. The models of Bohdal et al. [93] and Akers [91] show substantial deviations for both micro and mini-channels, and they cannot be used for multi-port channels.
In the condensation of the fluids in multi-port channels, the surface tension force plays a remarkable role in controlling the HTC, and its effect becomes more critical as channel diameter is reduced [94]. The effect of the surface tension force was explicitly considered in the new correlation given by Equation (12) using the Bond number. This dimensionless number (Bo) quantifies the relative importance of surface tension forces compared to the gravitational ones. Therefore, the Bo number is important for condensation HTC modeling of multi-port channels, especially for smaller diameters. However, in most of the previous correlations, there was no factor accounting for the effect of surface tension force.
It should be noted that some comparisons between the measured HTC reported by different sources were conducted at identical conditions in [94]. A high degree of consistency between analyzed experimental data was observed. It means that the accurate intelligent based models obtained for a broad range of operational and geometrical conditions capable of estimating HTC for new independent sources, properly. Moreover, the GP-based models outperform earlier models when examining data samples other than those used for developing GP models, as reported in the previous studies [56,94].  [57], (e) Crosser [92] and (f) Hosseini et al. [56]. The solid points represent the outcomes of different models. The solid line represents the best fitting between experimental and estimated data, and the dash lines represent the ±30% error bounds.
The accuracy evaluation of previous models and our model for estimating the HTC inside multi-port micro and mini-channels based on Kandlikar [19] classification was given in Table 5. This table verifies that GP correlation (Equation (13)) can predict the HTC in both micro and mini-channels reasonably well and exhibits the best agreements with the experimental data among all available models with AARD of 16.60% and 16.91% for micro-mini channels, respectively. In addition, our approach predicts almost 85% of data for both micro-mini channels with an error of less than 30%. The Kim and Mudawar correlation [52] predicts 63.44% of the data for mini-channels with an error of lower than 30%. However, its overall accuracy for mini-channels is not satisfactory, with an AARD of 36.94%. In addition, its deviations for micro-channels are much more significant, and its AARD for micro-channels is 59.10%.
The models developed by Dorao and Fernandino [49], Shah (2016) [90], Shah (2019) [57], and Crosser [92] do not provide acceptable results for multi-port micro and minichannels, and their AARD values for all cases are more than 35%. Furthermore, similar to  [57], (e) Crosser [92] and (f) Hosseini et al. [56]. The solid points represent the outcomes of different models. The solid line represents the best fitting between experimental and estimated data, and the dash lines represent the ±30% error bounds.
The physical trends of the established model and those of Kim and Mudawar [52], as the most accurate model among the previous correlations, are presented in Figures 4-7.
Here, the influence of the input parameters on the HTC in multi-port channels, are also discussed. Figure 4 illustrates the impacts of the vapor quality and flow mass flux on the HTC during the condensation of the R744 in a multi-port channel with a HD of 0.13 mm at the saturation temperature of 15 • C. The experimental data shows that the condensation HTC increased with increasing the vapor quality and mass flux. This was in accordance with Sikora et al. [95] findings for the HFE 7000 refrigerant in microchannels. According to Coleman and Garimella [96] flow patterns, the flow regime changes from slug and bubbly with thick liquid film to annular with thin liquid film by raising the vapor quality. In addition, increasing the mass flux can increase the shear stress at the two phases interface, and the liquid film thickness reduces. Therefore, increasing the mass flux and vapor quality decreases the liquid film thickness and increases the HTC. Figure 4 shows that the GP correlation shows good accordance with the experimental data. However, the correlation of Kim and Mudawar [52] underestimates the measured values. Here, the influence of the input parameters on the HTC in multi-port channels, are also discussed. Figure 4 illustrates the impacts of the vapor quality and flow mass flux on the HTC during the condensation of the R744 in a multi-port channel with a HD of 0.13 mm at the saturation temperature of 15 °C. The experimental data shows that the condensation HTC increased with increasing the vapor quality and mass flux. This was in accordance with Sikora et al. [95] findings for the HFE 7000 refrigerant in microchannels. According to Coleman and Garimella [96] flow patterns, the flow regime changes from slug and bubbly with thick liquid film to annular with thin liquid film by raising the vapor quality. In addition, increasing the mass flux can increase the shear stress at the two phases interface, and the liquid film thickness reduces. Therefore, increasing the mass flux and vapor quality decreases the liquid film thickness and increases the HTC. Figure 4 shows that the GP correlation shows good accordance with the experimental data. However, the correlation of Kim and Mudawar [52] underestimates the measured values.  (13)) and the model of Kim and Mudawar [52] with the experimental data [37] for R744 condensing at the saturation temperature of 15 ℃ inside a 0.13 mm channel under different vapor qualities and mass fluxes. Blue and red points represent the G = 600 kg. m . s and G = 800 kg. m . s , respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes, and Kim and Mudawar [52] outcomes.
The predicted HTCs for R134a are examined in contrast to the real data for saturation temperatures of 50 and 60 °C in Figure 5. This result showed that the HTCs for the saturation temperature of 50 °C were higher than those for 60 °C. Table 6 presents the thermophysical characteristics of the R134a at the saturation temperatures of 50 and 60  (13)) and the model of Kim and Mudawar [52] with the experimental data [37] for R744 condensing at the saturation temperature of 15°C inside a 0.13 mm channel under different vapor qualities and mass fluxes. Blue and red points represent the G = 600 kg.m −2 .s −1 and G = 800 kg.m −2 .s −1 , respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes, and Kim and Mudawar [52] outcomes.   (13)) and the model of Kim and Mudawar [52] with the experimental data [72] for R134a condensing at the mass flux of 800 kg. m . s inside a 0.33 mm channel under different saturation temperatures. Blue and red points represent the T = 50 °C and T = 60 °C, respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes and Kim and Mudawar [52] outcomes. Figure 6 compares the HTC of R744 at the saturation temperature of 0 ℃ for two different channel diameters. The experimental data and empirical correlations show that the HTCs of R744 in the channels with 0.68 mm hydraulic diameter was slightly higher than those for the channel with 0.78 mm hydraulic diameters. As was noted previously, the surface tension force plays a more critical role in the smaller channels, and the interactions of the phases increase. Therefore, the liquid film becomes thinner, and the HTC increases. Figure 6 shows that the new GP correlation is consistent with the experimental data, while the correlation of Kim and Mudawar [52] overestimated the HTC.  (13)) and the model of Kim and Mudawar [52] with the experimental data [72] for R134a condensing at the mass flux of 800 kg· m −2 · s −1 inside a 0.33 mm channel under different saturation temperatures. Blue and red points represent the T s = 50 • C and T s = 60 • C, respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes and Kim and Mudawar [52] outcomes.
Processes 2022, 10, x FOR PEER REVIEW 22 of 28 Figure 6. Comparisons of the predictions of the GP correlation (Equation (13)) and the model of Kim and Mudawar [52] with the experimental data [82] for R744 condensing at the mass flux of 800 kg · m · s and saturation temperature of 0 °C inside different channel diameters. Blue and red points represent the D = 68 mm and D = 78 mm, respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes and Kim and Mudawar [52] outcomes. Figure 7 shows the effect of the working fluids (R32 and R290) on the HTC. As shown in this figure, the condensation HTCs for R290 as the working fluid were higher than that for R32. The thermophysical properties of the R32 and R290 at the constant saturation temperature of 40 °C are compared in Table 7. Accordingly, R290 had lower liquid and vapor densities. Thus, the liquid and vapor velocities for R290 were higher than those for R32 at the same mass flux. Therefore, the shear stress at the interface increased, which resulted in a higher HTC. In addition, the latent heat of vaporization of R290 was higher Figure 6. Comparisons of the predictions of the GP correlation (Equation (13)) and the model of Kim and Mudawar [52] with the experimental data [82] for R744 condensing at the mass flux of 800 kg·m −2 ·s −1 and saturation temperature of 0 • C inside different channel diameters. Blue and red points represent the D = 68 mm and D = 78 mm, respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes and Kim and Mudawar [52] outcomes.  (13)) and the model of Kim and Mudawar [52] with the experimental data [73,76] for R290 and R32 condensing at the mass flux of 350 kg · m · s and saturation temperature of 40 ℃ inside a 1.16 mm channel. Blue and red points represent the R290 and R32, respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes, and Kim and Mudawar [52] outcomes, respectively.

Conclusions
A new general correlation for predicting the condensation HTC in multi-port mini and micro-channels was developed using the method of genetic programming. For validating the established model, 3503 experimental data points were collected from 21 published papers covering a broad range of working fluids and operating conditions. Key conclusions of this study are: 1. The effects of all affecting parameters ( , , , , ) on the HTC were considered in the GP correlation. Using these parameters led to a correlation that estimates the condensation HTC with reasonable accuracy. The new correlation estimated the HTC with a total AARD of 16.87% for a broad range of data samples.   (13)) and the model of Kim and Mudawar [52] with the experimental data [73,76] for R290 and R32 condensing at the mass flux of 350 kg·m −2 ·s −1 and saturation temperature of 40°C inside a 1.16 mm channel. Blue and red points represent the R290 and R32, respectively. The solid points, solid lines, and dash lines represent the experimental values, GP outcomes, and Kim and Mudawar [52] outcomes, respectively.
The predicted HTCs for R134a are examined in contrast to the real data for saturation temperatures of 50 and 60 • C in Figure 5. This result showed that the HTCs for the saturation temperature of 50 • C were higher than those for 60 • C. Table 6 presents the thermophysical characteristics of the R134a at the saturation temperatures of 50 and 60°C. Reducing the saturation temperature led to higher liquid density and lower vapor density. Therefore, the vapor velocity was increased, and the liquid velocity was reduced and led to a higher velocity difference at the two phases interface. Therefore, higher interfacial shear stress and thinner liquid film were obtained at the lower saturation temperature, and the HTC increased. In addition, R134a had a higher liquid thermal conductivity at the saturation temperature of 50 • C. Therefore, its thermal resistance in liquid film at this temperature was lower, which resulted in a higher HTC. Figure 5 shows that the GP and Kim and Mudawar [52] correlations showed the same trends for different saturation temperatures. However, the GP correlation led to a much better agreement with the experimental data.  Figure 6 compares the HTC of R744 at the saturation temperature of 0°C for two different channel diameters. The experimental data and empirical correlations show that the HTCs of R744 in the channels with 0.68 mm hydraulic diameter was slightly higher than those for the channel with 0.78 mm hydraulic diameters. As was noted previously, the surface tension force plays a more critical role in the smaller channels, and the interactions of the phases increase. Therefore, the liquid film becomes thinner, and the HTC increases. Figure 6 shows that the new GP correlation is consistent with the experimental data, while the correlation of Kim and Mudawar [52] overestimated the HTC. Figure 7 shows the effect of the working fluids (R32 and R290) on the HTC. As shown in this figure, the condensation HTCs for R290 as the working fluid were higher than that for R32. The thermophysical properties of the R32 and R290 at the constant saturation temperature of 40 • C are compared in Table 7. Accordingly, R290 had lower liquid and vapor densities. Thus, the liquid and vapor velocities for R290 were higher than those for R32 at the same mass flux. Therefore, the shear stress at the interface increased, which resulted in a higher HTC. In addition, the latent heat of vaporization of R290 was higher than that of R32. This information also justifies the higher HTCs for R290. Figure 7 also shows that the GP and Kim and Mudawar correlations predict the same trends for different working fluids. However, the GP correlation results were in closer accordance with the real data.

Conclusions
A new general correlation for predicting the condensation HTC in multi-port mini and micro-channels was developed using the method of genetic programming. For validating the established model, 3503 experimental data points were collected from 21 published papers covering a broad range of working fluids and operating conditions. Key conclusions of this study are:

1.
The effects of all affecting parameters (Re l , Pr l , X tt , Bo, x) on the HTC were considered in the GP correlation. Using these parameters led to a correlation that estimates the condensation HTC with reasonable accuracy. The new correlation estimated the HTC with a total AARD of 16.87% for a broad range of data samples. In addition, the percentages of all data with error lower than 20% and 30% for the new model were 70.04% and 84.43%, respectively. 2.
The previous models' predictions of the HTC were also compared to the measured data. The previous correlations showed significantly higher deviations from the experimental data compared to the new correlation. It was found that the previous models can not be considered as general correlations for estimating the condensation HTC in multi-port channels. However, the new GP correlation provides much-improved estimates for the HTCs. 3.
The new model and previous correlations were used for estimating the HTC in channels with different sizes. It was shown that the new correlation estimated the data for micro and mini-channels with the AARD of 16.60% and 16.91%, respectively, and was the most accurate for all cases among all correlations. In addition, the new GP correlation estimated more than 80% of all data for both micro and mini-channels with an error of lower than 30%. This is because using the Bond number in the new model explicitly includes the effects of surface tension in different channel sizes. The previous correlations, however, showed relatively large deviations for all cases. Furthermore, their deviations become larger when the channel's diameter decreases, and all correlations had AARD values of more than 50% for micro-channels. 4.
The new model developed by GP was shown to be suitable for estimating the condensation HTC in multi-port channels over a broad range of vapor qualities, mass velocities, saturation temperatures, channel diameters, and working fluids.
Knowledge of the parameters of greater effect on the target function is critical for developing further simple and precise correlations. Therefore, performing a Pearson's correlation coefficient analysis to select the most influential parameters in controlling the condensation of two-phase HTC in multi-port channels and also to determine the weight of each factor in the model is suggested as future work.