冷庫(kù)冷凝熱回收用三套管相變蓄熱器蓄熱特性模擬

羅小龍 袁艷平 曹曉玲 孫亮亮 向 波

?

冷庫(kù)冷凝熱回收用三套管相變蓄熱器蓄熱特性模擬

羅小龍 袁艷平 曹曉玲 孫亮亮 向 波

（西南交通大學(xué)機(jī)械工程學(xué)院 成都 610031）

為回收冷庫(kù)制冷系統(tǒng)產(chǎn)生的冷凝熱用于制取生產(chǎn)及生活熱水，在壓縮機(jī)和冷凝器之間接入三套管相變蓄熱器。通過(guò)焓法建立三套管相變蓄能換熱器的數(shù)學(xué)模型，分析制冷劑入口溫度、制冷劑流量及結(jié)構(gòu)參數(shù)對(duì)其蓄熱特性的影響。結(jié)果表明：入口溫度從110℃升高到130℃，熔化時(shí)間降低20.5%，且制冷劑入口溫度越高，總蓄熱量越大；制冷劑流量從0.008kg/s增加到0.024kg/s，熔化時(shí)間減小55.9%，改變流量只影響蓄熱快慢而不影響總蓄熱量大小，且隨著流量增大其對(duì)蓄熱性能的影響減弱；管徑對(duì)蓄熱性能的影響較小，較小的管徑比有利于減小流動(dòng)阻力而不影響蓄熱性能。

熱回收；三套管相變蓄熱器；蓄熱特性；焓法；模擬

0 引言

冷庫(kù)作為大型制冷場(chǎng)所，其制冷系統(tǒng)運(yùn)行時(shí)會(huì)向環(huán)境排放大量冷凝熱，若可以將排放的熱量回收用于制取生產(chǎn)及生活熱水，不僅節(jié)約能源，而且能減少對(duì)環(huán)境的熱污染。傳統(tǒng)的熱回收技術(shù)大多為顯熱回收，即直接加熱水或其他液體來(lái)蓄熱，這種方法的缺陷在于其蓄熱密度低、蓄熱裝置體積大、占用空間且成本較高。而潛熱蓄熱（相變蓄熱），利用相變材料（PCM）發(fā)生相變時(shí)會(huì)吸收或放出大量的潛熱來(lái)進(jìn)行蓄熱，其蓄熱密度高、放熱過(guò)程溫度波動(dòng)小[1]，能有效解決能量供需在時(shí)間、空間及強(qiáng)度上的不匹配現(xiàn)象[2]。

為有效回收制冷系統(tǒng)冷凝熱，本文提出在冷庫(kù)制冷系統(tǒng)壓縮機(jī)和冷凝器之間接入三套管相變蓄熱器。三套管相變蓄熱器是一種新型的集成節(jié)能設(shè)備，能夠?qū)⑵胀ㄌ坠軗Q熱器的換熱功能與相變材料的蓄熱功能有機(jī)地結(jié)合[3]，其在空調(diào)系統(tǒng)中的應(yīng)用已有大量文獻(xiàn)報(bào)道。文獻(xiàn)[4]在太陽(yáng)能熱泵系統(tǒng)供熱工況下，針對(duì)三套管相變蓄熱器穩(wěn)態(tài)傳熱模式建立數(shù)學(xué)模型，研究了相變材料層的溫度場(chǎng)變化，三套管外層為經(jīng)太陽(yáng)能集熱器換熱的低溫?zé)崴艧幔瑑?nèi)層為制冷劑吸熱。牛福新等[5,6]將三套管蓄能換熱器應(yīng)用于太陽(yáng)能-空氣源熱泵集成系統(tǒng)作為蒸發(fā)器，分別模擬研究了三套管蓄能換熱器的蓄熱和蓄冷動(dòng)態(tài)特性。Long等[7]將PCM應(yīng)用于空氣源熱泵熱水器，以三套管為模型模擬了相變單元蓄、放熱過(guò)程，研究了傳熱速率、相界面的移動(dòng)以及PCM的溫度場(chǎng)變化，并進(jìn)行了實(shí)驗(yàn)驗(yàn)證。文獻(xiàn)[8]基于空調(diào)或太陽(yáng)能領(lǐng)域的余熱回收研究了三套管相變蓄熱器的釋熱和蓄熱過(guò)程，得出了傳熱流體的入口溫度和流速對(duì)換熱特性的影響規(guī)律。上述針對(duì)三套管相變蓄熱器傳熱特性的研究其熱流體在外側(cè)管內(nèi)流動(dòng)、冷流體在內(nèi)側(cè)管內(nèi)流動(dòng)，而熱流體在內(nèi)側(cè)流動(dòng)、冷流體在外側(cè)流動(dòng)的模式缺乏相關(guān)研究。

本文在前人的研究基礎(chǔ)上，針對(duì)熱流體在內(nèi)側(cè)管內(nèi)流動(dòng)、冷流體在外側(cè)管內(nèi)流動(dòng)的三套管相變蓄熱器進(jìn)行研究。建立了基于焓法的三套管相變蓄熱單元數(shù)學(xué)模型，分析了制冷劑入口溫度、流量和結(jié)構(gòu)參數(shù)對(duì)其蓄熱特性的影響，為三套管相變蓄熱器在冷凝熱回收系統(tǒng)中的應(yīng)用打下基礎(chǔ)。

1 物理模型

本文以一直接蒸發(fā)式氨制冷系統(tǒng)為背景設(shè)計(jì)蓄熱器結(jié)構(gòu)，其高壓級(jí)壓縮機(jī)排氣量273.69m3/h，冷凝溫度38℃，PCM初始溫度為20℃。三套管相變蓄熱器由4組相變蓄能換熱單元組成，蓄能換熱單元的結(jié)構(gòu)如圖1所示，內(nèi)管為制冷劑管，中間管為相變蓄能管，外套管為冷水管道，各管管徑及管長(zhǎng)尺寸見(jiàn)表1，制冷劑與水流向相反。所用蓄熱材料為棕櫚蠟—膨脹石墨復(fù)合材料[9]，其相變溫度為78～82℃，相變潛熱為150.9kJ/kg，導(dǎo)熱系數(shù)為5.92W/(m·K)，密度為1000kg/m3。

圖1 三套管相變蓄能換熱單元示意圖

表1 三套管蓄能換熱器結(jié)構(gòu)參數(shù)

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

在建模之前進(jìn)行如下假設(shè)：（1）制冷劑在管內(nèi)作一維流動(dòng)；（2）兩相流在同一流動(dòng)截面上處于熱力平衡狀態(tài)；（3）制冷劑側(cè)能量方程中忽略動(dòng)能和勢(shì)能的影響；（4）忽略相變材料在液態(tài)時(shí)的自然對(duì)流；（5）忽略金屬管壁的熱阻和熱容；（6）忽略重力的影響。根據(jù)以上假設(shè)，傳熱過(guò)程中圓柱體各個(gè)半徑方向上溫度場(chǎng)相同，相變材料溫度場(chǎng)分布可視為是二維的，即沿徑向和軸向變化，取環(huán)形微元體，如圖2所示。

圖2 簡(jiǎn)化后的模型及微元體示意圖

2.1 制冷劑能量方程

制冷劑在管內(nèi)流動(dòng)過(guò)程中會(huì)發(fā)生相態(tài)的變化，在兩相區(qū)采用均相流模型，假設(shè)氣液兩相流速相等，其能量控制方程如下[10]：

m(h-h)=q（1）

q=πdΔxa(T-T) （2）

式中，r為制冷劑質(zhì)量流量，kg/s；j為制冷劑沿流動(dòng)方向第j點(diǎn)的焓，J/kg；r為熱流密度，W；0為內(nèi)層套管直徑，m；?為微元長(zhǎng)度，m；r為制冷劑側(cè)對(duì)流換熱系數(shù)，W/(m2·K)；r為制冷劑溫度，K；p為管壁溫度，K。

對(duì)于單相流（包括過(guò)熱氣態(tài)制冷劑、過(guò)冷液態(tài)制冷劑）管內(nèi)湍流換熱，r采用Dittus-Boelter關(guān)聯(lián)式表達(dá)[11]：

式中，為雷諾數(shù)；為普朗特?cái)?shù)；為制冷劑導(dǎo)熱系數(shù)，W/m·K；i為管內(nèi)徑，m。

對(duì)于兩相凝結(jié)換熱，換熱系數(shù)采用Shah關(guān)聯(lián)式[11]：

式中，L為蒸氣凝結(jié)成液相時(shí)的換熱系數(shù)，采用式（3）計(jì)算；為制冷劑的干度；為工作壓力，kPa；crit為臨界壓力，kPa。

2.2 相變材料能量方程

相變導(dǎo)熱問(wèn)題的求解采用焓法模型，焓法的主要思路是采用熱焓和溫度一起作為待求變量，在整個(gè)區(qū)域建立一個(gè)統(tǒng)一的能量方程，利用數(shù)值方法求出熱焓分布，然后確定兩相界面[12]。PCM微元體能量方程[10]如下：

能量守恒分析得出：

式中，為PCM的比焓，J/kg，這里取0℃時(shí)焓的基準(zhǔn)值為零，即PCM為固態(tài)時(shí)，=p；為PCM的導(dǎo)熱系數(shù)，W/m·K；為PCM節(jié)點(diǎn)溫度，K；i-1，i，j-1，j分別為微元體左、右、上、下表面的面積，m2。

求解過(guò)程需將溫度轉(zhuǎn)換為焓，假設(shè)相變過(guò)程發(fā)生在溫度區(qū)間[m-，m+]內(nèi)，溫度與焓的轉(zhuǎn)換關(guān)系如下：

式中，s為相變材料的固相比熱，J/kg·℃；l為相變材料的液相比熱，J/kg·℃；m為相變材料的相變潛熱，J/kg；m為相變材料的相變溫度，℃。

蓄熱過(guò)程內(nèi)表面與制冷劑對(duì)流換熱：

其他表面視為絕熱：

2.3 數(shù)值求解方法

將求解區(qū)域劃分為均勻網(wǎng)格，再將控制微分方程在以網(wǎng)格節(jié)點(diǎn)為中心的控制體積內(nèi)用有限差分法進(jìn)行數(shù)值離散。PCM能量方程的離散采用隱式格式，用Gauss-Seidel迭代法求解方程。對(duì)制冷劑方程，可根據(jù)初始條件直接求解代數(shù)方程，得到制冷劑沿流動(dòng)方向上的溫度場(chǎng)。編寫(xiě)MATLAB數(shù)值計(jì)算程序，制冷劑在不同狀態(tài)下的物性參數(shù)通過(guò)REFPROP軟件調(diào)取。在每個(gè)時(shí)間步長(zhǎng)內(nèi)，分別求解制冷劑和PCM能量方程，采用預(yù)估-校正法耦合計(jì)算，迭代直至收斂，收斂標(biāo)準(zhǔn)設(shè)置為10-6。

2.4 模型驗(yàn)證

為驗(yàn)證上述模型和計(jì)算方法的正確性，以已有文獻(xiàn)的實(shí)驗(yàn)結(jié)果對(duì)比驗(yàn)證。首先，利用Atal等[13]針對(duì)雙套管相變蓄熱器的實(shí)驗(yàn)研究，驗(yàn)證單相流體與相變材料換熱過(guò)程，對(duì)比蓄、放熱過(guò)程中相變材料平均溫度隨時(shí)間的變化，如圖3所示，模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)接近，誤差小于10%。制冷劑的換熱涉及相態(tài)變化，其原理較為復(fù)雜，為驗(yàn)證制冷劑能量方程，對(duì)比文獻(xiàn)[14]中冷凝器在環(huán)境溫度26.4℃和17.6℃兩種工況下的實(shí)驗(yàn)結(jié)果，如圖4所示，結(jié)果基本吻合，最大誤差為12.3%。綜上，本文所用三套管相變蓄熱器的數(shù)學(xué)模型及求解方法合理，可用于研究相變材料與制冷劑及水的換熱過(guò)程。

圖3 相變材料平均溫度對(duì)比

圖4 制冷劑冷凝過(guò)程溫度對(duì)比

3 結(jié)果分析

為研究相變蓄熱單元的在不同運(yùn)行條件下的蓄熱性能，通過(guò)數(shù)值模擬不同制冷劑入口溫度、質(zhì)量流量和管徑的影響。標(biāo)準(zhǔn)工況為制冷劑入口溫度r,in=120 ℃，質(zhì)量流量r=0.012 kg/s，管徑比1/0=2.4。在分析某一參數(shù)對(duì)蓄熱特性的影響時(shí)，其他參數(shù)固定。

3.1 制冷劑入口溫度影響

圖5顯示了制冷劑入口溫度對(duì)液相分?jǐn)?shù)的影響。由圖可知，整個(gè)蓄熱過(guò)程可分為三個(gè)階段：（1）固相顯熱蓄熱階段，PCM溫度低于相變溫度，液相分?jǐn)?shù)為0；（2）潛熱蓄熱階段，PCM溫度達(dá)到相變溫度并保持穩(wěn)定，液相分?jǐn)?shù)從0增加到1；（3）液相顯熱蓄熱階段，PCM全部熔化為液態(tài)，其溫度高于相變溫度，液相分?jǐn)?shù)為1。在第一階段，隨著制冷劑入口溫度升高，制冷劑與PCM之間的溫差增大，導(dǎo)致傳熱速率增加。傳入PCM的熱量越多，PCM溫度越快升高到相變溫度。所以隨著入口溫度升高，顯熱蓄熱階段縮短。PCM溫度達(dá)到相變溫度后進(jìn)入潛熱蓄熱階段，在這一階段，隨著入口溫度升高，液相分?jǐn)?shù)曲線斜率增加即熔化速率更大，直到液相分?jǐn)?shù)達(dá)到1，潛熱蓄熱過(guò)程結(jié)束。隨著入口溫度升高，固相顯熱蓄熱階段縮短，PCM熔化速率增加，所以蓄熱時(shí)間明顯減少。

圖5 制冷劑入口溫度對(duì)液相分?jǐn)?shù)的影響

圖6 制冷劑入口溫度對(duì)固-液相界面的影響

圖6顯示了制冷劑入口溫度對(duì)=20min時(shí)固-液相界面分布的影響，縱坐標(biāo)表示固-液相界面到內(nèi)層套管表面的距離。固-液相界面可反映不同位置的熔化速率，對(duì)某一特定的入口溫度，越靠近內(nèi)層管壁，PCM熔化越快；不同入口溫度下，隨著入口溫度升高，蓄熱速率增大，導(dǎo)致PCM熔化速率增加。所以隨著入口溫度升高，PCM在相同時(shí)間內(nèi)的熔化分?jǐn)?shù)增加，固-液相界面向制冷劑出口側(cè)移動(dòng)，且其變化是均勻的。

圖7 制冷劑入口溫度對(duì)蓄熱量的影響

不同制冷劑入口溫度下PCM蓄熱量隨時(shí)間變化情況如圖7所示。在固相顯熱蓄熱階段，由于傳熱溫差較大，蓄熱量迅速上升。然而在這一階段，能量以顯熱形式存儲(chǔ)，導(dǎo)致PCM溫度快速上升，溫差減小，所以傳熱速率隨時(shí)間逐漸降低，蓄熱量曲線斜率減小。不同的入口溫度下，其蓄熱量差別不大。在潛熱蓄熱階段，PCM溫度處于相變溫度區(qū)間，對(duì)于給定的入口溫度，傳熱溫差幾乎為定值，所以蓄熱量直線上升；且入口溫度越高，蓄熱量越大。在液相顯熱蓄熱階段，PCM溫度快速升高，溫差隨時(shí)間減小，所以傳熱速率隨時(shí)間減小，直到減小到0，蓄熱量達(dá)到最大值；且入口溫度越高，蓄熱過(guò)程越快達(dá)到穩(wěn)定，最終的蓄熱量也越大。

圖8 制冷劑入口溫度對(duì)蓄熱時(shí)間的影響

圖8給出了PCM蓄熱量達(dá)到穩(wěn)定時(shí)所用時(shí)間與制冷劑入口溫度的關(guān)系，由圖可知在制冷劑流量一定的情況下，蓄熱時(shí)間與入口溫度在一定范圍內(nèi)近似呈線性關(guān)系，這是因?yàn)槿肟跍囟仍礁?，制冷劑與PCM之間的平均溫差越大，所以完成蓄熱所用時(shí)間就越短。當(dāng)入口溫度從110℃升高到130℃，熔化時(shí)間從78min降到62min，降低了約20.5%。

3.2 制冷劑流量影響

圖9 制冷劑流量對(duì)液相分?jǐn)?shù)的影響

圖9顯示了制冷劑流量對(duì)液相分?jǐn)?shù)的影響。在固相顯熱蓄熱階段，隨著制冷劑流量增加，傳熱性能增強(qiáng)，導(dǎo)致PCM溫度快速上升到相變溫度，所以隨著流量增加顯熱蓄熱階段縮短。在潛熱蓄熱階段，隨著流量增加，傳熱性能增強(qiáng)，導(dǎo)致PCM熔化速率增大，液相分?jǐn)?shù)曲線斜率隨流量增加而增大，即潛熱蓄熱階段縮短。所以蓄熱時(shí)間隨流量增加快速減小。

圖10 制冷劑流量對(duì)固-液相界面的影響

圖10顯示了制冷劑流量對(duì)=20 min時(shí)固-液相界面的影響。隨著制冷劑流量增加，傳熱速率增大，PCM熔化速率加快，所以相同時(shí)間內(nèi)PCM的熔化分?jǐn)?shù)增加，固-液相界面向制冷劑出口側(cè)移動(dòng)。此外，越靠近制冷劑管壁面，流量對(duì)固-液相界面的影響越顯著，意味著增加制冷劑流量導(dǎo)致固-液相界面分布和PCM熔化速率的非均勻性增加。

不同流量下PCM蓄熱量隨時(shí)間變化情況如圖11所示。在固相顯熱蓄熱階段和潛熱蓄熱階段，隨著制冷劑流量增大，制冷劑與PCM之間的對(duì)流換熱系數(shù)增加，導(dǎo)致傳熱量增大。在液相顯熱蓄熱階段，傳熱速率隨時(shí)間降低，最終趨于0，蓄熱過(guò)程結(jié)束，且制冷劑流量越大蓄熱速率減小得越快，意味著PCM完成蓄熱所用時(shí)間越短。各流量下PCM蓄熱量最終趨于同一穩(wěn)定值，這是因?yàn)樵谥评鋭┤肟跍囟纫欢ǖ臈l件下，PCM溫度場(chǎng)最終趨于相同，故所蓄的熱量值相等。

圖11 制冷劑流量對(duì)蓄熱量的影響

圖12 制冷劑流量對(duì)蓄熱時(shí)間的影響

圖12給出了PCM蓄熱達(dá)到穩(wěn)定所需時(shí)間與制冷劑流量的關(guān)系。隨著制冷劑流量增加，對(duì)流換熱系數(shù)增大，導(dǎo)致制冷劑和PCM之間的傳熱熱阻減小，所以隨著制冷劑流量增加，PCM熔化速率增大，熔化時(shí)間減小。隨著流量從0.008kg/s增加到0.024kg/s，熔化時(shí)間從93min減小到41min，減小約55.9%，所以增加制冷劑流量可顯著減小熔化時(shí)間，但隨著流量繼續(xù)增加，其影響越來(lái)越弱。這是由于制冷劑流量增大時(shí)，雖然其與PCM的對(duì)流換熱熱阻減小，但PCM的導(dǎo)熱熱阻不變，所以對(duì)流換熱熱阻占總熱阻的比例逐漸下降，因而對(duì)蓄熱時(shí)間的影響也會(huì)減小。

3.3 管徑影響

在相同的制冷劑入口溫度和質(zhì)量流量下研究管徑對(duì)蓄熱特性的影響。為了保持PCM體積恒定，應(yīng)保持橫截面PCM面積為定值。不同結(jié)構(gòu)下管徑尺寸見(jiàn)表2。

表2 不同結(jié)構(gòu)下套管單元尺寸

圖13顯示了管徑比對(duì)液相分?jǐn)?shù)的影響。隨著管徑比增大，相同時(shí)間內(nèi)熔化分?jǐn)?shù)增加，但相比制冷劑入口溫度和流量，管徑的影響很弱。因?yàn)楣軓奖仍龃笠馕吨?減小，制冷劑流速增大，但傳熱面積相應(yīng)減小?？梢?jiàn)，隨著管徑比增大，流動(dòng)阻力會(huì)大幅增加，但蓄熱性能并不會(huì)得到很大提升。

圖13 管徑對(duì)液相分?jǐn)?shù)的影響

圖14 管徑對(duì)固-液相界面的影響

管徑比對(duì)固-液相界面的影響如圖14所示，縱坐標(biāo)表示固-液界面距制冷劑管的距離。隨著管徑比增大，制冷劑管管徑減小，制冷劑流速增大，導(dǎo)致PCM熔化速率加快，所以固-液界面向出口側(cè)移動(dòng)，隨著管徑比進(jìn)一步增大，影響逐漸減弱。此外，越靠近制冷劑管，管徑比對(duì)-液相界面分布的非均勻性增加。

圖15顯示了管徑比對(duì)蓄熱量的影響。固相顯熱蓄熱階段，不同管徑下PCM蓄熱量基本相同。潛熱蓄熱階段，隨著管徑比增大，蓄熱量增大，但增幅較小。液相顯熱蓄熱階段，不同結(jié)構(gòu)下蓄熱量逐漸趨于相等，可見(jiàn)管徑比對(duì)三套管相變蓄熱器蓄熱量的影響很小。

圖15 管徑對(duì)蓄熱量的影響

4 結(jié)論

三套管相變蓄能換熱器作為冷庫(kù)冷凝熱回收系統(tǒng)的核心部件，其蓄熱特性直接決定系統(tǒng)的性能，為此本文基于焓法模型建立了三套管相變蓄能換熱單元的數(shù)學(xué)模型，并通過(guò)數(shù)值模擬，研究了不同運(yùn)行條件和結(jié)構(gòu)參數(shù)下三套管相變蓄熱單元的蓄熱特性，結(jié)果表明：

（1）制冷劑入口溫度對(duì)蓄熱性能影響較大，入口溫度越高，PCM完成相變用時(shí)越短，在相同時(shí)間內(nèi)蓄得的熱量越多，且入口溫度對(duì)固-液相界面的影響是均勻的。

（2）制冷劑流量對(duì)蓄熱性能影響較大，流量越大，蓄熱速率越大，完成相變用時(shí)越短，且隨著流量的增大其對(duì)蓄熱性能的影響減弱。此外，流量增大，固-液相界面分布的非均勻性增加。

（3）管徑對(duì)蓄熱性能的影響較小，隨著管徑比增大，蓄熱速率加快，但蓄熱總量并不受影響。較小的管徑比有利于在蓄熱性能幾乎不受影響的情況下減小流動(dòng)阻力。

[1] 李偉,李新.翅片管簇式相變蓄熱器的實(shí)驗(yàn)研究[J].煤氣與熱力,2011,31(9):9-11.

[2] 袁艷平,向波,曹曉玲,等.建筑相變儲(chǔ)能技術(shù)研究現(xiàn)狀與發(fā)展[J].西南交通大學(xué)學(xué)報(bào),2016,51(3):585-598.

[3] 楊靈艷,倪龍,姚楊,等.蓄能型太陽(yáng)能與空氣源熱泵集成系統(tǒng)運(yùn)行模式[J].煤氣與熱力,2009,29(1):1-4.

[4] 楊靈艷,倪龍,姚楊,等.三套管相變蓄能換熱器供熱工況穩(wěn)態(tài)模擬分析[J].煤氣與熱力,2010,30(11):4-7.

[5] 牛福新,倪龍,姚楊,等.三套管蓄能換熱器的蓄熱特性研究[J].湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,38(7): 69-72.

[6] 牛福新,倪龍,姚楊,等.直接蒸發(fā)式三套管蓄能換熱器的蓄冷特性[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào),2011,43(6): 56-60.

[7] Long Jianyou, Zhu Dongsheng. Numerical and exper- imental study on heat pump water heater with PCM for thermal storage [J]. Energy and Buildings, 2008,40(4): 666-672.

[8] Long Jianyou. Numerical and experimental investig- ation for heat transfer in triplex concentric tube with phase change material for thermal energy storage[J]. Solar Energy, 2008,82(11):977-985.

[9] Xia Mingzhu, Yuan Yanping, Zhao Xudong, et al. Cold storage condensation heat recovery system with a novel composite phase change material[J]. Applied Energy, 2016,175:259-268.

[10] 黃挺.空氣源熱泵除霜用相變蓄能換熱器的模擬研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2007.

[11] 張春路.制冷空調(diào)系統(tǒng)仿真原理與技術(shù)[M].北京:化學(xué)工業(yè)出版社,2012.

[12] 袁修干.相變蓄熱技術(shù)的數(shù)值仿真及應(yīng)用[M].北京:國(guó)防工業(yè)出版社,2013.

[13] Atal A, Wang Y, Harsha M, etc. Effect of porosity of conducting matrix on a phase change energy storage device[J]. International Journal of Heat and Mass Transfer,2016,93:9-16.

Heat Storage Characteristics of Triplex Tube Heat Exchanger Used for Cold Storage Condensation Heat Recovery

Luo Xiaolong Yuan Yanping Cao Xiaoling Sun Liangliang Xiang Bo

( School of Mechanical Engineering, Southwest Jiaotong University, Chengdu, 610031 )

In order to recover the condensation heat of cold storage refrigeration system and use for domestic hot water preparation or industrial hot water supply, a triplex tube heat exchanger (TTHX) based on phase change materials was set between the compressor and the condenser. By establishing a mathematical model of TTHX according to enthalpy method, the dynamic characteristics was simulated, and the impact of refrigerant inlet temperature, mass flow rate and tube geometric parameters on the storage characteristics was investigated. The results show that with refrigerant inlet temperature increasing from 110℃ to 130℃, the melting time reduces 20.5%. The higher the refrigerant inlet temperature is, the greater the total heat storage capacity. When mass flow rate of refrigerant increases from 0.008kg/s to 0.024kg/s, the melting time reduces 55.9%. Changing the mass flow rate only affects heat storage rate without affecting the total heat storage capacity, and with the increase of mass flow rate, its effect on heat storage performance is weakened. Tube diameter has little effect on performance, and smaller diameter ratio helps to reduce flow resistance without affecting heat storage performance.

heat recovery; triplex tube heat exchanger; heat storage characteristics; enthalpy method; simulation

1671-6612（2018）06-579-07

TK02

A

建筑環(huán)境與能源高效利用四川省青年科技創(chuàng)新研究團(tuán)隊(duì)項(xiàng)目（2015TD0015）

羅小龍（1992-），男，在讀碩士研究生，E-mail：xiaolongluo@my.swjtu.edu.cn

袁艷平（1973-），男，博士，教授，E-mail：ypyuan@home.swjtu.edu.cn

2018-04-26