materials
Article Enhanced Thermoelectric Cooling through Introduction of Material Anisotropy in Transverse Thermoelectric Composites
Bosen Qian 1,* , Fei Ren 2, Yao Zhao 2, Fan Wu 3 and Tiantian Wang 4 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 * Correspondence: [email protected]
Received: 2 June 2019; Accepted: 21 June 2019; Published: 26 June 2019
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.
Keywords: thermoelectric cooling; transverse thermoelectricity; figure of merit; composite material; material anisotropy
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
Materials 2019, 12, 2049; doi:10.3390/ma12132049 www.mdpi.com/journal/materials Materials 2019, 12, 2049 2 of 13 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–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–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 Figure1, 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. Materials 2019, 12, x FOR PEER REVIEW 3 of 14 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 alignMaterials parallel2019, 12 to, 2049 each other, and a, w, z1, and z2 align parallel to each other. The cross-sectional plane3 of 13 of fiber (p–p) is parallel with the u–v plane of the matrix material.
(a) (b)
(c)
Figure 1. SchematicSchematic of of the the unit unit cell cell model model in in the the fibr fibrousous thermoelectric thermoelectric composite composite which which consisted consisted of aof fibrous a fibrous (F) (F)and and matrix matrix (M) (M) phase. phase. (a) (Materiala) Material property property coordinate coordinate system system of the of thefibrous fibrous phase; phase; (b) (materialb) material property property coordinate coordinate system system of the of thematrix matrix phase; phase; (c) (schematicc) schematic of the ofthe unit unit cell cell structure. structure.
TheWhen effective anisotropic properties material of propertiesthe unit cell are can introduced be derived into in thetwo component steps as illustrated phase (i.e., in fibrousFigure phase1c. In theand first matrix step, phase) the square of a compositefiber (F) is co material,mbined thewith alignment part of the of matrix the material phase (M local1) to propertyform a rectangular axis with compositerespect to theC1, compositewhere the material matrix block coordinate has the system same plays width a keyand roleheight in theas etheffective fibrous properties phase. In of the secondcomposite. step, Many the composite 1D inclusions, C1 is such combined as carbon with nanotubes the rest [of23 ,the24], matrix have di tofferent form material the unit properties cell. The effectivealong and thermal perpendicular and electrical to the fiberconducting axis direction. proper Therefore,ties of the inunit this cell study, can we be used calculated the subscript using “p”Kirchhoff’s to represent law, theand propertiesthe effective on Seebeck the cross-sectional coefficient can plane be ofcalculated the fiber using and subscript Thevenin’s “a” theorem to represent [25]. propertiesThe material along properties fiber axis of (Figurethe composite1a). As forC1 can the matrixbe derived phase, according we used to the Equation material (1): local axis system “uvw” to describe the anisotropic material properties (Figure1b). There are four sets of coordinate 1 systems = in Figure1, which are the local material coordinate system of the fibrous and matrix phase, = = 1 the coordinate system of composite C1(F + M1), and the coordinate system of unit cell C2(F + M1 + M2). In Figure1, the coordinate axes u, x 1, and x2 align parallel to each other; v, y1, and y2 align parallel to 1 each other,= and a, w, z1 , and z 2 align= parallel to each other. = The cross-sectional plane of fiber (p–p)(1) is 1 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 Figure1c. In the 1 first step, the square fiber (F) is combined with part of the matrix phase (M1) to form a rectangular = = = 1 composite C 1, where the matrix block has the same width and height as the fibrous phase. In the secondThe step, symbols the composite ρ, λ, S C, 1andis combined f stand withfor electrical the rest of resistivity, the matrix tothermal form the conductivity, unit cell. The eSeebeckffective coefficient,thermal and and electrical fiber volume conducting fraction properties. The subscripts of the unit p, a, cell u,can v, w, be x calculatedi, yi, and zi using standKirchho for componentff’s law, material’sand the eff propertiesective Seebeck along coe eachfficient material’s can be local calculated coordinate using systems. Thevenin’s The theorem material [ 25properties]. The material of the propertiesunit cell can of be the derived composite according C1 can to be Equation derived according(2): to Equation (1):
nρMu +ρFp λFp λMu (1+n) SFp λMu +nSMu λFp ρx = λx = Sx = 1 1+n 1 nλFp +λMu 1 nλFp +λMu ρMv ρFp (1+n) λFp +nλMv SFp ρMv +nSMv λFp ρy = λy = Sy = (1) 1 ρMv +nρFp 1 1+n 1 nρFp +ρMv ρMw ρFa (1+n) λFa +nλMw SFa ρMw +nSMw λFa ρz = λz = Sz = 1 ρMw +nρFa 1 1+n 1 nρFa +ρMw
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 Materials 2019, 12, x FOR PEER REVIEW 4 of 14
1 1 Materials 2019 12 , , 2049 1 4 of 13 = = 1 = 1 1 properties along each material’s local coordinate systems. The material properties of the unit cell can be derived according to Equation (2): 1 1 ! = ! ! S + 1 S = x1 ρMu 1 Mu ρx1 = ρx ρM u 1 p √ f − (2) 1 1 1 ρx2 = λx2 = λx1 + 1 λMu f Sx2 1= ! √ f (ρMu ρx )+ρx √ f − 1 1 − 1 1 1 ρx + ρMu √ f − 1 ! 1 Sy λMv + 1 SMv λy p λy λ 1 √ f 1 = + = 1 Mv = 1 − ρy2 f ρy1 ρv ρMv λy2 ! ! S y2 1 ! (2) − 1 1 1√ f 1 λy +λM 1 λy +λM = √ f − 1 v √ f − 1 v = 1 = ! 1 ! ! 1 S ρ + 1 S ρ 1 z1 Mw Mw z1 ρz1 ρMw 1 p √ f − ρz2 = λz2 = λz1 + 1 λMw f Sz2 = ! √ f (ρMw ρz )+ρz √ f − 1 − 1 1 1 ρz +ρMw √ f − 1 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 In Figure1, the fiber axis aligns parallel to the z 2-axis of the composite. When fibers are tilted transformation.aligned in the composite,For example, the if e theffective fibers properties in Figure of1 are the rotated composite by an can angle be calculated of θ around through x2-axis matrix into the configuration in Figure 2, the effective properties of the composite can be calculated as [26]: transformation. For example, if the fibers in Figure1 are rotated by an angle of θ around x2-axis into the configuration in Figure2 , the e ffective properties00 of the composite can be calculated as [26]: 1 0 2 = Px 0 0 2 2 (3) 2 2 1 Pxyz = 0 Py1cos (θ) + Pz sin (θ) Pz Py sin(2θ) (3) 0 2 2 2 2 2 − 2 1 2 2 0 2 Pz Py sin(2θ) Pz cos (θ) + Py sin (θ) 2 2 − 2 2 2
Figure 2. Schematic of a fibrous transverse thermoelectric composite with tilted fibers. Figure 2. Schematic of a fibrous transverse thermoelectric composite with tilted fibers. In a transverse thermoelectric material, heat flux and electrical current flow perpendicular to each other.In Baseda transverse on Figure thermoelectric2, we assume material, the electrical heat currentflux and flows electrical in the current y-direction flow and perpendicular heat flux flows to 2 each other. Based on Figure 2, we assume the electrical current flows in the y-direction andSzy heatT flux in the z-direction. The transverse thermoelectric figure of merit is defined as Z T = [10 ,27], trans ρyyλzz flows in the z-direction. The transverse thermoelectric figure of merit is defined as Z T= where Szy is the transverse Seebeck coefficient, ρyy is electrical resistivity in the y-direction, λzz is [10,27],thermal where conductivity Szy is the intransverse the z-direction, Seebeck T coefficient, is the operating ρyy is electrical temperature. resistivity The in term the Sy-direction,zy, ρyy, and λλzzzz isvalues thermal can conductivity be calculated in using the z-direction, Equation (1) T tois Equationthe operating (3). temperature. The term Szy, ρyy, and λzz valuesFigure can be2 calculatedcan also represent using Equation a fibrous (1) transverse to Equation thermoelectric (3). device operating under cooling mode,Figure where 2 can the topalso surface represent serves a fibrous as a cooling transverse surface thermoelectric with a temperature device of operating Tc and the under bottom cooling surface mode,serves where as a heat the sinktop surface with a temperatureserves as a cooling of Th. When surface the with electrical a temperature current flows of Tc in and the the y-direction, bottom surfacethe transverse serves as Peltier a heat eff ectsink will with trigger a temperature heat flux along of Th the. When z-direction. the electrical The maximum current flows cooling in capacity the y- direction,of the system the transverse is related to Peltier the Ztrans effectT value, will accordingtrigger heat to previousflux along studies the z-direction. [10,28]. In thisThe study, maximum the Th coolingwas set capacity to be the of measuring the system temperature is related to of the material ZtransT value, properties, according which to are previous shown instudies Table 1[10,28].. In this study, the Th was set to be the measuring temperature of material properties, which are shown in Table 1.
Materials 2019, 12, 2049 5 of 13
Table 1. Material properties used in this study.
Single Crystal Polycrystal Single Crystal Single Crystal Material Copper Copper Copper Bi2Te3 Bi2Te3 In4Se2.25 SnSe Measuring 300 K 300 K 600 K 700 K 300 K 600 K 700 K Temperature Sw (µV/K) 375 210 187 − 540 2.83 3.34 3.84 Su/Sv (µV/K) − − 313 − 5 5 − 4 2 ρw (Ω m) 4.5 10− 2.4 10− 2.0 10− 1.0 10− 1.7 4.0 5.0 · × 5 × 5 × 4 × 3 ×8 ×8 ×8 ρu ρv (Ω m) 1.5 10− 1.0 10− 1.0 10− 1.1 10− 10− 10− 10− \ 1· 1 × × × × λw (W m− K− ) 1.00 0.78 0.80 0.25 · ·1 1 400 386 377 λu λv (W m K ) 1.55 1.17 1.15 0.35 \ · − · − Reference [17][5][29][18][30,31][30,31][30,31] T: Temperature; S: Seebeck coefficient; ρ: electrical resistivity; λ: thermal conductivity.
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 shownMaterials in2019 Table, 12, x1 .FOR According PEER REVIEW to Table 1, the Copper phase was isotropic and Bi 2Te3 was anisotropic7 of 14 with planar (u–v plane) isotropy. WhenWhen anisotropicanisotropic materialmaterial propertiesproperties areare implementedimplemented inin thethe componentcomponent phasephase ofof aa composite,composite, thethe alignmentalignment of of component component material’s material’s local local property property coordinate coordinate system, system, with respect with to respect the composite to the coordinatecomposite system,coordinate will system, affect the will eff ectiveaffect propertiesthe effective of properties the composite. of th Ine composite. Figure1, the In u–v Figure plane 1, ofthe the u– matrixv plane material of the alignsmatrix perpendicular material aligns to the perpendicular fiber axis. If theto materialthe fiber property axis. If coordinatethe material in Figureproperty1b coordinate in Figure 1b is rotated around the u-axis 90 degrees, the anisotropy plane of Bi2Te3 will be is rotated around the u-axis 90 degrees, the anisotropy plane of Bi2Te3 will be parallel to the fiber axis direction.parallel to For the the fiber convenience axis direction. of analysis, For the we shallconvenience refer to theof twoanalysis, configurations we shall mentionedrefer to the above two asconfigurations configuration mentioned I and configuration above as II. configuration In configuration I and I, theconfiguration material property II. In axesconfiguration u, v, and wI, arethe material property axes u, v, and w are parallel to the composite axes x2, y2, and z2, respectively. In parallel to the composite axes x2, y2, and z2, respectively. In configuration II, the material property configuration II, the material property axes u, v, and w are parallel to the composite axes x2, z2, and axes u, v, and w are parallel to the composite axes x2, z2, and y2, respectively. The corresponding y2, respectively. The corresponding ZtransT and maximum cooling capacities (∆Tmax) based on these ZtransT and maximum cooling capacities (∆Tmax) based on these two configurations were calculated usingtwo configurations the derived mathematical were calculated model using with respect the derived to fiber mathematical rotation angle, modelθ, and fiberwith volumerespect fraction,to fiber frotation. The results angle, are θ, shown and fiber in Figure volume3. fraction, f. The results are shown in Figure 3.
(a) (b)
FigureFigure 3.3. MaximumMaximum cooling cooling capacity capacity (∆ (T∆maxTmax) of) the of theBi2Te Bi3 2matrix/CuTe3 matrix fiber/Cu fibercomposite composite at 300 at K 300under K underdifferent diff erentmaterial material local localcoordinate coordinate system system alignment alignment configurations. configurations. (a) Configuration (a) Configuration I; (b) I; (Configurationb) Configuration II. II.
BasedBased onon FigureFigure3 ,3, the the Z ZtranstransT and ∆Tmax valuesvalues exhibited exhibited similar similar trends trends for both configurations. configurations. TheThe maximummaximum valuesvalues appearredappearred atat rotationrotation anglesangles betweenbetween 7575 andand 8585 degrees.degrees. TheThe changechange inin fiberfiber 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–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