謝云云 阮如君

(桂林航天工業(yè)學(xué)院 能源與建筑環(huán)境學(xué)院，廣西 桂林 541004)

在我國經(jīng)濟(jì)穩(wěn)步發(fā)展過程中，能源作為主要推動力，是保證我國經(jīng)濟(jì)發(fā)展規(guī)模以及速度的基石[1]。目前世界能源儲備并不能滿足全世界人口的要求，傳統(tǒng)能源的開發(fā)與利用對環(huán)境造成不可逆的破壞，由于環(huán)境問題的日益突出與嚴(yán)重， 人們環(huán)保和節(jié)能意識的增強(qiáng)，研究開發(fā)新型環(huán)保的能源利用技術(shù)成為當(dāng)下各國研究者關(guān)注的焦點。目前，我國在積極開發(fā)利用太陽能、風(fēng)能、水能等綠色能源的同時，也在積極探索高效環(huán)保的能源利用技術(shù)，其中，相變儲熱是目前較為受關(guān)注的一項技術(shù)。相變儲熱材料在發(fā)生相變時具有非常大的相變潛熱，同時體積和溫度變化非常小，通過介質(zhì)在特定場合、特定時間以及特定溫度下發(fā)生相變，存儲或釋放大量的潛熱來完成相應(yīng)的工作。不同相變材料的相變溫度從低溫到高溫均有分布，因此已經(jīng)在很多領(lǐng)域得到應(yīng)用，比如在太陽能發(fā)電、工業(yè)余熱的回收、供熱供電、建筑材料中提供溫度調(diào)節(jié)等。關(guān)于冰—水相變換熱的研究，在凍土熱融方面有著很大的應(yīng)用價值[2]。對相變儲熱過程以及相變換熱器的研究分析，有利于提高能源利用率，具有重要的工程實際意義[3]。

換熱器作為相變儲熱裝置的主要組成部分，可以在供需的兩端進(jìn)行能量的傳遞和轉(zhuǎn)移，使需求一方的使用可以得到保障，隨著相變材料研究不斷深化以及其在眾多領(lǐng)域中得到廣泛利用，換熱器在許多相變儲能的項目中發(fā)揮了重要的樞紐作用[4]。目前，將相變儲能技術(shù)與換熱技術(shù)相結(jié)合，不斷開發(fā)出一種新型的儲能相變換熱器已成為研究熱點[5]。李洋等[6]對管殼式、板式和熱管式等相變換熱的優(yōu)缺點進(jìn)行了定性分析，隨后通過數(shù)值模擬方式，對板式和管式相變換熱器的二維相變?nèi)刍Ｐ瓦M(jìn)行了定量比較。發(fā)現(xiàn)管式相變換熱器的融化時間為6 h，板式相變換熱器則為8.5 h，這主要是由于兩種換熱器的換熱管/板的排布方式的差異導(dǎo)致，考慮到板式換熱器造型緊湊、加工簡單的優(yōu)點，因此具有巨大的優(yōu)勢。郭茶秀等[7]則是對板式相變換熱器的強(qiáng)化換熱進(jìn)行了相應(yīng)的研究，通過數(shù)值模擬，采用石蠟作為相變材料，通過在換熱板添加翅片的方式對板式換熱器進(jìn)行模擬計算，得到了不同時間的相界面位置、凝固總時間、壁面熱流、翅片溫度分布等參數(shù)。此外，還研究了翅片在長寬比參數(shù)不同的蓄熱器中的傳熱作用。尹點等[8]則是采用實驗和數(shù)值模擬的方式，設(shè)計了一種板式換熱器，通過非穩(wěn)態(tài)儲熱模擬計算以及實驗數(shù)據(jù)進(jìn)行參照對比，對該板式換熱器的性能進(jìn)行了分析研究，同時對相變材料層以及熱流層結(jié)構(gòu)參數(shù)的變化對于相變過程的影響也進(jìn)行了探究，發(fā)現(xiàn)相變材料的儲能速率與儲能量會受到相變材料層和熱流層厚度的影響，同時，儲能速率與儲能量也是呈反比。KUMAR等[9]對板式換熱器的水力和熱力約束進(jìn)行了優(yōu)化，以提高換熱器的靈敏度。 在所提出的方法中，多目標(biāo)列優(yōu)化(MOWO)用于優(yōu)化板式換熱器的參數(shù)。在MOWO中，板式換熱器的水平端口中心距、垂直端口中心距、放大系數(shù)、端口直徑、板厚、熱板數(shù)量和板間距等參數(shù)進(jìn)行了優(yōu)化，以獲得更好的靈敏度。以熱壓降和水壓降為目標(biāo)函數(shù)，使用MOWO增強(qiáng)傳熱，使壓降最小化。然后利用matlab實現(xiàn)了所提出的系統(tǒng)并對其性能進(jìn)行了分析發(fā)現(xiàn)，所獲得的性能證明了所提出系統(tǒng)的有效性。BOBIC等[10]在為現(xiàn)代節(jié)能區(qū)域供熱和制冷應(yīng)用開發(fā)智能控制算法時，了解板式換熱器的動態(tài)特性起著至關(guān)重要的作用。在這項工作中，研究了逆流板式換熱器的動態(tài)行為，以了解其在不同流體流動配置下由于入口溫度擾動引起的瞬態(tài)響應(yīng)。為此，提出了一個合適的理論模型。通過求解代表所研究板式換熱器的派生一維模型的集總參數(shù)系統(tǒng)，對溫度瞬態(tài)進(jìn)行數(shù)值評估。數(shù)值結(jié)果表明，當(dāng)流體從板式換熱器的入口流向熱的出口時，不同的流體流動配置會顯著影響瞬態(tài)溫度響應(yīng)以及整體溫降和傳熱速率。流體側(cè)通過使用紅外熱成像技術(shù)進(jìn)行一系列系統(tǒng)測試，對預(yù)測的瞬態(tài)進(jìn)行了實驗驗證。分析表明，所提出的模型與實驗吻合較好。熱成像測量的結(jié)果還可以更深入地了解溫度分布及其沿流體流動通道的瞬態(tài)前沿傳播。在未來的研究中，通過紅外熱成像這種實驗方法對于在預(yù)測釬焊板式換熱器的熱疲勞壽命時識別臨界溫度區(qū)以及擴(kuò)展所提出的理論模型以包括獲得的二維局部傳熱系數(shù)分布很有價值。BHATTAD等[11]使用不同濃度的混合納米流體(Al2O3+ MWCNT/水)對板式換熱器進(jìn)行了數(shù)值和實驗研究，通過改變流體濃度研究其對傳熱和壓降特性的影響。在數(shù)值模擬研究中，使用離散相位模型，并將結(jié)果與實驗結(jié)果以及同質(zhì)模型的結(jié)果進(jìn)行了比較。研究了不同操作參數(shù)(納米流體入口溫度、流速和體積濃度)對冷卻劑出口溫度、傳熱率、對流和總傳熱系數(shù)、努塞爾數(shù)、摩擦系數(shù)、壓降、泵送功率、效率的影響和績效指數(shù)。除此之外還研究了基礎(chǔ)流體、納米流體和混合流體的速度和溫度分布。通過使用混合納米流體，傳熱系數(shù)提高了39.16%(優(yōu)點)，而泵送功率的增加可以忽略不計1.23%(缺點)。在使用混合納米流體代替基礎(chǔ)流體時，由于傳熱和壓降特性的增強(qiáng)，這也體現(xiàn)了板式換熱器的有效性。季中敏等[12]將相變換熱混合工質(zhì)低溫板翅式換熱器的表面?zhèn)鳠崤c流阻特性進(jìn)行數(shù)值模擬。實驗在給定溫度下沿縱向模擬得到了傳熱系數(shù)和壓力梯度的變化曲線，并將其與MUSE換熱器設(shè)計仿真軟件的計算結(jié)果進(jìn)行了比較。

綜上所述，板式相變換熱器儲熱模塊具有比表面積大、厚度方向?qū)嵝院?、加工過程簡單、組配靈活、結(jié)構(gòu)緊湊等優(yōu)點，為后期維護(hù)提供了極大便利[13]。目前，板式換熱器在電子零部件熱管理、建筑溫控節(jié)能、凍土熱融管理以及太陽能熱利用等領(lǐng)域具有非常好的應(yīng)用前景。由于相變換熱器儲熱結(jié)構(gòu)特點對換熱過程以及換熱流體產(chǎn)生影響，使得換熱器的儲熱性能也受到影響，性能不穩(wěn)定。本文采用控制變量法，控制相變換熱器的體積不變，通過改變換熱器入口換熱流體流速以及溫度的方法，利用Fluent軟件分別進(jìn)行數(shù)值模擬計算，并研究不同工況下?lián)Q熱器的儲熱性能。通過本文的研究結(jié)果，為后期對設(shè)計和優(yōu)化相變換熱裝置、促進(jìn)儲熱技術(shù)的發(fā)展提供一定的指導(dǎo)作用[14]。

1 數(shù)理模型與驗證

1.1 物理模型

與普通的換熱器相比，相變換熱器的特點在于換熱器中需要同時布置換熱流體的管道和相變材料[15]。本文設(shè)計了一種單板式相變換熱器模型，如圖1所示。換熱器模型長度l=280 mm，寬度w=40 mm。存儲相變材料單元長度=168.56 mm，寬度=20.96 mm。模型利用熱空氣作為換熱流體，冰作為相變材料。左端為換熱流體的入口，右端為換熱流體的出口，其中換熱流體熱空氣流經(jīng)管道與儲存在方形管內(nèi)的相變材料冰進(jìn)行熱交換。

圖1 單板式相變換熱器模型示意圖(單位：mm)

考慮到相變換熱過程的復(fù)雜性，本文在對板式換熱器進(jìn)行模擬計算時，進(jìn)行了相應(yīng)的簡化，并做出相應(yīng)假設(shè)：

1)換熱器內(nèi)的流體流動過程中，不存在黏性耗散；

2)相變材料為各向同性；

3)忽略儲熱單元壁厚；

4)相變過程形成的糊狀區(qū)內(nèi)溫度呈線性變化；

5)不考慮相變材料熔化前后的體積變化；

6)相變材料的密度使用Boussinesq假設(shè)[16]。

本文所采用的相變材料(PCM)冰和換熱流體(HTF)空氣的物性參數(shù)如表1所示。

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

1.2 數(shù)學(xué)模型

換熱流體區(qū)的控制方程：

1)質(zhì)量守恒方程

(1)

2)動量守恒方程

x方向：

(2)

y方向：

(3)

能量守恒方程：

(4)

式中：

uf——空氣在x方向的速度，m/s；

vf——空氣在y方向的速度，m/s；

ρf——空氣密度，kg/m3；

cf——空氣的比熱，kJ/(kg·K)；

μf——空氣的動力黏度，kg·s/m2；

g——重力加速度，m/s2；

kf——空氣的導(dǎo)熱系數(shù)，W/(m·K)；

Tf——空氣溫度，K；

t——時間，s。

相變區(qū)的守恒方程：

1)質(zhì)量守恒方程

(5)

2)動量守恒方程

x方向：

(6)

y方向：

(7)

能量守恒方程：

(8)

式中：

u——液態(tài)PCM在x方向的速度，m/s；

v——液態(tài)PCM在y方向的速度，m/s；

ρ——PCM的密度，kg/m3；

μ——液態(tài)PCM的動力黏度，kg·s/m2；

β——PCM的熱膨脹系數(shù)，K-1；

k——PCM的導(dǎo)熱系數(shù)，W/(m·K)；

T0——Boussinesq假設(shè)的工作溫度，K；

H——金屬圓管或PCM的焓值，J。

模型中使用Boussinesq假設(shè)對相變材料內(nèi)部自然對流進(jìn)行模擬研究[16]?；赑ressure-Velocity Coupling算法，壓力模型為PRESTO！，動量方程為Second Order Upwind，能量方程為Second Order Upwind。

初始溫度：

T(x，y，t)=T0=265 K

(9)

入口邊界條件(x=0)：

Tair(x=0，t)=Tin=315K

(10)

u(x=0，t)=uin

(11)

v(x=0，t)=0

(12)

出口邊界條件(x=280 mm)：

(13)

(14)

(15)

對稱邊界條件(y=0，y=40 mm)：

(16)

(17)

(18)

1.3 網(wǎng)格無關(guān)性驗證

在網(wǎng)格無關(guān)性驗證中，劃分了三組不同網(wǎng)格數(shù)的模型進(jìn)行模擬計算，網(wǎng)格數(shù)分別為67 295、70 732和125 918。如圖2所示，三種網(wǎng)格數(shù)所模擬的出口溫度隨時間變化的曲線基本重合，為節(jié)省計算資源和時間，所以本文采用的換熱模型所劃分的網(wǎng)格數(shù)為70 732。

圖2 不同網(wǎng)格數(shù)模型出口溫度隨時間變化曲線

2 研究結(jié)果與分析

2.1 相變區(qū)溫度隨時間的變化

基于模型中提出的假設(shè)，在平板換熱器入口段、中間換熱區(qū)以及出口段選取三個特征點進(jìn)行研究，分別為T1(65，20)、T2(168.92，20)和T3(203.56，20)。圖3給出了單板式換熱器相變區(qū)溫度隨時間的變化曲線。三條溫度變化曲線都有一個共同特征：在某一時間段，溫度幾乎保持水平不發(fā)生變化，這也表明此時材料正處于相變階段，而相變過程的特點就是溫度保持不變，狀態(tài)發(fā)生變化。相變材料的融化分三個階段，第一階段是相變材料在固體狀態(tài)下吸熱升溫，直到達(dá)到相變材料的相變溫度后，進(jìn)入第二階段，此時，相變材料具有非常大的相變潛熱，材料發(fā)生狀態(tài)變化時，從外界吸收大量的熱量但是溫度保持不變；相變材料融化完成，繼續(xù)對其加熱，此時相變材料吸熱升溫。以T1為例，在初始階段(0～880 s)，相變材料的溫度隨著時間不斷升高，直至T1點開始融化；第二階段(880～3 270 s)，該過程為T1點的融化過程，T1點相變材料的溫度始終保持在273.15 K左右，直至完全融化；第三階段(3 270～13 040 s)，在該過程中，相變材料完全融化以后，溫度繼續(xù)隨著時間的變化而升高，直至達(dá)到和換熱流體接近的溫度，315 K。從圖中該可以看出，T1融化所消耗的時間相對于T2和T3的較短，T1點相變材料融化所需的時間為3 270 s，T2點相變材料融化所需的時間為5 190 s，T3點相變材料融化所需的時間為4 850 s，由于T2和T3點較為接近，所以融化時間也相對接近。由此可知，換熱流體熱空氣的溫度沿著流動的方向不斷降低，更靠近換熱流體熱空氣入口處的相變材料所需的融化時間較短，傳熱速率也相對較快。

圖3 單板式相變換熱器相變區(qū)溫度隨時間變化曲線

2.2 液相分?jǐn)?shù)隨時間的變化

圖4為單板式相變換熱器相變區(qū)不同點處液相分?jǐn)?shù)隨時間變化曲線，坐標(biāo)分別為h1(65，20)、h2(168.92，20)和h3(203.56，20)，與溫度監(jiān)測點位置一致。由圖看出，隨著時間的變化，液相分?jǐn)?shù)也逐漸增大。該過程分為三個階段，以h1點為例，第一階段(0～3 180 s)，為相變材料開始升溫的過程，液相分?jǐn)?shù)不發(fā)生改變，為0，直至h1點開始融化；第二階段(3 180～3 310 s)，該過程為相變材料的融化階段，h1到達(dá)相變的臨界點，發(fā)生相變，隨之液相分?jǐn)?shù)呈上升趨勢，直到完全融化，液相分?jǐn)?shù)為1；第三階段(3 310～13 040 s)，相變材料完全融化，為升溫過程，液相分?jǐn)?shù)也不再發(fā)生改變。從圖還可以看出，h1點的融化時間比h2點和h3點所需消耗的時間要短，表明了越靠近換熱流體熱空氣入口處的相變材料完全融化所需的時間更短，傳熱效率更高，而靠近換熱流體出口處的相變材料融化所消耗的時間較長。

圖4 液相分?jǐn)?shù)隨時間變化曲線

2.3 固液界面和溫度云圖的變化

圖5為單板式相變換熱器五個不同時刻相變材料融化過程中固液分界和溫度的云圖。由圖5(a)固液分界云圖可知，沿著空氣流速的方向，隨著換熱流體熱空氣的溫度逐漸降低，相變區(qū)的相變材料冰的融化速率也明顯降低。由圖5(b)溫度分布云圖可知，沿著換熱流體流速的方向，溫度不斷降低，換熱速率也逐漸減慢。這表明了，隨著時間的變化，靠近換熱流體熱空氣入口側(cè)的相變材料傳熱速率較快，融化所需消耗的時間也較短；沿著換熱流體流動的方向，溫度逐漸降低，傳熱速率也隨之減慢，融化所消耗的時間較長。

圖5 固液分界和溫度云圖

2.4 入口熱空氣流速和溫度對冰融化過程的影響

不同入口邊界條件下，選取監(jiān)測點依次為P1(65，20)，P2(168.92，20)和P3(203.56，20)，圖6中(a)、(b)、(c)分別為三個監(jiān)測點在不同入口邊界條件下溫度隨時間的變化曲線。由圖可知，當(dāng)入口溫度和入口流速都較大時，管內(nèi)的相變材料融化所需要的時間明顯比其他條件下的融化時間短；同時，當(dāng)入口溫度保持不變時，流速越大時，相變材料的融化速度越快；同理，當(dāng)流速保持不變時，入口溫度越高也能加速相變材料的融化。在入口邊界條件為入口溫度T=315 K，入口流速v=2.5 m/s和入口溫度T=315 K，入口流速v=2 m/s的換熱模型中，兩者入口溫度相同，但入口流速不同，前者P2點的相變材料總?cè)诨臅r間為5 190 s，后者P2點的相變材料總?cè)诨臅r間為5 640 s。此時在入口邊界條件為入口溫度T=315 K，入口流速v=2.5 m/s和入口溫度T=305 K，入口流速v=2.5 m/s的換熱模型中，兩者入口流速相同，但入口溫度不同，前者P2點的相變材料總?cè)诨臅r間為5 190 s，后者P2點的相變材料總?cè)诨臅r間為6 600 s。

圖6 不同入口邊界條件的相變區(qū)溫度隨時間的變化曲線

從圖6(a)和(c)可知，P1和P3的相變區(qū)的溫度隨時間變化趨勢也和P2的變化趨勢大致相同。這表明了在相同幾何參數(shù)的換熱模型中，入口溫度一定時，通過提高入口流速，可以提高傳熱效率；入口流速一定時，通過提高入口溫度，可以更進(jìn)一步加快換熱速率；當(dāng)同時提高入口溫度和入口流速時，模型中的相變材料融化所需要消耗的時間會越來越短，傳熱速率變得越來越快。由圖6的變化曲線可發(fā)現(xiàn)，當(dāng)換熱流體溫度為315 K，流速為2.5 m/s時，相變換熱器內(nèi)的相變材料融化最快，同時，對于換熱流體溫度較高的換熱器模型，相變材料的融化明顯快于換熱流體溫度低的模型，這也表明入口溫度相對入口流速而言，通過提高入口溫度的換熱模型，換熱效果更好。

3 結(jié)論

采用空氣作為換熱流體、冰作為相變材料，對單板式相變換熱器進(jìn)行數(shù)值模擬計算。研究分析了換熱器內(nèi)相變過程的變化規(guī)律；并通過改變相變換熱器的入口溫度和入口流速，設(shè)定了四組不同的入口邊界條件，研究了傳熱過程中相變區(qū)同一點的溫度變化、液相分?jǐn)?shù)變化，以及溫度云圖變化等，分析相變材料的變化規(guī)律以及換熱流體與相變材料的換熱過程。研究發(fā)現(xiàn)：

1)在單板式相變換熱器傳熱過程中，靠近換熱流體入口處的相變材料溫度較高，融化較快。沿著換熱流體空氣流動的方向，溫度降低，傳熱速率減慢；

2)入口邊界條件對傳熱速率有一定的影響，在換熱器幾何參數(shù)相同的情況下，提高換熱流體的入口溫度或流速，傳熱速率可以得到相應(yīng)的提高，其中提高入口溫度，換熱效果更明顯；

3)在換熱器幾何參數(shù)相同的情況下，同時提高換熱流體的入口溫度和入口流速可以有效地提高換熱器的傳熱速率，并很大程度地縮短融化的總時間。