熱電發(fā)電裝置的動態(tài)模型及分析

袁義生，熊 晶

(華東交通大學(xué)電氣學(xué)院，江西 南昌330013)

1 引言

能源緊張及環(huán)境污染的加劇使得廢熱回收技術(shù)成為提高能源利用率的一個重要研究方向。以汽車為例，傳統(tǒng)汽車石油燃料產(chǎn)生的能量約30%轉(zhuǎn)換成機械能來驅(qū)動汽車行駛，另外約40%以廢氣方式排放，還有約30%在引擎冷卻過程中流失，其能源利用率極低［1-2］。因此，研究如何將這些廢熱回收利用是一項極其有益的研究工作。

熱電發(fā)電裝置由熱端導(dǎo)熱片、熱電發(fā)電器件和冷端散熱片構(gòu)成，之后再經(jīng)后級直流變換器將熱電發(fā)電器件發(fā)出的電能轉(zhuǎn)換后送給蓄電池儲存。

但熱電發(fā)電及變流裝置要得到廣泛應(yīng)用，還有一些基礎(chǔ)性問題需要解決，其中重要一點就是熱電發(fā)電裝置的數(shù)學(xué)模型的建立。這是因為后級直流變流器的控制特性，包括做最大功率跟蹤和功率匹配［3］時的動態(tài)特性都與前級熱電發(fā)電裝置的動靜態(tài)模型密切相關(guān)。然而，現(xiàn)有的熱電發(fā)電器件的模型有兩種，一種是設(shè)計模型，另外一種是電等效模型。設(shè)計模型是基于器件物理特性而建立的微觀參數(shù)模型，主要為評估器件Seebeck效應(yīng)建立，此類模型主要是不考慮器件工作條件的物理靜態(tài)模型［4］，也有少數(shù)考慮了外部工作條件的準動態(tài)模型［5］。電等效模型則是將熱電發(fā)電器件的輸出等效為電壓源U和內(nèi)阻Rin串聯(lián)電路。此類模型可以由廠家給出特性參數(shù)。但問題是，器件廠家標示的器件參數(shù)是在器件兩端溫差ΔT和負載固定的靜態(tài)條件下測量得到的，并不反映熱電發(fā)電器件的動態(tài)特性，尤其是大信號特性。而且，熱電發(fā)電器件處于熱端尾氣管和冷端散熱器之間，這兩者的特性對熱電器件的輸出電壓動態(tài)特性都有影響;還有就是后級直流變換器的等效輸入電阻對熱電發(fā)電器件的輸出電壓也有影響。這些都是之前的研究沒有系統(tǒng)考慮的問題。

所以，本文建立熱端導(dǎo)熱片、熱電發(fā)電器件和冷端散熱器構(gòu)成的完整熱電發(fā)電裝置的動態(tài)模型。分析熱電發(fā)電器件輸出電壓的熱源響應(yīng)特性，以及輸出負載響應(yīng)特性。建立等效仿真電路模型，并用小裝置試驗測試驗證了模型的正確性。

2 熱電發(fā)電裝置結(jié)構(gòu)及原理

熱電發(fā)電裝置由熱端導(dǎo)熱片(散熱片)、熱電發(fā)電器件和冷端散熱器構(gòu)成，如圖1所示。其中，因為熱電發(fā)電器件為雙層半導(dǎo)體串聯(lián)結(jié)構(gòu)，表面不絕緣，所以在熱冷端散熱器與熱電發(fā)電器件之間還有絕緣陶瓷片。而虛線框的后級直流變換器部分可用一個等效負載電阻Ro來代替。

圖1 熱電發(fā)電裝置結(jié)構(gòu)Fig.1 Structure of thermoelectric generator equipment

熱電發(fā)電器件Seebeck效應(yīng)的輸出電壓表示為:

式中，α是 Seebeck系數(shù);ΔT是熱電發(fā)電器件熱端和冷端間的溫差。α值由材料特性決定。但ΔT卻不僅取決于熱電發(fā)電器件的傳熱特性，由圖1可知也與熱源變化、冷熱端散熱片傳熱特性和等效負載電阻Ro有關(guān)。所以，建立完整的熱電發(fā)電裝置模型需要逐個建立上述部件數(shù)學(xué)模型并連接構(gòu)成。

3 散熱片動態(tài)模型

對任意微小封閉容積，據(jù)熱力學(xué)第一定律［8］有:

式中，ρ、v、c分別為此封閉容積的密度、體積、比熱容;q表示能量流動的速率，即熱流量。此式表示控制容積中的能量速率的增大，必定等于進入控制容積的能量速率，減去離開控制容積的能量速率，再加上在同一瞬間控制容積內(nèi)產(chǎn)生熱量的速率。

由于導(dǎo)熱系數(shù)恒定，溫度沿厚度方向線性分布。如果散熱片很薄，可以使用兩端平均溫度代替散熱片溫度。定義其兩端溫度分別為T1和T2，則有:

式中，m為封閉容積的質(zhì)量。再結(jié)合基本的靜態(tài)熱阻定律:

式中，θ表示散熱片的熱阻?？梢缘玫綗崂涠松崞魍暾撵o動態(tài)特性等效電路如圖2所示。

4 熱電發(fā)電器件的動態(tài)等效模型

圖2 平板散熱器等效電路圖Fig.2 Equivalent model of flat heatsink

為了提高熱電轉(zhuǎn)換效率，熱電發(fā)電器件通常制作成平板狀。當熱電發(fā)電器件輸出接負載電阻Ro時，體內(nèi)的熱能就轉(zhuǎn)換成電能，有:

式中，Rm是熱電發(fā)電器件本身電阻。

對于半導(dǎo)體熱電模塊，考慮塞貝克效應(yīng)及帕爾帖效應(yīng)，單位時間內(nèi)從熱接觸層吸收的熱量為:

熱電器件單位時間內(nèi)放給冷端接觸層的熱量為［9］:

式中，Ta和Tc分別為熱電模塊熱面和冷面的溫度;θm為熱電模塊的熱阻。并且:

根據(jù)熱力學(xué)第一定律可得:

可以建立熱電模塊的等效模型，如圖3所示。

圖3 熱電發(fā)電器件等效模型Fig.3 Equivalent model of thermoelectric generator

5 完整模型及分析

5.1 完整模型

在上述分析的基礎(chǔ)上，據(jù)圖1裝置結(jié)構(gòu)可以得到完整的熱電發(fā)電裝置的等效電路模型如圖4所示。

圖4 熱電發(fā)電裝置等效電路模型Fig.4 Equivalent circuit model of thermoelectric equipment

在圖4中，T1、T2分別為熱端和冷端溫度;θ1為熱端散熱片的熱阻;C1為熱端散熱片的熱容，其值為熱端散熱片的質(zhì)量與比熱容的乘積的1/2;C3為熱電發(fā)電器件的熱容;θ4為冷端散熱片的總熱阻;C2為冷端散熱片的熱容;θ2和 θ3分別為熱端和冷端導(dǎo)熱陶瓷墊片熱阻和接觸熱阻之和;Rm和Ro分別為熱電發(fā)電器件的內(nèi)阻和等效負載電阻。因為絕緣陶瓷片熱容小，因此圖中忽略。

5.2 模型分析

(1)源響應(yīng)

源響應(yīng)是指熱源波動對Seebeck輸出電壓的影響。假定熱源信號階躍變化，而熱電發(fā)電器件輸出開路，即在圖4中設(shè)定I=0，Ro斷開，就可以分析模型的源響應(yīng)特性。這是一個多階的熱阻容θC慣性電路。熱電發(fā)電器件Seebeck效應(yīng)輸出電壓U與ΔT(Ta－Tc)成正比。當 qin是一個初始為零的階躍函數(shù)時，ΔT的穩(wěn)態(tài)值就等于qinθm，但其暫態(tài)分量則取決于各串聯(lián)的熱阻容乘積。在實際的系統(tǒng)中，熱端散熱片熱容大;熱電發(fā)電器件的熱容則最小;冷端散熱片的熱容較大。所以，ΔT的暫態(tài)分量主要取決于θ1C1和θ4C2兩個時間常數(shù)。即 Seebeck效應(yīng)輸出電壓U的源階躍響應(yīng)特性曲線是一個簡化的二階慣性環(huán)節(jié)。

(2)負載響應(yīng)

負載響應(yīng)是指等效負載電阻Ro(后級直流變換器電路等效輸入電阻)階躍變化時，Seebeck效應(yīng)輸出電壓U的響應(yīng)特性。此時可以假定熱源穩(wěn)定。等效電路如圖5所示。

分析負載由開路(R=∞)階躍變化至 Ro時的輸出電壓響應(yīng)。由于Ro和Rm的分壓，此時負載電壓Uo會瞬間降低。之后，因為IU/2的分流作用，流過 θm的熱量 qm下降;且(αTa－IRm/2)Iθm的存在，使得熱冷端溫差ΔT下降。所以，負載電壓Uo會繼續(xù)隨Seebeck電壓U下降而下降。由于熱電發(fā)電器件的效率很低，IU/2遠小于qa，所以負載電壓Uo繼續(xù)下降的過程中可以忽略熱端散熱器的影響qa不變。此時，負載電壓Uo基本按θmC3時間常數(shù)規(guī)律下降至穩(wěn)定值。

圖5 負載響應(yīng)等效電路Fig.5 Equivalent circuit of load response

綜上所述，可以看出:熱電發(fā)電裝置輸出電壓的變化在其物理特性基礎(chǔ)上，不僅與熱源的改變，也與后級負載的變化有關(guān)。

6 仿真與實驗

在實驗室搭建了一個小型熱電發(fā)電裝置，其中的熱電發(fā)電器件采用了一片Hi-Z公司生產(chǎn)的HZ-14(可見該公司網(wǎng)站)，熱源采用熱風槍來模擬，冷端采用強迫風冷使冷端散熱器外側(cè)溫度與空氣溫度接近。

按照圖4在電路仿真軟件SABER中建立仿真電路模型。其中的溫度，熱阻和熱容分別用電壓，電阻和電容來代替。

系統(tǒng)各部分的熱阻和熱容按照文獻［8］中公式計算，熱電發(fā)電器件的熱阻在說明書中查。參數(shù)見表1。

表1 仿真和試驗參數(shù)Tab.1 Parameters of simulation and experiment

實驗首先測試了裝置的負載響應(yīng)特性。在熱源及溫度穩(wěn)定的條件下，投卸等效輸出電阻 Ro，測試其輸出電壓響應(yīng)特性。其中，圖6(a)和6(b)分別為投載下仿真和實驗的負載電壓Uo波形圖。圖7(a)和7(b)分別為卸載下仿真和實驗的負載電壓Uo波形圖。由圖可見，投載瞬間負載電壓Uo先迅速下降，而后沿曲線下降約5s后穩(wěn)定。卸載時相反。實驗和仿真結(jié)果吻合。

圖6 投載下Uo仿真和實驗波形Fig.6 Simulation and experimental waveforms ofvoltage Uoin load-on transient

圖7 卸載下Uo仿真和實驗波形Fig.7 Simulation and experimental waveforms of voltage Uoin load-off transient

其次測試了裝置的源響應(yīng)特性。在熱源加熱溫度穩(wěn)定后斷開熱源，分別考察空載和帶載情況下的負載電壓Uo響應(yīng)。圖8(a)和8(b)分別是空載下仿真和實驗波形。圖9(a)和9(b)分別是帶載下仿真和實驗波形。可見，源響應(yīng)負載電壓Uo下降時間達到1800s，遠遠超過負載響應(yīng)時間，這與第4節(jié)的分析吻合。另外，帶載下的源響應(yīng)速度略快于空載下的源響應(yīng)速度，但并不顯著，這是因為熱電發(fā)電器件的效率低，輸出電功率遠小于傳遞的熱量。

圖8 空載下斷熱源時Uo仿真和實驗波形Fig.8 Simulation and experimental waveform of voltage Uoin hot source cut-off transient without load

圖9 帶載下斷熱源時Uo仿真和實驗波形Fig.9 Simulation and experimental waveform of voltage Uoin hot source cut-off transient with load

7 結(jié)論

建立了熱電發(fā)電裝置的動態(tài)模型和等效電路模型，分析了其源響應(yīng)和負載響應(yīng)特性。指出了熱電發(fā)電器件的輸出電壓源響應(yīng)是一個由散熱器決定的大時間慣性特性曲線，而負載響應(yīng)則是由熱電發(fā)電器件決定的小時間慣性特性曲線。這對后級直流變換器主電路和控制電路的進一步研究具有重要參考價值。建立的小型實驗裝置的實驗和仿真結(jié)果吻合，證明了模型的有效性。

［1］Yang Jihui.Potential applications of thermoelectric waste heat recovery in the automotive industry［A］.2005 International Conference on Thermoelectronics［C］.New York，USA，Aug.2005.19-23.

［2］Yang Jihui.Development of thermoelectric technology for automotive waste heat recovery［R］.Washington:DOE USA，2008.

［3］Rae Young Kim，Jih-Sheng Lai.A seamless mode transfer maximum power point tracking controller for thermoelectric generator applications［J］.IEEE Trans.on Power Electronics，2008，23(5):2310-2317.

［4］宋瑞銀，李偉，楊燦軍，等 (Song Ruiyin，Li Wei，Yang Canjun，et al.).微小型熱電發(fā)電器建模及優(yōu)化設(shè)計研究(Modeling and optimal design for mic-thermoelectric generator)［J］.太陽能學(xué)報 (Acta Energiae Solaris Sinica)，2006，27(6):554-558.

［5］Rasit Ahiske.New method for investigation of dynamic parameters of thermoelectric modules［J］.Electric Engin，2007，15(1):51-65.

［6］Simon Lineykin，Sam Ben-Yaakov.Modeling and analysis of thermoelectric modules［A］.IEEE APEC 2005［C］.New York，USA，2005.2019-2023.

［7］Chavez J，Ortega J，Salazar J，et al.Spice model of thermoelectric elements including thermal effects［A］.Proc.Instrumentation and Measurement Tech.Conference［C］.New York，USA，2000.1019-1023.

［8］Incropera F P，葛新石譯 (translated by Ge Xinshi).傳熱和傳質(zhì)基本原理 (Fundamentals of heat and mass transfer)［M］.北京:化學(xué)工業(yè)出版社 (Beijing:Chemical and Industry Press)，2009.

［9］Omer S A，Infield D G.Design optimization of thermoelectric devices for solar power generation［J］.Solar Energy Materials and Solar Cells，1998，(53):67-82.