鄧志安姜晨薇張雪婷梁　逍楊　帆

1.西安石油大學(xué) 2.中國石化銷售有限公司陜西咸陽石油分公司

?

超級開架式氣化器新型傳熱管內(nèi)流場及對流換熱的數(shù)值模擬

鄧志安1姜晨薇1張雪婷1梁逍1楊帆2

1.西安石油大學(xué) 2.中國石化銷售有限公司陜西咸陽石油分公司

鄧志安等. 超級開架式氣化器新型傳熱管內(nèi)流場及對流換熱的數(shù)值模擬. 天然氣工業(yè), 2016,36(4):90-95.

摘 要超級開架式氣化器（SuperORV）新型傳熱管主要是在普通傳熱管的內(nèi)管中加入了十字螺旋擾流桿，以強化傳熱管的換熱效果。為了研究該裝置的換熱效率，建立了與新型傳熱管和普通傳熱管相對應(yīng)的數(shù)學(xué)和模型物理模型，采用數(shù)值模擬的方法對兩種傳熱管內(nèi)流場及對流換熱性能進行了對比分析，得到了傳熱管內(nèi)流道的努賽爾數(shù)隨入口速度的變化圖，以及傳熱管不同位置處溫度和傳熱系數(shù)的分布規(guī)律。通過對現(xiàn)有實驗進行數(shù)值模擬，對比分析了模擬結(jié)果與實驗結(jié)果。結(jié)果表明：①FLUENT數(shù)值模擬方法能準(zhǔn)確描述傳熱管的傳熱特性；②十字螺旋擾流桿的存在，不僅能夠加強流體的湍流強度，而且還能有效減小邊界層厚度，產(chǎn)生強烈的二次流，加強了流體在徑向上的熱交換，提高了傳熱管的整體換熱能力；③流體入口速度與努賽爾數(shù)的增長幅度呈正比，傳熱系數(shù)隨著溫度的升高而增加，在相同的雷諾數(shù)下，新型傳熱管內(nèi)流體的平均溫度明顯高于普通傳熱管內(nèi)流體溫度。該研究成果能夠為超級開架式氣化器國產(chǎn)化進程提供參考。

關(guān)鍵詞超級開架式氣化器 螺旋擾流桿 強化傳熱 數(shù)值模擬 傳熱系數(shù) 溫度分布 換熱效能

開架式氣化器（ORV）是采用海水作為熱源的氣化器，因其具有操作簡單、運行成本低、生產(chǎn)安全、工作負(fù)荷可調(diào)等優(yōu)點而被廣泛應(yīng)用。ORV在運行過程中，由于LNG溫度過低（－163 ℃），會使得海水溫度急劇下降，為了避免板型管束下部特別是集液管外表面結(jié)冰，提高氣化器的換熱性能，在傳統(tǒng)ORV基礎(chǔ)上，開發(fā)出了SuperORV。相關(guān)資料表明[1]：SuperORV單根換熱管的氣化能力可以提高3倍左右，而海水用量、建造成本和安裝所需空間則可分別減少15％、10％和40％。

自2007年SuperORV技術(shù)首次應(yīng)用于日本大阪的Senboku LNG接收終端并取得成功后，近年來對SuperORV的應(yīng)用正逐步增多，Sekiguchi等[2-3]開發(fā)了一種高性能低成本的ORV設(shè)計方法，并采用該方法進行了低成本的緊湊型ORV設(shè)計。Hisada等[4]通過數(shù)值分析對ORV傳熱管進行了傳熱設(shè)計，得出了ORV傳熱管的最佳型式，建立了完整的ORV優(yōu)化設(shè)計方法并開發(fā)了一種可靠的具有低制造成本和良好運行性能的ORV。張尚文[5]對ORV的結(jié)構(gòu)和工作原理、LNG組分及物性、換熱管結(jié)構(gòu)及工藝設(shè)計等進行了分析介紹。Morimoto等[6]和Yamazaki等[7]分別對SuperORV的結(jié)構(gòu)、工作原理、傳熱特性和發(fā)展現(xiàn)狀進行了深入分析。盡管對SuperORV的應(yīng)用逐步增多，但是對SuperORV的研究仍然不多，大多數(shù)研究工作仍是針對ORV展開的。筆者建立了一種改進的SuperORV新型傳熱管模型，并采用數(shù)值模擬方法分析、探討了SuperORV新型傳熱管的換熱效果。

1 物理模型的建立

超級開架式氣化器新型傳熱管（以下簡稱為B管），是在超級開架式氣化器普通傳熱管（以下簡稱為A管）的基礎(chǔ)上，對內(nèi)套管流道加以改進，在其內(nèi)增加了一個十字螺旋擾流桿（圖1），圖2為其結(jié)構(gòu)示意圖，圖2中a為37.6 mm，b1為56.2 mm，b2 為42.8 mm，b3為53.6 mm。

圖1　螺旋擾流桿示意圖

圖2　SuperORV傳熱管的結(jié)構(gòu)示意圖

2 計算模型

2.1 控制方程

各控制方程[8]的表述如下（i，j = 1，2，3）。

2.1.1 連續(xù)方程

式中ρ表示流體的密度；t表示時間；xi表示i方向；ui表示流體速度沿i方向的分量。

2.1.2 動量守恒方程

式中xj表示j方向；uj表示流體速度沿j方向的分量；p表示靜壓力；τij表示應(yīng)力矢量；ρgi表示i方向的重力分量；Fi表示能源及阻力引起的其他能源項。

2.1.3 能量守恒方程

式中h表示傳熱系數(shù)；k表示湍動能；kτ表示湍流傳遞引起的傳導(dǎo)率；T表示熱力學(xué)溫度；Sh表示定義的體積源。

2.2 湍流模型

由于新型傳熱管內(nèi)加裝了十字螺旋擾流桿，甲烷流體在超級開架式氣化器新型傳熱管內(nèi)的流動會發(fā)生復(fù)雜的旋流以及多尺度二次流。因此，采用Realizable k-ε模型作為湍流模型。

2.2.1 湍動能k方程

式中σk表示k方程的湍流Prandtl數(shù)；μ表示層流黏性系數(shù)；μτ表示湍流黏性系數(shù)；K表示分子熱導(dǎo)率；Gk表示由層流速度梯度而產(chǎn)生的湍流動能；Gb表示由浮力而產(chǎn)生的湍流動能；ε表示耗散率；Ym表示由擴散產(chǎn)生的波動；Sk表示控制方程的源項。

2.2.2 湍流耗散率ε方程

式中τ表示應(yīng)力；σε表示ε方程的湍流Prandtl數(shù)；C2ε表示常量、C1ε表示常量、C3ε表示常量；Sε表示控制方程的源項。

其中

2.3 網(wǎng)格劃分

運用ANSYS ICEM CFD對模型進行網(wǎng)格劃分，采用有限容積法對計算區(qū)域進行離散化，采用正四面體和六面體結(jié)合的方式進行分塊劃分?？紤]到邊界層的影響，對近壁面處進行了網(wǎng)格加密處理，能更準(zhǔn)確地反映近壁面處流體的流動及傳熱特性（圖3）。通過對網(wǎng)格的不斷測試及改進，最終得到了新型傳熱管和普通傳熱管的網(wǎng)格數(shù)分別為150萬個和170萬個。

圖3　模型剖面網(wǎng)格示意圖

2.4 求解及邊界條件的設(shè)置

為了研究氣化器內(nèi)超臨界甲烷沿豎直管向上流動的相變傳熱過程，應(yīng)用壓力基求解器和Realizable k-ε模型，壓力和速度耦合采用SIMPLE方法[9]，動量方程、能量方程、湍流耗散率方程以及湍動能方程均采用二階迎分差格式[10]，甲烷在超臨界條件下的物性參數(shù)采用分段多項式擬合的方法處理，通過FLUENT的UDF[11]給定。UDF根據(jù)物性參數(shù)擬合方程自編的程序，采用C語言編寫，使用DEFINE宏來定義，加載不同工況下甲烷的不同物性參數(shù)，使其動態(tài)地與FLUENT求解器連接進行求解。

邊界條件的設(shè)置如下：①管內(nèi)介質(zhì)為液體甲烷[12-13]（質(zhì)量分?jǐn)?shù)為95％，故作為替代工質(zhì)）；②管材為鋁合金；③進口溫度為130 K，速度入口邊界條件為0.2～2.0 m/s；④壓力出口邊界條件為5～8 MPa；⑤無滑移固定壁面，外翅海水側(cè)液膜的換熱系數(shù)為5 800 W/（m2·K），該數(shù)據(jù)來自與該類氣化器有關(guān)的參考文獻[14]。

3 模擬結(jié)果及分析

3.1 管內(nèi)流場分析

流體進口速度為1 m/s時，A管和B管在沿傳熱管軸線0.2 m截面處的速度矢量局部分布見圖4。

從圖4可以看出，B管中花瓣狀流道內(nèi)的流體速度較A管分布均勻，且中心流域和近壁面流體流速的差異較A管小。由于螺旋擾流桿的離心作用，B管中的徑向運動明顯增強，加大了流體的擾動，從而改善了A管中內(nèi)管壁面處形成的較厚黏性底層，沒有形成滯留區(qū)，使B管整個管道的換熱效率得到大幅度增強。

雷諾系數(shù)為10 000時兩種流道出口處截面的速度場局部云圖見圖5。

圖4　流道中間截面處速度矢量局部分布圖

圖5　傳熱管出口處速度場局部云圖

與圖5-a中A管花瓣狀流道和內(nèi)管流道主流速度場呈現(xiàn)的同心圓相比較，B管內(nèi)由于受到擾流桿的擾流作用，產(chǎn)生旋流以及垂直于主流的二次流，增強了流體在徑向的混合力度及湍流強度，從而削弱了溫度邊界層厚度，使換熱效率得到提高。此外，A管中內(nèi)管壁面處的滯留現(xiàn)象在B管中也得到了很好的改善。

3.2 努賽爾數(shù)與流速的關(guān)系分析

圖6是相同進口速度下，A、B管內(nèi)管流道的努賽爾數(shù)（Nu）隨進口速度（v）的變化圖。從圖6可以看出，B管內(nèi)管的Nu明顯大于A管，且其增加的幅度隨著入口速度的增大而增大。B管內(nèi)的十字螺旋擾流桿促使流道內(nèi)溫度分布呈現(xiàn)較均勻狀態(tài)，同時提高了近壁面處邊界層的溫度梯度。同時，擾流產(chǎn)生得二次流使得徑向擾動程度增強，增大了流體的湍流度，進一步破壞了邊界層，近壁面近乎滯留態(tài)的流體逐步過渡為核心流擾動，使更多流體參與換熱，提高了對流換熱效果，增強了管內(nèi)流體的換熱效率。

圖6　兩種流道的努賽爾數(shù)隨進口速度變化圖

3.3 綜合換熱性能分析

圖7所示為A、B管在雷諾系數(shù)為10 000時管中不同位置流體加熱過程中溫度沿管長方向的分布情況。從圖7可以看出，A、B管溫度沿管長方向升高的程度有所差異，B管內(nèi)管流道中流體的平均溫度較A管高，而內(nèi)管流道壁面平均溫度較A管有所下降，且偏向管內(nèi)低溫流體。說明B管內(nèi)擾流桿使得花瓣狀流道內(nèi)的熱量能夠更好地傳入到內(nèi)管流道中，改善了A管中管壁處形成較厚黏性底層阻礙了熱量傳遞的現(xiàn)象。此外，當(dāng)流體流至管長0.1 m后，花瓣狀流道內(nèi)的流體溫度基本保持在一個比較穩(wěn)定的溫度上，驗證了SuperORV利用環(huán)狀間隙中的流體保溫以緩解管道外壁因溫度驟降而結(jié)冰的狀況。

圖7　A、B管不同位置的平均溫度分布圖

壓力為5 MPa、雷諾系數(shù)為10 000時，超級開架式氣化器新型傳熱管與普通傳熱管傳熱系數(shù)的分布情況見圖8。隨著管內(nèi)主流溫度的升高，B管十字螺旋擾流管和A管內(nèi)管的換熱系數(shù)均升高，且B管十字螺旋擾流管的傳熱系數(shù)較A管內(nèi)管的傳熱系數(shù)大，說明B管的十字螺旋擾流管引起的擾流作用加強了管內(nèi)兩側(cè)流體的換熱，提高了傳熱管的換熱效率。同時，當(dāng)溫度上升到160 K后，傳熱系數(shù)隨之上升并保持在一個較穩(wěn)定的狀態(tài)，說明此后管道內(nèi)流體的溫度變化較為穩(wěn)定，這也與圖5-b中溫度隨管長的分布規(guī)律一致。

圖8　A、B管的傳熱系數(shù)分布規(guī)律圖

從傳熱系數(shù)（h）沿管長（L）的變化規(guī)律（圖9）可以看出，A管和B管換熱系數(shù)均沿管長方向逐漸減小，結(jié)合圖7可以看出，沿管長方向換熱管的傳熱量雖然增加了，但是傳熱系數(shù)降低。因此換熱效率也隨之降低，故并不是換熱管越長越好。

圖9　傳熱系數(shù)沿管長的變化圖

4 數(shù)值模擬結(jié)果正確性驗證

為了驗證新型傳熱管強化傳熱數(shù)值模擬方法的正確性和可靠性，對本文參考文獻[15]中的壓力為23 MPa、質(zhì)量流速為458.8 kg/（m2·s）的工況進行了數(shù)值模擬，通過比對數(shù)值計算結(jié)果與實驗結(jié)果來驗證數(shù)值模擬手段的可靠性。

根據(jù)本文參考文獻[15]建立內(nèi)徑為10 mm、壁厚為3 mm、長為2.5 m的垂直圓管模型，采用非結(jié)構(gòu)自動體網(wǎng)格進行網(wǎng)格劃分，并在壁面處進行網(wǎng)格加密處理。圖10為圓管內(nèi)不同熱流密度情況下，壓力為23 MPa時內(nèi)壁溫度及傳熱系數(shù)隨主流溫度變化的數(shù)值計算結(jié)果，與本文參考文獻[15]中的實驗結(jié)果相近。其中，tb表示主流溫度，tw表示內(nèi)壁溫度。

從圖10可以看出，質(zhì)量流速為458.8 kg/（m2·s）的工況下，在入口不遠(yuǎn)處便出現(xiàn)溫度急劇升高、傳熱惡化的現(xiàn)象。這主要是由于熱流密度較大，近壁面溫度較快地達到擬臨界溫度，此時在壁面處流體的密度較大，形成了阻礙熱量傳遞的黏性底層且流速慢，甚至形成滯留區(qū)包覆在壁面上，從而導(dǎo)致傳熱惡化、壁面溫度突升。當(dāng)?shù)竭_擬臨界點時出現(xiàn)傳熱強化現(xiàn)象，且其強化傳熱效果隨熱流密度的增加逐漸減小，傳熱系數(shù)降低。但傳熱惡化和強化現(xiàn)象均在發(fā)展一段時間后恢復(fù)?？傊?，隨著熱流密度升高，管內(nèi)壁溫度升高，傳熱系數(shù)減小，傳熱特性逐漸變差。因此，減小熱流密度能改善傳熱效果，降低壁面溫度。

以上結(jié)論與本文參考文獻[15]中的實驗結(jié)果相符，表明采用FLUENT數(shù)值模擬的方法能夠準(zhǔn)確地描述傳熱管的傳熱特性。

5 結(jié)論

1）超級開架式氣化器新型傳熱管內(nèi)流道增加了十字螺旋擾流桿，不僅增大了換熱面積，還具有一定的擾流作用，提高了傳熱管的換熱效果。

2）流體入口速度越大，超級開架式氣化器新型傳熱管努賽爾數(shù)的增長幅度越大，其換熱效果比普通傳熱管的換熱效果越好。

3）在相同雷諾數(shù)的條件下，新型傳熱管內(nèi)流道的流體溫度明顯高于普通傳熱管內(nèi)管的流體溫度，且管壁溫度趨于內(nèi)流道溫度?；ò隊盍鞯纼?nèi)的溫度升至一定位置后保持較穩(wěn)定的狀態(tài)，驗證了SuperORV利用環(huán)狀間隙中的流體保溫以緩解管道外壁因溫度驟降而結(jié)冰的狀況。

4）換熱管并非越長越好，沿著管長方向雖然傳熱量增加，但是換熱系數(shù)減小，換熱效率降低。

5）隨著溫度的升高，超級開架式氣化器傳熱管傳熱系數(shù)逐漸增加，且新型傳熱管的傳熱系數(shù)大于普通傳熱管的傳熱系數(shù)；當(dāng)溫度上升至165 K后，管道內(nèi)流體的溫度將保持在一個較為穩(wěn)定的范圍內(nèi)。

參 考 文 獻

[1] 陳永東, 陳學(xué)冬. LNG成套裝置換熱器關(guān)鍵技術(shù)分析[J]. 天然氣工業(yè), 2010, 30(1): 96-100. Chen Yongdong, Chen Xuedong. Technology analysis of heat exchanger in large LNG plant and terminal[J]. Natural Gas Industry, 2010, 30(1): 96-100.

[2] Sekiguchi M, Hisada N. Development of optimum design method and high performance ORV for cost down[C]//The 13thInternational Conference & Exhibition on LNG, LNG13 Poster Session, 2001.

[3] Sekignchi M, Hisada N. Achievement of drastic downsizing and cost reduction for LNG open rack vaporizers by using optimum design methods[C]//The 22ndWorld Gas Conference, 2003.

[4] Hisada N, Sekiguchi M. Design and analysis of open rack LNG vaporizer[C]//Proceedings of Design and Analysis of Pressure Vessels, Heat Exchangers, and Piping Components, 25-29 July 2004, San Diego, California, USA.

[5] 張尚文. 液化天然氣開架式氣化器工藝研究和設(shè)計[J]. 石油化工設(shè)備, 2012, 41(3): 25-29. Zhang Shangwen. Process study and design liquefed natural gas open rack vaporizer[J]. Petro-Chemical Equipment, 2012, 41(3): 25-29.

[6] Morimoto N, Yamamoto S, Yamasaki Y, Shimokawatoko T, Shinkai K, Egashira S, et al. Development and practical application of a high performance open-rack LNG vaporizer (SuperORV)[R]. Osaka: Osaka Gas Co., Ltd., 2003.

[7] Yamazaki K, Shimokawatoko T, Yamasaki Y. Development of a new type of open-rack LNG vaporizer[R]. Osaka: Osaka Gas Co., Ltd. & Kobe Steel, Ltd., 2003.

[8] 孫幫成, 李明高. ANSYS FLUENT 14.0仿真分析與優(yōu)化設(shè)計[M]. 北京: 機械工業(yè)出版社，2014. Sun Bangcheng, Li Minggao. ANSYS FLUENT 14.0 simulation analysis and optimization design[M]. Beijing: China Machine Press, 2014.

[9] 李軍, 王晨, 桑芝富, 史永紅. 強化熱管管內(nèi)傳熱及阻力性能的數(shù)值研究[J]. 石油機械, 2009, 37(2): 37-40. Li Jun, Wang Chen, Sang Zhifu, Shi Yonghong. Numerical studies of heat transfer and resistance performance in strengthening heat pipe tube[J]. China Petroleum Machinery, 2009, 37(2): 37-40.

[10] Wang CC, Chi KY, Chang YJ. An experimental study of heat transfer and friction characteristics of typical lower fn-and-tube heat exchangers[J]. International Journal of Heat & Mass Transfer, 1998, 41(4): 817-822.

[11] 周俊杰，徐國權(quán)，張華俊. FLUENT工程技術(shù)與實例分析[M].北京: 中國水利水電出版社，2010. Zhou Junjie, Xu Guoquan，Zhang Huajun. FLUENT software engineering techniques and case analysis[M]. Beijing: China Water Power Press, 2010.

[12] Miana M, Del Hoyo R, Rodrigalvarez V, Valdes JR, Llorens R. Calculation models for prediction of liquefed natural gas ageing during ship transportation[J]. Applied Energy, 2010, 87(5): 1687-1700.

[13] 余美玲. 新型開架式內(nèi)螺旋管汽化器的數(shù)值模擬和優(yōu)化[D].廣州: 華南理工大學(xué), 2012. Yu Meiling. Numerical simulation and optimization of a new open-frame vaporizer with inner spiral tubes[D]. Guangzhou: South China University of Technology, 2012.

[14] 顧安忠. 液化天然氣技術(shù)[M]. 北京: 機械工業(yè)出版社, 2004. Gu Anzhong. The liquefed natural gas technology［M］. Beijing：China Machine Press, 2004.

[15] 王飛, 楊鈺, 顧漢洋, 趙萌, 李虹波, 盧東華. 垂直管內(nèi)超臨界水傳熱實驗研究[J]. 原子能科學(xué)技術(shù), 2013, 47(6): 933-939. Wang Fei, Yang Yu, Gu Hanyang, Zhao Meng, Li Hongbo, Lu Donghua. Experimental research on heat transfer performance of supercritical water in vertical tube[J]. Atomic Energy Science and Technology, 2013, 47(6): 933-939.

(修改回稿日期 2016-01-10 編 輯 何 明)

Numerical simulation analysis on the fow feld and convection heat transfer in a new heat transfer tube of SuperORV

Deng Zhi'an1, Jiang Chenwei1, Zhang Xueting1, Liang Xiao1, Yang Fan2
(1. Xi’an Shiyou University， Xi’an, Shaanxi 710065, China; 2. Shaanxi Xianyang Petroleum Company of Sinopec Sales Co., Ltd.， Xi’an, Shaanxi 710065, China)

NATUR. GAS IND. VOLUME 36, ISSUE 4, pp.90-95, 4/25/2016. (ISSN 1000-0976; In Chinese)

Abstract:A new heat transfer tube of SuperORV (super open rack vaporizer) is developed by adding a Phillips screw spoiler lever to the inner tube of the traditional heat transfer tube to strengthen the heat transfer efficiency. Mathematical and physical models corresponding to the new and the traditional heat transfer tubes were built to study the heat transfer efficiency of this new device. Two types of heat transfer tubes were compared in terms of their flow fields and convection heat transfer performances by means of numerical simulation. As a result, the relationship diagram of inlet velocity vs. Nusselt Number of flow path in heat transfer tubes was worked out, and the distribution laws of temperature and heat transfer coefficient at different positions of heat transfer tubes were drawn. After numerical simulation was carried out on the existing experiments, simulation results and experimental data were compared and analyzed. It is shown that based on the FLUENT numerical simulation method, the heat transfer characteristics can be described accurately. By virtue of screw spoiler levers, not only the turbulence intensity of the fluid is enhanced, but the boundary layer thickness is reduced effectively with a strong secondary flow, so the radial heat exchange of the fluid is strengthened and the overall heat transfer capacity of heat transfer tubes is improved. Inlet flow velocity is proportional to the increasing rate of Nusselt Number, and heat transfer coefficient increases with the rise of temperature. With the same Reynolds number, the average temperature of the fluid inside the new heat transfer tube is significantly higher than that in the traditional transfer tube. These results provide reference for the localization of SuperORV.

Keywords:SuperORV; Screw spoiler lever; Enhanced heat transfer; Numerical simulation; Heat transfer coefficient; Temperature distribution; Heat transfer efficiency

DOI:10.3787/j.issn.1000-0976.2016.04.014

基金項目：國家自然科學(xué)基金項目“開架式氣化器新型傳熱管強化傳熱機理研究及結(jié)構(gòu)優(yōu)化”（編號：51304160）。

作者簡介：鄧志安，1962年生，教授，博士；主要從事油氣儲運方面的教學(xué)與科研工作。地址：（710065）陜西省西安市西安石油大學(xué)石油工程學(xué)院。電話：（029）88382673，18629351539。ORCID:0000-0003-0459-3088。E-mail:1064514507@qq.com