蒲 亮 周 洋 齊 迪

(西安交通大學能源與動力工程學院制冷與低溫工程系 西安 710049)

星型翅片管換熱器結(jié)霜過程的數(shù)值模擬

蒲 亮 周 洋 齊 迪

(西安交通大學能源與動力工程學院制冷與低溫工程系 西安 710049)

采用數(shù)值模擬的方法對星型翅片管換熱器的結(jié)霜過程進行了研究。得到了霜層厚度和霜層密度等物性參數(shù)的變化規(guī)律，分析了其空氣側(cè)結(jié)霜過程中的傳熱傳質(zhì)特性，研究了翅片高度和翅片個數(shù)兩個參數(shù)對結(jié)霜過程的影響。結(jié)果表明，結(jié)霜現(xiàn)象會造成換熱器換熱性能的惡化,并且翅片個數(shù)越多、翅片高度越小，傳熱惡化更嚴重。

星型翅片管換熱器 數(shù)值模擬 結(jié)霜過程 傳質(zhì)傳熱特性

1 引 言

氣化器是低溫換熱器的一種，其主要應(yīng)用在液氮、液氧、LNG等各種低溫介質(zhì)的再氣化過程中?？諟厥綒饣魇堑蜏負Q熱器的一種常見形式，其主要結(jié)構(gòu)形式為星型翅片管。在氣化過程中，其從空氣中吸收熱量并將其傳遞給管內(nèi)的低溫介質(zhì)使相變而氣化。與其他氣化器相比，其結(jié)構(gòu)簡單、經(jīng)濟性好、符合節(jié)能環(huán)保的要求，因此在多個領(lǐng)域中得到了廣泛的應(yīng)用[1]。但星型翅片管氣化器管內(nèi)流動的是各種低溫液體介質(zhì)，空氣中的水蒸氣很容易在管壁外側(cè)凝華而結(jié)霜，這會造成星型翅片管換熱器傳熱性能的嚴重惡化。

目前，國內(nèi)外不少學者對換熱器結(jié)霜進行了大量的研究[2-4]。但在這些研究中，結(jié)霜模型多是在-18 —-25 ℃的普冷工況提出的，而低溫星型翅片管的管內(nèi)介質(zhì)的溫度多在-160 ℃以下，低溫下的結(jié)霜現(xiàn)象還會出現(xiàn)許多新問題，這限制了現(xiàn)有結(jié)霜模型的應(yīng)用；同時，目前對結(jié)霜工況下的星型翅片管的換熱性能的研究不多，且這些研究多關(guān)注于濕空氣溫濕度等環(huán)境參數(shù)對換熱特性和結(jié)霜特性的影響，針對結(jié)構(gòu)參數(shù)對結(jié)霜性能的研究很少。本文基于webb[5]等提出霜層表面濕空氣過飽和假設(shè)建立了星型翅片管換熱器的結(jié)霜模型，并通過UDF導入Fluent多相流模型的控制方程中，實現(xiàn)了對霜層生長過程的數(shù)值模擬，得到了霜層各物性參數(shù)的變化，并對結(jié)構(gòu)參數(shù)對結(jié)霜現(xiàn)象的影響進行了研究，這對于指導星型翅片管換熱器的設(shè)計和優(yōu)化具有一定的指導價值。

2 數(shù)學模型的建立與分析

2.1 物理模型

星型翅片管換熱器的結(jié)構(gòu)簡單，由星型翅片管按照特定的排列方式組成的。其基本結(jié)構(gòu)如圖1所示。翅片管由豎直圓管和在圓管外壁面輻射狀分布的豎直平直翅片組成，其主要結(jié)構(gòu)參數(shù)有翅片管長度L，翅片管內(nèi)徑D，翅片高度H，翅片厚度δ，翅片夾角θ等。

圖1 星型翅片管結(jié)構(gòu)示意圖Fig.1 Diagram of star-shaped finned-tube

為了降低計算量和計算難度，在進行數(shù)值模擬時選取了整個翅片管換熱器中部的一段作為研究對象(如圖2)，并將預先計算出的速度場和溫度場等作為其入口邊界條件。鑒于星型翅片管結(jié)構(gòu)的規(guī)則性，本文采用六面體結(jié)構(gòu)化網(wǎng)格進行了網(wǎng)格劃分，并進行了網(wǎng)格無關(guān)性驗證，最終選取網(wǎng)格數(shù)目為5萬。

圖2 計算域的幾何結(jié)構(gòu)Fig.2 Structure of calculation area

2.2 基本假設(shè)

對實際結(jié)霜過程進行數(shù)值模擬是十分復雜的，因此本章在進行數(shù)值模擬時做了如下假設(shè)：

(1)結(jié)霜過程簡化為水蒸氣在翅片管外壁面上的凝華過程，忽略掉水蒸氣的凝結(jié)過程[6]；

(2)忽略結(jié)霜后期出現(xiàn)的周期性融霜過程[7]；

(3)翅片管內(nèi)壁溫度恒定不變，等同于低溫介質(zhì)的蒸發(fā)溫度；

(4)忽略掉高度和重力的影響，不考慮霜層在豎直方向上的變化。

2.3 模擬方法

本文采用歐拉多相流模型結(jié)合UDF來對結(jié)霜過程進行模擬。在模擬過程中，主相為濕空氣，濕空氣由干空氣和水蒸氣組成，利用組分輸運模型來實現(xiàn)混合過程；冰設(shè)為次相。在結(jié)霜過程中，濕空氣中的水蒸氣會向冰相進行轉(zhuǎn)移而發(fā)生相變，相變過程還存在著相間質(zhì)量源項等的傳遞。

由結(jié)霜機理可知，質(zhì)量源項、能量源項等各個源項都與濕空氣中的水蒸氣到冰相的質(zhì)量傳遞率息息相關(guān)，因此，選取正確的傳質(zhì)模型是十分重要的。webb[5]過飽和水蒸氣假設(shè)指出霜層表面的濕空氣是處于過飽和狀態(tài)的；。基于此假設(shè)，計算出了相變過程中的相間質(zhì)量傳遞率，如式(1)，并將相對濕度和水的冰點溫度作為控制相變的空間判據(jù)。

(1)

式中：A為溫度系數(shù)；β為轉(zhuǎn)換系數(shù),1/m；ρa為濕空氣的密度,kg/m3；wa為濕空氣中水蒸氣的質(zhì)量分數(shù),%；wfs為濕空氣中水蒸氣的飽和質(zhì)量分數(shù),%。

在模型中，霜層的物理結(jié)構(gòu)是濕空氣相和冰相組成的多孔介質(zhì)，結(jié)構(gòu)如圖3，霜層的密度可以通過冰相體積分數(shù)間接求出。

(2)

式中：αi為冰相體積分數(shù)；ρfr為霜層密度,kg/m3。

圖3 網(wǎng)格劃分情況Fig.3 Mesh structure

圖4 霜層結(jié)構(gòu)示意圖Fig.4 Diagram of frost structure

2.4 模擬工況

本文主要研究了翅片高度和翅片厚度對結(jié)霜過程中傳熱傳質(zhì)特性的影響，翅片高度的變化范圍為30—120 mm，翅片個數(shù)的研究范圍為3—24個。翅片管內(nèi)壁選用定壁溫條件，數(shù)值為液化天然氣的沸點溫度111.65 K，空氣溫度為284.15 K，空氣濕度為54%。

3 模擬結(jié)果

3.1 實驗驗證與分析

為了對數(shù)值模擬結(jié)果進行驗證以確保其可靠性，本文先采用趙鵬[8]等人的實驗中的尺寸結(jié)構(gòu)(如表1)進行了數(shù)值模擬，并將數(shù)值模擬得出的霜層厚度與實驗結(jié)果進行了比較。

表1 星型翅片管換熱器結(jié)構(gòu)參數(shù)

圖5為霜層厚度的模擬值和實驗值的對比。由圖可知，兩者的變化趨勢是基本一致的。在結(jié)霜初期霜層生長較快，霜層厚度增加明顯；而隨著結(jié)霜過程的進行，霜層高度的增加量逐漸減小，霜層生長進入穩(wěn)定期。

圖5 霜層厚度的實驗驗證Fig.5 Experimental test of frost formation process

此外，在結(jié)霜初期，模擬所得的霜層厚度要略高于實驗值。這是由于結(jié)霜初期霜層較為稀疏，霜層易被吹走，導致霜層厚度降低；而數(shù)值模擬中忽略了這一過程，因此結(jié)霜初期霜層偏高?？偟膩碚f，模擬結(jié)果和實驗數(shù)據(jù)相比吻合性較好，霜層厚度的最大誤差為18.75%，平均誤差為9.33%。

圖6為不同時刻星型翅片管外側(cè)霜層的變化圖。由圖可知，翅片管外側(cè)的霜層高度和霜層密度都在逐漸增大。在初始階段，霜層厚度和密度增加較快；而隨著結(jié)霜過程的不斷進行，霜層厚度和霜層密度的增加量都在逐漸變小。

圖6 翅片管表面霜層生長示意圖Fig.6 Simulation results of frost formation process

此外，在結(jié)霜初始階段，結(jié)霜現(xiàn)象主要發(fā)生在圓管外側(cè)和翅片翅根處，這是由于在初始階段圓管外側(cè)和翅根處的溫度要相對較低，該處的相變驅(qū)動力相對較大。而隨著結(jié)霜過程的進行，圓管外側(cè)和翅根處和霜層率先進入霜層穩(wěn)定期，此時，結(jié)霜現(xiàn)象主要發(fā)生在翅片表面和翅尖處。這是因為在結(jié)霜后期，翅尖處的厚度相對較小，空隙率也相對較大，在霜層表面和霜層內(nèi)部發(fā)生相變的水蒸氣也較多，因此結(jié)霜后期翅尖處的結(jié)霜現(xiàn)象較為強烈。

3.2 翅片個數(shù)對結(jié)霜過程的影響

3.2.1 翅片個數(shù)對霜層表面溫度的影響

圖7為不同翅片個數(shù)的星型翅片管換熱器的霜層表面溫度隨時間的變化趨勢圖。由圖可知，隨著翅片個數(shù)的增加，在結(jié)霜后期，霜層表面溫度與翅片個數(shù)成正比關(guān)系，即翅片個數(shù)越多，霜層表面溫度越大。但在結(jié)霜過程前期，霜層表面溫度卻與翅片個數(shù)成反比關(guān)系，即翅片個數(shù)越多，霜層表面溫度越小。

圖7 霜層表面溫度隨翅片個數(shù)的變化Fig.7 Temperature of frost with number of fins

這是因為結(jié)霜前期，翅片個數(shù)越多，霜層密度和霜層導熱系數(shù)越大，兩個翅片之間的空氣的溫度越低，因而霜層表面的溫度更低。但隨著結(jié)霜過程進行，霜層的熱阻越來越大，翅片個數(shù)越多，結(jié)霜現(xiàn)象更嚴重，其霜層熱阻也就越大，這導致翅片個數(shù)較多時霜層表面溫度更接近于環(huán)境溫度；因而此時，翅片個數(shù)越多，霜層表面溫度越高。

3.2.2 翅片個數(shù)對結(jié)霜速度的影響

本文中采用傳質(zhì)系數(shù)來代表霜層的生長速率，其計算關(guān)系式如下：

式中：ΔM為濕空氣單位時間內(nèi)的總析濕量,kg/s，即總傳質(zhì)量；LMCD為對數(shù)平均濃度差；A為冷表面的面積,m2。

圖8不同翅片個數(shù)的星型翅片管結(jié)霜過程中的傳質(zhì)系數(shù)隨時間的變化趨勢圖。由圖可知，隨著翅片個數(shù)的增加，霜層的傳質(zhì)系數(shù)逐漸增加，這說明翅片個數(shù)較大時，結(jié)霜現(xiàn)象更為嚴重；但當翅片個數(shù)較多時增加趨勢逐漸減小。這是由于翅片個數(shù)越多，兩個翅片間的空氣的溫度就越低，達到相變條件的水蒸氣就越多，因此傳質(zhì)系數(shù)較大；但隨著翅片個數(shù)的增加，相鄰兩個翅片中的濕空氣的量逐漸減小，這限制了霜層的生長，因此翅片個數(shù)較多時，傳質(zhì)系數(shù)增加較小。

圖8 傳質(zhì)系數(shù)隨翅片的個數(shù)的變化Fig.8 Mass transfer coefficient with number of fins

3.2.3 翅片個數(shù)對空氣側(cè)換熱過程的影響

在結(jié)霜過程中，霜層的積累會造成換熱器流道的堵塞，引起換熱器性能的惡化；但結(jié)霜過程中存在著水蒸氣的凝華潛熱，這又對換熱過程產(chǎn)生了有利影響。在考慮了相變潛熱的前提下，本文對結(jié)霜工況下翅片管外側(cè)的換熱過程進行了研究。

圖9為考慮相變潛熱情況下的空氣側(cè)的換熱系數(shù)的變化趨勢圖。當翅片個數(shù)較多時，空氣側(cè)的換熱系數(shù)在結(jié)霜初期先增加而后減少。這是因為翅片個數(shù)越多，傳質(zhì)系數(shù)越高，相變潛熱的數(shù)值就越大；同時，此時霜層厚度較薄，霜層熱阻的影響較小，這兩者共同作用導致結(jié)霜前期的換熱系數(shù)呈上升趨勢。

圖9 空氣側(cè)換熱系數(shù)隨時間和翅片個數(shù)的變化Fig.9 Heat transfer coefficient with number of fins

如圖10所示，在霜層成熟期，隨著翅片個數(shù)的增加，空氣側(cè)的總換熱系數(shù)不斷減小。由圖9和圖10可知，當翅片個數(shù)為3時，空氣側(cè)總換熱系數(shù)下降了63.1%；而翅片個數(shù)為12時，下降的數(shù)值為79.3%，綜上所述，翅片個數(shù)不宜過多，最多不要超過12個，以防止霜層堵塞流道，造成嚴重的傳熱惡化。

圖10 時間為6 h時空氣側(cè)換熱系數(shù)隨翅片個數(shù)的變化趨勢Fig.10 Heat transfer coefficient with number of fins at 6 h

3.3 翅片高度對結(jié)霜過程的影響

3.3.1 翅片高度對霜層表面溫度的影響

圖11給出了霜層表面溫度隨時間和翅片高度的變化趨勢。由圖可知，在結(jié)霜初期，翅片高度越小，霜層表面溫度越低；而在結(jié)霜后期，翅片高度較小，霜層表面溫度越大。這是因為在結(jié)霜初期，翅片高度越小，結(jié)霜速率越大，霜層密度和導熱系數(shù)也越大，導致其霜層表面的溫度較低； 而在結(jié)霜后期，霜層熱阻越

圖11 霜層表面溫度翅片高度的變化Fig.11 Temperature with number of height of fins

來越大并成為主要熱阻，導致霜表面的溫度更接近空氣溫度，而翅片高度較小時，結(jié)霜現(xiàn)象更為嚴重，因此結(jié)霜末期翅片高度較小時霜層表面溫度更高。

3.3.2 翅片個數(shù)對結(jié)霜速度的影響

圖12為傳質(zhì)系數(shù)隨翅片高度的變化趨勢圖。由圖可知，翅片高度越高，傳質(zhì)系數(shù)越小；同時，隨著翅片高度的增加，傳質(zhì)系數(shù)的下降趨勢也逐漸趨于平緩，這說明增加翅片高度對霜層生長具有抑制作用。產(chǎn)生這種現(xiàn)象的原因是，翅片高度越低，相鄰翅片間的空氣溫度越低，進入翅片間進行相變的水蒸氣的量也越多，這導致結(jié)霜初期低翅片高度的翅片管傳質(zhì)系數(shù)更大。而隨著結(jié)霜過程的進行，發(fā)生相變的水蒸氣的量逐漸減少，因此傳質(zhì)系數(shù)逐漸減小。而翅片高度越小，結(jié)霜現(xiàn)象越嚴重，霜層對傳質(zhì)現(xiàn)象造成的阻礙作用也就越大，因此，低翅片高度的翅片管傳質(zhì)系數(shù)的衰減更大。

圖12 傳質(zhì)系數(shù)隨翅片高度的變化趨勢圖Fig.12 Mass transfer coefficient with height of fins

3.3.3 翅片高度對空氣側(cè)換熱過程的影響

圖13為考慮了相變潛熱的情況下的空氣側(cè)傳熱系數(shù)的變化趨勢。如圖所示，在結(jié)霜初期，空氣側(cè)傳熱系數(shù)先上升而后下降，這是因為結(jié)霜初期，傳質(zhì)系數(shù)較大，相變潛熱數(shù)值也相對較大，此時，相變潛熱的影響要大于霜層積累帶來的傳熱惡化的影響，因此傳熱系數(shù)上升；而隨著結(jié)霜過程的進行，霜層熱阻的影響逐漸增大，空氣側(cè)傳熱系數(shù)隨之變小。在霜層穩(wěn)定期，翅片高度越低，空氣側(cè)傳熱系數(shù)的下降幅度更大，換熱系數(shù)的數(shù)值也越小，如圖14，在翅片高度為30 mm時，空氣側(cè)傳熱系數(shù)降低了72.7%；而翅片高度為90 mm時，該數(shù)值為48.9%。綜上所述，翅片高度較大時，翅片管在結(jié)霜工況下有著更好的傳熱性能。

圖13 空氣側(cè)換熱系數(shù)隨時間和翅片高度的變化Fig.13 Heat transfer coefficient with height of fins

圖14 時間為6 h時空氣側(cè)換熱系數(shù)隨翅片高度的變化趨勢Fig.14 Heat transfer coefficient with height of fins at 6 h

4 結(jié) 論

本文采用對星型翅片管換熱器表面的結(jié)霜過程進行數(shù)值模擬，并對結(jié)構(gòu)參數(shù)對結(jié)霜過程中的傳質(zhì)傳熱特性的影響進行了研究，結(jié)果顯示：

(1)模擬結(jié)果與實驗結(jié)果吻合性良好，霜層厚度的最大誤差為18.75 %，平均誤差為9.33 %。

(2)結(jié)霜初期，翅片管的基管處結(jié)霜速度較快，結(jié)霜末期，結(jié)霜過程主要發(fā)生在翅尖處。

(3)霜層的生長速度和霜層的表面溫度的增加量隨著霜層的生長逐漸變小。

(4)結(jié)霜現(xiàn)象會造成傳熱惡化現(xiàn)象，翅片個數(shù)越多，翅片高度越小，傳熱惡化更嚴重。

1 張書堂，張 潔．液化天然氣空溫式氣化器傳熱過程的數(shù)值分析[J].上海煤氣， 2014(5)：1-4．

Zhang Shutang, Zhang Jie. Numerical Analysis of Heat Transfer Process in LNG Ambient Air Vaporizer [J]. Shanghai Gas, 2014(5)：1-4.

2 陳 彥．冷面結(jié)霜研究及風力機葉片覆冰的數(shù)值模擬[D]．北京：清華大學，2012．

Chen Yan. Study of frost formation and icing simulation of wind turbine blade [D]. Beijing: Tsinghua University, 2012.

3 Chen S, Yao S, Xie F. Analysis of Thermal Conductivity of Frost on Cryogenic Finned-tube Vaporizer Using Fractal Method[J]. Energy & Power Engineering, 2013 (4):109-115.

4 崔 靜. 結(jié)霜與抑霜機理研究及數(shù)值模擬[D]. 大連：大連理工大學, 2011.

Cui Jing. An Investigation of Mechanism for Frosting and Frost Restraint as well as Their Simulations [D]. Dalian: Dalian University of Technology, 2011.

5 Na B, Webb R L. New model for frost growth rate[J]. International Journal of Heat & Mass Transfer, 2004, 47(47):925-936.

6 Hayashi Y, Aoki A, Adachi S, et al. Study of Frost Properties Correlating With Frost Formation Types[J]. Journal of Heat Transfer, 1977, 99(2):239-245.

7 李瑞霞，吳曉敏，王維城．鋁表面結(jié)霜現(xiàn)象的實驗研究[J]．熱科學與技術(shù)，2005(4)：315-318．

8 Li Ruixia, Wu Xiaomin, Wang Weicheng. Experimental study of frost phenomena on aluminum plate [J]. Journal of Thermal Science and Technology, 2005(4)：315-318.

9 趙 鵬，李祥東，汪榮順,等．豎直星型翅片管結(jié)霜的實驗研究[J]．低溫與超導，2009，37(8)：9-13.

Zhao Peng, Li Xiangdong, Wang Rongshun, et al. Experimental study of frosting process on a vertical star-shaped finned-tube [J]. Cryogenics and Superconductivity, 2009,37(8):9-13.

Simulation of frost formation process of star-shaped finned-tube heat exchanger

Pu Liang Zhou Yang Qi Di

(Department of Refrigeration and Cryogenic Engineering，Xi’an Jiaotong University, Xi’an 710049,China)

The frost formation process of star-shaped finned-tube heat exchanger is studied by numerical simulation. In this paper, the frost property is obtained, and the heat and mass transfer characteristics during frost formation process is analyzed, and influences of fin number and fin height on frost formation process is researched as well. The results show that the heat transfer performance will deteriorate seriously due to the frost formation, and the shorter fin height and larger fin number will result in a more severer heat transfer deterioration.

star-shaped finned-tube heat exchanger; numerical simulation; frost formation process; heat and mass transfer characteristics

2016-06-21；

2016-09-29

中央高?；究蒲袠I(yè)務(wù)費專項資金資助。

蒲 亮，男，40歲，副教授、博士生導師。

TB657

A

1000-6516(2016)05-0033-06