韓昌亮，任婧杰，王焱慶，董文平，畢明樹(shù)

?

SCV耦合傳熱特性實(shí)驗(yàn)研究與數(shù)值模擬

韓昌亮，任婧杰，王焱慶，董文平，畢明樹(shù)

（大連理工大學(xué)化工機(jī)械學(xué)院，遼寧大連 116024）

浸沒(méi)燃燒式汽化器（SCV）是液化天然氣（LNG）接收站中一種必不可少的換熱設(shè)備，主要通過(guò)水浴系統(tǒng)作為中間介質(zhì)實(shí)現(xiàn)煙氣與LNG之間的熱量傳遞。搭建一套完整的SCV流動(dòng)換熱實(shí)驗(yàn)平臺(tái)，對(duì)其內(nèi)部復(fù)雜的傳熱特性進(jìn)行研究?？梢暬瘜?shí)驗(yàn)結(jié)果揭示了汽化器內(nèi)部一些獨(dú)特的流體動(dòng)力學(xué)現(xiàn)象（局部水浴結(jié)冰等），同時(shí)通過(guò)建立的氣液兩相混合物與跨臨界LNG耦合傳熱計(jì)算模型得到了換熱管束內(nèi)外局部流體溫度和局部傳熱系數(shù)分布曲線(xiàn)，并分析了LNG進(jìn)口壓力、LNG入口速度、初始水位高度以及煙氣進(jìn)氣量對(duì)NG出口溫度和水浴溫度的影響規(guī)律。研究成果能夠?yàn)镾CV國(guó)產(chǎn)化設(shè)計(jì)提供重要參考。

浸沒(méi)燃燒式汽化器；LNG；實(shí)驗(yàn)；流動(dòng)；傳熱；水浴溫度；數(shù)值模擬

引 言

隨著我國(guó)液化天然氣(LNG)進(jìn)口量逐年增加，人們對(duì)于接收站中的汽化設(shè)備愈發(fā)關(guān)注。據(jù)了解，國(guó)內(nèi)接收站在役LNG汽化器主要有開(kāi)架式汽化器(ORV)、浸沒(méi)燃燒式汽化器（SCV）和中間介質(zhì)型汽化器（IFV）等[1]。其中，SCV是一種利用浸沒(méi)燃燒技術(shù)輔以多相流換熱和跨臨界換熱技術(shù)開(kāi)發(fā)的新型高效換熱設(shè)備，具有結(jié)構(gòu)緊湊、操作靈活和換熱效率高等優(yōu)點(diǎn)[2]，廣泛用于接收站調(diào)峰系統(tǒng)。

SCV運(yùn)行時(shí)，燃料氣和空氣的混合氣在燃燒室內(nèi)被點(diǎn)燃，產(chǎn)生的高溫?zé)煔馔ㄟ^(guò)氣體分布器排入到水箱中，與水浴實(shí)現(xiàn)直接接觸式換熱，高速湍動(dòng)的水浴將煙氣熱量傳遞給換熱管束中的LNG。由于接收站工藝要求，通常需要對(duì)LNG進(jìn)行加壓處理，導(dǎo)致SCV管程操作壓力一般大于其臨界壓力。LNG在盤(pán)管內(nèi)經(jīng)歷了由液態(tài)到超臨界態(tài)的工藝轉(zhuǎn)化[3]。由于煙氣攜帶的能量基本能滿(mǎn)足LNG汽化所需熱量，水浴溫度基本保持恒定[4]。圖1給出了典型SCV內(nèi)部流體溫-熵圖。SCV傳熱過(guò)程主要包括兩個(gè)環(huán)節(jié)，即管外氣液兩相流橫掠管束換熱和管內(nèi)跨臨界LNG強(qiáng)制對(duì)流換熱。本研究?jī)?nèi)容為二者耦合傳熱過(guò)程。

目前，國(guó)內(nèi)外學(xué)者針對(duì)SCV流體傳熱性能報(bào)道較少。Park等[5]利用最小熵增方法在恒定壁面溫度和恒定熱通量假設(shè)條件下給出了汽化器優(yōu)化設(shè)計(jì)方法。齊超等[6]利用簡(jiǎn)化的一維傳熱模型分析了某SCV不同工況下的運(yùn)行特性。李仲珍等[7]、張康等[8]和靳書(shū)武等[9]利用CFD數(shù)值模擬方法對(duì)換熱管內(nèi)跨臨界LNG傳熱特性進(jìn)行探討，研究了不同操作參數(shù)對(duì)流體傳熱系數(shù)的影響規(guī)律。李泓鈺[10]以常溫水為管程換熱介質(zhì)，對(duì)水浴及排煙特性進(jìn)行了實(shí)驗(yàn)研究。然而上述研究均基于一定假設(shè)條件，得到的結(jié)論往往無(wú)法直接適用于實(shí)際SCV。同時(shí)，國(guó)內(nèi)使用的SCV均是從日本和德國(guó)進(jìn)口[11]，因此研究其耦合傳熱過(guò)程對(duì)我國(guó)LNG事業(yè)已經(jīng)迫在眉睫。

本工作首先搭建一套完整的可視化SCV傳熱性能測(cè)試平臺(tái)，針對(duì)其內(nèi)部復(fù)雜的換熱特性進(jìn)行實(shí)驗(yàn)研究，并且基于觀察到的現(xiàn)象建立耦合求解煙氣與水浴兩相混合物和跨臨界LNG之間耦合傳熱模型，對(duì)不同操作工況下SCV運(yùn)行特性進(jìn)行了分析?？蔀镾CV工程設(shè)計(jì)提供大量基礎(chǔ)性數(shù)據(jù)。

1 SCV可視化實(shí)驗(yàn)系統(tǒng)

1.1 實(shí)驗(yàn)系統(tǒng)

如圖2所示，本實(shí)驗(yàn)系統(tǒng)主要由浸沒(méi)燃燒器、水箱、煙氣分布器、換熱管束、低溫增壓泵、杜瓦罐和測(cè)量系統(tǒng)等組成。有關(guān)實(shí)驗(yàn)臺(tái)具體尺寸見(jiàn)文獻(xiàn)[12]。需要指出的是，出于安全考慮，實(shí)驗(yàn)過(guò)程中利用跨臨界液氮（LN2）替代LNG作為管程介質(zhì)。煙氣進(jìn)出口溫度、管程流體進(jìn)出口溫度、水浴溫度和換熱管壁面溫度采用PT100熱電阻進(jìn)行測(cè)量。管程操作壓力由壓力變送器進(jìn)行實(shí)時(shí)記錄。LN2流量由高壓旋進(jìn)渦旋流量計(jì)測(cè)得。另外，在水箱和圍堰前后側(cè)開(kāi)有透明視窗，采用高速攝像機(jī)(FASTCAM SA4)觀察殼程流體流動(dòng)狀態(tài)。

1—air blower; 2—volumetric flow meter; 3—ignition system; 4—needle valve; 5—flue tank (CH4); 6—submerged burner; 7—exhaust stack; 8—computer; 9—relief valve; 10—high-pressure ball valve; 11—vortex flowmeter; 12—one-way valve; 13—water tank; 14—weir; 15—flue gas distributor; 16—ball valve; 17—booster pump; 18—dewar; 19—high-speed camera

1.2 實(shí)驗(yàn)結(jié)果分析

圖3顯示了SCV殼程流場(chǎng)。從拍攝結(jié)果來(lái)看，SCV在開(kāi)車(chē)啟動(dòng)階段局部水浴出現(xiàn)了結(jié)冰現(xiàn)象，結(jié)冰部位主要位于底部換熱管區(qū)域和進(jìn)液總管處。這一方面是因?yàn)樵谠搮^(qū)域內(nèi)管內(nèi)流體溫度極低（大約110 K），管外煙氣熱量無(wú)法滿(mǎn)足管內(nèi)流體升溫需求，導(dǎo)致局部水浴溫度急劇下降至冰點(diǎn)以下。另一方面是由于煙氣對(duì)水浴擾動(dòng)作用較弱，水浴出現(xiàn)了局部流動(dòng)“死區(qū)”，此時(shí)應(yīng)調(diào)整系統(tǒng)參數(shù)來(lái)消除結(jié)冰現(xiàn)象。因?yàn)槿绻麚Q熱管壁面長(zhǎng)期存在絕熱邊界條件，將導(dǎo)致系統(tǒng)換熱量大大降低。同時(shí)可以看出水浴中存在有大量換熱氣泡，煙氣與水浴之間以及兩相混合物與管壁之間均發(fā)生直接接觸式傳熱，水浴換熱效率極高[13]。這種狀態(tài)保證了煙氣能將其攜帶的顯熱和水蒸氣冷凝的潛熱全部傳遞給水浴。最終殼程形成了氣液兩相流橫掠管束流動(dòng)形態(tài)。

表1給出了典型SCV運(yùn)行實(shí)驗(yàn)數(shù)據(jù)?？梢钥闯?，管程入口壓力大于LN2臨界壓力，而且管程流體出口溫度遠(yuǎn)大于LN2臨界溫度，因此LN2實(shí)現(xiàn)了跨臨界汽化工藝。在給定燃料量和助燃空氣量條件下，管程出口溫度和殼程對(duì)流傳熱系數(shù)隨初始水位高度增加呈現(xiàn)先升高后降低趨勢(shì)。當(dāng)水位高度從300 mm增至400 mm時(shí)，管程流體出口溫度升高了6 K左右。但是，當(dāng)水位高度繼續(xù)增加到450 mm時(shí)，出口流體溫度反而降低了4 K。因此，在該進(jìn)氣量條件下最佳匹配初始水位高度為400 mm。水浴溫度隨初始水位高度增加單調(diào)降低。

表1 SCV運(yùn)行實(shí)驗(yàn)數(shù)據(jù)

2 SCV耦合傳熱特性數(shù)值模擬

2.1 跨臨界LNG物性分析

LNG為多組分混合物，其主要成分是CH4。本研究使用的LNG組分含量見(jiàn)表2。LNG臨界壓力和臨界溫度分別為4.55 MPa和190.4 K，管程操作壓力一般為5～10 MPa，進(jìn)出口溫度為118～274 K，因此掌握跨臨界LNG物性對(duì)分析SCV換熱過(guò)程至關(guān)重要。圖4給出了利用REFPROP軟件[14]得到的LNG物性曲線(xiàn)。可以看出，在一定壓力下比熱容隨溫度升高先增大后減小，在臨界溫度附近達(dá)到最大值。并且隨壓力增大比熱容峰值變小。密度、熱導(dǎo)率和黏度均隨溫度升高而下降，在臨界點(diǎn)附近下降斜率最大。

表2 LNG各組分含量

2.2 物理模型

為了更好地反映SCV管束特征，本研究在構(gòu)建物理模型時(shí)選取兩根蛇形換熱管為管程計(jì)算區(qū)域，建立了如圖5所示的三維物理模型，此外還包括水箱、煙氣分布器和圍堰等。模型尺寸見(jiàn)表3。

2.3 數(shù)學(xué)模型與網(wǎng)格劃分

采用VOF多相流模型與組分輸運(yùn)模型相結(jié)合方法描述煙氣和水浴兩相混合物與換熱管束之間的傳熱。限于篇幅原因，相關(guān)的控制方程不再贅述，具體數(shù)學(xué)方程可以參見(jiàn)本題課組之前的工作[15-16]。

表3 計(jì)算區(qū)域主要幾何參數(shù)

采用Gambit軟件對(duì)圖5計(jì)算區(qū)域進(jìn)行劃分。為保證計(jì)算精度，首先對(duì)其進(jìn)行分塊處理，除分布器排氣孔附近和彎管處采用非結(jié)構(gòu)網(wǎng)格外，其余部分均采用結(jié)構(gòu)化網(wǎng)格。劃分情況如圖6所示。以0.001、0.002、0.004 m 3種不同網(wǎng)格尺度對(duì)管程局部LNG升溫情況進(jìn)行考察，以此驗(yàn)證網(wǎng)格獨(dú)立性。發(fā)現(xiàn)采用后兩種網(wǎng)格時(shí)模擬計(jì)算結(jié)果十分接近，最終采用0.002 m的網(wǎng)格尺寸，整個(gè)模型網(wǎng)格約為660萬(wàn)個(gè)。

2.4 邊界條件與計(jì)算方法

采用商業(yè)軟件Fluent14.5對(duì)非穩(wěn)態(tài)傳熱過(guò)程進(jìn)行數(shù)值計(jì)算。管程入口為速度入口邊界(t,in=120 K)，管程出口為壓力出口邊界。煙氣入口為速度入口邊界 (s,in=973 K)，煙氣出口為壓力出口邊界。換熱管定義為Coupled邊界條件。其余各面為無(wú)滑移絕熱條件。計(jì)算區(qū)域前、后面均為Symmetry條件，以便該模型可以代表整個(gè)SCV管束區(qū)域。

選取PISO算法進(jìn)行壓力速度耦合求解、Geo-Reconstruct方法追蹤兩相流體界面，壓力項(xiàng)選擇Body Force Weighted 進(jìn)行離散。松弛因子采用默認(rèn)值，動(dòng)量方程和能量方程采用二階迎風(fēng)格式進(jìn)行迭代，時(shí)間步長(zhǎng)選為5×10-5s。當(dāng)NG出口溫度和水浴溫度不隨時(shí)間變化時(shí)，近似認(rèn)為計(jì)算工況收斂。

本研究涉及的管內(nèi)局部傳熱系數(shù)和管外局部傳熱系數(shù)計(jì)算公式如下

(2)

式中，wi和wo分別為作用在管束內(nèi)、外壁面的熱通量，wi和wo分別為管束內(nèi)、外壁面溫度，L和b分別為局部LNG溫度和水浴溫度。

2.5 模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比驗(yàn)證

圖7給出了初始水位高度為400 mm、煙氣量為60 m3·h-1、管程入口壓力為4.5 MPa時(shí)模擬得到的管程流體出口溫度與實(shí)驗(yàn)值對(duì)比結(jié)果?？梢钥闯龆咦兓厔?shì)基本相同。實(shí)驗(yàn)中管程入口流量由45kg·h-1增至108kg·h-1時(shí)，管程流體出口溫度從304.0 K下降到272.7 K，相應(yīng)的數(shù)值模擬結(jié)果為從298.6 K下降到261.2 K。數(shù)值模擬結(jié)果低于實(shí)驗(yàn)值，兩者誤差在±20%以?xún)?nèi)。分析原因，主要是數(shù)值模擬中忽略了熱損失和煙氣熱輻射傳熱，并且忽略了管程壓力降對(duì)流體物性的影響?？紤]到SCV耦合傳熱過(guò)程的復(fù)雜性，該模擬精度滿(mǎn)足工程需求。因此可以利用以上模型開(kāi)展進(jìn)一步耦合傳熱分析。

表4 SCV數(shù)值模擬主要參數(shù)

3 結(jié)果與討論

3.1 SCV基本傳熱特性分析

本研究首先以表4工況為例對(duì)SCV基本傳熱特性進(jìn)行分析。圖8～圖10分別顯示了SCV管程溫度云圖、殼程流場(chǎng)和殼程溫度云圖。

從圖8可以看出管程LNG進(jìn)、出口溫度分別是120 K和273 K，出口溫度高于汽化器基本需求。這是因?yàn)楦咝∠到y(tǒng)將煙氣攜帶的熱量傳遞給管內(nèi)LNG，使其完成了從液態(tài)到超臨界態(tài)的轉(zhuǎn)化。由于LNG在換熱管內(nèi)不斷吸熱膨脹，其速度逐漸增加，該工況下進(jìn)、出口速度分別為0.1 m·s-1和2.0 m·s-1。圖9顯示了殼程流場(chǎng)，可見(jiàn)在煙氣噴入水浴后圍堰內(nèi)產(chǎn)生了很多尺寸不一的換熱氣泡（藍(lán)色和紅色分別代表氣體和液態(tài)水），氣泡群對(duì)水浴的擾動(dòng)作用非常強(qiáng)烈。過(guò)熱氣泡的存在不僅加大了氣液兩相之間的接觸面積，也增強(qiáng)了近壁面處兩相混合物與管壁對(duì)流傳熱效果，這一流體動(dòng)力學(xué)現(xiàn)象是SCV殼程熱交換機(jī)理所在。圖10顯示了殼程溫度云圖，高溫?zé)煔鈬娙胨『鬁囟妊杆傧陆担魵庥隼溽尫诺臐摕嶂饕性谠搮^(qū)域。水浴吸收了煙氣的熱量，又很快將熱量傳遞給管束內(nèi)LNG，最終導(dǎo)致水浴除了在煙氣分布器處附近存在一定高溫區(qū)外其余區(qū)域溫度基本恒定。該結(jié)果與文獻(xiàn)[3]中的結(jié)論基本一致。

從圖11可知，管內(nèi)局部傳熱系數(shù)沿管長(zhǎng)方向先增大后減小，傳熱系數(shù)在準(zhǔn)臨界點(diǎn)處達(dá)到最大值，原因是在擬臨界區(qū)域LNG比定壓熱容達(dá)到最大值，對(duì)傳熱起到了加強(qiáng)作用[17]。而殼程傳熱系數(shù)呈現(xiàn)出一定的波動(dòng)趨勢(shì)，平均殼程傳熱系數(shù)可以達(dá)到4950 W·m-2·K-1，這是因?yàn)榇罅繗馀莶粩鄶_動(dòng)水浴，同時(shí)氣液兩相混合物以較高速度沖刷管束，使得SCV殼程換熱效率很高，同一截面處管外傳熱系數(shù)幾乎是管內(nèi)傳熱系數(shù)的2倍，因此SCV換熱熱阻主要存在于管內(nèi)，該工況下SCV整體傳熱系數(shù)可以達(dá)到850 W·m-2·K-1。此外還可以發(fā)現(xiàn)管內(nèi)傳熱系數(shù)在彎管區(qū)域均會(huì)出現(xiàn)突變，這是因?yàn)長(zhǎng)NG在流經(jīng)彎管時(shí)在離心力作用下出現(xiàn)了二次流現(xiàn)象[18]，該現(xiàn)象有助于減薄流體熱邊界層厚度，提高流體傳熱系數(shù)，所以彎管區(qū)域流體傳熱系數(shù)大于直管段。

圖12顯示LNG溫度、管內(nèi)外壁溫度均沿流動(dòng)方向單調(diào)升高。特別地，在<1.15 m時(shí)，LNG為純液態(tài)，比定壓熱容較小，熱導(dǎo)率較大，因此在該區(qū)域內(nèi)流體升溫較快。在1.15 m<<3.5 m時(shí)，LNG處于超臨界區(qū)域，其熱導(dǎo)率和比定壓熱容均較小，升溫速率相比液態(tài)時(shí)慢一些。在整個(gè)換熱管區(qū)域內(nèi)，對(duì)于管壁溫度而言，沒(méi)有明顯的階躍存在，這主要?dú)w功于管外相對(duì)均勻的氣泡分布和管內(nèi)跨臨界LNG的良好流動(dòng)換熱性能。換熱管內(nèi)外壁溫差以及內(nèi)壁與LNG之間溫差沿管長(zhǎng)方向逐漸減小，這是因?yàn)樵浇咏艹坛隹贚NG升溫所需熱量越小，作用于管壁的熱通量也越小。除此之外，殼程局部水浴溫度基本保持不變。在該工況下煙氣出口溫度為295 K，所以SCV內(nèi)部流體溫度滿(mǎn)足以下順序：煙氣出口溫度>水浴溫度>NG出口溫度。

3.2 LNG進(jìn)口壓力對(duì)傳熱的影響

不同LNG進(jìn)口壓力主要影響管內(nèi)流體熱物性。本研究對(duì)5.43、6.93、8.38、9.55 MPa 4種不同進(jìn)口壓力下LNG升溫特性及水浴溫度進(jìn)行了分析。

圖13顯示了不同LNG進(jìn)口壓力對(duì)SCV傳熱特性的影響。從圖13可知，在<0.75 m和2.5 m<<3.5m時(shí)，LNG溫度基本不隨LNG進(jìn)口壓力改變而發(fā)生變化。而在擬臨界區(qū)域，壓力越高，LNG比定壓熱容峰值越小，局部傳熱系數(shù)也越小，因此LNG升溫更快。壓力升高，NG出口溫度略微升高。當(dāng)其他參數(shù)保持不變時(shí)，水浴溫度基本不隨LNG進(jìn)口壓力變化而改變。當(dāng)LNG進(jìn)口壓力從5.65 MPa增加到9.35 MPa時(shí)，水浴溫度改變了0.35 K。

3.3 LNG進(jìn)口速度對(duì)傳熱的影響

LNG進(jìn)口速度不同將產(chǎn)生不同的管內(nèi)流動(dòng)狀態(tài)。本研究對(duì)0.05、0.10、0.15、0.20 m·s-14種不同LNG進(jìn)口速度下管內(nèi)LNG溫度分布及水浴溫度進(jìn)行了分析。

圖14顯示了不同LNG進(jìn)口速度對(duì)SCV性能的影響?？梢钥闯觯S著LNG進(jìn)口速度的增加，在同一截面上LNG溫度變小。這是因?yàn)椋琇NG進(jìn)口速度增加，換熱管內(nèi)流體湍流作用加劇，有助于減薄流體熱邊界層厚度，增大表面?zhèn)鳠嵯禂?shù)[19]，因此LNG進(jìn)口速度增加可以有效地加強(qiáng)傳熱。NG出口溫度隨LNG進(jìn)口速度增加呈現(xiàn)逐漸降低趨勢(shì)。另外，LNG進(jìn)口速度增加意味著LNG將從殼程水浴吸收更多熱量，所以管外水浴溫度逐漸下降。經(jīng)計(jì)算，LNG進(jìn)口速度從0.05 m·s-1增加到2 m·s-1時(shí)，水浴溫度下降了8 K。

3.4 初始水位高度對(duì)傳熱的影響

在給定煙氣量情況下，不同初始水位高度主要影響殼程傳熱面積。本研究對(duì)300、350、400、450 mm 4種不同初始水位高度下管內(nèi)LNG升溫情況和水浴溫度進(jìn)行了分析。

圖15顯示了不同初始水位高度對(duì)SCV性能的影響?？梢钥闯?，隨著初始水位高度的增加 (300～400 mm)，管內(nèi)LNG升溫加快。這是因?yàn)闅こ虃鳠崦娣e增加導(dǎo)致LNG更容易從水浴系統(tǒng)中吸收來(lái)自煙氣的熱量。當(dāng)初始水位高度增加到450 mm時(shí)，水浴系統(tǒng)的體積過(guò)大，使其受到煙氣的擾動(dòng)程度減弱，LNG吸收到的熱量減少，最終導(dǎo)致管程N(yùn)G出口溫度下降。由于水浴體積的增加，同樣煙氣熱量進(jìn)入水浴后，水浴溫度呈現(xiàn)單調(diào)遞減趨勢(shì)。

3.5 煙氣量對(duì)傳熱的影響

在給定初始水位高度情況下，不同煙氣進(jìn)氣量主要影響水浴中含氣率和水浴流動(dòng)傳熱狀態(tài)。本研究對(duì)45、65、88、110 m3·h-14種不同煙氣量下管程LNG升溫情況和水浴溫度進(jìn)行了分析。

圖16顯示了不同煙氣量對(duì)系統(tǒng)性能的影響。在一定范圍內(nèi)，隨著煙氣量的增加，管內(nèi)LNG升溫更迅速。這是因?yàn)?，進(jìn)氣量增加，水浴內(nèi)部氣含率加大，有助于加劇水浴湍流程度和減薄流體邊界層，煙氣的熱量更容易傳遞給管內(nèi)LNG。繼續(xù)增加煙氣量，導(dǎo)致殼程氣含率和不凝性氣體含量過(guò)大，水浴與管壁接觸面積急劇減小，殼程傳熱過(guò)程主要受氣體支配，而氣體與管壁之間傳熱能力遠(yuǎn)小于液體與管壁之間[20]，水浴出現(xiàn)傳熱惡化現(xiàn)象，導(dǎo)致NG出口溫度下降。而煙氣量增加意味更多熱量被輸入到水浴中，因此水浴溫度逐漸升高。

4 結(jié) 論

采用實(shí)驗(yàn)與數(shù)值模擬相結(jié)合方法對(duì)SCV內(nèi)部流體耦合傳熱特性進(jìn)行研究，利用實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了模型的可靠性，得到的主要結(jié)論如下。

（1）由于管程流體超低溫特性，SCV進(jìn)液總管和底部換熱管附近水浴極易形成冰層，實(shí)際SCV操作中應(yīng)調(diào)整好系統(tǒng)參數(shù)來(lái)消除結(jié)冰現(xiàn)象。這是因?yàn)槿绻L(zhǎng)期存在冰層會(huì)導(dǎo)致系統(tǒng)換熱量大大降低。

（2）管內(nèi)LNG局部傳熱系數(shù)沿管長(zhǎng)方向先增加后減小，在擬臨界點(diǎn)附近達(dá)到最大值。管外局部傳熱系數(shù)由于大量氣泡的擾動(dòng)作用呈現(xiàn)波動(dòng)趨勢(shì)。同一截面上，管外傳熱系數(shù)遠(yuǎn)大于管內(nèi)傳熱系數(shù)，換熱熱阻主要存在于管束內(nèi)部。SCV換熱管束壁溫沒(méi)有明顯的階躍現(xiàn)象。

（3）水浴溫度是SCV運(yùn)行過(guò)程中一個(gè)重要參數(shù)，需要將其控制在合理范圍內(nèi)。模擬結(jié)果顯示，水浴溫度隨初始水位高度和進(jìn)口LNG速度增加而降低，隨煙氣量增加而升高。在一定參數(shù)范圍內(nèi)，管程入口壓力對(duì)水浴溫度影響較小。優(yōu)化后的煙氣量和水位高度可以使得NG出口溫度和水浴溫度滿(mǎn)足運(yùn)行要求。

符 號(hào) 說(shuō) 明

cp——比定壓熱容，kJ·kg-1·K-1 H ——初始水位高度，mm K——傳熱系數(shù)，W·m-2·K-1 L——管長(zhǎng)，m p ——壓力，MPa Q ——質(zhì)量流量，kg·h-1 q ——壁面熱通量，kW·m-2 s——熵，J·g-1·K-1 T ——溫度，K V ——體積流量，m3·h-1 v ——速度，m·s-1 l——熱導(dǎo)率，W·m-1·K-1 m——黏度，Pa·s r——密度，kg·m-3 下角標(biāo) a——空氣 b——水浴 cr——臨界 f ——燃料 fg ——煙氣 in ——進(jìn)口 L——局部 out ——出口 pc——擬臨界 s ——?dú)こ?t ——管程 wi ——內(nèi)壁面 wo ——外壁面

References

[1] 梅鵬程, 鄧春鋒, 鄧欣. LNG氣化器的分類(lèi)及選型設(shè)計(jì)[J]. 化學(xué)工程與裝備, 2016, (5): 65-70. MEI P C, DENG C F, DENG X. Classification and selection design of LNG vaporizer[J]. Chemical Engineering & Equipment, 2016, (5): 65-70.

[2] 康風(fēng)立, 孫海峰, 熊亞選, 等. 浸沒(méi)燃燒式LNG氣化器水浴氣化傳熱計(jì)算[J]. 油氣儲(chǔ)運(yùn), 2016, 35(4): 406-411. KANG F L, SUN H F, XIONG Y X,. Computation of heat transfer in water bath of LNG SCV[J]. Oil & Gas Storage and Transpotation, 2016, 35(4): 406-411.

[3] HAN D Y, XU Q Q, ZHOU D, et al. Design of heat transfer in submerged combustion vaporizer[J]. Journal of Natural Gas Science and Engineering,2016, 31: 76-85.

[4] MENG H Y, WANG S Z, ZHOU L, et al. Influence on inlet mass flow rate on heat transfer of supercritical liquefied natural gas in horizontal tubes[J]. Advanced Materials Research, 2014, 960: 433-437.

[5] PARK Y C, Kim J. Submerged combustion vaporizer optimization using entropy minimization method[J]. Applied Thermal Engineering, 2016, 103: 1071-1076.

[6] 齊超, 王博杰, 易沖沖, 等. 浸沒(méi)燃燒式氣化器的運(yùn)行特性及優(yōu)化[J]. 化工學(xué)報(bào), 2015, 66(S2): 199-205. QI C, WANG B J, YI C C, et al. Operation characteristics and optimization of submerged combustion vaporizer[J]. CIESC Journal, 2015, 66(S2): 199-205.

[7] 李仲珍, 郭少龍, 陶文銓. 超臨界LNG管內(nèi)流動(dòng)與換熱特性研究[J]. 工程熱物理學(xué)報(bào), 2013, 34(12): 2314-2317. LI Z Z, GUO S L, TAO W Q. Studies of supercritical convective heat transfer of LNG in tube[J]. Journal of Engineering Thermophysics, 2013, 34(12): 2314-2317.

[8] 張康, 韓昌亮, 任婧杰, 等. SCV蛇形換熱管內(nèi)超臨界LNG傳熱特性數(shù)值模擬[J]. 化工學(xué)報(bào), 2015, 66(12): 4788-4795. ZHANG K, HAN C L, REN J J, et al. Numerical simulation on heat transfer of supercritical LNG in coil tubes of submerged combustion vaporizer[J]. CIESC Journal, 2015, 66(12): 4788-4795.

[9] 靳書(shū)武, 武錦濤, 銀建中. 水平圓管內(nèi)超臨界甲烷對(duì)流換熱數(shù)值模擬[J]. 應(yīng)用科技, 2015, 42(5): 67-71. JIN S W, WU J T, YIN J Z. Numerical simulation of convective heat transfer of supercritical methane in a horizontal tube[J]. Applied Science and Technology, 2015, 42(5): 67-71.

[10] 李泓鈺. 浸沒(méi)燃燒式LNG氣化器模擬實(shí)驗(yàn)研究[D]. 北京: 北京建筑大學(xué), 2015. LI H Y. The study of basic characteristics of submerged combustion vaporizer[D]. Beijing: Beijing Architecture University, 2015.

[11] 粘權(quán)鑫, 郭少龍, 方文振, 等. 液化天然氣浸沒(méi)燃燒式氣化器數(shù)值模擬方法研究[J]. 西安交通大學(xué)學(xué)報(bào), 2016, 50(1): 68-72. NIAN Q X, GUO S L, FANG W Z, et al. Numerical simulation of liquefied natural gas submerged combustion vaporizer[J]. Journal of Xi’an Jiaot ong University, 2016, 50(1): 68-72.

[12] HAN C L, REN J J, WANG Y Q, et al. Experimental studies of shell-side fluid flow and heat transfer characteristics in a submerged combustion vaporizer[J].International Journal of Heat and Mass Transfer, 2016, 101: 436-444.

[13] RIBEIRO C P, LAGE P L C. Direct-contact evaporation in the homogeneous and heterogeneous bubbling regimes(Ⅰ): Experimental analysis[J]. International Journal of Heat and Mass Transfer, 2004, 47(17/18): 3825-3840.

[14] MCLINDEN M O, KLEIN S A, LEMMON E W,. NIST standard reference database 23: NIST thermodynamic and transport properties of refrigerants and refrigerant mixtures—REFPROP, Version 6.0 [DB/CD]. USA: National Institute of Standards and Technology, 1998.

[15] HAN C L, REN J J, DONG W P, et al. Numerical investigation of supercritical LNG convective heat transfer in a horizontal serpentine tube[J]. Cryogenics, 2016, 78: 1-13.

[16] 韓昌亮, 任婧杰, 董文平, 等. 沉浸式汽化器殼程流體傳熱實(shí)驗(yàn)與數(shù)值模擬[J]. 化工學(xué)報(bào), 2016, 67(10): 4095-4103. HAN C L, REN J J, DONG W P,. Experimental and numerical simulation on shell side fluid heat transfer in submerged combustion vaporizer[J].CIESC Journal, 2016. 67(10): 4095-4103.

[17] YU S Q, LI H X, LEI X L,. Influence of buoyance on heat transfer to water flowing in a horizontal tubes under supercritical pressure[J]. Applied Thermal Engineering, 2013, 59:380-388.

[18] CHENG X, KUANG B, YANG Y H. Numerical analysis of heat transfer in supercritical water cooled flow channels[J]. Nuclear Engineering and Design, 2007, 237(3): 240-252.

[19] 王珂, 謝金, 劉遵超, 等. 超臨界二氧化碳在微細(xì)管內(nèi)的換熱特性[J]. 化工學(xué)報(bào), 2014, 65(S1):323-327. WANG K, XIE J, LIU Z C, et al. Heat transfer characteristics of supercritical carbon dioxide in a micro-capillary tube[J].CIESC Journal, 2014, 65(S1): 323-327.

[20] 王樹(shù)立. 鼓泡塔氣液兩相流的數(shù)值模擬與實(shí)驗(yàn)研究[D]. 大連: 大連理工大學(xué), 2002. WANG S L. Numerical and experimental investigation on gas-liquid two-phase flow in a bubble column[D]. Dalian: Dalian University of Technology, 2002.

Experimental study and numerical simulation on coupled heat transfer characteristics of submerged combustion vaporizer

HAN ChangliangREN JingjieWANG YanqingDONG WenpingBI Mingshu

(School of Chemical Machinery, Dalian University of Technology, Dalian 116024, Liaoning, China)

Submerged combustion vaporizer is indispensable heat exchanger equipment for liquefied natural gas (LNG) receiving terminals, it mainly utilizes water bath system as intermediate media to achieve heat transfer between flue gas and LNG. In this work, an unabridged experimental apparatus is built to study the complex heat transfer characteristics of it. Experimental results reveal some unique fluid dynamic phenomena inside SCV system (local water bath freezes). Meanwhile, the simulated model is established for coupled flow and heat transfer process between two-phase mixture and trans-critical LNG. The influences of inlet LNG pressure, inlet LNG velocity, static water height, flue gas flux on the NG outlet temperature and water bath temperature are analyzed. The outcomes can provide some important guidance to localization design of SCV.

submerged combustion vaporizer; LNG; experiment; flow; heat transfer; water bath temperature; numerical simulation

10.11949/j.issn.0438-1157.20160973

TE 088

A

0438—1157（2017）03—0854—10

中央高?；究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金項(xiàng)目（DUT16QY29）。

2016-07-11收到初稿，2017-01-27收到修改稿。

聯(lián)系人：畢明樹(shù)。第一作者：韓昌亮（1987—），男，博士研究生。

2016-07-11.

Prof. BI Mingshu, bimsh@dlut.edu.cn

supported by the Fundamental Research Funds for the Central Universities (DUT16QY29).