Investigation of Solar Hybrid Electric / Thermal System with Radiation Concentrator and Thermoelectric Generator

An experimental study of a solar-concentrating system based on thermoelectric generators (TEGs) was performed. The system included an electrical generating unit with 6 serially connected TEGs using a traditional semiconductormaterial, Bi 2 Te 3 , which was illuminated by concentrated solar radiation on one side and cooled by running water on the other side. A sun-tracking concentrator with amosaic set of mirrors was used; its orientation towards the sun was achieved with two pairs of radiation sensors, a differential amplifier, and two servomotors. The hot side of the TEGs at midday has a temperature of around 200C, and the cold side is approximately 50C. The thermosiphon cooling system was designed to absorb the heat passing through the TEGs and provide optimal working conditions. The system generates 20W of electrical energy and 200W of thermal energy stored in water with a temperature of around 50C. The hybrid system studied can be considered as an alternative to photovoltaic/thermal systems, especially in countries with abundant solar radiation, such as Mexico, China, and India.


Introduction
Solar hybrid electric/thermal systems using photovoltaic (PV) panels combined with a water/air-filled heat extracting unit were designed and studied in many laboratories during the last three decades [1][2][3][4][5][6][7][8][9][10] and now are widely used throughout the world (England, Canada, China, Greece, India, and so forth).Some investigations were made [11][12][13][14][15][16] into the possibilities of using thermoelectric generators (TEGs) in solar hybrid systems, with the conclusion that TEGs can be successfully used in these systems instead of PV panels or together with them.An essential increase in thermoelectric conversion efficiency was achieved during the last decade [17][18][19], which is quite favorable for this kind of TEG's applications.With the traditional thermoelectric material Bi 2 Te 3 , the peak electric efficiency that could be obtained in such a system is 5% [16].
Chávez-Urbiola et al. [14] investigated different options of the construction of hybrid solar energy conversion systems using TEGs.They showed that these systems can be efficient (and economic in case of industrial production) even with the use of material and devices that are already available on the market, especially in countries with high solar insolation (Mexico, China, India, etc.).Below we describe the construction and detailed experimental investigation of one of the hybrid systems analyzed in the above-mentioned paper [14], namely, the system with a solar radiation concentrator, TEG, and water-filled heat extracting unit.Circulation of water was achieved by thermosiphon effect.The experiments were performed in Queretaro, Mexico, at 20 ∘ of northern latitude, in March 2012.

Description of the Hybrid System
A schematic of the system is shown in Figure 1 where the solar radiation flux (1) is concentrated by the mosaic mirror (2) onto the electric/thermal generating unit (3) (details of the TEG are shown in Figure 2) consisting of a radiation absorber (hot plate), TEG array, and a cooling plate that is in direct contact with water-circulating copper tubes.The thermosiphon water loop includes a water storage thermal tank (4) with tubes for water entrance and output.
The radiation-concentrating block consisted of 55 plane mirrors each having a size equal to that of the TEG array (8 × 12 cm 2 ), providing a concentration ratio (the number of mirrors focused on the heating plate multiplied by the mirror reflecting efficiency) of ∼52, and considering a reflection efficiency of 0.95.The mirrors were positioned in a parabolic curve, with the focal point over the heating plate of the TEG assembly; the angle of the inclination of each mirror was calculated to achieve this effect.The block (mirror holder) was attached to the 2-axis sun-tracking system (see [20] for details) equipped with 2 pairs of radiation sensors positioned in such a way that the difference in photo response in each couple is zero if the mirror holder is orientated towards the sun, giving the highest illumination of the absorber hot plate.The difference in photo response increases with disorientation (disorientation signal).This difference signal is applied to a PIC16F877 microcontroller, which monitors the system using two geared servomotors.
The TEG array includes 6 generating elements of the type TGM-127-1.4-2.5 based on Bi 2 Te 3 (made by Kryotherm, Saint Petersburg, Russia; each element is 4 × 4 × 0.5 cm 3 ).The electrical characteristics of the elements at different temperatures of the operation were given in a previous publication [14].

Calculation of Thermosiphon Loop
For the thermosiphon solar water heaters, the flow rate of the circulating water is conventionally calculated by equating the pressure head and the friction head.Pressure head is caused by density gradients in the loop, and the friction head is caused by friction in the plumbing arrangement.
The pressure head in the thermosiphon causes flow to occur.This flow in the collector is driven by the weight difference between the hot water column in the return pipe passing through the collector and the cold fluid column in the inlet pipe.The temperature conditions are given by the inlet temperature of the fluid   = 25 ∘ C and the inner surface temperature of the hot pipe   = 40 ∘ C; the density variations in the water along the collector are assumed to be linear for the calculations [21].The desired maximum temperature in the cooling plate should be around 50 ∘ C.
Imagine an opened thermosiphon loop as a U-tube containing a fluid with one column filled with hot fluid and the other with cold fluid.A height difference, ℎ, results due to the density differences.If instead of U-tube one has a closed loop, this ℎ leads to a driving force that produces the flow in the loop.
The continuity equations under static equilibrium in case of U-tube can be expressed by and the corresponding pressure head which is a function of the temperature and the total height of the columns.We can rewrite (2) as a function of the cold-and hot-side densities and considering a total length : To determine the thermal driving forces, it is necessary to take into account the values of   and   .Using the desired values of   = 25 ∘ C and   = 40 ∘ C, ℎ = 2.51 mm is obtained.
The friction head, flow rate, and convective coefficient are interrelated, but they also depend on several physical parameters that must be defined, such as piping type, materials, and pipe length, among others.
Using the Bernoulli equation, an energy conservation analysis can be made.For a pipe system [21] where  1 and  2 are inlet and outlet pressures,  1 and  2 are heights, and  1 and  2 are the corresponding flow velocities, the following can be written in terms of energy: where ℎ  is the added energy, ℎ  is the subtracted energy, ℎ  is the energy loss (friction head),  is the acceleration of gravity, and  is the specific weight of the fluid.On the other hand, it is necessary to include the Darcy equation for friction head ℎ  : where  is the piping length and  is its diameter that depends on the flow rate .As the friction factor  depends on the Reynolds number Re for the laminar flow [21] for a thermosiphon system, the pressure head is equalized to friction head causing the energy loss ℎ = ℎ  .
As a consequence, it is necessary to take into account the energy losses due to friction (major losses due to friction and minor losses due to changes in the size and direction of the flow path) in the loop.
The friction head can then be expressed for this case in terms of the friction factor and the flow rate The first term corresponds to the sum of losses in the inlet and outlet, where it is common to use the estimation friction coefficient   = ( eq /) = 1.5 for systems of this kind [21].
Taking in consideration the laminar flow and equalizing the thermal driving head with the friction head, we get Solving ( 9) for ℎ = 2.51 mm, a flow rate of   = 0.0436 m/s is obtained.
Once the flow rate is defined, the convective coefficient can be calculated [21].For the laminar region in the circular pipe and the temperature of 25 ∘ C, the corresponding Nusselt average number is Nu = 3.66 + 0.065 (/) Re Pr For a thermal length   ≈ 0.05 Re Pr  = 5.047 m, the condition   <  is satisfied, and the convection coefficient can be calculated as Once ℎ  and   are determined, we can take them as initial values for the design of the heat exchanger, which starts with computer simulations.

Design of the Heat Exchanger for Electric/Thermal Generating Unit
In order to determine the optimal configuration of the heat exchanger, several configurations were proposed and evaluated using commercial finite element method (FEM) software (COMSOL Multiphysics 4.2a).For the flow rate, the value obtained earlier was used:   = 0.0436 (m/s).The cooling plate temperature must not exceed 50 ∘ C, and the hot plate must be around 200 ∘ C. The solar power transformed into heat in the hot plate is around 200 W in an area of 0.08 × 0.12 m 2 , in correspondence with the 2 × 3 array of TEG elements.
The heat exchanger was designed to be as simple as possible, a flat plate attached to the commercial pipes.In Figure 3, the modeling results for several configurations are presented, changing parameters like hot plate location, pipe diameter, piping array, and welding material, among others.In this same figure, simulations from Figure 3(a) to 3(c) are for 1-inch-diameter and from Figure 3(d) to 3(f) are for 3/4inch-diameter type K copper pipe according to ASTM B-88 standard.The red areas are the hottest and the blue ones the coldest, in accordance with the reference bar in the right side of each model.
After evaluating a wide range of configurations, two options that best meet the conditions were selected and evaluated, and the results are shown in Figure 4. One-inch type K cooper pipe [21] was used in (a), obtaining a maximum value of 42.4 ∘ C in the center of the surface (red zone).3/4-inch type K pipe was chosen for case (b), leading to a maximum value of 39.4 ∘ C distributed in a more uniform way along the center of the surface.Thus, option (b) was chosen for the experiment.

Experimental Results
The actual system studied is shown in the photograph in Figure 5.The positions of the thermocouples are indicated by the red and blue points.The red points also indicate the location of the ink injection, used to give idea of the actual water flow rate.As can be seen in the image, the thermoelectric assembly is illuminated by the concentration block.
The results of the system's electrical and thermal characterization are presented in Figures 6 and 7. To estimate the system's efficiency (both electrical and thermal), the intensity of solar radiation was taken as 950 W/m, which corresponds to the direct normal irradiance (DNI) in Queretaro, Mexico, at 20 ∘ of northern latitude, at the equinox time of the year.First, the electric power generated by the system during daytime is shown in Figure 6.The measurements were taken with a matched load, so the data shown gives the maximum power available.One can see that the average power was 20 W, thus producing 120 Wh of electric energy between 10 am and 4 pm (the total energy obtained during the day was 175 Wh).These results correspond to a maximum electric efficiency of the system of 5%, which agrees well with the estimations made in [14] and with the results of modeling [16].
Figure 7 presents the thermal characteristics of the hybrid system.The average hot water tank temperature was approximately 45 ∘ C, which is sufficient for domestic applications.The variations in the thermic efficiency,  therm = ( therm / sol ) × 100, observed during the time of the experiment give in average of 50%, which is higher than that an traditional PV/thermal systems.The corresponding thermal power is 200 W, giving the 1.2 kWh of energy in the 6 h between 10 am and 4 pm.It can also be seen that the flow rate correlates with the thermal efficiency-higher flow is accompanied by higher efficiency-although the water temperature in the tank is lower.
The cost-efficiency estimation made in our previous publication for the hybrid system studied [14] showed that at industrial production (in quantities of hundreds of systems), the cost of the electric energy generation could be 3-4 US dollars per watt peak, which is almost the same as the cost of the energy generated by PV panels.This is in contrast with the typical case where the cost of energy production in hybrid systems is usually 50% higher than that in the individual devices.The cost of thermal energy in our system with TEGs is lower than that in the traditional PV/thermal systems because of the larger thermal power.

Conclusions
Performance of the designed hybrid system in the conditions at Queretaro, Mexico at the equinox time of the year has revealed that a systems electrical efficiency of 5% and thermal efficiency of 50% with the estimated cost of the electric energy production are practically equal to those of the traditional photovoltaic/thermal systems.Thus, we conclude that the solar hybrid system with the concentrator and the thermoelectric generator, even with the existing components, can be considered as a reasonable alternative to the traditional electric/thermal solar hybrid system.Taking into account the rapid progress in the development of new nanostructured and highly efficient thermoelectric materials, we can expect that in the near future performance of the TEG-based systems can surpass that of the traditional solar hybrid systems, in particular, in the solar-rich regions having relatively low latitude.

Figure 3 :
Figure 3: Computer simulations for different configurations (see text); the colors show the temperature distribution.

𝐴
: Effective flow area of the piping, m   : Specific heat capacity, J/Kg⋅ ∘ C : Piping diameter, m ℎ: Pressure head (thermal driving head), m:Friction factor :Gravity, m/s ℎ:Height,m ℎ  :Energy added to the fluid, J ℎ  :Cold fluid column height, m ℎ  :Average convection coefficient, W/m⋅ ∘ C ℎ ℎ :Hot fluid column height, m ℎ  :Friction head inside the piping, m ℎ  : Energyremovedtothefluid,J   :Friction coefficient :Thermalconductivity ,W/mK :Piping length, m  eq : Equivalent pipig length of the minor losses, m   :Thermic inlet length, m ṁ:Mas flow, Kg/s Nu:Average Nusselt number  1 :Inlet pressure, N/m  2 :Outlet pressure, N/m Pr:Prandtl number Q:Heatflux,W  sun : Solar heat input W  out : Bottomingheat transfer W  therm : Heat transfer to running water W Re: Reynolds number   : Outlet temperature in the fluid of the heating pipe, ∘ C   : Inlet temperature in the fluid of the heating pipe, ∘ C   : Inner surface temperature of the hot pipe, ∘ C : Flow velocity inside the piping, m/s  1 : Inletvelocityofthefluid,m/s  2 : Outlet velocity of the fluid, m/s  1 : Height at the inlet point, m  2 : Height at the outlet point, m  therm : Thermic efficiency : Fluid density, kg/m   : Cold fluid density, kg/m  ℎ : Hot fluid density, kg/m : Specific weight of the fluid, N/m : Dynamic viscosity, Kg/m⋅s.

Figure 6 :
Figure 6: Electric power generation as a function of time of day (a) and as a function of the temperature difference between the TEG plates (b).

Figure 7 :
Figure 7: Hot water tank temperature (a) and calculated thermal efficiency of the system and the flow rate in the thermosiphon (b).