Next Article in Journal
Liquid Phase Assisted Superplastic Deformation of TiO2-Doped ZTA Ceramics
Next Article in Special Issue
TiO2–SrTiO3 Biphase Nanoceramics as Advanced Thermoelectric Materials
Previous Article in Journal
Effects of Cell Temperature and Reactant Humidification on Anion Exchange Membrane Fuel Cells
Previous Article in Special Issue
p-n Control of AlMgB14-Based Thermoelectric Materials by Metal Site Occupancy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhanced Thermoelectric Cooling through Introduction of Material Anisotropy in Transverse Thermoelectric Composites

1
Key Laboratory of Traffic Safety on Track, Ministry of Education, School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China
2
Department of Mechanical Engineering, Temple University, Philadelphia, PA 19122, USA
3
Joint International Research Laboratory of Key Technology for Rail Traffic Safety, School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China
4
National and Local Joint Engineering Research Center of Safety Technology for Rail Vehicle, School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(13), 2049; https://doi.org/10.3390/ma12132049
Submission received: 2 June 2019 / Revised: 18 June 2019 / Accepted: 21 June 2019 / Published: 26 June 2019
(This article belongs to the Special Issue Advanced Thermoelectric Materials)

Abstract

:
Transverse thermoelectric materials can achieve appreciable cooling power with minimal space requirement. Among all types of material candidates for transverse thermoelectric applications, composite materials have the best cooling performance. In this study, anisotropic material properties were applied to the component phase of transverse thermoelectric composites. A mathematical model was established for predicting the performance of fibrous transverse thermoelectric composites with anisotropic components. The mathematical model was then validated by finite element analysis. The thermoelectric performance of three types of composites are presented, each with the same set of component materials. For each type of component, both anisotropic single-crystal and isotropic polycrystal material properties were applied. The results showed that the cooling capacity of the system was improved by introducing material anisotropy in the component phase of composite. The results also indicated that the orientation of the anisotropic component’s property axis, the anisotropic characteristic of a material, will significantly influence the thermoelectric performance of the composite. For a composite material consisting of Copper fiber and Bi2Te3 matrix, the maximum cooling capacity can vary as much as 50% at 300 K depending on the property axis alignment of Bi2Te3 in the composite. The composite with Copper and anisotropic SnSe single crystal had a 51% improvement in the maximum cooling capacity compared to the composite made of Copper and isotropic SnSe polycrystals.

1. Introduction

Thermoelectric material can transform thermal energy into electrical energy and vice versa. Its applications can be found in various fields such as temperature measurement, power generation, temperature control, etc. Conventional thermoelectric devices consist of both n-type and p-type thermoelectric legs. Within each leg, the heat flux and electrical current flow are parallel to each other [1,2]. Transverse thermoelectric devices make use of the transverse Seebeck effect so that the electrical current and heat flux flow perpendicular to each other [3]. There are four main types of transverse thermoelectric materials, which are anisotropic single-crystal material, polycrystal material with engineered anisotropy, anisotropic organic thin-film thermoelectrics, and anisotropic thermoelectric composites. The single-crystal material anisotropy is the result of an unsymmetrical lattice structure [4]; the material anisotropy in polycrystal is contributed to by both the lattice structure of the grain crystal and the grain geometry [5]. The organic materials anisotropic property is caused by the preferred alignment of polymer chains [6,7,8]. Studies have shown that the anisotropic thermoelectric composites can provide the best system efficiency among all of the candidates [9].
The transverse thermoelectric composite can further be categorized into layered and fibrous thermoelectric composites. These two types of composites can also be seen as two-dimensional (2D) and one-dimensional (1D) inclusion composites. For layered composites with isotropic components, Babin et al. [10] established a mathematical model for fast prediction of dimensionless transverse thermoelectric figure of merit (ZtransT), while other researchers performed finite element simulations and experimental tests and validated such a mathematical model [11,12,13,14,15]. Similar to layered composites, fibrous composites can also provide appreciable ZtransT values. A mathematical model was established by Qian [16] to study the thermoelectric performance of fibrous composites.
Until now, most studies on transverse thermoelectric composites assumed isotropic material properties in the component phase. However, many studies have already shown that thermoelectric materials can exhibit anisotropic properties, such as Bi2Te3 [5,17], SnSe [18], and organic thermoelectric PEDOT:PSS (poly(3,4-ethylenedioxythiophene) polystyrene sulfonate) [19,20], etc. As a result, anisotropic material properties were applied into the component phase of layered transverse thermoelectric composites, and the results showed that the maximum ZtransT could be improved by introducing material anisotropy in a polycrystal [21]. However, the effect of material anisotropy on fibrous composite still remains unknown.
In this study, anisotropic material properties were applied to the component phase of fibrous transverse thermoelectric composites. A mathematical model was established for predicting the effective material properties and thermoelectric cooling capacities of the composite. Finite element simulations were conducted to validate the mathematical model. A case study was conducted on a fibrous composite material containing Copper fibers and an anisotropic Bi2Te3 matrix, where the possible aggregate (anisotropic) properties for Bi2Te3 were calculated and applied into the component phase. The influence of the material anisotropy and the material property axis alignment in the composite were then discussed with respect to the maximum ZtransT and maximum cooling capacity of the composites. In the next step, the cooling performances of both the layered and the fibrous composites were compared. The thermoelectric properties of three representative composites, which were Bi2Te3/Copper, In4Se2.25/Copper, and SnSe/Copper, were applied into the component phase of the composites under different material property axis alignment configurations and extreme aggregate properties of the anisotropic component. The influence of material anisotropy on the enhancement of the maximum cooling capacity and maximum ZtransT of the transverse thermoelectric composites was demonstrated. The cooling capacities of the layered and fibrous composites were also compared. The derived mathematical model can serve as an efficient tool for selecting high-performance fibrous transverse thermoelectric composites, while the comparison between fibrous and layered thermoelectric composites can inspire and assist thermoelectric researchers in designing higher performance thermoelectric devices.

2. Mathematical Model for Fibrous Thermoelectric Composites with Anisotropic Components

In this section, analytical equations are derived for the effective properties of a transverse thermoelectric composite. In the mathematical model, unit cell structure is used to reduce the complexity of this problem. This unit cell model has been proven to be a convenient tool in studies on the effective properties of a material with periodical structures [16,22]. The schematic of a fibrous thermoelectric composite unit cell is shown in Figure 1, where a fiber (F) with a square cross-section was placed at the corner of a cubic matrix (M). The square fiber was used in the unit cell model for mathematical simplicity. The contact between the matrix and fibrous material was assumed to be ideal so that electrical and thermal contact resistance were not considered.
When anisotropic material properties are introduced into the component phase (i.e., fibrous phase and matrix phase) of a composite material, the alignment of the material local property axis with respect to the composite material coordinate system plays a key role in the effective properties of the composite. Many 1D inclusions, such as carbon nanotubes [23,24], have different material properties along and perpendicular to the fiber axis direction. Therefore, in this study, we used the subscript “p” to represent the properties on the cross-sectional plane of the fiber and subscript “a” to represent properties along fiber axis (Figure 1a). As for the matrix phase, we used the material local axis system “uvw” to describe the anisotropic material properties (Figure 1b). There are four sets of coordinate systems in Figure 1, which are the local material coordinate system of the fibrous and matrix phase, the coordinate system of composite C1(F + M1), and the coordinate system of unit cell C2(F + M1 + M2). In Figure 1, the coordinate axes u, x1, and x2 align parallel to each other; v, y1, and y2 align parallel to each other, and a, w, z1, and z2 align parallel to each other. The cross-sectional plane of fiber (p–p) is parallel with the u–v plane of the matrix material.
The effective properties of the unit cell can be derived in two steps as illustrated in Figure 1c. In the first step, the square fiber (F) is combined with part of the matrix phase (M1) to form a rectangular composite C1, where the matrix block has the same width and height as the fibrous phase. In the second step, the composite C1 is combined with the rest of the matrix to form the unit cell. The effective thermal and electrical conducting properties of the unit cell can be calculated using Kirchhoff’s law, and the effective Seebeck coefficient can be calculated using Thevenin’s theorem [25]. The material properties of the composite C1 can be derived according to Equation (1):
ρ x 1 = n ρ M u + ρ F p 1 + n λ x 1 = λ F p λ M u ( 1 + n ) n λ F p + λ M u S x 1 = S F p λ M u + n S M u λ F p n λ F p + λ M u ρ y 1 = ρ M v ρ F p ( 1 + n ) ρ M v + n ρ F p λ y 1 = λ F p + n λ M v 1 + n S y 1 = S F p ρ M v + n S M v λ F p n ρ F p + ρ M v ρ z 1 = ρ M w ρ F a ( 1 + n ) ρ M w + n ρ F a λ z 1 = λ F a + n λ M w 1 + n S z 1 = S F a ρ M w + n S M w λ F a n ρ F a + ρ M w
The symbols ρ, λ, S, and f stand for electrical resistivity, thermal conductivity, Seebeck coefficient, and fiber volume fraction. The subscripts p, a, u, v, w, xi, yi, and zi stand for component material’s properties along each material’s local coordinate systems. The material properties of the unit cell can be derived according to Equation (2):
ρ x 2 = ρ x 1 ρ M u f ( ρ M u ρ x 1 ) + ρ x 1 λ x 2 = ( λ x 1 + ( 1 f 1 )   λ M u ) f S x 2 = S x 1 ρ M u + ( 1 f 1 ) S M u ρ x 1 ( 1 f 1 ) ρ x 1 + ρ M u ρ y 2 = f ( ρ y 1 ρ v ) + ρ M v λ y 2 = λ y 1 λ M v f ( ( 1 f 1 ) λ y 1 + λ M v ) S y 2 = S y 1 λ M v + ( 1 f 1 ) S M v λ y 1 ( 1 f 1 ) λ y 1 + λ M v ρ z 2 = ρ z 1 ρ M w f ( ρ M w ρ z 1 ) + ρ z 1 λ z 2 = ( λ z 1 + ( 1 f 1 )   λ M w ) f S z 2 = S z 1 ρ M w + ( 1 f 1 ) S M w ρ z 1 ( 1 f 1 ) ρ z 1 + ρ M w
In Figure 1, the fiber axis aligns parallel to the z2-axis of the composite. When fibers are tilted aligned in the composite, the effective properties of the composite can be calculated through matrix transformation. For example, if the fibers in Figure 1 are rotated by an angle of θ around x2-axis into the configuration in Figure 2, the effective properties of the composite can be calculated as [26]:
P x y z = [ P x 2 0 0 0 P y 2 c o s 2 ( θ ) + P z 2 s i n 2 ( θ ) 1 2 ( P z 2 P y 2 ) s i n ( 2 θ ) 0 1 2 ( P z 2 P y 2 ) s i n ( 2 θ ) P z 2 c o s 2 ( θ ) + P y 2 s i n 2 ( θ ) ]
In a transverse thermoelectric material, heat flux and electrical current flow perpendicular to each other. Based on Figure 2, we assume the electrical current flows in the y-direction and heat flux flows in the z-direction. The transverse thermoelectric figure of merit is defined as Z trans T = S z y 2 T ρ y y λ z z [10,27], where Szy is the transverse Seebeck coefficient, ρyy is electrical resistivity in the y-direction, λzz is thermal conductivity in the z-direction, T is the operating temperature. The term Szy, ρyy, and λzz values can be calculated using Equation (1) to Equation (3).
Figure 2 can also represent a fibrous transverse thermoelectric device operating under cooling mode, where the top surface serves as a cooling surface with a temperature of Tc and the bottom surface serves as a heat sink with a temperature of Th. When the electrical current flows in the y-direction, the transverse Peltier effect will trigger heat flux along the z-direction. The maximum cooling capacity of the system is related to the ZtransT value, according to previous studies [10,28]. In this study, the Th was set to be the measuring temperature of material properties, which are shown in Table 1.
In the next step, experimentally measured material properties are implemented into the derived equations. Previous studies [9] have shown that high-performance transverse thermoelectric composites usually consist of one semiconducting thermoelectric phase and one highly conductive phase, where the thermoelectric phase provides high Seebeck coefficient and low thermal conductivity, the conducting phase provides low electrical resistivity. Therefore, we chose anisotropic Bi2Te3 polycrystal as the matrix phase and Copper as the fibrous phase for the composite. The properties of these materials are shown in Table 1. According to Table 1, the Copper phase was isotropic and Bi2Te3 was anisotropic with planar (u–v plane) isotropy.
When anisotropic material properties are implemented in the component phase of a composite, the alignment of component material’s local property coordinate system, with respect to the composite coordinate system, will affect the effective properties of the composite. In Figure 1, the u–v plane of the matrix material aligns perpendicular to the fiber axis. If the material property coordinate in Figure 1b is rotated around the u-axis 90 degrees, the anisotropy plane of Bi2Te3 will be parallel to the fiber axis direction. For the convenience of analysis, we shall refer to the two configurations mentioned above as configuration I and configuration II. In configuration I, the material property axes u, v, and w are parallel to the composite axes x2, y2, and z2, respectively. In configuration II, the material property axes u, v, and w are parallel to the composite axes x2, z2, and y2, respectively. The corresponding ZtransT and maximum cooling capacities (∆Tmax) based on these two configurations were calculated using the derived mathematical model with respect to fiber rotation angle, θ, and fiber volume fraction, f. The results are shown in Figure 3.
Based on Figure 3, the ZtransT and ∆Tmax values exhibited similar trends for both configurations. The maximum values appearred at rotation angles between 75 and 85 degrees. The change in fiber volume fraction under a constant rotation angle had a minor influence on the ZtransT and ∆Tmax values. The maximum ZtransT and ∆Tmax values in configuration I were 0.29 and 34 K, while the values of maximum ZtransT and ∆Tmax in configuration II were 0.22 and 27 K. These results indicate a 26% difference in the peak thermoelectric performance among the two configurations. Hence, the alignment of the material’s local property coordinates with respect to the composite coordinate system had significant influence on the composite’s cooling performance.
Many thermoelectric materials are polycrystals. These polycrystals may exhibit certain degrees of texture, and thus, anisotropic properties as a result of the fabrication processes. There exists theoretical models that can correlate the aggregate properties of polycrystals with its single-crystal material property. Among these theoretical models, the Voigt model and Reuss model can provide upper and lower bounds for the aggregate properties of polycrystals, respectively [32,33,34]. In this study, we take experimentally measured values from a single-crystal Bi2Te3 (Table 1) and calculate possible aggregate material properties using Reuss and Voigt models. The Bi2Te3 polycrystal has isotropic Seebeck coefficients and exhibits planar isotropy (u–v plane) in electrical resistivity and thermal conductivity based on existing experimental results [5,18,35,36]. The anisotropy ratio terms r λ and   r ρ are used to relate the material properties in different axis directions, such that r λ = λ u / λ w ,   r ρ = ρ w / ρ u . By importing these hypothetical aggregate properties into the Bi2Te3 phase, the maximum cooling capacities were calculated for the Bi2Te3 matrix/Cu fiber composite. The results are shown in Figure 4. According to the definition of r λ and   r ρ , the bottom-left corner of Figure 4a,b represents the maximum cooling capacity of a polycrystal Bi2Te3 with isotropic material properties. The results in Figure 4 show that by introducing material anisotropy in a fibrous transverse thermoelectric composite, the cooling capacity of the composite can be improved by 7% (Figure 4b, Voigt model) and 21% (Figure 4a, Reuss model).

3. Finite Element Simulation for Fibrous Thermoelectric Composite with Anisotropic Components

In order to validate the effectiveness of the mathematical model, finite element simulations were carried out using COMSOL (COMSOL Multiphysics, COMSOL Inc., Burlington, MA, USA). A unit cell model based on Figure 1c was constructed using polycrystal Bi2Te3 as the matrix phase and Copper as the fiber phase assuming the isotropy plane of Bi2Te3 was perpendicular to the fiber axis. For all finite element models involved in this study, mesh convergence studies were carried out and the convergence error was limited to below 3%.
According to the unit cell model shown in Figure 1c, square fiber was used for mathematical simplicity. In the finite element simulations, both cylindrical fiber and square fiber were investigated so as to explore the effect of fiber shapes. The effective properties of the unit cells are shown in Figure 5.
According to Figure 5, the finite element simulation results and mathematical model agreed well with each other. The conductivities of the unit cell increased as the volume fraction of the fibers increased, since the fiber was more conductive than the matrix material. In Figure 5f, the effective Seebeck coefficient had negative values at a low fiber volume fraction, but became positive as the fiber volume fraction increased. This was because the Bi2Te3 phase had a negative Seebeck coefficient and Copper had a positive Seebeck coefficient. The effective Seebeck coefficient was dominated by the matrix phase at low fiber volume fractions, and by the fibrous phase at high fiber volume fractions.
Some minor deviations were observed between the mathematical model and the finite element simulation model. This was because the mathematical model assumes a one-dimensional flux flow, i.e., the electrical and heat flux only flow along the applied field, whereas fluxes in other directions are not considered. These fluxes are caused by the material inhomogeneity and the shape feature of fibers. When the fiber volume fraction is very small or very large, homogeneity of the composite will inevitably be compromised. Moreover, when the direction of the applied potential is not perpendicular to the interface between the matrix and fibrous phase, the secondary dimensional flux will occur. Both situations above are the cause of deviation between the mathematical model and finite element simulation results.
Finite element models with tilted aligned fibers in a matrix were also constructed to study the cooling capacity of the composite. The properties of polycrystal Bi2Te3 and Copper were applied. The fiber rotation angle was fixed at 80 degrees. According to Figure 2, the composite can be seen as fiber arrays periodically aligned in the x-direction. Therefore, to simplify the finite element model, only one fiber array was constructed inside the matrix block, and periodic boundary conditions were applied on the surfaces parallel to the y-z plane. The geometry of the model was 2 mm × 500 mm × 10 mm, while the fiber diameter was adjusted to suit different fiber volume fractions. The bottom surface of the device was fixed at 300 K to serve as a heat sink, while the top surface was subjected to natural convective heat transfer. Two 1 mm thick Copper blocks were added on both ends of the y-direction of the composite to serve as electrodes. Twenty cases were studied with different fiber volume fractions, fiber shapes, and with respect to the two configurations mentioned above. Within each case, the input electrical current density was adjusted until the maximum temperature difference between the top and bottom surface was reached.
Figure 6 provides temperature distributions on the y–z plane for both cylindrical fiber model and square fiber model at 30% fiber volume fraction. Enlarged views of the top surface are shown in the insets. For both composites in Figure 6, a temperature gradient caused by the transverse Peltier effect can be clearly observed in the vertical direction. In the horizontal direction, the temperature was uniformly distributed except at the ends of the electrodes. In the insets of Figure 6, the abrupt temperature change at the interfaces of the fiber and matrix phases was the result of the Peltier effect.
The surface temperature was averaged on the top surface of each finite element model. The maximum temperature difference in each case was calculated by subtracting the heat sink temperature by the averaged cooling surface temperature. The results are presented in Table 2. In general, the mathematical model results agreed well with the finite element simulation results with at most a 10% difference. These small discrepancies were mainly caused by the secondary dimensional flux that was previously discussed. The only exception in Table 2 was the case for the cylindrical fiber with f = 0.1 under configuration I, where the ∆Tmax calculated using the mathematical model was more than 10% higher than the finite element simulation result. This was because under small or large fiber volume fraction cases, the homogeneity of the composite was very poor, which conflicts with the assumption in the mathematical model. The differences between the two configurations also agreed with the results shown in Figure 3, which emphasized the significant effect of the component material’s local property axes alignment in the composite coordinate system. During the finite element analysis, it was found that the size of the fibers has a slight influence on the TE performance. In anisotropic thermoelectric composites, smaller fibers provide the composite with better homogeneity. Although, under the same volume fraction of fibers the cooling capacity of the composite may not vary much with the size of the fibers, large fiber geometry is likely to cause higher temperature fluctuations in between the boundaries of the fibrous phase and the matrix phase.

4. Cooling Capacity Comparison for 1-Dimensional and 2-Demensional Inclusion Thermoelectric Composites with Anisotropic Components

The cooling performance of fibrous transverse thermoelectric composites were investigated with both a mathematical approach and a simulation approach. While the fibrous transverse thermoelectric composite can be treated as a 1D inclusion composite, the layered transverse thermoelectric composite can be treated as a 2D inclusion composite. The cooling capacity of the layered transverse thermoelectric composites with anisotropic components have been investigated thoroughly in previous studies [21]. Hence, it is informative to compare the cooling performances of transverse thermoelectric composites between 1D and 2D inclusion composites. In this section, three types of anisotropic thermoelectric crystals were chosen as the semiconducting phase of the composite, while Copper was chosen as the conducting phase of the composite. The material properties are listed in Table 1. The thermoelectric properties of Bi2Te3 single crystal were measured at 300 K, the thermoelectric properties of In4Se2.25 single crystal were measured at 600 K, and the thermoelectric properties of SnSe single crystal were measured at 700 K. The Copper served as the fibrous phase in the fibrous composite thermoelectric materials; its thermoelectric properties at 600 K and 700 K are also given in Table 1. It should be noticed that, in transverse thermoelectric composites, the transverse ZT matrix is anisotropic. The study made use of engineered anisotropy in composite materials and optimized the transverse ZT with respect to one specific element in the anisotropic ZT matrix. The anisotropic ZT values for other elements in the matrix are not discussed.
The semiconducting materials used in this comparison, i.e., Bi2Te3, In4Se2.25, and SnSe, all belong to planar isotropic material. Therefore, for both 1D and 2D inclusion composites consisting of semiconducting crystal and isotropic Copper, two representative configurations were used. Assuming the local material property axis as uiviwi (where the subscript i differentiates the type of materials), and the isotropy plane of material lies in the uivi plane. For fibrous composite materials, the definition of the two configurations were explained in previous sections of this study. For layered composite materials, in configuration III, the uivi plane aligns parallel to the layered plane of materials; in configuration IV, the uivi plane aligns perpendicular to the planar surface of the layers. Besides the four configurations, the conducting properties of isotropic polycrystal Bi2Te3 and SnSe, listed in Table 1, can be calculated using Reuss and Voigt models. Since Bi2Te3 and SnSe single crystals have isotropic Seebeck properties, their isotropic polycrystal Seebeck properties are the same as anisotropic single crystals. The conducting properties of the isotropic polycrystals are the average of the results calculated by the Reuss and Voigt models, which is referred as ‘Average RV’ in Table 3, respectively.
The results in Table 3 present both maximum ZtransT as well as the maximum cooling capacity (∆Tmax) for composites with different anisotropic components and configurations using the derived mathematical model. It can be seen that, for composites which consist of the same combination of component materials but under different configurations, the alignment of component material’s property axis in the composite has a major influence on ZtransT as well as the ∆Tmax values. For Bi2Te3/Copper composite, the difference in thermoelectric performances caused by the alignment of the anisotropic component’s property axis was 50%. For the SnSe/Copper composite, this difference exceeded 400%. This difference was caused by the variance of properties among different material property axes in the anisotropic component phase, where larger variance usually leads to bigger variance in the ZtransT and ∆Tmax values of a composite.
When the cooling performance of the transverse thermoelectric was compared under the same component combination but with different inclusion types, it was found that for the composite with the same component volume fraction, the fibrous and layered composites yielded very similar ∆Tmax as well as ZtransT values. However, it should be noted that the peak performances of the layered and fibrous composites under the same volume fraction were not identical to each other. The ∆Tmax values were rounded up to the nearest integer, and the maximum ZtransT values were rounded up to the second decimal in Table 3. Although according to Table 3, the ZtransT and ∆Tmax values in between different configurations have same values, it should be noted that none of the results in Table 3 are exactly identical to one another.
By comparing the results between the composite with anisotropic component material and isotropic component material, it was found that the anisotropic component in single-crystal form in both the layered transverse thermoelectric composite and the fibrous transverse thermoelectric composite can provide a better cooling performance compared to their isotropic polycrystal counterpart. For the Bi2Te3/Copper composite at 300 K, the composite with the anisotropic Bi2Te3 single crystal improved the maximum ZtransT and ∆Tmax by as much as 15% compared to the composite with the isotropic Bi2Te3 polycrystal. For the SnSe/Copper composite at 700 K, the improvement in the maximum ZtransT and ∆Tmax for composite with SnSe single crystal component was as much as 51% compared to composite with isotropic SnSe polycrystal components. For the In4Se2.25/Copper composite, the single-crystal form of In4Se2.25 had anisotropic Seebeck properties, but there are no existing studies regarding the effective properties of polycrystals with anisotropic Seebeck properties. Therefore, the comparison on the cooling performances for the In4Se2.25/Copper composite between composites with isotropic and anisotropic In4Se2.25 phase properties was not presented.

5. Conclusions

This study investigated the transverse thermoelectric properties of fibrous composites with anisotropic component materials. A mathematical model was built for predicting the transverse thermoelectric figure of merit (ZtransT) and maximum cooling capacity (∆Tmax) of the composite. The effectiveness of the mathematical model was verified by finite element simulations. The agreement among the two approaches indicated that the mathematical model can serve as an efficient tool for selecting and screening prospecting candidates for fibrous transverse thermoelectric composites with anisotropic material components. A case study using anisotropic Bi2Te3 polycrystal as the matrix and isotropic Copper as the fiber was performed and a ∆Tmax of 34 K was reached.
The variance in the maximum ZtransT and ∆Tmax values for both the layered and fibrous transverse thermoelectric composites was investigated with respect to the possible anisotropic properties of component phase. The results showed that by enhancing the anisotropic profile of polycrystal, the maximum ZtransT and ∆Tmax for both the layered and fibrous transverse thermoelectric composites can be improved. For the SnSe/Copper composite at 700 K, the improvement of maximum ZtransT and ∆Tmax values for the composite with anisotropic SnSe single-crystal components can be as much as 51% compared to the composites with isotropic SnSe polycrystals.
This study also showed that the alignment of the component’s local property axes in the composite can lead to significant variance in the ZtransT and ∆Tmax values of the composite. For the Bi2Te3/Copper composite, the difference in the thermoelectric performance caused by the alignment of the anisotropic component’s property axis was 50%. For the SnSe/Copper composite, this difference exceeded 400%. This discrepancy was mainly caused by the variance in the anisotropic properties among different material property axis directions in the component phase.

Author Contributions

Conceptualization, B.Q. and F.R.; methodology, B.Q.; software, F.W.; validation, B.Q., Y.Z. and T.W.; writing—original draft preparation, B.Q.; writing—review and editing, F.R., Y.Z. and T.W.; funding acquisition, B.Q. and F.R.

Funding

This research was funded by the Central South University and Temple University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ren, P.; Liu, Y.; He, J.; Lv, T.; Gao, J.; Xu, G. Recent advances in inorganic material thermoelectrics. Inorg. Chem. Front. 2018, 5, 2380–2398. [Google Scholar] [CrossRef]
  2. Twaha, S.; Zhu, J.; Yan, Y.; Li, B. A comprehensive review of thermoelectric technology: Materials, applications, modelling and performance improvement. Renew. Sustain. Energy Rev. 2016, 65, 698–726. [Google Scholar] [CrossRef]
  3. Goldsmid, H.J. Introduction to Thermoelectricity; Springer: Berlin, Germany, 2016. [Google Scholar]
  4. He, J.; Tritt, T.M. Advances in thermoelectric materials research: Looking back and moving forward. Science 2017, 357, eaak9997. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Yan, X.; Poudel, B.; Ma, Y.; Liu, W.; Joshi, G.; Wang, H.; Lan, Y.; Wang, D.; Chen, G.; Ren, Z. Experimental studies on anisotropic thermoelectric properties and structures of n-type Bi2Te2. 7Se0. 3. Nano Lett. 2010, 10, 3373–3378. [Google Scholar] [CrossRef] [PubMed]
  6. Wu, J.; Sun, Y.; Pei, W.B.; Huang, L.; Xu, W.; Zhang, Q. Polypyrrole nanotube film for flexible thermoelectric application. Synth. Met. 2014, 196, 173–177. [Google Scholar] [CrossRef]
  7. Yao, C.J.; Zhang, H.L.; Zhang, Q. Recent progress in thermoelectric materials based on conjugated polymers. Polymers 2019, 11, 107. [Google Scholar] [CrossRef] [PubMed]
  8. Du, Y.; Xu, J.; Paul, B.; Eklund, P. Flexible thermoelectric materials and devices. Appl. Mater. Today 2018, 12, 366–388. [Google Scholar] [CrossRef]
  9. Goldsmid, H.J. Application of the transverse thermoelectric effects. J. Electron. Mater. 2011, 40, 1254–1259. [Google Scholar] [CrossRef]
  10. Babin, V.P.; Gudkin, T.S.; Dashevskii, Z.M.; Dudkin, L.D.; Iordanishvili, E.K.; Kaidanov, V.I.; Kolomoets, N.V.; Narva, O.M.; Stil’bans, L.S. Artificial anisotropic thermoelements and their ultimate opportunities. Sov. Phys. Semicond. 1974, 8, 478. [Google Scholar]
  11. Sakai, A.; Kanno, T.; Takahashi, K.; Tamaki, H.; Adachi, H.; Yamada, Y. Enhancement in performance of the tubular thermoelectric generator (TTEG). J. Electron. Mater. 2013, 42, 1612–1616. [Google Scholar] [CrossRef]
  12. Kyarad, A.; Lengfellner, H. Transverse peltier effect in tilted Pb–Bi2Te3 multilayer structures. Appl. Phys. Lett. 2006, 89, 192103. [Google Scholar] [CrossRef]
  13. Kyarad, A.; Lengfellner, H. Al-Si multilayers: A synthetic material with large thermoelectric anisotropy. Appl. Phys. Lett. 2004, 85, 5613–5615. [Google Scholar] [CrossRef]
  14. Reitmaier, C.; Walther, F.; Lengfellner, H. Transverse thermoelectric devices. Appl. Phys. A 2010, 99, 717–722. [Google Scholar] [CrossRef]
  15. Kanno, T.; Sakai, A.; Takahashi, K.; Omote, A.; Adachi, H.; Yamada, Y. Tailoring effective thermoelectric tensors and high-density power generation in a tubular Bi0.5Sb1.5Te3/Ni composite with cylindrical anisotropy. Appl. Phys. Lett. 2012, 101, 011906. [Google Scholar] [CrossRef]
  16. Qian, B.; Ren, F. Transverse thermoelectricity in fibrous composite materials. Energies 2017, 10, 1006. [Google Scholar]
  17. Kaibe, H.; Tanaka, Y.; Sakata, M.; Nishida, I. Anisotropic galvanomagnetic and thermoelectric properties of n-type Bi2Te3 single crystal with the composition of a useful thermoelectric cooling material. J. Phys. Chem. Solids 1989, 50, 945–950. [Google Scholar] [CrossRef]
  18. Zhao, L.D.; Lo, S.H.; Zhang, Y.; Sun, H.; Tan, G.; Uher, C.; Wolverton, C.; Dravid, V.P.; Kanatzidis, M.G. Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals. Nature 2014, 508, 373–377. [Google Scholar] [CrossRef]
  19. Wei, Q.; Mukaida, M.; Kirihara, K.; Naitoh, Y.; Ishida, T. Recent progress on PEDOT-based thermoelectric materials. Materials 2015, 8, 732–750. [Google Scholar] [CrossRef] [PubMed]
  20. Wei, Q.; Mukaida, M.; Kirihara, K.; Ishida, T. Experimental studies on the anisotropic thermoelectric properties of conducting polymer films. ACS Macro Lett. 2014, 3, 948–952. [Google Scholar] [CrossRef]
  21. Qian, B.; Zhao, Y.; Ren, F. Effect of material anisotropy on the transverse thermoelectricity of layered composites. Int. J. Energy Res. 2018, 43, 181–188. [Google Scholar] [CrossRef] [Green Version]
  22. Nan, C.W. Physics of inhomogeneous inorganic materials. Prog. Mater. Sci. 1993, 37, 1–116. [Google Scholar] [CrossRef]
  23. Hone, J.; Llaguno, M.C.; Nemes, N.M.; Johnson, A.T.; Fischer, J.E.; Walters, D.A.; Casavant, M.J.; Schmidt, J.; Smalley, R.E. Electrical and thermal transport properties of magnetically aligned single wall carbon nanotube films. Appl. Phys. Lett. 2000, 77, 666–668. [Google Scholar] [CrossRef]
  24. Deng, F.; Zheng, Q.S.; Wang, L.F.; Nan, C.W. Effects of anisotropy, aspect ratio, and nonstraightness of carbon nanotubes on thermal conductivity of carbon nanotube composites. Appl. Phys. Lett. 2007, 90, 021914. [Google Scholar] [CrossRef]
  25. MacDonald, D.K.C. Thermoelectricity: An Introduction to The Principles; Courier Corporation: New York, NY, USA, 2006. [Google Scholar]
  26. Nye, J.F. Physical Properties of Crystals. Mater. Today 2007, 10, 391–397. [Google Scholar]
  27. Goldsmid, H.J. Porous materials in synthetic transverse thermoelements. J. Electron. Mater. 2010, 39, 1987–1989. [Google Scholar] [CrossRef]
  28. Qian, B.; Ren, F. Cooling performance of transverse thermoelectric devices. Int. J. Heat Mass Transf. 2016, 95, 787–794. [Google Scholar] [CrossRef]
  29. Rhyee, J.S.; Cho, E.; Lee, K.H.; Lee, S.M.; Kim, S.I.; Kim, H.-S.; Kwon, Y.S.; Kim, S.J. Thermoelectric properties and anisotropic electronic band structure on the In4Se3−x compounds. Appl. Phys. Lett. 2009, 95, 212106. [Google Scholar] [CrossRef]
  30. Burkov, A.T.; Vedernikov, M.V. Thermoelectric properties of metallic materials. In CRC Handbook of Thermoelectrics; CRC Press: Boca Raton, FL, USA, 1995. [Google Scholar]
  31. ASM Handbook, Volume 2—Properties of Pure Metals, Properties and Selection: Nonferrous Alloys and Special-Purpose Materials; ASM Handbook: Material Park, OH, USA, 1990; pp. 1099–1201.
  32. Voigt, W. Ueber die beziehung zwischen den beiden elasticitätsconstanten isotroper körper. Ann. Phys. Berlin 1889, 274, 573–587. [Google Scholar] [CrossRef]
  33. Willemse, M.W.M.; Caspers, W.J. Electrical conductivity of polycrystalline materials. J. Math. Phys. 1979, 20, 1824–1831. [Google Scholar] [CrossRef]
  34. Reuss, A. Berechnung der fließgrenze von mischkristallen auf grund der plastizitätsbedingung für einkristalle. ZAMM J. Appl. Math. Mech./Z. Angew. Math. Mech. 1929, 9, 49–58. [Google Scholar] [CrossRef]
  35. Kim, Y.M.; Lydia, R.; Kim, J.-H.; Lin, C.-C.; Ahn, K.; Rhyee, J.-S. Enhancement of thermoelectric properties in liquid-phase sintered Te-excess bismuth antimony tellurides prepared by hot-press sintering. Acta Mater. 2017, 135, 297–303. [Google Scholar] [CrossRef]
  36. Chen, C.L.; Wang, H.; Chen, Y.-Y.; Day, T.; Snyder, G.J. Thermoelectric properties of p-type polycrystalline SnSe doped with Ag. J. Mater. Chem. A 2014, 2, 11171–11176. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of the unit cell model in the fibrous thermoelectric composite which consisted of a fibrous (F) and matrix (M) phase. (a) Material property coordinate system of the fibrous phase; (b) material property coordinate system of the matrix phase; (c) schematic of the unit cell structure.
Figure 1. Schematic of the unit cell model in the fibrous thermoelectric composite which consisted of a fibrous (F) and matrix (M) phase. (a) Material property coordinate system of the fibrous phase; (b) material property coordinate system of the matrix phase; (c) schematic of the unit cell structure.
Materials 12 02049 g001
Figure 2. Schematic of a fibrous transverse thermoelectric composite with tilted fibers.
Figure 2. Schematic of a fibrous transverse thermoelectric composite with tilted fibers.
Materials 12 02049 g002
Figure 3. Maximum cooling capacity (∆Tmax) of the Bi2Te3 matrix/Cu fiber composite at 300 K under different material local coordinate system alignment configurations. (a) Configuration I; (b) Configuration II.
Figure 3. Maximum cooling capacity (∆Tmax) of the Bi2Te3 matrix/Cu fiber composite at 300 K under different material local coordinate system alignment configurations. (a) Configuration I; (b) Configuration II.
Materials 12 02049 g003
Figure 4. Maximum cooling capacity (∆Tmax) of Copper fiber/Bi2Te3 matrix transverse thermoelectric composite as a function of material anisotropy ratio ( r λ ,   r ρ ) of the Bi2Te3 phase. The aggregated values of the Bi2Te3 phase were calculated from single-crystal data (Table 1) using the (a) Reuss model and (b) Voigt model.
Figure 4. Maximum cooling capacity (∆Tmax) of Copper fiber/Bi2Te3 matrix transverse thermoelectric composite as a function of material anisotropy ratio ( r λ ,   r ρ ) of the Bi2Te3 phase. The aggregated values of the Bi2Te3 phase were calculated from single-crystal data (Table 1) using the (a) Reuss model and (b) Voigt model.
Materials 12 02049 g004
Figure 5. Finite element analysis results on effective properties of the unit cell (Figure 1c) in directions parallel and perpendicular to the fiber axis. (a) Thermal conductivity-Perpendicular to fiber axis; (b) Thermal conductivity-Along fiber axis; (c) Electrical Resistivity-Perpendicular to fiber axis; (d) Electrical Resistivity -Along fiber axis; (e) Seebeck coefficient-Perpendicular to fiber axis; (f) Seebeck coefficient-Along fiber axis.
Figure 5. Finite element analysis results on effective properties of the unit cell (Figure 1c) in directions parallel and perpendicular to the fiber axis. (a) Thermal conductivity-Perpendicular to fiber axis; (b) Thermal conductivity-Along fiber axis; (c) Electrical Resistivity-Perpendicular to fiber axis; (d) Electrical Resistivity -Along fiber axis; (e) Seebeck coefficient-Perpendicular to fiber axis; (f) Seebeck coefficient-Along fiber axis.
Materials 12 02049 g005aMaterials 12 02049 g005b
Figure 6. Temperature contour of fibrous transverse thermoelectric under cooling application. Upper: composite with square fibers. Lower: Composite with cylindrical fibers. Insect: Enlarged view of the composite’s cooling surface.
Figure 6. Temperature contour of fibrous transverse thermoelectric under cooling application. Upper: composite with square fibers. Lower: Composite with cylindrical fibers. Insect: Enlarged view of the composite’s cooling surface.
Materials 12 02049 g006
Table 1. Material properties used in this study.
Table 1. Material properties used in this study.
MaterialSingle Crystal
Bi2Te3
Polycrystal
Bi2Te3
Single Crystal
In4Se2.25
Single Crystal
SnSe
CopperCopperCopper
Measuring Temperature300 K300 K600 K700 K300 K600 K700 K
Sw (μV/K)−210−187−375−5402.833.343.84
Su/Sv (μV/K)−313
ρw (Ω·m)4.5 × 10−52.4 × 10−52.0 × 10−41.0 × 10−21.7 × 10−84.0 × 10−85.0 × 10−8
ρu\ρv (Ω·m)1.5 × 10−51.0 × 10−51.0 × 10−41.1 × 10−3
λw (W·m−1·K−1)1.000.780.800.25400386377
λu\λv (W·m−1·K−1)1.551.171.150.35
Reference[17][5][29][18][30,31][30,31][30,31]
T: Temperature; S: Seebeck coefficient; ρ: electrical resistivity; λ: thermal conductivity.
Table 2. Comparison of ∆Tmax between the mathematical model and the finite element simulation model.
Table 2. Comparison of ∆Tmax between the mathematical model and the finite element simulation model.
Volume Fraction of CopperConfiguration IConfiguration II
Mathematical Model
(K)
Square Fiber
(K)
Cylindrical Fiber
(K)
Mathematical Model
(K)
Square Fiber
(K)
Cylindrical Fiber
(K)
0.132.5929.9525.6723.2820.7421.37
0.230.6328.7628.3919.6919.3118.91
0.329.0326.5927.3617.9318.0217.86
0.428.4026.4127.2517.3318.4118.22
0.528.5728.1627.9817.5519.0419.22
Table 3. ∆Tmax and maximum ZtransT for different types of transverse thermoelectric composites.
Table 3. ∆Tmax and maximum ZtransT for different types of transverse thermoelectric composites.
Component 1Bi2Te3
Single Crystal
In4Se2.25
Single Crystal
SnSe
Single Crystal
Component 2Copper
Operating Temperature300 K600 K700 K
-∆TmaxZtransT∆TmaxZtransT∆TmaxZtransT
Fibrous
Composite
Configuration I300.24790.35920.35
Configuration II200.16650.28220.07
Average RV (Reuss and Voigt models)260.21N/AN/A610.21
Layered
Composite
Configuration III200.16650.28220.07
Configuration IV300.24800.35920.35
Average RV260.21N/AN/A610.21

Share and Cite

MDPI and ACS Style

Qian, B.; Ren, F.; Zhao, Y.; Wu, F.; Wang, T. Enhanced Thermoelectric Cooling through Introduction of Material Anisotropy in Transverse Thermoelectric Composites. Materials 2019, 12, 2049. https://doi.org/10.3390/ma12132049

AMA Style

Qian B, Ren F, Zhao Y, Wu F, Wang T. Enhanced Thermoelectric Cooling through Introduction of Material Anisotropy in Transverse Thermoelectric Composites. Materials. 2019; 12(13):2049. https://doi.org/10.3390/ma12132049

Chicago/Turabian Style

Qian, Bosen, Fei Ren, Yao Zhao, Fan Wu, and Tiantian Wang. 2019. "Enhanced Thermoelectric Cooling through Introduction of Material Anisotropy in Transverse Thermoelectric Composites" Materials 12, no. 13: 2049. https://doi.org/10.3390/ma12132049

APA Style

Qian, B., Ren, F., Zhao, Y., Wu, F., & Wang, T. (2019). Enhanced Thermoelectric Cooling through Introduction of Material Anisotropy in Transverse Thermoelectric Composites. Materials, 12(13), 2049. https://doi.org/10.3390/ma12132049

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