黑迎曉，曲明璐，黃 萍，張童瑤，李 釗

（上海理工大學(xué) 環(huán)境與建筑學(xué)院，上海 200093）

0 引言

相變蓄能具有蓄熱密度大等諸多優(yōu)點，受到廣泛關(guān)注和研究。相變蓄熱器可應(yīng)用于復(fù)疊式空氣源熱泵中進行蓄能除霜［1-3］，在常規(guī)復(fù)疊式空氣源熱泵中加入相變蓄熱器，可將制熱模式下低溫級多余的熱能儲存到其中，除霜期間蓄熱器蓄存熱量用于低溫級除霜和高溫級供熱［4-5］。但是高低溫級制冷劑分別為R134a和R410A，而且各級的制冷劑流量、溫度也因各自不同的運行工況出現(xiàn)不同，因此相變蓄熱器會出現(xiàn)多溫位釋熱情況。QU等［6］針對此問題采取了5種不同的室外工況的實驗研究分析，結(jié)果表明高溫級的熱量供給占相變蓄熱器總釋熱量的57.4%～35.0%，低溫級占剩余的42.6%～65%，同時得出制冷劑流量是影響蓄熱器高低溫級釋熱量分配的主要因素。然而對于蓄熱器性能的研究，需要利用適當(dāng)數(shù)值模型對其分析。目前，對相變蓄熱的數(shù)值模擬研究主要包括蓄熱器形態(tài)結(jié)構(gòu)［7-9］，相變材料蓄放熱規(guī)律［10-11］，蓄熱器的優(yōu)化設(shè)計［12-14］，但未見針對多溫位蓄熱器釋熱特性的數(shù)值研究。因此，本文利用ANSYS FLUENT 軟件對簡化的三維蓄熱器單元模塊進行建模，研究多溫位釋熱特性。

1 試驗與模型建立

1.1 蓄熱器結(jié)構(gòu)

為驗證物理模型的準確性，需要利用試驗結(jié)果驗證，本文采用復(fù)疊式空氣源熱泵系統(tǒng)實驗進行實驗操作，實驗臺的結(jié)構(gòu)如圖1所示，系統(tǒng)各部件參數(shù)詳見文獻［5］。本研究中關(guān)鍵部件蓄熱器由外殼體和換熱器內(nèi)芯組成，其中換熱器內(nèi)芯為翅片管式換熱器，由高溫級盤管及低溫級盤管分排交錯排列構(gòu)成。圖2和表1分別示出蓄熱器的設(shè)計結(jié)構(gòu)以及蓄熱器的結(jié)構(gòu)參數(shù)。蓄熱器的外殼均包有保溫棉。

圖1 實驗臺結(jié)構(gòu)示意Fig.1 Schematic structural diagram of the lab bench

圖2 蓄熱器結(jié)構(gòu)Fig.2 Structural diagram of thermal energy storage

表1 蓄熱器結(jié)構(gòu)參數(shù)Tab.1 Structure parameters of thermal energy storage

相變材料RT10充滿在換熱器內(nèi)芯與外殼體之間。表2和圖3分別示出功率補償器（DSC）測得的RT10主要物性參數(shù)和相應(yīng)測試結(jié)果。試驗中，各層相變材料和高低溫級制冷劑溫度利用T型熱電偶（0.1F.S）進行測量。高低溫級制冷劑流量由轉(zhuǎn)子式流量計（精度1.6%）測量。

表2 RT10主要物性參數(shù)Tab.2 The main physical parameters of RT10

圖3 RT10的功率補償器（DSC）測試結(jié)果Fig.3 The DSC test result of RT10

1.2 試驗方案與結(jié)果

本試驗中室內(nèi)外環(huán)境溫濕度在蓄熱和放熱操作期間維持室內(nèi)側(cè)模擬工況干球溫度：（22±0.1）℃，相對濕度 50%±2%；室外側(cè)模擬工況干球溫度分別為-9，-15 ℃，相對濕度均75%±3%，除霜量控制在2.0 kg。根據(jù)試驗以及數(shù)據(jù)處理得到的試驗結(jié)果見表3，試驗步驟詳見文獻［3］。

表3 釋熱過程中熱量耦合特性參數(shù)及試驗結(jié)果Tab.3 Thermal coupling characteristic parameters and experimental results in the heat release process

1.3 數(shù)值模型建立

1.3.1 物理模型

本文采用ANSYS FLUENT建立該蓄熱器三維模型，尺寸選取與實驗系統(tǒng)中蓄熱器單元實際尺寸一致，蓄熱器三維模型如圖4所示。

圖4 蓄熱器三維模型Fig.4 3D model of the thermal energy storage

為了方便計算，對該物理模型做如下假設(shè)和設(shè)置：

（1）設(shè)置材料各物性參數(shù)為常數(shù)，僅考慮自然對流而引起的流體密度的變化；

（2）液相區(qū)中的流體為不可壓縮牛頓流體；

（3）蓄熱套管關(guān)于制冷劑軸線與兩管中線軸對稱；

（4）除了與制冷劑的耦合表面，相變材料的其他外表面均設(shè)置為絕熱；

（5）不考慮流體的黏性耗熱；

（6）其他物性參數(shù)為常數(shù)，相變材料密度滿足Boussinesq假設(shè)。

1.3.2 模型驗證

根據(jù)之前的研究發(fā)現(xiàn)，由于影響高低溫級釋熱的主導(dǎo)因素為制冷劑流量［6］，因此，本文換熱邊界條件設(shè)置為不同熱流密度邊界條件。

在室外干球溫度-9 ℃，相對濕度為75%，結(jié)霜量為2 kg（工況一）的試驗工況下對蓄熱器作為釋熱過程作模擬分析。設(shè)置相變材料初始溫度為21 ℃。由表3中的試驗結(jié)果可得，低溫級釋熱量為1 398.7 kJ，高溫級為1 343.9 kJ，該工況除霜時間為530 s。由表1中參數(shù)計算可得蓄熱器高、低溫級換熱面積均為3.8 m2。由式（1）計算得到試驗過程中高、低溫級平均熱流密度分別為667.3，694.4 W/m2。

式中 q ——平均熱流密度，W/m2；

Q ——釋熱量，J；

S ——換熱面積，m2；

t ——除霜時間，s。

利用ANSYS FLUENT數(shù)值模擬蓄熱器高溫級相變材料溫度與工況1下所測得的相變材料平均溫度進行對比，如圖5所示，該模擬工況下相變材料初始溫度為294 K（21 ℃），圖中可以看出相變材料溫度曲線較為吻合，試驗所測相變材料平均溫度為281.5 K（8.4 ℃），模擬相變材料平均溫度為281.1 K（8.0 ℃），計算得誤差為4.8%。因此，該模型及相應(yīng)參數(shù)設(shè)置相對準確，可用于接下來的分析研究。

圖5 蓄熱器相變材料溫度試驗、模擬結(jié)果對比Fig.5 Simulation comparison with experiments chart of PCM temperature on thermal energy storage

2 模擬結(jié)果分析

將驗證后的該數(shù)值模型應(yīng)用于該蓄熱器在不同熱流密度邊界條件下的釋熱過程的模擬，分析對應(yīng)的釋熱過程中蓄熱器高、低溫級相變材料液相率變化及其溫度變化情況，具體邊界條件見表4。第1組工況即為驗證模擬設(shè)定的工況一，第2組工況根據(jù)表3試驗工況二計算可得，低溫級釋熱量約為高溫級釋熱量的2.1倍，故在此設(shè)置模擬邊界條件為低溫級恒定熱流密度值約為高溫級的2.1倍。為和第2組進行對照，第3組工況設(shè)置高溫級恒定熱流密度值約為低溫級的2.0倍。各組相變材料與紫銅管的初始溫度均設(shè)為21 ℃。

表4 蓄熱器邊界條件Tab.4 Thermal energy storage boundary conditions

2.1 液相率分析

3組不同工況下蓄熱器高低溫級相變材料在200 s時刻液相率如圖6所示。200 s時蓄熱器高溫級相變材料液相率分別為0.71，0.82，0.46，低溫級分別為0.68，0.57，0.85。同一時刻，第1組工況下高低溫級相變材料液相率值很接近，高溫級相變材料液相率略高于低溫級液相率值，說明高溫級相變材料相變速率略低于低溫級相變速率。第2組工況下，同一時刻，高溫級相變材料液相率值遠高于低溫級的數(shù)值，說明高溫級相變材料相變速率遠低于低溫級相變速率。第3組工況下，同一時刻，高溫級相變材料液相率值低于低溫級的數(shù)值，說明高溫級相變材料相變速率高于低溫級相變速率。圖7示出了3組工況下蓄熱器在除霜期內(nèi)高低溫級液相率隨時間變化情況。

圖6 各工況200 s時蓄熱器內(nèi)相變材料液相率Fig.6 Liquid fraction diagram of PCM in the thermal energy storage at 200 seconds for each operating condition

圖7 3組工況相變材料液相率變化Fig.7 Change in liquid fraction of PCM in the three groups of conditions

從圖7可以看出，第1組工況下，高低溫級液相率曲線幾乎重合，原因是此工況下高低溫級熱流密度差別很小。第2組工況下，蓄熱器高低溫級相變材料液相率差值隨時間變化逐漸增大，低溫級相變材料液相率值很快由1降至接近0，而高溫級相變材料液相率值在500 s時只降低至0.5左右，高低溫級相變差距明顯，說明該工況下，蓄熱器低溫級很快完成了相變，而高溫級相變速率緩慢。第3組工況，蓄熱器高低溫級相變材料液相率差值隨時間變化逐漸增大，低溫級液相率500 s時降至0.4，高溫級液相率值在400 s時已基本完成相變。該工況下，蓄熱器高溫級相變材料很快完成了相變，而低溫級相變速率則很緩慢。由圖可見，第2組工況低溫級率先完成相變，向低溫級釋熱更快，即可以更快完成低溫級除霜。

2.2 溫度分析

3組不同工況下蓄熱器內(nèi)高低溫級相變材料200 s時刻溫度對比如圖8所示。第1組工況下，蓄熱器高溫級相變材料溫度在200 s為7.6 ℃，低溫級相變材料溫度為11.1 ℃。同一時刻，蓄熱器高低溫級相變材料溫度接近。第2組工況下，蓄熱器高溫級相變材料溫度200 s為，7.8 ℃，低溫級相變材料溫度為9.1 ℃，高溫級釋熱速率低于低溫級釋熱速率，當(dāng)?shù)蜏丶壔就瓿上嘧儠r，高溫級相變材料大多仍處于相變過程中。第3組工況下，蓄熱器高溫級相變材料溫度在200 s為6.9 ℃，低溫級相變材料溫度為7.9 ℃，大約40 s后，高溫級相變材料到達凝固點，大約100 s后，低溫級相變材料到達凝固點。第3組工況則與第2組工況相反，高溫級釋熱速率高于低溫級釋熱速率，當(dāng)?shù)蜏丶壢栽谙嘧儠r，高溫級已基本完成相變。

圖8 各工況200 s時蓄熱器內(nèi)相變材料液溫度Fig.8 Temperature diagram of PCM in the thermal energy storage at 200 s for each operating condition

圖9示出3組工況下蓄熱器在500 s內(nèi)高低溫級相變材料溫度隨時間變化情況圖。

圖9 3組工況相變材料溫度變化Fig.9 Change in temperature of PCM in the three groups of conditions

從圖9中看出，第1組工況下，蓄熱器高低溫級相變材料溫度曲線幾乎重合。第2組工況下，蓄熱器高低溫級相變材料溫度差值隨時間變化而增大，高溫級相變材料溫度在140 s后開始緩慢降低，低溫級相變材料釋熱，0～60 s內(nèi)迅速降低，隨后進入相變過程，溫度緩慢降低，210 s相變結(jié)束后，溫度再次快速下降。第3組工況下，蓄熱器低溫級相變材料溫度在110 s后隨時間變化緩慢降低，高溫級相變材料溫度在0～40 s內(nèi)迅速降低，隨后進入相變過程，溫度緩慢降低，200 s后相變結(jié)束，溫度再次快速下降。該工況下，蓄熱器低溫級釋熱速度緩慢，高溫級釋熱則比較迅速?？梢?，第1組工況高低溫級相變材料釋熱時間基本一致，而第2組工況低溫級相變材料釋熱快，第3組高溫級相變材料釋熱快，這與液相率分析的結(jié)果相一致。

3 結(jié)論

（1）第1組高低溫恒定熱流密度接近工況下，高低溫級液相率以及溫度變化曲線在釋熱過程最為重合，表明工況下高低溫級液相率變化速率最為一致、溫度變化最為均勻。

（2）第2組低溫級熱流密度約為高溫級的2.1倍工況下，低溫級的液相率與溫度曲線變化相對較大，表明該工況下低溫級液相率和溫度下降更迅速。

（3）第3組高溫級熱流密度值為低溫級的2.0倍工況下，高溫級的液相率與溫度曲線變化相對較大，表明該工況下高溫級液相率和溫度下降更迅速。

本文所建立的模型能夠針對相變蓄熱器多溫位釋熱情況進行預(yù)測和分析，并可以進一步應(yīng)用于多溫位釋熱相變蓄熱器的結(jié)構(gòu)優(yōu)化研究。