# Analysis of a Hybrid PV/Thermal Solar-Assisted Heat Pump System for Sports Center Water Heating Application

Authors: Y. Bai; T. T. Chow; C. Ménézo; P. Dupeyrat
Journal: International Journal of Photoenergy
 (2012)

---

## Abstract

The application of solar energy provides an alternative way to replace the primary source of energy, especially for large-scale installations. Heat pump technology is also an effective means to reduce the consumption of fossil fuels. This paper presents a practical case study of combined hybrid PV/T solar assisted heat pump (SAHP) system for sports center hot water production. The initial design procedure was first presented. The entire system was then modeled with the TRNSYS 16 computation environment and the energy performance was evaluated based on year round simulation results. The results show that the system COP can reach 4.1 under the subtropical climate of Hong Kong, and as compared to the conventional heating system, a high fractional factor of energy saving at 67% can be obtained. The energy performances of the same system under different climatic conditions, that include three other cities in France, were analyzed and compared. Economic implications were also considered in this study.


---

## Body

## 1. Introduction

There are various applications of solar energy for medium-grade heating, including water heating, space heating, distillation, and the like. The utilization of solar energy for hot water production is one most popular application. However, large rooftop or other space in buildings is needed for installing enough solar collector arrays to satisfy the full services requirements in the large-scale and/or multiple-demand cases. From this point of view, the combined solar-assisted heat-pump system appears a suitable alternative, which not only saves building space but also reduces the reliance on utilities electricity supply. The combined system is then able to work more efficiently.Unlike the conventional solar thermal collectors, the hybrid photovoltaic/thermal (PV/T) collectors make possible higher energy outputs per unit surface area because the absorbed solar radiation is converted into electricity and usable heat simultaneously. The studies on this topic were initiated in the 1970s and since then many innovative systems have been introduced. Several recent technical reviews [1–3] cover the latest developments of PV/T collector designs and their performances in terms of electrical/thermal outputs as well as application potentials. The increased attention and research outputs in this area indicate that hybrid PV/T technologies present many attractive features for wide applications and industrialization opportunities for large scale production. Specifically in the research area of hot water co-generation, a series of experimental and numerical studies have been undertaken on the stand-alone and building-integrated photovoltaic/water-heating (PV/W) systems for warm climate application [4–6]. The year round thermal and cell conversion efficiencies were found 37.5% and 9.4% respectively for the case of BiPV/W application in Hong Kong. The overall heat transmission through the PVW wall was reduced to 38% of the normal building façade. Santbergen et al. [7] took a detailed analysis of the energy yield of solar domestic hot water systems with covered sheet-and-tube PV/T collectors. A detailed quantitative analysis of all loss mechanisms inherent to the PV/T collectors was performed on top of those related to PV modules and conventional thermal collectors. The annual electrical efficiencies of the PV/T systems they investigated were found lower than the plain PV systems (up to 14% relatively) and the annual thermal efficiencies were also lower than the conventional thermal collector systems (up to 19% relatively). With the aim to improve the overall system performance, many parametric studies have been done [8–11]. Charalambous et al. [8] carried out an optimization study based on the “low-flow” concepts, of which the advantages include improved system performance, smaller pump, reduced size and thickness of tubing and insulation, less construction work and time for the optimum absorber configuration, and thus also cost saving. The optimized flow rate for the header and riser in serpentine PV/T collectors is determined, respectively, by using the EES code. A similar study was taken by Cristofari et al. [9] on copolymer PV/T collectors.To improve the global system efficiency in large-scale installations of hybrid PV/T technology, the combined use with other efficient energy systems (such as heat pump systems) can be promising.A number of research investigations have been conducted in the design, modeling, and testing of solar-assisted heat pump systems. Ozgener and Hepbasli [12, 13] reviewed the reported works on energy and exergy analysis of SAHP systems. Bridgeman and Harrison [14] conducted a preliminary experimental evaluation of indirect SAHP system for domestic hot water applications; tests were conducted with a range of evaporator supply temperatures and the results indicated that the COP can span from 2.8 to 3.3, depending on the evaporator and condenser temperatures. Dikici et al. [15] performed an energy-exergy analysis of solar-assisted heat pump system for space heating in a test room of 60 m2. The system COP was determined as 3.08 while the exergy loss of the solar collector was found to be 1.92 kW. The authors concluded that the COP increased when the exergy loss of evaporator decreased. Through system simulation  Li and Yang [16] studied a parallel SAHP system for supplying hot water to a hypothetical residential building; the author concluded that the solar collector area affects the optimum flow rate significantly. Wang et al. [17] developed a novel indirect-expansion solar-assisted multifunctional heat pump (IX-SAMHP). The proposed system not only works in operation modes included in the two household appliances, but also operates in four new energy-saving operation modes for space cooling, space heating, and water heating. The experimental results indicated that the IX-SAMHP on cloudy days can produce hot water with considerable reduction in electricity consumption than a solar water heater and in cold winter can operate in much higher COP (3.5–4.2) than a domestic heat pump. Sterling and Collins [18] carried out a feasibility analysis of an indirect SAHP system for domestic water heating as compared to: (i) a traditional solar domestic hot water system and (ii) an electric domestic hot water system. It was found that the electrical consumption and operating cost were most favorable with the indirect SAHP system.From our literature review, we found that the experimental and theoretical analysis of hybrid SAHP system for large-scale water heating application is very limited. In this study, an indirect hybrid PV/T solar-assisted heat pump system has been investigated from the initial design phase to the detailed analysis of the annual performance through numerical simulation. The cost investment on the proposed system was also estimated.


## 2. System Description

A hypothetical sports center was taken as a reference case. A perspective view of the building is shown in Figure1. The floor area of the sports center is 3200 m2, corresponding to 80 m (L) by 40 m (W). The provisions include an indoor swimming pool, one general-games sport hall, several general minigames rooms, one gym, and the relevant services spaces like changing rooms, canteen, lobby, and office. The central plant room is located at the ground level. In this case study, the designed energy system was used for hot water production to cater for bathing and hand washing in the changing rooms.Figure 1
Prospective view of the sports center model.A simplified schematic diagram of the proposed PVT-SAHP system is presented in Figure2. The solar heat source is for water preheating. An indirect solar water heating mode was chosen in this case. The solar preheating system mainly consists of an array of hybrid PV/T collectors, a plate heat exchanger, and solar storage tanks. On the hot side of the heat exchanger, the fluid in the thermal collector circulates according to the temperature differential between the collector outlet T1 and the bottom of the storage tank T7, which was monitored by an on-off differential controller R1. Upper and lower dead bands were set at 10°C and 2°C, respectively. On the cold side of the exchanger, cold water is taken from the bottom of the tank and delivered to the heat exchanger connected to the collector where it gains solar energy and returns to the tank at a higher temperature. The water circulation is controlled by the on-off differential controller R2 according to the tem differential between the entrance of the heat exchanger on the hot side T2 and the bottom of the tank T7. The dead bands values are the same as inR1. The feed water stream is divided into two parts. On part enters from the bottom of the solar tank, and another part is mixed with the water exiting from the heat pump system. The proportion is controlled by a temperature-controlled liquid flow diverter according to the designed load water temperature. After preheating by the solar collector system and the small back-up electric heater in the tank, water then passes through the heat pump for reheating as the temperature is still not reaching the desired delivery water temperature. The internal heater in the solar tank is activated only during the winter season. The external auxiliary heater is functioned to keep the desired pool water temperature, of which the set point temperature should be higher than the required delivery temperature in order to cover the heat losses in the pipe distribution system.Figure 2
Schematic overview of the indirect PV/T solar-assisted heat pump system.


## 3. Design of Solar Water Preheating System

### 3.1.  Calculation of the Daily Heating Demand of the Sports Center

According to the practical design guides of sports center hot water heating, the water load can be determined according to: (i) the daily average hot water consumption per person or (ii) the hourly water consumption per shower faucet and lavabo.In our study, the hot water load is calculated by using the daily water loadQh, which depends on the number of showers faucets and lavabos in the provision, hence
(1)m˙dis=∑qhn0b,
where qh  is the hourly water consumption of shower faucets and lavabos, in kg/h; n0 is the number of the shower faucets and lavabos, b is the simultaneous using factor of the shower facets and lavabos within one hour, which is usually taken as 100% for the sports center cases.The hot water demand is then calculated by:(2)mh˙=krm˙dis,
where kr is the hot water mixing factor. This is determined according to the designed hot water temperature at the outlet of the heat pump system and the delivery temperature, as shown (3) that follows:(3)kr=Tdis-TfwTh-Tfw,
where Tdis is the delivery water temperature after mixing, in °C; Th is the hot water temperature at the outlet of the tank, in °C; Tfw is the feed water temperature, in °C.The required heating load is then calculated on the basis of:(4)Qh=mh˙Cp(Th-Tfw).


### 3.2. Estimation of the Required Solar Collector Area

The required thermal collector area is dependent on several factors including the daily water heating load of the sports center, the collector characteristics, and the climatic conditions. The area of solar collector is determined by the following relation in the direct heating mode:(5)Acd=24Qf[JTηcd(1-ηL)],
where Ac is the total collector area when the indirect heating mode is used, in m2; JT is the average daily solar radiation level, in kJ/m2·d; ηL the heat loss coefficient of the hydraulic piping system, normally between 0.2 and 0.3 for the well-insulated case; ηcd is the collector efficiency determined by practical test, in general this value varies between 0.25 and 0.5.For the indirect heating mode, as in this study, the required collector area is determined on considering the heat exchanger characteristics, as shown in (6):
(6)Ac=Acd(1+δAcdKhxAhx),
where Acd is the total collector area when the direct heating mode is used, in m2; Ahx is the heat exchanger surface area, in m2; f is the general solar fraction according to the different climatic zone; δ is the collector heat loss coefficient, in W/(m2·K); Khx is the heat transfer coefficient of the heat exchanger, in W/(m2·K).In this case study, arbitrarilyηL=0.2 and   ηcd=0.4.


### 3.3. Modeling of the Heat Pump System

In the heat pump unit, the heating capacity of the condenser can be calculated by(7)Qcond=mrf(hout-hin),
where mrf is mass flow rate of the refrigerant (R410A), kg/s; hout and hin are the specific enthalpy of R410A at the inlet and outlet of the condenser, kJ/kg.The heat transfer to water is determined by the following relation:(8)Qdhw=mdhwCp(Tdhw_in-Tdhw_out),
where mdhw is mass flow rate of hot water, kg/s; Tdhw_in and Tdhw_out are the inlet and outlet hot water temperatures, °C.The power consumption of the heat pump unit is the total energy consumed by the compressor, water pumps, and the controller, in kWh.(9)Qhp=∑τ=0t(Pcomp+Pcontroller+Ppump),
wherePcomp,Pcontroller, and Ppump represent the power of the compressor, the controller, and the water pumps, respectively, in kW.The coefficient of performance of the overall heating system can be defined as(10)COPsys=QdhwQhp.


## 3.1.  Calculation of the Daily Heating Demand of the Sports Center

According to the practical design guides of sports center hot water heating, the water load can be determined according to: (i) the daily average hot water consumption per person or (ii) the hourly water consumption per shower faucet and lavabo.In our study, the hot water load is calculated by using the daily water loadQh, which depends on the number of showers faucets and lavabos in the provision, hence
(1)m˙dis=∑qhn0b,
where qh  is the hourly water consumption of shower faucets and lavabos, in kg/h; n0 is the number of the shower faucets and lavabos, b is the simultaneous using factor of the shower facets and lavabos within one hour, which is usually taken as 100% for the sports center cases.The hot water demand is then calculated by:(2)mh˙=krm˙dis,
where kr is the hot water mixing factor. This is determined according to the designed hot water temperature at the outlet of the heat pump system and the delivery temperature, as shown (3) that follows:(3)kr=Tdis-TfwTh-Tfw,
where Tdis is the delivery water temperature after mixing, in °C; Th is the hot water temperature at the outlet of the tank, in °C; Tfw is the feed water temperature, in °C.The required heating load is then calculated on the basis of:(4)Qh=mh˙Cp(Th-Tfw).


## 3.2. Estimation of the Required Solar Collector Area

The required thermal collector area is dependent on several factors including the daily water heating load of the sports center, the collector characteristics, and the climatic conditions. The area of solar collector is determined by the following relation in the direct heating mode:(5)Acd=24Qf[JTηcd(1-ηL)],
where Ac is the total collector area when the indirect heating mode is used, in m2; JT is the average daily solar radiation level, in kJ/m2·d; ηL the heat loss coefficient of the hydraulic piping system, normally between 0.2 and 0.3 for the well-insulated case; ηcd is the collector efficiency determined by practical test, in general this value varies between 0.25 and 0.5.For the indirect heating mode, as in this study, the required collector area is determined on considering the heat exchanger characteristics, as shown in (6):
(6)Ac=Acd(1+δAcdKhxAhx),
where Acd is the total collector area when the direct heating mode is used, in m2; Ahx is the heat exchanger surface area, in m2; f is the general solar fraction according to the different climatic zone; δ is the collector heat loss coefficient, in W/(m2·K); Khx is the heat transfer coefficient of the heat exchanger, in W/(m2·K).In this case study, arbitrarilyηL=0.2 and   ηcd=0.4.


## 3.3. Modeling of the Heat Pump System

In the heat pump unit, the heating capacity of the condenser can be calculated by(7)Qcond=mrf(hout-hin),
where mrf is mass flow rate of the refrigerant (R410A), kg/s; hout and hin are the specific enthalpy of R410A at the inlet and outlet of the condenser, kJ/kg.The heat transfer to water is determined by the following relation:(8)Qdhw=mdhwCp(Tdhw_in-Tdhw_out),
where mdhw is mass flow rate of hot water, kg/s; Tdhw_in and Tdhw_out are the inlet and outlet hot water temperatures, °C.The power consumption of the heat pump unit is the total energy consumed by the compressor, water pumps, and the controller, in kWh.(9)Qhp=∑τ=0t(Pcomp+Pcontroller+Ppump),
wherePcomp,Pcontroller, and Ppump represent the power of the compressor, the controller, and the water pumps, respectively, in kW.The coefficient of performance of the overall heating system can be defined as(10)COPsys=QdhwQhp.


## 4. System Simulation

The entire system simulation was carried out based on the initial design parameters under the TRNSYS transient simulation environment [19]. This computer program is popularly used by the international scientific community in the energy and thermal engineering field. It is based on interconnecting elementary modules called TYPE, which are either components of the studied system such as storage tanks or particular functions such as the weather data reader which allows connecting the selected area weather database to the system.The numerical model of each component used in this study was either available in the existing library or developed as a new component. The hybrid PV/T collector model was developed based on a prototype from the collaborative research of the Fraunhofer Institute for Solar Energy Systems (Fraunhofer ISE) and the National Institute of Applied Science in Lyon, France (INSA Lyon) and with the support of EDF R&D. In this prototype, the physical configurations of the absorber, that include the water channels and the solar cells, are shown in Figure3. Four strings of eight pseudosquare sc-Si PV cells were connected in series and laminated with encapsulant and a polymer film onto the surface of a flat Rollbond[20] aluminum heat exchanger with an electric insulating coating. In spite of slightly lower thermal properties than copper, the lower price of aluminum enables a higher plate thickness to be considered (1 mm instead of 0.2 mm) in order to provide a better stiffness to the solar cells in the absence of glazing [21, 22]. The absorber was inserted into a collector frame. The collector front cover was a 4 mm AR glass with a transmission above 0.93 and the distance between the glass cover and the absorber was 20 mm. Thermal insulation material was then applied on the back side of the absorber. The main design data of this prototype is given in Table 1.Table 1
Key parameters of the PVT collector unit.
Front glazing
Thickness of glass
0.004 m
Emissivity of glass
0.9
Thermal capacity
840 kJ/(kg·K)
Density
2500 kg/m3
Depth of air gap underneath
0.02 mm
PV cell
Length of the PV cell
0.156 m
Width of the PV cell
0.156 m
Thickness
0.0002 m
Thermal conductivity
148 W/(m·K)
Emissivity
0.91
Absorptivity
0.93
Reference electrical efficiency of the sc-Si module
14%
Packing factor
0.79
Thermal absorber
Length of the module
1.350 m
Width of the module
0.73 m
Thermal capacity
910 kJ/(kg·K)
Density
2700 kg/m3
Absorptivity
0.94
Number of tubes
16
Thickness of insulation layer
0.05 m
Insulation thermal conductivity
0.04 W/(m·K)Description of the PV-T absorber plate rear side (a) and front side (b) [23].
(a)
(b)The simulation model is based on a nodal approach originally developed by Fraisse et al. [24, 25] and later modified by Dupeyrat et al. [26] according to the actual physical arrangements. 12 temperature nodes have been considered for the different interfaces from the cover to the backside of the collector including the collector fluid temperature. The nodal temperatures then can be obtained by solving the group of equations derived from the energy balance of each interface. This “TYPE” is different from the other models of solar thermal collectors in TRNSYS because it takes into account several important dynamic properties such as:(i)
the long-wave radiation exchange between the cover and the environment (exchanges between the ground and sky are separated),(ii)
the thermophysical characteristics of the collector (instead of the empirical correlation obtained from experiments).For the evaluation of the convective exchange in the sealed-air channel between the PV cell layer and the cover, the air properties are determined as a function of the air temperature (instead of using constant values). The coefficients of correlation are obtained from the experimental results.(i)
Energy balance of the cover(11)CCTCF-TCIΔt=K1(TAF-TCF)+KC1(Tamb-TCF)+KR1(Ts-TCF)+BcGAC,with:
(12)CC=ρC⋅eC⋅AC⋅CpC,K1=AC(hci+hri),KC1=AC(hc+hri),KR1=AChrs,TCI and TCF are the cover temperatures at the beginning and the end of every time step, K; ρC is the cover density, kg/m3; eC  is the thickness of the cover, m; CpCis the specific heat of the cover, kJ/(kg·K);hci is the convective heat transfer coefficient between the absorber and the cover, W/(m2·K);hri  is the radiative heat transfer coefficient between the absorber and the cover, W/(m2·K); hc is the convective heat transfer coefficient between the cover and ambience, W/(m2·K); hrs is the radiative heat transfer coefficient between the cover and the sky, W/(m2·K).(ii)
Energy balance of the absorber(13)CATAF-TAIΔt=K1(TCF-TAF)+K2(Tbs-TAF)+K3(Tf-TAF)+BGACwith:
(14)CA=ρA⋅eA⋅AC⋅CpA+ρT⋅CpT⋅π4⋅[(DT+2⋅eT)2-DT2]⋅LT⋅NT.
The heating capacity of the absorber including the tubes can be written as
(15)K2=ACeins/λins+1/hbs,K3=h0πDTLTNT,
where, TAI and TAF are the absorber temperatures at the beginning and the end of the time step respectively, K;B is the optical factor of the absorber; ρA is the absorber density, kg/m3; eA is the thickness of the absorber, m; CpA is the specific heat of the absorber, kJ/(kg·K); ρT is the tube density, kg/m3; eT is the thickness of the tube, m; CpT  is the specific heat of the tube, kJ/(kg·K); DT is the diameter of the tubes, m; LT is the length of the tubes, m; NT is the number of tubes; eins is the thickness of the insulation, m; λins is the thermal conductivity of the insulation, W/(m·K); Tbs is the collector backside temperature, °C; h0 is the convective heat transfer coefficient of the fluid, W/(m2·K).(iii)
Energy balance of the fluid(16)CfTfF-TfIΔt=K4(Tif-TfF)+K3(TAF-TfF),where,
(17)Cf=ρf⋅Cpf⋅πDT24⋅LT⋅NT,K4=mf˙⋅Cpf(1/(1-e-αLT)-1/αLT),α=h0πDT(mf˙⋅Cpf⋅NT).TfI and TfF are the fluid temperatures at the beginning and the end of the time step, K; ρf is the fluid density, kg/m; Cpfis the specific heat of the cover, kJ/kg·K); m˙f is the fluid mass flow rate, kg/s.The electrical efficiency is determined as follows:(18)EffPV=EffPV_ref⋅(1+βr(TPV-TPV_ref))⋅γ⋅τ1,
where γ is the PV cell packing factor, τ1 is the transmittance of the covers, EffPV_ref is the reference cell efficiency at the reference operating temperature TPV_ref=298.15 K; TPV is the PV cell surface temperature, K; βr is the temperature coefficient, %/K.At every time step, the exchange coefficients are to be calculated by using the initial temperatures of the nodes. The three differential equations, namely, (11), (13), and (16), which are defined according to the energy balance of every nodes, are solved numerically.The storage tank (Type 60) was modeled as a component of the vertical stratified cylindrical tank including internal heat exchanger. The thermal stratification can be modeled by assuming that the tank consists ofN (N⩽100) fully mixed equal volume segments. The degree of stratification is determined by the value of N. If N is equal to 1, the storage tank is modeled as a fully mixed tank and no stratification effects are possible.Listed in Tables2, 3  and 4 are the system parameters and input data considered in the simulation process.Table 2
Simulation parameters of sanitary equipment and required temperature.
Mass flow rate per shower equipment
300 L/h
Mass flow rate per lavabo
30 L/h
Number of the shower equipment
150
Number of the lavabo
80
Bathing water temperature
35°CTable 3
Simulation parameters of the solar water preheating system.
Collectors area
600 m2
Number of  collectors in series
10
Number of the array group in parallel
12
Mass flow rate
40 L/h·m2
Auxiliary heater power
5 kW
Tank model with internal coil exchanger
Horizontal cylinder tank
Orientation/tilt angle
South-facing/23°Table 4
Water source heat pump unit characteristics under the standard condition.
Type of refrigerant
R410A
Heat capacity
14.61 kW (40°C/45°C hot water, constant flow)
C.O.P
4.9
Cooling capacity
10.37 (12°C/7°C chilled water, constant flow)
E.E.R
4.5


## 5. Results and Discussion

Both the system energy and exergy performance were determined. The year round simulation results were obtained using the typical meteorological year (TMY) hourly weather data of Hong Kong. The simulation time step was 12 minutes.


### 5.1. Energy Evaluation Criteria

The rate of useful thermal energy produced from the PV/T collector can be evaluated as:(19)Qu=m˙fCp(Toc-Tic).
The electrical energy generated in kWh is given by:(20)Eelect=EffPV⋅Ac⋅G.
Annual thermal equivalent of electrical energy produced can be evaluated by using:
(21)Qelecannual=Eelecannual0.38.
An efficiency of 0.38 was arbitrarily used for the conversion of thermal energy into electrical energy for thermal power plants.The overall annual thermal output can be evaluated by:(22)Qt=Qelecannual+Quannual.
Collector thermal efficiency, which is the ratio of the useful energy gain Qu to the absorbed solar energy by the collectors, can be evaluated using (23), in that(23)Effthermal=Qu(ACG).
The fraction of energy saving is another evaluation criterion to estimate the system performance as compared to the conventional system by using the primary source energy. It can be determined by:
(24)Qsav=Qconvention-Qaux,
where, Qaux is the sum of the auxiliary consumption including the energy consumption of the auxiliary heater, the components (compressor, pumps, and, etc.) in the heat pump system:(25)Fsav=QsavQconvention.


### 5.2. Energy Evaluation Criteria

The exergy output is critical to define the real performance of PV/T system [25–27]. This can be calculated from the expression given by Fraisse et al.  [28], as follows:(26)Qexth=∑k=112∑j=1N∑i=1nQu˙(1-TaTsun),
where Ta is the ambient temperature in Kevin, and Tsun is temperature of the sun (5777 K).As the electrical output of a PV/W collector is a form of exergy, the total annual exergy output of the collector can be obtained by(27)Qexann=Eelect+Qexth.
The exergy from the solar radiation can be given by:
(28)Qexsolar=GA(1-TaTsun).
The exergy efficiency of the collector then can be calculated by
(29)Effexcol=QexannQexsolar.
The exergy efficiency of the overall system is then given as:
(30)Effexsys=m˙mCp(Tdis-Tfw)-Tam˙mln⁡(Tdis/Tfw)+EelectQexsolar+Qaux.


### 5.3. Daily and Monthly System Energy and Exergy Performance in Hong Kong

Figures4-5 show the variations of the collector inlet and outlet water temperatures (hot side of the heat exchanger), cold side inlet, and outlet temperatures of the heat exchanger connected to tank, delivery water temperature, and the outdoor air temperature of two typical sunny days, one in winter (mid February) and the other in summer (mid August) of Hong Kong. The mass flow rates in the solar loop and in the heat exchanger-tank loop are 2400 l/h and 2000 l/h, respectively. The fluid circulation on both sides is regulated by the corresponding temperature differential controller. During the system operating period, the temperature difference between the inlet and outlet of the solar collectors can be close to 10°C when there is water flow in the solar loop. As expected, the temperature of load water Tdis could be maintained around 40°C, which is slight higher than the expected pool water temperature on considering the heat losses in pipe distribution. The collector temperature is found higher in summer than in winter for the same solar radiation level.Figure 4
Daily variation of temperatures and the incident solar radiation, mi February.Figure 5
Daily variation of temperatures and the incident solar radiation, mi August.The monthly energy gain of the hybrid collectors is shown in Figure6. Qu represents the useful thermal output of the collectors. Epv represents the electricity generation from the collectors. It can be observed that the variation of the heat gain of the solar collectors agrees well with the incident solar radiation. The heat gain is low during the first quarter of the year as the solar radiation level is also relatively low in this period. The highest useful heat gain of 5.08E+04 kWh and the electricity generation of 9.57E+03 kWh as well are obtained in July. As shown in Figure 7, the collector thermal efficiency varies between 0.37 and 0.57 over the year. But the PV electrical efficiency is relatively constant (fluctuating between 10.1% and 10.7%). The year round thermal and electrical efficiencies are 49.3% and 10.3%, respectively. The combined energy efficiency is then 76.3%, with 0.38 as the conversion factor.Figure 6
Monthly thermal and electrical energy gain of arrays of solar collectors and the incident solar radiation.Figure 7
Monthly collector efficiencies.Figure8 shows the solar energy delivered to the tank Qhxc and the total energy auxiliary energy consumption of the system (auxiliary heater, heat pump components). It can be observed that the transferred solar energy via the heat exchanger is slightly lower than the useful thermal energy collected in the collector. This is because of the heat losses induced by the heat exchanger. The overall auxiliary energy consumption is more important for the winter season. This is when the small internal auxiliary heater is activated to keep the desired water temperature before the reheating by the heat pump system. Compared to the conventional electrical heating system, the monthly fraction of energy saving was around 65% over the winter period. This fraction became higher in the warm season and was around 90%.Figure 8
Monthly system auxiliary energy consumption and the energy saving.The monthly thermal performance of the heat pump series and the system COP are given in Figure9. Qhw represents the heat transfer to hot water; Qhp is the sum of the energy consumption of the compressor and the controller in the heat pump system. The energy consumption of the heat pump system is more important for the winter season when the water temperature exiting the storage tank is relatively low. More reheat energy is needed to push up to the desired delivery temperature. The system COP in the warm season is relatively low. This is so especially during the last three summer months when the temperature of heat source turns high, and this affects the heat pump performance.Figure 9
Monthly water source heat pump energy consumption and the system COP.


### 5.4. Comparison of Annual Energy and Exergy Performance under Different Climates

The annual energy and exergy performance of the designed system were evaluated under different climates, that is, in Hong Kong under subtropical climate and in three other cities of France. There are three climate zones in France, namely, Oceanic, Continental, and Mediterranean. The three cities Paris, Lyon and Nice, each could represent the typical climatic condition of each zone, were indicated on the map in Figure10. The selected four cities can represent also most of the similar region climate of the world. The numerical computation of the annual performance was performed based on the hourly TMY weather data of individual cities. All cities were assumed to have the same load profile as in Hong Kong and thereby creating a common basis for comparison. The feed water temperature was determined according to the fresh water temperature and the return load water temperature as shown in Figure 11.Figure 10
Climatic zones in France.Figure 11
Feed water temperature of different cities.Table5 lists the estimation of annual energy and exergy performance in the four cities. Regarding the total incident solar radiation, Nice has the highest level of solar radiation which yields also the highest outputs (useful thermal energy at 4.71E+05kWh and electricity at 1.02E+05kWh). Paris has the least energy gains as the solar radiation level and ambient temperature in this zone are relatively low. Consequently, more auxiliary energy was consumed in Paris. In terms of the system thermal output that affected by the freshwater temperature, Hong Kong has the lowest thermal output, at 1.21E+06kWh. Concerning the exergy output, Hong Kong and Paris are positioned at the same low level, whereas Lyon has the best exergy performance of 1.89E+05kWh.Table 5
Annual energy and exergy assessment of the designed system under different climates.
City
Solar radiation (kWh)
Collector useful thermal energy (kWh)
Collector electrical output (kWh)
Overall thermal output (kWh)
Exergy output (kWh)
Auxiliary energy consumption (kWh)
Hong Kong
8.74E + 05
4.31E + 05
9.00E + 04
1.21E + 06
1.29E + 05
4.19E + 05
Paris
6.74E + 05
2.93E + 05
7.24E + 04
1.63E + 06
1.29E + 05
6.26E + 05
Lyon
7.82E + 05
3.46E + 05
8.28E + 04
1.76E + 06
1.89E + 05
5.29E + 05
Nice
9.78E + 05
4.71E + 05
1.02E + 05
1.79E + 06
1.83E + 05
4.67E + 05The annual system efficiencies are given in Table6. The electrical efficiency for these four cities vary slightly between 10.29% and 10.75%. The collector thermal efficiency and the overall efficiency are relatively low for Paris and Lyon and are around 5% lower than Hong Kong—the city with the highest collector overall efficiency. Concerning the performance of the heat pump series, the system COP remains almost the same for all the cities, and the fraction of energy saving is all around 78%. Overall speaking, the system performs less desirable in exergy efficiency, at around 10% for all climatic conditions.Table 6
Annual efficiencies of the designed system.
City
Col. electrical efficiency (%)
Col. thermal efficiency (%)
Col. overall efficiency (%)
Exergy efficiency (%)
C.O.P of HP system (%)
Fraction of energy saving (%)
Hong Kong
10.29
49.24
76.32
10.06
4.11
66.55
Paris
10.75
43.45
71.73
9.85
4.34
67.98
Lyon
10.59
44.29
72.16
9.13
4.52
71.05
Nice
10.42
48.11
75.54
11.93
4.34
75.63


## 5.1. Energy Evaluation Criteria

The rate of useful thermal energy produced from the PV/T collector can be evaluated as:(19)Qu=m˙fCp(Toc-Tic).
The electrical energy generated in kWh is given by:(20)Eelect=EffPV⋅Ac⋅G.
Annual thermal equivalent of electrical energy produced can be evaluated by using:
(21)Qelecannual=Eelecannual0.38.
An efficiency of 0.38 was arbitrarily used for the conversion of thermal energy into electrical energy for thermal power plants.The overall annual thermal output can be evaluated by:(22)Qt=Qelecannual+Quannual.
Collector thermal efficiency, which is the ratio of the useful energy gain Qu to the absorbed solar energy by the collectors, can be evaluated using (23), in that(23)Effthermal=Qu(ACG).
The fraction of energy saving is another evaluation criterion to estimate the system performance as compared to the conventional system by using the primary source energy. It can be determined by:
(24)Qsav=Qconvention-Qaux,
where, Qaux is the sum of the auxiliary consumption including the energy consumption of the auxiliary heater, the components (compressor, pumps, and, etc.) in the heat pump system:(25)Fsav=QsavQconvention.


## 5.2. Energy Evaluation Criteria

The exergy output is critical to define the real performance of PV/T system [25–27]. This can be calculated from the expression given by Fraisse et al.  [28], as follows:(26)Qexth=∑k=112∑j=1N∑i=1nQu˙(1-TaTsun),
where Ta is the ambient temperature in Kevin, and Tsun is temperature of the sun (5777 K).As the electrical output of a PV/W collector is a form of exergy, the total annual exergy output of the collector can be obtained by(27)Qexann=Eelect+Qexth.
The exergy from the solar radiation can be given by:
(28)Qexsolar=GA(1-TaTsun).
The exergy efficiency of the collector then can be calculated by
(29)Effexcol=QexannQexsolar.
The exergy efficiency of the overall system is then given as:
(30)Effexsys=m˙mCp(Tdis-Tfw)-Tam˙mln⁡(Tdis/Tfw)+EelectQexsolar+Qaux.


## 5.3. Daily and Monthly System Energy and Exergy Performance in Hong Kong

Figures4-5 show the variations of the collector inlet and outlet water temperatures (hot side of the heat exchanger), cold side inlet, and outlet temperatures of the heat exchanger connected to tank, delivery water temperature, and the outdoor air temperature of two typical sunny days, one in winter (mid February) and the other in summer (mid August) of Hong Kong. The mass flow rates in the solar loop and in the heat exchanger-tank loop are 2400 l/h and 2000 l/h, respectively. The fluid circulation on both sides is regulated by the corresponding temperature differential controller. During the system operating period, the temperature difference between the inlet and outlet of the solar collectors can be close to 10°C when there is water flow in the solar loop. As expected, the temperature of load water Tdis could be maintained around 40°C, which is slight higher than the expected pool water temperature on considering the heat losses in pipe distribution. The collector temperature is found higher in summer than in winter for the same solar radiation level.Figure 4
Daily variation of temperatures and the incident solar radiation, mi February.Figure 5
Daily variation of temperatures and the incident solar radiation, mi August.The monthly energy gain of the hybrid collectors is shown in Figure6. Qu represents the useful thermal output of the collectors. Epv represents the electricity generation from the collectors. It can be observed that the variation of the heat gain of the solar collectors agrees well with the incident solar radiation. The heat gain is low during the first quarter of the year as the solar radiation level is also relatively low in this period. The highest useful heat gain of 5.08E+04 kWh and the electricity generation of 9.57E+03 kWh as well are obtained in July. As shown in Figure 7, the collector thermal efficiency varies between 0.37 and 0.57 over the year. But the PV electrical efficiency is relatively constant (fluctuating between 10.1% and 10.7%). The year round thermal and electrical efficiencies are 49.3% and 10.3%, respectively. The combined energy efficiency is then 76.3%, with 0.38 as the conversion factor.Figure 6
Monthly thermal and electrical energy gain of arrays of solar collectors and the incident solar radiation.Figure 7
Monthly collector efficiencies.Figure8 shows the solar energy delivered to the tank Qhxc and the total energy auxiliary energy consumption of the system (auxiliary heater, heat pump components). It can be observed that the transferred solar energy via the heat exchanger is slightly lower than the useful thermal energy collected in the collector. This is because of the heat losses induced by the heat exchanger. The overall auxiliary energy consumption is more important for the winter season. This is when the small internal auxiliary heater is activated to keep the desired water temperature before the reheating by the heat pump system. Compared to the conventional electrical heating system, the monthly fraction of energy saving was around 65% over the winter period. This fraction became higher in the warm season and was around 90%.Figure 8
Monthly system auxiliary energy consumption and the energy saving.The monthly thermal performance of the heat pump series and the system COP are given in Figure9. Qhw represents the heat transfer to hot water; Qhp is the sum of the energy consumption of the compressor and the controller in the heat pump system. The energy consumption of the heat pump system is more important for the winter season when the water temperature exiting the storage tank is relatively low. More reheat energy is needed to push up to the desired delivery temperature. The system COP in the warm season is relatively low. This is so especially during the last three summer months when the temperature of heat source turns high, and this affects the heat pump performance.Figure 9
Monthly water source heat pump energy consumption and the system COP.


## 5.4. Comparison of Annual Energy and Exergy Performance under Different Climates

The annual energy and exergy performance of the designed system were evaluated under different climates, that is, in Hong Kong under subtropical climate and in three other cities of France. There are three climate zones in France, namely, Oceanic, Continental, and Mediterranean. The three cities Paris, Lyon and Nice, each could represent the typical climatic condition of each zone, were indicated on the map in Figure10. The selected four cities can represent also most of the similar region climate of the world. The numerical computation of the annual performance was performed based on the hourly TMY weather data of individual cities. All cities were assumed to have the same load profile as in Hong Kong and thereby creating a common basis for comparison. The feed water temperature was determined according to the fresh water temperature and the return load water temperature as shown in Figure 11.Figure 10
Climatic zones in France.Figure 11
Feed water temperature of different cities.Table5 lists the estimation of annual energy and exergy performance in the four cities. Regarding the total incident solar radiation, Nice has the highest level of solar radiation which yields also the highest outputs (useful thermal energy at 4.71E+05kWh and electricity at 1.02E+05kWh). Paris has the least energy gains as the solar radiation level and ambient temperature in this zone are relatively low. Consequently, more auxiliary energy was consumed in Paris. In terms of the system thermal output that affected by the freshwater temperature, Hong Kong has the lowest thermal output, at 1.21E+06kWh. Concerning the exergy output, Hong Kong and Paris are positioned at the same low level, whereas Lyon has the best exergy performance of 1.89E+05kWh.Table 5
Annual energy and exergy assessment of the designed system under different climates.
City
Solar radiation (kWh)
Collector useful thermal energy (kWh)
Collector electrical output (kWh)
Overall thermal output (kWh)
Exergy output (kWh)
Auxiliary energy consumption (kWh)
Hong Kong
8.74E + 05
4.31E + 05
9.00E + 04
1.21E + 06
1.29E + 05
4.19E + 05
Paris
6.74E + 05
2.93E + 05
7.24E + 04
1.63E + 06
1.29E + 05
6.26E + 05
Lyon
7.82E + 05
3.46E + 05
8.28E + 04
1.76E + 06
1.89E + 05
5.29E + 05
Nice
9.78E + 05
4.71E + 05
1.02E + 05
1.79E + 06
1.83E + 05
4.67E + 05The annual system efficiencies are given in Table6. The electrical efficiency for these four cities vary slightly between 10.29% and 10.75%. The collector thermal efficiency and the overall efficiency are relatively low for Paris and Lyon and are around 5% lower than Hong Kong—the city with the highest collector overall efficiency. Concerning the performance of the heat pump series, the system COP remains almost the same for all the cities, and the fraction of energy saving is all around 78%. Overall speaking, the system performs less desirable in exergy efficiency, at around 10% for all climatic conditions.Table 6
Annual efficiencies of the designed system.
City
Col. electrical efficiency (%)
Col. thermal efficiency (%)
Col. overall efficiency (%)
Exergy efficiency (%)
C.O.P of HP system (%)
Fraction of energy saving (%)
Hong Kong
10.29
49.24
76.32
10.06
4.11
66.55
Paris
10.75
43.45
71.73
9.85
4.34
67.98
Lyon
10.59
44.29
72.16
9.13
4.52
71.05
Nice
10.42
48.11
75.54
11.93
4.34
75.63


## 6. Economic Analysis

The investment on a combined solar-assisted heat pump system includes the material and labour costs of the solar collectors, thermal storage tank, heat pump series, accessories (such as heat exchanger, pumps, controller, valves and, etc.) and the relevant system test and transportation. By making reference to the expenses of other similar SAHP systems in Hong Kong and in France, the costs of the entire system were evaluated and shown in Table7. The most important investment item lies in the solar collectors. In France, according to the low carbon promotion policy, every installation of photovoltaic system could benefit by 11% tax credit [29] based on the total system investment. In addition, there exist some kinds of local subvention depending on the region and relevant department. In this case study, we took on average 2% of the total investment as the sum of the subvention. The average cost per kWh of electricity in Hong Kong is 1 HK$ and this is 0.1074 € (equivalent 1.1 HK$) [30] in France. The cost payback period was then estimated for each city. The results are listed in Table 7. Hong Kong has the longest payback period of 10.52 years even though the cost of material is less expensive compared to the other French cities. So far in Hong Kong, there is no subvention formulated for the solar thermal or electrical installations. The payback periods for the French cities are all around 9 years. In Nice, the city which possesses the richest annual solar radiation, this could be shortened to 8.4 years. However, the French government lowers the tax reduction for photovoltaic system year by year (25% for 2010, 22% for 2011 and 11% for 2012). The payback period will be probably longer in future, if the technical advancements are not to be considered.Table 7
Investment of the designed system.
System configuration
Component
Hong Kong
Paris
Lyon
Nice
Collector area (m2)
600
600
600
600
Heat pump heating capacity (kW)
175
270
240
220
Power of auxiliary heater unit (kW)
15
15
50
50
Storage tank (m3)
60
60
60
60
Cost (HK$)
Solar collectors
5.10E + 06
6.00E + 06
6.00E + 06
6.00E + 06
Water tank
4.02E + 05
2.40E + 06
2.40E + 06
2.40E + 06
Auxiliary heater
5.94E + 03
1.05E + 04
1.05E + 04
1.05E + 04
Heat pump systems
2.78E + 05
3.20E + 06
2.56E + 06
2.40E + 06
Accessories
5.38E + 05
9.20E + 05
8.56E + 05
8.40E + 05
Transportation
2.18E + 04
2.24E + 05
1.98E + 05
1.92E + 05
Cost of Balance of system plus site installation and testing
1.50E + 04
4.60E + 05
4.28E + 05
4.20E + 05
Total investment
6.36E + 06
1.32E + 07
1.25E + 07
1.23E + 07
Tax reduction and governmental subvention (percentage of total investment)
NA
13%
13%
13%
Energy saving (kWh)
6.72E + 05
1.19E + 06
1.12E + 06
1.16E + 06
Payback period
10.52
8.82
8.83
8.40


## 7. Conclusions

This study was intended to investigate the potential application of the hybrid PV/T solar-assisted heat pump system for indoor sports centre water heating under the subtropical climatic condition of Hong Kong, as well as the multiple climates in France. An initial design procedure was first presented, and the energy performance was evaluated through the use of the TRNSYS transient simulation tool. The numerical results showed that the designed system can well satisfy the system energy demands. The hot water supply temperature can be reheated to 40°C during the operating period. The mean heating COP of the SAHP system is found to be 4.3. The overall efficiency of PV/T collectors can reach 76% in Hong Kong. The global fractional energy saving factor can achieve a maximum of 75% in Nice. The payback period of 10.52 years for Hong Kong is comparatively long. The stake of the further research led on this topic should be on the improvement of system exergy efficiency by optimizing the system configuration and parametric studies. The system life cycle assessment should also be in place to assess the environmental impacts.



---

*Source: 265838-2012-08-09.xml*