Model of R134a Liquid–Vapor Two-Phase Heat Transfer Coe ﬃ cient for Pulsating Flow Boiling in an Evaporator Using Response Surface Methodology

: The present study proposes a model to predict the heat transfer coe ﬃ cient in R134a liquid–vapor two-phase pulsating ﬂow boiling in an evaporator using the experimental data and response surface methodology (RSM). The model is based on the current existing empirical correlation for R134a liquid–vapor two-phase continuous ﬂow with an imposed modiﬁcation factor. The model for the imposed modiﬁcation factor is the function of the pulsating period and inlet / outlet vapor quality, which is obtained using the limited experimental data. An analysis of variance (ANOVA) is carried out to test the signiﬁcance of the model and normal probability of residuals is analyzed as well. Results show that the regression model produces a mean error of − 4.3% and a standard deviation of 15.4%, compared to experimental results. Of the data 95.1% is contained inside a ± 50% error window, which indicates that the proposed model could predict the heat transfer coe ﬃ cient of R134a liquid–vapor two-phase pulsating ﬂow boiling well.


Introduction
Numerous strategies, including passive and active techniques such as a fin-tube, coiled wires, and swirl/pulsating flow generator are widely used to improve the thermal performance of heat exchangers [1][2][3]. Pulsating flow (pulsed flow), characterized as the periodic change of one or some flow parameters, occurs in numerous systems (refrigeration, chemical engineering, etc.), and has a special significance in the operation of modern industrial equipment technologies. Studies on pulsating flow in a heat exchanger are essential because they are beneficial to understanding flow characteristics and improve heat transfer performance. To date, heat transfer characteristics in pulsating flow have been investigated previously by many researchers with conflicting conclusions, especially for the liquid-vapor two-phase flow [4][5][6].
Different from the single-phase flow, the two-phase heat transfer coefficient (HTC) predictions require not only a broad, accurate experimental database, but also a good knowledge of the heat transfer mechanism. For this purpose, many research proposed different empirical correlations to predict the HTC in the two-phase flow. One classical type is based on the HTC of liquid flow with Dittus and Bolter correlation [7] and on some dimensionless numbers (boiling number Bo, convective number Co, Froude number Fr, Martinelli parameter X tt , etc.), which can be usually written as follows: h tp = h l · f (Bo, Fr, X tt , . . .) (1) Energies 2020, 13 Shah presented a correlation by considering the nucleate boiling component and the forced convection component, and selected the larger of the two components [8]. Saiz Jabardo et al. proposed a correlation that is related to Bo, Fr, and X tt [9]. The models of Chaddock and Buzard [10], Wattelet et al. [11], and Panek [12] mainly are the functions of X tt . The second category is the power-type multiplication model for two boiling components that is analogous to the Dittus-Boelter correlation, and also can be written as: where A is the area of the cross section. Tran et al. carried out R12 flow boiling transfer experiments in small channels and proposed a correlation with a function of Bo, We l , and ρ l /ρ g [13]. Kew and Cornwell [14] observed an increase of HTC with the vapor quality in larger tubes and modified the Lazarek and Black's correlation [15]. In term of application, these two categories of correlations are simple as they do not require the identification of the flow pattern. However, Silva Lima et al. conducted an experimental study on flow boiling in a horizontal smooth tube for R134a, and found the correlation above mentioned tended to underpredict the experimental results [16]. Mikielewicz presented a new method for flow boiling HTC in both conventional channels and minichannels by the hypothesis that energy dissipation resulting from the friction is linked to the heat transfer in the flow boiling process [17]. By the analyses for various experimental investigations in flow boiling, it is found that both the convective and nucleate boiling heat transfer mechanism played import roles in the two-phase heat transfer process [18]. The third format usually can be written as: where n is the power of number. The factor S (suppression factor, <1) reflects a lower effective superheat available in forced convection, which is opposed to pool boiling [19], because of the thinner boundary layer. The factor F (>1) indicates that much higher velocities could enhance forced convection heat transfer for the two-phase flow compared to the single-phase fluid flow. Gungor and Winterton assumed that the convective and nucleate boiling heat transfer have the same effect, and presented a general correlation for forced convection boiling in tubes and annuli with n = 1 [20]. However, they pointed that the HTC was overpredicted for the high vapor quality and underpredicted for the low quality [21]. It indicates that the simple superposition of the two heat transfer mechanism is not suitable. Therefore, Liu et al. modified the superposition relationship of nucleate boiling heat transfer and forced convection heat transfer in the form of the square root [22]. Based on superposition methods and considering the flow pattern characteristics, Kattan et al. proposed a flow pattern based method to develop the correlation [23]. The assumption of the linear variation of the dry angle without the consideration of flow patterns leads to the underprediction at low vapor quality. Subsequently, Wojtan updated the calculation of dry angle for slug flow, stratified-wavy flow and slug-stratified wavy flow to avoid any jump in HTC at any transition boundary and modified Kattan's correlation [24]. According to the research results from the literatures, heat transfer in two-phase flow boiling is very complex, and related to many parameters such as the mass flow rate, heat flux, vapor quality, flow pattern, even fluid properties, and so on. Results show the different HTC correlations in flow boiling have different accuracy for different working conditions. It should be noted that all the correlations above are developed for the liquid-vapor two-phase continuous flow. The introduction of pulsating flow would further increase the complexity of two phase flow boiling as the flow, the liquid-vapor interface and flow pattern is always dynamic instead of a fully developed state in continuous flow. In our previous study, the heat transfer in two-phase pulsating flow is strongly affected by the pulsation period, mass flux, and vapor quality. Additionally, the heat transfer enhancement mechanism of heat transfer is the change of the liquid-vapor interface and wetted wall caused by the flow regime transition [25,26]. To the best of our knowledge, few literatures are devoted to obtain the HTC model for the two-phase pulsating flow. Besides, the correlation of the heat transfer Energies 2020, 13, 3540 3 of 19 coefficient in the two-phase flow needs a large database with a huge number of experiments. Response surface methodology (RSM) can obtain the regression equation of factors and responses by limited experimental designs and both considering the interaction relationship between factors. Therefore, the main objective of this study is to develop an HTC model for the two-phase pulsating flow based on the current experimental data and RSM.
In this work, a model has been proposed to predict HTC for the R134a liquid-vapor two-phase pulsating flow. In this model, the HTC in the two-phase pulsating flow is predicted using the HTC in continuous flow multiplied by a ratio factor. The ratio factor is the function of the pulsating period, mass flux, and vapor quality. One hundred and sixty experimental data under a variable mass flux, pulsating period, and vapor quality are used as the data base and processed with RSM to obtain the ratio factor of HTC of pulsating flow to that of continuous flow. Although the model is based on limited experimental data, it provides a simple tool to predict heat transfer of the pulsating flow, which could be extremely valuable for the pulsating system design.

Pulsating Flow Boiling Facility
Pulsating flow is always generated by the nature of system or flow instability, which is not controlled. In our experiment, the flow is pulsed with a setting period controlled by a solenoid valve. For the comparison of both continuous flow and pulsating flow, two identical loops with two identical evaporators (HX1 and HX2, see Figure 1) were created. For the continuous flow, one of solenoid valves keeps open and other is closed, a continuous flow enters one evaporator exchanging heat with ambient air. For the pulsating flow, the two solenoid valves alternately open and close with a periodic time setting (when one solenoid valve opens the other was closed, and vice versa). In this way, the flow through the evaporators is pulsed, and is steady outside the evaporators. Therefore, the two loops are never closed simultaneously when pulsation occurs, which could ensure a steady state inlet pressure and a continuous mass flux before the evaporator. heat transfer coefficient in the two-phase flow needs a large database with a huge number of experiments. Response surface methodology (RSM) can obtain the regression equation of factors and responses by limited experimental designs and both considering the interaction relationship between factors. Therefore, the main objective of this study is to develop an HTC model for the two-phase pulsating flow based on the current experimental data and RSM.
In this work, a model has been proposed to predict HTC for the R134a liquid-vapor two-phase pulsating flow. In this model, the HTC in the two-phase pulsating flow is predicted using the HTC in continuous flow multiplied by a ratio factor. The ratio factor is the function of the pulsating period, mass flux, and vapor quality. One hundred and sixty experimental data under a variable mass flux, pulsating period, and vapor quality are used as the data base and processed with RSM to obtain the ratio factor of HTC of pulsating flow to that of continuous flow. Although the model is based on limited experimental data, it provides a simple tool to predict heat transfer of the pulsating flow, which could be extremely valuable for the pulsating system design.

Pulsating Flow Boiling Facility
Pulsating flow is always generated by the nature of system or flow instability, which is not controlled. In our experiment, the flow is pulsed with a setting period controlled by a solenoid valve. For the comparison of both continuous flow and pulsating flow, two identical loops with two identical evaporators (HX1 and HX2, see Figure 1) were created. For the continuous flow, one of solenoid valves keeps open and other is closed, a continuous flow enters one evaporator exchanging heat with ambient air. For the pulsating flow, the two solenoid valves alternately open and close with a periodic time setting (when one solenoid valve opens the other was closed, and vice versa). In this way, the flow through the evaporators is pulsed, and is steady outside the evaporators. Therefore, the two loops are never closed simultaneously when pulsation occurs, which could ensure a steady state inlet pressure and a continuous mass flux before the evaporator.    Figure 1 shows a schematic of the two-phase pulsating flow boiling experimental system [25,26]. There are three loops in the system: refrigerant loop (red line), air loop (orange line), and chilled loop (blue line). In the refrigerant loop, the subcooled refrigerant R134a is pumped by a gear pump from a subcooler. Then the fluid passes through a flow meter and is heated by a pre-heater to control the vapor quality before the evaporators. After that, the two-phase R134a enters the two same sub-loops, starting from the solenoid valves, then the heat exchangers (HX1 and HX2), where it exchanges heat with air flow in a wind tunnel. After the heat exchanger, R134a is heated into superheated vapor by the superheaters. Finally, the two sub-loops meet and the vapor passes the condenser and is cooled by a chilled system. The heat exchanger used in our experiment was the fin-tube heat exchanger with 385 mm × 300 mm × 30 mm (length × width × thickness). The inner diameter of copper tube was 6.2 mm. The plain aluminum fins had a fin density of 15 fpi and fin thickness of 0.3 mm. The heat exchanger is shown in Figure 2.
Energies 2020, 13, x FOR PEER REVIEW 4 of 19 Figure 1 shows a schematic of the two-phase pulsating flow boiling experimental system [25,26]. There are three loops in the system: refrigerant loop (red line), air loop (orange line), and chilled loop (blue line). In the refrigerant loop, the subcooled refrigerant R134a is pumped by a gear pump from a subcooler. Then the fluid passes through a flow meter and is heated by a pre-heater to control the vapor quality before the evaporators. After that, the two-phase R134a enters the two same sub-loops, starting from the solenoid valves, then the heat exchangers (HX1 and HX2), where it exchanges heat with air flow in a wind tunnel. After the heat exchanger, R134a is heated into superheated vapor by the superheaters. Finally, the two sub-loops meet and the vapor passes the condenser and is cooled by a chilled system. The heat exchanger used in our experiment was the fin-tube heat exchanger with 385 mm × 300 mm × 30 mm (length × width × thickness). The inner diameter of copper tube was 6.2 mm. The plain aluminum fins had a fin density of 15 fpi and fin thickness of 0.3 mm. The heat exchanger is shown in Figure 2. The chilled loop is used to cool the refrigerant vapor in the condenser and subcooler. The chilled fluid is the mixture of water and ethylene glycol with a volume ratio of 1:1. The temperature of chilled fluid is measured before and after the condenser. The air loop provides air flow to exchange heat with two-phase R134a in the heat exchangers. The air flow rate is measured by a nozzle system. The measurement method and design of the nozzle refer to the American Society of Mechanical Engineers (ASME) Standard MFC-3M [27] and ASME Standard PTC 19.5 [28]. The air flow is circulated by a blower controlled by a speed controller. The temperature of the air flow is measured before and after the heat exchangers.

Test Conditions and Data Reduction
The experiments were designed to ensure that HTC could be evaluated only under conditions of the liquid-vapor two-phase flow in the heat exchanger tubes. To be specific, the vapor quality at the inlet of the evaporator is controlled by adjusting the input power of the preheater according to the mass flow rate. The product of the refrigerant mass flow rate multiplying the enthalpy difference before and after the preheater is used to estimate the input power of the preheater. The outlet vapor quality is controlled by adjusting the air flow rate and the input power of the superheater. Additionally the product of the refrigerant mass flow rate multiplying the enthalpy difference before and after the superheater was used to estimate the input power of the superheater. The matrix of the experimental design is given in Table 1. The chilled loop is used to cool the refrigerant vapor in the condenser and subcooler. The chilled fluid is the mixture of water and ethylene glycol with a volume ratio of 1:1. The temperature of chilled fluid is measured before and after the condenser. The air loop provides air flow to exchange heat with two-phase R134a in the heat exchangers. The air flow rate is measured by a nozzle system. The measurement method and design of the nozzle refer to the American Society of Mechanical Engineers (ASME) Standard MFC-3M [27] and ASME Standard PTC 19.5 [28]. The air flow is circulated by a blower controlled by a speed controller. The temperature of the air flow is measured before and after the heat exchangers.

Test Conditions and Data Reduction
The experiments were designed to ensure that HTC could be evaluated only under conditions of the liquid-vapor two-phase flow in the heat exchanger tubes. To be specific, the vapor quality at the inlet of the evaporator is controlled by adjusting the input power of the preheater according to the mass flow rate. The product of the refrigerant mass flow rate multiplying the enthalpy difference before and after the preheater is used to estimate the input power of the preheater. The outlet vapor quality is controlled by adjusting the air flow rate and the input power of the superheater. Additionally the product of the refrigerant mass flow rate multiplying the enthalpy difference before and after the superheater was used to estimate the input power of the superheater. The matrix of the experimental design is given in Table 1. The meausred parameters before and after the evaporator (HX) is shown in Figure 3. The air inlet/outlet temperature and flow rate were measured to calculate the heat exchange (q a ) at the air side as shown in Equation (4). The refrigerant mass flux (m a ), subcooling temperature (T ri ) before the preheater, the power of the preheater (q hin ) and the super heater (q hsup ), and the superheat temperature (T rsup ) after the superheater were measured to calculate the heat exchange (q r ) at the refrigerant side as shown in Equation (5). The overall heat transfer coefficient U was calculated based on the average heat transfer rate q ave using the LMTD method, as shown in Equation (7). Finally the refrigerant side heat transfer coefficient h r was derived by subtracting the thermal resistance of air side and wall heat conduction.
where h ri and h rsup are the enthalpy of the refrigerant before the preheater and after the superheater. h r and h a are the refrigerant and air side HTC. r o and r i are the inner and outer diameter of the tube.
A o is the heat transfer ares of air side. η is the fin efficiency. Additionally, h a is calculated using correlation Kim-Youn-Webb's correlation [29].
Energies 2020, 13, x FOR PEER REVIEW 5 of 19 The meausred parameters before and after the evaporator (HX) is shown in Figure 3. The air inlet/outlet temperature and flow rate were measured to calculate the heat exchange (qa) at the air side as shown in Equation (4). The refrigerant mass flux (ma), subcooling temperature (Tri) before the preheater, the power of the preheater (qhin) and the super heater (qhsup), and the superheat temperature (Trsup) after the superheater were measured to calculate the heat exchange (qr) at the refrigerant side as shown in Equation (5). The overall heat transfer coefficient U was calculated based on the average heat transfer rate qave using the LMTD method, as shown in Equation (7). Finally the refrigerant side heat transfer coefficient hr was derived by subtracting the thermal resistance of air side and wall heat conduction.
where hri and hrsup are the enthalpy of the refrigerant before the preheater and after the superheater. hr and ha are the refrigerant and air side HTC. ro and ri are the inner and outer diameter of the tube. Ao is the heat transfer ares of air side. η is the fin efficiency. Additionally, ha is calculated using correlation Kim-Youn-Webb's correlation [29]. The sensors, ranges, and accuracy of the measured parameters are listed in Table 2. The measurement accuracy caused an averaged error of ±10.4% for the overall HTC, and ±12% -±16.5% for the refrigerant side HTC. More detail about the experimental approach, test conditions, and data source can be found in our prior work [25,26]. The sensors, ranges, and accuracy of the measured parameters are listed in Table 2. The measurement accuracy caused an averaged error of ±10.4% for the overall HTC, and ±12%-±16.5% for the refrigerant side HTC. More detail about the experimental approach, test conditions, and data source can be found in our prior work [25,26].

HTC Correlations in the Two-Phase Continuous Flow
In this section, the experimental HTC data in the two-phase continuous flow were compared to several existing types of correlations, which are listed in Table 3. The mean error and standard deviation of predicted HTC by different correlations were analyzed. The suitable correlation was selected to predict the results in the two-phase continuous flow, which will be one important component of the HTC in the two-phase pulsating flow. Table 3. Heat transfer coefficient (HTC) correlations in previous literature.

Author Correlation
Mikielewicz [17] h tp 25 For conventional channels, m = 0 For small-diameter channels and mini channels, m = −1 Figure 4 presents the comparisons of the experimental HTC data in the two-phase continuous flow to the different types of existing empirical correlations. It can be seen that the correlations proposed by Shah, Gungor, Liu, and Kattan tended to overpredict the experimental results with mean errors of 42.7%, 37%, 31%, and 40.9%, respectively. Only about 50% of the data fell inside ±30% error window for these four correlations. The dispersion occurred mainly in the higher inlet vapor quality region with a narrower vapor quality range (mainly intermittent and annular flow pattern regions). Compared with the four correlations above, the correlations developed by Wojtan and Mikielewicz had higher accuracy to predict the experimental data with small mean errors of −4.6% and 10.4%. These two correlations predicted 86.7% and 79.4% of the data inside ±30% error windows, respectively.
The summary of all the statistical analysis from the comparisons between experimental data and predicted values is given in Table 4. Although the mean absolute error in Wojtan's model was larger than that in Mikielewicz's correlation, Wojtan's model had the most data points inside an error window of ±30% and ±50%. Therefore, the Wojtan's model was selected to predict HTC in the R134a two-phase continuous flow.
Energies 2020, 13  The summary of all the statistical analysis from the comparisons between experimental data and predicted values is given in Table 4. Although the mean absolute error in Wojtan's model was larger than that in Mikielewicz's correlation, Wojtan's model had the most data points inside an error window of ±30% and ±50%. Therefore, the Wojtan's model was selected to predict HTC in the R134a two-phase continuous flow.

Description about RSM
The RSM, firstly proposed by Box and Wilson [30], is a method to accurately predict the input-response relationships by fully considering the parameter interaction in engineering systems [31][32][33]. It has been popularized in many fields for the design, development, and promotion of new products, as well as in the optimization of existing product designs. With the aid of limited experimental designs, RSM can achieve the regression model of parameters and responses, and intuitively show the correlation curve between parameters and responses using the contour plot and response surface. Therefore, the RSM is used to rearrange the experimental design and further develop the correlation of HTC ratio in the two-phase pulsating flow.
Usually, a second-order model in RSM can be written as follows: where y is the response, x i and x j represent independent variables; n is the number of variables; β 0 , β i , β ii , and β ij (i = 0, 1, 2, . . . , n; j = 0, 1, 2, . . . , n) are coefficients for the intercept, linear, quadratic, and interaction terms in the regression model respectively; and ε is the statistical error. The generic form of HTC correlation is usually written as follows: where x 1 , x 2 , . . . , x n are general variables closed to HTC. a 1 , a 2 , and a n are the power of the number. Additionally, C is the constant term and n is the number of variables. Taking the logarithm of both sides in Equation (9), one can get Neglecting the statistical error and rewriting Equation (8) In order to declare the reason why the second-order RSM model is used, the comparisons of Equations (10) and (11) are listed in Table 5.
It can be seen from the Table 5 that Equations (10) and (11) have a similar format, indicating that the RSM is feasible to obtain the HTC correlation.  (10) and (11).

Rearrangement of Experimental Data
In the two-phase pulsating flow, the pulsation period, mass flux, and vapor quality all have remarkable effects on the heat transfer performance with a strong interaction. In order to facilitate the analysis, a dimensionless variable called the Strouhal number (St), which represents the ratio of the inertial forces due to the instability of the flow to the inertial forces caused by velocity change, is defined as: where The St number considers the effects of both the pulsation period and mass flux comprehensively. Besides this factor, inlet and outlet vapor quality are also chosen as the designing factors because they determine the flow patterns combined with mass flux. The HTC ratio defined as the ratio of HTC in the pulsating flow (h pulsating ) to that in the continuous flow (h continuous ) is selected as the response. It can be written as follows: In this paper, the number of variables is 3, that is, n = 3. x 1 , x 2 , and x 3 represents the St number, inlet vapor quality x in , and outlet vapor quality x out . The experimental variables and varying ranges are given in Table 6. First the factors and varying ranges are the input and the experimental matrix is generated by the RMS design. Then taking the logarithm of both variables and responses, the experimental matrix can be regenerated. Eventually, taking exponentials for both sides of the RSM regression model, the HTC ratio correlation was achieved. According to the design by RSM, 160 of 264 experimental data were selected to obtain the HTC ratio. The experimental designing parameters and obtained HTC ratio were taken as the Napierian logarithm, and input for the matrix designing by RSM. The RSM design matrix in the form of the Napierian logarithm is shown in Table 6. The factors were tabulated in the middle three columns and the response (the ratio h*) was given in the far right column of Table A1, please see the Appendix A.

Analysis of Variance (ANOVA)
An ANOVA was performed to observe the significance of the model and designing parameters statistically. Table 7 shows the ANOVA results of the regression model for ln(h*). In ANOVA, the sum of squares was adopted to assess the square of deviation from the grand mean. Dividing the sum of squares by the degrees of freedom, one can obtain the mean squares. The F value was calculated using the ratio of the mean square for each term to the mean square for the residual. The value of F was greater than the F-table value or the p-value was less than 0.05, indicating that the model was significant or the corresponding variables had significant influences on the responses. As Table 7 shows, the main effects of ln(St), ln(x in ), and ln(x out ), as well as the interaction terms of ln(St) and ln(x in ) and ln(x in ) and ln(x out ) and the quadratic terms of (ln(St)) 2 and (ln(x in )) 2 were statistically significant because their p-values were smaller than 0.05. "Adeq Precision" measures the signal to noise ratio. An "Adeq Precision" greater than 4 is desirable. The ratio of 34.433 indicates an adequate signal. Therefore, this model can be used to predict the results.

RSM Regression Model for the HTC Ratio
After the ANOVA, the regression model developed by RSM is finally given as follows: Taking the exponent for both sides, the Equation (15) can be transformed as: where The normal probability plot for the residuals obtained from RSM for the regression model is displayed in Figure 5. It can be seen that all the residuals generally fell on or near the diagonal line, indicating that the errors satisfied normal distribution. In addition, almost all internally studentized residuals were close to zero, indicating that the regression model of the response had a better prediction precision.
The normal probability plot for the residuals obtained from RSM for the regression model is displayed in Figure 5. It can be seen that all the residuals generally fell on or near the diagonal line, indicating that the errors satisfied normal distribution. In addition, almost all internally studentized residuals were close to zero, indicating that the regression model of the response had a better prediction precision.  Figure 6 presents the comparisons between RSM predicted results and 246 experimental data of the HTC ratio. It can be seen that the RSM model tended to predict the experimental data accurately with a mean error of 0.76% and a standard deviation of 5.9%. Of the data 91.4% was contained inside a ± 10% error window while 100% of the data lay in a ±20% error window, indicating that the prediction of the RSM model of HTC ratio was satisfactory. Next, the RSM model of the HTC ratio and Wojtan's model in the continuous flow were combined to develop a model to predict the HTC in the two-phase pulsating flow.  Figure 6 presents the comparisons between RSM predicted results and 246 experimental data of the HTC ratio. It can be seen that the RSM model tended to predict the experimental data accurately with a mean error of 0.76% and a standard deviation of 5.9%. Of the data 91.4% was contained inside a ± 10% error window while 100% of the data lay in a ±20% error window, indicating that the prediction of the RSM model of HTC ratio was satisfactory. Next, the RSM model of the HTC ratio and Wojtan's model in the continuous flow were combined to develop a model to predict the HTC in the two-phase pulsating flow.

HTC model for R134a Two-Phase Pulsating Flow Boiling
The heat transfer characteristic in the R134a two-phase pulsating flow is much more complicated compared with the two-phase continuous flow. It can be affected by many parameters, such as the pulsation period, mass flow rate, vapor quality, and so on. After obtaining the HTC ratio model by RSM, the final prediction of HTC in the two-phase pulsating flow can be written as Wojtan's model multiplying by the RSM model (Equation (16)).
The comparison between the correlation predicted and experimental HTC in the two-phase pulsating flow is depicted in Figure 7. The results obtained with the new model shows a low mean error of −4.3% and a low standard deviation of 15.4%. Of the data 81.2% and 95.1% were contained

HTC model for R134a Two-Phase Pulsating Flow Boiling
The heat transfer characteristic in the R134a two-phase pulsating flow is much more complicated compared with the two-phase continuous flow. It can be affected by many parameters, such as the pulsation period, mass flow rate, vapor quality, and so on. After obtaining the HTC ratio model by RSM, the final prediction of HTC in the two-phase pulsating flow can be written as Wojtan's model multiplying by the RSM model (Equation (16)).
The comparison between the correlation predicted and experimental HTC in the two-phase pulsating flow is depicted in Figure 7. The results obtained with the new model shows a low mean error of −4.3% and a low standard deviation of 15.4%. Of the data 81.2% and 95.1% were contained inside the ±30% and ±50% error window, indicating that the correlation could predict the HTC in the two-phase pulsating flow well. It should be noted that a large dispersion between prediction and experimental values occurred at a high mass flow rate with a high inlet vapor quality. The error was mainly due to the overprediction of the annular zone in the two-phase continuous flow, as shown in Figure 4e.
two phase pulsating flow.

HTC model for R134a Two-Phase Pulsating Flow Boiling
The heat transfer characteristic in the R134a two-phase pulsating flow is much more complicated compared with the two-phase continuous flow. It can be affected by many parameters, such as the pulsation period, mass flow rate, vapor quality, and so on. After obtaining the HTC ratio model by RSM, the final prediction of HTC in the two-phase pulsating flow can be written as Wojtan's model multiplying by the RSM model (Equation (16)).
The comparison between the correlation predicted and experimental HTC in the two-phase pulsating flow is depicted in Figure 7. The results obtained with the new model shows a low mean error of −4.3% and a low standard deviation of 15.4%. Of the data 81.2% and 95.1% were contained inside the ±30% and ±50% error window, indicating that the correlation could predict the HTC in the two-phase pulsating flow well. It should be noted that a large dispersion between prediction and experimental values occurred at a high mass flow rate with a high inlet vapor quality. The error was mainly due to the overprediction of the annular zone in the two-phase continuous flow, as shown in Figure 4e.

Conclusions
In this study, a new model for predicting HTC in the R134a liquid-vapor two-phase pulsating flow was developed based on the heat transfer model in the two-phase continuous flow and response surface methodology. The conclusions can be summarized as follows: (1) Six existing models for heat transfer in the two-phase continuous flow were compared with our experimental data. Results show that the Wojtan's model had the smallest mean error and the most data lying in the ±50% error window. Therefore, the Wojtan's model was selected as a basis for predicting HTC in the R134a two-phase continuous flow. (2) RSM was carried out to rearrange the experimental data to obtain the regression model for the HTC ratio. ANOVA was performed to test the statistical significance of the model. The small error between RSM predicted data and experimental results shows that the model was satisfactory.