張 猛，耿堯辰（中國天辰工程有限公司，天津 300400）

圓塊孔式石墨換熱器總傳熱系數(shù)的模擬計算

張 猛，耿堯辰
（中國天辰工程有限公司，天津 300400）

介紹了一種利用 HTRI軟件中的套管式換熱器模型對圓塊孔式石墨換熱器進行計算的模擬方法，并將模擬結(jié)果與廠家計算值進行比較，大部分結(jié)果偏差在 10%以內(nèi)，表明此種計算方法具有較好的工程指導(dǎo)意義。

圓塊孔式石墨換熱器；總傳熱系數(shù)；HTRI

1 概述

不透性石墨制換熱設(shè)備具有優(yōu)良的耐腐蝕性和傳熱性能，被廣泛應(yīng)用于氯堿行業(yè)、磷酸濃縮以及鋼鐵酸洗等涉及腐蝕性介質(zhì)傳熱的工藝過程中。石墨換熱器按其結(jié)構(gòu)形式主要分為列管式和圓塊孔式兩種，其中圓塊孔式換熱器因其具有密封性能良好、使用壓力高以及結(jié)構(gòu)緊湊等優(yōu)點常被應(yīng)用于操作條件苛刻或有熱沖擊及震動的場合。

圓塊孔式石墨換熱器由若干圓形截面的浸漬石墨塊組成，石墨塊之間用O型圈密封。 在石墨塊上有許多平行于軸線和垂直于軸線的圓孔流道，前者稱為縱向流道，后者稱為橫向流道，其中，橫向流道位于各縱向流道的間隔中間。一種介質(zhì)（一般為工藝介質(zhì)）沿軸向流道從上往下或從下往上流動；而另一種介質(zhì)（一般為服務(wù)介質(zhì)）沿橫向流道流動，相鄰兩塊換熱塊橫向流動方向相逆。兩種介質(zhì)利用換熱塊石墨實體進行 間壁式換熱[1]。

盡管圓塊孔式石墨換熱器已得到廣泛的應(yīng)用，其傳熱計算方法卻鮮有報道，通用的換熱器計算軟件也沒有其計算模型。南京工業(yè)大學(xué)的潘虎自上世紀80年代開始對圓塊孔式石墨換熱器的傳熱計算進行研究，并結(jié)合實驗數(shù)據(jù)提出了縱向、橫向流道側(cè)傳熱膜系數(shù)以及總傳熱系數(shù)的計算公式[2-3]。 然而，此公式僅適用于簡單的非相變傳熱計算，并不能滿足工藝設(shè)計中多組分冷凝、蒸發(fā)等復(fù)雜的相變傳熱計算。

目前，圓塊孔式石墨換熱器的傳熱計算被掌握在少數(shù)石墨設(shè)備制造廠家手中，設(shè)計單位和用戶無法有效的校核其換熱器設(shè)計的合理性，因此在換熱器的設(shè)計和訂貨過程中處于不利位置。本文介紹了一種利用 HTRI軟件對多組分相變傳熱過程進行計算的方法，并將其應(yīng)用到中國天辰公司承接項目中圓塊孔式石墨換熱器的校核計算，大部分計算的總傳熱系數(shù) K 縱（以縱向孔道表面積為基準）與廠家計算值的偏差小于 10%，提供了一種較可靠的圓塊孔式石墨換熱器的校核計算方法。

2 總傳熱系數(shù)K縱的計算公式

圓塊孔式石墨換熱器的傳熱機理及過程與列管式的相同，都 包括以下 3 個 過程[4]。

（1）熱介質(zhì)將熱量傳給高溫壁面；

（2）熱量由高溫壁面通過石墨間壁傳給低溫壁面；

（3）熱量由低溫壁面?zhèn)鹘o冷介質(zhì)。

則K縱應(yīng)由3部分組成：

式中：K縱—總傳熱系 數(shù)，W/（m2·℃）；

α縱—工藝側(cè) 傳熱膜系數(shù)，W/（m2·℃）；

α橫—服務(wù)側(cè) 傳熱膜系數(shù)，W/（m2·℃）；

A縱—縱向孔道總面積，m2；

A橫—橫向孔道總面積，m2；

A均—平均傳熱表面積，m2，A

k—石墨導(dǎo)熱系數(shù)，W/（m·℃）；

V′—石墨塊鉆去孔道后的剩余體積，m3。

考慮工藝側(cè)及服務(wù)側(cè)的污垢熱阻R縱和R橫，總污垢熱阻為則

其中，A縱、A橫、A均以及b′由石墨塊的結(jié)構(gòu)決定，R橫和R縱由工藝條件或者工程經(jīng)驗給出，k 根據(jù)文獻取 100 W/（m·℃）。 根據(jù)圓 塊孔式石墨換 熱 器的結(jié)構(gòu)，可將工藝介質(zhì)和服務(wù)介質(zhì)分別按照在垂直和水平圓管內(nèi)作強制對流或發(fā)生相變傳熱來考慮，通過HTRI 軟件模擬計算或通過 Dittus－Boelter 關(guān)聯(lián)式等半經(jīng)驗公式來求得α縱和α橫。

3 計算舉例

以項目中某汽提塔冷凝器為例進行校核計算，換熱器的結(jié)構(gòu)性質(zhì)見表1。

換熱器的流股性質(zhì)見表2。

工藝介質(zhì)的質(zhì)量組成見表3。

石墨塊的結(jié)構(gòu)參數(shù)見圖1。

表1 換熱器的結(jié)構(gòu)性質(zhì)

表2 換熱器的流股性質(zhì)

表3 工藝介質(zhì)的質(zhì)量組成

圖1 石墨塊的結(jié)構(gòu)參數(shù)

3.1 服務(wù)側(cè)傳熱膜系數(shù)α橫的計算

此例中服務(wù)側(cè)為水的無相變強制對流傳熱，α橫可 由 Dittus－Boelter 關(guān) 聯(lián) 式 ， 即計算得出。

3.2 工藝側(cè)傳熱膜系數(shù) α縱的計算

此例中工藝側(cè)為多組分蒸汽冷凝的復(fù)雜相變傳熱，利用HTRI軟件中的套管式換熱器模型（Jacketed Pipe Exchanger）來模擬此傳熱過程，設(shè)定工藝熱物流走管內(nèi)（Hot Pipe），計算結(jié)果中的管內(nèi)傳熱膜系數(shù)（Pipe h）即為所求α縱。 下面簡要介紹一下程序中各參數(shù)的設(shè)置。

在 Input Summary選項中 ， 將 Case Mode 選為Rating。 在 Jacket Geometry 中輸入夾套管的參數(shù)，其中夾套管內(nèi)徑選為 55 mm （后文會討論夾套管內(nèi)徑對模擬結(jié)果的影響），熱流股走管內(nèi)，安裝方向為垂直。 在 Process Conditions 中輸入冷熱物流的流股性質(zhì)，熱物流的流量為 13.58 kg/h，冷凝水由于溫升已確定故不必輸入流量。 在 Pipe Geometry 中輸入垂直圓管的參數(shù)，其中管長 L=0.34×8=2.72 m （忽略各換熱塊間的間距），壁厚選取換熱塊的當(dāng)量壁厚 7mm，管外徑 OD=10+7×2=24（mm），管子類型為光滑管。

在 Geometry 選項中 ，僅將 Pipe Material中的管子材質(zhì)選為石墨， 導(dǎo)熱系數(shù)設(shè)為 100 W/（m·℃），其余選項設(shè)為默認值。 設(shè)置完成后在 Nozzle Location選項中確認模擬換熱器的安裝方式及熱物流流動方向是否與石墨換熱器相同，見圖 2。

圖2 模擬換熱器的安裝方式及物流流動方向

在 Process 選項中，定義冷熱物流的名稱并輸入兩側(cè)的污垢熱阻，根據(jù)工程經(jīng)驗水側(cè)及工藝側(cè)的污垢熱阻分別設(shè)為0.000 396 和 0.000 198 （m2·℃ ）/W。至此，換熱器結(jié)構(gòu)性質(zhì)參數(shù)已輸入完畢，下面進行流股性質(zhì)及組成的輸入。

在 Hot Fluid Properties 選項中 ， 打開 Property Generator，將 Property package 選為 VMG Thermo，氣液相模型選為默認的 Advanced Peng-Robinson。 在Composition選項中，將 Composition Basis 選 為 Mass后按照表 2.3 輸入工藝介質(zhì)的組成及質(zhì)量分數(shù)。 之后 在 Conditions 選 項 中 輸 入 壓 力 500 kPa 及 溫 度 范圍 42 ℃~154 ℃，生成流股物性并將結(jié)果全部導(dǎo)入。

在 Cold Fluid Properties 選項中，將 Physical Property Input Option 一 項設(shè)定為 Component by Component。 在 Components 選項中，選擇 Package 為VMG Thermo，在數(shù)據(jù)庫中選擇組分為水。

至此所有參數(shù)輸入完畢，運行模擬。報告中Pipeh 一項即為所求的 α縱，結(jié)果為 260.7W/（m2·℃）。

3.3 夾套管內(nèi)徑對 α縱模擬結(jié)果的影響

一般認為一側(cè)的傳熱膜系數(shù)只與該側(cè)流體的流動情況有關(guān)，因為冷卻水流量已確定，所以改變夾套管的內(nèi)徑會改變冷卻水的流速進而改變環(huán)側(cè)的傳熱膜系數(shù)，但卻不會對管內(nèi)傳熱膜系數(shù)的大小造成影響。為了驗證夾套管內(nèi)徑對管內(nèi)傳熱膜系數(shù)的影響，本文進行了靈敏度分析，結(jié)果見圖 3。

圖3 靈敏度分析結(jié)果

由靈敏度分析的結(jié)果可以看出，夾套管內(nèi)徑的取值對于環(huán)側(cè)傳熱膜系數(shù)影響很大，而對管內(nèi)傳熱膜系數(shù)基本沒有影響。因此，在模擬過程中夾套管內(nèi)徑的取值只需大于圓管外徑即可，在本例中將夾套管內(nèi)徑設(shè)為 55 mm。

3.4 總傳熱系數(shù) K縱的計算

考慮污垢熱阻時：

不考慮污垢熱阻時：

模擬結(jié)果與廠商計算值的比較見表4所示。

表4 模擬結(jié)果與廠商計算值的比較

所有樣本的偏差分析見表5。

表5 樣本的偏差分析

由表5可以看出偏差多為負偏差，即模擬的總傳熱系數(shù)K值小于廠家計算值。在設(shè)置相關(guān)參數(shù)進行模擬時， 管長為石墨塊的組合高度，HTRI軟件將縱向孔道認為是一根連續(xù)的圓管進行傳熱計算。而真實情況是相鄰塊體之間留有 2~4 mm 的間隙 （組裝O型密封圈產(chǎn)生）， 此空間相當(dāng)于縱向流道中介質(zhì)再分布的湍流室，可提高縱向孔道的傳熱膜系數(shù)。因此模擬計算值普遍小于廠商計算值，可將模擬計算值乘以合適的系數(shù)予以校正。

4 總結(jié)

根據(jù)圓塊孔式石墨換熱器的結(jié)構(gòu)， 借助 HTRI軟件中套管式換熱器的模型，本文提出了一種針對圓塊孔式石墨換熱器的校核計算方法，可將其應(yīng)用于圓塊孔式石墨換熱器校核計算中，改變設(shè)計單位和用戶在換熱器的設(shè)計與訂貨過程中的不利位置，具有較好的工程指導(dǎo)意義。

［1］許志遠，等．石墨制化工設(shè)備，北京：化學(xué)工業(yè)出版社，2003.

［2］潘 虎．圓塊式石墨換熱器傳熱系數(shù)的計算，化工設(shè)計通訊，1988，14（2）：60～63.

［3］潘 虎．圓塊式石墨換熱器的傳熱系數(shù)，石油化工設(shè)備，1978，16（11）：18～22.

［4］錢頌文．換熱器設(shè)計手冊，北京：化學(xué)工業(yè)出版社，2006.

Total heat transfer coefficient simulation calculation of graphite block heat exchanger

ZHANG Meng， GENG Yao-chen
（China Tianchen Engineering Corporation,Tianjin 300400,China）

On the basis of applying double-pipe heat exchanger model from HTRI，a heat-transfer calculation method for graphite block heat exchanger is studied in the paper.Most biases between calculation results and manufacture parameters are less than 10%，indicating the method is applicable to engineering design.

graphite block heat exchanger；total heat transfer coefficient；HTRI

TQ051.5

：B

：1009－1785(2017)02-0030-04

2017-01-27