1. Introduction
The use of air conditioning systems is expanding rapidly around the world. An estimated amount of 700 million air conditioners will be operating in the world by 2030 [
1]. This growing demand for air conditioning systems has enormous impacts on the environment.
Currently, a number of present regulations have been applied worldwide to control the use of harmful refrigerants [
2,
3,
4]. The key implications of the use of conventional refrigerants include the depletion of the ozone layer and global warming. Based on the Montreal Protocol, a complete abolishment of chlorofluorocarbons (CFCs) was decided, due to their high ozone depletion potential (ODP) [
5]. In addition, the phase out of hydrochlorofluorocarbon (HCFC) refrigerants was implemented [
6]. On the other hand, the F-gas Regulation, first issued in the European Union in 2006, aimed to introduce measures for the reduction of fluorinated gases—hydrofluorocarbons (HFCs) and perfluorocarbons (PFCs)—in form of a phase-down, due to their high global warming potential (GWP) [
7].
Instead, natural refrigerants are proposed as the substitute for the harmful refrigerants. Carbon dioxide (CO2) is a natural, low cost, non-flammable, non-toxic refrigerant. Subsequently, it has emerged as a credible natural refrigerant to replace HFCs and HCFCs. However, its unique critical point, high critical pressure of 73.8 bar, and low critical temperature of 30.98 °C, remarkably affects the performance of CO2 refrigeration systems, as well as imposes special design and control challenges. When the ambient temperature is higher than the critical temperature of CO2, the system operates at supercritical conditions, and the heat rejection process occurs at a supercritical regime. In consequence, a phase change does not take place, and the heat exchanger in which this change of state occurs is called the gas cooler.
The impact of the gas cooler on the CO
2 refrigeration systems plays an important role, due to its high exergy loss. Therefore, it is considered vital to be further investigated and designed properly [
8]. The finned-tube type for gas coolers is well established in the heating, ventilation, and air conditioning (HVAC) and refrigeration industries, due to its compactness and manufacturing flexibility. The design of the finned-tube heat exchangers affects considerably the overall heat transfer performance and system efficiency. Particularly, fin and tube thickness and the respective materials, spacing, and dimensions of the tubes and fins are crucial parameters of the design [
9]. Fundamental studies about the heat transfer characteristics during the heat rejection process in tubes have been performed theoretically and experimentally by many researchers since Lorentzen and Pettersen [
10] proposed the transcritical CO
2 cycle for mobile air conditioning systems.
Pitla et al. [
11] conducted an investigation about heat transfer phenomena and pressure losses of CO
2 at supercritical conditions into a tube. They found that the majority of the deviations between the numerical and experimental values are within ± 20%, and a new heat transfer correlation was presented. Son and Park [
12] carried out an experiment in order to investigate the gas cooling process of CO
2 in terms of heat transfer coefficient and pressure drop of the refrigerant. They described the variations of local heat transfer coefficient in the cooling process in the direction of the flow and proposed a more accurate heat transfer correlation. Zhang et al. [
13] evaluated the performance of a printed circuit heat exchanger for cooling CO
2 with water. The analysis concluded that rapid variations in the thermodynamic properties of supercritical CO
2 increase entropy generation and therefore, to optimize the second law efficiency of the investigated heat exchanger, higher CO
2 mass flow rates should be used. Jadhav et al. [
14] evaluated, using simulations, CO
2 gas coolers for air conditioning applications. For their investigation, a counter crossflow plain fin and staggered tube configurations were considered. According to the study, transverse tube spacing, gas cooler width, and air volumetric flow were the most influential parameters in the heat transfer mechanisms of the gas cooler.
Liu et al. [
15] investigated experimentally the supercritical CO
2 characteristics in horizontal tubes with inner diameter of 4, 6, and 10.7 mm in terms of heat transfer phenomena and pressure losses. The authors concluded that the tube diameter significantly affects the heat transfer performance, and they proposed a new heat transfer correlation for the large diameter. Jiang et al. [
16] investigated the convection heat transfer of CO
2 at supercritical pressures in a vertical small tube with inner diameter of 2.0 mm, experimentally and numerically. They studied the effects of various operational parameters and buoyancy on convection heat transfer in a small diameter. They concluded that when the CO
2 bulk temperatures are in the near-critical region, the local heat transfer coefficients vary significantly along the tube. Chai et al. [
17] investigated, using simulations, the performance of finned-tube supercritical CO
2 gas coolers, combining a distributed modeling approach with the ε-NTU method. The results indicated that the performance of the gas coolers was enhanced by higher mass flow rates and lower tube diameters at the expense of higher pressure drops.
Other researchers have also investigated the performance of the air-cooled CO
2 gas coolers. Cheng et al. [
18] presented an analysis of heat transfer and pressure drop experimental data and correlations for supercritical CO
2 cooling in macro- and micro-channels. Ge and Cropper [
19] presented a detailed mathematical model for air-cooled finned-tube CO
2 gas coolers. They used a distributed method in order to obtain more accurate refrigerant thermophysical properties and local heat transfer coefficients during cooling processes. The model was compared with published test results. The comparison showed that the approach temperature and the heat capacity are simultaneously improved with the increase of heat exchanger circuit numbers. Marcinichen et al. [
20] conducted simulations to optimize the working fluid charge of the gas cooler. The optimal design of the study reduced CO
2 charge by 14%, compared to a reference design. Moreover, the analysis revealed the importance of the oil concentrations in the CO
2 pressure drop, which is up to 2.65 times higher for oil concentrations of up to 3%. Zilio et al. [
21] experimentally evaluated two different gas coolers, one with continuous, and one with separated fins, and on two different circuit arrangements for a transcritical CO
2 cycle. Using a coil with fins, a heat flux improvement of up to 5.6% was identified, which corresponded to a coefficient of performance (COP) increase of up to 6.6% for a conventional CO
2 refrigeration cycle.
Gupta and Dasgupta [
22] applied a similar modelling method to the one from Ge and Cropper, [
19] in order to evaluate the performance of the heat exchanger being affected by the airflow velocity. Here, a higher gas cooler performance is achieved at a higher air flow velocity as it decreases the refrigerant’s approach temperature, and thus the heating capacity of the gas cooler is increased. Santosa et al. [
23] built two CO
2 finned-tube gas coolers with different structural designs and controls, connected with a test rig of a CO
2 booster refrigeration system. They carried out experiments at different operating conditions while they developed models of the finned-tube CO
2 gas cooler. The analysis was conducted based on the distributed and lumped methods. They concluded that the heat exchanger design can affect the performance of both the component and the integrated system.
Although the heat transfer and pressure drop characteristics of the supercritical CO2 in tubes have been investigated extensively using experimental and theoretical methods, research on the air-cooled finned-tube CO2 gas coolers is still limited.
In this paper, mathematical calculations of the finned-tube CO2 gas cooler are conducted, in order to establish a reliable design procedure. With focus on the heat transfer characteristics of the air side, the developed model is validated with an experimental setup using water as working fluid. Investigations of the effects of fan frequency, water inlet temperature, and water mass flow on the overall heat transfer coefficient are conducted, while deviations between the model and the test results are extracted according to the fan frequency. In addition, potential heat transfer correlations for the air- and refrigerant-side heat transfer coefficients have been studied. Finally, the model was applied to identify a reliable and efficient finned-tube CO2 gas cooler design, as well as to evaluate its performance in different off-design conditions under varying ambient temperatures.
2. Materials and Methods
This study is part of a larger project of CO
2 air conditioning systems for residential applications and focuses on the gas cooler. A scheme of the considered CO
2 air conditioning system is depicted in
Figure 1. Particularly, an efficient and reliable air finned-tube gas cooler was designed based on the boundary conditions, which are given in
Table 1.
In order to design the CO
2 gas cooler, a script in MATLAB R2019a [
24] was developed. The simulation model calculated the overall heat transfer coefficient
U of the gas cooler, based on the mass flows, inlet and outlet temperatures, and pressures of the medium. Subsequently, the required exchange area
AR was determined. Both parameters were based on the following equations:
where
and
Ai represent the outer and inner surface of the tube, respectively.
The U-value consists of three parts: air convection, refrigerant convection, and the conduction. Compared to the other two heat transfer contributions, the conduction plays a minor role. The heat transfer coefficients on the air- and refrigerant-side are crucial for the overall heat transfer coefficient, and the validation for them are considered necessary. In order to validate the calculations, the overall heat transfer coefficient of a defined air-cooled heat exchanger was investigated experimentally.
The experimental part was based on a test rig using water and employed a specific design of an air-cooled heat exchanger (HEX). The HEX type was a finned tube with a fan air cooling system. The U-value was investigated experimentally for the entire HEX for different conditions. The second part of the validation consisted of modelling the heat exchanger to simulate the finned tube of the HEX.
2.1. Theoretical Model
For the model calculations, a script was created in MATLAB R2019a [
24], modelling the defined gas cooler. The model overall heat transfer coefficient of the HEX was calculated from Equation (1).
2.1.1. Air-Side Heat Transfer
In-line arrangement and circular finned tubes were assumed, in order to model the HEX. Based on the assumption of crossflow type, the air- and refrigerant-side heat transfer coefficients can be calculated. The proposed correlation from VDI-Heat Atlas [
25] was used in order to calculate the Nusselt number, using the following equation:
with
C = 0.22 for in-line arrangement.
/
is the ratio of the finned surface to the surface of the base tube, and for circular fins was calculated from the following equation:
The Reynolds number was calculated by the equation:
where
is the velocity in the smallest cross-section and was calculated from the following equation:
The air-side heat transfer coefficient was calculated from its definition, as the following equation shows.
However, the air-side heat transfer coefficient was affected by the fins. The fins should be taken into consideration, thus the following equation was used [
25]:
The fin efficiency is defined as the ratio of the heat removed by the fin to the heat removed by the fin at wall temperature. The efficiency of the fin was calculated from the following equation:
with [
25]:
and
for circular fins [
25]. So, the equation of the overall heat transfer coefficient used the updated air-side heat transfer coefficient as follows:
2.1.2. Refrigerant-Side Heat Transfer
On the refrigerant-side, the Gnielinski correlation [
26] was used to calculate the Nusselt number:
which is valid in the range of
.
In addition, the friction factor
was calculated by:
Here, the Reynolds number was defined as
and
G was defined as the mass velocity, and was calculated from the following equation:
Finally, heat transfer coefficient at the refrigerant side was calculated from:
Further investigation of potential heat transfer correlations was conducted. More specifically, comparisons between the correlations of Gnielinski [
26] and Dittus–Boelter [
27], and the correlations proposed by VDI-Heat Atlas [
25] and by Schmidt [
28] were made for the refrigerant- and air-side, respectively. The equations below illustrate the Dittus–Boelter’s [
27] and Schmidt’s [
28] correlations, respectively:
where
n = 0.3 for the fluid being cooled.
where
C = 0.3 for in-line arrangement.
2.2. Experimental Set Up
As it is referred, the U-value consists of three parts: the air convection, refrigerant convection, and the conduction. In the present case, the air-side heat transfer represents the main thermal resistance. Thus, the experimental set up aimed to identify a suitable heat transfer correlation with focus on the air side. In order to validate the calculations, especially the air-side heat transfer, a well-known working medium should be used for the experiments on the refrigerant-side. Here, water with well-known thermophysical properties and reliable heat transfer correlations at single-phase regime was selected as a working medium.
The test rig consisted of the heater, the gas cooler, the measurement equipment, and controls. To enable the information to be read and recorded, the instrumentations were connected to a data logging system. The test rig is shown in
Figure 2 below.
The air-cooled HEX used for the experiments was the CU-713CX2 from Panasonic. The heater was from the Single® company, model STW 150/1-18-45-KS7. The K-type thermocouples and pressure transducers used were from OMEGA company with uncertainties of ±0.2 °C and ±1.5%, respectively, while the mass flow valve and meter with uncertainties of ±0.5% were manufactured by Bürkert.
The experiments were carried out for different water mass flows, fan frequencies, and inlet water temperatures, as the following
Table 2 shows.
Figure 3a,b illustrate the air mass flow rate and fan power as function of the fan frequency.
Using Equations (20)–(24), the experimental overall heat transfer coefficient was calculated by:
2.3. Model Application
The validated model was finally applied to define an efficient and reliable gas cooler design, as well as to evaluate its performance for different off-design conditions. Within scope, the on-design analysis investigated different designs of heat exchangers, such as the number of rows and length of the tube. The target of the off-design analysis was to evaluate the operation of the gas cooler in different ambient temperatures.
2.3.1. On-Design
The developed model was utilized to design an air-finned CO2 gas cooler. Therefore, three and four numbers of rows and finned tubes with inner diameters of 6.85, 16, and 22.5 mm were considered. In order to identify an efficient gas cooler design, a compromise between the bundle area and the refrigerant pressure drop of the gas cooler was aimed for. The model was provided with the inlet, outlet temperatures and pressures, medium mass flows, and the duty of the heat exchanger as input variables. In addition, the design properties of the tube were specified, so the total area of the tube was calculated. The air and refrigerant heat transfer coefficients, which were initially unknown, were assumed, and then the overall heat transfer coefficient was calculated from Equation (1). The required exchanged area was calculated from Equation (2), and so the required number of tubes was obtained. In the iteration process, the air and refrigerant heat transfer coefficients were updated using Equations (3)–(17). The iteration was continued until the relative tolerance of the overall heat transfer coefficient for two continuous iterations was equal to 0.001.
The pressure drop was calculated by the following equation:
which Filonenko [
27] applies for
.
In case of
, Blasius [
27] correlations was applied:
The thermophysical properties of air and refrigerant like density, viscosity, specific heat capacity, and thermal conductivity were obtained from REFPROP version 10 database. The model used the mean temperature and pressure of the mediums in order to calculate the heat transfer coefficients.
2.3.2. Off-Design
The calculations were used to evaluate the overall performance of the gas cooler in different conditions. The defined gas cooler was investigated in different ambient temperatures. The boundary conditions and operational parameters for the off-design analysis were defined according to
Table 3.
Based on the off-design data, the overall heat transfer coefficient was calculated using Equations (3)–(17). In addition, investigation of different potential heat exchanger designs was conducted by comparing their overall heat transfer coefficient and the refrigerant pressure drop. The potential heat exchanger’s designs were based on the on-design analysis, and were chosen in terms of pressure losses and bundle area. The selected air-cooled HEXs were designed with identical finned-tube properties.
4. Conclusions
A design procedure for a CO2 finned-tube gas cooler was developed and validated experimentally. The experimental focus was laid on the validation of the air-side heat transfer. Absolute deviations between the model and the experiments were extracted and prove the reliability of the selected heat transfer correlations. The developed simulation model was used to design an efficient gas cooler and evaluate its performance in different ambient temperatures. Based on these results, the following conclusions are justified:
The deviations between the calculations and the experiments are highly affected by the fan frequency, the water mass flow, and the water inlet temperature. It was found that the absolute average deviations for 60–80 Hz are less than 10%. The increase in the fan frequency, the water mass flow, and the water inlet temperature caused an improvement of the overall heat transfer coefficient of the heat exchanger. The comparison between potential heat transfer correlations showed that the combination of the correlations proposed by VDI-Heat Atlas [
25] and by Gnielinski [
26] for the air- and refrigerant-side heat transfer coefficient, respectively, approached the experimental results better.
The heat exchanger with four rows has a smaller bundle area than with three rows, while the investigation of different tubes showed that the heat exchanger designed with the smallest diameter has the smallest bundle area and the highest refrigerant pressure drop. As a compromise between the bundle area and the refrigerant pressure drop, a gas cooler of 2.11 m2 and refrigerant pressure drop of 0.23 bar was defined. Comparison between potential gas coolers was made, and the results show that the defined gas cooler operates more efficiently for different ambient conditions, compared to other potential heat exchangers.
Due to the pressure limitations of the equipment, the experimental validation of the model was conducted with water as working fluid. This procedure enables a reliable validation of the applied heat transfer correlation at the air side. However, the use of CO2 as working fluid in the tubes would improve validation approach. Therefore, experiments with CO2 are suggested for further work, next to the off-design analysis for the entire air conditioning system including the developed gas cooler model.