DOI：10.19912/j.0254-0096.tynxb.2021-1322 文章編號：0254-0096（2023）03-0525-08

摘 要：通過焓-孔隙法建立數(shù)值模型，研究石蠟RT50在管殼式太陽能儲熱器的傳熱特性，并分析儲熱器安裝形式、不同管徑和管外肋片結(jié)構(gòu)對相變儲熱過程的影響規(guī)律。結(jié)果表明，立式儲熱器在管徑較大時（相變材料區(qū)域厚度較薄時）熔化速度比臥式儲熱器快；而在管徑較小時（相變材料區(qū)域較厚時）熔化速度比臥式儲熱器慢。臥式儲熱器相變材料熔化過程的初始時間段和末尾時間段熱傳導占主導地位，而中間時間段自然對流換熱占主導地位，且熔化速度主要由自然對流換熱決定。此外，管外側(cè)安裝肋片能夠顯著提高傳熱系數(shù)，縮短熔化時間。相較于無肋片管道，臥式儲熱器采用環(huán)肋管道可縮短熔化時長31.6%；采用直肋管道可縮短熔化時長42.1%。直肋比環(huán)肋具有更好的傳熱效果，主要原因為直肋能夠強化熔化過程中的自然對流。

關(guān)鍵詞：太陽能；相變材料；熔化；肋片；自然對流；熱能儲存

中圖分類號：TK513.5 " "文獻標志碼：A

0 引 言

近年來太陽能的開發(fā)利用成為國內(nèi)外的研究熱點。由于太陽能分布不均、能流密度低以及間歇性日照的特點，儲能環(huán)節(jié)成大多數(shù)太陽能應(yīng)用過程中的關(guān)鍵［1］。國家能源局在2016年啟動的中國第一批太陽能熱發(fā)電示范項目中明確要求太陽能儲熱設(shè)備必須具備最少4 h的儲熱能力?，F(xiàn)有的儲能方式主要分為3種：顯熱儲能、潛熱儲能、熱化學儲能。其中，潛熱儲能是利用相變材料（phase change material， PCM）發(fā)生相變時吸熱或放熱的特性儲存或釋放能量，也叫相變儲熱。相比之下，潛熱儲能能量儲存密度較大、化學穩(wěn)定性好、應(yīng)用溫度范圍較廣，被認為是最有應(yīng)用前景的儲能方式之一［2］。

管殼式太陽能相變儲熱系統(tǒng)由于結(jié)構(gòu)簡單，應(yīng)用十分廣泛。近年來，眾多學者通過理論、實驗和數(shù)值等方法對管殼式相變儲熱系統(tǒng)傳熱性能進行了大量研究。Lacroix［3］建立了一個理論模型來預(yù)測帶有PCM管殼式儲熱系統(tǒng)的瞬態(tài)行為，結(jié)果表明PCM的導熱系數(shù)低是潛熱儲能系統(tǒng)的一個基本缺點，這導致PCM系統(tǒng)的低傳熱速率和延遲響應(yīng)。劉洋等［4］通過實驗研究了不同熱源溫度與環(huán)境溫度下儲能散熱過程中儲能單元的熱性能，結(jié)果發(fā)現(xiàn)添加泡沫銅對儲能單元的儲散熱具有較好的優(yōu)化效果。

對于相變儲熱來說，通過數(shù)值建模的方式，易于優(yōu)化容器的尺寸和形狀以及評估儲熱系統(tǒng)的性能，且無需大量的研究經(jīng)費［5］。張仲彬等［6］建立相變蓄熱膠囊的三維有序及無序堆積數(shù)值模型，發(fā)現(xiàn)當采用有序堆積時系統(tǒng)的各項性能指標最優(yōu)。李椿等［7］建立殼管式相變儲能換熱器仿真模型，研究了蓄熱過程中熱源管壁溫、熱管口徑規(guī)格對傳熱效率的影響。

與此同時，也有一些學者研究了PCM流動及肋片強化對相變儲熱過程的影響。Hosseini等［8］采用實驗和數(shù)值相結(jié)合的方法研究了管殼式換熱器內(nèi)石蠟RT50熔化過程中浮力驅(qū)動對流的作用，結(jié)果表明適當提高進水溫度可有效減少傳熱時間。Pahamli等［9］通過數(shù)值模擬的方法研究了偏心度、流量大小以及進口水溫等參數(shù)對PCM熔化過程的影響。田楊等［10］數(shù)值模擬研究不同轉(zhuǎn)速下加肋旋轉(zhuǎn)的蓄熱器的蓄熱特性，研究表明加肋并旋轉(zhuǎn)可顯著提升相變蓄熱器的蓄熱效果，縮短蓄熱時間。Sciacovelli等［11］使用不同的樹形肋片來提高管殼式潛熱儲能裝置的性能，同時提出一種改進瞬態(tài)工況系統(tǒng)的形狀優(yōu)化策略。

本文通過焓-孔隙法建立數(shù)值模型，研究不同形式太陽能儲熱器（立式和臥式）、不同管徑及管外加裝不同肋片等參數(shù)對傳熱過程PCM傳熱熔化行為的影響，進而闡明相變儲熱過程流動傳熱強化機理，為相變儲熱過程強化設(shè)計提供參考。

1 物理及數(shù)值模型

1.1 物理模型

1.1.1 立式、臥式換熱器儲熱單元

目前大多數(shù)國內(nèi)外對管殼式儲能換熱器的研究主要集中在設(shè)計及運行參數(shù)，對于儲熱器形式的影響研究缺乏對比。因此本文針對不同形式的管殼式儲能換熱器對相變儲熱過程的影響進行探究。管殼式換熱器儲熱單元形狀如圖1所示，其整體長度為600 mm，根據(jù)擺放形式分為臥式（圖1a）和立式（圖1b）。重力作用為豎直向下。傳熱流體（heat transfer fluid，HTF）為水工質(zhì)，走管程，其中立式為下口進，上口出。截面尺寸如圖1c所示。殼側(cè)填充的PCM材料為石蠟RT50。表1為RT50的各項熱物性參數(shù)。

RT50黏度[μ]與其溫度[T]的關(guān)系為式（1）：

式中：[A=0.819]，[B=-0.01546]。

1.1.2 翅片形狀和尺寸

對于管殼式換熱器，安裝肋片增加換熱面積以達到強化傳熱的目的是以往研究中常用的方法。尤其對于無機鹽類以及有機相變材料來說，其導熱系數(shù)通常較小，這也是相變儲熱應(yīng)用中的短板之一。以往研究表明，肋片管能有效提高相變過程中的傳熱效率。有學者對肋片的高度間距進行探究，也有學者設(shè)計了各類肋片結(jié)構(gòu)以提高傳熱效率。但大多新型肋片設(shè)計較為復(fù)雜，難以得到廣泛應(yīng)用；而對于常見的一些工程應(yīng)用中的肋片類型，其傳熱強化效果缺乏定性定量的認識。本文在維持可力片表面積與體積不變的前提下，分別研究直肋與環(huán)肋對相變材料熔化過程的強化作用，并進行對比分析。為了便于比較，在圖1a臥式儲熱器管上加裝直肋與環(huán)肋，其形狀尺寸分別如圖2、圖3所示。

1.2 數(shù)值模型

1.2.1 焓-孔隙率法模型

已有研究表明，在相變儲熱過程中，對流傳熱效果不可忽略。在確保結(jié)果準確性的基礎(chǔ)上，為了簡化數(shù)值模型，本文提出以下假設(shè)：1）考慮儲熱器保溫效果，儲熱器殼壁面設(shè)置為絕熱，并考慮管壁厚度；2）PCM熔化與凝固過程中考慮自然對流效應(yīng)；3）水工質(zhì)和相變材料的熱物理性質(zhì)被認為與溫度無關(guān)，均勻且各向同性。

對于相變過程的求解，應(yīng)用較為廣泛的是焓-孔隙法模型。其通過對式（4）～式（7）焓值的分布進行求解，從而確定相界面。控制方程包括連續(xù)性方程、動量方程和能量方程。

連續(xù)性方程如式（2）所示：

式中：[▽]——矢量微分算子；[V]——速度矢量。

動量方程如式（3）所示：

式中：[t]——時間，s；[ρ]——密度，kg/m3；p——壓強，Pa；[μ]——動力黏度，Pa·s；[β]——熱膨脹系數(shù)，1/K；[g]——重力加速度，m/s2；[T]——溫度，K；[Tref]——參考溫度，K；[S]——源項。

能量方程如式（4）所示：

式中：[cp]——比熱容，J/（kg·K）；[k]——導熱系數(shù)，W/（m·K）；[hsens]——顯熱焓，J/kg；[hlat]——潛熱焓，J/kg。

顯熱焓可表示為式（5）：

潛熱焓可表示為式（6）：

式中：[L]——完全熔化潛熱；[λ]——液相分數(shù)，[λ]可表示為式（7）：

式中：[Ts]——固體溫度，K；[Tliq]——液體溫度，K。動量方程中的源項[S]為達西定律的阻尼項，其表達式為式（8）：

式中：[Amush]——糊狀區(qū)常數(shù)，其值通常在104～107之間。本文取106計算求解。

1.2.2 邊界及初始條件

PCM的初始溫度[T0]為[25 ℃]；HTF進口溫度[Tf，i]設(shè)為70 ℃，入口流速為0.044 m/s；儲能單元殼壁面設(shè)置為絕熱。

2 結(jié)果與討論

2.1 立、臥式儲熱單元傳熱過程

對于立式和臥式儲熱單元，采用相同數(shù)量的結(jié)構(gòu)化網(wǎng)格及時間步長計算，其結(jié)果如圖4所示。圖4中立式儲熱單元內(nèi)PCM熔化時間比臥式更短，但并不意味立式儲熱單元的絕對優(yōu)越性。當采用管內(nèi)徑[di]為6 mm、外徑[do]為10 mm，而其他條件不變的儲熱單元做實驗對比，此時PCM熔化的過程如圖5所示。對比圖4、圖5可看到，當采用更小的管徑時（相當于增加相變材料區(qū)域厚度），臥式儲熱單元內(nèi)RT50熔化速度更快。這說明立式儲熱單元內(nèi)相變材料熔化過程受相變材料區(qū)域厚度影響較大。對于立式儲熱單元來講，由于自然對流發(fā)生在豎直方向，取豎直方向截面流動及液相云圖如圖6所示。由圖6可看到，整個相變?nèi)刍^程中流動區(qū)域、流動趨勢未發(fā)生明顯變化。但由于對流傳熱主要發(fā)生在固液兩相區(qū)及其附近，而導熱主要將熱量通過管外壁導向固液兩相區(qū)及固相區(qū)，熔化初始階段固液兩相區(qū)及固相區(qū)較大，導熱作用較強，而固液兩相區(qū)通過導熱快速產(chǎn)生后，在熔化過程中間階段兩相區(qū)域大小并未明顯增大。因此在熔化前期，導熱效果減弱，對流有所增強，但是熔化速度并未顯著放緩，而在60 min后，固相區(qū)消失，固液兩相區(qū)逐漸減小，但導熱與對流作用同時減弱，因而熔化速度迅速減小。而對于臥式儲熱單元，其PCM熔化速度經(jīng)歷了復(fù)雜的變化過程。結(jié)合圖7中心截面熔化及流動情況分析，發(fā)現(xiàn)熔化過程按傳熱情況可分為3個階段：初始時間段導熱主導區(qū)、中間時間段自然對流主導區(qū)及末尾時間段導熱主導區(qū)。

1）初始階段，直接接觸HTF的PCM薄層由于熱傳導作用快速升溫熔化，這個過程持續(xù)時間很短。由于PCM的導熱系數(shù)較小，當液相層到達一定厚度時，其熱傳導過程發(fā)生急劇惡化，液相變化趨于平緩。

2）中間階段，隨著傳熱過程的進行，液相區(qū)逐漸擴大。由于溫度差導致的熱膨脹效果，液相區(qū)及部分固液相混合區(qū)發(fā)生自然對流效應(yīng)。在密度差的驅(qū)動下，液相沿著溫度較高的管壁面一側(cè)爬升，沿著溫度較低的固相區(qū)域一側(cè)下降，形成環(huán)形對流，如圖7所示。在固液混合區(qū)邊界處的自然對流極大地增加了RT50的熔化速度，同時固液兩相區(qū)持續(xù)擴大。

3）末尾階段（60 min以后），隨著液相區(qū)的擴大，液相線的位置持續(xù)降低到下半?yún)^(qū)，固液混合區(qū)與液相之間邊界減小，同時液相區(qū)內(nèi)密度差減小，對流流動速度減小，且此時腔體上半部液相區(qū)內(nèi)的流動對傳熱失去強化效果，因此熔化速率急劇減小。最終，隨著固相區(qū)減薄，熱傳導效果增強，熔化速率略微增大直至完全熔化。

2.2 直肋、環(huán)肋傳熱強化效果

對于管內(nèi)徑[di]為10 mm、外徑[do]為14 mm的臥式儲熱單元，在管上加裝肋片可極大強化相變儲熱過程，提高傳熱效果與PCM熔化速率。從圖8可看到，用直肋管代替無肋管后PCM熔化速率提高了42.1%，而環(huán)肋管取代無肋管后PCM熔化速率提高了31.6%。顯然，直肋管的效果更好。圖9為直肋管儲熱單元中心截面處的熔化以及流動情況示意圖，可看到，加裝直肋之后其內(nèi)部PCM流動區(qū)域明顯擴大，基本覆蓋整個區(qū)域，對流傳熱得到顯著增強。因此PCM熔化速度大幅度提高，極大減少了熔化時間。圖10為直肋儲熱單元中心截面溫度分布云圖。

對于加裝環(huán)肋的情況，考慮到三維對稱性，在離中心位置0、3、6、9 mm處取4個截面進行觀察，結(jié)果如圖11所示。圖12為各個截面在傳熱進行到15 min時的液相分布及流動情況?？煽吹?，取肋片所在的A截面，無流動跡線。B截面液相流動區(qū)域最大。但隨著距離增大，液相流動區(qū)域有一定程度縮小，這與環(huán)肋肋片的附近導熱效果更好有關(guān)。取B截面觀察其不同時刻的流動情況，如圖13所示?？煽吹剑湓诓煌瑫r刻的流動狀況與無肋管的類似，環(huán)肋肋片的存在并沒有擴大其對流流動域，因此對于對流傳熱的強化十分有限，可得出，環(huán)肋對于傳熱過程傳熱強化主要在熱傳導方面，而直肋可強化熱對流及熱傳導過程。

2.3 儲熱器內(nèi)部自然對流形成機理

對于直肋管儲熱單元來說，其截面流動情況可簡化為圖14所示，開始傳熱后，熱的肋片和管壁與冷的殼壁形成溫度勢差，PCM熔化后沿熱壁面爬升，而遇到冷的壁面下降，在相對封閉的空間形成一個完整的流動循環(huán)，從而大幅度強化對流傳熱。而對于環(huán)肋管儲熱單元來說，其內(nèi)部PCM流動情況可簡化為圖15所示。由于溫度高的肋片彼此平行，因此

無法形成溫度勢差驅(qū)動。熱流沿著肋片爬升后，遇到頂部冷的殼壁，只能被迫改變方向，無法在當前截面內(nèi)完成循環(huán)流動，因此與無肋管內(nèi)PCM流動狀態(tài)類似。所以在管上加裝環(huán)肋并不能有效強化對流傳熱過程。

綜上所述，對于常見的導熱系數(shù)較小的相變材料，加裝肋片對熱傳導強化效果有限，應(yīng)當考慮其對自然對流傳熱過程的強化。設(shè)計原則應(yīng)當包括，在封閉空間內(nèi)，考慮熱壁面與冷壁面豎直相對，以便于形成流動循環(huán)。同時，當冷、熱壁面有一定高度時，對于受熱膨脹的流體來說，熱壁面與冷壁面的溫度梯度應(yīng)當豎直向上（如圖16所示），以便熱流可以沿著壁面爬升。

3 結(jié) 論

具有儲熱系統(tǒng)的太陽能熱發(fā)電技術(shù)是未來可再生能源中最具應(yīng)用前景的發(fā)電技術(shù)之一，因此對于儲熱系統(tǒng)的研究和改進具有重要意義。本文通過焓-孔隙法建立管殼式太陽能儲熱器數(shù)值模型，在考慮PCM密度隨溫度變化的基礎(chǔ)上，研究RT50在不同管殼式儲熱器中的熔化行為，并探究相變儲能傳熱過程的結(jié)構(gòu)優(yōu)化方法。得到的結(jié)論如下：

1） 對比研究立式和臥式換熱器儲熱單元內(nèi)PCM傳熱熔化過程，發(fā)現(xiàn)臥式儲熱單元內(nèi)PCM熔化過程可分為3個階段：初始時間段導熱主導，中間時間段自然對流主導及末尾時間段導熱主導。而立式儲熱單元內(nèi)PCM熔化速度對相變材料區(qū)域的厚度更加敏感。

2） 對于翅片強化過程研究表明，在臥式管殼換熱器內(nèi)管上加裝直肋后PCM完全熔化時間縮短了42.1%，而加裝環(huán)肋完全熔化時間縮短了31.6%，直肋管的效果比環(huán)肋管更好，這是由于直肋強化了PCM熔化過程中的自然對流傳熱。

3） 對于一些常見的相變儲熱材料，其導熱系數(shù)通常很小，肋片等結(jié)構(gòu)對傳熱導熱過程的強化效果有限，應(yīng)當考慮其自然對流傳熱強化傳熱過程。

［參考文獻］

［1］ 韓宗偉， 王一茹， 阿不來提·依米提， 等. 太陽能熱泵相變蓄熱供暖系統(tǒng)參數(shù)影響研究［J］. 太陽能學報， 2015，36（8）： 2028-2035.

HAN Z W， WANG Y R， YIMIT A， et al. Study on influence of parameters of solar energy heat pump system on" phase" change" heat" storage［J］." Acta" energiae" solaris sinica， 2015， 36（8）： 2028-2035.

［2］ AGYENIM F， HEWITT N， EAMES P， et al. A review of materials， heat transfer and phase change problem formulation for latent heat thermal energy storage systems （LHTESS）［J］." " Renewable" "and" "sustainable" "energy reviews， 2010， 14（2）： 615-628.

［3］ LACROIX M. Numerical simulation of a shell-and-tube latent heat thermal energy storage unit［J］. Solar energy， 1993， 50（4）： 357-367.

［4］ 劉洋， 段建國， 賀秀芬， 等. 低溫相變儲能單元強化傳熱特性實驗［J］. 太陽能學報， 2020， 41（6）： 335-340.

LIU Y， DUAN J G， HE X F， et al. Experiments on enhanced heat transfer performance of low temperature latent thermal energy storage unit［J］. Acta energiae solaris sinica， 2020， 41（6）： 335-340.

［5］ VISKANTA R. Phase-change heat transfer//Solar heat storage： latent heat materials［M］. CRC Press， 2018： 153-222.

［6］ 張仲彬， 孟雨欣， 劉永強. 相變膠囊堆積形態(tài)對蓄熱系統(tǒng)性能的影響分析［J］. 太陽能學報， 2021， 42（2）： 37-42.

ZHANG Z B，MENG Y X，LIU Y Q. Influence of phase change capsule packing configuration on performance of" "heat" storage" system［J］." Acta" energiae" solaris" sinica， 2021， 42（2）： 37-42.

［7］ 李椿， 王志華， 王建春， 等. 殼管式相變儲能換熱器性能研究與場協(xié)同效應(yīng)分析［J］. 太陽能學報， 2020， 41（3）： 226-233.

LI C， WANG Z H， WANG J C， et al. Performance study and field synergy analysis of shell and tube phase change energy" storage" heat" exchanger［J］. Acta" energiae" solaris sinica， 2020， 41（3）： 226-233.

［8］ HOSSEINI M J， RANJBAR A A， SEDIGHI K， et al. A combined experimental and computational study on the melting behavior of a medium temperature phase change storage material inside shell and tube heat exchanger［J］. International communications in heat and mass transfer， 2012， 39（9）： 1416-1424.

［9］ PAHAMLI Y， HOSSEINI M J， RANJBAR A A， et al. Analysis of the effect of eccentricity and operational parameters in PCM-filled single-pass shell and tube heat exchangers［J］. Renewable energy， 2016， 97： 344-357.

［10］ 田揚， 趙明， 胡明禹， 等. 加肋旋轉(zhuǎn)對相變蓄熱器蓄熱性能的影響及場協(xié)同分析［J］. 太陽能學報， 2021， 42（3）： 395-400.

TIAN Y， ZHAO M， HU M Y， et al. Effect of ribbed rotation on heat storage performance of phase change thermal storage unit and field synergy analysis［J］. Acta energiae solaris sinica， 2021， 42（3）： 395-400.

［11］ SCIACOVELLI A， GAGLIARDI F， VERDA V. Maximization of performance of a PCM latent heat storage system with innovative fins［J］. Applied energy， 2015， 137： 707-715.

NUMERICAL SIMULATION STUDY OF HEAT TRANSFER CHARACTERISTICS ON SOLAR TUBE-AND-SHELL PHASE

CHANGE HEAT STORAGE UNIT

Han Tao1，Ma Yanhua2，F(xiàn)ang Jiabin1，Yan Haobing1，Wei Jinjia1，Dong Siyuan2

（1. School of Chemical Engineering and Technology， Xi’an Jiaotong University， Xi’an 710049， China；

2. Faculty of Mechanical and Electrical Engineering， Yunnan Agricultural University， Kunming 650201， China）

Abstract：A numerical model is established using the enthalpy-porosity approach to study the heat transfer characteristics of a tube-and-shell phase-change heat exchanger filled with paraffin wax RT50. The influence of exchanger placement forms， tube diameters and fin structures on the phase change process of RT50 is analyzed. The results depicts that the vertical heat exchanger has a faster melting rate than the horizontal one as the tube diameter is large. However， the opposite results are obtained in case the tube diameter is small. For the horizontal exchanger， the heat conduction is dominant at the beginning and end of the melting process， while the natural convection plays a more important role at the intermediate stage of melting. Besides， the duration of the melting is mainly determined by the natural convection. In addition， the installation of fins on the outside of the tube can significantly increase the heat transfer coefficient and shorten the melting time. Compared with finless tube， the utilization of ring-fin tube in the horizontal heat storage device can shorten the melting time by 31.6%， and the utilization of straight-fin tube can shorten the melting time by 42.1%. Therefore， the straight fins have a better effect than the ring fins. The main reason is that straight fins can strengthen the natural convection during the melting process.

Keywords：solar energy； phase change material； melting； fins； natural convection； thermal energy storage