Analysis of the Effect of the CaCl 2 Mass Fraction on the Efficiency of a Heat Pump Integrated Heat-Source Tower Using an Artificial Neural Network Model

An existing idle cooling tower can be reversibly used as a heat-source tower (HST) to drive a heat pump (HP) in cold seasons, with calcium chloride (CaCl2) aqueous solution commonly selected as the secondary working fluid in an indirect system due to its good thermo-physical properties. This study analyzed the effect of CaCl2 mass fraction on the effectiveness (ε) of a closed HST and the coefficient of performance (COP) of a HP heating system using an artificial neural network (ANN) technique. CaCl2 aqueous solutions with five different mass fractions, viz. 3%, 9%, 15%, 21%, and 27%, were chosen as the secondary working fluids for the HSTHP heating system. In order to collect enough measured data, extensive field tests were conducted on an experimental test rig in Changsha, China which experiences hot summer and cold winter weather. After back-propagation (BP) training, the three-layer (4-9-2) ANN model with a tangent sigmoid transfer function at the hidden layer and a linear transfer function at the output layer was developed for predicting the tower effectiveness and the COP of the HP under different inlet air dry-/wet-bulb temperatures, hot water inlet temperatures and CaCl2 mass fractions. The correlation coefficient (R), mean relative error (MRE) and root mean squared error (RMSE) were adopted to evaluate the prediction accuracy of the ANN model. The results showed that the R, MRE, and RMSE between the training values and the experimental values of ε (COP) were 0.995 (0.996), 2.09% (1.89%), and 0.005 (0.060), respectively, which indicated that the ANN model was reliable and robust in predicting the performance of the HP. The findings of this paper indicated that in order to guarantee normal operation of the system, the freezing point temperature of the CaCl2 aqueous solution should be sufficiently (3–5 K) below its lowest operating temperature or lower than the normal operating temperature by about 10 K. The tower effectiveness increased with increasing CaCl2 mass fraction from 0 to 27%, while the COP of the HP decreased. A tradeoff between the tower effectiveness and the COP of the HP should be considered to further determine the suitable mass fraction of CaCl2 aqueous solution for the HSTHP heating system. The outputs of this study are expected to provide guidelines for selecting brine with an appropriate mass fraction for a closed HSTHP heating system for actual applications, which would be a reasonable solution to improve the system performance.


Introduction
A mechanical draft cooling tower is usually employed to release the heat withdrawn from buildings [1][2][3][4], but it also can be fully and reversibly used as a heat-source-tower (HST) to drive a heat pump (HP) heating system in cold seasons [5]. Heat-source tower heat pumps (HSTHPs), having the same advantage as common water-cooling air conditioners and water source heat pumps, can utilize sensible and latent heat of air to efficiently carry out three-combined supply of air-conditioning, heating, and hot water for buildings [6]. Recently, a number of HSTHP systems were installed in China, and have been operating satisfactorily for several years [7,8].
Many researchers have carried out a series of experimental studies and theoretical analyses on the performance characteristics of HST. Liang et al. [9] studied the various rules of solution regeneration rate of an open HST installed in a sub-tropical region of China via experimental and theoretical analysis. It was found that increasing the auxiliary heat greatly improved the solution regeneration rate. Wen et al. [10] conducted an experimental study on the heat transfer coefficient between air and liquid, such as water and ethylene glycol, with 40% mass concentration in a cross-flow HST. An empirical correlation of the heat transfer coefficient was developed by a regression method. Cheng et al. [11,12] analyzed the merits and demerits of closed HSTHP systems in hot-summer and cold-winter regions of China. It was reported that HSTHP can still maintain high energy efficiency under frost prevention conditions. Zhang et al. [13] developed an analytical model for the coupled heat and mass transfer processes in a counter flow HST based on operating conditions. Wu et al. [14] conducted an experimental investigation on the influence factors of the Lewis number in a cross flow HST and predicted the performance characteristics of the HST heat pump heating system in winter using an artificial neural network (ANN) technique [15] and an adaptive neuro-fuzzy inference approach [16]. Tan and Deng [17] found that the use of an HST would achieve higher energy efficiency than the use of electrical heating as the backup heat provisions when the building space cooling load is reduced. Subsequently, a method was developed for evaluating the heat and mass transfer characteristics in an HST for heat recovery [18]. A complete method was also provided for analyzing the heat and mass transfer characteristics within and at the boundaries of the HST together with the analytical method developed previously [19].
Although the performance characteristics of cooling towers integrated with cooling water systems have been widely reported in the literature, there are few researches focusing on HSTHP systems [20][21][22][23]. The operational characteristics of HSTHP systems are expected to be significantly different from that of common cooling water systems. These differences mainly include reduction of latent heat exchange and the relatively low secondary working fluid temperature, which is usually below 0˝C. Brine antifreezes instead of water are usually selected as the circulation medium due to their good thermo-physical properties. From the present research status, very few researchers have studied the calcium chloride (CaCl 2 ) mass fraction influence on the efficiency of an HSTHP heating system. In fact, this influence is non-linear, and the factors are interactional and inter-coupling. ANN modeling is good at preferably catching nonlinear rule instead of classical mathematic modeling. Meanwhile, ANN models were successfully used to predict the cooling towers' performance [24][25][26]. In this paper, ANN technology is also adopted to simulate the performance of the HSTHP system.
This study aims to analyze the effect of the CaCl 2 mass fraction on the effectiveness (ε) of an HST and the coefficient of performance (COP) of an HP by means of experimental testing and ANN technology. A series of experiments were conducted on a test rig in Hunan University of China, which experiences hot summer and cold winter weather. CaCl 2 aqueous solutions with five different mass fractions, viz. 3%, 9%, 15%, 21%, and 27%, were chosen as the secondary working fluids in this study. Extensive measured data, including inlet air dry-/wet-bulb temperatures and hot water inlet/outlet temperatures, were obtained under different CaCl 2 mass fractions and constant input powers of fan, circulation pump, and compressor. The ANN approach was adopted to deal with the experimental data. In the ANN model, inlet air dry-/wet-bulb temperatures, hot water inlet temperature, and CaCl 2 mass fraction were input parameters, and the tower effectiveness and the COP of the HP unit were output parameters. The findings should facilitate engineers to select brine antifreezes with a suitable mass fraction for an HSTHP system.

Description of the Closed HSTHP System
An experimental closed HSTHP system [12] was installed in Hunan University, China, which experiences hot summer and cold winter weather. Extensive field experimental work on the performance characteristics of the HSTHP system was carried out in the experimental test rig. The tested HST, which was an induced draft counter-flow type and specially designed, was utilized fully and efficiently as a heat source to extract heat from ambient air to drive an HP heating system. The HP system was an integrated water source HP unit with a rated power of 4.1 kW and a total heating capacity of 16 kW. Figure 1 presents the schematic diagram of the experimental system. The heating system consisted of five working fluids, viz. moist air, antifreeze solution, circulating solution, refrigerant (R22), and water. CaCl 2 aqueous solution, which was an alternative fluid for application in low temperature conditions (below 0˝C), was selected as the secondary working fluid due to its good thermo-physical properties [27].

Description of the Closed HSTHP System
An experimental closed HSTHP system [12] was installed in Hunan University, China, which experiences hot summer and cold winter weather. Extensive field experimental work on the performance characteristics of the HSTHP system was carried out in the experimental test rig. The tested HST, which was an induced draft counter-flow type and specially designed, was utilized fully and efficiently as a heat source to extract heat from ambient air to drive an HP heating system. The HP system was an integrated water source HP unit with a rated power of 4.1 kW and a total heating capacity of 16 kW. Figure 1 presents the schematic diagram of the experimental system. The heating system consisted of five working fluids, viz. moist air, antifreeze solution, circulating solution, refrigerant (R22), and water. CaCl2 aqueous solution, which was an alternative fluid for application in low temperature conditions (below 0 °C), was selected as the secondary working fluid due to its good thermo-physical properties [27]. Ambient air is sucked into the HST unit by the fan with a rated power of 0.55 kW, and then indirectly cooled by the CaCl2 aqueous solution in the finned heat exchange coil inside the HST. Under frosting conditions, the antifreeze solution, sprayed on the top of the finned coil, can protect the surface of the heat exchanger from frosting. The cold aqueous solution is pumped into the HST by a circulating solution pump with a rated power of 0.75 kW to absorb heat from the ambient air, and then pumped back into the evaporator of the HP, to exchange heat with the refrigerant of R22. The water entering the condenser is heated at the same time for a hot water supply. Therefore, the HST unit, replacing a boiler or other heating systems, can provide a low temperature heat source for HP units for hot water production. Fan, antifreeze solution pump, circulating solution pump, and compressor are all controlled by frequency converter speed control systems in order to improve the system's operation stability as well as to save energy.

System Performance Indexes
The tower effectiveness and the COP of the HP unit were selected as the indexes to evaluate the system performance under non-frosting operating conditions and constant input power of fan, circulating pump, and compressor. Similar to the definition of cooling tower effectiveness, the effectiveness of HST can be defined as the ratio of the actual inlet temperature rise of the CaCl2 aqueous solution to the ideal maximum temperature rise. Ambient air is sucked into the HST unit by the fan with a rated power of 0.55 kW, and then indirectly cooled by the CaCl 2 aqueous solution in the finned heat exchange coil inside the HST. Under frosting conditions, the antifreeze solution, sprayed on the top of the finned coil, can protect the surface of the heat exchanger from frosting. The cold aqueous solution is pumped into the HST by a circulating solution pump with a rated power of 0.75 kW to absorb heat from the ambient air, and then pumped back into the evaporator of the HP, to exchange heat with the refrigerant of R22. The water entering the condenser is heated at the same time for a hot water supply. Therefore, the HST unit, replacing a boiler or other heating systems, can provide a low temperature heat source for HP units for hot water production. Fan, antifreeze solution pump, circulating solution pump, and compressor are all controlled by frequency converter speed control systems in order to improve the system's operation stability as well as to save energy.

System Performance Indexes
The tower effectiveness and the COP of the HP unit were selected as the indexes to evaluate the system performance under non-frosting operating conditions and constant input power of fan, circulating pump, and compressor. Similar to the definition of cooling tower effectiveness, the effectiveness of HST can be defined as the ratio of the actual inlet temperature rise of the CaCl 2 aqueous solution to the ideal maximum temperature rise.
The COP of HP is expressed as follows: COP " c w m w pT w,o´Tw,i q P comp (2) The correlation equation between these two indexes was established as follows: εc aq ρ aq v aq pT wb,i´Taq,i q " pCOP´ηqˆP comp As shown in the correlation equation above, the relationship between the tower effectiveness and the COP of the HP is non-linear. The influence factors included the specific heat capacity, the density, the volume flow rate of CaCl 2 aqueous solution, the CaCl 2 solution temperature and the ambient wet-bulb temperature at the inlet of the HST, as well as the input power and efficiency of the compressor.

Testing Procedure
A set of instruments were installed to measure the related operating parameters, including ambient air dry-/wet-bulb temperatures, dry-/wet-bulb temperatures of the air outlet, air velocity, flow rates of the aqueous solution and the hot water, inlet and outlet temperatures of the aqueous solutions, as well as the input power of fan, antifreeze solution pump, circulating solution pump and compressor. The detailed arrangement of the instruments is presented in Figure 1. Table 1 shows the specification and accuracy of the instruments. Table 1. Specifications of instruments. Figure 1)

Instruments (Type) Accuracy
Dry-and wet-bulb temperature of inlet (T1) and outlet (T2) air Platinum resistance thermometer (PT100)˘0 .2˝C Discharge temperature of condensate water(T3) Inlet (T4) and outlet (T5) aqueous solution temperature Inlet (T6) and outlet (T7) water temperatures In this study, the input power of the fan, circulating solution pump, and the compressor were kept constant and were 0.51 kW, 0.75 kW, and 4.1 kW, respectively. During the whole experimental process, the input power of the antifreeze solution pump was zero because the closed HST was operated under non-frosting conditions. The hot water flow rate was 1.31 m 3 /h. CaCl 2 aqueous solutions with five different mass fractions, viz. 3%, 9%, 15%, 21%, and 27%, were chosen as the secondary working fluid, respectively. The other relevant parameters were recorded by color paperless recorders every five minutes after the system achieved stability for an hour. Each set of data was chosen as the system achieved steady state, which required that the temperatures at all test points fluctuated within˘0.2˝C for longer than 20 min.

Uncertainty Analysis
Uncertainty analyses for the experimental results were expressed by the following equation [28]: Obviously, the uncertainty of the tower effectiveness was influenced by the measurement accuracy of the ambient wet-bulb temperature, and the CaCl 2 aqueous solution temperatures at the inlet and outlet of the HST. The uncertainty analyses for the COP focused on the temperature difference between the inlet and outlets of the hot water, the input power of the compressor, and the hot water mass flow rate. For example, when T w,o = 42.8˝C and T w,i = 33.6˝C, the uncertainty of the COP, calculated based on the given operating parameters, was 4.7%, in which the uncertainties caused by the temperature difference between the inlet and outlets of the hot water and the input power of the compressor accounted for 85.3% and 10.2% of the total uncertainty, respectively, while the uncertainty caused by the hot water mass flow rate was only 0.5%. Therefore, the key point of improving the test accuracy is to improve the temperature measurement accuracy.

Three-Layer BP Network
The HSTHP is a complex system, mainly consisting of the processes of heat and mass transfer, electro-thermal conversion, and electro-kinetic mass transport. It is quite difficult to accurately describe these processes using classical modelling techniques. Thus, an ANN model, which can be used to resolve complicated and non-linear problems of internal heat transfer, was adopted in this paper [15,25,[29][30][31]. The model with one hidden layer can meet the simulative requirements [32]. Compared with other ANN models, it is quite appropriate by using the BP network to solve complicated problems. Thus, a three-layer BP network with a complicated non-linear mapping function was adopted to realize the mapping from input to output. As shown in Figure 2, the input layer has four nodes, including inlet air dry-bulb temperature T a,i , inlet air wet-bulb temperature T wb,i , hot water inlet temperature T w,i , and the mass concentration of CaCl 2 solution w. Also the output layer has two nodes, viz. the tower effectiveness and the COP of the HP. In addition, the value ranges of the input parameters are shown in Table 2.
Obviously, the uncertainty of the tower effectiveness was influenced by the measurement accuracy of the ambient wet-bulb temperature, and the CaCl2 aqueous solution temperatures at the inlet and outlet of the HST. The uncertainty analyses for the COP focused on the temperature difference between the inlet and outlets of the hot water, the input power of the compressor, and the hot water mass flow rate. For example, when Tw,o = 42.8 °C and Tw,i = 33.6 °C, the uncertainty of the COP, calculated based on the given operating parameters, was 4.7%, in which the uncertainties caused by the temperature difference between the inlet and outlets of the hot water and the input power of the compressor accounted for 85.3% and 10.2% of the total uncertainty, respectively, while the uncertainty caused by the hot water mass flow rate was only 0.5%. Therefore, the key point of improving the test accuracy is to improve the temperature measurement accuracy.

Three-Layer BP Network
The HSTHP is a complex system, mainly consisting of the processes of heat and mass transfer, electro-thermal conversion, and electro-kinetic mass transport. It is quite difficult to accurately describe these processes using classical modelling techniques. Thus, an ANN model, which can be used to resolve complicated and non-linear problems of internal heat transfer, was adopted in this paper [15,25,[29][30][31]. The model with one hidden layer can meet the simulative requirements [32]. Compared with other ANN models, it is quite appropriate by using the BP network to solve complicated problems. Thus, a three-layer BP network with a complicated non-linear mapping function was adopted to realize the mapping from input to output. As shown in Figure 2, the input layer has four nodes, including inlet air dry-bulb temperature Ta,i, inlet air wet-bulb temperature Twb,i, hot water inlet temperature Tw,i, and the mass concentration of CaCl2 solution w. Also the output layer has two nodes, viz. the tower effectiveness and the COP of the HP. In addition, the value ranges of the input parameters are shown in Table 2.    The BP network part was implemented under the Matlab environment, and the tangent sigmoid function and the linear transfer function were chosen at the hidden layer and the output layer, respectively. The data set consisted of 319 input-output pairs, therein, 70% of the data set were randomly assigned as the training set, the remaining 30% were employed for testing the network. All the data were normalized into´1 to 1 in order to improve the prediction agreement. The normalization and anti-normalization functions for the training data were premnmx and postmnmx functions, and the normalization and anti-normalization functions for testing data were tramnmx and postmnx functions, respectively. During the course of training, the traingdx function was used to act as the training function in this paper, because it has higher stability and a faster convergence rate [33].

Performance Analysis of ANN
In order to assess the accuracy of the BP network model, the mean relative error (MRE), correlation coefficient (R) and root mean square error (RMSE) were used as the characteristic parameters to evaluate the consistency between training and prediction results [15,25,34].
The MRE, which shows the mean ratio between the error and the network output values, is given by [15,25,34] MREp%q 00ˆa i´bi a iˇ (  5) The correlation coefficient (R) is a measurement of how well the variation in the predicted outputs is explained by the experimental values, and the R value between the experimental values and the predicted outputs is given by [15,25,34], where cov(a,b) is the covariance between a and p sets which represent the experimental and network predicted output sets, respectively, and is given by In addition, cov(a,a) and cov(b,b) are the auto covariance of a and b sets, respectively, and are expressed by: covpa, aq " Epa´µ a q 2 (8) covpb, bq " Epb´µ b q 2 The correlation coefficient values close to +1 indicate a stronger agreement between training and predicted values, while the values close to´1 indicate a stronger negative relationship between training and predicted values.
The root mean square error (RMSE), which could reveal the accuracy of the model, is calculated as: [15,34],

Determination of the Node Number in the Hidden Layer
Three formulas, which can receive the node number in the hidden layer, can be determined by Equations (10)-(12) [35], n 1 " ? n`m`c 0 (12) For the three-layer BP network, the calculated formula of the node number in the hidden layer is defined by [33], Another empirical formula [36] was also developed to determine the node number in the hidden layer, which is expressed as follows: According to the above formulas, the node number ranged from 2-13 for the hidden layer. Figure 3 illustrates the dependence between the neuron number and mean squared error (MSE) for the variable learning rate backpropagation algorithm. As shown in Figure 3, the MSE of the network is much higher for the second (MSE 0.0155), the third (MSE 0.0070), and the fourth (MSE 0.0052) hidden neurons than for the others. For the other neurons in the hidden layer, the corresponding MSE fluctuates around 0.001. With nine hidden neurons, the MSE reached its minimum value of 0.0007. Hence, the neural network containing nine hidden neurons was chosen as the best case. log n n  (13) For the three-layer BP network, the calculated formula of the node number in the hidden layer is defined by [33], 1 n nm  (14) Another empirical formula [36] was also developed to determine the node number in the hidden layer, which is expressed as follows: According to the above formulas, the node number ranged from 2-13 for the hidden layer. Figure 3 illustrates the dependence between the neuron number and mean squared error (MSE) for the variable learning rate backpropagation algorithm. As shown in Figure 3, the MSE of the network is much higher for the second (MSE 0.0155), the third (MSE 0.0070), and the fourth (MSE 0.0052) hidden neurons than for the others. For the other neurons in the hidden layer, the corresponding MSE fluctuates around 0.001. With nine hidden neurons, the MSE reached its minimum value of 0.0007. Hence, the neural network containing nine hidden neurons was chosen as the best case.  Figure 4 illustrates the training process, validation, and test mean squared errors for the variable learning rate backpropagation algorithm. The best validation performance is about 0.0007 at 131 iterations, after that the training was stopped because the training error did not change. These results indicate that the ANN network model is reliable.

Thermo-Physical Properties of CaCl2 Aqueous Solution
Water is an excellent secondary working fluid for air conditioning when the lowest temperature is not less than +3 °C, while CaCl2 aqueous solution is an alternative for applications below 0 °C due to the freezing point depressor and the high water affinity [27]. CaCl2 aqueous solution was earlier used as a brine in various air conditioning systems. In this work, CaCl2 aqueous solutions with different mass fractions, viz. 3%, 9%, 15%, 21%, and 27%, were chosen as the secondary working fluids in the experiments. The mass fractions and the corresponding temperatures of these aqueous solutions could be used to calculate their thermo-physical properties according to the methods described by Melinder [27] and Conder [37]. Figures 5-8 show the basic thermo-physical properties, viz. freezing point, density, viscosity, specific heat, and thermal conductivity. When the mass fraction varied from 0%-27% (0% ≤ w ≤ 27%), the freezing point and the relative density (to saturated liquid water at the same temperature) may be represented by a single function of the mass fraction, as shown in Figure 5. At the freezing point temperature, ice crystals begin to form in equilibrium if there is no sub-cooling. The relative density increases with increasing mass fraction. The specific heat capacity, dynamic viscosity, and thermal conductivity of CaCl2 aqueous solution may be represented by a function of the mass fraction and the temperature. Figures 6 and 7 show that the specific heat capacity and the thermal conductivity increase with increasing temperature and decrease with increasing mass fraction. As shown in Figure 8, the dynamic viscosity increases with increasing mass fraction and decreases with increasing temperature. Overall, CaCl2 aqueous solution, which provides an alternative option for application below 0 °C, could give sufficient freezing protection, small pressure drop, and good heat transfer performance.

Thermo-Physical Properties of CaCl 2 Aqueous Solution
Water is an excellent secondary working fluid for air conditioning when the lowest temperature is not less than +3˝C, while CaCl 2 aqueous solution is an alternative for applications below 0˝C due to the freezing point depressor and the high water affinity [27]. CaCl 2 aqueous solution was earlier used as a brine in various air conditioning systems. In this work, CaCl 2 aqueous solutions with different mass fractions, viz. 3%, 9%, 15%, 21%, and 27%, were chosen as the secondary working fluids in the experiments. The mass fractions and the corresponding temperatures of these aqueous solutions could be used to calculate their thermo-physical properties according to the methods described by Melinder [27] and Conder [37]. Figures 5-8 show the basic thermo-physical properties, viz. freezing point, density, viscosity, specific heat, and thermal conductivity. When the mass fraction varied from 0%-27% (0% ď w ď 27%), the freezing point and the relative density (to saturated liquid water at the same temperature) may be represented by a single function of the mass fraction, as shown in Figure 5. At the freezing point temperature, ice crystals begin to form in equilibrium if there is no sub-cooling. The relative density increases with increasing mass fraction. The specific heat capacity, dynamic viscosity, and thermal conductivity of CaCl 2 aqueous solution may be represented by a function of the mass fraction and the temperature. Figures 6 and 7 show that the specific heat capacity and the thermal conductivity increase with increasing temperature and decrease with increasing mass fraction. As shown in Figure 8, the dynamic viscosity increases with increasing mass fraction and decreases with increasing temperature. Overall, CaCl 2 aqueous solution, which provides an alternative option for application below 0˝C, could give sufficient freezing protection, small pressure drop, and good heat transfer performance.

Training Results and Performance Analysis of the ANN Model
Over many trials, the MRE reached the minimum, the R of training results was very close to +1 and the results had the least RMSE when the node number in the hidden layer was 9. As shown in Figure 9a, the R, MRE, and RMSE between training values and experimental values for ε are 0.995%, 2.09%, and 0.005%, respectively, while Figure 9b shows those for COP are 0.996%, 1.89%, and 0.060%, respectively. These results demonstrate that the predicted values of ε and COP are quite accurate and the training results have a very good agreement according to the characteristic parameters MRE, R, and RMSE. Thus, the ANN model can effectively predict the tower effectiveness and the COP of the HP under normal operating conditions.

Effects of CaCl2 Mass Fraction on the Efficiency of the HSTHP
In this paper, we mainly focus on studying the effect of CaCl2 mass fraction on the tower effectiveness and the COP of the HP. In order to visualize the effect, the ε and COP predicted by the ANN model under different CaCl2 mass fractions are indicated in Figure 10. This figure depicts the same variation trends of the predicted and the experimental values with the CaCl2 mass fraction changing from 0% to 27%, and the other three input parameters kept constant (viz. the inlet temperature of the hot water, inlet air dry-and wet-bulb temperatures are equal to 35.0 °C, 10.0 °C,

Training Results and Performance Analysis of the ANN Model
Over many trials, the MRE reached the minimum, the R of training results was very close to +1 and the results had the least RMSE when the node number in the hidden layer was 9. As shown in Figure 9a, the R, MRE, and RMSE between training values and experimental values for ε are 0.995%, 2.09%, and 0.005%, respectively, while Figure 9b shows those for COP are 0.996%, 1.89%, and 0.060%, respectively. These results demonstrate that the predicted values of ε and COP are quite accurate and the training results have a very good agreement according to the characteristic parameters MRE, R, and RMSE. Thus, the ANN model can effectively predict the tower effectiveness and the COP of the HP under normal operating conditions.

Training Results and Performance Analysis of the ANN Model
Over many trials, the MRE reached the minimum, the R of training results was very close to +1 and the results had the least RMSE when the node number in the hidden layer was 9. As shown in Figure 9a, the R, MRE, and RMSE between training values and experimental values for ε are 0.995%, 2.09%, and 0.005%, respectively, while Figure 9b shows those for COP are 0.996%, 1.89%, and 0.060%, respectively. These results demonstrate that the predicted values of ε and COP are quite accurate and the training results have a very good agreement according to the characteristic parameters MRE, R, and RMSE. Thus, the ANN model can effectively predict the tower effectiveness and the COP of the HP under normal operating conditions.

Effects of CaCl2 Mass Fraction on the Efficiency of the HSTHP
In this paper, we mainly focus on studying the effect of CaCl2 mass fraction on the tower effectiveness and the COP of the HP. In order to visualize the effect, the ε and COP predicted by the ANN model under different CaCl2 mass fractions are indicated in Figure 10. This figure depicts the same variation trends of the predicted and the experimental values with the CaCl2 mass fraction changing from 0% to 27%, and the other three input parameters kept constant (viz. the inlet temperature of the hot water, inlet air dry-and wet-bulb temperatures are equal to 35.0 °C, 10.0 °C,

Effects of CaCl 2 Mass Fraction on the Efficiency of the HSTHP
In this paper, we mainly focus on studying the effect of CaCl 2 mass fraction on the tower effectiveness and the COP of the HP. In order to visualize the effect, the ε and COP predicted by the ANN model under different CaCl 2 mass fractions are indicated in Figure 10. This figure depicts the same variation trends of the predicted and the experimental values with the CaCl 2 mass fraction changing from 0% to 27%, and the other three input parameters kept constant (viz. the inlet temperature of the hot water, inlet air dry-and wet-bulb temperatures are equal to 35.0˝C, 10.0˝C, and 7.0˝C, respectively). It can be seen that the predicted values are in good agreement with the experimental values. The tower effectiveness increases with increasing CaCl 2 mass fraction. This may be explained that with increasing CaCl 2 mass fraction, the thermal conductivity of the CaCl 2 aqueous solution decreases which in turn gives a negative effect to the heat transfer in the HST. Because the density and the dynamic viscosity of the CaCl 2 aqueous solution increase with increasing CaCl 2 mass fraction from 0% to 27%, a decline in the mass flow rate of the CaCl 2 aqueous solution can lead to provision of adequate retention time of the heat transfer fluid to ensure positive heat transfer. Meanwhile, the decreasing specific heat capacity requires less thermal energy for each degree rising. Therefore, the actual inlet temperature rise rate of CaCl 2 aqueous solution is faster than that of the ideal maximum. The COP of the HP decreases with increasing CaCl 2 mass fraction. The results indicate that the lower the CaCl 2 mass fraction, the higher the COP of the HP. The reasons can be explained in that the heat absorbed by the HST and transported by the pipeline network decreases with increasing CaCl 2 mass fraction due to the decrease of the heat transfer coefficient and the CaCl 2 aqueous solution flow rate.  A previous study [27] showed that for a suitable secondary working fluid selected in an indirect system, its freezing point temperature should be lower than the lowest expected temperature and the normal operating temperature by about 10 K such that the fluid can be pumped through the system easily. Another suggestion given in ref. [38] was that the freezing point of the secondary fluid should be sufficiently (3-5 K) below its lowest operating temperature. The main findings of this paper are that the tower effectiveness increased with increasing CaCl2 mass fraction from 0% to 27%, while the COP of the heat pump decreased. Therefore, a tradeoff between the tower effectiveness and the COP of the heat pump should be considered to further determine the suitable mass fraction of CaCl2 aqueous solution for the HSTHP heating system.

Conclusions
This study analyzed the effect of CaCl2 mass fraction on the COP of a heat pump (HP) and the effectiveness of a closed heat-source tower (HST) using artificial neural network (ANN) technology. An ANN model based on a variable learning rate back-propagation algorithm was developed. Sufficient field experimental data were obtained from an existing heat-source tower heat pump (HSTHP) heating system under non-frosting operating conditions. After back-propagation training with the data, the three-layer (4-9-2) ANN model with a tangent sigmoid transfer function (tansig) at the hidden layer and a linear transfer function (purelin) at the output layer was used to efficiently predict the tower effectiveness and the COP of the HP.
Correlation coefficient (R), mean relative error (MRE) and root mean squared error (RMSE) were adopted to evaluate the prediction accuracy of the ANN model. The results demonstrated that the R, MRE, and RMSE between training values and experimental values for ε (COP) were 0.995 A previous study [27] showed that for a suitable secondary working fluid selected in an indirect system, its freezing point temperature should be lower than the lowest expected temperature and the normal operating temperature by about 10 K such that the fluid can be pumped through the system easily. Another suggestion given in ref. [38] was that the freezing point of the secondary fluid should be sufficiently (3-5 K) below its lowest operating temperature. The main findings of this paper are that the tower effectiveness increased with increasing CaCl 2 mass fraction from 0% to 27%, while the COP of the heat pump decreased. Therefore, a tradeoff between the tower effectiveness and the COP of the heat pump should be considered to further determine the suitable mass fraction of CaCl 2 aqueous solution for the HSTHP heating system.

Conclusions
This study analyzed the effect of CaCl 2 mass fraction on the COP of a heat pump (HP) and the effectiveness of a closed heat-source tower (HST) using artificial neural network (ANN) technology. An ANN model based on a variable learning rate back-propagation algorithm was developed. Sufficient field experimental data were obtained from an existing heat-source tower heat pump (HSTHP) heating system under non-frosting operating conditions. After back-propagation training with the data, the three-layer (4-9-2) ANN model with a tangent sigmoid transfer function (tansig) at the hidden layer and a linear transfer function (purelin) at the output layer was used to efficiently predict the tower effectiveness and the COP of the HP.
Correlation coefficient (R), mean relative error (MRE) and root mean squared error (RMSE) were adopted to evaluate the prediction accuracy of the ANN model. The results demonstrated that the R, MRE, and RMSE between training values and experimental values for ε (COP) were 0.995 (0.996), 2.09% (1.89%), and 0.005 (0.060), respectively, which indicated that the ANN model was reliable and could be used to predict both performance indexes of the HSTHP system with a high degree of accuracy.
Results showed that in order to guarantee the normal operation of the system, the freezing point temperature of the CaCl 2 aqueous solution should be sufficiently (3-5 K) below its lowest operating temperature or lower than the normal operating temperature by about 10 K. The tower effectiveness increased with increasing CaCl 2 mass fraction from 0% to 27%, while the COP of the HP decreased. Therefore, a tradeoff between the tower effectiveness and the COP of the HP should be considered to further determine the appropriate CaCl 2 mass fraction for an HSTHP heating system. The conclusions are expected to provide guidelines for selecting brine with a suitable mass fraction in actual applications, which is a reasonable solution to improve the efficiency of a closed HSTHP system.