吳麟麟 胡迎香 汪若塵 羅丁

DOI： 10.3969/j.issn.1671-7775.2024.03.003

開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

摘要： 為了預(yù)測(cè)溫差發(fā)電（thermoelectric generator，TEG）系統(tǒng)的動(dòng)態(tài)特性，基于COMSOL Multiphy-sics建立了用于求解溫差發(fā)電系統(tǒng)溫度場(chǎng)分布的瞬態(tài)計(jì)算流體力學(xué)（computational fluid dynamics，CFD）模型和用于研究溫差發(fā)電模塊瞬態(tài)響應(yīng)特性的分析模型，提出了混合瞬態(tài)CFD-分析模型，并經(jīng)過(guò)瞬態(tài)試驗(yàn)驗(yàn)證.結(jié)果表明：由于熱慣性的影響，TEG系統(tǒng)的轉(zhuǎn)化效率會(huì)出現(xiàn)一個(gè)瞬時(shí)的較高值；相較于尾氣溫度和質(zhì)量流量的瞬態(tài)波動(dòng)，熱電半導(dǎo)體的熱端溫度和冷端溫度會(huì)存在時(shí)滯；在美國(guó)環(huán)保局的高速公路燃油經(jīng)濟(jì)性測(cè)試（highway fuel economy test，HWFET）模式循環(huán)工況下，瞬態(tài)模型求解得到整個(gè)溫差發(fā)電系統(tǒng)的平均輸出功率、平均轉(zhuǎn)化效率分別為35.63 W和3.40%，瞬態(tài)模型的輸出電壓平均誤差為6.41%；該模型能夠以較高的精度及較短的計(jì)算時(shí)間預(yù)測(cè)溫差發(fā)電系統(tǒng)在瞬態(tài)熱源激勵(lì)下的瞬態(tài)響應(yīng)特性.

關(guān)鍵詞：? 溫差發(fā)電系統(tǒng)； 尾氣余熱回收； 混合瞬態(tài)CFD-分析模型； 瞬態(tài)響應(yīng)特性； 熱慣性

中圖分類號(hào)： U466; TK11? 文獻(xiàn)標(biāo)志碼：? A? 文章編號(hào)：?? 1671-7775（2024）03-0265-08

引文格式：? 吳麟麟，胡迎香，汪若塵，等. 汽車尾氣溫差發(fā)電系統(tǒng)瞬態(tài)回收性能分析［J］.江蘇大學(xué)學(xué)報(bào)（自然科學(xué)版），2024，45（3）：265-272.

收稿日期：?? 2022-02-24

基金項(xiàng)目：? 國(guó)家自然科學(xué)基金資助項(xiàng)目（51977100）

作者簡(jiǎn)介：?? 吳麟麟（1970—），男，江蘇泰興人，教授 （llwu@ujs.edu.cn），主要從事汽車尾氣溫差發(fā)電研究.

胡迎香（1996—），女，河南信陽(yáng)人，碩士研究生（ujs_cj@outlook.com），主要從事汽車尾氣溫差發(fā)電研究.

Transient recovery performance investigation of thermoelectric

generation system for vehicle exhaust

WU Linlin， HU Yingxiang， WANG Ruochen， LUO Ding

（School of Automotive and Traffic Engineering， Jiangsu University， Zhenjiang， Jiangsu 212013， China）

Abstract： To predict the dynamic characteristics of the thermoelectric generator （TEG） system， a transient computational fluid dynamics （CFD） model for solving the temperature field distribution of TEG system and an analysis model for investigating the transient response characteristics of TEG modules were established based on COMSOL Multiphysics. A hybrid transient CFD-analysis model was proposed and validated through transient experiments. The results show that due to the thermal inertia effect， the conversion efficiency of TEG system exhibits momentary peak value. Compared to the transient fluctuations of exhaust gas temperature and mass flow rate， there is time lag in the hot and cold side temperatures of the thermoelectric semiconductor. Under the highway fuel economy test （HWFET） cycle conditions defined by the U.S. Environmental Protection Agency， by the transient model， the average output power and the average conversion efficiency of the entire TEG system are solved with 35.63 W and 3.40%， respectively. The average error of the transient model′s output voltage is 6.41%. By the proposed model， the transient response characteristics of the TEG system under transient thermal source excitation can be predicted with high accuracy and short computation time.

Key words：? thermoelectric generator system; exhaust waste heat recovery; hybrid transient CFD-analysis model; transient response characteristics; thermal inertia

發(fā)動(dòng)機(jī)在工作過(guò)程中，尾氣攜帶大量熱量以廢氣的形式排放到空氣中，降低了發(fā)動(dòng)機(jī)的工作效率.YANG J.等［1］研究表明，汽車發(fā)動(dòng)機(jī)排氣損耗占總能量的比例約為40%，造成了巨大的能源浪費(fèi)和環(huán)境污染.將發(fā)動(dòng)機(jī)尾氣余熱進(jìn)行回收利用是提高發(fā)動(dòng)機(jī)整體工作效率的有效手段，也是實(shí)現(xiàn)汽車綠色節(jié)能的有效途徑.溫差發(fā)電（thermoelectric generator，TEG）系統(tǒng)對(duì)于回收尾氣熱能具有無(wú)運(yùn)動(dòng)部件、無(wú)噪聲和無(wú)排放［2］等優(yōu)點(diǎn)，市場(chǎng)前景廣闊.

為了研究TEG系統(tǒng)在汽車尾氣廢熱回收領(lǐng)域中的應(yīng)用前景，國(guó)內(nèi)外學(xué)者制造了大量樣機(jī)進(jìn)行試驗(yàn).S. EZZITOUNI等［3］研究了TEG樣機(jī)集成到柴油機(jī)上的輸出性能，結(jié)果表明，在與外界絕熱的情況下，TEG系統(tǒng)的輸出功率可達(dá)43 W.WANG Y. P.等［4］將TEG樣機(jī)安裝到越野車上進(jìn)行道路試驗(yàn)，結(jié)果表明，在70 Ω的外部電阻和125 km/h的車速下，TEG系統(tǒng)可產(chǎn)生609 W的輸出功率.

由于試驗(yàn)設(shè)備、條件及成本的限制，進(jìn)行試驗(yàn)難度較大，測(cè)試精度較低，為了指導(dǎo)試驗(yàn)樣機(jī)進(jìn)行結(jié)構(gòu)優(yōu)化與設(shè)計(jì)，學(xué)者們建立了理論模型對(duì)TEG系統(tǒng)進(jìn)行分析預(yù)測(cè).在目前的理論模型中，計(jì)算流體力學(xué)（computational fluid dynamics，CFD）模型和分析模型應(yīng)用較為廣泛.CFD模型可以求解出TEG系統(tǒng)的物理場(chǎng)分布，且具有較高的精度，是使用廣泛的理論分析模型.N. KEMPF等［5］使用CFD模型對(duì)穩(wěn)態(tài)的TEG系統(tǒng)換熱器進(jìn)行數(shù)值分析，結(jié)果表明，在靠近換熱器入口處增加溫差發(fā)電模塊（thermoelectric module，TEM）個(gè)數(shù)并在靠近換熱器出口處使用熱導(dǎo)率更高的換熱器材料，可以最大限度提高燃油效率.A. MASSAGUER等［6］針對(duì)TEM建立了三維穩(wěn)態(tài)的CFD模型，將恒定溫度邊界分別加載在溫差發(fā)電模塊的熱端和冷端，求解得到溫差發(fā)電模塊的電壓分布云圖，并分析了溫差發(fā)電模塊在不同溫差下的輸出性能.以上CFD模型都基于穩(wěn)態(tài)工況［7］，瞬態(tài)的CFD模型大多限于對(duì)溫差發(fā)電模塊的研究［8-9］.

雖然CFD模型可用于分析TEG系統(tǒng)的熱力學(xué)性能，但用于整個(gè)TEG系統(tǒng)時(shí)會(huì)消耗大量的計(jì)算機(jī)資源及計(jì)算時(shí)間.而分析模型可以快速獲得理論分析結(jié)果，指導(dǎo)溫差發(fā)電裝置的初步性能分析［10］.S. KUMAR等［11］對(duì)整個(gè)TEG系統(tǒng)建立穩(wěn)態(tài)分析模型，研究結(jié)果表明，在尾氣溫度、質(zhì)量流量分別為550 ℃和35 g/s的條件下，TEG系統(tǒng)的輸出功率和轉(zhuǎn)化效率分別為552 W和3.33%.WANG Y. C.等［12］通過(guò)建立穩(wěn)態(tài)分析模型對(duì)TEM進(jìn)行優(yōu)化，結(jié)果表明，根據(jù)溫度分布選擇TEM相應(yīng)的材料可以顯著提高輸出功率和輸出效率.以上分析模型都是穩(wěn)態(tài)的，雖然有針對(duì)整個(gè)TEG系統(tǒng)的瞬態(tài)分析模型，但是TEG系統(tǒng)內(nèi)部的傳熱極其復(fù)雜，主要包括流體的流動(dòng)、流體與固體之間的傳熱和翅片的傳熱等，難以準(zhǔn)確刻畫它們之間的熱量傳遞，導(dǎo)致誤差偏大.然而，針對(duì)TEM的分析模型忽略了與外界的熱損失，基于能量守恒計(jì)算其輸出性能，誤差較低.

筆者針對(duì)上述研究不足，完成從穩(wěn)態(tài)到瞬態(tài)，從單一溫差發(fā)電模塊到整個(gè)TEG系統(tǒng)的跨越.首先，針對(duì)換熱器內(nèi)部的流場(chǎng)和熱場(chǎng)區(qū)域建立CFD模型；其次，針對(duì)TEM的熱場(chǎng)和電場(chǎng)區(qū)域建立分析模型，不僅可以得到TEG系統(tǒng)在瞬態(tài)工況下較為精確的溫度場(chǎng)分布結(jié)果，還可以較快計(jì)算TEG系統(tǒng)的輸出性能；最后，進(jìn)行瞬態(tài)試驗(yàn)驗(yàn)證，分析尾氣溫度和尾氣質(zhì)量流量對(duì)TEG系統(tǒng)瞬態(tài)余熱回收性能的影響.

1? 混合瞬態(tài)CFD-分析模型

1.1? TEG系統(tǒng)結(jié)構(gòu)

TEG系統(tǒng)的基本工作原理是利用TEM的塞貝克效應(yīng)將尾氣的熱能轉(zhuǎn)化為電能.本研究采用的是具有穩(wěn)定性能的平板式TEG系統(tǒng)，如圖1所示.

TEG系統(tǒng)由換熱器、TEM及冷卻裝置組成，其中，單個(gè)TEM包含128對(duì)P型和N型熱電半導(dǎo)體、256個(gè)銅電極片和2個(gè)陶瓷板.TEG系統(tǒng)的主要傳熱過(guò)程是尾氣攜帶熱量從換熱器入口流入，以對(duì)流換熱形式將熱量傳遞到換熱器，由換熱器通過(guò)熱傳導(dǎo)將吸收的熱量傳遞到TEM，最后，TEM通過(guò)熱電效應(yīng)將熱能轉(zhuǎn)化為電能.在換熱器內(nèi)部插入集熱翅片，其作用是增大換熱面積從而強(qiáng)化換熱器的傳熱效果，提高TEG系統(tǒng)的輸出效率.采用水冷循環(huán)進(jìn)行冷卻，原因是其較風(fēng)冷散熱具有更高的比熱容及更好的冷卻性能.

選擇重型卡車為研究對(duì)象，將TEG系統(tǒng)安裝在排氣管與消聲器之間.為了匹配排氣管的尺寸，換熱器的進(jìn)出口直徑都為50.0 mm，換熱器主體的長(zhǎng)、寬、高分別為212.0、110.0、70.0 mm，換熱器內(nèi)部包含20個(gè)長(zhǎng)、寬、高分別為200.0、2.0、20.0 mm的集熱翅片.基于Bi2Te3制成的TEM的長(zhǎng)、寬、高分別為40.0、40.0、0.8 mm.冷卻裝置的進(jìn)出口直徑都為5.5 mm，長(zhǎng)、寬、高分別為200.0、40.0、12.0 mm.TEG系統(tǒng)的材料屬性如表1所示，其中T為溫度.

1.2? TEG系統(tǒng)的瞬態(tài)CFD模型

1.2.1? 幾何簡(jiǎn)化

幾何簡(jiǎn)化包括TEM的結(jié)構(gòu)簡(jiǎn)化和TEG系統(tǒng)的結(jié)構(gòu)簡(jiǎn)化.TEG系統(tǒng)的簡(jiǎn)化流程如圖2所示.

為了減少仿真時(shí)間，將由陶瓷板、銅電極片及熱電半導(dǎo)體組成的TEM簡(jiǎn)化成一個(gè)導(dǎo)熱板.為了避免簡(jiǎn)化后的TEM產(chǎn)生誤差，簡(jiǎn)化后的TEM與簡(jiǎn)化前的TEM熱導(dǎo)率保持相等.由于整個(gè)TEG系統(tǒng)完全對(duì)稱，所以采用1/4簡(jiǎn)化模型進(jìn)行仿真，同時(shí)保持仿真精度并節(jié)省仿真時(shí)間.由于1/4 TEG系統(tǒng)4個(gè)TEM是串聯(lián)的，因此，電流相同，從尾氣入口到出口方向分別定義為TEM1、TEM2、TEM3、TEM4.熱電半導(dǎo)體結(jié)構(gòu)如圖3所示，其中：hce、hco、hte分別為陶瓷板、銅電極片和熱電半導(dǎo)體的高度.

簡(jiǎn)化后TEM與未簡(jiǎn)化前TEM的熱導(dǎo)率相等，熱導(dǎo)率計(jì)算公式為

λTEM=2hce+2hco+hteAce2hceλceAce+2hcoλcoAco+128hteλteAte Th（t）+Tc（t）Th（t）-Th_co（t），（1）

式中：Ace、Aco、Ate分別為陶瓷板、銅電極片和熱電半導(dǎo)體的橫截面積；λce、λco、λte分別為陶瓷板、銅電極片和熱電半導(dǎo)體的熱導(dǎo)率；Th（t）、Tc（t）分別為TEM的熱端溫度和冷端溫度；Th_co（t）為銅電極片的熱端溫度.

將表1數(shù)據(jù)代入式（1），可得簡(jiǎn)化后TEM熱導(dǎo)率與溫度之間的關(guān)系：

λTEM=-2.91×10-9T4+4.33×10-6T3-

0.002 4T2+0.58T-51.06.（2）

1.2.2? 控制方程

COMSOL Multiphysics和ANSYS是2個(gè)常用于計(jì)算CFD模型的有限元分析軟件，選擇COMSOL Multiphysics軟件建立瞬態(tài)CFD模型來(lái)研究TEG系統(tǒng)的熱力學(xué)性能.

換熱器和冷卻裝置中的流體流動(dòng)滿足如下質(zhì)量、動(dòng)量和能量守恒方程：

ρt+·（ρv）=0，（3）

t（ρv）+·（ρvv）=·μ（v+vT）-p，（4）

t（ρT）+ρvT=·λcT，（5）

式中：ρ為流體密度；t為時(shí)間；v為速度矢量；μ為動(dòng)力黏度；p為流體壓力；λ為熱導(dǎo)率；c為比熱容.

尾氣和水的流動(dòng)遵循k-ε湍流模型方程：

t（ρk）+xi（ρkui）=xjμ+μtσkkxj+

Gk+Gb-ρε-YM，（6）

t（ρε）+xi（ρεui）=xjμ+μtσεεxj-

Rε+C1εεk（Gk+C3εGb）-C2ερε2k ，（7）

式中： k為湍流動(dòng)能；ui為雷諾平均速度；μt為湍動(dòng)黏度；σk=1.0；Gk為由于平均速度梯度產(chǎn)生的湍流動(dòng)能；Gb為由于浮力產(chǎn)生的湍流動(dòng)能；YM為可壓縮介質(zhì)中波動(dòng)膨脹對(duì)總耗散率的貢獻(xiàn)；ε為湍流能量耗散率；σε=1.3；R=0；C1ε=1.44；C2ε=1.92；C3ε=0.

對(duì)于換熱器和冷卻裝置的固體區(qū)域，也滿足能量守恒，即

Tt（ρc）=·（λT）.（8）

1.2.3? 邊界條件

TEG系統(tǒng)的輸出性能受到尾氣質(zhì)量流量和溫度的影響較大，而尾氣質(zhì)量流量和尾氣溫度與發(fā)動(dòng)機(jī)、車輛質(zhì)量、負(fù)載量和循環(huán)工況等相關(guān).基于ADVISOR車輛仿真軟件平臺(tái)對(duì)重型卡車在高速公路燃油經(jīng)濟(jì)性測(cè)試（highway fuel economy test，HWFET）工況下進(jìn)行實(shí)況動(dòng)態(tài)仿真.重型卡車的型號(hào)為VEH_ralphs_grocery，質(zhì)量為7 068 kg，載質(zhì)量為2 000 kg，其發(fā)動(dòng)機(jī)型號(hào)為diesel engine，排量為7.2 L，最大功率為206 kW.由于排氣系統(tǒng)是從環(huán)境溫度慢慢上升，經(jīng)歷一段時(shí)間的循環(huán)后，尾氣溫度和尾氣質(zhì)量流量才達(dá)到平衡，所以選擇在第20個(gè)循環(huán)周期輸出重型卡車的瞬態(tài)尾氣溫度和尾氣質(zhì)量流量.經(jīng)仿真可得，重型卡車在第20個(gè)循環(huán)工況下，尾氣質(zhì)量流量和尾氣溫度隨時(shí)間的變化如圖4所示，尾氣質(zhì)量流量在0到40 g/s之間振蕩，尾氣溫度在590到630 K之間變化，且尾氣溫度的變化較平緩.

在瞬態(tài)CFD模型中，TEG系統(tǒng)與環(huán)境接觸的壁面上設(shè)置了對(duì)流換熱邊界，其中，自然對(duì)流換熱系數(shù)設(shè)置為15 W/（m2·K），環(huán)境溫度設(shè)置為300.00 K.在換熱器的入口處，采用ADVISOR得到的瞬態(tài)尾氣質(zhì)量流量和瞬態(tài)尾氣溫度作為尾氣流道入口的邊界條件.在冷卻裝置的入口處，采用穩(wěn)態(tài)溫度為363.15 K和穩(wěn)態(tài)流速為1 m/s的冷卻水作為冷卻水流道的入口邊界條件.在換熱器和冷卻裝置的出口處，均設(shè)置為壓力出口，且出口壓力為標(biāo)準(zhǔn)大氣壓.

1.3? TEM的瞬態(tài)分析模型

1.3.1? TEM的結(jié)構(gòu)

TEM是TEG系統(tǒng)的核心發(fā)電元件，單個(gè)TEM結(jié)構(gòu)如圖5所示.在TEM的實(shí)際工作過(guò)程中，尾氣熱量通過(guò)換熱器傳遞到TEM的熱端，冷卻水在TEM的冷端進(jìn)行散熱，TEM兩端產(chǎn)生溫差，導(dǎo)致內(nèi)部電子從熱端向冷端發(fā)生轉(zhuǎn)移形成電動(dòng)勢(shì).整個(gè)過(guò)程產(chǎn)生了焦耳熱、帕爾貼熱及湯姆遜熱.實(shí)際瞬態(tài)分析模型中，在TEM的兩端連接一個(gè)17 Ω的負(fù)載電阻以計(jì)算TEM的輸出功率.

1.3.2? 控制方程

在TEM的工作過(guò)程中，其熱端吸收的熱量和冷端散失的熱量分別為

Qh（t）=αI（t）Th_te（t）+

λteAtehteTh_te（t）-Tc_te（t）-12I2（t）Rin，（9）

Qc（t）=αI（t）Tc_te（t）+

λteAtehteTh_te（t）-Tc_te（t）+12I2（t）Rin，（10）

式中：α為熱電半導(dǎo)體的塞貝克系數(shù)；I（t）為電流；Rin為內(nèi)阻；Th_te（t）、Tc_te（t）分別為熱電半導(dǎo)體熱端和冷端的溫度.

I（t）=α（Th_te（t）-Tc_te（t））Rin+RL，（11）

式中：RL為負(fù)載電阻.

TEM熱端吸收的熱量和冷端散失的熱量還可以表示為

Qh（t）=Th（t）-Th_te（t）Rh_co+ce，（12）

Qc（t）=Tc_te（t）-Tc（t）Rc_co+ce，（13）

式中：Rh_co+ce為熱端銅電極片和陶瓷板的總熱阻； Rc_co+ce為冷端銅電極片和陶瓷板的總熱阻.

TEM的輸出功率和轉(zhuǎn)化效率分別為

PTEM（t）=Qh（t）-Qc（t）=I2（t）RL，（14）

ηTEM（t）=PTEM（t）Qh（t）.（15）

上述瞬態(tài)分析模型僅計(jì)算了單個(gè)TEM的輸出功率，將4個(gè)TEM的輸出功率相加可計(jì)算得到1/4 TEG系統(tǒng)的輸出功率，故1/4 TEG系統(tǒng)的輸出功率、吸收的尾氣熱量和轉(zhuǎn)化效率分別為

PTEG（t）=PTEM1（t）+PTEM2（t）+

PTEM3（t）+PTEM4（t），（16）

QTEG（t）=caqm（t）（Tin（t）-Tout（t）），（17）

ηTEG（t）=PTEG（t）QTEG（t），（18）

式中：PTEM1（t）、PTEM2（t）、PTEM3（t）、PTEM4（t）分別為TEM1、TEM2、TEM3、TEM4的輸出功率；ca為尾氣比熱容；qm（t）為尾氣質(zhì)量流量；Tin（t）為進(jìn)氣口尾氣溫度；Tout（t）為出氣口尾氣溫度.

1.3.3? 邊界條件

瞬態(tài)CFD模型得到1/4 TEG系統(tǒng)TEM1、TEM2、TEM3、TEM4的平均熱端溫度和平均冷端溫度隨時(shí)間的變化如圖6所示，其中：Th1（t）、Th2（t）、Th3（t）、Th4（t）分別為TEM1、TEM2、TEM3、TEM4的平均熱端溫度；Tc1（t）、Tc2（t）、Tc3（t）、Tc4（t） 分別為TEM1、TEM2、TEM3、TEM4的平均冷端溫度.將Th1（t）、Th2（t）、Th3（t）、Th4（t）、Tc1（t）、Tc2（t）、Tc3（t）、Tc4（t）作為瞬態(tài)分析模型的邊界條件.

從圖6可以看出，由于水的比熱容較高，相比于TEM熱端溫度的變化幅度，冷端溫度的變化較為平緩.從圖6a可以看出，時(shí)間為0～100 s，TEM3的熱端溫度高于TEM4，這是由于TEM4接近換熱器出口，換熱器出口是從大截面到小截面，從而產(chǎn)生回流，導(dǎo)致溫度較高.

2? 混合多物理場(chǎng)數(shù)學(xué)模型的建立

2.1? 網(wǎng)格獨(dú)立性驗(yàn)證

1/4 TEG系統(tǒng)的網(wǎng)格系統(tǒng)如圖7所示，在對(duì)混合瞬態(tài)CFD模型進(jìn)行有限元仿真分析時(shí)，網(wǎng)格靈敏度對(duì)于模型精度影響較大，需要對(duì)網(wǎng)格進(jìn)行獨(dú)立性分析.選擇TEM1熱端溫度的誤差作為判別網(wǎng)格精度的標(biāo)準(zhǔn)，選擇網(wǎng)格IV作為參照標(biāo)準(zhǔn)進(jìn)行誤差計(jì)算.混合瞬態(tài)CFD-分析模型的網(wǎng)格分析如表2所示，網(wǎng)格II的計(jì)算時(shí)間為270 min，網(wǎng)格數(shù)為751 672個(gè)，與網(wǎng)格IV的誤差僅差0.05%，在可接受范圍內(nèi).TEM1的熱端溫度隨著網(wǎng)格數(shù)的增加而降低，在減少仿真時(shí)間和滿足網(wǎng)格精度的前提下，選擇網(wǎng)格II作為混合瞬態(tài)CFD-分析模型的網(wǎng)格劃分標(biāo)準(zhǔn).

2.2? 瞬態(tài)試驗(yàn)驗(yàn)證

2.2.1? 試驗(yàn)臺(tái)架

為了驗(yàn)證混合瞬態(tài)CFD-分析模型結(jié)果的準(zhǔn)確性，設(shè)計(jì)了TEG系統(tǒng)的實(shí)物裝置，如圖8所示.搭建了TEG系統(tǒng)的瞬態(tài)試驗(yàn)臺(tái)架，如圖9所示，進(jìn)行了試驗(yàn)結(jié)果與模型結(jié)果的對(duì)比驗(yàn)證.在實(shí)際的瞬態(tài)性能測(cè)試中，使用熱風(fēng)機(jī)（型號(hào)為F1-R1055）提供熱空氣來(lái)模擬廢氣，熱風(fēng)機(jī)的最大功率為5 kW，空氣溫度和空氣流速可分別在0～350 ℃和0～5.5 m3/min的范圍內(nèi)進(jìn)行調(diào)節(jié).考慮到熱量的流失，實(shí)際的進(jìn)氣口溫度低于熱風(fēng)機(jī)設(shè)定值.選用熱電偶傳感器（型號(hào)為WRNT）測(cè)量入氣口處的熱空氣溫度并用溫度數(shù)據(jù)記錄儀（型號(hào)為RDXL4SD）進(jìn)行記錄，在出氣口處使用熱線風(fēng)速計(jì)（型號(hào)為HHF-SD1）測(cè)量實(shí)時(shí)的熱空氣流速.由于熱線風(fēng)速儀的工作溫度不能高于50 ℃，所以在TEG實(shí)物裝置后連接了一個(gè)中冷器，其作用是保護(hù)熱線風(fēng)速計(jì).其中熱風(fēng)機(jī)、TEG系統(tǒng)、中冷器和熱線風(fēng)速計(jì)之間由帶法蘭的鋁管連接.選用水冷方式進(jìn)行散熱，水的流速約為10 L/min.為了測(cè)量電壓，將電子負(fù)載（型號(hào)為IT8500+，ITECH）與TEG實(shí)物裝置相連形成閉合回路，并將阻值設(shè)置為17 Ω.在實(shí)際的瞬態(tài)試驗(yàn)中，穩(wěn)定工作條件下的水溫約為284.85 K.在仿真時(shí)采用了1/4簡(jiǎn)化模型，只測(cè)試了1/4 TEG系統(tǒng)的輸出電壓.

2.2.2? 試驗(yàn)結(jié)果對(duì)比

通過(guò)試驗(yàn)得到600 s內(nèi)的熱空氣溫度和熱空氣流速，將其輸入到瞬態(tài)CFD-分析模型中進(jìn)行求解，可得到1/4 TEG系統(tǒng)的輸出電壓.TEG系統(tǒng)瞬態(tài)模型預(yù)測(cè)與瞬態(tài)試驗(yàn)的輸出電壓對(duì)比如圖10所示，兩者的變化趨勢(shì)大致相同，時(shí)間為0～600 s，模型預(yù)測(cè)電壓比試驗(yàn)輸出電壓平均高6.41%，這主要是由于瞬態(tài)試驗(yàn)測(cè)量誤差造成的；其次，瞬態(tài)模型設(shè)置了假設(shè)條件，忽略了實(shí)際的熱量損失，模型預(yù)測(cè)電壓比試驗(yàn)輸出電壓有更大的響應(yīng)滯后，這是由于模擬中不能考慮儀器的熱慣性，即當(dāng)空氣溫度變化時(shí)，溫度傳感器受熱慣性的影響不會(huì)立即響應(yīng)，只能將延遲的信號(hào)記錄下來(lái)，由此在進(jìn)行數(shù)值計(jì)算時(shí)存在一定時(shí)滯.

3? 瞬態(tài)性能分析

3.1? 溫度分布

利用上述瞬態(tài)CFD模型研究了在HWFET循環(huán)工況下1/4 TEG系統(tǒng)的溫度場(chǎng)分布.當(dāng)時(shí)間t分別為0、200、400、600 s時(shí)，1/4TEM模塊和1/4 TEG系統(tǒng)的溫度分布如圖11所示.

從圖11可以看出：t=400 s時(shí)換熱器內(nèi)壁溫度最高，尾氣質(zhì)量流量和尾氣溫度最大，TEM兩端溫差也最大；t=200 s與t=600 s時(shí)的尾氣質(zhì)量流量幾乎相同，在尾氣溫度卻偏高的情況下，t=600 s時(shí)，熱電半導(dǎo)體的熱端溫度反而更高.由圖6b可以看出：在t=200 s時(shí)，尾氣溫度從低到高；在t=600 s時(shí)，尾氣溫度從高到低，溫度受熱慣性的影響沒有立即響應(yīng)，在較短時(shí)間內(nèi)仍然保持原來(lái)較高溫度狀態(tài)，從而導(dǎo)致溫度升高.

1/4TEG系統(tǒng)的溫度分布如表3所示.

［8］? LUO D， YAN Y Y， WANG R C， et al. Numerical investigation on the dynamic response characteristics of a thermoelectric generator module under transient temperature excitations［J］. Renewable Energy， 2021，170：811-823.

［9］? WANG S L， LIU H B， GAO Y W， et al. Transient supercooling performance of thermoelectric coolers with a continuous double current pulse［J］. Journal of the Taiwan Institute of Chemical Engineers， 2021，120：127-135.

［10］? LUO D， WANG R C， YU W， et al. Theoretical analysis of energy recovery potential for different types of con-ventional vehicles with a thermoelectric generator［J］. Energy Procedia， 2019，158：142-147.

［11］? KUMAR S， HEISTER S D， XU X， et al. Thermoelectric generators for automotive waste heat recovery systems part II： parametric evaluation and topological stu-dies［J］. Journal of Electron Mater， 2013，42（6）：944-955.

［12］? WANG Y C， DAI C S， WANG S X. Theoretical analysis of a thermoelectric generator using exhaust gas of vehicles as heat source［J］. Applied Energy， 2013，112：1171-1180.

[JY，2][HT5SS]（責(zé)任編輯? 賈國(guó)方）