Experimental parameters of the symmetrical and asymmetrical leg thermoelectric modules.
In this chapter, the impact of the shape of thermoelectric legs and parasitic contact resistances from metal electrodes and device wiring on thermoelectric figure of merit ZT is addressed. First section deals with the influence of the legs geometry on ZT. The shape of the legs is crucial in the thermoelectric performance of the thermoelectric devices. Unlike to conventional geometry thermoelectric legs, non-constant cross-section legs could help by lowering the overall thermal conductance of the device so as to increase the temperature gradient along legs, hence harnessing the Thomson effect, which is generally neglected in constant square cross-section thermoelectric legs. The final section is devoted to the electrical contact engineering of the device. Parasitic contact and wiring resistances play an important role in the performance of the device because they increase the isothermal resistance of the device. As the isothermal resistance of the device increases, the ZT decreases.
- asymmetrical legs
- energy harvesting
- heat recovery
- heat recycling
- thermoelectric module
Waste heat occurs in many areas of daily life (natural and industrial processes). In general, generation of 1 W of power using the state-of-the art heat engines requires about 3 W of energy input, which results in dumping into the environment the equivalent of about 2 W of power in the form of heat. Thereby, waste heat recovery to create a source of energy is an important technology for a stable supply of electricity and environmental protection. In particular, the heat loss due to low temperature, which is below 450 K, represents a large portion of waste heat emitted from an automobile and industry. For instance, in the US manufacturing sector alone, more than 3000 TW of waste heat energy is lost each year, an amount equivalent to more than 1.72 billion barrels of oil . Hence, efficiently reclaiming even a small portion of such waste heat would itself nearly satisfy the electricity needs of the planet . Thermoelectricity enabling directly converting heat into electricity using thermoelectric converters is a promising energy technology, which is under intense research to provide a solution to thermal energy recovery and management. However, wide-scale applications of this technology are limited due to the relative low thermoelectric efficiency of current materials and devices. For this reason, thermoelectric technology is only used in niche applications where its solid-state nature outweighs its poor efficiency. Therefore, an attractive and sustainable solution to the energy problem would be the development of solid-state thermoelectric devices which could help to recover this waste heat efficiently. The performance of a thermoelectric material is evaluated by the dimensionless figure of merit,
On the other hand, at device level evidently the transfer from material to device affects significantly the intrinsic figure of merit of the material; hence the different materials and substrates constituting the thermoelectric device as well as its architecture play an important role in the thermoelectric figure of merit of the device, as consequence, an effective thermoelectric figure of merit lower than the intrinsic one could be presented in the thermoelectric module . In this sense, significant enhancements in
Evidently, the performance of thermoelectric devices depends on many factors, such as the temperature difference between hot and cold plates, thermoelectric legs and device material properties, as well as the configuration and arrangement of thermoelectric legs. Among these factors, the geometry and configuration of the thermoelectric legs are crucial. In this sense, several studies under a purely theoretical context have predicted the effects of varying the geometry of the legs on the thermoelectric performance of the device. For instance, it has been reported that a smaller number of shorter legs have the potential to achieve the same power per unit module area as a greater number of longer legs . It has also been studied that the effect of the number of legs and their heights on the maximum output power and efficiency of the thermoelectric generator in order to find the optimum arrangements of legs . Moreover, it has been analyzed that the thermoelectric performance of the device by variations in the cross-section of the legs (legs with non-constant cross sections), and predictions show that legs with trapezoid legs (linear variation cross-section) have higher nominal power density than quadratic and exponential variations in cross-sectional legs [34, 35]. In addition, the tailoring of the geometry configuration of the legs in line with the device operating conditions has also been analyzed via thermoelectric legs tapering , as well as the thermoelectric performance of thermoelectric generators and coolers in relation to the geometry of non-constant cross-section legs has been formulated thermodynamically , and via genetic algorithms [38, 39, 40] for pyramid, cylindrical, and cuboid shapes with the aim of accomplishing device geometric optimization.
As mentioned above, the design and optimization of thermoelectric legs have previously been studied theoretically. However, experimental work confirming these theoretical predictions has not been carried out so far, mainly because of the difficulties involving in fabrication of thermoelectric legs with complex geometrical structures. Section 2 of the present chapter deals with the influence of the legs geometry on ZT. Unlike conventional geometry thermoelectric legs, non-constant cross-section legs could help by lowering the overall thermal conductance of the device so as to increase the temperature gradient along legs, hence harnessing the Thomson effect, which is generally neglected in constant square cross-section thermoelectric legs as it shows a complete theoretical and experimental analysis of asymmetrical legs. Section 3 is devoted to the electrical contact engineering of the device. Parasitic contact and wiring resistances play an important role in the performance of the device because they increase the isothermal resistance of the device. Finally, in Section 4, conclusions about of the results presented in this chapter are shown.
2. Theoretical and experimental analysis of symmetrical and asymmetrical legs
According to previous theoretical studies, the investigation of the geometric structure of thermoelectric legs is essential, as their geometry affects the performance of devices. It has been reported that asymmetrical shape of thermoelectric legs can lead to a decrease in the thermal and electrical conductance, which in turn improve the Seebeck voltage due to an increase in the temperature difference . In the present work, the temperature differences of the proposed asymmetrical legs are calculated by using the Fourier law. Figure 1a shows the proposed asymmetrical thermoelectric leg. For the sake of simplicity in the analysis, it is considered that such leg is constituted of four simpler geometrical legs connected thermally in parallel. By analyzing this simpler leg as shown in Figure 1b, it is possible to determine the temperature profile of the asymmetrical leg. Based on the Fourier law, the heat conduction along the simpler leg is defined by
Where, is the heat flux through the simpler leg, is the thermal conductivity, represents the slope of the pyramid, is the half-length of the smaller cross section end, and
Similarly, for a rectangular thermoelectric leg, the temperature profile is given by
where is the half-length of the bigger cross-section end.
The thermoelectric modules experimentally analyzed in this chapter have rectangular thermoelectric legs based on Bi2Te3 with a typical dimension of around 1.7 × 1.7 × 2.1 mm. Based on this dimensions and using Eqs. (2) and (3), Figure 3a shows the temperature profiles of an asymmetrical thermoelectric leg with a slant angle of 10° and a symmetrical thermoelectric leg under different heat fluxes, respectively. The upper solid lines (color online) represent the temperature profiles of the asymmetrical legs, whilst the lower dashed lines (color online) correspond to the temperature profiles of the symmetrical thermoelectric legs. As expected, as the heat flux decreases (), the temperature difference across the leg also decreases. Besides, a larger temperature difference is obtained across the asymmetrical leg as compared with the symmetrical one for a given heat flux. Therefore, by lowering the overall thermal conductance of the device via asymmetrical legs, the temperature gradient in the legs is increased, thus Seebeck voltage across terminals must be significantly increased. Besides, for simplicity it is worth to mention that the thermal conductivity
Figure 3b shows the temperature profiles of the asymmetrical thermoelectric leg as a function of the position for different slant angles of the pyramid for a given heat flux. As expected, as the slant angle increases, the temperature difference across the leg increases; however, such increase is limited to a critical slant angle (see inset), at which a maximum temperature difference across the leg can be achieved without affect the leg length; such limit angle is given by
Wherefrom, the limit angle for a thermoelectric leg with the length of 2.1 mm is around 22°, for this reason in Figure 2b as the slant angle achieves this value, the temperature rise increase drastically. It is worth to mention that in such modeling, convective effects were not taken into consideration because experimentally the device was tested under vacuum in order to avoid heat losses. Clearly, asymmetrical legs could help in two ways, by lowering the overall thermal conductance of the device so as to increase the temperature gradient in the legs, and by harnessing the Thomson effect, that depends on the temperature gradient in the legs and the temperature variation of Seebeck coefficient of the material in the operating temperature range, which is generally neglected in conventional rectangular thermoelectric legs. Thomson coefficient is given by , when a temperature gradient is imposed on a thermoelectric material
Although in this chapter the influences of the temperature dependence of the thermal conductivity and Seebeck coefficient are not considered in the model given by Eqs. (2) and (3) for the sake of simplicity; however, the main conclusions arising from the model are still valid. In addition, a simulation study has been carried out using a finite element simulation on 3D geometries for a thermoelectric device consisting of nine pairs of pyramidal legs and rectangular legs respectively using the dedicated software “Thermoelectric module of COMSOL Multiphysics in Seebeck mode”. Figure 4a, b show the simulation results of the temperature profiles and terminal voltages in the asymmetrical device respectively due to the presence of a given heat flux. For comparison, Figure 4c, d show the temperature profile and the terminal voltage in the symmetrical device respectively. Clearly, asymmetrical legs based thermoelectric device presents a higher open circuit terminal voltage than its rectangular counterpart as a consequence of the larger temperature difference generated in the legs because of their asymmetry.
The modules were made using p-type and n-Type Bi2Te3 thermoelectric materials available from Thermal Electronics Corp with ZT~1. For clarity, the full process so as to fabricate the thermoelectric legs is shown in Figure 5. The raw thermoelectric material in the form of a rod is covered by a layer of wax (wax-70) so as to hold the material in the cutting base during the cutting process, see Figure 5a. Then, it is cut into slices of 2.1 mm in thickness, which is the length of the thermoelectric leg as shown in Figure 5b. The wax is cleaned up using a warm solution of water and aquaclean-900 at a concentration of 10 mg/ml after slices cutting process as shown in Figure 5c. Next, a layer of Ni ranging from 0.5 and 1 μm is coated on both surfaces of the slide by electroplating. The Ni layer works as a diffusion barrier between the solder (Sn/Pb 60/40) and the thermoelectric legs, see Figure 5d. Prior to cutting, the Ni-electroplated slide is again fixed on a graphite plate using wax (wax-70) with the aim of keeping the legs during cutting, as you can see in Figure 5e. Subsequently, thermoelectric legs with the regular geometry of 1.7 × 1.7 × 2.1 mm are obtained using a circular saw cutting machine (Accutom-100), see Figure 5f. Finally, the wax is removed from legs as previously indicated in Figure 5c.
In order to obtain the asymmetrical legs, a similar process as described above is employed. However, in this case, a tilted base is used to hold the graphite plate during cuttings. Hence, so as to obtain the asymmetrical geometry (truncated square pyramid), it was necessary to do a cutting in every face of the rectangular leg by rotating the graphite tilted base 90° during each cutting. The slope of the base depends on the desired slant angle in the thermoelectric legs as shown in Figure 6a, b, clearly, the graphite plate is tilted at an angle of 10°. Besides, it has been attempted to fabricate the legs with a slant angle close to the critical angle of 22° so as to maximize the thermal gradient. However, it has been observed that angles higher than 10° produce legs with a small cross-section in the thin end. Such legs tend to be very fragile because of the mechanical properties of Bi2Te3, and they fracture during the assembling process due to the pressure and thermal treatments applied. Figure 6c, d show the side and top view of the asymmetrical legs used for the fabrication of the thermoelectric device.
Figure 7a, b show the images of the fabricated modules with symmetrical and asymmetrical legs, respectively, wherefrom differences in the geometry of the thermoelectric legs can be observed by comparing both images. At first instance, the performance of the devices has been evaluated by means of the hot-plate method. In such method, the thermoelectric modules are installed between heating and cooling plates. Next, a dc voltage with regular increments is applied to the 130 Ω resistive heater of the heating plate by using the BK Precision 9184 dc power supply. This action induces a constant heat flow
Parallel, the open circuit thermal voltage is measured along with the temperature rise across the device, such voltage is detected by using the Metrohm Autolab B. V. system; then, it is also plotted
On the other hand, when a thermoelectric device is connected to any load it is desirable that such device be able to transfer the greatest amount of power to the load. In this sense, applying the Theory of Maximum Power Transfer assumes a simple electrical circuit with a voltage source
The current is
Then is considered a variation in power when the load resistance is changed, but internal resistance of thermoelectric device
In this sense, in order to evaluate the maximum power given by the modules the maximum output power between the asymmetrical and the symmetrical modules has been evaluated by way of the maximum power transfer theorem. In this case, an identical heat flux of around 4 mW/mm2 was supplied to both modules by applying an electrical current of 90 mA to the Ohmic heater; once modules reached the steady state different load resistances were connected into the module and the voltage across the resistance was recorded, by using those voltages and resistances the output power was estimated. Figure 9 shows the obtained results, it can be observed that asymmetrical module delivers more power than the symmetrical one once the load resistance equals the device internal resistance. In Fact, the asymmetrical thermoelectric module shows to have almost twofold the maximum delivered power as compared to conventional one with a constant square cross-section. Besides, by estimating the maximum available power per unit amount of material (mass of the legs) it has been obtained 433 μW/gram and 1.57 mW/gram for the symmetrical and the asymmetrical modules, respectively.
Evidently, these modules apparently present low output power, however by comparing these modules with several commercially available they have very competitive output power values . For instance, by extrapolating the data shown in Ref.  to ΔT = 20°C a TEG module based on Bi2Te3 model FERROTEG 9501/71/040B with 71 pairs, and 22 mm × 22 mm generates a maximum output power around 1.5 mW. In our case, for modules with only nine pairs, we obtain 0.3 and 0.5 mW for symmetrical and asymmetrical modules, respectively. Nevertheless, by the projection of our modules to 71 pairs we would obtain 2.36 and 3.94 mW, respectively. Besides, if we compare our module against TEG-FERROTEG 9500/127/100B module based on Bi2Te3 with 127 pairs, and 40 × 40 mm under ΔT = 20°C, which delivers an output power around 2.5 mW, we would obtain by a similar extrapolation 4.23 and 7.05 mW for symmetrical and asymmetrical modules, respectively, under ΔT = 20°C. It is worth to mention that ΔT scale in Ref.  is in logarithmic scale, so it can be closely compared to ΔT = 20°C and ΔT = 30°C, for the asymmetrical module.
The thermoelectric figure of merit of the fabricated modules was also evaluated by using impedance spectroscopy technique . In this method, the thermoelectric figure of merit is determined by measuring the adiabatic and isothermal responses of the module under electrical excitation. Under the adiabatic condition at steady state (i.e., ω = 0), the total impedance of the module can be written as 
Where is the thermoelectric resistance, and it refers to the resistance of the device which is the result of the temperature difference induced between the ends of the sample due to the Peltier effect, is the isothermal resistance, and it is the resistance of the device excluding thermal effects but including parasitic resistances such as contact and wire resistances, .
According to the Harman method, the thermoelectric figure of merit is given by the ratio between the thermoelectric resistance and the isothermal resistance of the system ; hence, by dividing Eq. (6) by ,
As well as, in terms of the thermoelectric resistance
Therefore, by using Eqs. (7) or (8), the effective thermoelectric figure of merit of a device can be accomplished. The adiabatic and isothermal resistances can be easily accessed via electrical impedance measurements . Likewise, parasitic resistances from module with symmetrical and asymmetrical legs are evaluated by applying Transmission Line Method (TLM)  and so they can be removed from
During the measurements, the samples were isolated and suspended to provide adiabatic conditions in a similar way as required in the Harman method . Figure 10a, b show the experimental electrical impedance curves obtained for the symmetrical and asymmetrical leg modules, respectively. In both curves, the thermoelectric, adiabatic, and isothermal resistances are indicated in order to access to their respective values. Nevertheless, it is well known that a material has more than one contribution to its impedance response, which is often the case of thermoelectric materials where thermoelectric impedance, isothermal impedance, and contact impedance have distinct contributions. Hence, one can witness more than one semi-circle, often overlapping each other which makes impossible to distinguish them. One of the ways to model such a behavior in a simple model can be using in three series–parallel RC elements circuit. In Figure 10a, b, the solid line corresponds to the obtained fitting results. For clarity, such resistance results, as well as the effective thermoelectric figure of merit of the symmetrical and asymmetrical modules are shown in Table 1. Evidently, the thermoelectric figure of merit of the asymmetrical module is almost two-fold the thermoelectric figure of merit of the symmetrical module, such result confirms the enhanced thermoelectric performance of the asymmetrical module as a consequence of the larger temperature rise generated in the legs because of their asymmetry. Hence, harnessing of the Thomson coefficient via asymmetrical legs could be an important strategy in order to accomplish thermoelectric devices with enhanced performance.
On the other hand, it is worth to mention that the present experimental research is mainly focused on the development of devices for applications at room temperature (i.e. 300 K), in that case, it is not necessary to measure the temperature dependence of ZT. Besides, our devices are based on P and N-type Bi2Te3, it is well known that such materials present an optimal thermoelectric performance at around room temperature; hence, operation of such materials must be well below 100°C, so an operation condition above this temperature will damage the device because by applying an excessive heat flux it could damage the device due to the melting of the weld joining the thermoelectric legs. In this sense, it is not possible to operate such device under a high-temperature rise away from room temperature would affect seriously their performance.
Moreover, according to the values shown in Table 1, evidently, the parasitic electrical resistances play an important role in the performance of the device. For instance, if we take into consideration only parasitic contact effects (i.e. parasitic electrical contact resistance between legs and ceramic plates) and neglect the effect of parasitic resistances generated by cable wiring, a value of 0.43 and 0.73 on
3. Impact of parasitic contact electrical resistances on
ZT of the thermoelectric device
Thermoelectric device engineering involves the formation of several intrinsic parasitic resistances that affect the thermoelectric module performance. In this sense, the TLM has been applied to discard the parasitic resistances and demonstrate that the increase on
Figure 11a shows different lengths in the symmetrical and asymmetric thermoelectric legs as well as their respective electrical resistances as a function of length. The total measured resistance consists of several components:
Where RW1 y RW2 are wiring resistances, RC1 y RC3 are contact resistances due to metal contacts, RC2 is associated with the metallic contact between the junction of the P-type and N-type thermoelectric legs, and RP-TE and RN-TE define the internal resistance of the P-type and N-type thermoelectric legs, respectively.
Therefore, the total parasitic electrical resistance Rp is given by the contact resistance Rc and the wiring resistance Rw, then Eq. (10) can be rewritten as:
By way of the TLM, it is possible to measure the total parasitic electrical resistance Rp. In this sense, thermoelectric legs have been fabricated with 2, 3, and 4 mm in length. The TLM is a technique used to determine the contact resistance between a metal and a semiconductor. First, the electrical resistance is measured for each length and then each resistance is presented as a function of length as shown in Figure 12. In the limit of a zero-length resistor, the residual resistance would be just the contact resistance. Then can be found from the graph by extrapolating back to L = 0. Then, the parasitic resistance of the P-N junction is the sum of such interceptions; hence, the total parasitic resistance of the device is estimated by multiplying this value by the number of P-N junctions in the device (in this case, 9 P-N pairs).
The total parasitic resistance
Where, y values are defined by slope interception with y-axes, as mentioned above. Such values are shown as inset table of Figure 12. This result is the total parasitic resistance of thermoelectric module with nine pairs that includes the wiring resistance of wire used during measurements as well as contact resistance.
In addition, by way of the four-probe AC method, it is possible to measure the wiring resistance
Now, contact resistance Rc of the thermoelectric module with symmetrical legs is calculated as:
By applying a similar procedure, it is possible to measure the parasitic resistance Rp using the TLM in the asymmetric device, in this case, the resistance as a function of the length is shown in Figure 14.
The total parasitic resistance
Where, y values are defined by slope interception with
Now, contact resistance
In this particular case, has been obtained around 116 mΩ y 183 mΩ for the symmetrical and asymmetrical device, respectively. Using this information,
In the present chapter, it has been done an experimental demonstration of the influence of the device legs geometry as well as parasitic electrical contact resistance on ZT. Results prove that asymmetrical thermoelectric module shows to have almost twofold the thermoelectric figure of merit as compared to conventional one with a constant square cross-section. Thermal analysis of the device via analytical, as well as numerical modeling unveils an increment in the temperature gradient and Seebeck voltage across the device with asymmetrical thermoelectric legs. Such result confirms that the thermoelectric enhancement is due to the harnessing of Thompson effect which is normally neglected in rectangular legs devices. Additionally, the impact of parasitic electrical contact and wiring resistances on the thermoelectric module performance is shown. In this sense, a significant decrement on ZT due to parasitic effects is observed. Thereby, the general results of the present chapter experimentally prove that geometrical configuration of the device legs can improve significantly the thermoelectric performance of the device opening a new route to the development of enhanced performance thermoelectric modules via device engineering.
This work was supported by the National Council for Science and Technology-Conacyt Mexico, through the Grant for fundamental research No. 241597 and national issues No.1358. A.F.M. thanks to Conacyt Mexico for fellowship, as well as the thermoelectric laboratory at the Cardiff University for facilities.
Conflict of interest
Authors declare no conflicts of interest.
Tritt TM, Subramanian MA. Thermoelectric materials, phenomena, and applications: A bird’s eye view. MRS Bulletin. 2016; 31(3):188. DOI: 10.1557/mrs2006.44
Pop E. Energy dissipation and transport in nanoscale devices. Nano Research. 2010; 3:147-169
Wang S, Fu F, She X, Zheng G, Li H, Tang X. Optimizing thermoelectric performance of Cd-doped β-Zn4Sb3 through self-adjusting carrier concentration. Intermetallics. 2011; 19:1823-1830
Dresselhaus MS, Chen G, Tang MY, Yang RG, Lee H, Wang DZ, et al. New directions for low-dimensional thermoelectric materials. Advanced Materials. 2007; 19:1043-1053
Hicks LD, Dresselhaus MS. Thermoelectric figure of merit of a one-dimensional conductor. Physical Review B. 1993; 47:16631-16634
Hicks LD, Dresselhaus MS. Effect of quantum-well structures on the thermoelectric figure of merit. Physical Review B. 1993; 47:12727-12731
Chen G, Dresselhaus MS, Dresselhaus G, Fleurial J-P, Caillat T. Recent developments in thermoelectric materials. International Materials Review. 2003; 48:45-66
Snyder GJ, Toberer ES. Complex thermoelectric materials. Nature Materials. 2008; 7:105-114
Venkatasubramanian R, Siivola E, Colpitts T, O’Quinn B. Thin-film thermoelectric devices with high room-temperature figures of merit. Nature. 2001; 413:597-602
Sun X, Zhang Z, Dresselhaus MS. Theoretical modeling of thermoelectricity in Bi nanowires. Applied Physics Letters. 2016; 4005:26-29
Zhu T, Ertekin E. Phonon transport on two-dimensional graphene/boron nitride superlattices. Physical Review B. 2014; 90:195209
Chen Z-G, Han G, Yang L, Cheng L, Zou J. Nanostructured thermoelectric materials: Current research and future challenge. Progress in Natural Science: Materials International. 2012; 22:535-549
Poudel B, Hao Q, Ma Y, Lan Y, Minnich A, Yu B, et al. High-thermoelectric performance of nanostructured bismuth antimony telluride bulk alloys. Science. 2008; 320:634-638
Oh T-S. Thermoelectric characteristics of p-type (Bi,Sb)2Te3/(Pb,Sn)Te functional gradient materials with variation of the segment ratio. Journal of Electronic Materials. 2009; 38(7):1041
Vining CB. An inconvenient truth about thermoelectrics. Nature Materials. 2009; 8:83-85
Hsu FK. Cubic AgPbmSbTe2+m: Bulk thermoelectric materials with high figure of merit. Science. 2004; 303:818-821
Zhao XB, Ji XH, Zhang YH, Zhu TJ, Tu JP, Zhang XB, et al. Bismuth telluride nanotubes and the effects on the thermoelectric properties of nanotube-containing nanocomposites. Applied Physics Letters. 2016:062111
Ni HL, Zhao XB, Zhu TJ, Ji XH, Tu JP. Synthesis and thermoelectric properties of Bi2Te3 based nanocomposites. Journal of Alloys and Compounds. 2005; 397:317-321
Slack GA, Hussain MA. The maximum possible conversion efficiency of silicon-germanium thermoelectric generators. Journal of Applied Physics. 2016; 2694
Kim W, Zide J, Gossard A, Klenov D, Stemmer S, Shakouri A, et al. Thermal conductivity reduction and thermoelectric figure of merit increase by embedding nanoparticles in crystalline semiconductors. Physical Review Letters. 2006; 96:45901
Alvarez-Quintana J. Impact of the substrate on the efficiency of thin film thermoelectric technology. Applied Thermal Engineering. 2015; 84:206-210
Salvador JR, Cho JY, Ye Z, Moczygemba JE, Thompson AJ, Sharp JW, et al. Conversion efficiency of skutterudite-based thermoelectric modules. Physical Chemistry. 2014; 16:12510-12520
Taguchi K, Terakado K, Ogusu M, Matumoto A, Kayamoto T, Okura K, et al. Linear shaped Si-Ge thermoelectric module. In: Reference Number F2000A045 in Proceedings of Seul 2000 FISITA World Automotive Congress; 2000. pp 1-5
Joshi G, He R, Engber M, Samsonidze G, Pantha T, Dahal E, et al. NbFeSb-based p-type half-Heuslers for power generation applications. Energy & Environmental Science. 2014; 7:4070-4076
Fujisaka T, Sui H, Suzuki RO. Design and numerical evaluation of Cascade-type thermoelectric modules. Journal of Electronic Materials. 2013; 42:1688-1696
Sun T, Peavey JL, David Shelby M, Ferguson S, O’Connor BT. Heat shrink formation of a corrugated thin film thermoelectric generator. Energy Conversion and Management. 2015; 103:674-680
Madan D, Chen A, Wright PK, Evans JW. Printed Se-doped MA n-type Bi2Te3 thick-film thermoelectric generators. Journal of Electronic Materials. 2012; 41(6):1481
Chen A, Madan D, Wright PK, Evans JW. Dispenser-printed planar thick-film thermoelectric energy generators. J Micromechanics Microengineering. 2011; 104006:21
Zheng XF, Liu CX, Yan YY, Wang Q. A review of thermoelectrics research – Recent developments and potentials for sustainable and renewable energy applications. Renewable and Sustainable Energy Reviews. 2014; 32:486-503
Aswal DK, Basu R, Singh A. Key issues in development of thermoelectric power generators: High figure-of-merit materials and their highly conducting interfaces with metallic interconnects. Energy Conversion and Management. 2016; 114:50-67
Hamid Elsheikh M, Shnawah DA, Sabri MFM, Said SBM, Haji Hassan M, Ali Bashir MB, et al. A review on thermoelectric renewable energy: Principle parameters that affect their performance. Renewable and Sustainable Energy Reviews. 2014; 30:337-355
Fateh H, Baker CA, Hall MJ, Shi L. High fidelity finite difference model for exploring multi-parameter thermoelectric generator design space. Applied Energy. 2014; 129:373-383
Hodes M. Optimal pellet geometries for thermoelectric power generation. IEEE Trans Components Packag Technol. 2010; 33:307-318
Sahin AZ, Yilbas BS. The thermoelement as thermoelectric power generator: Effect of leg geometry on the efficiency and power generation. Energy Conversion and Management. 2013; 65:26-32
Shi Y, Mei D, Yao Z, Wang Y, Liu H, Chen Z. Nominal power density analysis of thermoelectric pins with non-constant cross sections. Energy Conversion and Management. 2015; 97:1-6
Yilbas BS, Ali H. Thermoelectric generator performance analysis: Influence of pin tapering on the first and second law efficiencies. Energy Conversion and Management. 2015; 100:138-146
Ali H, Sahin AZ, Yilbas BS. Thermodynamic analysis of a thermoelectric power generator in relation to geometric configuration device pins. Energy Conversion and Management. 2014; 78:634-640
Xuan XC, Ng KC, Yap C, Chua HT. Optimization of two-stage thermoelectric coolers with two design configurations. Energy Conversion and Management. 2002; 43:2041-2052
Cheng Y-H, Lin W-K. Geometric optimization of thermoelectric coolers in a confined volume using genetic algorithms. Applied Thermal Engineering. 2005; 25:2983-2997
Cheng Y-H, Shih C. Maximizing the cooling capacity and COP of two-stage thermoelectric coolers through genetic algorithm. Applied Thermal Engineering. 2006; 26:937-947
Mu Y, Chen G, Yu R, Li G, Zhai P, Li P. Effect of geometric dimensions on thermoelectric and mechanical performance for Mg2Si-based thermoelectric unicouple. Materials Science in Semiconductor Processing. 2014; 17:21-26
Ioffe AF. Semiconductor thermoelements and thermoelectric cooling. Infosearch London. 1957
Yamashita O. Resultant Seebeck coefficient formulated by combining the Thomson effect with the intrinsic Seebeck coefficient of a thermoelectric element. Energy Conversion and Management. 2009; 50(9):2394
Salerno D. Ultralow voltage energy harvester uses thermoelectric generator for battery-free wireless sensor. Journal of Analog Innovation. 2010; 20(3):1
García-cañadas J, Min G. Impedance spectroscopy models for the complete characterization of thermoelectric materials. Journal of Applied Physics. 2014; 116:174510
Iwasaki H, Koyano M, Hori H. Evaluation of the figure of merit on thermoelectric materials by Harman method. Japanese Journal of Applied Physics. 2002; 41:6606
Reeves GK, Harrison HB. Obtaining the specific contact resistance from transmission line model measurements. IEEE Electron Device Letters. 1982; 3:111-113