Optimum Thermal Concentration of Solar Thermoelectric Generators (STEG) in Realistic Meteorological Condition

Global warming and air pollution concerns make renewable energies inevitable. Thermoelectric (TE) generators—solid-state devices which can convert thermal energy into electricity—are one of the candidates to capture the energy of the sun’s rays. Impact of high thermal on flat panel Solar Thermoelectric Generator (STEG) performance is known. In this research, a method to optimize thermal concentration in realistic terrestrial condition is introduced. To this end, a Simulink model of the STEG was developed, and module performance curves are determined. According to the results, Thermal concentration in realistic condition is more than double, compared to standard condition. The efficiency of the STEG was 4.5%, 6.8%, and 7.7% when the module figure of merit (ZT) was set to 0.8, 1.2, and 1.5, respectively, in locations with a higher ratio of diffused radiation (e.g., Aalborg and Denmark). These values corresponded to 70%, 106%, and 121% of the electrical power produced by parabolic troughs under the same meteorological condition. Furthermore, the possibility of controlling the ratio of heat and electricity in the cogeneration system is possible by controlling the heating flow or electric current. Heating flow can be controlled by the electrical current in STEG over 17 percent of its value in optimum condition.


Introduction
The thermoelectric effect has been known since 1821 when Thomas Seebeck observed the electric potential by applying a temperature gradient. However, low efficiency restricted their applications to thermocouple sensors and niche places, such as space exploration [1]. The first Solar Thermoelectric Generator (STEG) with the efficiency of 0.63% was introduced by Telkes in 1954 [2], who predicted that solar absorbers could increase efficiency to 1.05%. Promising theoretical [3] and experimental [4] studies in the 1990s revealed thermoelectric (TE) properties could significantly enhanced, paving the way for further studies. On the other hand, green energy has become increasingly important during the past few decades, mostly due to global warming and air pollution concerns [5][6][7]. There are three major types of solar plants, including photovoltaic, flat, and concentrated thermal absorbers.  [25].
Based on the above argument, this work aims to develop a method to optimize the thermal concentration of STEG with respect to realistic terrestrial conditions. To end this goal, a 1D model of STEG was developed in MATLAB Simulink where TEG performance was computed over a wide range of solar radiation and cold side temperature. The model simulates STEG behavior, and it can control current, voltage, and cold side temperature. To optimize TC over a year, the optimum values of each day are averaged. This is carried out in a loop to check the peak hot side temperature and adjust TC to protect module.

Solar TEG Model
Electrical equivalent model for thermal studies, which is based on the analogy between thermal and electrical parameters (Table 1), was used to simulate thermoelectric generators [26][27][28][29]. In the present study, MATLAB Simulink Simscape library was exploited to develop the model. The model includes the solar absorber, TEG, cold side heat exchanger, and electrical load. In Figure 2, the schematic of the two conceptual designs of the STEGs, which use a selective absorber, as well as vacuum insulated chamber and transparent glass are presented.   [25].
Based on the above argument, this work aims to develop a method to optimize the thermal concentration of STEG with respect to realistic terrestrial conditions. To end this goal, a 1D model of STEG was developed in MATLAB Simulink where TEG performance was computed over a wide range of solar radiation and cold side temperature. The model simulates STEG behavior, and it can control current, voltage, and cold side temperature. To optimize TC over a year, the optimum values of each day are averaged. This is carried out in a loop to check the peak hot side temperature and adjust TC to protect module.

Solar TEG Model
Electrical equivalent model for thermal studies, which is based on the analogy between thermal and electrical parameters (Table 1), was used to simulate thermoelectric generators [26][27][28][29]. In the present study, MATLAB Simulink Simscape library was exploited to develop the model. The model includes the solar absorber, TEG, cold side heat exchanger, and electrical load. In Figure 2, the schematic of the two conceptual designs of the STEGs, which use a selective absorber, as well as vacuum insulated chamber and transparent glass are presented.  The energy flow is shown in Figure 2a. The sun irradiance ( ) is absorbed by the selective absorber, passes through the collector plate, and then flows through the thermoelectric legs ( Figure  2b). The thermal concentration, the ratio between the collector ( ), and TEG ( ) surface areas are above 100, which causes high temperature and radiative losses , . Vacuum insolation can minimize conductive losses( , ). Moreover, the thermal conductivity in a typical vacuumed solar collector is in the order of 10 W mK ⁄ . In Figure 2b, the dissipated heat ( ) rejects to ambient.
Meanwhile, the rejected heat is an input of the heat for the CHP system, presented in Figure 2c. Energy balance equations of the system are, as follows: In STEGs, efficiency of the system is the product of thermal efficiency and electrical efficiency where: Thermal efficiency (3) depends on absorber properties, insulation quality, and temperature difference, while electrical efficiency Equation (4) is the function of TEG properties and the temperature difference across the TEG. The temperature must be optimized due to the inverse effect of the temperature difference on the thermal and electrical efficiencies. A STEG system includes a solar absorber, a TEG device, cold side, and an electrical load, as shown in Figure 3. Each part of the The energy flow is shown in Figure 2a. The sun irradiance ( . Q Solar ) is absorbed by the selective absorber, passes through the collector plate, and then flows through the thermoelectric legs ( Figure 2b). The thermal concentration, the ratio between the collector (A Coll ), and TEG (A TEG ) surface areas are above 100, which causes high temperature and radiative losses ( . Q L, Rad ). Vacuum insolation can minimize conductive losses ( . Q L, Cond ). Moreover, the thermal conductivity in a typical vacuumed solar collector is in the order of 10 −3 W/mK. In Figure 2b, the dissipated heat ( . Q C ) rejects to ambient. Meanwhile, the rejected heat is an input of the heat for the CHP system, presented in Figure 2c. Energy balance equations of the system are, as follows: .
In STEGs, efficiency of the system is the product of thermal efficiency and electrical efficiency where: Thermal efficiency (3) depends on absorber properties, insulation quality, and temperature difference, while electrical efficiency Equation (4) is the function of TEG properties and the temperature difference across the TEG. The temperature must be optimized due to the inverse effect of the temperature difference on the thermal and electrical efficiencies. A STEG system includes a solar absorber, a TEG device, cold side, and an electrical load, as shown in Figure 3. Each part of the system is masked to sweep across variables by assigning the thermal concentration and cold side temperature from the MATLAB workspace. system is masked to sweep across variables by assigning the thermal concentration and cold side temperature from the MATLAB workspace.

Solar Absorber Model
The input heat to solar absorber is the product of the collector surface and radiation intensity. The input energy and the absorbed heat are expressed as, = . .
The thermal losses can be expressed as: In the Simulink view, the solar absorber acts as a current source, controlled by absorber temperature, radiation intensity, and the system dynamic. The cold side was considered as a controllable voltage (temperature) source, and the electrical load block as a controllable resistor. In the case of maximum electrical efficiency, the load's resistance should be equal to the TEG's internal resistance.

TEG Model
Peltier's, Seebeck's, and Ohm's laws are governing equations for TEG devices. By dividing TEG into stacked blocks, temperature-dependent material properties can be used. In this approach, the TEG is divided into several thin TEGs, which are thermally and electrically connected in series. The output of each block is an input for the next block. Due to only small temperature differences along each block, it is possible to assume that temperature is constant and use locally constant temperaturedependent materials' properties. Figure 4 shows a schematic of this method.

Solar Absorber Model
The input heat to solar absorber is the product of the collector surface and radiation intensity. The input energy and the absorbed heat are expressed as, . .
The thermal losses can be expressed as: . .
In the Simulink view, the solar absorber acts as a current source, controlled by absorber temperature, radiation intensity, and the system dynamic. The cold side was considered as a controllable voltage (temperature) source, and the electrical load block as a controllable resistor. In the case of maximum electrical efficiency, the load's resistance should be equal to the TEG's internal resistance.

TEG Model
Peltier's, Seebeck's, and Ohm's laws are governing equations for TEG devices. By dividing TEG into stacked blocks, temperature-dependent material properties can be used. In this approach, the TEG is divided into several thin TEGs, which are thermally and electrically connected in series. The output of each block is an input for the next block. Due to only small temperature differences along each block, it is possible to assume that temperature is constant and use locally constant temperature-dependent materials' properties. Figure 4 shows a schematic of this method. The heat transfer and electrical equations in each block are, Properties of the blocks are computed based on average temperature of the block, number of the blocks and temperature dependent properties of the module can be expressed as: TEG properties for different ZT values are presented in Table 2. In Case #1 (ZT = 0.8), the coefficients were extracted by fitting curves for experimental data. Afterward, the coefficients were modified to make maximum ZT equal to 1.2 and 1.5 for Cases #2 and #3, respectively. In the present study, the value of the figure of merit refers to an effective ZT module [30]. Table 2. TEG temperature dependent properties in one verified (ZT = 0.8) and two (Z = 1.2 and ZT =1.5) modified conditions.  Each block has a Seebeck coefficient, and an electrical and thermal resistance corresponding to its average temperature.

Module Figure of Merit
The heat transfer and electrical equations in each block are, Properties of the blocks are computed based on average temperature of the block, number of the blocks and temperature dependent properties of the module can be expressed as: TEG properties for different ZT values are presented in Table 2. In Case #1 (ZT = 0.8), the coefficients were extracted by fitting curves for experimental data. Afterward, the coefficients were modified to make maximum ZT equal to 1.2 and 1.5 for Cases #2 and #3, respectively. In the present study, the value of the figure of merit refers to an effective ZT module [30]. Table 2. TEG temperature dependent properties in one verified (ZT = 0.8) and two (Z = 1.2 and ZT = 1.5) modified conditions.

Meteorological Data
The Joint Research Centre, which is the commission's science and knowledge service, provides an online database of environmental data related to renewable energies [25] exploited in this study. Two places in Europe with different insolation and temperature regimes were selected: Aalborg (57 • 2 11 N, 9 • 55 21 E) in the North of Denmark, which is mostly cloudy with diffused isolation and 7-18 daylight hours. The second location is Gibraltar (36 • 19 36 N, 5 • 37 29 E), which has more direct radiation and 10-14 daylight hours. This selection was based on a diffused to direct ratio of sun isolation (see Figure 1) in order to compare the STEG and concentrated CHP. The values of direct, diffused, and total irradiance in Aalborg and Gibraltar are shown in Table 3, according to which the diffused irradiance in Aalborg was 19% higher compared to Gibraltar. However, global irradiance was 21% lower.

Optimization Process
The optimization process was mainly carried out to find the desired value of TC so that the maximum STEG performance could be obtained in a specific location in a year. The inputs are meteorological data, STEG performances curves (at different cold side temperatures and incident radiation), and maximum hot side temperature of the TEG, which is considered a limiting factor during the process. The pairs of maximum efficiency and correlated TC were calculated based on air temperature, solar insolation, and STEG curves. The weighted mean of TC and η were the yearly optimal TC and efficiency of the system. The hot side temperature limitation was applied in the final step. TC increased discriminately to ensure that temperature is in the safe zone for the module. The maximum tolerable temperature is based on solder durability, For instance, it is considered between 400 and 450 K.

Experimental Setup
To obtain the module properties and validate the model, a sample module was tested by Tegeta at the Aalborg University [31]. This setup was designed initially to characterize high temperature oxide materials [31], and it was modified for the mid temperature module studied in this paper. Figure 5 shows the modified setup in which the ceramic heater is replaced by a ULTRAMIC 600 supplied with Keysight B2900 to maintain the hot side at the desired temperature. T-type thermocouples, with a diameter of 50 micrometers and a nominal response time of 20 milliseconds sense temperatures. The Keysight 34972 logged data every 10 s. The second channel of the Keysight B2900 was used to run constant reverse voltage and emulate the electrical load. Silicon thermal past is used between layers to ensure low thermal resistance between the hot side and TE module as well as the module and heat sink. The measurement was run over different temperatures, and all data was stored in a text file then imported to a Microsoft Excel spreadsheet for further analysis. The standard procedure to characterize TE modules was: set constant temperature and electrical load until thermal and electrical stability was achieved then log the corresponding data for each electrical load. The data was processed later to fit the output power versus electrical load. The corresponded electrical load of maximum output power is the internal electrical resistance of module. was processed later to fit the output power versus electrical load. The corresponded electrical load of maximum output power is the internal electrical resistance of module.

Experimental Validation of Simulink Model
Zmeter is a known method to characterize TE materials and modules; details can be found in previous papers [32,33]. The module properties were extracted by a series of experiments at different temperatures. Later, the obtained values of temperature, electrical voltage, current, and heat flow were used to fit curves, presented in Table 2. The nominal figure of merit for the tested module is 0.8; higher values can be produced by adjusting the data. Afterwards, the setup was used to validate the model. Predicted and measured values for the maximum power and corresponding electrical properties are presented in Table 4. The fill factor of the module was 0.7 and the number of thermocouples was 100. This module is designed to harvest energy in sensory systems.  Figure 6 shows the measured and calculated power module at two temperature differences between hot and cold sides. The mismatching between calculated and measured electrical resistance is due to sidewall thermal losses, and it is increased by increasing the average temperature, comparing between Figure 6a,b. However, the difference is acceptably less than 10 percent of the total power [34]. This proves that the model can represent a real system with an acceptable error.

Experimental Validation of Simulink Model
Zmeter is a known method to characterize TE materials and modules; details can be found in previous papers [32,33]. The module properties were extracted by a series of experiments at different temperatures. Later, the obtained values of temperature, electrical voltage, current, and heat flow were used to fit curves, presented in Table 2. The nominal figure of merit for the tested module is 0.8; higher values can be produced by adjusting the data. Afterwards, the setup was used to validate the model. Predicted and measured values for the maximum power and corresponding electrical properties are presented in Table 4. The fill factor of the module was 0.7 and the number of thermocouples was 100. This module is designed to harvest energy in sensory systems.  Figure 6 shows the measured and calculated power module at two temperature differences between hot and cold sides. The mismatching between calculated and measured electrical resistance is due to sidewall thermal losses, and it is increased by increasing the average temperature, comparing between Figure 6a,b. However, the difference is acceptably less than 10 percent of the total power [34]. This proves that the model can represent a real system with an acceptable error.

Performance Curves of TE Module
The temperature distribution along the heat conduction direction in the TEG module is shown in Figure 7. Due to different heat transfer mechanisms, the temperature distribution along the TE elements is not linear. In order to have a faster simulation, the thermal mass of the solar collector was eliminated from the simulations. It is worth noting that the long response time of the TEGs is an advantage that causes higher stability in power generation systems, thereby diminishing the necessity of fast MPPT [35] algorithms and providing an intrinsic storage capability for the system. According to Tables 2 and 5, the STEG time constant mostly depends on the solar thermal collector. For instance, a 1 ℃ change in 500 μm thick copper takes 17 s when radiation is 1000 W m ⁄ . Table 5. Thermal properties of thermoelectric materials and solar collector.

Performance Curves of TE Module
The temperature distribution along the heat conduction direction in the TEG module is shown in Figure 7. Due to different heat transfer mechanisms, the temperature distribution along the TE elements is not linear. In order to have a faster simulation, the thermal mass of the solar collector was eliminated from the simulations. It is worth noting that the long response time of the TEGs is an advantage that causes higher stability in power generation systems, thereby diminishing the necessity of fast MPPT [35] algorithms and providing an intrinsic storage capability for the system. According to Tables 2 and 5, the STEG time constant mostly depends on the solar thermal collector. For instance, a 1 change in 500 µm thick copper takes 17 s when radiation is 1000 W/m 2 .

Performance Curves of TE Module
The temperature distribution along the heat conduction direction in the TEG module is shown in Figure 7. Due to different heat transfer mechanisms, the temperature distribution along the TE elements is not linear. In order to have a faster simulation, the thermal mass of the solar collector was eliminated from the simulations. It is worth noting that the long response time of the TEGs is an advantage that causes higher stability in power generation systems, thereby diminishing the necessity of fast MPPT [35] algorithms and providing an intrinsic storage capability for the system. According to Tables 2 and 5, the STEG time constant mostly depends on the solar thermal collector. For instance, a 1 ℃ change in 500 μm thick copper takes 17 s when radiation is 1000 W m ⁄ .   Higher radiation intensity at constant cold side temperature increases the efficiency and decreases the corresponding TC, whereas an increase in the cold side temperature leads to a reduction in both efficiency and the corresponding TC, which is not a desirable condition. The opposite effects of air temperature and solar radiation intensity cause the relative efficiency in winter and summer days. For example, two different conditions in summer (600 W/m 2 and 300 K) and winter (400 W/m 2 and 280 K), have an equal efficiency of 4.8% for ZT = 0.8 and different values of TCs. Therefore, it is necessary to perform long-term optimization instead of only STC condition optimization. While the module with higher ZT increases efficiency, it simultaneously decreases TC. Incensement in ZT is due to the lower thermal conductivity of TE materials, which may lead to the manifestation of optimum temperature gradient at lower TC. For example, increasing ZT from 0.8 to 1.2 leads to the decrease of TCs at 200 W/m 2 and 500 W/m 2 by order of 200 and 110, respectively.
Energies 2018, 11, x FOR PEER REVIEW 10 of 16 Figure 8 exhibits the maximum STEG efficiency (a-c) and optimum TC (d-f) contours at different terrestrial conditions for various ZT values. Higher radiation intensity at constant cold side temperature increases the efficiency and decreases the corresponding TC, whereas an increase in the cold side temperature leads to a reduction in both efficiency and the corresponding TC, which is not a desirable condition. The opposite effects of air temperature and solar radiation intensity cause the relative efficiency in winter and summer days. For example, two different conditions in summer (600 W m and 300 K ⁄ ) and winter (400 W m and 280 K ⁄ ), have an equal efficiency of 4.8% for ZT = 0.8 and different values of TCs. Therefore, it is necessary to perform long-term optimization instead of only STC condition optimization. While the module with higher ZT increases efficiency, it simultaneously decreases TC. Incensement in ZT is due to the lower thermal conductivity of TE materials, which may lead to the manifestation of optimum temperature gradient at lower TC. For example, increasing ZT from 0.8 to 1.  (d-f) Optimum thermal concentration for different TEG properties according to cold side temperature and solar radiation intensity. In the case where ZT, cold side temperature, and radiation are 0.8, 300 K, and 300 W m ⁄ , respectively, the optimum TC is 750 and the maximum efficiency is 3.8%. Figures 9 and 10 show how the optimum efficiency and TC can vary according to the terrestrial parameters, which confirms the necessity of yearly design optimization. Each blue dot in the figures refers to one data point. As observed, points with high values of TC were in early mornings or late evenings. Table 6 shows the daily energy production expectation for different values of ZT in comparison to the parabolic trough, which is a mature CHP technology with an average efficiency of 12.5% that can only capture DNI. Due to the higher amount of diffused radiation in Aalborg (49%), the increased ZT had more of an impact on energy production, compared to Gibraltar where diffused radiation is 33%. The parabolic trough can capture 67%of sun rays. STEG with ZT = 0.8 can generate 70% of the expected electrical power of the parabolic trough, and it could increase to 106% and 121% with ZT of 1.2 and 1.5, respectively. (d-f) Optimum thermal concentration for different TEG properties according to cold side temperature and solar radiation intensity. In the case where ZT, cold side temperature, and radiation are 0.8, 300 K, and 300 W/m 2 , respectively, the optimum TC is 750 and the maximum efficiency is 3.8%. Figures 9 and 10 show how the optimum efficiency and TC can vary according to the terrestrial parameters, which confirms the necessity of yearly design optimization. Each blue dot in the figures refers to one data point. As observed, points with high values of TC were in early mornings or late evenings. Table 6 shows the daily energy production expectation for different values of ZT in comparison to the parabolic trough, which is a mature CHP technology with an average efficiency of 12.5% that can only capture DNI. Due to the higher amount of diffused radiation in Aalborg (49%), the increased ZT had more of an impact on energy production, compared to Gibraltar where diffused radiation is 33%. The parabolic trough can capture 67%of sun rays. STEG with ZT = 0.8 can generate 70% of the expected electrical power of the parabolic trough, and it could increase to 106% and 121% with ZT of 1.2 and 1.5, respectively.    In the second scenario, the dissipated heat on the cold side was the input to the heating system. This might cause a higher cold side temperature and lower total efficiency. Figure 11 shows the maximum efficiency and the corresponding TC for different cold side temperatures. The cold side temperature is assumed to be constant throughout a year, which is possible by control of the liquid flow and TEG current. The desirable water temperature of flat solar heat collectors was 320-350 K, which obtained efficiency between 3.5% and 2.5% for ZT = 0.8 in Aalborg. By increasing the ZT to 1.2, the efficiency increased to 6 and 4.5, respectively. The STEG efficiency in Gibraltar in the coproduction scenario was relatively higher compared to Aalborg. At lower water temperatures, this difference became more significant. By considering constant temperature for both locations, Aalborg lost the lower ambient air temperature merit. Compared to similar studies [19,36], where optimum TC varies between 120 and 270, the proposed model shows the significantly different values for TC. In the original study on STEG [19] it is between 180 and 220. Studying different materials by 3D simulations [36] shows the similar range for optimum TC. This is mainly due to the constant STC condition. The seasonal effect was studied previously [36] for two scenarios, winter and summer, and it is pointed out that the terrestrial condition significantly degrades STEG performance. Kossyvakis et al. [36] proposed to change the Fresnel lens or the focal point. A developed model represents an effort to overcome the seasonal effect by applying realistic data in the design step. In the second scenario, the dissipated heat on the cold side was the input to the heating system. This might cause a higher cold side temperature and lower total efficiency. Figure 11 shows the maximum efficiency and the corresponding TC for different cold side temperatures. The cold side temperature is assumed to be constant throughout a year, which is possible by control of the liquid flow and TEG current. The desirable water temperature of flat solar heat collectors was 320-350 K, which obtained efficiency between 3.5% and 2.5% for ZT = 0.8 in Aalborg. By increasing the ZT to 1.2, the efficiency increased to 6 and 4.5, respectively. The STEG efficiency in Gibraltar in the coproduction scenario was relatively higher compared to Aalborg. At lower water temperatures, this difference became more significant. By considering constant temperature for both locations, Aalborg lost the lower ambient air temperature merit. Compared to similar studies [19,36], where optimum TC varies between 120 and 270, the proposed model shows the significantly different values for TC. In the original study on STEG [19] it is between 180 and 220. Studying different materials by 3D simulations [36] shows the similar range for optimum TC. This is mainly due to the constant STC condition. The seasonal effect was studied previously [36] for two scenarios, winter and summer, and it is pointed out that the terrestrial condition significantly degrades STEG performance. Kossyvakis et al. [36] proposed to change the Fresnel lens or the focal point. A developed model represents an effort to overcome the seasonal effect by applying realistic data in the design step.

Controlling by Electric Current
According to the results, higher TC had better efficiency in the early morning and late noon due to higher energy flux. Meanwhile, lower TC provided more efficiency at noon because of lower heat losses. The results showed that the optimum point of TC was between these two values. There are two parameters with the ability to regulate the electricity generator and rejected heat at cold side: (i) electric current, which can be referred to as electrical load resistance, and (ii) cold side thermal resistance, which can be recognized as a volumetric flow of heat exchanger fluid. Due to the Peltier effect, the electric current variation exerts an impact on the thermal conduction. Figure 12 exhibits the variation of the load resistance from 0 (short circuit) to 80 ohms (open circuit), causing a 17.5% change in dissipated heat at the cold side of TEG. It means that when the solar radiation variation is below 17.5%, the water temperature can be regulated. Reversely, the output power can be kept constant by controlling the water flow of the heat exchanger.

Controlling by Electric Current
According to the results, higher TC had better efficiency in the early morning and late noon due to higher energy flux. Meanwhile, lower TC provided more efficiency at noon because of lower heat losses. The results showed that the optimum point of TC was between these two values. There are two parameters with the ability to regulate the electricity generator and rejected heat at cold side: (i) electric current, which can be referred to as electrical load resistance, and (ii) cold side thermal resistance, which can be recognized as a volumetric flow of heat exchanger fluid. Due to the Peltier effect, the electric current variation exerts an impact on the thermal conduction. Figure 12 exhibits the variation of the load resistance from 0 (short circuit) to 80 ohms (open circuit), causing a 17.5% change in dissipated heat at the cold side of TEG. It means that when the solar radiation variation is below 17.5%, the water temperature can be regulated. Reversely, the output power can be kept constant by controlling the water flow of the heat exchanger. resistance, which can be recognized as a volumetric flow of heat exchanger fluid. Due to the Peltier effect, the electric current variation exerts an impact on the thermal conduction. Figure 12 exhibits the variation of the load resistance from 0 (short circuit) to 80 ohms (open circuit), causing a 17.5% change in dissipated heat at the cold side of TEG. It means that when the solar radiation variation is below 17.5%, the water temperature can be regulated. Reversely, the output power can be kept constant by controlling the water flow of the heat exchanger.

Conclusions
There are several challenges in designing of thermoelectric generator systems mainly due to the several heat and electricity conversion mechanisms. The design of thermoelectric generator systems gets more complicated when dealing with an optimal design for terrestrial applications, which imposes more dynamic parameters to the system. STEGs have been designed and tested under standard test conditions based on the thermal concentration (TC) concept, showing the potential to reach energy conversion efficiency from sun rays to electricity above 4%. In this study, the new Simulink model to simulate STEGs under different conditions was presented and validated by experimental tests. The optimization results confirmed the necessity of long-term optimization. The optimum TC and maximum efficiency of STEG at two different realistic locations of Aalborg and Gibraltar were determined and the results showed that, in locations with a higher ratio of diffuse to

Conclusions
There are several challenges in designing of thermoelectric generator systems mainly due to the several heat and electricity conversion mechanisms. The design of thermoelectric generator systems gets more complicated when dealing with an optimal design for terrestrial applications, which imposes more dynamic parameters to the system. STEGs have been designed and tested under standard test conditions based on the thermal concentration (TC) concept, showing the potential to reach energy conversion efficiency from sun rays to electricity above 4%. In this study, the new Simulink model to simulate STEGs under different conditions was presented and validated by experimental tests. The optimization results confirmed the necessity of long-term optimization. The optimum TC and maximum efficiency of STEG at two different realistic locations of Aalborg and Gibraltar were determined and the results showed that, in locations with a higher ratio of diffuse to direct radiation, STEGs could generate up to 70% of the energy of parabolic trough systems. This can increase by 106 to 121% using recently developed thermoelectric materials or cascade TEGs. The results demonstrated that the coupling between electrical and thermal flows provided the opportunity to control the temperature of cooling fluid and electrical power production. The flat-panel STEGs showed a greater number of advantages in locations with a higher ratio of diffuse to direct radiation. Considering promising new materials, which had higher ZTs and lower prices, STEGs could be more efficient even in single generation systems. In addition to the development of the materials, there are engineering challenges in this field, including the geometry of STEGs to provide uniform temperature at the hot and cold sides of the TEG, manufacturing challenges to provide maximum insolation between the hot and cold sides of the TEG, and the optimal STEG dimensions to achieve the greatest economic benefit.