Next Article in Journal
Multi-Objective Optimization Design of an Electrohydrostatic Actuator Based on a Particle Swarm Optimization Algorithm and an Analytic Hierarchy Process
Previous Article in Journal
A Stochastic Approach to Energy Policy and Management: A Case Study of the Pakistan Energy Crisis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Meysam Karami Rad
1,2,*,
Mahmoud Omid
1,*,
Ali Rajabipour
1,
Fariba Tajabadi
3,
Lasse Aistrup Rosendahl
4 and
Alireza Rezaniakolaei
4
1
Department of Agricultural Machinery Engineering, Faculty of Agricultural Engineering and Technology, University of Tehran, Karaj 31787, Iran
2
Laboratory of Organic Electronics, Department of Science and Technology, Linköping University, 60174 Norrköping, Sweden
3
Materials and Energy Research Center, Imam Khomeini Blvd, Meshkin-Dasht, Karaj 31787, Iran
4
Department of Energy Technology, Aalborg University, Pontoppidanstraede 111, 9220 Aalborg, Denmark
*
Authors to whom correspondence should be addressed.
Energies 2018, 11(9), 2425; https://doi.org/10.3390/en11092425
Submission received: 24 August 2018 / Revised: 6 September 2018 / Accepted: 10 September 2018 / Published: 13 September 2018
(This article belongs to the Section A: Sustainable Energy)

Abstract

:
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.

1. 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. Photovoltaic can directly convert photons energy to electricity, where based on semiconductor band gap theory efficiency is limited to 32.8% for single-junction GaAs cell. However, practical obstacles make the real efficiency a fraction of the theoretical one [8]. Solar thermal collectors absorb the energy of the sun’s rays, which is typically used to provide warm water [9]. Having a large emissivity in the solar wavelengths and a small emissivity in the long wavelengths, selective solar absorbers can convert more than 95% of solar energy into heat [10], which leads to efficiency within the range of 40 to 60% with working fluid temperature below 100 [11]. Parabolic troughs were present in 96.3% of installed concentrated solar power plants (CSPs) until 2011, with the peak efficiency of 21%. Meanwhile, the electricity efficiency for achieved system level of sun rays is in the range of 10 to 15% with 18% mid-term perspective [12,13]. A major drawback is that CSPs only capture DNI rays and their performance is degraded under cloudy conditions [13].
TEGs require a high temperature difference to show good performance [2,14] while the solar flux is constant. Different types of STEGs, hybrid photovoltaic–thermoelectric generators, and cogeneration of heat and electricity by concentrated photovoltaic systems were reviewed in previous studies [15,16,17,18,19,20,21,22,23,24,25,26,27,28]. Moreover, concentrated optical systems [15,16] and flat-panel STEGs, which use the thermal concentration (TC) [2] to conduct heat in the absorber laterally to concentrate heat where the TEG is located, have been evaluated before.
In 2017, Olsen et al. [15] designed a high-temperature STEG with segmented n-type leg consisting of skutterudite and L a 3 T e 4 , and p-type leg including skutterudite and Y b 14 M n 21 S b 11 in the designed thermal/optical cavity, so that the whole complement of light wavelengths could be captured and an efficiency of 15% could be achieved. In the aforementioned research, the efficiency of 18% was predicted by higher optical concentration (200–300 Suns and T h = 1000   ) and active cooling on the cold side of the TEG. This shows a 2-fold gain, compared to experimental results [16] in 2016 with efficiency of 7.4% (TC = 211, selective absorber, vacuum insulators, segmented legs, and T h = 600   ).
Among several researchers who applied TC concepts [17,18,19,20,21], Kraemer et al. [19] carried out an interesting study, reporting the highest STEG efficiency of 5% without any optical concentrator. They designed a flat-panel STEG with a selective absorber and nanostructured thermoelectric materials. According to the aforementioned research, efficiency was in the range of 4.3 to 5.2% at AM1.5 while the ZT material was equal to 1.04 at an optimum temperature of 100   [15,22]. Kraemer et al. also predicted an efficiency of 14% by a more efficient material with ZT = 2 and an optical concentration of 10. Clearly, the concentrated STEGs show better performance compared to flat panels. Nonetheless, STEGs have some disadvantages, including the complexity of the system and the fact they can only capture DNI [16]. Generally, concentrating systems are only economical for locations with DNI above 1800   kWkm 2 year 1 [13], this restricts the application of concentrator systems in low-DNI locations.
The theoretical studies [23] on flat-panel STEGs include several simplifications, such as temperature independent material properties, geometry-independent radiation heat losses, and neglecting thermal/electrical contact resistances. A more detailed model for STEGs was proposed by Kraemer et al. [24], who considered all of the aforementioned effects and compared the model with previous studies in order to highlight the importance of each assumption. Furthermore, some modified simplifications were introduced to maintain the accuracy of the model and reduce the computational time. By considering a geometric ratio ( L / A T E G ) to determine thermal resistance of the STEG for cold side temperatures between 25 and 100 °C, the optimum geometric ratio was determined to be between 0.32 and 0.41 m. These researchers also calculated the output energy for the period of one day, concluding that the results were not exactly sensitive to the geometric ratio. However, their model was unable to generalize this conclusion throughout the year due to the variation of the temperature and radiation.
Figure 1 shows solar irradiation of two locations with different solar incident specifications, including Aalborg in the North of Denmark with high diffused irradiance and low global irradiance and Gibraltar city in Gibraltar with lower diffused irradiance and higher global irradiance. In contrast to Gibraltar, annual direct irradiance in Aalborg is below the minimum DNI required for concentrated systems ( 1800   kWkm 2 year 1 ) . With respect to the actual terrestrial conditions (e.g., air temperature and irradiance in a period of one year), and high dependency of STEG performance on thermal concentration value and cold side temperature, the optimization objective of the STEG must be achieved to maximize the efficiency for a year instead of a specific standard test condition (STC) in order to avoid obtaining yearly performance below the expected rate.
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.

2. Materials and Methods

2.1. 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 ( Q ˙ S o l a r ) 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 C o l l ) , and TEG ( A T E G ) surface areas are above 100, which causes high temperature and radiative losses ( Q ˙ L ,   R a d ) . Vacuum insolation can minimize conductive losses ( Q ˙ L ,   C o n d ) . 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:
Q ˙ S o l a r = Q ˙ L ,   R a d + Q ˙ L , C o n d + P E + Q ˙ C
Q ˙ C = Q ˙ P e l t i e r + Q ˙ C o n d
In STEGs, efficiency of the system is the product of thermal efficiency and electrical efficiency where:
η R H = Q ˙ S o l a r Q ˙ L ,   R a d Q ˙ L , C o n d Q ˙ S o l a r
η H E = P E Q ˙ S o l a r Q ˙ L ,   R a d Q ˙ L , C o n d = P E P E + Q ˙ C
η = ( η R H ) × ( η H E )
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.

2.2. 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,
Q ˙ H = Q ˙ A b s o r b e r Q ˙ L , R a d Q ˙ L , C o n d
Q ˙ A b s o r b e = τ · α · Q ˙ S o l a r
Q ˙ S o l a r = q ˙ S o l a r · T C · A T E G
The thermal losses can be expressed as:
Q ˙ L , R a d = A c o l l ε · σ S B · [ T a b s o r b e r 4 T a m b 4 ]
Q ˙ L , c o n d = K e n c l o s e [ T a b s o r b e r T a m b ]
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.

2.3. 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,
Q ˙ h i = [ S . I + K t h i ] [ T h i T c i ] V i I
V i = S [ T h i T c i ]
I = V R e + R L o a d
V = i = 1 n V i
R e = i = 1 n R i
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:
R e i = R e ( T i ) i
K t h i = [ K t h ( T i ) ] · i
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].

2.4. 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.

2.5. 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.
T C o p t = η 1 T C 1 q s o l , 1 + η 2 T C 2 q s o l , 2 + η n T C n q s o l , n η 1 q s o l , 1 + η 2 q s o l , 2 + + η n q s o l , n
η o p t = η 1 q s o l , 1 + η 2 q s o l , 2 + η n q s o l , n q s o l , 1 + q s o l , 2 + + q s o l , n
Each point refers to 15 min of daylight. The weighted average (Equations (18) and (19)) gives the optimum values of TC and η efficiency in the period of one year.

2.6. 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.

3. Results and Discussion

3.1. 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.

3.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 Table 2 and Table 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 .
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 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.

3.3. Optimization Results

Figure 9 and Figure 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.

3.4. 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.

4. 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.

Author Contributions

Software, M.K.R.; Validation, M.K.R. and A.R.; Formal Analysis, M.K.R. and A.R.; Data Curation, M.K.R.; Writing-Original Draft Preparation, M.K.R.; Writing-Review & Editing, M.K.R.; Visualization, M.K.R.; Supervision, M.O. and F.T.; Funding Acquisition, L.A.R. and M.K.R.

Acknowledgments

The authors would like to acknowledge the research center of the University of Tehran, Iran, for financial support of this project, and the NanoCaTe project granted by the FP7 Program under Grant No. 604647 for providing the thermoelectric module and test facilities. The first author (M.K.R.) has finished the paper while he is working at Linköping University. Hence, He would kindly acknowledge Laboratory of Organic Electronics (LOE) and LiU Biblioteket for the financial support.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

κ Thermal conductivity [ W / mK ] Rectrical resistance [ Ω ]
σ Electrical conductivity [ 1 / Ω . m ] l Thickness of vacuum layer
S Seebeck coefficient [ μ V / K ] VVtage [ V ]
R t h Thermal resistance [ K W ] Electrical current [ A ]
K Thermal conductivity [ W K ] Q ˙ Heat flux [ W ]
C t h Thermal mass [ J K ] PPower [ W ]
T Temperature [ K ] n Data point index
L Length [m] τ glass transparency
ASurface Area [ m 2 ] α absorptance
η efficiency ε effective emittance
q ˙ S o l a r solar intensity [ W / m 2 ] σSBStefan Boltzmann constant
Abbreviation
FEMFinite Element MethodSTEGSolar Thermoelectric Generator
DNIDirect Normal Incident [ W / m 2 ] TEermoelectric
CPHCombined heat and power TEGThermoelectric Generator
CSPConcentrated Solar power PlantsTCthermal concentration
SolSolar
Subscribes
ambAmbientLLoss
cCold sideTEGThermoelectric Generator
collCollectorthThermal
EElectrical R H Ray to heat
hHot side H E Heat to electricity
iSection index
nNumber of section

References

  1. Rowe, D. Applications of nuclear-powered thermoelectric generators in space. Appl. Energy 1991, 40, 241–271. [Google Scholar] [CrossRef]
  2. Telkes, M. Solar thermoelectric generators. J. Appl. Phys. 1954, 25, 765–777. [Google Scholar] [CrossRef]
  3. Hicks, L.; Dresselhaus, M. Effect of quantum-well structures on the thermoelectric figure of merit. Phys. Rev. B 1993, 47, 12727. [Google Scholar] [CrossRef]
  4. Hicks, L.; Dresselhaus, M. Thermoelectric figure of merit of a one-dimensional conductor. Phys. Rev. B 1993, 47, 16631. [Google Scholar] [CrossRef]
  5. Sokka, L.; Sinkko, T.; Holma, A.; Manninen, K.; Pasanen, K.; Rantala, M.; Leskinen, P. Environmental impacts of the national renewable energy targets—A case study from finland. Renew. Sustain. Energy Rev. 2016, 59, 1599–1610. [Google Scholar] [CrossRef]
  6. Noppers, E.H.; Keizer, K.; Milovanovic, M.; Steg, L. The importance of instrumental, symbolic, and environmental attributes for the adoption of smart energy systems. Energy Policy 2016, 98, 12–18. [Google Scholar] [CrossRef] [Green Version]
  7. Inglesi-Lotz, R. The impact of renewable energy consumption to economic growth: A panel data application. Energy Econ. 2016, 53, 58–63. [Google Scholar] [CrossRef] [Green Version]
  8. Luque, A.; Hegedus, S. Handbook of Photovoltaic Science and Engineering; John Wiley & Sons: ‎Hoboken, NJ, USA, 2011. [Google Scholar]
  9. Tian, Y.; Zhao, C.-Y. A review of solar collectors and thermal energy storage in solar thermal applications. Appl. Energy 2013, 104, 538–553. [Google Scholar] [CrossRef] [Green Version]
  10. Cao, F.; McEnaney, K.; Chen, G.; Ren, Z. A review of cermet-based spectrally selective solar absorbers. Energy Environ. Sci. 2014, 7, 1615–1627. [Google Scholar] [CrossRef] [Green Version]
  11. Kumar, R.; Rosen, M.A. A critical review of photovoltaic–thermal solar collectors for air heating. Appl. Energy 2011, 88, 3603–3614. [Google Scholar] [CrossRef]
  12. Zhang, H.; Baeyens, J.; Degrève, J.; Cacères, G. Concentrated solar power plants: Review and design methodology. Renew. Sustain. Energy Rev. 2013, 22, 466–481. [Google Scholar] [CrossRef]
  13. Müller-Steinhagen, H.; Trieb, F. Concentrating solar power. Ingenia Inform. QR Acad. Eng. 2004, 18, 43–50. [Google Scholar]
  14. Chen, W.-H.; Wu, P.-H.; Lin, Y.-L. Performance optimization of thermoelectric generators designed by multi-objective genetic algorithm. Appl. Energy 2018, 209, 211–223. [Google Scholar] [CrossRef]
  15. Olsen, M.; Warren, E.; Parilla, P.; Toberer, E.; Kennedy, C.; Snyder, G.; Firdosy, S.; Nesmith, B.; Zakutayev, A.; Goodrich, A. A high-temperature, high-efficiency solar thermoelectric generator prototype. Energy Procedia 2014, 49, 1460–1469. [Google Scholar] [CrossRef]
  16. Kraemer, D.; Jie, Q.; McEnaney, K.; Cao, F.; Liu, W.; Weinstein, L.A.; Loomis, J.; Ren, Z.; Chen, G. Concentrating solar thermoelectric generators with a peak efficiency of 7.4%. Nat. Energy 2016, 1, 16153. [Google Scholar] [CrossRef] [Green Version]
  17. Fuschillo, N.; Gibson, R.; Eggleston, F.; Epstein, J. Solar thermoelectric generator for near-earth space applications. IEEE Trans. Electron Devices 1966, 13, 426–432. [Google Scholar] [CrossRef]
  18. Goldsmid, H.; Giutronich, J.; Kaila, M. Solar thermoelectric generation using bismuth telluride alloys. Sol. Energy 1980, 24, 435–440. [Google Scholar] [CrossRef]
  19. Kraemer, D.; Poudel, B.; Feng, H.-P.; Caylor, J.C.; Yu, B.; Yan, X.; Ma, Y.; Wang, X.; Wang, D.; Muto, A. High-performance flat-panel solar thermoelectric generators with high thermal concentration. Nat. Mater. 2011, 10, 532–538. [Google Scholar] [CrossRef] [PubMed]
  20. Lv, S.; He, W.; Hu, D.; Zhu, J.; Li, G.; Chen, H.; Liu, M. Study on a high-performance solar thermoelectric system for combined heat and power. Energy Convers. Manag. 2017, 143, 459–469. [Google Scholar] [CrossRef]
  21. Sun, D.; Shen, L.; Yao, Y.; Chen, H.; Jin, S.; He, H. The real-time study of solar thermoelectric generator. Appl. Therm. Eng. 2017, 119, 347–359. [Google Scholar] [CrossRef]
  22. Joshi, G.; Lee, H.; Lan, Y.; Wang, X.; Zhu, G.; Wang, D.; Gould, R.W.; Cuff, D.C.; Tang, M.Y.; Dresselhaus, M.S. Enhanced thermoelectric figure-of-merit in nanostructured p-type silicon germanium bulk alloys. Nano Lett. 2008, 8, 4670–4674. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, G. Theoretical efficiency of solar thermoelectric energy generators. J. Appl. Phys. 2011, 109, 104908. [Google Scholar] [CrossRef]
  24. Kraemer, D.; McEnaney, K.; Chiesa, M.; Chen, G. Modeling and optimization of solar thermoelectric generators for terrestrial applications. Sol. Energy 2012, 86, 1338–1350. [Google Scholar] [CrossRef] [Green Version]
  25. Photovoltaic Geographical Information System—Interactive Maps. Available online: http://re.jrc.ec.europa.eu (accessed on 14 June 2017).
  26. Mitrani, D.; Tomé, J.A.; Salazar, J.; Turó, A.; García, M.J.; Chávez, J.A. Methodology for extracting thermoelectric module parameters. IEEE Trans. Instrum. Meas. 2005, 54, 1548–1552. [Google Scholar] [CrossRef]
  27. Dalola, S.; Ferrari, M.; Ferrari, V.; Guizzetti, M.; Marioli, D.; Taroni, A. Characterization of thermoelectric modules for powering autonomous sensors. IEEE Trans. Instrum. Meas. 2009, 58, 99–107. [Google Scholar] [CrossRef]
  28. Phillip, N.; Maganga, O.; Burnham, K.J.; Dunn, J.; Rouaud, C.; Ellis, M.A.; Robinson, S. Modelling and Simulation of a Thermoelectric Generator for Waste Heat Energy Recovery in Low Carbon Vehicles. In Proceedings of the 2012 2nd International Symposium on Environment Friendly Energies and Applications, Newcastle upon Tyne, UK, 25–27 June 2012. [Google Scholar]
  29. Kane, A.; Verma, V.; Singh, B. Temperature Dependent Analysis of Thermoelectric Module Using Matlab/Simulink. In Proceedings of the 2012 IEEE International Conference on Power and Energy (PECon), Kota Kinabalu, Malaysia, 2–5 December 2012. [Google Scholar]
  30. Kim, H.S.; Liu, W.; Chen, G.; Chu, C.-W.; Ren, Z. Relationship between thermoelectric figure of merit and energy conversion efficiency. Proc. Natl. Acad. Sci. USA 2015, 112, 8205–8210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Man, E.A.; Schaltz, E.; Rosendahl, L.; Kolaei, A.R.; Platzek, D. Experimental characterization procedure of thermoelectric generator modules under non-equilibrium thermal conditions. In Proceedings of the 2014 International Conference on Thermoelectrics, Nashville, TN, USA, 6–10 July 2014. [Google Scholar]
  32. Amatya, R.; Mayer, P.; Ram, R. High temperature z-meter setup for characterizing thermoelectric material under large temperature gradient. Rev. Sci. Instrum. 2012, 83, 075117. [Google Scholar] [CrossRef] [PubMed]
  33. Ando Junior, O.; Calderon, N.; de Souza, S. Characterization of a thermoelectric generator (teg) system for waste heat recovery. Energies 2018, 11, 1555. [Google Scholar] [CrossRef]
  34. Shakouri, A. Recent developments in semiconductor thermoelectric physics and materials. Annu. Rev. Mater. Res. 2011, 41, 399–431. [Google Scholar] [CrossRef]
  35. Dalala, Z.; Saadeh, O.; Bdour, M.; Zahid, Z. A new maximum power point tracking (mppt) algorithm for thermoelectric generators with reduced voltage sensors count control. Energies 2018, 11, 1826. [Google Scholar] [CrossRef]
  36. Kossyvakis, D.; Vossou, C.; Provatidis, C.; Hristoforou, E. Computational analysis and performance optimization of a solar thermoelectric generator. Renew. Energy 2015, 81, 150–161. [Google Scholar] [CrossRef]
Figure 1. Direct normal irradiation and the ratio between diffused and global irradiation of two different places in Europe with different annual direct irradiance ( 2046   kWm 2 year 1 in Gibraltar and 1157   kWm 2 year 1 in Aalborg) [25].
Figure 1. Direct normal irradiation and the ratio between diffused and global irradiation of two different places in Europe with different annual direct irradiance ( 2046   kWm 2 year 1 in Gibraltar and 1157   kWm 2 year 1 in Aalborg) [25].
Energies 11 02425 g001
Figure 2. Schematic of the studied Solar Thermoelectric Generator (STEG) designs. (a) Energy flow diagram in the system; (b) electric power generation; (c) heat and electricity cogeneration.
Figure 2. Schematic of the studied Solar Thermoelectric Generator (STEG) designs. (a) Energy flow diagram in the system; (b) electric power generation; (c) heat and electricity cogeneration.
Energies 11 02425 g002
Figure 3. First level of STEG Simulink. Solar absorber, TEG, cold side, and load are four masked parts of a STEG system. The TC and cold side temperature can be fed from the Matlab workspace while other parameters are set inside the simulink.
Figure 3. First level of STEG Simulink. Solar absorber, TEG, cold side, and load are four masked parts of a STEG system. The TC and cold side temperature can be fed from the Matlab workspace while other parameters are set inside the simulink.
Energies 11 02425 g003
Figure 4. A schematic of the dividing thermoelectric generators (TEG) into stacked layers. Each block has a Seebeck coefficient, and an electrical and thermal resistance corresponding to its average temperature.
Figure 4. A schematic of the dividing thermoelectric generators (TEG) into stacked layers. Each block has a Seebeck coefficient, and an electrical and thermal resistance corresponding to its average temperature.
Energies 11 02425 g004
Figure 5. Experimental setup to study the transient behavior of the thermoelectric generator.
Figure 5. Experimental setup to study the transient behavior of the thermoelectric generator.
Energies 11 02425 g005
Figure 6. Experimental results and model prediction for module power output. (a) Small temperature difference and (b) high temperature difference.
Figure 6. Experimental results and model prediction for module power output. (a) Small temperature difference and (b) high temperature difference.
Energies 11 02425 g006
Figure 7. The temperature distribution along TEG legs.
Figure 7. The temperature distribution along TEG legs.
Energies 11 02425 g007
Figure 8. (ac) Efficiency contours for different TEG properties according to cold side temperature and solar radiation intensity. (df) 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%.
Figure 8. (ac) Efficiency contours for different TEG properties according to cold side temperature and solar radiation intensity. (df) 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%.
Energies 11 02425 g008
Figure 9. Maximum efficiency and optimum TC in Aalborg. Black lines are the weighted averages for each year and the blue dots are the value for each point. X axes are data points, each point represents 15 min of sunshine.
Figure 9. Maximum efficiency and optimum TC in Aalborg. Black lines are the weighted averages for each year and the blue dots are the value for each point. X axes are data points, each point represents 15 min of sunshine.
Energies 11 02425 g009
Figure 10. Maximum efficiency and optimum TC in Gibraltar. Black lines are the weighted average for a year and blue dots are the value for each point. X axes are data points, each point represents 15 min of sunshine.
Figure 10. Maximum efficiency and optimum TC in Gibraltar. Black lines are the weighted average for a year and blue dots are the value for each point. X axes are data points, each point represents 15 min of sunshine.
Energies 11 02425 g010
Figure 11. Maximum efficiency and corresponded TC for different values of ZT (ac) in the case where cold side temperature is constant.
Figure 11. Maximum efficiency and corresponded TC for different values of ZT (ac) in the case where cold side temperature is constant.
Energies 11 02425 g011
Figure 12. The effect of electrical load and thermal resistance of cold side on dissipated heat at cold side.
Figure 12. The effect of electrical load and thermal resistance of cold side on dissipated heat at cold side.
Energies 11 02425 g012
Table 1. Analogy between thermal and electrical parameters in this study.
Table 1. Analogy between thermal and electrical parameters in this study.
Thermal QuantityElectrical Equivalent
Temperature, [ K ] Voltage, [V]
Heat flow, Q ˙ = ( Δ T R t h )   [ W ] Current, I = ( Δ V R )   [ A ]
Thermal conductivity, κ   [ W mK ] Electrical conductivity, σ   [ 1 Ω m ]
Thermal mass, C t h   [ J K ] Electrical capacitance, [F]
ermal resistance, R t h   [ K W ] Electrical resistance, R   [ Ω ]
Table 2. TEG temperature dependent properties in one verified (ZT = 0.8) and two (Z = 1.2 and ZT = 1.5) modified conditions.
Table 2. TEG temperature dependent properties in one verified (ZT = 0.8) and two (Z = 1.2 and ZT = 1.5) modified conditions.
Module Figure of Merit A s ( 10 7 ) B s ( 10 4 ) C s ( 10 4 ) A R ( 10 5 ) B R ( 10 4 ) C R ( 10 4 ) A t h ( 10 4 ) B t h ( 10 4 ) C t h
Z T m a x = 0.8 2 2 4.4 5 141 4418 2 3.78
Z T m a x = 1.2 2 2 5 5 136 4241 2 4.736
ZTmax = 1.5 2 2 5.3 5 127 3976 2 4.903
S ( T ) = A s T 2 + B s T + C s
R e ( T ) = A R T 2 + B R T + C R
R t h ( T ) = A t h T 2 + B t h T + C t h
TEG dimensions (L × W × H) are 8.7 × 8.7 × 1.4   ( mm 3 ) .
Table 3. Yearly average of direct, diffused, and global irradiance in Aalborg and Gibraltar.
Table 3. Yearly average of direct, diffused, and global irradiance in Aalborg and Gibraltar.
Location Direct   [ W . h / m 2 d a y ] Diffused   [ W . h / m 2 d a y ] Global   [ W . h / m 2 d a y ]
Aalborg322030936313
Gibraltar505026017651
Table 4. Comparison between measured and calculated values of matched electrical load and respected output power.
Table 4. Comparison between measured and calculated values of matched electrical load and respected output power.
Cold SideTemperature DifferenceR [ohm]P [mW/cm2]
MeasuredCalculatedMeasuredCalculated
TC = 25 °CdT = 15 °C12.8512.141.161.14
TC = 25 °CdT = 20 °C13.212.141.391.21
TC = 40 °CdT = 15 °C13.1412.141.231.17
TC = 40 °CdT = 20 °C13.712.141.451.29
Table 5. Thermal properties of thermoelectric materials and solar collector.
Table 5. Thermal properties of thermoelectric materials and solar collector.
-- Thermal   Conductivity   [ W / m K ] Heat   Capacity   [ J / g K ] Density   [ g / c m 3 ]
Copper3850.3858.9
B i 2 T e 3 *1.50.5447.29
* The average value for n and p-type materials.
Table 6. The daily energy production in Aalborg and Gibraltar for different values of ZT compared with the parabolic trough.
Table 6. The daily energy production in Aalborg and Gibraltar for different values of ZT compared with the parabolic trough.
LocationParabolic trough
[W. h/m2day]
ZT = 0.8
[W. h/m2day]
ZT = 1.2
[ W . h / m 2 d a y ]
ZT = 1.5
[ W . h / m 2 d a y ]
Aalborg402284430486
South of Spain631360564614

Share and Cite

MDPI and ACS Style

Karami Rad, M.; Omid, M.; Rajabipour, A.; Tajabadi, F.; Aistrup Rosendahl, L.; Rezaniakolaei, A. Optimum Thermal Concentration of Solar Thermoelectric Generators (STEG) in Realistic Meteorological Condition. Energies 2018, 11, 2425. https://doi.org/10.3390/en11092425

AMA Style

Karami Rad M, Omid M, Rajabipour A, Tajabadi F, Aistrup Rosendahl L, Rezaniakolaei A. Optimum Thermal Concentration of Solar Thermoelectric Generators (STEG) in Realistic Meteorological Condition. Energies. 2018; 11(9):2425. https://doi.org/10.3390/en11092425

Chicago/Turabian Style

Karami Rad, Meysam, Mahmoud Omid, Ali Rajabipour, Fariba Tajabadi, Lasse Aistrup Rosendahl, and Alireza Rezaniakolaei. 2018. "Optimum Thermal Concentration of Solar Thermoelectric Generators (STEG) in Realistic Meteorological Condition" Energies 11, no. 9: 2425. https://doi.org/10.3390/en11092425

APA Style

Karami Rad, M., Omid, M., Rajabipour, A., Tajabadi, F., Aistrup Rosendahl, L., & Rezaniakolaei, A. (2018). Optimum Thermal Concentration of Solar Thermoelectric Generators (STEG) in Realistic Meteorological Condition. Energies, 11(9), 2425. https://doi.org/10.3390/en11092425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop