Study on the performance of a solar assisted air source heat pump system for building heating
Cai-hua Liang,Xiao-song Zhang,Xiu-wei Li,Xia Zhu
DOI: https://doi.org/10.1016/j.enbuild.2011.04.028
IF: 7.201
2011-01-01
Energy and Buildings
Abstract:This paper proposed a new solar assisted air source heat pump system with flexible operational modes to improve the performance of the heating system. A mathematical model was established on the solar assisted air source heat pump system for building heating with a heating capacity of 10 kW, and an air source heat pump unit was developed to validate the model. The effect of the solar collector area on the performance of the system running in Nanjing was studied. The results showed that the COP of the heat pump unit was enhanced with the increase of the solar radiation density during the typical sunny day in the heating season. In addition, the COP also increased in proportion to the solar collector area. Compared with the case when the solar collector area was 0 m 2 , the COP increase of the heat pump and the energy-saving rate were 11.22% and 24% respectively when the solar collector area was 40 m 2 . Meanwhile, the solar equivalent generation power efficiency could reach 11.8%. Keywords Building heating system Solar energy Air source heat pump Solar equivalent generation power efficiency Nomenclature A v Circulation area of electronic expansion valve (m 2 ) C v Flow coefficient of electronic expansion valve (dimensionless) d b Diameter of root of fins (m) d i Internal diameter of pipe (m) d o External diameter of pipe (m) f Impact factor (dimensionless) G com Flow of refrigerant in compressor (kg/s) h Enthalpy (kJ/kg) I Solar radiation intensity (W/m 2 ) m v Flow of refrigerant in electronic expansion valve (kg/s) N Row number of finned tube (dimensionless) n Rotational speed of compressor (r/min) Nu Nusselt number (dimensionless) p Pressure (Pa) P cr Critical pressure of refrigerant (Pa) Pr Prandtl number (dimensionless) q Heating capacity (kW) Re Reynolds number (dimensionless) s Fin interval (m) s 2 Interval between tubes along the air direction (m) T a Environment temperature (°C) T i Hot water temperature inside collector tube (°C) V th Theoretical capacity of compressor (m 3 ) V com Working capacity of compressor (m 3 ) W Power (kW) X tt Martinelli Number (dimensionless) x Dryness of refrigerant (dimensionless) Greek symbols a Heat transfer coefficient (W/(m 2 °C)) λ Heat conduction coefficient (W/m K) μ Dynamic viscosity of refrigerant (Pa s) ρ Intensity (kg/m 3 ) η Instant collecting efficiency of solar collector (dimensionless) η com Capacity efficiency of compressor (dimensionless) η k Shaft efficiency of compressor (dimensionless) ν Absorption specific volume of compressor (m 3 /kg) Subscripts a Environment ashp Air source heat pump com Compressor dis Compressor discharge i Inlet, inner face l Liquid o Outlet, outer face ow-hp Output water of the air source heat pump ow-hs Output water of the heating system rw-hp Return water of the air source heat pump rw-hs Return water of the heating system segp Solar equivalent generation power sh Superheated zone sc Subcooled zone suc Compressor suction tp Two-phase zone v Vapor, electronic expansion valve w Water 1 Introduction Conventional building heating schemes mainly include boiler heating, heat pump heating and solar energy heating. With the deterioration of energy shortage and environment pollution, solar energy can be used as the substitution of conventional energy for building heating and plays an important role in reducing the dependence on conventional energy resources. Solar energy has the characteristics of intermittence and low density, which largely restrict the application of solar heating. Heat pump can transit the heat from low-grade heat sources to high-grade heat sources. So the solar assisted heat pump heating system, which combines heat pump technology with solar heating technology, can solve the intermittent problem of solar energy. Since the 1950s, many scholars have studied on the various solar assisted heat pumps [1–4] . Kaygusuz theoretically and experimentally studied a solar assisted heat pump with a phase transition heat conservation facility [5] . Han et al. built up a mathematical model of solar assisted ground source heat pump with a potential energy storage device and obtained the effect of the heat conservation device on the performance of solar assisted ground source heat pump [6] . Ozgener and Hepbasli set up a greenhouse heating system based on the solar assisted ground source heat pump and studied the operating characteristics of the system. Meanwhile, he got the effect of the arrangement of pipes on the performance of solar assisted ground source heat pump heating system [7,8] . Kuang et al. built up an experimental testbed of solar assisted heat pump heating system and carried out experiments. The effects of environmental parameters and structural parameters on the performance of the system were obtained [9] . Meanwhile, he developed a direct-expansion solar assisted heat pump system for domestic use, which could offer space heating in winter, air conditioning in summer and hot water during the whole year, and studied the long-term performance of the system [10] . Xu et al. proposed a new type of solar-air source heat pump water heater, whose flat-plate heat collector/evaporator with spiral-finned tubes can obtain energy from both solar irradiation and ambient air for hot water heating, and simulated the operating performance of the system [11] . Li et al. built up a direct-expansion solar assisted heat pump water heater, studied the system performance and proposed some suggestions for the system optimization [12] . Dikici and Akbulut set up an experimental device of solar assisted heat pump with a flat plate collector and studied the effect of energy loss on the performance of the system [13] . Badescu designed a model of solar assisted air source heat pump with a sensible heat storage device and mainly studied the effect of heat storage device on the performance of the system [14] . Georgiev carried out an experimental study with solar collectors as the heat source and got the effect of the collector parameters on the performance of the system [15] . Wang and Qi studied the geothermal heat characteristics of solar assisted ground source heat pump, and the results showed that the geothermal heat characteristics depended strongly on the intensity of solar radiation and the matching degree between the solar collector area and the capacity of water tank [16] . The solar assisted heat pumps mentioned above can be classified into direct-expansion solar assisted heat pump [9–12] and indirect-expansion solar assisted heat pump [4–8,13–16] . According to the different connection styles between the solar collector and the heat pump, indirect-expansion solar assisted heat pump can also be divided into three styles: series style [13,15] , parallel style [4] and dual heat source style [5–8,14,16] . Direct-expansion solar assisted heat pump, series style and dual heat source style utilize solar energy as their low-grade heat source. The utilization of solar energy needs to consume electrical energy by the heat pump to improve the energy grade. For parallel style, the utilization of solar energy is independent of the heat pump and the performance of heat pump cannot be improved by solar energy assistance. Therefore, we proposed a new solar assisted heat pump system for building heating. In this system, solar energy can be used to improve the performance of heat pump. The matching degree between the solar collector area and the system integral capacity is one of the most important factors which affect the performance of the system. So by setting up the model of this heating system, the paper further studied the effects of the solar collector area on the operating characteristic and the performance of the heating system, which could be helpful for the design optimization and cost-effective evaluation of the solar assisted heat pump system. 2 New solar assisted air source heat pump system for building heating Because of the flexibility, convenience and low investment of air source heat pump, we combined solar energy with air source heat pump to form a new solar assisted air source heat pump system (SAASHPS) for building heating. The whole heating system consisted of two parts: air source heat pump unit and solar collection unit. The schematic diagram of the new SAASHPS was shown in Fig. 1 . The output water temperature for heat supply of the system ( T ow-hs ) was 45 °C and the return water temperature of the system ( T rw-hs ) was 40 °C. In the loop of hot water, the condenser of air source heat pump was in series with solar collector, and the air source heat pump was in front of the solar collector. The system worked by three modes: (1) single air source heat pump mode, which is operated when the sun is not available or solar radiant intensity is low, then the available heat of radiation is smaller than or equal to the heat dissipation of solar collector with the water temperature of 45 °C (the available heat of radiation refers to the heat which can be absorbed by water in the solar collector); (2) combination mode, which is operated when the available heat of radiation is more than the heat dissipation of solar collector with the water temperature of 45 °C, and the available heat of radiation is not adequate to increase the water temperature by 5 °C; (3) single solar collector mode, which is operated when the solar collector can provide adequate heat to increase the water temperature by 5 °C. The specific process of each mode is as follow: 1) Single air source heat pump mode: the solenoid valves 2 and 4 are closed, and the solenoid valves 1 and 3 are opened; the return water passes through the solenoid valve 1 and enters the shell and tube heat exchanger, then the water cools down refrigerant and absorbs the heat, whose temperature increases to 45 °C; after it goes out of the shell and tube heat exchanger, it gets through the solenoid valve 3 to reach the hot water container, and flows out of system after being pressurized by the pump to the air handling unit. 2) Combination mode: the solenoid valves 2 and 3 are closed, and the solenoid valves 1 and 4 are opened; the return water passes through the solenoid valve 1 and enters the shell and tube heat exchanger, then the water cools down refrigerant and absorbs the heat, whose temperature increases; after it goes out of the shell and tube heat exchanger, it gets through the solenoid valve 4 to reach the solar collector and absorbs the heat; when the temperature increases to 45 °C, it enters the hot water container and flows out the system after being pressurized by the pump to the air handling unit. 3) Solar collector direct heat supply mode: the solenoid valves 1 and 3 are closed, and the solenoid valves 2 and 4 are opened; the return water passes through the solenoid valves 2 and 4, then the water directly enters the solar collector to absorb the heat; when the temperature rises to 45 °C, it enters the hot water container and flows out the system after being pressurized by the pump to the air handling unit. The compressor in the air source heat pump is a variable capacity compressor. Under the combination mode, the compressor runs at variable capacity to maintain the T ow-hs at 45 °C. The higher the solar radiant intensity is, the lower capacity the compressor runs with, and vice versa. The testing results by solar radiant intensity instrument are used as the reference for the switching between the modes. Under the combination mode, the temperature of output water from the condenser of the air source heat pump ( T ow-hp ) declines with the increase of solar radiant intensity. In this way, the condensing pressure and condensing temperature of the heat pump system will be decreased, which make the improvement of heat pump performance. 3 Establishment of the heating system model 3.1 Model of the air source heat pump unit As we aim to obtain the operational characteristic and the performance rules of the ASASHPHS at steady state, a steady-state model of the heating system was built up. The SAASHPS was designed without heat storage, and its rated heating capacity was 10 kW. So the heating capacity of the air source heat pump was 10 kW under the rated heating condition. The compressor of the heat pump was a frequency conversion compressor, and the refrigerant of heat pump was R22. The finned tube heat exchanger worked as evaporator, the shell and tube heat exchanger worked as condenser, and the throttle was an electronic expansion valve. In the air source heat pump model, the following assumptions about heat exchangers were made: (1) it is thermal isolation between the shell and tube heat exchanger and the environment; (2) the pressure inside heat exchanger and the temperature of refrigerant inside two-phase zone are constant; (3) only pure refrigerant is considered; (4) axial heat transfer is neglected; and (5) the refrigerant from each channel of the heat exchanger is fully mixed. 3.2 Model of the shell and tube heat exchanger 1) Heat transfer coefficients of refrigerant in single-phase zone The heat transfer coefficient of superheated zone was α sh , and the heat transfer coefficient of subcooled zone was α sc , which can both be calculated by the standard Dittus–Boelter formula: (1) α s h = 0.023 λ v d i Re v 0.8 Pr v 0.3 (2) α s c = 0.023 λ l d i Re l 0.8 Pr l 0.3 2) Heat transfer coefficient of refrigerant in two-phase zone Local heat transfer convection coefficient in the two-phase zone, a tp ( x ), was evaluated by Dobson formula: (3) α t p ( x ) = 0.023 λ l d i R e l 0.8 P r l 0.4 1 + 2.22 X t t 0.889 X t is Martinelli value, and the equation was as follow: (4) X t t = 0.551 p P c r 1 − x x 0.9 α tp can be got by (5) α t p = ∫ x i x o d x ∫ x i x o ( 1 / α t p ( x ) ) d x 3) Heat transfer coefficient of external water The shell side of the shell and tube heat exchanger was equipped with baffle plates. Heat exchange happened by the water flowing vertically and horizontally outside the pipes. The heat transfer coefficient of external water, α o , was (6) α o = 0.25 λ d o R e 0.6 P r 0.33 3.3 Model of the finned tube exchanger 1) Heat transfer coefficient of refrigerant in the superheated zone In the superheated zone, the refrigerant was in the state of superheated vapor and moved in the form of turbulent flow inside the pipes. Heat transfer coefficient of refrigerant, α sh , was evaluated by Petukhov–Popov equation: (7) α s h = λ v d i ⋅ ( f / 8 ) R e v P r v 1.07 + 12.7 ( f / 8 ) 0.5 ( Pr v 2 / 3 − 1 ) f is impact factor, and the equation was as follow: (8) f = ( 1.82 lg R e v − 1.64 ) − 2 2) Average heat transfer coefficient of refrigerant in two-phase zone The local heat transfer coefficient in the two-phase zone, a tp ( x ), which was related to the refrigerant dryness fraction, x , may be evaluated by (9) α t p ( x ) = 3.0 X t t 2 / 3 α l (10) X t t = μ l μ v 0.1 1 − x x 0.9 ρ v ρ l 0.5 The heat transfer coefficient of refrigerant in pure liquid phase, α l , can be calculated by the standard Dittus–Boelter formula, (11) α l = 0.023 R e l 0.8 P r l 0.4 λ l d i The average heat transfer coefficient of the two-phase zone, α tp , may be evaluated by the following formula: (12) α t p = ∫ x i x o d x ∫ x i x o ( 1 / α t p ( x ) ) d x 3) Heat transfer coefficient of external air The fin type of the finned tube exchanger was a straight one-piece fin. The heat transfer coefficient of external air, α f , can be obtained [17] : (13) N u = 0.982 R e 0.424 s d b − 0.0887 N ⋅ s 2 d b − 0.1590 (14) α f = λ d b N u 3.4 Model of the compressor The suction pressure and discharge pressure of compressor were assumed as the vaporizing pressure and condensing pressure of the system respectively. The flow of refrigerant in compressor, G com , depended on the theoretical capacity of compressor V th , the capacity efficiency of compressor η com and the specific volume of refrigerant at the entrance of compressor ν com : (15) G c o m = η c o m V t h v c o m (16) V t h = V c o m ⋅ n 60 η com was provided by the compressor manufacturer as a constant. The input power of compressor, W in , can be calculated by the following formulas: (17) W o u t = ( h d i s − h s u c ) ⋅ G c o m (18) W i n = W o u t η k 3.5 Model of the electronic expansion valve Under steady working condition, the flow of the electronic expansion valve, m v , was (19) m v = C v A v ρ v , i ( p v , i − p v , o ) 3.6 Model of the solar collection unit The collector of the solar collection unit was the heat pipe evacuated tube solar thermal collector of model BTZ22. The instantaneous efficiency of collector, η , can be obtained (Neglecting the loss in the heat exchange between the pipes and the environment): (20) η = 0.682 − 2.32 T i − T a I 3.7 Validation of the mathematical model To validate the mathematical model of the system, an air source heat pump unit was developed. The compressor of the unit was a frequency conversion rotary compressor of model THS20MC6-Y, and the compressor frequency can be adjusted from 15 Hz to 110 Hz by a general frequency converter, the main parameters of the compressor are shown in Table 1 . The shell and tube heat exchanger of the unit was model RER-20, whose main parameters are shown in Table 2 in detail. The shape of the finned tube heat exchanger of the unit was U-shaped with flat fin, and the row number was 2. The main parameters of the finned tube heat exchanger are shown in Table 3 . The layout of measuring spots is shown in Fig. 1 . The temperature was measured by the platinum-resistance temperature sensors with an accuracy of 0.1 °C and the measuring range from −50 to 150 °C. The measuring range of the pressure sensors was from 0 to 2.5 MPa with an accuracy of 0.25%. The water flow rate was measured by the turbine flow meter of model LWY-15C with an accuracy of 1% and the measuring range from 0 to 6 m 3 /h. The power of the compressor was measured by the power meter of model WT230 with an accuracy of 0.1%. All the data for every spot were collected in real-time by a data collector Agilent of model 34970a. In order to get a steady experimental environment for the unit, the unit was placed in an environmental chamber as shown in Fig. 2 , which was used to obtain constant temperature and relative humidity of the air. The temperature and relative humidity of the air in the chamber could be adjusted; the temperature could be manually regulated from −5 °C to 43 °C with an accuracy of 0.2 °C; the relative humidity could be regulated from 5% to 95% with an accuracy of 2%. The simulations and tests of the unit have been done under the environmental temperature of 7 °C, the return water temperature of the air source heat pump ( T rw-hp ) of 40 °C and a variety of the output water temperature of the air source heat pump ( T ow-hp ). Fig. 3 shows the effect of T ow-hp on the relative COP of the air source heat pump, which took the COP of air source heat pump under the rated condition as a benchmark. In this figure, it can be seen that the variation trend of the simulation results was in good agreement with the experimental results, and the maximum deviation was less than 8%. 4 Results and discussion 4.1 Performance of the heating system in a typical sunny day in the heating season One of the main factors which affect the performance of the solar assisted air source heat pump is the matching degree between the solar collector area and the system. So in order to obtain the effect of the solar collector area on the performance of the system, the study was carried out on the new SAASHPS whose heating capacity was 10 kW under the rated condition. The paper studied five different solar collector areas of 0 m 2 , 10 m 2 , 20 m 2 , 30 m 2 and 40 m 2 . It was equal to the sole working of air source heat pump when the solar collector area was 0 m 2 . In the simulation process, we did not consider the change in the thermal loads of buildings and keep the heating capacity of 10 kW constant. Taking the weather condition of Nanjing as an example, the paper chose the weather parameters of typical 24 h (a typical sunny day) in the heating season in winter as the operating condition of the heating system. The weather parameters of a typical sunny day are shown in Fig. 4 . In this day, the max environment temperature was 9.5 °C, the min value was −1.2 °C, and the max intensity of solar radiation was 537 W/m 2 . The simulations were done based on the model of the heating system. In the typical sunny day, the real-time curves of T ow-hp are shown in Fig. 5 . As can be seen from the figure, the heating system operated under the single air source heat pump mode and T ow-hp was equal to T ow-hs before 8:00 when there was no sunshine; with the increase of solar radiation intensity, the heating system ran at the combination mode, and T ow-hs was kept at 45 °C, while T ow-hp declined with the increase of solar radiation intensity; T ow-hp climbed up when the radiation intensity decreased. For the solar collector area of 20 m 2 , the lowest T ow-hp reached 42.15 °C. With the increase of solar collector area, the lowest T ow-hp became lower and lower; when the solar collector area was 40 m 2 , the system ran at the single solar collector mode, the air source heat pump shut down, and T ow-hp was 40 °C between 12:00 and 14:00. Fig. 6 shows the power curves of the air source heat pump system ( W ashp ) in the typical sunny day. In this figure, the heating system ran at the single air source heat pump mode and W ashp did not change with the increase of solar collector area when the sun was not available; with the increase of solar radiation intensity, the heating system operated under the combination mode and W ashp decreased; when the solar collector area was 0 m 2 , W ashp also changed with the increase of solar radiation intensity because of the change of environment temperature; W ashp became less and less with the increase of solar collector area when solar radiation was available; the heating system ran at the single solar collector mode from 12:00 to 14:00 when the solar collector area was 40 m 2 , at the same time, the air source heat pump shut down. In the typical sunny day, the heating capacity curves of the air source heat pump system ( q ashp ) are shown in Fig. 7 . From the figure, the heating system ran at the single air source heat pump mode and q ashp was constant for different solar collector areas when there was not sunshine. The heating system operated under the combination mode and q ashp decreased with the increase of solar radiation intensity. q ashp was gradually declining with the increasing of solar collector area; when the solar collector area was 40 m 2 , q ashp was zero and the air source heat pump shut down from 12:00 to 14:00, at the same time, the heating system worked at the single solar collector mode. Fig. 8 shows the COP curves of the air source heat pump in the typical sunny day when the sun was available, and the heating system ran at the combination mode. In this figure, the COP gradually declined under the combination mode with the increase of solar radiation intensity; when the solar collector area was 20 m 2 , the solar radiation intensity reached the max value at 13:00, compared with the solar collector area of 0 m 2 , the COP of the air source heat pump was improved by 17%; the COP of air source heat pump increased by 8.3% on average under the combination mode from 8:00 to 17:00. Thus it can be seen that it is able to greatly ameliorate the working conditions, enhance the performance and reduce the energy consumption while the new heating system can make full use of solar energy. With the increase of solar collector area, the COP of the air source heat pump increased at the same time. When the solar collector area was 40 m 2 , the heating system ran at the single solar collector mode, and the air source heat pump shut down between 12:00 and 14:00, and thereby the curve was discontinuous. Table 4 shows the performance of the new SAASHPS at different solar collector areas in the typical sunny day. As can be seen from this table, the total solar energy and available solar energy increased in a similarly linear way with the increase of solar collector area. Meanwhile, the capacity of the heating system did not change with the increase of solar collector area, but the heating capacity of the air source heat pump unit decreased. The heating system realized the direct application of solar energy. The application of solar energy reduced the heating capacity of air source heat pump unit, ameliorated the working conditions and improved the efficiency of the air source heat pump unit. The amount of energy saving and energy-saving rate were also improved with the increase of solar collector area. Compared with the case when the solar collector area was 0 m 2 , the saved power of the air source heat pump unit was 15.69 kW h and the energy-saving rate was 24% when the solar collector area was 40 m 2 . Compared with the air source heat pump sole heating system, the new system realized the solar generation. The solar equivalent generation power efficiency ( e segp ) was got by comparing the saved power consumption of air source heat pump to the total solar energy. Fig. 9 indicates the curve of e segp with the solar collector area in the typical sunny day. In this figure, it can be seen that e segp of the new heating system could reach 11.8% when the solar collector area was 40 m 2 . With the increase of solar collector area, e segp gradually decreased. Though initial investment was increased for the additional sets of solar collection unit compared with the air source heat pump sole heating system, the new SAASHPS was able to generate electricity indirectly by solar energy (the saved power consumption by solar energy was equivalent to the indirect solar generation). The increased investment of system was far less than that of the solar photovoltaic generation equipment. Therefore, in this situation, the thermal utilization of solar energy was superior to the solar photovoltaic generation. 4.2 Performance of the heating system in the heating season The heating season of Nanjing lasts from December to February. The paper further studied the characteristics of the new heating system in the whole heating season. According to the model of the heating system, we made a simulation in the heating season calculated by 92 days with hour as the basic unit and 2208 h in total. In the simulation process, we selected a building with the thermal load of 10 kW under the rated condition, and the thermal load was affected by the weather. In addition, the heating capacity of the heating system was also affected by the weather. By calculating the performance of the heating system in the whole heating season hour by hour, we made the integral performance of the new SAASHPS in Table 5 . The effect of the solar collector area on e segp in the heating season is shown in Fig. 10 . As can be seen from Table 5 , in the whole heating season, the trends of all the performance parameters of the heating system were the same with those in the typical sunny day. The total solar energy and the available solar energy in the heating system increased with the increase of solar collector area in the similarly linear way. The heating capacity of the air source heat pump unit decreased with the increase of solar collector area. The amount of energy saving and energy-saving rate of the system increased with the increase of solar collector area, but the absolute values of both changes were different. In the whole heating season, compared with the case when the solar collector area was 0 m 2 , the heating system could save 453.43 kW h with the energy-saving rate of 9.90% when the solar collector area was 40 m 2 . Meanwhile, in Fig. 10 , it can be seen that e segp was 5%. In the whole heating season, the energy-saving rate and e segp were lower than those in the typical sunny day. The main reasons were that there were changes caused by rainy and cloudy days and the variation of solar radiation intensity in the whole heating season. 5 Conclusions This paper proposed a new solar assisted air source heat pump system for building heating, and studied the performance of the system. Meanwhile, aiming to resolve the problem about the optimization between the solar collector area and the system, the effect of the solar collector area on the performance of this new heating system was studied. The following conclusions can be drawn: 1) Taking the weather conditions of Nanjing as an example, the COP of the air source heat pump increased with the increase of solar radiation intensity in the typical sunny day in the heating season. At the same time, the COP of the air source heat pump increased with the increase of solar collector area. When the solar collector area was 20 m 2 , the solar radiation intensity reached the max value at 13:00, compared with the case when the solar collector area was 0 m 2 , the COP of the air source heat pump unit could be improved by 17%; the COP of air source heat pump is increased by 8.3% on average from 8:00 to 17:00. 2) Meanwhile, in the typical sunny day in the heating season, when the solar collector area was 40 m 2 , compared with the solar collector area of 0 m 2 , the heating system with the heating capacity of 10 kW could save 15.69 kW h with the energy-saving rate of 24% and e segp of 11.8%; the COP of the air source heat pump was improved by 11.22% on average from 8:00 to 17:00. 3) In the whole heating season in winter, compared with the solar collector area of 0 m 2 , the heating system could save 453.43 kW h with the energy-saving rate of 9.9% and e segp of 5% when the solar collector area was 40 m 2 . Acknowledgements This research was supported by the grants from the Natural Science Foundation of China (no. 50676018 ), the 11th Five Year Science and Technology Support Key Project of China (no. 2008BAJ12B02 ), the Natural Science Foundation of Jiangsu Province (key project) (no. BK2010029 ) and the Science Foundation of Southeast University ( KJ2010397 ). References [1] P. Sporn E.R. Ambrose The heat pump and solar energy Proceedings of the World Symposium on, Applied Solar Energy Phoenix, Ariz 1955 [2] T.L. Freeman J.W. Mitchell T.E. Audit Performance of combined solar-heat pump systems Solar Energy 22 1979 125 135 [3] A.R. Day T.G. Karayiannis Solar-assisted heat pump research and development Building Service Engineering Research & Technology 15 1994 71 80 [4] P.J. Martínez A. Velázquez A. Viedma Performance analysis of a solar energy driven heating system Energy and Buildings 37 2005 1028 1034 [5] K. Kaygusuz Experimental and theoretical investigation of a solar heating system with heat pump Renewable Energy 21 2000 79 102 [6] Z. Han M. Zheng F. Kong F. Wang Z. Li T. Bai Numerical simulation of solar assisted ground-source heat pump heating system with latent heat energy storage in severely cold area Applied Thermal Engineering 28 2008 1427 1436 [7] O. Ozgener A. Hepbasli Experimental performance analysis of a solar assisted ground-source heat pump greenhouse heating system Energy and Buildings 37 2005 101 110 [8] O. Ozgener A. Hepbasli Modeling and performance evaluation of ground source (geothermal) heat pump systems Energy and Buildings 39 2007 66 75 [9] Y. Kuang R. Wang Y. Xu Studies on the performance of a direct-expansion solar-assisted heat pump water heating system Journal of Engineering Thermophysics 25 2004 737 740 (in Chinese) [10] Y.H. Kuang R.Z. Wang Performance of a multi-functional direct-expansion solar assisted heat pump system Solar Energy 80 2006 795 803 [11] G.Y. Xu X.S. Zhang S.M. Deng A simulation study on the operating performance of a solar-air source heat pump water heater Applied Thermal Engineering 26 2006 1257 1265 [12] Y.W. Li R.Z. Wang J.Y. Wu Y.X. Xu Experimental performance analysis and optimization of a direct expansion solar-assisted heat pump water heater Energy 32 2007 1361 1374 [13] A. Dikici A. Akbulut Performance characteristics and energy–exergy analysis of solar-assisted heat pump system Building and Environment 43 2008 1961 1972 [14] V. Badescu Model of a thermal energy storage device integrated into a solar assisted heat pump system for space heating Energy Conversion Manage 44 2003 1589 1604 [15] A. Georgiev Testing solar collectors as an energy source for a heat pump Renewable Energy 33 2008 832 838 [16] H. Wang C. Qi Performance study of underground thermal storage in a solar-ground coupled heat pump system for residential buildings Energy and Buildings 40 2008 1278 1286 [17] W. Li W. Tao H. Kang Analysis of the elasto-plastic stress in screwed-pipe joint Chinese Journal of Mechanical Engineering 33 1997 81 86 (in Chinese)