Next Article in Journal
Design of Uninterruptible Power Supply Inverters for Different Modulation Techniques Using Pareto Front for Cost and Efficiency Optimization
Next Article in Special Issue
Review of Oil–Water Flow Characteristics of Emptying by Water Displacing Oil in Mobile Pipelines
Previous Article in Journal
Experimental Research on Performance Comparison of Compressed Air Engine under Different Operation Modes
Previous Article in Special Issue
Numerical Modeling of Laser Heating and Evaporation of a Single Droplet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Study on Boiling Heat Transfer Characteristics and Performance of Low-Temperature Heating System of R744 and Its Azeotropic Refrigerant

1
School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
2
School of Energy and Power Engineering, Northeast Electric Power University, Jilin 132000, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(3), 1313; https://doi.org/10.3390/en16031313
Submission received: 6 December 2022 / Revised: 10 January 2023 / Accepted: 19 January 2023 / Published: 26 January 2023

Abstract

:
R744 is the most competitive and ideal natural refrigerant when flammability and toxicity are strictly limited. However, there are still some problems when it is applied to a heating system. For example, the discharge pressure of the system exceeds 10 MPa, it increases the cost of the system, and the cycle efficiency is also low. To solve these problems, this paper proposes to replace R744 by mixing R744 and ethane at a ratio of (77.6/22.4) to form an azeotropic refrigerant. At present, there is little research on R744 azeotropic refrigerant. Therefore, this paper first establishes the CFD model and compiles the UDF program to focus on flow boiling heat transfer characteristics, and then, it analyzes the performance of R744 and its azeotropic refrigerant in a low-temperature heating system. The results show that the heat transfer coefficient of R744 and its azeotropic refrigerant decreases with an increase in mass flux and increases with an increase in heat flux and saturation temperature; the heat transfer coefficient of azeotropic refrigerant is greater than R744; and there is no dryness under the same conditions. Under a given operating condition, there is a critical point that makes the performance of azeotropic refrigerant better than R744, and this critical point is related to the outlet temperature of a gas cooler, and the system discharge temperature of azeotropic refrigerant is significantly lower than that of R744. In conclusion, azeotropic refrigerant has certain advantages in heat transfer and system performance compared with R744, which will also play an important role in promoting the replacement of refrigerant in the future.

1. Introduction

At present, most of the refrigerant used in refrigeration and heating products is synthetic refrigerant, which has greenhouse effect potential (ODP) and ozone layer destruction potential (GWP), which have serious impacts on the environment. The natural working medium R744 receives wide concern from the industry thanks to its good environmental characteristics; however, there are still some problems when it is applied to heating systems. For example, the discharge pressure of the system exceeds 10 MPa, it increases the cost of the system, and the cycle efficiency is also low. To solve these problems, many scholars mixed R744 with other refrigerants. The research on refrigerant involves mainly the mechanism of flow boiling and the performance of the whole system. There are many experiments but few simulations on the mechanism of flow boiling in tubes with R744, and the research on R744 azeotropic refrigerant has been even more rarely reported.
In terms of experiments, Grauso et al. [1] studied the flow boiling heat transfer characteristics of R744/R290 with mass ratios of 83.2/16.8 and 70/30, respectively, through a horizontal smooth tube. The result showed that the flow boiling heat transfer characteristics of R744/R290 were evaluated. Afroz et al. [2] studied the heat transfer and pressure drop characteristics of binary mixtures (mass ratios: 39/61 and 21/79, respectively) in a horizontal smooth tube. The result showed that the heat transfer performance of R744/dimethyl ether was better than that of pure R744, and with an increase in the R744 mass ratio, the calculation formula was modified according to the experimental results. Onaka et al. [3] studied the evaporation heat transfer coefficients of pure dimethyl ether and a mixture of R744 and dimethyl ether (two mass ratios: 10/90 and 25/75, respectively) in a horizontal smooth tube with an inner diameter of 4.2 mm. The results showed that the heat transfer coefficient can be reduced by 20% for 10/90 and 48% for 25/75 by adding R744. Zhu et al. [4] studied the evaporation heat transfer characteristics of R744/R290 (75/25, 50/50, and 25/75) in a horizontal smooth microchannel pipe (inner diameter of 2 mm). The results showed that the heat transfer coefficient of the mixture is between R744 and R290, which is directly proportional to the ratio of R744 and inversely proportional to the evaporation temperature. Yun et al. [5] studied the flow boiling heat transfer of R744 in a horizontal tube with an inner diameter of 6 mm. The results showed that the heat transfer coefficient of R744 is on average 47% higher than that of R134a under the same conditions. Zhao et al. [6] studied the boiling heat transfer of R744 at low temperatures in a horizontal tube with an inner diameter of 4.57 mm. The results showed that the boiling heat transfer coefficient of R744 is much lower at low temperatures but that it increases with vapor quality. Oh et al. [7] studied the flow boiling heat transfer and pressure drop characteristics of R744 in a horizontal tube with an inner diameter of 4.57 mm. The results showed that the boiling pressure drop of R744 increased with an increase in the mass flux and decreased with an increase in the saturation temperature. Yun et al. [8] studied the convective boiling heat transfer coefficients and dry-out phenomena of R744 in rectangular microchannels whose hydraulic diameters ranged from 1.08 to 1.54 mm. The result showed that the hydraulic diameter decreases from 1.54 to 1.27 mm and from 1.27 to 1.08 mm at a heat flux of 15 kW/m2 and a mass flux of 300 kg/m2·s, respectively, and that the heat transfer coefficients increase by 5% and 31%, respectively. Ozawa et al. [9] studied the flow pattern and boiling heat transfer of R744 in horizontal small-bore tubes. The result showed that the low surface tension and latent heat also have significant influences on two-phase flow patterns and heat transfer. Huai et al. [10] studied the flow boiling characteristics of R744 in multiport minichannels in an experiment. The result showed that the pressure drop along the test section is very small, and the two-phase R744 flow exhibits a higher heat transfer coefficient than that of the single-phase liquid or vapor flow. Jiang et al. [11] studied the characteristics of heat transfer for the R744 flow boiling at a low temperature in a minichannel. The result showed that the heat transfer coefficient increases with an increase in the mass flux rate but decreases with an increase in the tubes’ inner diameter and saturation temperature, and the boiling heat transfer of R744 has a greater effect on nucleate boiling and highly depends on the heat flux. Zhang et al. [12] studied the flow boiling heat transfer characteristics of R744 in a horizontal minitube in an experiment. The result showed that an increase in the heat flux has a significant effect on enhancing the nucleate boiling heat transfer, while it can also speed up the dry-out process and decrease the starting vapor quality of dry-out. Kim et al. [13] studied the evaporative heat transfer and pressure drop of R744 flowing upward in vertical smooth and microfin tubes with a diameter of 5 mm in an experiment. The result showed that the average evaporation heat transfer coefficients for the microfin tubes were approximately 111–207% higher than those for the smooth tube under the same test conditions, and the PF was increased from 106% to 123%.
Regarding the simulation of R744, there are many studies on supercritical flow heat transfer, but few on boiling heat transfer. Wei et al. [14] used FLUENT software to simulate the flow boiling heat transfer of R744/R290 (50/50) in smooth tubes and microfin tubes with an inner diameter of 4 mm and a length of 1.4 m. The results showed that the convective heat transfer coefficient of microfin tubes is higher than that of smooth tubes. Wu et al. [15] studied the flow boiling heat transfer characteristics of R744/R290 in a horizontal smooth tube with an inner diameter of 2 mm. The initial evaporation temperature range was 0–10 °C, the heat flux was 5–15 kW/m2, and the mass flux rate was 200–400 kg/(m2·s). The results showed that the heat transfer coefficient of the mixture was between R744 and R290 and decreased with an increase in the R290 mass ratio, which was directly proportional to the initial evaporation temperature, heat flux, and mass flux rate, and the critical dryness was also proportional to the mass flux rate of R290.
On system performance, Kim et al. [16] carried out theoretical and experimental studies on the application of R744/R134a and R744/R290 in an auto-cascade refrigeration system. The results showed that with an increase in the R744 content, the COP increases and the system high pressure decreases. Sarkar et al. [17] proposed that R744/R600 and R744/R600a should be used for medium- and high-temperature heat pumps and that they should heat and cool at the same time. The results showed that R744/R600a instead of R114 for high-temperature heat pump was a good choice. Niu et al. [18] tested the performance of R744/R290 (71/29 mole fraction) for cascade refrigeration. The results showed that the COP and the cooling capacity were higher than those in the R13 system. Onaka et al. [19] analyzed the application of R744/RE170 in a heat pump water heater system. The results showed that the COP of the system was higher than that of pure R744. In addition, the research on the R744 transcritical cycle system is increasing day by day, and the research focuses mainly on the cooling performance of the R744 transcritical cycle, the R744 compressor, the expander [20,21,22,23,24,25], and the supercritical R744 [26,27,28,29,30,31,32].
In terms of the heat transfer phenomenon, Akram et al. [33] studied the results of the double-diffusion convection of Sisko nanofluids creeping along an asymmetric channel with an inclined magnetic field. The results showed that with an increase in the Dufour and thermophoresis variables, the concentration, temperature, and nanoparticle fraction significantly increase. Khan et al. [34] studied the influence factors on the creeping pumping of fourth-stage fluids in nonuniform channels. The results showed that as the influence of Brownian motion increased, the density of nanoparticles increased, which led to an increase in the fraction of nanoparticles. Athar et al. [35] studied the influencing factors of the creeping propulsion of oldroyd-4 constant nanofluids in asymmetric channels. The results showed that the temperature increased with an increase in Brownian motion and the thermophoresis constraint, while the concentration decreased. Noor et al. [36] studied the forced convection of magnetohydrodynamics (MHD) flow and its heat transfer under forced convection caused by a moving lid in a trapezoidal enclosure. The results showed that with an increase in the Reynolds number and heated length, the forced convection phenomenon became more obvious, and when the simultaneous effects of the moving lid force and heated sources moved away from the surface, the local Nusselt number decreased. Wong et al. [37] evaluated the viscous heating mechanism of Newtonian fluid filled in the cavity under the external force on the top lid by means of numerical simulation. The results showed that compared with the nonslip condition, the lower Reynolds number and larger Prandtl number (with the free-slip effect) can reduce the temperature distribution in the cavity at a faster speed, and the free-slip acted like a lubricant.
Because the multicomponent mixture can form azeotropic only in a certain proportion, none of several mixtures can form azeotropic refrigerant in a specific proportion. According to REFPROP 9.0, it is found that R744 and ethane can form azeotropic refrigerant at a ratio of 77.6/22.4. Because ethane is a flammable and explosive working medium, even if R744 can be used as a good fire-extinguishing agent to inhibit the combustion of ethane, some accidents will inevitably be encountered owing to equipment and operation problems in the experimental process, so some accidents can be avoided by using the numerical simulation method.
To sum up, this paper will conduct a numerical study on the flow boiling heat transfer characteristics of R744 and its azeotropic refrigerant. In addition, the performance of azeotropic refrigerants and R744 used in low-temperature heating systems is analyzed, so as to provide a theoretical basis for the future application of azeotropic refrigerant in system and refrigerant replacement.

2. Numerical Approach

2.1. Geometric Model

In this paper, the flow boiling heat transfer of R744 and its azeotropic refrigerant in a horizontal circular tube is studied. The geometric model is in Figure 1, because of the large amount of computation for the complex boiling flow in the three-dimensional tube grid; in order to save time and reduce the amount of computation, the geometric model is simplified as a two-dimensional rectangular fine tube with a length of 100 mm and a diameter of 1 mm.

2.2. Mathematical Model

The flow boiling problem of refrigerant in a horizontal pipe belongs to the category of a gas–liquid two-phase flow, which involves complex problems such as flow and heat transfer. The VOF multiphase flow model simulates the gas–liquid two-phase flow and heat transfer of refrigerant by solving the momentum equation, the continuity equation of each phase volume fraction, and the energy conservation equation, and the VOF model can also more accurately track the phase interface, so it is convenient to observe the flow pattern in the pipe [38]. The VOF model is often used to simulate the transient process of bubble formation during refrigerant flow boiling in channels. However, in order to obtain more-accurate flow patterns in a steady-state simulation, VOF is also applicable. Therefore, in this paper, the VOF model is chosen as the multiphase flow model for the steady-state of the flow boiling heat transfer of R744 and its azeotropic refrigerant mixture in a horizontal tube. All governing equations of the VOF model can be described as follows.

2.2.1. Basic Equation

Volume fraction equation:
q = 1 n α q = 1
where α q is the volume fraction of the q phase of the working fluid.
Continuity equation:
( ρ q v q ) = 0
where ρ q is the density of the q phase of the working fluid and v q is the velocity of the q phase of the working fluid.
Momentum equation:
( ρ v v ) = p + [ μ ( v + v T ) ] + ρ g + F v o l
F v o l = σ α q ρ q κ q α q 0.5 ( ρ 1 + ρ 2 )
κ 1 = Δ α 1 | α 1 | . κ v = Δ α v | α v |
where p is the pressure; g is the acceleration vector of gravity; F v o l is the volume force; ρ and μ are the density and the viscosity of each calculation unit, respectively; σ is the surface tension coefficient; and κ q is the interface curvature of the q phase of the working fluid. The surface tension is the surface force, which can be transformed into volume force by using the CFS (continuum surface force) model in [39].
Energy equation:
t ( ρ E ) + [ v ( ρ E + p ) ] = ( λ T ) + S E
where E is the energy; T is the temperature; λ is the thermal conductivity of the working fluid; and S E is the energy source term.

2.2.2. Turbulence Model

The RNG κ-ε model is adopted because it is suitable for simulating the complex flow and the high Re (Reynolds number), and the combination of it and the wall function method offers better performance for the near wall flow. The RNG k-ε model is given in Equation (7). The Reynolds number of R744 is 3571.4–7633.6, and the Reynolds number of the azeotropic refrigerant is 3571.4–8333.3.
( ρ k ) t + ( ρ k u i ) x i = x j [ α k μ e f f k x j ] + G k + ρ ε ( ρ ε ) t + ( ρ ε u i ) x i = x j [ α ε μ e f f k x j ] + C 2 ε * k G k C 2 ε ρ ε 2 k
where k is the turbulent kinetic energy; ε is the turbulent dissipation rate; ui is the instantaneous velocity in the rectangular coordinate system; xi, xj are the different directions in the rectangular coordinates; μ e f f is the effective dynamic viscosity coefficient; G k is the turbulent kinetic energy due to the time averaged velocity gradient; and α k = α ε = 1.39, C 2 ε * = 1.68, and C 2 ε = 1.42.

2.2.3. Interphase Transition Model

On the basis of the phase change evaporation condensation model proposed by Lee in [40], a user-defined function (UDF) program is compiled, which includes the source terms of the mass transfer and the energy transfer. The interphase mass transfer source term m can be described as follows:
m l v = r 1 α 1 ρ 1 T T s a t T s a t , T T s a t m v l = r v α v ρ v T T s a t T s a t , T < T s a t
where r 1 and r v are the interphase heat transfer coefficients, which indicate the heat exchange intensity. In the paper, both r 1 and r v are set to be 0.1. Too much coefficient will lead to the divergence of numerical calculation, and too little coefficient will cause more deviation in the gas–liquid interface temperature from the saturation temperature.
The interphase energy transfer source term Q can be described as follows:
Q = h f g m
where hfg is the vaporization latent heat.

2.3. Boundary Condition

In this paper, the velocity inlet boundary condition is selected at the entrance of the pipe, the flow direction is along the positive direction of the x-axis, and the velocity of inlet is calculated by G/ ρ . The outlet boundary condition is the outflow boundary condition, and the nonslip for the wall is set as the wall boundary condition. Keep the external boundary of the working fluid flow pipe as a fixed heat flux. The test conditions include the tube diameter (D), the mass flux (G), the heat flux (q), the saturation pressure (Psat), the saturation temperature (Tsat), and the specific simulation conditions, which are shown in Table 1.

2.4. Solution Strategy

In this paper, fluent 16.0 is used for the numerical simulation, the VOF multiphase flow model is used to track the gas–liquid interface, the energy equation is opened, the RNG K-ε turbulence model is selected, and the wall enhanced wall treatment is checked. Because this paper is based on the steady boiling process to compare the heat transfer characteristics of R744 and its azeotropic refrigerant, the common PISO algorithm has better results for the transient boiling process [40,41], while for the steady boiling process, through continuous calculation, it is found that the flow pattern and the results of the coupled algorithm are closer to the actual situation. The coupled algorithm is used for pressure and velocity coupling; the check coupled with the volume fraction gradient, in discrete options, is selected as the least-squares cell, the pressure is selected as the body force weighted, the volume fraction is selected as compressive, and the rest are selected as second-order upwind. The relaxation factor is selected as the default value.

2.5. The Thermophysical Properties of Working Fluids

Table 2 and Table 3 show the thermal properties of R744 and its azeotropic refrigerant at different saturation temperatures; all physical parameters are checked from the software of refprop9.0; and thanks to the fact that the different working fluids have different saturation pressures at the same saturation temperature, this paper only controls the inlet of the R744 and its azeotropic refrigerant, which are saturated liquids with the same saturation temperature.

2.6. Mesh Independent and Method Validation

A quadrilateral grid is used for the computational domain, and it is necessary to study the influence of the grid number on the results. Table 4 shows the local boundary layer grid of three grids, the number of the grid ranges are from 156,104 to 226,228, and the quality of the grid is 1. For the independent verification of the grid size, the working conditions (G = kg/m2·s, q = kW/m2, Tsat = K) are applied to three kinds of grids. Through the simulation, the temperature values at different positions are obtained. Taking the temperature value obtained by the grid number (a) as the benchmark (max relative error of T = 0), the maximum temperature error is compared with the other two groups of grids. Table 5 shows the relative error of grids (b) and (c) relative to (a), which are 0.053% and 0.081%, respectively. This shows that when the number of grids is larger, there is no substantial impact on the results. Finally, grid (a) is selected for the following simulation, and the grid division is shown in Figure 2.
In order to prove the reliability of the simulation, it is necessary to select the existing experimental data to compare them with the simulation values, but there is no experiment on the flow boiling heat transfer of the R744 azeotropic refrigerant at present. Therefore, the experimental conditions of pure R744 in [5] are selected to compare them with the simulated values of pure R744. Figure 3 shows the change in the flow boiling heat transfer coefficient with dryness. It can be found that the overall trend of the flow boiling heat transfer coefficient between the simulation value and the experimental value is the same, and most of the errors are within the acceptable range. However, there is a large error between the simulation data and the experimental data in the high dryness region, which is caused by the turbulence model, the multiphase flow boiling model, and the algorithm. The simulated values tend to decrease in the high dry area, while there are some special jumping points in the experimental values. In addition, the numerical calculation of the boiling heat transfer under different conditions will have certain deviations, and many research results have shown that the simulation value of the boiling heat transfer coefficient has a large error compared with the experimental value, but this paper focuses on a comparison of the heat transfer characteristics of R744 and its azeotropic refrigerant. Therefore, as long as the variation law obtained by numerical calculation is consistent with that obtained by the experiment, the simulation result is relatively reliable.

3. Result and Discussion

3.1. Pattern Distribution of Liquid Volume Fraction

Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 show the flow pattern of the flow boiling heat transfer of R744 and its azeotropic refrigerant under different working conditions. The evolution process of bubbles from the single-phase flow to the bubble flow and then to the slug flow can be clearly seen. Because of the working conditions, the overall dryness is not high, so the flow pattern reaches only the slug flow.
It can be seen from Figure 4 and Figure 5 that when the mass flux increases, the number of bubbles in the R744 tube decreases. This is because nuclear boiling dominates when the mass flux is small, and the dryness at the same position is larger, so more bubbles can be generated. With the increase in the mass flux, the number of bubbles in the azeotropic refrigerant library seems to first decrease and then increase, but in fact, the bubble flow transitions to the slug flow, and the number of bubbles does not increase.
It can be seen from Figure 6 and Figure 7 that the number of bubbles in R744 and its azeotropic refrigerant tube increases with the increase in the heat flux, which indicates that when the mass flux and the saturation temperature are constant, the boiling phenomenon becomes more intense with the increase in the heat flux.
It can be seen from Figure 8 and Figure 9 that the number of bubbles in the R744 tube first increases and then decreases with the increase in the saturation temperature. This is because when the temperature is between 253 K and 263 K, the number of bubbles increases, which means that the heat transfer effect is better. However, when the bubble decreases from 263 K to 273 K, it means that the heat transfer performance is further improved because there are enough bubbles at 263 K. If the bubbles continue to increase, the tube wall will be dried up. Therefore, the decrease in bubbles represents a further improvement in heat transfer performance. But this phenomenon is not so obvious in the azeotropic refrigerant, which may be related to its physical properties.

3.2. Effect of Mass Flux on Heat Transfer Coefficient of R744 and Azeotropic Refrigerant

Figure 10 and Figure 11 show the change in the heat transfer coefficient with the quality of R744 and its azeotropic refrigerant under different mass fluxes. It can be seen that the heat transfer coefficient of the azeotropic refrigerant is greater than that of R744 under the same working conditions. This is because from the point of view of the thermophysical properties, at the same saturation temperature, the viscosity of the azeotropic refrigerant is lower than that of R744, and the strong nucleate boiling makes the liquid film on the inner wall of the tube thinner and easier to break. Therefore, the enhancement effect of the low viscosity nucleate boiling on the heat transfer is greater than that of surface tension, and the heat transfer coefficient increases with the increase in viscosity. At the same time, low viscosity is easier to wet the wall with; therefore, the heat transfer performance of the azeotropic refrigerant is better than that of R744.
From Figure 10, we find that the heat transfer coefficient of R744 decreases with the increase in the mass flux. This is because when the mass flux is low, nucleate boiling dominates and the effect is greater than that of the mass flux. However, with the increase in quality, the heat transfer coefficient of the low mass flux sharply decreases. The reason may be that the wall surface is not timely wetted, owing to the low mass flux, while there is no drying phenomenon at the high mass flux, owing to its low dryness. However, with the increase in quality, the heat transfer coefficient of the high mass flux tends to exceed that of the low mass flux, which is also because the high mass flux can keep wetting the wall and prevent drying up. However, the azeotropic refrigerant did not dry out under the same conditions, which indicated that the occurrence of the dryness of the azeotropic refrigerant was higher than that of R744, which also showed that the heat transfer performance of the azeotropic refrigerant was better than that of R744.
In addition, it can be found from Figure 10 and Figure 11 that within the given quality range, when the mass fluxes of R744 are 50 kg/m2·s, 75 kg/m2·s, and 100 kg/m2·s, the corresponding average heat transfer coefficients are 5137.58 W/m2K, 4678.85 W/m2K, and 4656.93 W/m2K, respectively, and when the mass flux is 50 kg/m2·s, the average heat transfer coefficients are 9.8% and 10.3% higher than when the mass fluxes are 75 kg/m2·s and 100 kg/m2·s, respectively. Similarly, the corresponding average heat transfer coefficients of the azeotropic refrigerants are 6527.55 W/m2K, 5881.89 W/m2K, and 5406.93 W/m2K; when the mass flux is 50 kg/m2·s, its average heat transfer coefficient is 10.9% and 20.7% higher than when the mass flux is 75 kg/m2·s and 100 kg/m2·s, respectively. The results show that the heat transfer performance of the azeotropic refrigerant is better than R744, and the increase in its heat transfer coefficient has certain advantages.
Finally, under the given working conditions, when the mass flux of R744 is 50 kg/m2·s, 75 kg/m2·s, and 100 kg/m2·s, its heat transfer coefficient decreases by 25.8%, 7.1%, and 1.7%, respectively, with the increase in quality. The heat transfer coefficient of the azeotropic refrigerant decreases by 17.1%, 7.9%, and 2.5%, respectively.

3.3. Effect of Heat Flux on Heat Transfer Coefficient of R744 and Azeotropic Refrigerant

From Figure 12 and Figure 13, we can find the heat transfer coefficient of R744 and its azeotropic increases with the increase in the heat flux, but it decreases with the increase in quality. This is because the heat flux increases and enhances the nucleate boiling and convective heat transfer, which makes the heat transfer coefficient increase. However, because of the low mass flux, the effect of the convection heat transfer is still less than that of the nucleate boiling, which causes the wall to not be wetted in time, so the heat transfer coefficient decreases. In addition, it is found that the decrease in the heat transfer coefficient at the high heat flux is larger than that at the low heat flux, which is also because the high heat flux intensifies the nucleate boiling, increases the film on the wall, and aggravates the heat transfer deterioration.
In addition, it can be found from Figure 12 and Figure 13 that within the given quality range, when the heat flux of R744 is 10 kW/m2, 20 kW/m2, and 30 kW/m2, the corresponding average heat transfer coefficients are 5137.58 W/m2K, 9648.19 W/m2K, and 12,688.94 W/m2K, respectively, and when the heat flux is 30 kW/m2, the average heat transfer coefficients are 146.9% and 31.5% higher than when the heat fluxes are 10 kW/m2 and 20 kW/m2, respectively. Similarly, the corresponding average heat transfer coefficients of the azeotropic refrigerants are 6527.55 W/m2K, 12,176.86 W/m2K, and 16,017.91 W/m2K; when the heat flux is 30 kW/m2, its average heat transfer coefficient is 145.3% and 31.5% higher than when the heat fluxes are 10 kW/m2 and 20 kW/m2, respectively. It can be seen that under the influence of different heat fluxes, the increases in their heat transfer coefficients are basically the same.
Finally, under the given working conditions, when the heat flux of R744 is 10 kW/m2, 20 kW/m2, and 30 kW/m2, its heat transfer coefficient decreases by 25.8%, 40.8%, and 61.2%, respectively, with the increase in quality. The heat transfer coefficient of the azeotropic refrigerant decreases by 17.1%, 41.5%, and 61.7%, respectively.

3.4. Effect of Tsat on Heat Transfer Coefficient of R744 and Azeotropic Refrigerant

From Figure 14 and Figure 15, we can find that the heat transfer coefficient of R744 and its azeotropic increases with the increase in Tsat because the viscosity of the fluid decreases with the increase in the saturation temperature, and the decrease in the viscosity of the fluid increases the influence of the convective heat transfer, which is conducive to wetting the wall. At the same time, the nuclear boiling makes the liquid film on the inner wall of the tube thinner and easier to break, so the heat transfer coefficient increases when the viscosity is low. In addition, it is found that when R744 is 253 K, the heat transfer coefficient sharply decreases after the quality reaches a certain degree, which is also due to the high viscosity and the earlier drying, resulting in the deterioration of the heat transfer. However, because of the low viscosity of the azeotropic refrigerant increases, there is no sharp decrease in the heat transfer coefficient thanks to the lack of dry-out.
In addition, it can be found from Figure 14 and Figure 15 that within the given quality range, when Tsat of R744 are 253 K, 263 K, and 273 K, the corresponding average heat transfer coefficients are 5137.58 W/m2K, 6289.71 W/m2K, and 8004.55 W/m2K respectively, and when the Tsat is 273 K, the average heat transfer coefficients are 55.8% and 27.3% higher than when the Tsat values are 253 K and 263 K, respectively. Similarly, the corresponding average heat transfer coefficients of the azeotropic refrigerants are 6527.55 W/m2K, 7992.84 W/m2K, and 10,303.96 W/m2K; when the Tsat is 273 K, its average heat transfer coefficient is 61.9% and 28.9% higher than when the Tsat values are 253 K and 263 K, respectively. It can be seen that under the influence of different Tsat values, the increase in the heat transfer coefficient of the azeotropic refrigerants has certain advantages.
Finally, under the given working conditions, when the Tsat of R744 is 253 K, 263 K, and 273 K, its heat transfer coefficient decreases by 25.8%, 19.5%, and 24.6%, respectively, with the increase in quality. The heat transfer coefficient of the azeotropic refrigerant decreases by 17.1%, 20.3%, and 25.3%, respectively.

3.5. Comparison of Heat Transfer Performance of R744 and Its Azeotropic Refrigerant

Figure 16, Figure 17 and Figure 18 show the change in the heat transfer coefficient with the quality of R744 and its azeotropic refrigerant under different working conditions. It can be seen that the greater the heat transfer coefficient of the two, the greater the difference between them, indicating that the heat transfer performance of the azeotropic refrigerant is better than that of R744. On this basis, if the heat transfer is enhanced by changing the working conditions, the heat transfer performance of the azeotropic refrigerant will be even stronger in R744. In addition, it can be seen from Figure 16, Figure 17 and Figure 18 that the boiling heat transfer coefficient of the azeotropic refrigerant is 24.37% higher than that of R744 on average when the heat flux is 10 kW, the Tsat is 253–273 K, and the mass flux is 50–100 kg/m2·s. Similarly, when the mass flux is 50 kg/m2·s, the Tsat is 253 K, and the heat flux is 10–30 kW; the boiling heat transfer coefficient of the azeotropic refrigerant is 26.61% higher than that of R744 on average. When the mass flux is 50 kg/m2·s, the heat flux is 10 kW and the Tsat is 253–273 K; the boiling heat transfer coefficient of the azeotropic refrigerant is 27.73% higher than that of R744 on average.

3.6. Comparison of Low-Temperature Heating System Performance of R744 and Its Azeotropic Refrigerant

In Figure 19, the principle of the subcritical and the transcritical heating cycle systems is demonstrated. The supercritical cycle system uses the large temperature drop in the working medium to exchange heat with the outside world, so it uses a gas cooler instead of a condenser. Figure 20 shows a simple subcritical low-temperature heating cycle: the gas of refrigerant enters the compressor from the evaporator (1); the working fluid is compressed into a high-temperature gas state by the compressor, cooled by sensible heat of 2–3, and then condensed and releases heat by condensation of 3–4 (2); then through 4–5, the condensed liquid phase is led through the pressure relief expansion valve, which reduces the pressure of the liquid and partly evaporates it (3); and finally, the fluid in this state is at its saturation point when it arrives in the evaporator, where it receives heat at low temperatures until it has completely evaporated (4). After this step, the gas is again compressed in the compressor, and the whole cycle is repeated. Figure 21 shows a simple transcritical low-temperature heating cycle: the cycle process is roughly the same as that of the subcritical cycle, but the difference is that the subcritical 2–4 process changes to the 2–3 cooling process with a large temperature drop in the gas cooler, and the state point 3 is unsaturated. Table 6 shows the critical pressure and temperature of R744 and its azeotropic refrigerant. The simulation uses EES software as the calculation tool.
Figure 22 and Figure 23 show the variation in the COP and the discharge temperature with the evaporation temperature at a certain condenser outlet temperature in a subcritical system, respectively. The discharge pressure is constant to the saturated pressure of the working medium at the outlet temperature of the condenser.
As can be seen from Figure 22, the COP increases with the increase in the evaporation temperature, the COP decreases with the increase in the condensation temperature, which is due to the increase in the evaporation temperature, the increase in suction pressure, and the constant discharge pressure, which leads to the decrease in power consumption and the increase in the COP. It is worth noting that the COP of R744 is lower than that of the azeotropic refrigerant when the outlet temperature of the condenser is low, as shown in Figure 23a,b, and the increase is very small. When the outlet temperature of the condenser is high, as shown in Figure 23c, the COP of R744 is higher than that of the azeotropic refrigerant, and the increase in the COP of the azeotropic refrigerant becomes slow. It can be concluded that the outlet temperature of the condenser is the key factor in determining whether the COP of the azeotropic refrigerant exceeds that of R744.
It can be seen from Figure 23 that no matter how the condenser outlet temperature changes, the discharge temperature decreases with the increase in the evaporation temperature, and the discharge temperature of the azeotropic refrigerant is lower than R744, which indicates that the discharge temperature can be significantly reduced when using the azeotropic refrigerant compared with R744, making the system more stable.
Figure 24 and Figure 25 show the variation in the COP and the discharge temperature with the evaporation temperature at a certain condenser outlet temperature in a transcritical system, respectively. The discharge pressure is 7.5 MPa.
The change in the COP with the evaporation temperature is shown in Figure 24 when the discharge pressure is certain. Compared with Figure 22, it is obvious that there is an intersection point between R744 and azeotropic COP (corresponding to the evaporation temperature): the COP of R744 on the left side of intersection will be higher than that of the azeotropic refrigeration, while the COP of R744 on the right side will be lower than that of the azeotropic refrigeration. We call this intersection the critical point for determining whether the azeotropic refrigerant COP is higher than that of R744. In addition, the position of the critical point is determined by the condenser temperature. The position of the critical point moves to the right with the increase in the condenser temperature, which means that the COP of R744 is higher than that of the azeotropic refrigerant in a given evaporation temperature range. Therefore, if we want to ensure that the COP of the azeotropic refrigerant is higher than that of R744, we need to appropriately reduce the condenser temperature.
It can be concluded from Figure 25 that the discharge temperature of the azeotropic refrigerant is lower than that of R744 in a subcritical system or a transcritical system, which provides a strong basis for the azeotropic refrigerant to replace R744.

4. Conclusions

In this paper, in order to develop new mixed refrigerant to replace R744, a 2D numerical simulation was established to simulate the flow boiling heat transfer of R744 and its azeotropic refrigerant in a horizontal pipe and then analyze the performance of R744 and its azeotropic refrigerant in the low-temperature heating system. The following conclusions are available:
  • At the low mass flux, the heat transfer coefficients of R744 and its azeotropic refrigerant decrease with the increase in the mass flux.
  • Within a given working condition, the heat transfer coefficients of R744 and its azeotropic refrigerant increase with the increase in the heat flux and the Tsat.
  • Under the influence of different factors, the average boiling heat transfer coefficient of the azeotropic refrigerant is 24.37% (mass flux), 26.61% (heat flux), and 27.73% (Tsat) higher than R744, respectively.
  • Under the given working conditions, compared with R744, the azeotropic refrigerant is not found to be dry.
  • Compared with R744, the bubble change in the azeotropic refrigerant is more conducive to enhancing the heat transfer.
  • It is found that the critical evaporation temperature makes the COP of the azeotropic refrigerant higher than that of R744.
  • The critical evaporation temperature is determined by Tco, and it increases with the increase in the Tco.
In conclusion, compared with R744, the azeotropic refrigerant has certain advantages in the boiling heat transfer, and by controlling the condensation temperature, the system performance of the azeotropic refrigerant can also be better than that of R744. According to the conclusion of this paper, it has certain reference significance for the replacement research of R744. At the same time, relevant experimental research will be carried out in the future to better understand its heat transfer and system characteristics.

Author Contributions

Conceptualization, D.S.; methodology, D.S.; software, X.Z.; validation, Z.L.; formal analysis, D.S.; investigation, H.Z.; resources, Z.L.; data curation, D.S.; writing—original draft preparation, D.S.; writing—review and editing, D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is unavailable due to privacy.

Conflicts of Interest

We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work, there is no professional or other personal interest of any nature or kind in any product, service and/or company that could be construed as influencing the position presented in, or the review of, the manuscript entitled.

Nomenclature

c p specific heat at constant pressure (kJ/kg·K)
C 2 ε turbulence model’s constants
C 2 ε * turbulence model’s constants
Dtube diameter (m)
Eenergy
Fvolume force
gacceleration due to gravity (m/s2)
Gmass flux (kg/(m2·s))
Gkturbulent kinetic energy due to time averaged velocity gradient
hheat transfer coefficient (W/m2·K)
hgfvaporization latent heat (kJ/kg)
kturbulence kinetic energy (m2/s2)
mmass of working fluid (kg)
ppressure (MPa)
Qquantity of heat (kW)
qheat flux (kW/m2)
ReReynolds number
r interphase heat transfer coefficients
ssource term
Tthermodynamic temperature (k)
vvelocity (m/s)
xcartesian coordinates (m)
Xquality
Greek symbols
α volume fraction
ε rate of dissipation of k (m2/s3)
λ thermal conductivity (W/(m·K))
μ dynamic viscosity (Pa·s)
ρ density (kg/m3)
σ tension (N/m)
σ k turbulent Prandtl number for k
σ ε turbulent Prandtl number for ε
κ interface curvature
Subscripts
cooutlet temperature of condenser
ddischarge temperature
eevaporation
effeffective
lliquid phase
i,jgeneral spatial indices
qphase
satsaturation
vvapor phase
Acronyms
R744carbon dioxide
COPcoefficient of performance
GWPglobal warming potential
ODPozone-depleting potential

References

  1. Grauso, S.; Mastrullo, R.; Mauro, A.W.; Vanoli, G.P. R744 and propane blends: Experiments and assessment of predictive methods for flow boiling in horizontal tubes. Int. J. Refrig. 2011, 34, 1028–1039. [Google Scholar] [CrossRef]
  2. Afroz, H.M.; Miyara, A.; Tsubaki, K. Heat transfer coefficients and pressure drops during in-tube condensation of R744/DME mixture refrigerant. Int. J. Refrig. 2008, 31, 1458–1466. [Google Scholar] [CrossRef]
  3. Onaka, Y.; Miyara, A.; Tsubaki, K. Experimental study on evaporation heat transfer of R744/DME mixture refrigerant in a horizontal smooth tube. Int. J. Refrig. 2010, 33, 1277–1291. [Google Scholar] [CrossRef]
  4. Zhu, Y.; Wu, X.; Wei, Z. Heat transfer characteristics and correlation for R744/propane mixtures flow evaporation in a smooth mini tube. Appl. Therm. Eng. 2015, 81, 253–261. [Google Scholar] [CrossRef]
  5. Yun, R.; Kim, Y.; Kim, M.S.; Choi, Y. Boiling heat transfer and dryout phenomenon of R744 in a horizontal smooth tube. Int. J. Heat Mass Transf. 2003, 45, 2353–2361. [Google Scholar] [CrossRef]
  6. Zhao, X.; Bansal, P. Flow boiling heat transfer characteristics of CO2 at low temperatures. Int. J. Refrig. 2007, 30, 937–945. [Google Scholar] [CrossRef]
  7. Oh, H.-K.; Son, C.-H. Flow boiling heat transfer and pressure drop characteristics of CO2 in horizontal tube of 4.57-mm inner diameter. Appl. Therm. Eng. 2011, 31, 163–172. [Google Scholar] [CrossRef]
  8. Yun, R.; Kim, Y.; Kim, M.S. Convective boiling heat transfer characteristics of CO2 in microchannels. Int. J. Heat Mass Transf. 2005, 48, 235–242. [Google Scholar] [CrossRef]
  9. Ozawa, M.; Ami, T.; Ishihara, I.; Umekawa, H.; Matsumoto, R.; Tanaka, Y.; Yamamoto, T.; Ueda, Y. Flow pattern and boiling heat transfer of CO2 in horizontal small-bore tubes. Int. J. Multiph. Flow 2009, 35, 699–709. [Google Scholar] [CrossRef]
  10. Huai, X.; Koyama, S.; Zhao, T.; Shinmura, E.; Hidehiko, K.; Masaki, M. An experimental study of flow boiling characteristics of carbon dioxide in multiport mini channels. Appl. Therm. Eng. 2004, 24, 1143–1463. [Google Scholar] [CrossRef]
  11. Jiang, L.; Liu, J.; Zhang, L.; Liu, Q.; Xu, X. Characteristics of heat transfer for CO2 flow boiling at low temperature in mini-channel. Int. J. Heat Mass Transf. 2017, 108, 2120–2129. [Google Scholar] [CrossRef]
  12. Liang, Z.; Linlin, J.; Jianhua, L.; Yue, Z. Investigation of flow boiling heat transfer characteristics of CO2 in horizontal mini-tube. Int. J. Therm. Sci. 2019, 138, 109–115. [Google Scholar] [CrossRef]
  13. 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]
  14. Wei, Z.F.; Wu, X.M.; MO, S.J. Numerical simulation on flow boiling heat transfer of CO2/propane mixture in tubes with/without micro-fins. J. Eng. Thermophys. 2011, 32, 1019–1022. [Google Scholar]
  15. Wu, X.M.; Zhao, R.; Wei, Z.F.; Huang, X.J.; Wang, W.C. Experimental studies on flow heat transfer of R744/R290 mixture in a horizontal tube. J. Eng. Thermophys. 2013, 34, 706–709. [Google Scholar]
  16. Kim, S.G.; Kim, M.S. Experiment and simulation on the performance of an auto-cascade refrigeration system using carbon dioxide as a refrigerant. Int. J. Refrig. 2002, 25, 1093–1101. [Google Scholar] [CrossRef]
  17. Sarkar, J.; Bhattacharyya, S. Assessment of blends of R744 with butane and isobutane as working fluids for heating applications. Int. J. Therm. Sci. 2009, 48, 1460–1465. [Google Scholar] [CrossRef]
  18. Niu, B.; Zhang, Y. Experimental study of the refrigeration cycle performance for the R744/R290 mixtures. Int. J. Refrig. 2007, 30, 37–42. [Google Scholar] [CrossRef]
  19. Fukuta, M. Performance and characteristics of compressor/expander combination for R744 cycle. In Proceedings of the 7th IIR Gustav Lorentzen Conference on Natural Working Fluids, Trondheim, Norway, 29–31 May 2006. [Google Scholar]
  20. Nickl, J.; Will, G.; Quack, H.; Kraus, W.E. Integration of a three-stage expander into a R744 refrigeration system. In Proceedings of the 6th IIR Gustav-Lorentzen Natural Working Fluids Conf, Glasgow, UK, 29 August 2004. 4/A/12.20. [Google Scholar]
  21. Bo, Z.; Xueyuan, P.; Bei, G.; Ziwen, X.; Pengcheng, S. Design and experimental validation of the slider-based free piston expander for trancritical R744 refrigeration cycle. In Proceedings of the 7th IIR Gustav Lorentzen Conference on Natural Working Fluids, Trondheim, Norway, 29–31 May 2006. [Google Scholar]
  22. Huff, H.J.; Radermacher, R.; Preissner, M. Experimental investigation of a scroll expander in a carbon dioxide air-conditioning system. In Proceedings of the 21th International Congress of Refrigeration Washington, Washington, DC, USA, 17–22 August 2003. ICR485. [Google Scholar]
  23. Hua, T.; Yitai, M.; Minxia, L.; Wei, W. Study on Expansion Power Recovery in R744 Trans-Critical Cycle. Energy Convers. Manag. 2010, 51, 2516–2522. [Google Scholar]
  24. Hua, T.; Yitai, M.; Minxia, L.; Haiqing, G.; Zhongyan, L. Experimental research on impacts of non-condensation gas to expander performance used in carbon dioxide heating. Appl. Therm. Eng. 2011, 31, 1943–1949. [Google Scholar] [CrossRef]
  25. Pettersen, J.; Rieberer, R.; Leister, A. Heat transfer and pressure drop characteristics of supercritical carbon dioxide in micro-channel tubes under cooling. In Preliminary Proceedings of the 4th IIR-Gustav Lorentzen Conference on Natural Working Fluids at Purdue, West Lafayette, IN, USA, 25 July 2000; pp. 99–106. [Google Scholar]
  26. Olson, D.A. Heat transfer of supercritical carbon dioxide flowing in a cooled horizontal tube. In Preliminary Proceedings of the 4th IIR-Gustav Lorentzen Conference on Natural Working Fluids at Purdue, West Lafayette, IN, USA, 25 July 2000; pp. 251–258. [Google Scholar]
  27. Liao, S.M.; Zhao, T.S. Measurements of heat transfer coefficients from supercritical carbon dioxide flowing in horizontal mini/micro channels. Trans. ASME J. Heat Transf. 2002, 124, 413–420. [Google Scholar] [CrossRef]
  28. Pitla, S.S.; Groll, E.A.; Ramadhhyani, S. New correlation to predict the heat transfer coefficient during in-tube cooling of turbulent supercritical R744. Int. J. Refrig. 2002, 25, 887–895. [Google Scholar] [CrossRef]
  29. Pitla, S.S.; Bhatia, K.; Khetarpal, V.; Strikis, G. Numerical heat transfer analysis in heat exchangers for transcritical R744 system. In Preliminary Proceedings of the 4th IIR-Gustav Lorentzen Conference on Natural Working Fluids at Purdue, West Lafayette, IN, USA, 25 July 2000; pp. 307–314. [Google Scholar]
  30. Liao, S.M.; Zhao, T.S. A numerical investigation of laminar convection of supercritical carbon dioxide in vertical mini/micro tubes. Prog. Comput. Fluid Dyn. 2002, 2, 144–152. [Google Scholar] [CrossRef]
  31. Zingerli, A.; Groll, E.A. Influence of refrigeration oil on the heat transfer and pressure drop of supercritical R744 during in-tube cooling. In Preliminary Proceedings of the 4th IIR-Gustav Lorentzen Conference on Natural Working Fluids at Purdue, West Lafayette, IN, USA, 25 July 2000; pp. 269–278. [Google Scholar]
  32. Safia, A.; Maria, A.; Khalid, S.; Alia, R.; Metib, A.; Taseer, A. Impact of Partial Slip on Double Diffusion Convection of Sisko Nanofluids in Asymmetric Channel with Peristaltic Propulsion and Inclined Magnetic Field. Nanomaterials 2022, 12, 2736. [Google Scholar]
  33. Yasir, K.; Safia, A.; Alia, R.; Aewar, H.; Alsulaimani, H.A. Effects of Double Diffusive Convection and Inclined Magnetic Field on the Peristaltic Flow of Fourth Grade Nanofluids in a Non-Uniform Channel. Nanomaterials 2022, 12, 3037. [Google Scholar]
  34. Athar, M.; Khan, Y.; Akram, S.; Saeed, K.; Alameer, A.; Hussain, A. Consequence of Double-Diffusion Convection and Partial Slip on Magneto-Oldroyd-4 Constants Nanofluids with Peristaltic Propulsion in an Asymmetric Channel. Complexity 2022, 2022, 7634357. [Google Scholar] [CrossRef]
  35. Noor, N.F.M.; Haq, R.U.; Wong, H.F.; Alzahrani, A.K.; Ullah, M.Z. Flow and heat transfer due to partially heated moving lid in a trapezoidal cavity with different constraints at inner circular obstacle. Int. Commun. Heat Mass Transf. 2022, 135, 106111. [Google Scholar] [CrossRef]
  36. Wong, H.F.; Sohail, M.; Siri, Z.; Noor, N.F.M. Numerical Solutions for Heat Transfer of An Unsteady Cavity with Viscous Heating. Comput. Mater. Contin. 2021, 68, 319–336. [Google Scholar] [CrossRef]
  37. Xie, L.; Xie, Y.; Yu, J. Phase distributions of boiling flow in helical coils in high gravity. Int. J. Heat Mass Transf. 2015, 80, 7–15. [Google Scholar] [CrossRef]
  38. Brackbill, J.U.; Kothe, D.B.; Zemach, C. A continuum method of modeling surface tension. J. Comput. Phusis 1992, 10, 335–354. [Google Scholar] [CrossRef]
  39. Yang, Z.; Peng, X.F.; Ye, P. Numerical and experimental investigation of two phase flow during boiling in a coiled tube. Int. J. Heat Mass Transf. 2008, 51, 1003–1016. [Google Scholar] [CrossRef]
  40. Yang, Z.; Peng, X.F.; Ye, P. Numerical study of flow boiling flow patterns and pressure drop of R134a in small tubes under high flight acceleration. J. Braz. Soc. Mech. Sci. Eng. 2020, 42, 482. [Google Scholar]
  41. Niu, X.; Luo, S.; Fan, L.L.; Zhao, L. Numerical simulation on the flow and heat transfer characteristics in the one-side heating helically coiled tube. Appl. Therm. Eng. 2016, 106, 579–587. [Google Scholar] [CrossRef]
Figure 1. Geometric model.
Figure 1. Geometric model.
Energies 16 01313 g001
Figure 2. Grid diagram of computing domain.
Figure 2. Grid diagram of computing domain.
Energies 16 01313 g002
Figure 3. Change in heat transfer coefficient with dryness (the experimental conditions include D = 6 mm, G = 240 kg/m2·s, Tsat = 283 K, and q = 10–20 kW/m2 [7]).
Figure 3. Change in heat transfer coefficient with dryness (the experimental conditions include D = 6 mm, G = 240 kg/m2·s, Tsat = 283 K, and q = 10–20 kW/m2 [7]).
Energies 16 01313 g003
Figure 4. Volume fraction distribution of R744 liquid phase at different mass fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.28; (b) Tsat = 253 K, G = 75 kg/m2·s, q = 10 kW, X = 0.19; (c) Tsat = 253 K, G = 100 kg/m2·s, q = 10 kW, X = 0.14.
Figure 4. Volume fraction distribution of R744 liquid phase at different mass fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.28; (b) Tsat = 253 K, G = 75 kg/m2·s, q = 10 kW, X = 0.19; (c) Tsat = 253 K, G = 100 kg/m2·s, q = 10 kW, X = 0.14.
Energies 16 01313 g004
Figure 5. Volume fraction distribution of R744 liquid phase at different heat fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.28; (b) Tsat = 253 K, G = 50 kg/m2·s, q = 20 kW, X = 0.56; (c) Tsat = 253 K, G = 50 kg/m2·s, q = 30 kW, X = 0.84.
Figure 5. Volume fraction distribution of R744 liquid phase at different heat fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.28; (b) Tsat = 253 K, G = 50 kg/m2·s, q = 20 kW, X = 0.56; (c) Tsat = 253 K, G = 50 kg/m2·s, q = 30 kW, X = 0.84.
Energies 16 01313 g005
Figure 6. Volume fraction distribution of R744 liquid phase at different Tsat values: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.28; (b) Tsat = 263 K, G = 50 kg/m2·s, q = 10 kW, X = 0.31; (c) Tsat = 273 K, G = 50 kg/m2·s, q = 10 kW, X = 0.35.
Figure 6. Volume fraction distribution of R744 liquid phase at different Tsat values: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.28; (b) Tsat = 263 K, G = 50 kg/m2·s, q = 10 kW, X = 0.31; (c) Tsat = 273 K, G = 50 kg/m2·s, q = 10 kW, X = 0.35.
Energies 16 01313 g006
Figure 7. Volume fraction distribution of az-refrigerant liquid phase at different mass fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.26; (b) Tsat = 253 K, G = 75 kg/m2·s, q = 10 kW, X = 0.18; (c) Tsat = 253 K, G = 100 kg/m2·s, q = 10 kW, X = 0.13.
Figure 7. Volume fraction distribution of az-refrigerant liquid phase at different mass fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.26; (b) Tsat = 253 K, G = 75 kg/m2·s, q = 10 kW, X = 0.18; (c) Tsat = 253 K, G = 100 kg/m2·s, q = 10 kW, X = 0.13.
Energies 16 01313 g007
Figure 8. Volume fraction distribution of az-refrigerant liquid phase at different heat fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.26; (b) Tsat = 253 K, G = 50 kg/m2·s, q = 20 kW, X = 0.52; (c) Tsat = 253 K, G = 50 kg/m2·s, q = 30 kW, X = 0.78.
Figure 8. Volume fraction distribution of az-refrigerant liquid phase at different heat fluxes: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.26; (b) Tsat = 253 K, G = 50 kg/m2·s, q = 20 kW, X = 0.52; (c) Tsat = 253 K, G = 50 kg/m2·s, q = 30 kW, X = 0.78.
Energies 16 01313 g008
Figure 9. Volume fraction distribution of az-refrigerant liquid phase at different Tsat values: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.26; (b) Tsat = 263 K, G = 50 kg/m2·s, q = 10 kW, X = 0.29; (c) Tsat = 273 K, G = 50 kg/m2·s, q = 10 kW, X = 0.32.
Figure 9. Volume fraction distribution of az-refrigerant liquid phase at different Tsat values: (a) Tsat = 253 K, G = 50 kg/m2·s, q = 10 kW, X = 0.26; (b) Tsat = 263 K, G = 50 kg/m2·s, q = 10 kW, X = 0.29; (c) Tsat = 273 K, G = 50 kg/m2·s, q = 10 kW, X = 0.32.
Energies 16 01313 g009
Figure 10. Effect of mass flux on heat transfer coefficient of R744.
Figure 10. Effect of mass flux on heat transfer coefficient of R744.
Energies 16 01313 g010
Figure 11. Effect of mass flux on heat transfer coefficient of azeotropic refrigerant.
Figure 11. Effect of mass flux on heat transfer coefficient of azeotropic refrigerant.
Energies 16 01313 g011
Figure 12. Effect of heat flux on heat transfer coefficient of R744.
Figure 12. Effect of heat flux on heat transfer coefficient of R744.
Energies 16 01313 g012
Figure 13. Effect of heat flux on heat transfer coefficient of azeotropic refrigerant.
Figure 13. Effect of heat flux on heat transfer coefficient of azeotropic refrigerant.
Energies 16 01313 g013
Figure 14. Effect of Tsat on heat transfer coefficient of R744.
Figure 14. Effect of Tsat on heat transfer coefficient of R744.
Energies 16 01313 g014
Figure 15. Effect of Tsat on heat transfer coefficient of azeotropic refrigerant.
Figure 15. Effect of Tsat on heat transfer coefficient of azeotropic refrigerant.
Energies 16 01313 g015
Figure 16. Comparison of heat transfer coefficients between azeotropic refrigerant and R744 at different mass fluxes: (a) G = 100 kg/m2·s; (b) G = 75 kg/m2·s; (c) G = 50 kg/m2·s.
Figure 16. Comparison of heat transfer coefficients between azeotropic refrigerant and R744 at different mass fluxes: (a) G = 100 kg/m2·s; (b) G = 75 kg/m2·s; (c) G = 50 kg/m2·s.
Energies 16 01313 g016
Figure 17. Comparison of heat transfer coefficients between azeotropic refrigerant and R744 at different heat fluxes: (a) q = 10 kW; (b) q = 20 kW; (c) q = 30 kW.
Figure 17. Comparison of heat transfer coefficients between azeotropic refrigerant and R744 at different heat fluxes: (a) q = 10 kW; (b) q = 20 kW; (c) q = 30 kW.
Energies 16 01313 g017
Figure 18. Comparison of heat transfer coefficients between azeotropic refrigerant and R744 at different Tsat values: (a) Tsat = 253 K; (b) Tsat = 263 K; (c) Tsat = 273 K.
Figure 18. Comparison of heat transfer coefficients between azeotropic refrigerant and R744 at different Tsat values: (a) Tsat = 253 K; (b) Tsat = 263 K; (c) Tsat = 273 K.
Energies 16 01313 g018
Figure 19. Scheme of the subcritical vapor compression heating working cycle.
Figure 19. Scheme of the subcritical vapor compression heating working cycle.
Energies 16 01313 g019
Figure 20. T-S diagram of the subcritical vapor compression heating working cycle.
Figure 20. T-S diagram of the subcritical vapor compression heating working cycle.
Energies 16 01313 g020
Figure 21. T-S diagram of the transcritical vapor compression heating working cycle.
Figure 21. T-S diagram of the transcritical vapor compression heating working cycle.
Energies 16 01313 g021
Figure 22. Comparison of COP values of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser: (a) Tco = 273 K; (b) Tco = 278 K; (c) Tco = 283 K.
Figure 22. Comparison of COP values of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser: (a) Tco = 273 K; (b) Tco = 278 K; (c) Tco = 283 K.
Energies 16 01313 g022
Figure 23. Comparison of discharge temperature of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser: (a) Tco = 273 K; (b) Tco = 278 K; (c) Tco = 283 K.
Figure 23. Comparison of discharge temperature of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser: (a) Tco = 273 K; (b) Tco = 278 K; (c) Tco = 283 K.
Energies 16 01313 g023
Figure 24. Comparison of COP of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser: (a) Tco = 273 K; (b) Tco = 278 K; (c) Tco = 283 K.
Figure 24. Comparison of COP of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser: (a) Tco = 273 K; (b) Tco = 278 K; (c) Tco = 283 K.
Energies 16 01313 g024
Figure 25. Comparison of discharge temperature of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser.
Figure 25. Comparison of discharge temperature of azeotropic refrigerant and R744 varies with evaporation temperature at different outlet temperatures of condenser.
Energies 16 01313 g025
Table 1. Simulation conditions.
Table 1. Simulation conditions.
CaseD (mm)G (kg/m2·s)Q (kW/m2)Psat (MPa)Tsat (K)
R744150/75/10010/20/302.3/3.03/3.92253/263/273
azeotropic refrigerant150/75/10010/20/301.96/2.64/3.47253/263/273
Table 2. Thermophysical parameters of R744.
Table 2. Thermophysical parameters of R744.
ParameterR744R744R744
Tsat (K)253263273
ρ 1 (kg/m3)1032.4983.7928.33
ρ v (kg/m3)51.4570.8597.18
CPL (kJ/(kg·K))2.16362.30462.5377
CPV (kJ/(kg·K))1.28671.50501.8578
λ 1 (mW/(m·K))134.82122.72110.61
λ v (mW/(m·K))15.06916.92819.621
μ 1 (μPa·s)1.4 × 10−41.2 × 10−41.0 × 10−4
μ v (μPa·s)1.31 × 10−51.38 × 10−51.4 × 10−5
hgf (kJ/kg)282.78259231.35
σ (N/m)8.626.534.57
Psat (MPa)1.962.643.47
Table 3. Thermophysical parameters of the azeotropic refrigerant.
Table 3. Thermophysical parameters of the azeotropic refrigerant.
ParameterAz-RefrigerantAz-RefrigerantvAz-Refrigerant
Tsat (K)253263273
Liquid density (kg/m3)748.9704.14649.59
Vapor density (kg/m3)58.9981.49114.01
L of specific heat capacity (kJ/(kg·K))2.52782.82313.4261
V of specific heat capacity (kJ/(kg·K))1.76032.19783.1137
T-conductivity of l (mW/(m·K))109.5398.5687.47
T-conductivity of v (mW/(m·K))17.5420.5325.54
Liquid viscosity (μPa·s)0.9 × 10−40.78 × 10−40.65 × 10−4
Vapor viscosity (μPa·s)1.2 × 10−51.3 × 10−51.4 × 10−5
Latent heat of vaporization (kJ/kg)301.12276.34247.42
Tension (N/m)5.243.562.03
Psat (MPa)2.33.033.92
Table 4. The local boundary layer grid of three grids.
Table 4. The local boundary layer grid of three grids.
Number of Grid156,104182,104226,228
Grid Energies 16 01313 i001 Energies 16 01313 i002 Energies 16 01313 i003
(a)(b)(c)
Grid quality111
Table 5. Independent verification of grid size.
Table 5. Independent verification of grid size.
Number of GridsGrid QualityMax Relative Error of T
156,10410%
182,10410.053%
226,22810.081%
Table 6. Critical parameters of R744 and azeotropic refrigerant.
Table 6. Critical parameters of R744 and azeotropic refrigerant.
ParameterCritical Pressure (MPa)Critical Temperature (K)
R7447.37304.13
azeotropic refrigerant4.87305.32
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, D.; Zhang, X.; Liu, Z.; Zhang, H. Comparative Study on Boiling Heat Transfer Characteristics and Performance of Low-Temperature Heating System of R744 and Its Azeotropic Refrigerant. Energies 2023, 16, 1313. https://doi.org/10.3390/en16031313

AMA Style

Sun D, Zhang X, Liu Z, Zhang H. Comparative Study on Boiling Heat Transfer Characteristics and Performance of Low-Temperature Heating System of R744 and Its Azeotropic Refrigerant. Energies. 2023; 16(3):1313. https://doi.org/10.3390/en16031313

Chicago/Turabian Style

Sun, Dahan, Xin Zhang, Zhongyan Liu, and Hao Zhang. 2023. "Comparative Study on Boiling Heat Transfer Characteristics and Performance of Low-Temperature Heating System of R744 and Its Azeotropic Refrigerant" Energies 16, no. 3: 1313. https://doi.org/10.3390/en16031313

APA Style

Sun, D., Zhang, X., Liu, Z., & Zhang, H. (2023). Comparative Study on Boiling Heat Transfer Characteristics and Performance of Low-Temperature Heating System of R744 and Its Azeotropic Refrigerant. Energies, 16(3), 1313. https://doi.org/10.3390/en16031313

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

Article Metrics

Back to TopTop