許鵬飛 唐 豪

?

新型結(jié)構(gòu)熱電元件的多物理場耦合研究

許鵬飛 唐 豪

（南京航空航天大學(xué)能源與動力學(xué)院 南京 210016）

為了改善熱電元件的工作性能，提出了兩種新型結(jié)構(gòu)的熱電元件，分別是半傾斜型和雙傾斜型，傾斜角度57.99°。參考傳統(tǒng)的垂直型熱電元件，采用數(shù)值模擬的方法，通過控制電流、熱端溫度和冷熱端溫差等參量，對三種結(jié)構(gòu)的熱電元件進行了溫度場、電場耦合分析。結(jié)果顯示：雙傾斜型電熱元件能夠?qū)崿F(xiàn)最低的冷端溫度，且隨著電流的增大，不同結(jié)構(gòu)冷端溫度的差值先增大后減?。辉谛」β蕳l件下，雙傾斜型電熱元件制冷量最大；雙傾斜結(jié)構(gòu)能夠達到最高的制冷系數(shù)，但在大功率條件下制冷系數(shù)卻最低；傾角為鈍角處存在場強突變現(xiàn)象；最大電流密度均出現(xiàn)在不與半導(dǎo)體接觸的銅片區(qū)域中；雙傾斜結(jié)構(gòu)更適用于微小型低功耗熱電設(shè)備。

熱電元件；傾斜；溫度場；電場

0 引言

熱電制冷，又稱溫差電制冷，是一種基于帕爾貼效應(yīng)的主動制冷方式，具有重要的應(yīng)用價值和廣闊的發(fā)展前景。其原理如圖1所示，電源接通后，電流依次通過銅片、N型半導(dǎo)體、銅片、P型半導(dǎo)體和銅片，由于空穴和電子在半導(dǎo)體材料內(nèi)部發(fā)生遷移，在銅片和半導(dǎo)體的接觸面處產(chǎn)生冷熱效應(yīng)，一端制冷，溫度降低，另一端制熱，溫度升高，最終形成溫度梯度。因此，熱電制冷的實質(zhì)是溫度場和電場的耦合效應(yīng)[1]。熱電制冷同時具有明顯的優(yōu)點及缺點。

圖1 熱電制冷原理圖

相比于機械式制冷，熱電制冷優(yōu)點顯著，其無需制冷介質(zhì)，潔凈無污染，無需外界施加機械動力，可將發(fā)熱器件的工作溫度降低至環(huán)境溫度以下，無振動、無噪聲，制冷器具有模塊化結(jié)構(gòu)，質(zhì)量輕、體積小、使用壽命長。然而，熱電制冷的缺點也很明顯，制冷量小，制冷效率偏低，為了達到更好的制冷效果，熱電制冷器本身需要很大的電流，熱耗較大。

近年來，國內(nèi)外學(xué)者針對熱電制冷從不同角度進行了一系列深入研究，取得了碩果累累的成就。Perez-Aparicio[2]等人對三維、非線性、變物性參數(shù)的熱電元件進行了有限元模擬，給出了電勢、溫度和熱流量的分布狀況，所得結(jié)果與實驗及理論計算非常吻合，這對于優(yōu)化變物性參數(shù)熱電設(shè)備具有重要指導(dǎo)意義。毛佳妮[3]等人針對兩級熱電制冷器，在第三類邊界條件下，建立了冷熱端介質(zhì)環(huán)境擾動的物理模型，推導(dǎo)得到了兩級熱電組件的關(guān)鍵參數(shù)表達式，為多級熱電制冷器內(nèi)部結(jié)構(gòu)優(yōu)化提供了理論支持。李曉[4]等人分析了熱電制冷的物理模型，討論了不同因素對制冷性能的影響，為熱電空調(diào)的應(yīng)用提供了理論依據(jù)。杜芳莉[5]等人提出將熱電制冷技術(shù)應(yīng)用于筆記本散熱，并且有效回收利用制冷器熱端產(chǎn)生的熱量。王振雨[6]等人針對傳統(tǒng)風(fēng)冷式熱電空調(diào)器厚度較大的缺點，采用熱管橫向引熱的方法，提出了一種新型熱電空調(diào)器，通過數(shù)值模擬對比發(fā)現(xiàn)，改進后厚度減小的熱電空調(diào)器仍能保持原有的制冷性能，因此更加適用于超薄微小空間。Oliveira[7]等人提出了熱電元件二維模型，采用有限元法和TDMA算法推導(dǎo)出數(shù)學(xué)控制方程，并將模型計算結(jié)果與實驗結(jié)果作對比，發(fā)現(xiàn)誤差在10%以內(nèi)。

本文提出了兩種新型傾斜結(jié)構(gòu)的熱電元件，同時以常規(guī)的垂直型結(jié)構(gòu)作參考，基于ANSYS Workbench Thermal-Electric模塊分別對三種熱電元件進行多物理場耦合研究，分析探討了傾斜結(jié)構(gòu)元件的制冷性能，闡明傾斜結(jié)構(gòu)的性能特點，專門為微小型、小功耗熱電設(shè)備的性能優(yōu)化提供了新的思路和借鑒。

1 熱電元件分析模型

1.1 幾何模型

本文涉及三種不同結(jié)構(gòu)形式的熱電元件，如圖2所示，為別為（a）垂直型，即常規(guī)結(jié)構(gòu)[8]；（b）半傾斜型；（c）雙傾斜型。

圖2 熱電元件幾何模型

三種模型垂直于紙面的厚度都是5mm。圖2（a）垂直型即傳統(tǒng)常規(guī)熱電元件，同時也是其他兩種類型的參考結(jié)構(gòu)。圖2（b）半傾斜型，P型半導(dǎo)體保持不變，N型半導(dǎo)體設(shè)計成傾斜角度為57.99°的平行六面體結(jié)構(gòu)，但橫截面積和垂直高度與P型半導(dǎo)體相同。圖2（c）雙傾斜型，P型、N型半導(dǎo)體均設(shè)計成傾斜角度為57.99°的平行六面體結(jié)構(gòu)，同時橫截面積和垂直高度不變。三種熱電元件結(jié)構(gòu)形式不同，但無論是P型半導(dǎo)體還是N型半導(dǎo)體，體積均不發(fā)生變化，即三種結(jié)構(gòu)形式下的半導(dǎo)體材料用量相同，因此三種結(jié)構(gòu)具有邏輯上的可比性[8]。

1.2 物理模型

為了問題的簡化和計算的方便，本文提出以下五點假設(shè)：（1）熱電元件在定常條件下運行；（2）湯姆森效應(yīng)在熱電制冷過程中是一種二級效應(yīng)，因而可將其忽略不計[9]；（3）忽略冷、熱端與銅片的接觸熱阻；（4）熱電元件側(cè)面絕熱；（5）熱電材料為常物性。因為本文著重研究三種結(jié)構(gòu)熱電元件的制冷性能對比，所以常物性并不影響熱電元件性能的相對優(yōu)劣。材料的物性參數(shù)參考文獻[10]，取=303K時的值，如表1所示。

表1 材料物性參數(shù)

1.3 控制方程[11]

在熱電制冷過程中，熱通量三維控制方程可表述為：

式中，為塞貝克系數(shù)，V/K；為熱導(dǎo)率，W/(m·K)。將（2）式帶入（1）式，可得：

式中，代表電阻率，Ω·m；綜合（4）—（6）式，可得：

2 網(wǎng)格獨立性檢驗

本文借助ANSYS Workbench中的mesh功能對三種模型劃分結(jié)構(gòu)化網(wǎng)格，然后進行網(wǎng)格獨立性檢驗。設(shè)定不同的最小網(wǎng)格尺寸，分別為1mm、0.5mm、0.2mm，設(shè)定電流值為2.5A，熱端溫度為45℃，其余邊界條件保持默認。表2是不同最小網(wǎng)格尺寸所對應(yīng)的冷端溫度和相對偏差。

表2 垂直型網(wǎng)格獨立性檢驗

表3 半傾斜型網(wǎng)格獨立性檢驗

表4 雙傾斜型網(wǎng)格獨立性檢驗

這里以最小網(wǎng)格尺寸0.2mm為參考值，當(dāng)最小網(wǎng)格尺寸為1mm時，三種模型中最大相對偏差為0.112%；當(dāng)最小網(wǎng)格尺寸為0.5mm時，三種模型中最大相對偏差為0.019%。由于相對偏差均小于1%，所以為了計算的方便，三種模型網(wǎng)格劃分均可設(shè)定最小網(wǎng)格尺寸為1mm。

3 多物理場耦合分析

3.1 冷端溫度變化規(guī)律

對三種熱電元件施加50℃的熱端溫度，通以大小相等的電流，分析對比熱電元件冷端溫度的變化規(guī)律，結(jié)果如圖3所示。

圖3 冷端溫度隨電流的變化關(guān)系

隨著電流的增大，冷端溫度均逐漸降低，同時溫度下降的速率也逐漸減小。冷端溫度由低到高依次是雙傾斜型、半傾斜型、垂直型。因此，雙傾斜型可以獲取更低的制冷溫度。進一步觀察圖3還可以發(fā)現(xiàn)，伴隨電流的增大，不同結(jié)構(gòu)熱電元件冷端溫度之差先升高后降低。例如，當(dāng)電流=1A時，相比于垂直型，半傾斜型擁有1.647℃的差值，而雙傾斜型則對應(yīng)有3.445℃的差值；當(dāng)=5A時，同樣相比于垂直型，兩項差值分別是4.7363℃和9.3949℃；當(dāng)=10A，兩項差值又分別是4.004℃和7.594℃。

3.2 制冷量隨電流、溫差的變化規(guī)律

先后探討不同結(jié)構(gòu)熱電元件的制冷量Q關(guān)于電流、溫差Δ的變化規(guī)律，同時進行歸納總結(jié)。

控制冷熱端溫差為定值，改變電流大小，得到不同結(jié)構(gòu)熱電元件的制冷量與電流的變化規(guī)律?，F(xiàn)以Δ=20℃為例，結(jié)果如圖4所示。

雖然熱電元件結(jié)構(gòu)不同，但熱電制冷量均呈現(xiàn)相同的變化趨勢。隨著電流的增大，制冷量先增大后減小，所以每條變化曲線都擁有一個極值點。垂直型、半傾斜型、雙傾斜型元件的最大制冷量分別是1.1208W、0.99164W、0.88379W，同時最大制冷量所對應(yīng)的電流值又分別是22A、19A、17A。因此，得出結(jié)論：雙傾斜型能夠取得的制冷量峰值最小，但消耗的功率也最小，為2.2064W；垂直型可取得最大的制冷量，但消耗的功率也最大，為2.9099W。從圖4還可以發(fā)現(xiàn)，當(dāng)電流小于5A時，制冷量由大到小分別是雙傾斜型、半傾斜型、垂直型。進一步得出結(jié)論：在小功率工況下，雙傾斜型性能最優(yōu)，大功率工況垂直型最優(yōu)。

圖4 制冷量隨電流的變化關(guān)系

控制電流為定值，改變冷熱端溫差，得到不同結(jié)構(gòu)熱電元件的制冷量Q與溫差Δ的變化曲線。以小功率=4.5A（＜5A）為例，結(jié)果如圖5所示。

圖5 制冷量隨溫差的變化關(guān)系

通過多項式擬合，制冷量與溫差近似成線性關(guān)系，其中垂直型斜率最大，雙傾斜型斜率最小。當(dāng)Δ＞10℃，制冷量由大到小依次是雙傾斜型、半傾斜型、垂直型，且隨著溫差的增大，不同結(jié)構(gòu)熱電元件的制冷量之差越來越大。當(dāng)Δ=10℃時，相比于垂直型，半傾斜型的制冷量差值為0.0035W，而雙傾斜型的制冷量差值為0.0053W；當(dāng)Δ=50℃時，兩者制冷量差值又分別是0.0468W和0.0875W。因此，在較高溫差的前提下，垂直型性能最差，而雙傾斜型性能最佳，最具優(yōu)勢，且溫差越大，其優(yōu)勢越明顯。

3.3 制冷系數(shù)分析

制冷系數(shù)是評價熱電制冷性能的一項重要指標(biāo)。本文繼續(xù)以溫差Δ=20℃為例，研究對比不同結(jié)構(gòu)熱電元件的制冷系數(shù)隨電流的變化情況，變化曲線如圖6所示。

隨著電流的增大，制冷系數(shù)先增大后減小，每條變化曲線都有一個極值。垂直型、半傾斜型、雙傾斜型的最高制冷系數(shù)分別是1.9764、1.9708和1.9765，對應(yīng)的電流值又分別是3.7A、3.3A和3A。從中可以發(fā)現(xiàn)，雙傾斜型最高制冷系數(shù)最大，而半傾斜型最低，但前者只比后者高了0.29%。同時，雙傾斜型達最高制冷系數(shù)所需功率最小，為0.0907W，而垂直型所需功率最大，為0.1097W。進一步觀察圖6還可以發(fā)現(xiàn)，當(dāng)電流低于3A時，雙傾斜型制冷系數(shù)最高，而垂直型最低；當(dāng)電流高于3A時，垂直型制冷系數(shù)最高，而雙傾斜型最低。因此，雙傾斜型適用于小功率工況，而垂直型適用于大功率工況。

圖6 制冷系數(shù)隨電流的變化關(guān)系

3.4 電場分析

ANSYS Workbench中的Thermal-Electric模塊是一種熱電耦合功能模塊，既可以模擬溫度場的分布特點，還可以實現(xiàn)電場的仿真，從而能夠便于多物理場耦合分析?，F(xiàn)針對三種結(jié)構(gòu)的熱電元件，均施加電流=3A、冷熱端溫差Δ=30℃的邊界條件，模擬得到電勢分布圖，如圖7所示。經(jīng)對比分析發(fā)現(xiàn)，三種結(jié)構(gòu)的電勢從P型半導(dǎo)體到N型半導(dǎo)體逐漸減小至零，但電勢分布都是二維變化，在Y軸方向上電勢保持不變。同時，銅片中的電勢值處處相等，分布均勻。對于垂直型元件，電勢在半導(dǎo)體中沿Z軸線性分布。

對于半傾斜型元件，在垂直的P型半導(dǎo)體中，電勢仍關(guān)于Z軸線性變化，在傾斜的N型半導(dǎo)體中，電勢沿棱邊所在方向線性變化，同時在傾角為鈍角處存在電勢梯度，在傾角為銳角處則不存在，其對場強分布造成的影響將在下文進行分析。雙傾斜型元件的電勢分布與半傾斜型元件呈現(xiàn)相同的規(guī)律。

下面分析電場強度規(guī)律，如圖8所示為電場強度矢量分布圖。

對比矢量圖可以發(fā)現(xiàn)，半導(dǎo)體區(qū)域內(nèi)電場強度較大，遠遠高于銅片中的電場強度。這主要是因為銅片電阻率小，僅為1.7×10-8Ω·m，而P、N型半導(dǎo)體則達到8.98×10-6Ω·m、8.36×10-6Ω·m，是前者的500倍左右，所以在圖7中，銅片電勢下降幅度近似于零，電勢等值分布，無法形成電勢梯度。對于垂直型元件，半導(dǎo)體內(nèi)電場強度均勻分布，大小1.78V/m左右，方向豎直向上或向下。這是因為半導(dǎo)體內(nèi)電勢在Z軸方向上一維線性變化，即圖7（a）所示，繼而形成均勻的電勢梯度。對于半傾斜型元件，如圖8（b）所示，P型半導(dǎo)體保持垂直，因而電場強度分布同垂直型元件，但傾斜的N型半導(dǎo)體，電場強度在上下傾斜角處陡然增大，最高可達4.59V/m，中間區(qū)域電場強度仍能保持均勻分布，但大小遠不如傾斜區(qū)域。雙傾斜型熱電元件也有類似的現(xiàn)象。如圖8（c）所示，P、N半導(dǎo)體上下傾角處電場強度陡然增大，最高可達4.82V/m，而中間區(qū)域電場強度分布均勻，大小2V/m左右。同時結(jié)合半傾斜型場強分布圖，還可發(fā)現(xiàn)，場強突變均發(fā)生在傾角為鈍角處，而銳角處無突變現(xiàn)象。分析原因，根本在于電極的傾斜外形結(jié)構(gòu)。借助圖7（c）可知，在鈍角傾角處電勢變化幅度大，因此電場強度較大，而在銳角傾角處，電勢均勻分布，電場強度近似為零。

最后討論電流密度分布規(guī)律。如圖9所示為電流密度分布圖，將三種結(jié)構(gòu)進行對比分析，可以發(fā)現(xiàn)，電流密度分布規(guī)律彼此相似，其中最大值均在6.5×105A/m2以上，且出現(xiàn)在不與半導(dǎo)體接觸的銅片區(qū)域內(nèi)。

三種結(jié)構(gòu)的最小電流密度均出現(xiàn)在銅片頂邊區(qū)域，但數(shù)值存在差異。垂直型、半傾斜型、雙傾斜型的最小值分別是78.59A/m2、78.58A/m2、27299A/m2。垂直型和半傾斜型由于都包含垂直電極，因此電流密度最小值幾乎相等，僅相差0.01A/m2。而雙傾斜型則完全不同于前兩者，當(dāng)電流通過銅片與半導(dǎo)體的接觸面時，由于存在傾角，電流方向并非90°彎折，造成拐角處電流密度較大。

圖7 三種熱電元件的電勢分布云圖

圖8 三種熱電元件的電場強度矢量圖

圖9 三種熱電元件的電流密度分布云圖

4 總結(jié)與展望

本文針對三種不同結(jié)構(gòu)的熱電元件，采用數(shù)值模擬的方法對其進行了多物理場耦合分析，所得結(jié)論主要有：

（1）雙傾斜型熱電元件能夠?qū)崿F(xiàn)最低的冷端溫度，不同結(jié)構(gòu)的冷端溫度差值隨電流的增大而先增大后減??；

（2）三種結(jié)構(gòu)中，雙傾斜型制冷量峰值最低，但在小功率條件下制冷量卻最大；

（3）雙傾斜型能夠達到最高的制冷系數(shù)，且消耗的功率最?。淮蠊β蕳l件下垂直型制冷系數(shù)最高；

（4）傾角為鈍角處存在電場強度突變，而傾角為銳角或垂直型結(jié)構(gòu)并無此現(xiàn)象；

（5）最大電流密度出現(xiàn)在不與半導(dǎo)體接觸的銅片區(qū)域中；

（6）小功率設(shè)備宜采用雙傾斜型結(jié)構(gòu)，而大功率設(shè)備采用垂直型結(jié)構(gòu)。

[1] 胡仁喜,閆波,康士廷,等.ANSYS15.0多物理耦合場有限元分析從入門到精通[M].北京:機械工業(yè)出版社,2015.

[2] J L Pérez-Aparicio, Palma R, Taylor R L. Finite element analysis and material sensitivity of Peltier thermoelectric cells coolers[J]. International Journal of Heat & Mass Transfer, 2012,55(4):1363-1374.

[3] 毛佳妮,王世飛,江述帆,等.第三類熱邊界條件下兩級熱電制冷器的性能評價及優(yōu)化[J].低溫工程,2016,(2):1-7.

[4] 李曉,孫寶紅,彭慶龍.半導(dǎo)體空調(diào)的理論研究與性能分析[J].制冷與空調(diào),2013,27(2):110-113,129.

[5] 杜芳莉,樊麗娟,白大雨,等.熱電制冷在筆記本電腦散熱系統(tǒng)中的應(yīng)用[J].制冷與空調(diào),2014,28(3):310-313.

[6] 王振雨,唐豪.風(fēng)冷式熱電空調(diào)器的改進研究[J].制冷與空調(diào),2016,16(8):18-22.

[7] Oliveira K S M, Cardoso R P, Hermes C J L. Numerical assessment of the thermodynamic performance of ther- moelectric cells via two-dimensional modelling[J].App- lied Energy,2014,130(5):280-288.

[8] 許鵬飛,唐豪.傾斜結(jié)構(gòu)熱電元件的多物理場耦合研[J].南京航空航天大學(xué)學(xué)報,2018,50(4):494-500.

[9] 徐德勝.半導(dǎo)體制冷與應(yīng)用技術(shù)[M].上海:上海交通大學(xué)出版社,1999.

[10] Wang X D, Huang Y X, Cheng C H, et al. A three-dimensional numerical modeling of thermoelectric device with consideration of coupling of temperature field and electric potential field[J]. Energy, 2012,47(1):488-497.

[11] Chen W H, Liao C Y, Hung C I. A numerical study on the performance of miniature thermoelectric cooler affected by Thomson effect[J]. Applied Energy, 2012,89(1):464-473.

Multiphysics Coupling Research on Thermoelectric Element with New Structure

Xu Pengfei Tang Hao

( College of Energy and Power Engineering, Nanjing University of Aeronautics & Astronautics, Nanjing, 210016 )

Two kinds of thermoelectric elements with new structures are proposed to improve the working performance of thermoelectric element. They are semi-sloping type and double-sloping type respectively. The slope angle is 57.99°. Taking the traditional vertical thermoelectric element as a reference, temperature field and electric field coupling research is carried on the three kinds of thermoelectric elements with different structures by means of numerical simulation and controlling current, hot side temperature and temperature difference between cold side and hot side. The results show that: Double-sloping type can achieve the lowest cold side temperature. The difference between cold side temperature of different structures increases first and decreases second with the increase of current. The cooling capacity of double-sloping type is the largest in the condition of small power. The double-sloping structure can achieve the highest, but makes the lowestin the condition of large power. The electric field intensity will change suddenly where the slope angle is obtuse. The maximum of current density occurs in copper region where copper does not touch with semiconductor. Double-sloping structure is more suitable for the thermoelectric equipment with small size and small power.

thermoelectric element; sloping structure; temperature field; electric field

1671-6612（2018）06-567-07

TN248.41

A

南京321計劃重點項目資助（2013B03004）

許鵬飛（1992-），男，碩士研究生，主要從事熱電制冷方面的研究，E-mail：1660494071@qq.com

唐 豪（1960-），男，教授，博士生導(dǎo)師，主要從事空調(diào)及低溫工程方面的研究，E-mail：hao.tang@nuaa.edu.cn

2018-01-17