管狀陣列相變復合材料的傳熱特性

屈黎明，程志強，姚養(yǎng)庫，岳國富，楊文峰

(中國重型機械研究院股份公司，陜西 西安 710032)

0 引言

高效的熱存儲技術，特別是相變儲熱技術是合理有效利用現(xiàn)有能源、優(yōu)化使用可再生能源和提高能源利用效率的重要技術，也是近20年來世界節(jié)能領域一個非?；钴S的研究方向［1－6］。隨著工業(yè)余熱回收利用等眾多領域的迅猛發(fā)展，這為相變儲熱技術的進一步研究和發(fā)展提供了豐厚的條件。根據(jù)相變形式的不同，相變材料可以分為固－固相變、固－液相變、固－氣相變和液－氣相變4類。其中固－固相變、固－液相變是研究和實際中采用較多的相變類型。

目前，相變儲熱應用較多的相變蓄熱材料主要為有機物，但是有機相變儲熱材料一般導熱系數(shù)較低，致使在蓄熱系統(tǒng)的應用中傳熱性能差、蓄熱利用率低，降低了系統(tǒng)的效率。近年來，提出了復合相變儲能材料的概念［1－3］，它能有效克服單一的無機物或者有機物相變儲熱材料存在的導熱系數(shù)低的缺點，從而有效提高熱量存儲的速率，同時也改善了相變材料的應用效果及拓展其應用范圍。因而復合相變儲熱材料的結構優(yōu)化及相關傳熱機理、傳熱特性已成為相變儲能材料領域的研究熱點，其中由多孔金屬填充相變材料構成的相變復合材料受到廣泛關注［7－13］。但是多孔金屬本身的復雜結構帶來一定的理論和數(shù)值分析的困難，同時泡沫金屬的制作過程中的精度控制較難保證，且制作成本也較高。

陣列形式管狀結構內(nèi)填充相變材料構成的復合相變材料具有結構簡單、加工制作成本低的優(yōu)點，具有很好的工程應用價值。同時，作為其可以看做是多孔復合相變材料的一種簡化近似，從機理上也可以體現(xiàn)出多孔復合相變材料的基本傳熱特征。本文采用計算流體力學方法，對由不同管形的金屬管狀陣列構成的相變復合材料內(nèi)的非穩(wěn)態(tài)相變傳熱過程進行數(shù)值模擬，并分析管狀陣列構成的相變復合材料的傳熱特性。

1 計算模型

本文對兩種孔型的陣列結構，即圓孔型和方孔型所構成的相變復合材料進行研究?？仔椭兴畛涞南嘧儍δ懿牧线x取為石蠟，管狀陣列的基體材料取為銅，所研究的兩種管狀陣列復合結構的幾何形狀如圖1所示。

兩種不同孔型陣列的具體參數(shù)設置如表1所示?？椎慕Y構參數(shù)選取主要考慮孔隙率的不同及不同的管長對結果的影響。這里孔隙率定義為孔的體積占陣列總體積的份額。其中方形孔1和方形孔2用來研究不同孔隙率的影響，方孔型2和圓孔型1用來研究不同孔型結構的影響，圓孔型1和圓孔型2用來研究不同管長的影響。

圖1 本文所研究的兩種管狀陣列復合結構的外形示意圖

表1 管型結構的選擇

本文采用FLUENT軟件對金屬管狀陣列構成的相變復合材料內(nèi)的非穩(wěn)態(tài)相變傳熱過程進行數(shù)值模擬。該軟件提供了豐富的流動和傳熱求解模型，可對多種傳熱方式進行耦合求解。對本文的問題，由于相變只發(fā)生在陣列材料的孔隙中，在求解過程中金屬材料內(nèi)部的傳熱為無相變的熱傳導問題，而空隙中的石蠟會隨著溫度的升高而發(fā)生相變，需要考慮其中的相變換熱過程。

由于相變傳熱問題的溫度場在相界面附近一般存在大梯度，因而直接采用以溫度為控制變量的熱傳導方程存在數(shù)值穩(wěn)定性的問題。為了克服這一問題，相變換熱問題常采用焓法進行求解，F(xiàn)LUENT軟件在處理熔化/凝固過程時即采用了焓法。該方法引入一個液相率參量來表示液態(tài)物質在控制容積中所占的比例。不需要直接跟蹤相界面位置的變化，而是通過液相比例來間接的跟蹤相界面位置的變化，液體組分的計算基于焓的平衡來求解。

由于本文中不考慮流體速度，因此對于凝固/融化的問題，能量方程可以寫為

式中 ρ——材料的密度;

S——原項;

H——相變材料的焓。

可由顯焓h和潛熱ΔH計算為

其中，h=href+為參考焓，Tref為參考溫度，cp為定壓比熱。

焓法引入的液體分數(shù)β根據(jù)液化的程度進行

取值

其中Tsolidus為固化溫度，Tliquidus為液化溫度。借助液體分數(shù)β，(2)式中的潛熱量ΔH可以由查閱手冊得到的相變材料的相變潛熱L計算為

在研究周期性陣列結構時，同時兩個端面的溫度分別設定為500 K和288 K，初始溫度為常溫288 K。對于所研究的陣列型結構，考慮到傳熱過程的具有周期對稱性，為了減小計算量，將一個周期陣列的避面設為周期性邊界，這樣只需對一個周期陣列的模型進行三維模擬。

2 計算結果分析

2.1 整體三維溫度場分布特性

圖2(a)～(d)給出了計算得到的在t=2.78 h時刻時表1給出的4個管狀陣列相變復合材料內(nèi)的三維溫度場切片。從計算得到的溫度場等溫線可以看出，金屬銅內(nèi)的溫度場擴散速度均快于石蠟，這是由于金屬銅具有較大的導熱系數(shù)的原因。同時可以明顯注意到銅和石蠟內(nèi)部溫度分布的不均勻性，在x=－1.5 m和x=1.5 m位置金屬銅部分和石蠟部分等溫線移動速度比x=0 m位置要快，石蠟部分溫度和金屬銅部分溫度只差大于x=0 m位置石蠟部分和金屬銅部分溫差，說明兩側熱擴散速度比中間快。對比圖2(a)和圖2(b)，管長對三維溫度場的分布沒有明顯影響。

圖2 在t=2.78 h時刻分別在x=－1.5 m，x=0 m，x=1.5 m(從左至右)處的三維溫度場切片

由圖2(c)可以看出金屬銅部分溫度場擴散速度快于石蠟部分溫度場移動速度，并且位置為x=－1.5 m和x=1.5 m處溫度場與位置為x=0 m處溫度場相比，位置為x=0 m處切面上的等溫線分布更密集，石蠟部分和金屬銅部分溫差更大，說明兩側的熱傳導率高于中間部分。對比圖2(c)和圖2(d)可以看出，當孔隙率變大時軸向瞬態(tài)溫度場移動速度變慢了。這是由于隨著孔隙率增大，復合相變材料中多孔金屬骨架的材料金屬銅減少了，從而導致復合相變材料的導熱率降低，溫度場的移動速度變慢。

對比圖2(a)和2(c)可以得出相同孔隙率下不同元胞結構對復合相變材料傳熱性能的影響?？梢钥闯觯叫卧Y構中的軸向瞬態(tài)溫度場分布比圓孔型結構中的軸向瞬態(tài)溫度場分布要均勻，圓孔型結構中的軸向瞬態(tài)溫度場隨位置的變化較明顯，而位置變化對方形元胞結構中的軸向瞬態(tài)溫度場的影響較小。另外由圖2(c)和(d)可以看出，當孔隙率曾大時，金屬銅部分和石蠟部分溫差變大，等溫線趨于更不均勻。

2.2 石蠟和銅中的平均溫度分布特性

由前面的分析看出，銅和石蠟中的熱擴散存在一定速度差異。為更好的分析該復合材料的傳熱特性，下面對石蠟和銅內(nèi)部的平均溫度場進行對比分析。

圖3 在時刻t=2.78 h石蠟和金屬銅的平均溫度隨位置變化曲線

圖3給出了在時刻t=2.78 h石蠟和金屬銅內(nèi)部的平均溫度隨位置變化。從圖中可以看出，熱流方向為由右至左，即高溫端在右側。按照熱流方向，其中金屬銅部分的溫度隨位置變化曲線在靠右端邊界處存在斜率不變的區(qū)域，接著斜率開始減小，最終趨于定常。整個變化趨勢存在很明顯的非線性，這說明參與相變導熱的金屬銅，其內(nèi)部的導熱規(guī)律和僅為單一的銅材料的導熱規(guī)律有較大不同。

對于石蠟，其內(nèi)部平均溫度隨位置變化曲線在初期斜率與金屬銅相同，且它們的曲線接近重合。另外石蠟的溫度下降快于金屬銅部分溫度曲線，直至在相變點出其曲線出現(xiàn)轉折，重新下降后其曲線最終與金屬銅溫度曲線重合。由圖可以看出石蠟溫度曲線在相變前后具有很明顯的變化，這是由于液態(tài)石蠟的導熱率要高于固態(tài)石蠟，同時石蠟融化時吸收相變潛熱，導致經(jīng)過相變點后熱流量減少，石蠟的溫度曲線變平緩。對比圖3(c)和(d)可以看出當空隙率增大時，相變材料的相變點移動速度變慢了，相同孔隙率下，且圓孔型結構中相變材料在相變點處的溫差比方形管的小。

3 結論

本文對填充石蠟的圓形和方形的兩種不同管形的金屬銅管狀陣列所構成的相變復合材料內(nèi)的非穩(wěn)態(tài)相變傳熱過程進行了數(shù)值模擬，分析了管狀陣列構成的相變復合材料的傳熱特性。主要得到以下結論:

(1)對于該管狀陣列結構，金屬銅內(nèi)的溫度場擴散速度優(yōu)于石蠟，且銅和石蠟內(nèi)部溫度分布具有明顯的不均勻性。相同孔隙率條件下，方形管狀結構中的瞬態(tài)溫度場分布的均勻性優(yōu)于圓孔型結構中的軸向瞬態(tài)溫度場分布。

(2)當孔隙率增大時，相變材料中熱擴散速度明顯降低，這是由于孔隙率曾大，導致金屬銅材料減少，復合相變材料的有效導熱系數(shù)減小所致。

(3)數(shù)值模擬是研究相變復合材料內(nèi)傳熱過程的有效手段，可以應用于工程中對相變材料的傳熱特性進行預測和設計。

［1］陳愛英，汪學英，曹學增.相變儲能材料的研究進展與應用［J］.材料導報，2003，17(5):42－72.

［2］呂學文，考宏濤.復合相變儲能材料的傳熱性能研究進展［J］.節(jié)能，2009(28):9－11.

［3］方貴銀，李輝.復合相變蓄熱材料研制及性能分析［J］.現(xiàn)代化工，2003，23(12):30－33.

［4］沈學忠，張仁元.相變儲能材料的研究和應用［J］.節(jié)能技術，2006，24(5):460－463.

［5］周劍敏，張東，吳科如.相變儲能復合材料實驗研究與分析［J］.節(jié)能技術，2003，21(6):5－7.

［6］王俊，曲延濤，孟廣軍，等.相變儲能材料應用于基站設備初探［J］.節(jié)能技術，2008，26(2):151－154.

［7］彭東華，陳振乾，施明恒.泡沫金屬內(nèi)相變材料融化傳熱過程的數(shù)值模擬分析［J］.工程熱物理學報，2009，30 (6):1025－1028.

［8］張濤，余建祖.泡沫銅作為填充材料的相變儲熱實驗［J］.北京航空航天大學學報，2007，33(9):1021－1024.

［9］程文龍，韋文靜.高孔隙率泡沫金屬相變材料儲能、傳熱特性［J］.太陽能學報，2007，28(7):739－744.

［10］崔海亭，曹向茹，彭培英.泡沫金屬相變材料熔化傳熱過程的數(shù)值分析［J］.功能材料，2010，41(8):1439－1442.

［11］杲東彥，陳振乾.格子Boltzmann方法模擬泡沫金屬內(nèi)相變材料熱傳導融化傳熱過程［J］.熱科學與技術，2011，10(1):6－11.

［12］崔海亭，劉鳳青，朱金達，等.高孔隙率泡沫金屬對相變蓄熱的強化研究［J］.河北科技大學學報，2010，31(2): 93－96.

［13］遲蓬濤，謝永奇，余建祖，等.一種新型儲能裝置充冷過程的實驗及分析［J］.北京航空航天大學學報，2011，37 (9):1070－1075.