Theoretical Investigation of the Temperature Limits of an Actively Cooled High Concentration Photovoltaic System

: Concentrator photovoltaics have several advantages over ﬂat plate systems. However, the increase in solar concentration usually leads to an increase in the solar cell temperature, which decreases the performance of the system. Therefore, in this paper, we investigate the performance and temperature limits of a high concentration photovoltaic Thermal system (HCPVT) based on a 1 cm 2 multi-junction solar cell subjected to a concentration ratio from 500 × to 2000 × by using three di ﬀ erent types of cooling ﬂuids (water, ethylene glycol and water mixture (60:40), and syltherm oil 800). The results show that, for this conﬁguration, the maximum volumetric temperature of the solar cell did not exceed the manufacturer’s recommended limit for the tested ﬂuids. At 2000 × the lowest solar cell temperature obtained by using water was 93.5 ◦ C, while it reached as high as 109 ◦ C by using syltherm oil 800, which is almost equal to the maximum operating limit provided by the manufacturer (110 ◦ C). Overall, the best performance in terms of temperature distribution, thermal, and electrical e ﬃ ciency was achieved by using water, while the highest outlet temperature was obtained by using syltherm oil 800.


Introduction
Concentrator photovoltaics (CPV) has attracted a great deal of attention in the field of electricity generation due to its competitive cost in comparison with conventional PV systems for areas with high solar irradiance [1]. The key advantage of the CPV system is the use of high-efficiency concentrator cells, which can reach as high as 44.4% and 46% for three and four multi-junction (MJ) cells at 25 • C, respectively [2]. A key challenge, however, is the increase in the temperature of the solar cell, which adversely affects its electrical efficiency. In light of this challenge, the MJ solar cell manufacturer (AZUR SPACE) has set a maximum operating temperature for the CPV assembly to be 110 • C in order to avoid the damage of the solar cell [3]. Araki et al. [4] indicated that the solar cell temperature increases rapidly with the concentration ratio and can reach up to 1400 • C for 500× when left fully insulated.
Both passive and active cooling mechanisms can be deployed to maintain the solar cell temperature. However, an active cooling system offers a great advantage over the passive option as it allows the use of the waste heat from the system in other thermal applications, which increases the total efficiency of the system as well. Aldossary et al. [5] studied an active cooling system for a high concentrator photovoltaic (500×) for a 1 cm 2 MJ solar cell. Their system considered a rectangular cooling channel. The maximum solar cell temperature reached 65 • C at 0.01 m/s. Alsiyabi et al. [6] investigated the performance of the CPV system at 500× using a microchannel heat sink. They reported that the maximum cell temperature reached 54 • C using four layers of a microchannel heat sink, which adds to the cost of the system. Abo-Zahhad et al. [7] studied a CPV system subjected to a concentration ratio of 1000× using a water impingement jet. The researchers revealed that the solar cell temperature fluctuates between 55 and 67 • C through the use of water. In another study by the same researchers [8], a high concentration photovoltaic system (HCPV) under passive cooling conditions was considered. The maximum temperature of the solar cell reached 160 • C at 1500× and convective heat flux of 1600 W/m 2 K However, there is no reported research about the performance of an HCPV using the active cooling method for a concentration ratio higher than 500× for a 1 cm 2 MJ solar cell with the possibility of replacing water as a coolant with an ethylene glycol and water mixture or syltherm oil 800 for this system to avoid the mineral sedimentation issues of water on the fins at high temperatures [9]. Besides, they offer high boiling temperatures if compared with water, which gives them another advantage. Therefore, in this paper, we investigate the performance and temperature limits of HCPV using a mini-channel heat sink under different concentration ratios (500×-2000×). In addition, we present a comparison of solar cell temperatures when using water, ethylene glycol and water mixture (60:40), and syltherm oil 800 under a different range of mass flow rates. Performance characteristics are shown, categorized in terms of maximum solar cell temperature, outlet fluid temperature, and the thermal and electrical efficiencies of the system.

Mathematical Modeling
The 3D computational domain, shown in Figure 1, consists of: (1) 1 cm 2 MJ solar cell assembly; (2) ceramic layer; (3) copper layers; and (4) aluminum mini-channel heat sink. Their thermal properties and dimensions are presented in Tables 1 and 2. The germanium layer is considered the heat source of the system [10]. The total input energy (Q in ) is converted into electrical power (P sc,elec ) and the remainder is transformed into thermal heat Q sc,heat . The performance of the system is compared using different working media, water, ethylene glycol and water mixture (60:40), and syltherm oil 800 under different ranges of mass flow rates. The trial and error method is used to evaluate the temperature of the solar cell (T sc ) using the following relation [11,12]: where η sc and η ref are the solar cell efficiency at a certain temperature and reference cell efficiency at an ambient temperature of 25 • C, respectively. β ref is the temperature coefficient of the MJ solar cell and is equal to 0.0047 %/K as provided by the manufacturer [13] and T ref is the reference temperature of the solar cell, which is equal to the ambient temperature of 25 • C.
Energies 2020, 13, x FOR PEER REVIEW 2 of 10 reported that the maximum cell temperature reached 54 °C using four layers of a microchannel heat sink, which adds to the cost of the system. Abo-Zahhad et al. [7] studied a CPV system subjected to a concentration ratio of 1000× using a water impingement jet. The researchers revealed that the solar cell temperature fluctuates between 55 and 67 � through the use of water. In another study by the same researchers [8], a high concentration photovoltaic system (HCPV) under passive cooling conditions was considered. The maximum temperature of the solar cell reached 160 °C at 1500× and convective heat flux of 1600 W m 2 K ⁄ . However, there is no reported research about the performance of an HCPV using the active cooling method for a concentration ratio higher than 500× for a 1 cm 2 MJ solar cell with the possibility of replacing water as a coolant with an ethylene glycol and water mixture or syltherm oil 800 for this system to avoid the mineral sedimentation issues of water on the fins at high temperatures [9]. Besides, they offer high boiling temperatures if compared with water, which gives them another advantage. Therefore, in this paper, we investigate the performance and temperature limits of HCPV using a mini-channel heat sink under different concentration ratios (500×-2000×). In addition, we present a comparison of solar cell temperatures when using water, ethylene glycol and water mixture (60:40), and syltherm oil 800 under a different range of mass flow rates. Performance characteristics are shown, categorized in terms of maximum solar cell temperature, outlet fluid temperature, and the thermal and electrical efficiencies of the system.

Mathematical Modeling
The 3D computational domain, shown in Figure 1, consists of: (1) 1 cm 2 MJ solar cell assembly; (2) ceramic layer; (3) copper layers; and (4) aluminum mini-channel heat sink. Their thermal properties and dimensions are presented in Tables 1 and 2. The germanium layer is considered the heat source of the system [10]. The total input energy �Q in � is converted into electrical power �P sc,elec � and the remainder is transformed into thermal heat �Q sc,heat �. The performance of the system is compared using different working media, water, ethylene glycol and water mixture (60:40), and syltherm oil 800 under different ranges of mass flow rates. The trial and error method is used to evaluate the temperature of the solar cell (T sc ) using the following relation [11,12]: where η sc and η ref are the solar cell efficiency at a certain temperature and reference cell efficiency at an ambient temperature of 25 °C, respectively. β ref is the temperature coefficient of the MJ solar cell and is equal to 0.0047 %/K as provided by the manufacturer [13] and T ref is the reference temperature of the solar cell, which is equal to the ambient temperature of 25 °C. The total input energy on the surface of the solar cell can be represented as: The thermal heat and the electrical power produced from the cell are calculated as follows:   The total input energy on the surface of the solar cell can be represented as: The thermal heat and the electrical power produced from the cell are calculated as follows: P sc,elec = Q in η sc (4) where DNI, CR, and A cell are the direct solar irradiance (W/m 2 ), concentration ratio, and solar cell area (m 2 ), respectively. Hence, the thermal energy to be carried by the fluid: where · m, and c p,f are the mass flow rate (kg/s) and the specific heat capacity of the working medium (J/kg K). T f,in and T f,out are the inlet and outlet fluid temperatures ( • C), respectively. Therefore, the thermal efficiency, η th , can be calculated by the following relation: Then, the overall efficiency of the HCPVT system can be calculated as: where P elec,net and η opt are the net electrical power (W) and optical efficiency: P elec,net = P sc,elec − P pump (8) Energies 2020, 13,1902 4 of 10 The pumping power, P pump , can be calculated by the following relation: where · V and ∆P are the volume flow rate m 3 and pressure drop (Pa), respectively. The thermophysical properties of all the tested fluids are temperature dependent. The properties of water that were used during the computation were adapted from [14]: In the case of ethylene glycol and water mixture (60:40), the following equations were used. For syltherm oil 800, the Loikits data [15] were adapted: The computational domain ( Figure 2) was discretized and solved using the Finite Element Method in COMSOL-Multiphysics 5.3 software [16]. Using the Generalized Minimum Residual (GMR) method, the 3D conjugated heat transfer equations were solved iteratively under steady-state conditions. The top surface of the solar cell, copper, and ceramic layers are subjected to natural convection and radiation heat transfer. Therefore, the heat balance equation to describe the heat transfer on the top surface of the solar cell can be written as: The energy equation: The 3D computational domain was solved using the Finite Element Method (FEM) according to these assumptions:


The MJ solar cell is subjected to concentration ratios varying from 500× to 2000×.  The germanium layer is the heat source of the system.  The top surfaces of the solar cell, copper, and ceramic layers, as presented in Figure 2, are subjected to natural convection (h = 15 W (m 2 K) ⁄ ) [12] and surface to ambient radiation.  The sides of the module and its back surface are insulated.  The ambient temperature is fixed at 25 °C.  The flow inside the channel is steady, incompressible, and laminar.  The effects of the body force and viscous dissipation are ignored.  The inlet temperature of the fluid is fixed at 25 °C and the outlet pressure is atmospheric.
The governing equations were solved using for the following boundary conditions: At the channel inlet: At the channel outlet: The volumetric heat source in the germanium layer: where V is the volume of the germanium layer.

Grid Independence Study
The first step to ensure the accuracy of the modeling was to carry out a grid independence study. The software allows the use of unstructured mesh (free tetrahedral) to solve the computational domain. In Figure 3a, the variation of the volumetric maximum solar cell temperature is plotted against the grid size. The number of elements was varied from 0.2 to 1.9 million elements. The maximum temperature has an almost constant value from 1.5 to 1.9 million elements. Therefore, the best number of elements, 1.5 million, in terms of accuracy and computation time was selected and  The heat balance equations for both copper and ceramic layers are as follows: The convective and radiative heat transfer can be calculated from the following relations [17]: where h is the convective heat transfer coefficient (W / m 2 K ), T s is the exposed surface temperature ( • C), T ∞ is the ambient temperature and is assumed to be 25 • C, ε is the emissivity of the surfaces that are subjected to radiation, and σ is the Stefan-Boltzmann constant. For radiation heat transfer, the sky temperature was assumed to be equal to the ambient temperature as introduced in Equation (22) for simplification, as suggested in the literature [3,17]. For the MJ solar cell domain, the law of energy conversion for a steady-state with a heat source term · q can be written as follows: Further, for the heat sink, the heat conduction equation should be: For the fluid domain, the continuity equation: The momentum equation: The energy equation: The 3D computational domain was solved using the Finite Element Method (FEM) according to these assumptions: • The MJ solar cell is subjected to concentration ratios varying from 500× to 2000×.

•
The germanium layer is the heat source of the system.

•
The top surfaces of the solar cell, copper, and ceramic layers, as presented in Figure 2, are subjected to natural convection h = 15 W/ m 2 K [12] and surface to ambient radiation.
• The sides of the module and its back surface are insulated.

•
The ambient temperature is fixed at 25 • C.

•
The flow inside the channel is steady, incompressible, and laminar.

•
The effects of the body force and viscous dissipation are ignored.

•
The inlet temperature of the fluid is fixed at 25 • C and the outlet pressure is atmospheric.
The governing equations were solved using for the following boundary conditions: At the channel inlet: Energies 2020, 13,1902 6 of 10 u x=0 = u in (29) At the channel outlet: P x=l = P atm (30) The volumetric heat source in the germanium layer: where V is the volume of the germanium layer.

Grid Independence Study
The first step to ensure the accuracy of the modeling was to carry out a grid independence study. The software allows the use of unstructured mesh (free tetrahedral) to solve the computational domain. In Figure 3a, the variation of the volumetric maximum solar cell temperature is plotted against the grid size. The number of elements was varied from 0.2 to 1.9 million elements. The maximum temperature has an almost constant value from 1.5 to 1.9 million elements. Therefore, the best number of elements, 1.5 million, in terms of accuracy and computation time was selected and the generated mesh is presented in Figure 3b. After generating the mesh and running the simulation, the mass and energy residuals were calculated using the software to ensure that they have a value lower than the relative tolerance error set in the software (0.001).
Energies 2020, 13, x FOR PEER REVIEW 6 of 10 the generated mesh is presented in Figure 3b. After generating the mesh and running the simulation, the mass and energy residuals were calculated using the software to ensure that they have a value lower than the relative tolerance error set in the software (0.001).

Validation Study
To ensure the accuracy of the numerical modeling further, several validation studies were conducted. The results of the present model were compared with the previously published results of Al Siyabi et al. [12] and Theristis et al. [17] at the same dimension and operating conditions. A good agreement was found between the present results and the compared studies in terms of maximum solar cell temperature with a maximum error of 2.4% with Al Siyabi et al. [12] and 0.7% with Theristis et al. [17], as presented in Figure 4a

Validation Study
To ensure the accuracy of the numerical modeling further, several validation studies were conducted. The results of the present model were compared with the previously published results of Al Siyabi et al. [12] and Theristis et al. [17] at the same dimension and operating conditions. A good agreement was found between the present results and the compared studies in terms of maximum solar cell temperature with a maximum error of 2.4% with Al Siyabi et al. [12] and 0.7% with Theristis et al. [17], as presented in Figure 4a,b.
To ensure the accuracy of the numerical modeling further, several validation studies were conducted. The results of the present model were compared with the previously published results of Al Siyabi et al. [12] and Theristis et al. [17] at the same dimension and operating conditions. A good agreement was found between the present results and the compared studies in terms of maximum solar cell temperature with a maximum error of 2.4% with Al Siyabi et al. [12] and 0.7% with Theristis et al. [17], as presented in Figure 4a

Results and Discussion
This section investigates the impact of the concentration ratio and the mass flow rate of water, ethylene glycol and water mixture (60:40), and syltherm-oil 800 on the maximum volumetric solar cell temperature, temperature distribution on the solar cell surface, outlet fluid temperature, and thermal and electrical efficiency of the system.

Results and Discussion
This section investigates the impact of the concentration ratio and the mass flow rate of water, ethylene glycol and water mixture (60:40), and syltherm-oil 800 on the maximum volumetric solar cell temperature, temperature distribution on the solar cell surface, outlet fluid temperature, and thermal and electrical efficiency of the system.

Maximum Solar Cell Temperature and Temperature Distribution over the Surface of the Solar Cell
The impact of using water, ethylene glycol and water mixture (60:40), and syltherm oil 800 on the performance of the HCPVT system at different concentration ratios is presented in this section. Figure 5 compares the variation of the maximum volumetric solar cell temperature with the mass flow rate for water, ethylene glycol and water mixture, and syltherm oil. The highest temperature was observed by using oil (109.1 • C), while the lowest was by using water (93.5 • C) at the same mass flow rate and solar concentration ratio of 2000×. Increasing the mass flow rate had a positive effect on reducing the maximum volumetric solar cell temperature, reaching 34.7 • C at a solar concentration ratio of 500×. In general, all the calculated maximum temperatures were lower than the recommended limit provided by the manufacturer. In addition, the temperature distribution for the different working media is presented in Figure 6. Water offered the most uniform temperature distribution followed by ethylene glycol and water mixture (60:40), the average temperature varying from 82.16 to 88.48 • C, respectively.

Maximum Solar Cell Temperature and Temperature Distribution over the Surface of the Solar Cell
The impact of using water, ethylene glycol and water mixture (60:40), and syltherm oil 800 on the performance of the HCPVT system at different concentration ratios is presented in this section. Figure 5 compares the variation of the maximum volumetric solar cell temperature with the mass flow rate for water, ethylene glycol and water mixture, and syltherm oil. The highest temperature was observed by using oil (109.1 �), while the lowest was by using water (93.5 �) at the same mass flow rate and solar concentration ratio of 2000×. Increasing the mass flow rate had a positive effect on reducing the maximum volumetric solar cell temperature, reaching 34.7 °C at a solar concentration ratio of 500× . In general, all the calculated maximum temperatures were lower than the recommended limit provided by the manufacturer. In addition, the temperature distribution for the different working media is presented in Figure 6. Water offered the most uniform temperature distribution followed by ethylene glycol and water mixture (60:40), the average temperature varying from 82.16 to 88.48 °C, respectively.

Outlet Fluid Temperature
The variation of the outlet temperatures of the fluids with the concentration ratio and mass flow rate are presented in Figure 7. The outlet temperature decreased with increasing the mass flow rate, while it increased with the concentration ratio. The highest outlet temperature was observed in the case of using the syltherm oil (73.5 °C), while the lowest was in the case of using water (58.5 °C) at a solar concentration of 2000× and a mass flow rate of 6.74 × 10 −4 kg/s. Further, the outlet temperature by using ethylene glycol and water mixture reached as high as 65 °C at the lowest tested mass flow rate at the concentration ratio of 2000× and 34 °C at a concentration ratio of 500×. Therefore, increasing the concentration ratio had a positive impact on the outlet temperature for all the tested fluids. Syltherm oil 800 offered the highest outlet fluid temperature in comparison with other fluids especially for a concentration ratio above 1000×. Overall, deciding which fluid to use depends mainly on the designed concentration ratio and the available heat recovery application.

Outlet Fluid Temperature
The variation of the outlet temperatures of the fluids with the concentration ratio and mass flow rate are presented in Figure 7. The outlet temperature decreased with increasing the mass flow rate, while it increased with the concentration ratio. The highest outlet temperature was observed in the case of using the syltherm oil (73.5 • C), while the lowest was in the case of using water (58.5 • C) at a solar concentration of 2000× and a mass flow rate of 6.74×10 −4 kg/s. Further, the outlet temperature by using ethylene glycol and water mixture reached as high as 65 • C at the lowest tested mass flow rate at the concentration ratio of 2000× and 34 • C at a concentration ratio of 500×. Therefore, increasing the concentration ratio had a positive impact on the outlet temperature for all the tested fluids. Syltherm oil 800 offered the highest outlet fluid temperature in comparison with other fluids especially for a concentration ratio above 1000×. Overall, deciding which fluid to use depends mainly on the designed concentration ratio and the available heat recovery application.

Electrical and Thermal Efficiencies
The solar cell temperature and the electrical efficiency were calculated by equation (1) using the trial and error method [11,12] by inputting the reference efficiency as provided by the manufacturer. Therefore, the electrical efficiency depends on the concentration ratio as well as the solar cell temperature. The electrical efficiency of the HCPVT system increased significantly with an increase in the mass flow rate of the working media, as presented in Figure 8. Using syltherm oil decreased the electrical efficiency if compared with water by 2.2% at 6.74 × 10 −4 kg/s, especially at a concentration ratio of 2000×. This can be attributed to the high surface temperature of the solar cell that is detected in the case of using the syltherm oil if compared with the other fluids at the same mass flow rates. The highest electrical efficiency was 40.8% in the case of using water and a concentration ratio of 500× and mass flow rate of 0.006 kg/s, while the lowest efficiency was 31.2% in the case of using syltherm oil and a concentration ratio of 2000× and mass flow rate of 6.74 × 10 −4 kg/s. Further, increasing the concentration ratio had a negative impact on the calculated electrical efficiency as a large portion of the concentrated solar irradiance converted into heat. In total, among the tested fluids water showed the highest performance in terms of electrical efficiency especially at concentration ratios above 1000×. Figure 9 shows the effect of increasing the mass flow rate on the thermal efficiency of the HCPVT system at different concentration ratios. The thermal efficiency of the HCPVT system increased with an increase in the mass flow rate and concentration ratio. The thermal efficiency varied from 42% to 65.4% for the concentration ratio between 500× and 2000×.
The lowest values are observed in the case of using syltherm oil 800 due to its lower specific heat at

Electrical and Thermal Efficiencies
The solar cell temperature and the electrical efficiency were calculated by equation (1) using the trial and error method [11,12] by inputting the reference efficiency as provided by the manufacturer. Therefore, the electrical efficiency depends on the concentration ratio as well as the solar cell temperature. The electrical efficiency of the HCPVT system increased significantly with an increase in the mass flow rate of the working media, as presented in Figure 8. Using syltherm oil decreased the electrical efficiency if compared with water by 2.2% at 6.74×10 −4 kg/s, especially at a concentration ratio of 2000×. This can be attributed to the high surface temperature of the solar cell that is detected in the case of using the syltherm oil if compared with the other fluids at the same mass flow rates. The highest electrical efficiency was 40.8% in the case of using water and a concentration ratio of 500× and mass flow rate of 0.006 kg/s, while the lowest efficiency was 31.2% in the case of using syltherm oil and a concentration ratio of 2000× and mass flow rate of 6.74×10 −4 kg/s. Further, increasing the concentration ratio had a negative impact on the calculated electrical efficiency as a large portion of the concentrated solar irradiance converted into heat. In total, among the tested fluids water showed the highest performance in terms of electrical efficiency especially at concentration ratios above 1000×. Figure 9 shows the effect of increasing the mass flow rate on the thermal efficiency of the HCPVT system at different concentration ratios. The thermal efficiency of the HCPVT system increased with an increase in the mass flow rate and concentration ratio. The thermal efficiency varied from 42% to 65.4% for the concentration ratio between 500× and 2000×. The lowest values are observed in the case of using syltherm oil 800 due to its lower specific heat at constant pressure if compared with the other fluids. The thermal efficiency increased from 23.97% to 57.45% when using syltherm oil 800, while it varied from 34.65% to 62.8% when using the ethylene glycol and water mixture. In the case of thermal efficiency, increasing the concentration ratio had a positive impact as shown in the results. mass flow rates. The highest electrical efficiency was 40.8% in the case of using water and a concentration ratio of 500× and mass flow rate of 0.006 kg/s, while the lowest efficiency was 31.2% in the case of using syltherm oil and a concentration ratio of 2000× and mass flow rate of 6.74 × 10 −4 kg/s. Further, increasing the concentration ratio had a negative impact on the calculated electrical efficiency as a large portion of the concentrated solar irradiance converted into heat. In total, among the tested fluids water showed the highest performance in terms of electrical efficiency especially at concentration ratios above 1000×. Figure 9 shows the effect of increasing the mass flow rate on the thermal efficiency of the HCPVT system at different concentration ratios. The thermal efficiency of the HCPVT system increased with an increase in the mass flow rate and concentration ratio. The thermal efficiency varied from 42% to 65.4% for the concentration ratio between 500× and 2000×.
The lowest values are observed in the case of using syltherm oil 800 due to its lower specific heat at constant pressure if compared with the other fluids. The thermal efficiency increased from 23.97% to 57.45% when using syltherm oil 800, while it varied from 34.65% to 62.8% when using the ethylene glycol and water mixture. In the case of thermal efficiency, increasing the concentration ratio had a positive impact as shown in the results.

Conclusions
This study explored the performance of an HCPVT system using MJ solar cell under active cooling conditions and solar concentration ratios varying from 500× to 2000× . In addition, we investigated the possibility of replacing water with ethylene glycol and water mixture (60:40) and syltherm oil 800. In terms of solar cell temperature, water offered the most uniform temperature distribution and the lowest volumetric temperature at the examined flow rates and concentration ratios. Further, in terms of the outlet temperature, syltherm oil 800 showed the highest fluid temperature, especially for a concentration ratio above 1000×, which would be suitable for thermal applications that require a high temperature. At concentration ratios above 1000×, using syltherm oil reduced the electrical efficiency by more than 1% when compared with water, while at lower concentration ratios the reduction was insignificant. Water achieved the highest thermal efficiency due to its high specific heat followed by ethylene glycol and water mixture and lastly syltherm oil 800. Overall, deciding which fluid to be used depends mainly on the designed concentration ratio, the required outlet temperature, and thermal energy to be recovered. For future implementation, to reflect the outdoor conditions, a transient simulation would be better than the steady-state case. For more precise calculations, a coupled electrical and thermal simulation would be of importance.

Conclusions
This study explored the performance of an HCPVT system using MJ solar cell under active cooling conditions and solar concentration ratios varying from 500× to 2000×. In addition, we investigated the possibility of replacing water with ethylene glycol and water mixture (60:40) and syltherm oil 800. In terms of solar cell temperature, water offered the most uniform temperature distribution and the lowest volumetric temperature at the examined flow rates and concentration ratios. Further, in terms of the outlet temperature, syltherm oil 800 showed the highest fluid temperature, especially for a concentration ratio above 1000×, which would be suitable for thermal applications that require a high temperature. At concentration ratios above 1000×, using syltherm oil reduced the electrical efficiency by more than 1% when compared with water, while at lower concentration ratios the reduction was insignificant. Water achieved the highest thermal efficiency due to its high specific heat followed by ethylene glycol and water mixture and lastly syltherm oil 800. Overall, deciding which fluid to be used depends mainly on the designed concentration ratio, the required outlet temperature, and thermal energy to be recovered. For future implementation, to reflect the outdoor conditions, a transient simulation would be better than the steady-state case. For more precise calculations, a coupled electrical and thermal simulation would be of importance. Funding: The Newton-Mosharafa Fund (a UK-Egypt partnership program for postgraduate scholarships) is acknowledged for the financial support given to Asmaa Ahmed.

Conflicts of Interest:
The authors declare no conflict of interest.