Next Article in Journal
Experimental Studies on Thermophysical and Electrical Properties of Graphene–Transformer Oil Nanofluid
Next Article in Special Issue
Influence of Hydrodynamic Conditions on Micromixing in Microreactors with Free Impinging Jets
Previous Article in Journal
The Influence of Inflow Swirl on Cavitating and Mixing Processes in a Venturi Tube
Previous Article in Special Issue
Temperature Uniformity in Cross-Flow Double-Layered Microchannel Heat Sinks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Condensation of an Azeotropic Mixture inside 2.5 mm ID Minitubes

Department of Industrial Engineering, University of Padova, via Venezia 1, 35131 Padova, Italy
*
Author to whom correspondence should be addressed.
Fluids 2020, 5(4), 171; https://doi.org/10.3390/fluids5040171
Submission received: 31 August 2020 / Revised: 24 September 2020 / Accepted: 30 September 2020 / Published: 2 October 2020
(This article belongs to the Special Issue Recent Advances in Single and Multiphase Flows in Microchannels)

Abstract

:
The ongoing miniaturization of air conditioning and refrigeration systems, in order to limit, as much as possible, the refrigerant charge, calls for smaller and smaller heat exchangers. Besides, the new environmental regulations are calling for new pure refrigerants or refrigerants mixtures with lower values of global warming potentials (GWPs). In this context, this paper analyzes the possible implementation of minitubes during condensation of the azeotropic mixture R513A. Two minitubes are tested: a smooth tube with an inner diameter of 2.5 mm, and a microfin tube with an inner diameter at the fin tip of 2.4 mm. The effects of vapor quality (varied in the range 0.10–0.99), of mass velocity (varied in the range 200–1000 kg m−2 s−1), and of saturation temperature (30 °C and 40 °C) on the heat transfer coefficient are investigated. The experimental results indicate that the heat transfer coefficient increases as both vapor quality and mass velocity increase, both in the case of the smooth tube and of the microfin tube, but the slope of the heat transfer coefficient trend respect to vapor quality is higher in the case of the microfin tube. The microfin tube shows, on average, heat transfer coefficients are 79% higher than those of the smooth tube under the same working conditions. Since R513A is a possible substitute of R134a, some experimental data during condensation heat transfer are also compared against those for R134a. Finally, the experimental results are compared against values estimated by empirical correlations available in the open literature.

Graphical Abstract

1. Introduction

R513A is an azeotropic mixture made of R1234yf and R134a (0.56/0.44 by mass). R1234yf is an HydroFluroOlefin with a Global Warming Potential lower than 1 [1], and it is placed in the A2L flammability class. R134a is among the most used HydroFluoroCarbons, not flammable (A1 flammability class), widely implemented in existing air conditioning systems, but it must be phased out in a near future due to its high GWP, which is 1430. Therefore, the presence of R134a makes the mixture not flammable, whereas the presence of R1234yf lowers the warming impact of the mixture. Coupled with small sized tubes, R513A could be a viable solution for the next generation of air conditioning equipment.
R32 and R134a local heat transfer coefficients were measured by Matkovic et al. [2]. Experimental tests were carried out at 40 °C saturation temperature inside a 0.96 mm ID circular minitube for mass velocities in the range 100–1200 kg m−2 s−1. The heat transfer coefficient increased with mass velocity and vapor quality; therefore, condensation was shear stress dominated, except some data at 100 and 200 kg m−2 s−1, which overlapped in the graph.
Zhang et al. [3] experimentally investigated the condensation of R22, R410A, and R407C inside two stainless steel tubes with inner diameters of 1.088 mm and 1.289 mm. Tests were run for mass velocities in the range 300–600 kg m−2 s−1, for saturation temperatures of 30 °C and 40 °C. The experimental results indicated that the condensation heat transfer coefficient increases with mass velocity and vapor quality, increasing more rapidly in the high vapor quality region.
More recently, Guo et al. [4] obtained heat transfer coefficients during the condensation of R1234ze(E), R290, R161, and R41 in a smooth horizontal tube with an inner diameter of 2 mm. Mass velocity was varied from 200 to 400 kg m−2 s−1, and saturation temperature from 35 °C to 40 °C. The experimental results showed that the heat transfer coefficient increases as vapor quality and mass velocity increase, and as saturation temperature decreases. Other works related to condensation in minichannels can be found in [5,6,7].
Adding microfins on the inner circumference of the tube has proven to enhance the two-phase heat transfer. Small-sized microfin tubes were initially studied for applications that involve carbon dioxide, due to its high working pressure [8,9,10]. Considering microfin tubes classified as minitubes, i.e., with diameters ranging between 1 and 3 mm [11], just a few studies investigated refrigerants’ two-phase heat transfer. Experimental measurements of the heat transfer coefficient were carried out by Diani et al. [12], who studied R1234yf and R1234ze(E) condensing inside a 2.4 mm ID microfin tube. Mass velocities were varied from 300 to 1000 kg m−2 s−1, for saturation temperatures of 30 °C and 40 °C. R1234ze(E) and R134a showed similar values of heat transfer coefficients, whereas R1234yf showed lower values.
Bashar et al. [13] investigated the condensation heat transfer of R134a flowing inside a microfin tube having an equivalent diameter of 2.17 mm and inside a smooth tube with an inner diameter of 2.14 mm. Data were measured for mass velocities from 50 to 300 kg m−2 s−1, and saturation temperatures from 20 °C to 30 °C. The heat transfer coefficient was higher in the case of the microfin tube for all the investigated mass fluxes and, considering the mass velocity 100 kg m−2 s−1, it was 2–5 times greater than that of the smooth tube.
Diani et al. [14] conducted experiments during R1234yf flow boiling inside a microfin tube with an outer diameter of 3.0 mm, for mass fluxes from 375 to 940 kg m−2 s−1, heat fluxes from 10 kW m−2 to 50 kW m−2, with a saturation temperature at the inlet of the test section of 30 °C. The experimental results showed that the flow boiling phenomenon is controlled by nucleate boiling, especially at low vapor qualities, and by convective boiling, especially at high vapor qualities.
Jige and Inoue [15] investigated R32 flow boiling inside microfin tubes with equivalent diameters of 2.1 mm, 2.6 mm, and 3.1 mm. Mass velocity was varied from 50 to 200 kg m−2 s−1 and heat flux from 5 to 40 kW m−2. The saturation temperature was 15 °C. The effect of tube diameter on the flow boiling heat transfer coefficient was more appreciable at higher vapor quality and lower mass velocity (the lower the tube diameter, the higher the heat transfer coefficient). Frictional pressure drops were highly affected by tube diameter: the lower the tube diameter, the higher the pressure drop.
Hirose et al. [16] studied R1234ze(E) condensation inside three different microfin tubes with the same outer diameter of 4 mm. The effects of fin parameters, i.e., of helix angle, fin height, and fin number were analyzed for mass velocities in the range 50–400 kg m−2 s−1, at a saturation temperature of 35 °C.
There are many candidates to replace R134a. Some of them are pure fluids, such as R1234ze(E) and R1234yf; others are azeotropic mixtures, such as R513A, R515A and R516A; others are zeotropic mixtures, such as R450A, R445A, R430A, R436A, and R456A. Heredia-Aricapa et al. [17] overviewed possible refrigerants mixtures as replacements of common HFCs like R134a, R404A, and R410A. Among the replacements of R134a, they identified R445A, R430A, R515A, R436A, R456A, and R515A. R436A and R430A are the alternatives with the greatest differences in the thermophysical properties, due to the presence of hydrocarbons in the mixtures. The comparison among fluids was given in terms of COP and cooling capacity in a vapor compression system.
Kedzierski et al. [18] carried out experimental tests during pool boiling on a reentrant cavity surface for R134a, R1234yf, R513A, and R450A. At constant wall superheat and pressure, with a saturation temperature of 277.6 K, the heat flux for R1234yf and R513A was, on average, 16% and 19% less than that for R134a, respectively, for heat fluxes between 20 kW m−2 and 110 kW m−2, whereas the heat flux for R450A was, on average, 57% less than that of R134a for heat fluxes between 30 kW m−2 and 110 kW m−2. Other works related to R513A can be found in Diani et al. [19,20], but no comparisons with R134a were reported.
Being a relatively new mixture, to date, no work about R513A condensation inside mini microfin tubes with diameters lower than 3 mm can be found in the open literature. Two minitubes are tested in this research: the first one is a smooth tube with an inner diameter of 2.5 mm, and the second one is a microfin tube with an inner diameter at the fin tip of 2.4 mm. This work is aimed at enlarging the experimental database of heat transfer coefficients during R513A condensation inside smooth and microfinned minitubes. The wide range of operative conditions will permit us to understand how each working parameter, i.e., vapor quality, mass velocity, and saturation temperature, affects the two-phase heat transfer both in the case of the smooth tube and of the microfin tube, highlighting the differences among the two tubes. Finally, the experimental condensation heat transfer coefficients will be compared against values estimated by empirical correlations.

2. Experimental Facility and Test Sections

The experimental facility is in the Laboratory of Heat Transfer in MicroGeometries at the Department of Industrial Engineering of University of Padova, Italy. Figure 1 shows a schematic of the liquid pumped circuit. A magnetically coupled gear pump was chosen to avoid any possible oil traces in the refrigerant. The pump circulates the refrigerant into a Coriolis effect mass flowmeter, which presents a reading accuracy of ±0.1%. The refrigerant enters an evaporator in subcooled conditions, and it exits in superheated conditions. The evaporator is a BPHE where refrigerant flows in countercurrent with hot water, supplied in the hot water loop. An electric boiler, regulated with a PID controller, guarantees an inlet water temperature of 60 °C. A tube-in-tube pre-condenser is aimed at partially condensing the superheated vapor, which flows in the inner tube. Cold water, supplied by a water-cooled stabilized chiller, flows in the annulus. The chiller permits to set the temperature of the water at the inlet of the pre-condenser. In addition, a valve permits variation in the water flow rate, and therefore it is possible to adjust the heat flow rate exchanged in the pre-condenser to set the refrigerant conditions at the inlet of the test section. The refrigerant partially condenses in the test section, and it is later fully condensed and subcooled in the post-condenser, which is another BPHE. In this case, the refrigerant flows in counter current with a mixture of water/glycol supplied by an air-cooled chiller. The refrigerant passes through a drier filter prior to re-enter into the pump. The saturation conditions are controlled by means of a damper, which is connected to the air pressure line.
As depicted in Figure 1, temperatures (measured by T-type thermocouples with an accuracy of ±0.05 K) and pressures of the refrigerant (measured by absolute pressure transducers with an accuracy of ±1950 Pa) are measured in several locations throughout the circuit, so it is possible to determine its thermodynamic state. The refrigerant conditions at the inlet of the test section, i.e., at the exit of the pre-condenser, derive from an energy balance at the pre-condenser, as will be explained in the next section. The water flow rate in the cold-water loop at the pre-condenser, necessary for the energy balance, is measured by a magnetic flowmeter, with a reading accuracy of ±0.25%. The water temperature difference at the pre-condenser is measured by a T-type thermopile, which has an accuracy of ±0.03 K.
Two different minitubes made of copper were tested. The first one is a smooth tube with an internal diameter of 2.5 mm and an outer diameter of 4.5 mm. The second one is a microfin tube with an inner diameter at the fin tip of 2.4 mm, an outer diameter of 3.0 mm, having 40 fins, and each fin has an height of 0.12 mm and an apex angle of 43°; the helix angle is 7°. Both the test section for the smooth tube and the one for the microfin tube are composed of the tested tube, which is horizontally located, around which a 1.9 mm ID smooth tube is wrapped. Cold water, supplied by a thermostatic bath and measured by a magnetic flowmeter with an accuracy of ±0.50% of the reading, flows in the wrapped tube, to remove the heat flow rate necessary for the condensation of the refrigerant flowing inside the horizontally located tested tube. Two curves, one at the inlet and another at the outlet, permit to host two T-type thermocouples (accuracy ±0.05 K), which measure the water inlet and outlet temperatures in the test section. On the external wall of the tested tubes, T-type thermocouples were glued to measure the wall temperature. Figure 2 shows a schematic, as well as the geometrical dimensions, of the test section for the smooth tube. The blue dots represent the thermocouples which measure the water temperatures, whereas the red dots represent the thermocouples attached on the external wall of the tested tube. The tested tubes with the wrapped tubes were finally inserted inside an aluminum housing with a U-shape. This housing was filled with an alloy of tin/lead, which was casted inside. This alloy is needed to have a good thermal contact between the wrapped tube, inside which the cold water flows, and the tested tube, inside which the refrigerant condenses. AF/Armaflex was used as thermal insulation around the test section to limit the heat losses through the external ambient. Pressure ports were soldered downstream and upstream of the test sections: the upstream pressure port is connected to an absolute pressure transducer (accuracy of ±1950 Pa), whereas both upstream and downstream pressure ports are connected to a differential pressure transducer (accuracy of ±25 Pa).

3. Data Analysis

The vapor quality at the inlet of the test section, hTS,in, can be extrapolated from the specific enthalpy at the exit of the pre-condenser, hPC,out, calculated as:
h T S , i n = h P C , o u t = h v s m ˙ w , P C · c p , w · Δ T w , P C m ˙ r e f
where hvs is the specific enthalpy of the refrigerant in superheated vapor conditions at the exit of the evaporator, which is calculated from the knowledge of its temperature and pressure. m ˙ w , p c is the mass flow rate of the water flowing in the pre-condenser, cp,w is the water specific heat, ∆Tw,pc is the water temperature difference between outlet and inlet of the pre-condenser, and m ˙ r e f is the refrigerant mass flow rate. The inlet vapor quality can be calculated from the knowledge of hTS,in as:
x i n = h T S , i n h L h V h L
where hV and hL are the specific enthalpies of saturated vapor and liquid, respectively, calculated from the inlet pressure. REFPROP 10 (Lemmon et al. [21]) was used to calculate all the refrigerant thermodynamic properties, since it implements updated equations for R513A.
In a similar way, the refrigerant specific enthalpy at the outlet of the test section, hTS,out, can be derived from an energy balance, as:
h T S , o u t = h T S , i n q w , T S m ˙ r e f = h T S , i n m ˙ w , T S · c p , w · ( t w , T S , o u t t w , T S , i n ) m ˙ r e f
where qw,TS is the heat flow rate calculated in the water side, and, in this case, m ˙ w , T S is the mass flow rate of the water flowing in the test section, tw,TS,out and tw,TS,in are the water temperature at the outlet and inlet of the test section, respectively. Once hTS,out is known, the outlet vapor quality is calculated as:
x o u t = h T S , o u t h L h V h L
where hV and hL are the specific enthalpy of saturated vapor and liquid, respectively, calculated at the outlet pressure. Experimental results, discussed in the next section, will be referred to the mean vapor quality, which is the arithmetic average between inlet vapor quality and outlet vapor quality.
The experimental tests permitted to calculate, both for the smooth tube and for the microfin tube, the heat transfer coefficient (HTC), as:
H T C = q w , T S A D · ( t ¯ s a t t ¯ w a l l )
where t ¯ s a t and t ¯ w a l l are the mean saturation and wall temperatures. AD is equal to the inner surface area in the case of the smooth tube, whereas it is equal to the inner surface area of an equivalent smooth with the inner diameter at the fin tip in the case of the microfin tube.
Following the procedure suggested by Kline and McKlintock [22], considering a level of confidence of 95%, the mean uncertainty on the heat transfer coefficient is ±4.2%, whereas on the mean vapor quality is ±0.027.

4. Experimental Results

4.1. Smooth Tube

4.1.1. Single-Phase Results

The test section developed for the smooth tube was first verified during single phase flow involving the forced convection of subcooled liquid. The main purposes of these tests were to check the heat flow rate calculated on the water side and the one on the refrigerant side and to verify if the measuring technique was suitable by calculating the single phase heat transfer coefficients and by comparing these results against values predicted by well-known correlations from the literature.
These tests were run with an inlet subcooling level between 17 K and 21 K, at a pressure of about 11.5 bar, for mass velocities from 400 to 1000 kg m−2 s−1. Mass velocity is defined as:
G = 4 · m ˙ r e f π · D i 2
where Di is the inner diameter. The heat transfer coefficient can be calculated as in Equation (5) but in this case, t ¯ r e f is used instead of t ¯ s a t . Two T-type thermocouples, downstream and upstream of the test section, were used to measure the inlet and outlet refrigerant temperatures.
The experimental results demonstrate that the heat flow rates calculated on the water side and on the refrigerant side are, on average, within ±4%. Figure 3 reports the heat transfer coefficients plotted against the mass velocity. As expected, the heat transfer coefficient increases as the mass velocity increases. Figure 3 also reports the comparison between the experimental values and the values predicted by the correlations of Dittus and Boelter [23] and of Petukhov and Popov [24]. The correlation of Petukhov and Popov [24] aptly predicts the experimental values, with a relative deviation of 3.9%, an absolute deviation of 4.0% and a standard deviation of 4.0%. The correlation of Dittus and Boelter [23] shows a relative deviation of −0.9%, an absolute deviation of 7.8% and a standard deviation of 9.1%.
Therefore, since the percentage difference between heat flow rate on the refrigerant side and the one on the water side is within ±4%, and since empirical correlations taken from the literature are able to estimate the experimental values of liquid-phase heat transfer coefficients, the test section for the 2.5 mm ID smooth tube is deemed to be verified.

4.1.2. Condensation Results

The operating conditions during the condensing tests for the 2.5 mm ID smooth tube are the following: mass velocity in the range 200–1000 kg m−2 s−1 and saturation temperature at the inlet of 30 °C and 40 °C. Experimental tests involve partial condensation inside the test tube, with a vapor quality difference between inlet and outlet from 0.08 to 0.40, depending on the mass velocity. Therefore, in order to investigate vapor qualities from 0.1 to 0.99, tests were run by varying the inlet vapor quality.
The collected heat transfer coefficients during condensation inside the smooth tube at a saturation temperature at the inlet of the test section of 40 °C are reported in Figure 4. Only data which involve an exchanged heat flow rate higher than 30 W are reported, in order to avoid data with high uncertainties. Data are plotted against the mean vapor quality in the test section at different mass velocities. A clear effect of both mass velocity and vapor quality is visible: the higher the mass velocity and the vapor quality, the higher the heat transfer coefficient, due to the diminishing liquid film thickness. According to the flow regime map developed by Doretti et al. [25] and reported by Cavallini et al. [26], all the experimental data fall in the ∆T-independent zone, i.e., annular flow, except data at G = 200 kg m−2 s−1 up to a mean vapor quality of 0.5, which fall into the ∆T-dependent zone, i.e., they are in stratified flow regime. Data at higher mass velocity shows a higher slope of the heat transfer coefficient trend, highlighting a higher convective effect as the mass velocity increases.
Figure 5 reports the data collected during condensation inside the smooth tube for a saturation temperature at the inlet of the test section of 30 °C. All these data fall in the ∆T-independent zone and, as a result, the heat transfer coefficient increases as both the vapor quality and mass velocity increase. The comparison with Figure 4 reveals the effect of the saturation temperature. A higher saturation temperature implies a higher vapor density; thus, at constant mass velocity, it also implies a lower vapor velocity, with consequent lower convective effects. This aspect can be observed in the heat transfer coefficients reported in Figure 4 and Figure 5: the heat transfer coefficient tends to decrease as the saturation temperature increases, due to lower convective effects.

4.2. Microfin Tube

4.2.1. Single-Phase Results

In a similar way to the validation procedure depicted in Section 4.1.1, single-phase tests were carried out during R1234yf liquid flow to validate the test section for the 2.4 mm ID microfin tube, and the experimental results revealed that the heat flow rate calculated on the water side was within ±1.5 W the one on the refrigerant side, and that the experimental single-phase heat transfer coefficients were aptly estimated by empirical correlations from the literature. Therefore, the test section for the microfin tube was deemed verified. More details can be found in Diani et al. [12].

4.2.2. Condensation Results

The operating conditions during the condensing tests for the 2.4 mm ID microfin tube are the following: mass velocity in the range 200–1000 kg m−2 s−1, saturation temperature at the inlet of 30 °C and 40 °C. In this case, the mass velocity is calculated as reported in Equation (6) but the inner diameter at the fin tip is used as Di.
Figure 6 shows the condensation heat transfer coefficients plotted against the mean vapor quality for the microfin tube with a saturation temperature at the inlet of the test section of 40 °C. Again, only data which involve an exchanged heat flow rate higher than 30 W are reported. According to the flow regime map developed by Doretti et al. [25] and reported by Cavallini et al. [26], all the experimental data fall in the annular flow regime. Therefore, the higher the vapor quality and mass velocity, the higher the heat transfer coefficient will be. Compared to the results during condensation inside the smooth tube reported in Figure 4, the condensation heat transfer coefficients for the microfin tube are higher, and the comparison reveals that data at mass velocities of 200, 300, and 400 kg m−2 s−1 are particularly enhanced. Indeed, the slope of the heat transfer coefficients with mean vapor quality is now much higher for these low mass velocities: taking for instance G = 400 kg m−2 s−1, the heat transfer coefficient passes from approximately 2000 W m−2 K−1 at xmean = 0.3 to a value of about 5000 W m−2 K−1 at xmean = 0.8 for the smooth tube, i.e., it is 2.5 times higher, whereas it passes from approximately 3500 W m−2 K−1 at xmean = 0.3 to a value of about 12500 W m−2 K−1 at xmean = 0.8 for the microfin tube, i.e., it is 3.5 times higher.
Results with an imposed saturation temperature at the inlet of the test section of 30 °C for the microfin tube are reported in Figure 7. Again, all the experimental data fall in the annular flow regime, and the heat transfer coefficient increases with mass velocity and vapor quality. The comparison with Figure 6 reveals the effect of saturation temperature: as stated for the smooth tube, the higher the saturation temperature, the lower the heat transfer coefficient, due to lower convective effects that result from the lower vapor velocity.

4.3. Comparison among the Tubes

The heat transfer coefficients of smooth tube and microfin tube are here compared under the same working conditions. To better explain the comparison, the parameter called the Enhancement Factor (EF) is introduced:
E F = H T C m i c r o f i n H T C s m o o t h
where HTCmicrofin and HTCsmooth represent the heat transfer coefficient of the microfin tube and of the smooth tube, respectively, under the same operative conditions.
The Enhancement Factor is plotted versus the mean vapor quality in Figure 8, considering a saturation temperature of 40 °C. Adding microfins on the internal area of the microfin tube improves the thermal behavior during condensation; indeed, the EF is always higher than 1. In addition, the benefits of the microfins increase as mass velocity decreases, i.e., the highest EF is shown at the lowest mass velocities. At these lowest mass velocities, the EF increases along with the increasing vapor quality. Similar conclusions are made at tsat = 30 °C.
The area enhancement ratio, i.e., the ratio between the total inner surface area of the microfin tube and the inner surface of a smooth tube having the fin tip diameter as inner diameter, is 1.89. It is worth highlighting that, at G = 300 and 400 kg m−2 s−1 for mean vapor qualities higher than 0.4, the EF is higher than the mere increase of the heat transfer area, and thus the microfins have a positive effect on the heat transfer especially at low mass velocities, increasing the turbulence of the liquid film. Considering all the experimental data, the average Enhancement Factor is 1.79.

4.4. Comparison among Refrigerants

Being one of the possible direct drop in replacements of R134a in the near future, the performance of R513A is compared against that of R134a. Heat transfer coefficients during R134a condensation inside a microfin tube with an inner diameter at the fin tip of 2.4 mm are borrowed by Diani et al. [12]. The comparison can start by comparing the different thermophysical and transport properties of the two fluids. Table 1 reports the main properties of R134a and R513A.
Among the listed properties, R134a shows a lower vapor density than R513A, i.e., at constant mass velocity the vapor velocity is higher, with consequent higher convective effects. Furthermore, the liquid thermal conductivity of R134a is higher than that of R513A, with consequently better thermal characteristics of the liquid film thickness. Considering the pressure drop, the higher vapor velocity at constant mass velocity for R134a will also be reflected in the hydraulic behavior.
Figure 9 shows a comparison between R134a and R513A heat transfer coefficient during condensation inside the microfin tube. As it appears from the figure, the heat transfer coefficients of R134a are higher than those of R513A, especially in the high vapor quality region, due to the aforementioned motivations. Similar conclusions can be drawn at a saturation temperature of 40 °C.
Frictional pressure gradients can be evaluated from the total measured pressure drop by subtracting the momentum pressure gradient. The model of Rouhani and Axelsson [27] is applied to estimate the void fraction, needed for the calculation of the momentum term. Figure 10 shows the comparison of frictional pressure gradients at tsat = 30 °C. The frictional pressure gradient increases as the mass velocity increases and as the mean vapor quality increases, until a maximum value is reached at mean vapor qualities in the range 0.8–0.9, and then it slightly decreases. The comparison between the two fluids reveals that R134a shows a slightly higher frictional pressure gradient, especially at high vapor qualities, due to its lower vapor density with the consequent higher vapor velocity at constant mass velocity. A similar conclusion can be drawn at a saturation temperature of 40 °C.

5. Empirical Modelling

5.1. Smooth Tube

This section provides a comparison between the experimental condensation heat transfer coefficients collected for the smooth tube and the values estimated by empirical correlations available in the literature. Figure 11 shows a comparison between the experimental heat transfer coefficients for the smooth tube and the values estimated by the correlations developed by Dobson and Chato [28], Shah [29], Cavallini et al. [30], and Dorao and Fernandino [31]. The correlation of Dobson and Chato [28] is able to estimate the experimental values with a relative deviation of 29.6%, an absolute deviation of 29.6%, and a standard deviation of 16.2%; the correlation of Shah [29] with a relative deviation of 16.3%, an absolute deviation of 18.0%, and a standard deviation of 18.4%; the correlation of Cavallini et al. [30] with a relative deviation of 5.9%, an absolute deviation of 9.7%, and a standard deviation of 13.6%; the correlation of Dorao and Fernandino [31] with a relative deviation of 2.5%, an absolute deviation of 15.3%, and a standard deviation of 22.8%. Therefore, the comparison reveals that the correlation of Cavallini et al. [30], even if it was developed for internal diameters larger than 3 mm, as well as that by Dorao and Fernandino [31], are the ones that better predict the experimental values.

5.2. Microfin Tube

This section shows the comparison between the condensation heat transfer coefficients for the microfin tube against the values predicted by the correlation developed by Cavallini et al. [26]. A comparison between experimental and calculated condensation heat transfer coefficients is reported in Figure 12. The correlation can predict the experimental values with a relative, absolute, and standard deviation of 3.5%, 21.2%, and 31.1%, respectively. This correlation was developed from data obtained during condensation, for mass velocity in the range 80–890 kg m−2 s−1, saturation temperature in the range −15 °C–70 °C, inside microfin tubes with inner diameters at the fin tip from 5.9 mm and 14.2 mm. This aspect may partly explain the reason why the correlation tends to underestimate the experimental values at the lowest mass velocities in the high vapor quality zone, whose trend is peculiar of small-sized microfin tubes.

6. Conclusions

The paper presented experimental data during condensation of R513A inside a smooth minitube with an inner diameter of 2.5 mm and inside a mini microfin tube with an inner diameter at the fin tip of 2.4 mm. The following working conditions were investigated: mass velocities in the range 200–1000 kg m−2 s−1 and saturation temperatures at the inlet of the test section of 40 °C and 30 °C.
The test section for the smooth tube was verified during single-phase tests, involving liquid forced convection: once the heat flow rate calculated on the water side and the one on the refrigerant side were verified to be within ±4%, and once the experimental heat transfer coefficients were found to be well estimated by correlations from the literature, the new test section was deemed verified.
The condensation tests revealed that the heat transfer coefficient increases as both vapor quality and mass velocity increases, both for the smooth tube and for the microfin tube. However, the increase of the heat transfer coefficient with vapor quality is more evident with the microfin tube at mass velocities up to 400 kg m−2 s−1, highlighting its better convective properties, in particular at the lowest mass velocities. Concerning the saturation temperature, the higher the saturation temperature, the lower the heat transfer coefficient, for both the smooth tube and the microfin tube.
The average Enhancement Factor, i.e., the ratio between the heat transfer coefficient of the smooth tube and the one of the microfin tube under the same operative conditions, is 1.79, with the highest values shown at G = 300 kg m−2 s−1 at high vapor quality.
Experimental R513A heat transfer coefficients and frictional pressure drops obtained for the microfin tube were compared against the counterpart with R134a. R134a shows higher values of condensation heat transfer coefficient, particularly at high vapor qualities, but also slightly higher values of frictional pressure gradients.
The comparison against values predicted by empirical correlations revealed that the correlation of Cavallini et al. [30] and Dorao and Fernandino [31] are the most suitable to predict the condensation heat transfer coefficient inside the smooth tube, and the correlation of Cavallini et al. [26] is suitable to predict the condensation values for the microfin tube.

Author Contributions

Conceptualization, A.D.; validation, A.D., investigation, A.D.; writing—original draft preparation, A.D.; visualization, A.D.; supervision, L.R.; funding acquisition, L.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MIUR through PRIN Project 2017F7KZWS_005 and by Università degli Studi di Padova through Project CPDA107382 and Project BIRD172935/17.

Acknowledgments

Jacopo Carnio is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Aarea [m2]
cspecific heat [J kg−1 K−1]
Ddiameter [m]
Gmass velocity [kg m−2 s−1]
hspecific enthalpy [J kg−1]
HTCheat transfer coefficient [W m−2 K−1]
Llength [m]
m ˙ mass flow rate [kg s−1]
ppressure [bar]
qheat flow rate [W]
PrPrandtl number [-]
ttemperature [°C]
xthermodynamic vapor quality [-]
Greek symbols
ppressure drop [Pa]
Ttemperature difference [K]
λthermal conductivity [W m−1 K−1]
μdynamic viscosity [Pa s]
ρdensity [kg m−3]
σsurface tension [N m−1]
Subscripts
ffrictional
iinner
inat the inlet
Lof the saturated liquid
outat the outlet
PCat the pre-condenser
rreduced
refof the refrigerant
satsaturated
TSin the test section
Vof the saturated vapor
vsof the superheated vapor
wof the water
wallat the wall

References

  1. Hodnebrog, O.; Etminiam, M.; Fuglestvedt, J.S.; Marston, G.; Myhre, G.; Nielsen, C.J.; Shine, K.P.; Wallington, T.J. Global warming potentials and radiative efficiencies of halocarbons and related compounds: A comprehensive review. Rev. Geophys. 2013, 51, 300–378. [Google Scholar] [CrossRef] [Green Version]
  2. Matkovic, M.; Cavallini, A.; Del Col, D.; Rossetto, L. Experimental study on condensation heat transfer inside a single circular minichannel. Int. J. Heat Mass Transf. 2009, 52, 2311–2323. [Google Scholar] [CrossRef]
  3. Zhang, H.Y.; Li, M.J.; Liu, N.; Wang, B.X. Experimental investigation of condensation heat transfer and pressure drop of R22, R410A and R407C in mini.tubes. Int. J. Heat Mass Transf. 2012, 55, 3522–3532. [Google Scholar] [CrossRef]
  4. Guo, Q.; Li, M.; Gu, H. Condensation heat transfer characteristics of low-GWP refrigerants in a smooth horizontal mini tube. Int. J. Heat Mass Transf. 2018, 126, 26–38. [Google Scholar] [CrossRef]
  5. Azzolin, M.; Bortolin, S.; Del Col, D. Convective condensation at low mass flux: Effect of turbulence and tube orientation on the heat transfer. Int. J. Heat Mass Transf. 2019, 144, 118646. [Google Scholar] [CrossRef]
  6. Azzolin, M.; Berto, A.; Bortolin, S.; Moro, L.; Del Col, D. Condensation of ternary low GWP zeotropic mixtures inside channels. Int. J. Refrig. 2019, 103, 77–90. [Google Scholar] [CrossRef]
  7. Murphy, D.L.; Macdonald, M.P.; Mahvi, A.J.; Garimella, S. Condensation of propane in vertical minichannels. Int. J. Heat Mass Transf. 2019, 137, 1154–1166. [Google Scholar] [CrossRef]
  8. Gao, L.; Honda, T.; Koyama, S. Experiments on flow boiling heat transfer of almost pure CO2 and CO2-oil mixtures in horizontal smooth and microfin tube. HVAC&R 2007, 13, 415–425. [Google Scholar]
  9. Kim, Y.J.; Cho, J.M.; Kim, M.S. Experimental study on the evaporative heat transfer and pressure drop of CO2 flowing upward in vertical smooth and micro-fin tubes with the diameter of 5 mm. Int. J. Refrig. 2008, 31, 771–779. [Google Scholar] [CrossRef]
  10. Dang, C.; Haraguchi, N.; Hihara, E. Flow boiling heat transfer of carbon dioxide inside a small-sized microfin tube. Int. J. Refrig. 2010, 33, 655–663. [Google Scholar] [CrossRef]
  11. Qu, W.; Mudawar, I. Measurement and correlation of critical heat flux in two-phase micro.channel heat sinks. Int. J. Heat Mass Transf. 2004, 47, 2045–2059. [Google Scholar] [CrossRef]
  12. Diani, A.; Campanale, M.; Cavallini, A.; Rossetto, L. Low GWP refrigerants condensation inside a 2.4 mm ID microfin tube. Int. J. Refrig. 2018, 86, 312–321. [Google Scholar] [CrossRef]
  13. Bashar, M.K.; Nakamura, K.; Kariya, K.; Miyara, A. Experimental study of condensation heat transfer and pressure drops inside a small diameter microfin and smooth tube at low mass flux condition. Appl. Sci. 2018, 8, 2146. [Google Scholar] [CrossRef] [Green Version]
  14. Diani, A.; Cavallini, A.; Rossetto, L. R1234yf flow boiling heat transfer inside a 2.4 mm microfin tube. Heat Transf. Eng. 2017, 38, 303–312. [Google Scholar] [CrossRef]
  15. Jige, D.; Inoue, N. Flow boiling heat transfer and pressure drop of R32 inside 2.1 mm, 2.6 mm and 3.1 mm microfin tubes. Int. J. Heat Mass Transf. 2019, 123, 566–573. [Google Scholar] [CrossRef]
  16. Hirose, M.; Jige, D.; Inore, N. Optimum fin geometries on condensation heat transfer and pressure drop of R1234ze(E) in 4 mm outside diameter horizontal microfin tubes. Sci. Technol. Built Environ. 2019, 25, 1271–1280. [Google Scholar] [CrossRef]
  17. Heredia-Aricapa, Y.; Belman-Flores, J.M.; Mota-Babiloni, A.; Serrano-Arellano, J.; Garcìa-Pablon, J.J. Overview of low GWP mixtures for replacement of HFC refrigerants: R134a, R404A and R410A. Int. J. Refrig. 2020, 111, 113–123. [Google Scholar] [CrossRef]
  18. Kedzierski, M.A.; Lin, L.; Kang, D. Pool boiling of low-global warming potential replacements for R134a on a reentrant cavity surface. J. Heat Transf. 2018, 140, 121502. [Google Scholar] [CrossRef]
  19. Diani, A.; Rossetto, L. R513A flow boiling heat transfer inside horizontal smooth tube and microfin tube. Int. J. Refrig. 2019, 107, 301–314. [Google Scholar] [CrossRef]
  20. Diani, A.; Rossetto, L. Characteristics of R513A evaporation heat transfer inside small-diameter smooth and microfin tubes. Int. J. Heat Mass Transf. 2020, 162, 120402. [Google Scholar] [CrossRef]
  21. Lemmon, E.W.; Huber, M.L.; McLinden, M.O. NIST Standard Reference Database 23, NIST Reference Fluid Thermodynamic and Transport Properties, REFPROP, Version 10.0, Standard Reference Data Program; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2018.
  22. Kline, S.J.; McKlintock, F.A. Describing the uncertainties ins single-sample experiments. Mech. Eng. 1953, 75, 3–8. [Google Scholar]
  23. Dittus, F.W.; Boelter, L.M.K. Heat transfer in automotive radiators of the annular type. Univ. Calif. Publ. Eng. 1930, 2, 443–461. [Google Scholar]
  24. Petukhov, B.S.; Popov, V.N. Theoretical calculation of heat exchange and frictional resistance in turbulent flow in tubes of an incompressible fluid with variable physical-properties. High Temper. 1963, 1, 69–83. [Google Scholar]
  25. Doretti, L.; Fantini, F.; Zilio, C. Flow patterns during condensation of three refrigerants: Microfin tube vs. smooth tube. Proceeding of the IIR International Conference on Thermophysical Properties and Transfer Processes pf Refrigerant, Vicenza, Italy, 31 August–2 September 2005. [Google Scholar]
  26. Cavalini, A.; Del Col, D.; Mancin, S.; Rossetto, L. Condensation of pure and near-azeotropic refrigerants in microfin tubes: A new computational procedure. Int. J. Refrig. 2009, 32, 162–174. [Google Scholar] [CrossRef]
  27. Rouhani, S.Z.; Axelsson, E. Calculation of void fraction in the subcooled and boiling region. Int. J. Heat Mass Transf. 1970, 13, 383–393. [Google Scholar] [CrossRef]
  28. Dobson, M.K.; Chato, J.C. Condensation heat transfer inside smooth horizontal tubes. J. Heat Transf. ASME 1998, 120, 193–213. [Google Scholar] [CrossRef]
  29. Shah, M.M. An improved and extended general correlation for heat transfer during condensation in plain tubes. HVAC&R Res. 2009, 15, 889–913. [Google Scholar]
  30. Cavallini, A.; Del Col, D.; Doretti, L.; Matkovic, M.; Rossetto, L.; 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]
  31. Dorao, C.A.; Fernandino, M. On the heat transfer deterioration during condensation of binary mixtures. Appl. Phys. Lett. 2019, 114, 171902. [Google Scholar] [CrossRef]
Figure 1. Schematic of the test rig.
Figure 1. Schematic of the test rig.
Fluids 05 00171 g001
Figure 2. Schematic of the test section for the 2.5 mm ID smooth tube.
Figure 2. Schematic of the test section for the 2.5 mm ID smooth tube.
Fluids 05 00171 g002
Figure 3. Single-phase heat transfer coefficient for the smooth tube vs. mass velocity. Experimental data and estimations by Dittus and Boelter [23] and by Petukhov and Popov [24].
Figure 3. Single-phase heat transfer coefficient for the smooth tube vs. mass velocity. Experimental data and estimations by Dittus and Boelter [23] and by Petukhov and Popov [24].
Fluids 05 00171 g003
Figure 4. Heat transfer coefficient during condensation for the smooth tube vs. mean vapor quality at tsat = 40 °C. G expressed in [kg m−2 s−1].
Figure 4. Heat transfer coefficient during condensation for the smooth tube vs. mean vapor quality at tsat = 40 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g004
Figure 5. Heat transfer coefficient during condensation for the smooth tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Figure 5. Heat transfer coefficient during condensation for the smooth tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g005
Figure 6. Heat transfer coefficient during condensation for the microfin tube vs. mean vapor quality at tsat = 40 °C. G expressed in [kg m−2 s−1].
Figure 6. Heat transfer coefficient during condensation for the microfin tube vs. mean vapor quality at tsat = 40 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g006
Figure 7. Heat transfer coefficient during condensation for the microfin tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Figure 7. Heat transfer coefficient during condensation for the microfin tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g007
Figure 8. Enhancement Factor (EF) vs. mean vapor quality at tsat = 40 °C. G expressed in [kg m−2 s−1].
Figure 8. Enhancement Factor (EF) vs. mean vapor quality at tsat = 40 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g008
Figure 9. R134a and R513A heat transfer coefficient during condensation for the microfin tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Figure 9. R134a and R513A heat transfer coefficient during condensation for the microfin tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g009
Figure 10. R134a and R513A frictional pressure gradient for the microfin tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Figure 10. R134a and R513A frictional pressure gradient for the microfin tube vs. mean vapor quality at tsat = 30 °C. G expressed in [kg m−2 s−1].
Fluids 05 00171 g010
Figure 11. Experimental and calculated heat transfer coefficients for the 2.5 mm ID smooth tube. Models developed by Dobson and Chato [28], Shah [29], Cavallini et al. [30], and Dorao and Fernandino [31].
Figure 11. Experimental and calculated heat transfer coefficients for the 2.5 mm ID smooth tube. Models developed by Dobson and Chato [28], Shah [29], Cavallini et al. [30], and Dorao and Fernandino [31].
Fluids 05 00171 g011
Figure 12. Experimental and calculated heat transfer coefficients for the 2.4 mm ID microfin tube. Model developed by Cavallini et al. [26].
Figure 12. Experimental and calculated heat transfer coefficients for the 2.4 mm ID microfin tube. Model developed by Cavallini et al. [26].
Fluids 05 00171 g012
Table 1. Thermophysical and transport properties of R134a and R513A. Values from REFPROP 10 (Lemmon et al. [21]).
Table 1. Thermophysical and transport properties of R134a and R513A. Values from REFPROP 10 (Lemmon et al. [21]).
PropertiesR134aR513A
tsat [°C]3030
psat [bar]7.7028.211
ρL [kg m−3]11881115
ρV [kg m−3]37.5443.50
λL [W m−1 K−1]0.0790.068
λV [W m−1 K−1]0.0140.015
μL [μPa s]183156
μV [μPa s]11.911.8
cp,L [J kg−1 K−1]14471436
cp,V [J kg−1 K−1]10661092
PrL [-]3.353.29
PrV [-]0.8850.888
σ [N m−1]0.00740.0061
pr [-]0.1900.225

Share and Cite

MDPI and ACS Style

Diani, A.; Rossetto, L. Condensation of an Azeotropic Mixture inside 2.5 mm ID Minitubes. Fluids 2020, 5, 171. https://doi.org/10.3390/fluids5040171

AMA Style

Diani A, Rossetto L. Condensation of an Azeotropic Mixture inside 2.5 mm ID Minitubes. Fluids. 2020; 5(4):171. https://doi.org/10.3390/fluids5040171

Chicago/Turabian Style

Diani, Andrea, and Luisa Rossetto. 2020. "Condensation of an Azeotropic Mixture inside 2.5 mm ID Minitubes" Fluids 5, no. 4: 171. https://doi.org/10.3390/fluids5040171

APA Style

Diani, A., & Rossetto, L. (2020). Condensation of an Azeotropic Mixture inside 2.5 mm ID Minitubes. Fluids, 5(4), 171. https://doi.org/10.3390/fluids5040171

Article Metrics

Back to TopTop