孫悅桐, 孟凡凱, 周 林, 徐辰欣

(海軍工程大學(xué) 動(dòng)力工程學(xué)院,武漢 430033)

熱電制冷又稱溫差電制冷或半導(dǎo)體制冷,可通過帕爾貼效應(yīng)將電能直接轉(zhuǎn)化為溫差[1].與傳統(tǒng)的壓縮制冷相比,熱電制冷無運(yùn)動(dòng)部件,可靠性高,壽命長、易維護(hù);模塊化特點(diǎn)使其在制冷量增大時(shí)沒有顯著的性能損失;無工作流體,可避免泄漏危害環(huán)境[2-4].隨著熱電材料技術(shù)的迅速發(fā)展,熱電制冷技術(shù)成為制冷與低溫領(lǐng)域的研究熱點(diǎn),在微電子系統(tǒng)、激光通信、醫(yī)療等領(lǐng)域得到了廣泛應(yīng)用[5-7].

一些特殊領(lǐng)域的設(shè)備需要在低溫環(huán)境中維持運(yùn)行.如基于異質(zhì)結(jié)勢壘型結(jié)構(gòu)的中波紅外探測器工作溫度在150～195 K的范圍內(nèi)[8];各油庫和油站使用的測定儀器需要滿足300～200 K逐點(diǎn)降溫的要求;X射線太空望遠(yuǎn)鏡需要在低溫的狀態(tài)下工作,制冷溫度需達(dá)到170 K左右[9],進(jìn)而實(shí)現(xiàn)降低噪聲、減弱靈敏度的效果.單級熱電制冷器最大制冷溫差只能達(dá)到60 K左右[10],對于更大溫差、更大溫度跨度的場合,單級熱電制冷器無法滿足需求,此時(shí)需要采用多級制冷,通過增加級數(shù)獲得更大的制冷溫差.

多級熱電制冷即采用垂直級聯(lián)(金字塔狀)排列方式,低一級冷端為高一級熱端散熱,由高溫端至低溫端,制冷溫度逐漸降低,低溫級模塊可達(dá)到最低制冷溫度[11].隨著溫度跨度的增加,多級熱電制冷器的性能會(huì)迅速降低[12],對多熱熱電制冷器性能的分析與優(yōu)化具有重要意義.Karimi等[13]討論了金字塔型多級制冷器的熱流密度和性能系數(shù)最大值,數(shù)值模擬結(jié)果表明,熱端散熱器的熱阻是決定多級熱電制冷器整體性能的關(guān)鍵因素.Parashchuk等[14]開發(fā)了一種多級熱電制冷器,利用四級商用模塊和兩級超低溫模塊設(shè)計(jì)并制作了冷卻樣機(jī),工作溫度可達(dá)140 K.Phong和Shih[15]在建模結(jié)果與性能尺寸權(quán)衡的基礎(chǔ)上構(gòu)建了多級熱電制冷器的實(shí)驗(yàn)?zāi)Ｐ?實(shí)驗(yàn)結(jié)果顯示,在熱負(fù)荷為10 mW的情況下,六級熱電制冷器可達(dá)到的溫度差為137 K.Sofradir公司運(yùn)用四級熱電制冷器生產(chǎn)的處于非常溫運(yùn)行的面陣規(guī)格為320×256的MCT淺制冷型焦平面紅外探測器,該類探測器可在200 K和250 K兩個(gè)定點(diǎn)溫度下工作[16].

綜合有關(guān)多級熱電制冷的研究可見,一是研究模型中大多未考慮制冷器兩端換熱器件的影響,而外部熱阻有限溫差傳熱引起的不可逆性對制冷性能有很大影響;二是性能指標(biāo)多為制冷量和制冷系數(shù),缺少協(xié)調(diào)制冷能力和經(jīng)濟(jì)性能的綜合指標(biāo).有限時(shí)間熱力學(xué)是現(xiàn)代熱力學(xué)理論的一個(gè)重要分支,實(shí)現(xiàn)了熱力學(xué)、傳熱學(xué)、流體力學(xué)和化學(xué)反應(yīng)動(dòng)力學(xué)等多個(gè)基礎(chǔ)學(xué)科理論的交叉融合,著重研究在有限時(shí)間、有限尺寸、有限面積、有限速率傳熱等約束下的各類過程、循環(huán)和裝置的最優(yōu)性能,是分析能量轉(zhuǎn)換裝置的重要理論與有力工具[17-20],已廣泛應(yīng)用于化學(xué)熱機(jī)[21]、燃料電池[22]、布朗熱機(jī)[23]、熱電裝置[24]等能量轉(zhuǎn)換裝置的分析中,得到了一系列不同于傳統(tǒng)分析結(jié)果的新結(jié)論和一系列更接近實(shí)際、更具普適意義的新成果.本文基于有限時(shí)間熱力學(xué)理論,充分考慮外部熱源與熱電制冷器之間的不可逆性,建立空冷式多級熱電制冷器有限時(shí)間熱力學(xué)模型,提出協(xié)調(diào)性能系數(shù)評價(jià)指標(biāo),分析工作電流、熱電臂橫截面積和制冷溫差對裝置性能的影響,可為多級熱電制冷器的結(jié)構(gòu)設(shè)計(jì)與性能優(yōu)化提供參考.

1 多級熱電制冷器模型

圖1給出了多級熱電制冷器有限時(shí)間熱力學(xué)模型.圖中:I為電流.區(qū)別于非平衡熱力學(xué)模型,本模型考慮了熱電模塊端面與熱源之間的換熱器熱阻,體現(xiàn)出外部有限傳熱面積下的有限速率傳熱不可逆性對裝置性能的影響.裝置主要分為兩個(gè)部分:熱電制冷模塊和熱沉.熱電制冷模塊與環(huán)境空氣之間的熱量交換通過熱沉完成,并安裝風(fēng)扇加強(qiáng)換熱.裝置熱電偶總數(shù)為M,每個(gè)熱電偶由一個(gè)P型和一個(gè)N型熱電偶臂組成.由高溫端向低溫端各級熱電偶數(shù)分別為mi(1≤i≤n).TE為高溫?zé)嵩礈囟?TL為低溫?zé)嵩礈囟?熱電制冷器第1級熱端溫度(放熱溫度)和第n級冷端溫度(吸熱溫度)分別為Th和Tc,由第1級冷端至第n級熱端各級間溫度分別為Ti(1≤i≤n-1).高低溫?zé)嵩次艧崃繛镼E和QL,制冷模塊放、吸熱量分別為Qh和Qc,各節(jié)點(diǎn)熱流量分別為Qi(1 ≤i≤n-1).

圖1 多級熱電制冷器有限時(shí)間熱力學(xué)模型Fig.1 Finite time thermodynamic model of multistage thermoelectric cooler

忽略級間熱漏損失,兩級間熱流量等于高一級熱端的散熱量,同時(shí)等于低一級冷端的吸熱量.由非平衡熱力學(xué)理論,可得由高溫端至低溫端的各節(jié)點(diǎn)的熱流量,分別為

Qh=m1[αITh-K(Th-T1)+ 0.5I2R-0.5μI(Th-T1)]

(1)

Q1=m1[αIT1-K(Th-T1)- 0.5I2R+0.5μI(Th-T1)]=m2[αIT1-K(T1-T2)+ 0.5I2R-0.5μI(T1-T2)]

(2)

Qi=mi[αITi-K(Ti-1-Ti)- 0.5I2R+0.5μI(Ti-1-Ti)]=mi+1[αITi-K(Ti-Ti+1)+ 0.5I2R-0.5μI(Ti-Ti+1)] (2≤i≤n-2)

(3)

Qn-1=mn-1[αITn-1-K(Tn-2-Tn-1)- 0.5I2R+0.5μI(Tn-2-Tn-1)]=mn[αITn-1-K(Tn-1-Tc)+ 0.5I2R-0.5μI(Tn-1-Tc)]

(4)

Qc=mn[αITc-K(Tn-1-Tc)- 0.5I2R+0.5μI(Tn-1-Tc)]

(5)

式(2)～(4)可整理成關(guān)于未知節(jié)點(diǎn)溫度Ti(1≤i≤n-1)的由n-1個(gè)方程組成的方程組:

[(m1-m2)(αI-0.5μI)+(m1+m2)K]×T1-m2(K+0.5μI)T2=m1(K-0.5μI)Th+(m1+m2)0.5I2R

(6)

mi(0.5μI-K)Ti-1+[(mi-mi+1)(αI- 0.5μI)+(mi+mi+1)K]Ti-mi+1(K+ 0.5μI)Ti+1=(mi+mi+1)0.5I2R

(7)

mn-1(0.5μI-K)Tn-2+[(mn-1-mn)(αI-

0.5μI)+(mn-1+mn)K]Tn-1=

mn(K+0.5μI)Tc+(mn-1+mn)0.5I2R

(8)

由式(6)～(8)聯(lián)立解得中間連接層溫度Ti(1≤i≤n-1)為

(9)

(10)

(11)

考慮熱電模塊兩端與熱源不可逆?zhèn)鳠岬挠绊?根據(jù)有限時(shí)間熱力學(xué)理論,熱端和冷端的熱流量可分別表示為

(12)

式中:RE和RL分別為熱、冷端的傳熱熱阻.

由能量守恒方程式可得:

Qh=QE,Qc=QL

(13)

聯(lián)立式(1)～(5)、(9)～(13),得到熱電制冷器熱、冷端溫度.把Ti(1 ≤i≤n-1)、Th、Tc代入式(1)～(5)得流過各節(jié)點(diǎn)的熱流量,可得各級功率為

P1=Qh-Q1

(14)

Pi=Qi-1-Qi(2≤i≤n-1)

(15)

Pn=Qn-1-Qc

(16)

各級制冷系數(shù)為

(17)

εi=Qi/Pi(2≤i≤n-1)

(18)

εn=Qc/Pn

(19)

裝置整體的輸入功率和制冷系數(shù)分別為

(20)

(21)

2 熱阻分析

冷熱兩端的熱阻包括4個(gè)部分:陶瓷基板導(dǎo)熱熱阻Rcp、接觸熱阻Rc、熱沉基板熱阻Rex和熱沉對流換熱熱阻Rcv.式(12)中的RE和RL可表示為上述4個(gè)熱阻之和.圖2所示為相應(yīng)的熱阻網(wǎng)絡(luò)圖.圖中:RP和RN分別為P型和N型熱電偶的熱阻.

圖2 熱阻網(wǎng)絡(luò)Fig.2 Thermal resistance network

模塊的陶瓷基板熱阻、接觸熱阻和熱沉基板熱阻分別為

(22)

式中:δcp、δc和δex分別為陶瓷基板、導(dǎo)熱硅脂和熱沉基板的厚度;Acp、Ac和Aex分別為陶瓷基板、導(dǎo)熱硅脂和熱沉基板的面積;λcp、λc和λex分別為陶瓷、導(dǎo)熱硅脂和熱沉基板的導(dǎo)熱率.

熱電制冷器兩側(cè)通過氣體與固體壁面對流換熱,造成熱阻較大,需要采用熱沉和風(fēng)扇組合強(qiáng)化換熱.本文模型中采用等截面直肋熱沉,其結(jié)構(gòu)如圖3所示.圖中:δr為肋片厚度;br為肋片間距;H為肋片高度;W為肋片長度.

圖3 熱沉結(jié)構(gòu)Fig.3 Structure of heat sink

裝置冷熱兩側(cè)空氣為肋壁的強(qiáng)迫風(fēng)冷,對流換熱熱阻Rcv為[25]

(23)

式中:h1為對流換熱系數(shù);λr為肋片熱導(dǎo)率.

圖4給出了一定工況下,制冷量與制冷系數(shù)隨熱端散熱熱阻的變化趨勢.圖中:tE和tL分別為高溫?zé)嵩礈囟群偷蜏責(zé)嵩礈囟?當(dāng)其他參數(shù)設(shè)定,隨著冷熱兩端散熱熱阻的增大,裝置的制冷量和制冷系數(shù)逐漸降低,熱阻大小對制冷性能有著較大影響.因此,本文在熱力學(xué)模型中考慮了冷熱源與熱電模塊之間的有限速率傳熱,彌補(bǔ)了現(xiàn)有研究中未考慮制冷器兩端換熱器件外部熱阻有限溫差傳熱引起的不可逆性對制冷性能影響的不足.

圖4 制冷量、制冷系數(shù)與熱端熱阻的關(guān)系Fig.4 Relationship between cooling capacity, coefficient of performance and thermal resistance at hot end

3 評價(jià)指標(biāo)

3.1 熱力學(xué)完善度

制冷系數(shù)由設(shè)備運(yùn)行的環(huán)境溫度和制冷溫度決定,其數(shù)值無法反映出裝置的不可逆程度,在不同環(huán)境和制冷溫度條件下也沒有可比性.為了更好地對比反映裝置的性能,馬一太等[26]首次提出熱力學(xué)完善度ηre這一評價(jià)指標(biāo),并將其應(yīng)用于壓縮式制冷的計(jì)算中.熱力學(xué)完善度定義為裝置實(shí)際的制冷系數(shù)ε與同一溫限下逆向卡諾循環(huán)制冷系數(shù)εc之比,即

(24)

熱力學(xué)完善度反映了設(shè)備實(shí)際狀況與理想循環(huán)之間的差距以及不可逆性,可以將不同制冷溫差下的同一制冷裝置或者具有相同運(yùn)行模式的制冷裝置在相同的基礎(chǔ)上對比研究制冷特性,有助于更全面地反映裝置的制冷特性和不可逆程度.

3.2 協(xié)調(diào)性能系數(shù)

由熱電制冷的性能曲線可知,最大制冷量和最大制冷系數(shù)對應(yīng)不同的最佳電流,制冷能力與制冷的經(jīng)濟(jì)性無法同時(shí)達(dá)到最優(yōu).為了尋求制冷性能的協(xié)調(diào),文獻(xiàn)[27]中提出了將制冷量和制冷系數(shù)兩個(gè)性能指標(biāo)與各自最大值的比值再線性加權(quán),作為目標(biāo)函數(shù)的雙目標(biāo)優(yōu)化數(shù)學(xué)模型,即

(25)

式中:QLmax為所能達(dá)到的最大制冷量;εmax為所能達(dá)到的最大制冷系數(shù);x為0～1的比例系數(shù).

系數(shù)G體現(xiàn)了制冷量與制冷系數(shù)的協(xié)調(diào),x的取值大小反映了對這兩個(gè)性能指標(biāo)的偏好.當(dāng)性能系數(shù)G達(dá)到最大時(shí),表示一定偏好下制冷量與制冷系數(shù)協(xié)調(diào)達(dá)到最佳.

為了體現(xiàn)制冷量與耗功協(xié)調(diào)性能,定義協(xié)調(diào)性能系數(shù)F:

(26)

式中:Pmax為最大制冷量對應(yīng)的耗功.

系數(shù)F體現(xiàn)了制冷量與耗功的協(xié)調(diào),當(dāng)其取最大值時(shí),意味著取得了兩者的最佳折衷.耗功在一定程度上反映了經(jīng)濟(jì)性能的高低,系數(shù)F越大,說明制冷量越接近最大制冷量、同時(shí)耗功相對較小,表明了制冷性能和經(jīng)濟(jì)性能越理想.工程應(yīng)用中,當(dāng)追求裝置制冷量與經(jīng)濟(jì)性能的統(tǒng)一時(shí),推薦協(xié)調(diào)性能系數(shù)F作為衡量標(biāo)準(zhǔn);當(dāng)裝置對制冷量或制冷系數(shù)有所偏好,推薦采用協(xié)調(diào)性能系數(shù)G,可對x進(jìn)行不同的取值.

4 影響參數(shù)分析

熱電制冷器選用熱電材料碲化鉍(Bi2Te3),熱電臂橫截面積為A=1.4 mm×1.4 mm,單元長度為L=2 mm.選取熱電五級結(jié)構(gòu)制冷器進(jìn)行分析,第1級對數(shù)為m1=300,隨著級數(shù)的增加,熱帶單元對數(shù)逐級減半,總的熱電偶對數(shù)為582對.材料的物性參數(shù)會(huì)隨溫度的變化而變化,參數(shù)擬合公式[28]如下:

αP=-αN=(22 224.0+930.6T- 0.990 5T2)×10-9

(27)

ρP=ρN=(5 112.0+163.4T+ 0.627 9T2)×10-10

(28)

λP=λN=(62 605.0-277.7T+ 0.413 1T2)×10-4

創(chuàng)始人王勇的命運(yùn)轉(zhuǎn)折是改革開放40年中的標(biāo)志性事件，他從未避諱自己的出身，在多次采訪中，均表示自己本是鄒平電業(yè)局的一個(gè)臨時(shí)工?；蛟S是進(jìn)入體制無望，或許是時(shí)代創(chuàng)業(yè)浪潮的吸引，也或許是一個(gè)偶然的機(jī)遇，1986年，36歲的王勇離開電業(yè)局當(dāng)選西王村支部書記，開始在西王村創(chuàng)辦企業(yè)。

(29)

μP=-μN(yùn)=(930.6T-1.981T2)×10-9

(30)

式中:αj、ρj、λj和μj(j=P, N)分別表示材料的塞貝克系數(shù)、電阻率、熱導(dǎo)率和湯姆遜系數(shù).冷熱兩端換熱器一般為鋁制,文中討論溫度范圍內(nèi)熱導(dǎo)率變化較小,取為定值λex=222 W/(m·K);空氣比熱容、熱導(dǎo)率等物性參數(shù)取相應(yīng)熱源溫度下的取值.

數(shù)值計(jì)算中,熱電材料物性中的定性溫度采用熱端與冷端的平均溫度.熱端為300 K、冷端為213 K時(shí),熱電材料物性參數(shù)計(jì)算結(jié)果分別為:塞貝克系數(shù)為 1.957 6 V/K;電阻率為 8.833 5 Ω/m;熱導(dǎo)率為 1.855 4 W/(m·K);湯姆遜系數(shù)為 1.083 6 V/K.

4.1 工作電流的影響

圖5和圖6分別給出了不同低溫?zé)嵩礈囟认卵b置的制冷量QL和制冷系數(shù)ε隨工作電流的變化趨勢.由圖可見,當(dāng)電流在0.8～2.8 A的范圍內(nèi),隨著電流的增大,制冷量和制冷系數(shù)都呈現(xiàn)先增大后減小的趨勢,存在最佳電流,但是兩者的變化情況有所不同.當(dāng)?shù)蜏責(zé)嵩礈囟葹?60 ℃ 時(shí),制冷量最大為2.245 W,對應(yīng)的最佳工作電流為2.55 A;電流為1.3 A時(shí)對應(yīng)最大的ε為0.059.最大的QL和ε對應(yīng)的最佳電流相差較大,說明制冷能力與經(jīng)濟(jì)性能無法同時(shí)達(dá)到最優(yōu).隨著制冷溫度的降低,QL和ε均減小.相較于低溫?zé)嵩礈囟葹?50 ℃ 來說,溫度降低10 ℃,對應(yīng)的最大QL和ε分別降低了19.18%和29.37%,溫度降低20 ℃,對應(yīng)的最大QL和ε分別降低了38.21%和52.21%.小電流時(shí),焦耳熱效應(yīng)的影響很小可忽略不計(jì),但由于熱量與電流的平方呈正比,隨著電流的增加,對制冷量和制冷效率的影響逐漸凸顯,導(dǎo)致QL和ε存在最大值而后開始降低.隨著低溫?zé)嵩礈囟鹊慕档?QL取最大值所對應(yīng)的最佳工作電流基本保持不變,ε取最大值所對應(yīng)的最佳工作電流隨之增大.由圖6可以發(fā)現(xiàn),相同的ε對應(yīng)兩個(gè)不同的工作電流,考慮到制冷量可知,當(dāng)ε相同時(shí)選擇大電流可以獲得更好的制冷能力.

圖5 不同制冷溫度下制冷量與工作電流關(guān)系Fig.5 Relationship between cooling capacity and working current at different cooling temperatures

圖7給出了不同低溫?zé)嵩礈囟认聼崃W(xué)完善度ηre與工作電流I的關(guān)系.熱力學(xué)完善度隨電流的變化趨勢與ε相同,其值反映了裝置的不可逆程度,可以排除其他的影響因素直接進(jìn)行比較.當(dāng)?shù)蜏責(zé)嵩礈囟葹?50 ℃ 時(shí),熱力學(xué)完善度隨電流的增加先增大后減小,在I=1.15 A時(shí)取最大值為 0.028 8,對應(yīng)的工作電流稱為最佳工作電流.當(dāng)tL=-60 ℃時(shí),熱力學(xué)完善度在I=1.3 A時(shí)取最大值為 0.024 1,相較tL=-50 ℃ 降低了16.44%;當(dāng)tL=-70 ℃ 時(shí),熱力學(xué)完善度在I=1.5 A時(shí)取最大值為0.019 0,相較tL=-50 ℃ 降低了33.86%.隨著低溫?zé)嵩礈囟鹊慕档?ηre取最大值所對應(yīng)的最佳工作電流也隨之增大.與圖6的結(jié)果對比表明,不同制冷溫差下,εmax的變化相較ηremax的變化更大,因?yàn)棣莚e是個(gè)相對值,可以更直觀地反映出不同制冷溫差下裝置的性能,能夠直接進(jìn)行比較的優(yōu)點(diǎn)使得引入熱力學(xué)完善度這一評價(jià)指標(biāo)有著重要意義.圖8給出了最大制冷溫差與電流的關(guān)系.由圖可知,最大制冷溫差隨著工作電流的增大呈現(xiàn)先增大后減小的趨勢,存在最佳電流使得制冷溫差達(dá)到最大.當(dāng)電流為2.35 A時(shí),制冷溫差取最大值為129 ℃.結(jié)果表明,采用多級熱電制冷可以獲得較低的制冷溫度,實(shí)現(xiàn)大溫差制冷,在需要低溫工作環(huán)境的設(shè)備中具有廣泛應(yīng)用前景.

圖6 不同制冷溫度下制冷系數(shù)與工作電流關(guān)系Fig.6 Relationship between coefficient of performance and working current at different cooling temperatures

圖7 不同制冷溫度下熱力學(xué)完善度與工作電流關(guān)系Fig.7 Relationship between thermodynamic perfection and working current at different cooling temperatures

圖8 最大制冷溫差與工作電流的關(guān)系Fig.8 Relationship between maximum cooling temperature difference and working current

圖9所示為不同低溫?zé)嵩礈囟认聟f(xié)調(diào)性能系數(shù)與工作電流的關(guān)系.圖9(a)體現(xiàn)了制冷量與耗功的協(xié)調(diào),當(dāng)系數(shù)達(dá)到最大時(shí)表示一定工作條件下獲得最佳折衷.系數(shù)F隨著工作電流的增加呈現(xiàn)先增大后減小的趨勢,結(jié)果表明存在最佳工作電流使得F取最大值.當(dāng)tL=-50 ℃ 時(shí),在I=1.6 A時(shí)F取最大值為0.371,相較于tL=-60 ℃和tL=-70 ℃ 時(shí)的最大值分別增加了12.21%和31.64%.隨著低溫?zé)嵩礈囟鹊慕档?系數(shù)F取最大值所對應(yīng)的最佳工作電流也隨之增大.圖9(b)表示通過權(quán)衡尋找制冷量和制冷系數(shù)的全局最優(yōu),本文取x=0.5進(jìn)行計(jì)算.系數(shù)G隨著工作電流的增加呈現(xiàn)先迅速增大后緩慢減小的趨勢,存在最佳工作電流使得G取最大值.當(dāng)tL=-70 ℃時(shí),在I=1.9 A時(shí)G取最大值為0.885,相較于tL=-60 ℃ 和tL=-50 ℃ 時(shí)的最大值分別增加了4.37%和9.38%,結(jié)果表明低溫?zé)嵩礈囟葘ο禂?shù)G的影響較小,且當(dāng)電流到達(dá)一定大小后對系數(shù)的影響效果較小.通過與圖5和圖6對比可以看出,G取得最大值的工作電流位于QLmax和εmax對應(yīng)工作電流的中間,此時(shí)制冷量與制冷系數(shù)獲得最佳折衷,可以實(shí)現(xiàn)綜合效益最大化,充分體現(xiàn)了協(xié)調(diào)性能系數(shù)對雙目標(biāo)的協(xié)調(diào)意義,為熱電裝置運(yùn)行工況的選擇提供理論參考.

圖9 不同制冷溫度下協(xié)調(diào)性能系數(shù)與工作電流關(guān)系Fig.9 Relationship between coordination performance coefficient and working current at different cooling temperatures

4.2 熱電臂橫截面積的影響

圖10 不同工作電流下制冷量與橫截面積關(guān)系Fig.10 Relationship between cooling capacity and cross-sectional area at different working currents

圖11 不同工作電流下制冷系數(shù)與橫截面積關(guān)系Fig.11 Relationship between coefficient of performance and cross-sectional area at different working currents

圖12所示為不同工作電流下熱力學(xué)完善度ηre與橫截面積A的關(guān)系.熱力學(xué)完善度隨橫截面積的變化趨勢與制冷系數(shù)相同,但是在數(shù)值上相差較大,結(jié)果表明熱電裝置的不可逆程度較大.隨著工作電流的增大,最大ηre對應(yīng)的最佳橫截面積增大.工作電流I=1.5 A時(shí),橫截面積A=2.2 mm2時(shí)熱力學(xué)完善度ηre取最大值為0.024;當(dāng)I取2 A和2.5 A時(shí),最大熱力學(xué)完善度ηre對應(yīng)的最佳橫截面積分別為3.0 mm2和3.7 mm2.

圖12 不同工作電流下熱力學(xué)完善度與橫截面積關(guān)系Fig.12 Relationship between thermodynamic perfection and cross-sectional area at different working currents

圖13所示為不同工作電流下協(xié)調(diào)性能系數(shù)與橫截面積的關(guān)系.圖13(a)中系數(shù)F隨著橫截面積的增加呈現(xiàn)先迅速增大后較為緩慢減小的趨勢.工作電流I=1.5 A時(shí),橫截面積A=1.8 mm2時(shí)系數(shù)F取最大值為0.330.結(jié)果表明,通過協(xié)調(diào)存在最佳橫截面積使得制冷量與裝置耗功獲得最佳折衷.圖13(b)中系數(shù)G隨著橫截面積的增加呈現(xiàn)先增大后減小的趨勢,存在最佳橫截面積使G取最大值.工作電流I=1.5 A,橫截面積A=1.7 mm2時(shí)系數(shù)G取最大值為0.848;當(dāng)I取2 A和2.5 A時(shí),F和G的最大值對應(yīng)的最佳橫截面積分別為2.2 mm2和2.8 mm2.結(jié)果表明,通過賦予權(quán)重可以使制冷量與制冷系數(shù)之間有一個(gè)很好的權(quán)衡,有效協(xié)調(diào)兩者之間因?yàn)闄M截面積變化而出現(xiàn)此消彼長的矛盾.

圖13 不同工作電流下協(xié)調(diào)性能系數(shù)與橫截面積關(guān)系Fig.13 Relationship between coordination performance coefficient and cross-sectional area at different working currents

4.3 制冷溫差的影響

圖14和圖15分別給出了不同工作電流下裝置的制冷量QL和制冷系數(shù)ε隨制冷溫差Δt的變化趨勢.由圖可以發(fā)現(xiàn),當(dāng)制冷溫差在60～120 ℃的范圍內(nèi),當(dāng)I一定時(shí),制冷量和制冷系數(shù)隨著Δt的增大呈現(xiàn)單調(diào)減小的趨勢.由于溫差的不斷增大,平均溫度降低,熱電模塊的物性參數(shù)發(fā)生改變,造成制冷量的不斷下降.圖16所示為不同工作電流下熱力學(xué)完善度ηre與制冷溫差Δt的關(guān)系.工作電流I=2.5 A時(shí),制冷溫差Δt=78 ℃時(shí)熱力學(xué)完善度ηre取最大值為0.013;當(dāng)I取2 A和1.5 A時(shí),制冷溫差分別為78 ℃和70 ℃時(shí)對應(yīng)ηre的最大值,分別為0.019和0.026,相比I=2.5 A條件下的最大ηre分別增大了42.2%和96.6%.結(jié)果表明,熱力學(xué)完善度隨制冷溫差的增大呈現(xiàn)先緩慢增大后減小的趨勢,與制冷量和制冷系數(shù)的變化趨勢有所不同,存在最佳制冷溫差使得熱力學(xué)完善度最大,即不可逆程度達(dá)到最小.這是由于對于給定的環(huán)境溫度,逆卡諾循環(huán)制冷系數(shù)εc=TE/Δt-1與溫差呈反比關(guān)系,熱電制冷器制冷系數(shù)與溫差呈類拋物線關(guān)系,導(dǎo)致存在最佳溫差使得熱電制冷器制冷系數(shù)最接近逆卡諾循環(huán)制冷系數(shù),即熱力學(xué)完善度最大.這個(gè)指標(biāo)獨(dú)立于制冷量和制冷系數(shù),為系統(tǒng)的評價(jià)與優(yōu)化提供了更全面的理論基礎(chǔ).

圖14 不同工作電流下制冷量與制冷溫差關(guān)系Fig.14 Relationship between cooling capacity and temperature difference at different working currents

圖15 不同工作電流下制冷系數(shù)與制冷溫差關(guān)系Fig.15 Relationship between coefficient of performance and temperature difference at different working currents

圖16 不同工作電流下熱力學(xué)完善度與制冷溫差關(guān)系Fig.16 Relationship between thermodynamic perfection and temperature difference at different working currents

5 結(jié)論

基于有限時(shí)間熱力學(xué)理論,建立了空冷式多級熱電制冷器熱力學(xué)模型,提出了新的性能評價(jià)指標(biāo),分析了工作電流、熱電臂橫截面積和制冷溫差對裝置性能的影響,主要結(jié)論有:

(1) 當(dāng)制冷溫差為87 ℃時(shí),工作電流分別為2.55 A和1.30 A,制冷量和制冷系數(shù)可分別達(dá)到最大值2.245 W和0.059;協(xié)調(diào)性能系數(shù)能夠?qū)⒅评淞?、制冷系?shù)和耗功進(jìn)行協(xié)調(diào)后得到綜合性能的最佳折衷,獲得最佳運(yùn)行參數(shù).電流分別為1.70 A 和1.75 A時(shí),協(xié)調(diào)性能系數(shù)F和G分別取得最大值,獲得制冷量、耗功與制冷系數(shù)的最佳折衷.

(2) 一定工況下存在最佳熱電臂橫截面積對應(yīng)最優(yōu)的評價(jià)指標(biāo),超過最佳值后,由于模塊電阻和熱阻減小、模塊產(chǎn)熱量和熱漏量增大,制冷性能降低.綜合考慮制冷能力和經(jīng)濟(jì)性能的協(xié)調(diào),當(dāng)工作電流為2 A時(shí),最佳橫截面積為2～3 mm2.

(3) 制冷量和制冷系數(shù)隨制冷溫差增大逐漸減小;熱力學(xué)完善度隨制冷溫差增大呈現(xiàn)先增大后減小的趨勢,即存在最佳制冷溫差使得熱力學(xué)完善度最大,系統(tǒng)不可逆程度最小.