劉俊杰，王國清，張利軍，周先鋒

（中國石化 北京化工研究院，北京 100013）

含有扭曲片的裂解爐線性廢熱鍋爐的數(shù)值模擬和工藝計算

劉俊杰，王國清，張利軍，周先鋒

（中國石化 北京化工研究院，北京 100013）

利用數(shù)值模擬技術(shù)對紐帶和扭曲片強化傳熱管進行了模擬，將模擬結(jié)果應(yīng)用于商業(yè)傳熱軟件中，對含有扭曲片的裂解爐線性廢熱鍋爐進行了工藝計算。計算了兩種強化傳熱管的速度場、壓降場和溫度場，給出了不同雷諾數(shù)下的壓降、努賽爾數(shù)、強化傳熱因子、壓降比和努塞爾數(shù)比。數(shù)值模擬結(jié)果表明，扭曲片強化傳熱管的強化傳熱性能優(yōu)于紐帶強化傳熱管，前者的壓降和努塞爾數(shù)均小于后者，前者的強化傳熱因子大于后者。工藝計算結(jié)果表明，裂解原料分別為加氫尾油、石腦油及乙烷時，具有扭曲片強化傳熱管的線性廢熱鍋爐的傳熱量比普通線性廢熱鍋爐的傳熱量增加了8.2%，6.1%，2.4%。

紐帶；扭曲片；強化傳熱；裂解爐；線性廢熱鍋爐；數(shù)值模擬；乙烯裝置

管內(nèi)插入物是強化管內(nèi)單相流體傳熱的一種有效方法，紐帶和扭曲片是兩種比較重要的插入物結(jié)構(gòu)形式。Masoud等［1］對各種紐帶進行了研究，提出了一種帶有缺口的紐帶強化傳熱元件，獲得了較高的強化傳熱因子。Sarada等［2］以空氣為介質(zhì)對不同寬度的紐帶的強化傳熱效果進行了研究，擬合了努塞爾數(shù)和弗魯?shù)聰?shù)的經(jīng)驗公式。洪蒙納等［3-5］對扭曲片強化傳熱機理也進行了實驗和數(shù)值模擬研究，給出了合理的理論解釋。Smith等［6-7］對不同間距的扭曲片的強化傳熱效果進行了研究，擬合了經(jīng)驗公式。王國清等［8-14］對扭曲片進行了數(shù)值模擬研究，優(yōu)化了扭曲片的結(jié)構(gòu)，在裂解爐輻射段爐管中取得了很好的工業(yè)應(yīng)用效果。上述研究過程大多是在小管徑換熱管中進行，大管徑高流速換熱管的數(shù)據(jù)并不是很多。同時，如何將扭曲片強化傳熱管推廣應(yīng)用到石化領(lǐng)域的換熱器中也是一個值得研究的問題。線性廢熱鍋爐（廢鍋）是高氣速大直徑換熱器，因此線性廢鍋是扭曲片強化傳熱管一個很好的應(yīng)用場合。但有關(guān)具有扭曲片強化傳熱管的線性廢鍋的流動和傳熱的基礎(chǔ)數(shù)據(jù)比較匱乏。同時，由于商業(yè)傳熱軟件缺少扭曲片的壓降和傳熱模型，因此如何利用商業(yè)傳熱軟件進行含有扭曲片的線性廢鍋的工業(yè)計算也是一個亟待解決的問題。

本工作通過對線性廢鍋中紐帶和扭曲片強化傳熱管兩種不同強化傳熱構(gòu)件進行數(shù)值模擬，并將模擬結(jié)果應(yīng)用于商業(yè)傳熱軟件中，對具有扭曲片的線性廢鍋進行工藝計算。

1 含有強化傳熱元件的線性廢鍋的數(shù)值模型的建立

1.1 兩種強化傳熱管的實體建模

線性廢鍋的內(nèi)徑為81.6 mm，模擬區(qū)域為扭曲片管前端5倍直徑以內(nèi)和扭曲片管后端50倍直徑以內(nèi)的區(qū)域，扭曲比為2.5，利用GAMBIT軟件實現(xiàn)具有強化傳熱元件的線性廢鍋的實體建模。由于扭曲片強化傳熱管的結(jié)構(gòu)復(fù)雜，無法采用結(jié)構(gòu)化網(wǎng)格，利用四面體的非結(jié)構(gòu)網(wǎng)格進行劃分，線性廢鍋的其余部分采用六面體網(wǎng)格。紐帶貫穿于整個線性廢鍋，整體采用四面體網(wǎng)格。扭曲片強化傳熱管的實體建模見圖1，網(wǎng)格數(shù)為80萬左右。紐帶強化傳熱管的實體建模見圖2，網(wǎng)格數(shù)為110萬左右。

圖1 用于強化傳熱管的扭曲片的結(jié)構(gòu)Fig.1 Structure of a twisted tape used in strengthened heat transfer tube(twisted tape tube).

圖2 用于強化傳熱管的紐帶的結(jié)構(gòu)Fig.2 Structure of a bond used in strengthened heat transfer tube(bond tube).

1.2 數(shù)值計算模型

1.2.1 數(shù)學(xué)模型

通過比較各種數(shù)學(xué)模型的適用范圍和優(yōu)缺點以及分析構(gòu)件的具體結(jié)構(gòu)，在計算扭曲片強化傳熱管時，采用realize k-ε模型對流體流動過程進行模擬。計算中選取的連續(xù)性方程、動量方程、能量方程和湍動能方程如下：

連續(xù)性方程：

動量方程：

能量方程：

上述各方程中，μt=μm+μturd；λt=λm+λturd；Dt=Dm+Dturd。

湍動能方程：

其中，k方程為：

ε方程為：

1.2.2 求解方法和邊界條件

在數(shù)值模擬過程中，采用隱式求解器，對各偏微分方程采用二階迎風(fēng)格式進行離散，利用壓力速度耦合的SIMPLIC算法。線性廢鍋的進口速度為60，80，100，120 m/s，管壁溫度恒定為326 ℃，管出口的邊界條件為壓力出口。

1.3 強化傳熱元件的評價函數(shù)

通常利用強化傳熱因子來評價強化傳熱元件的優(yōu)劣，強化傳熱因子（nfactor）為：

2 結(jié)果與討論

2.1 兩種強化傳熱管的速度場分析

扭曲片和紐帶強化傳熱管的流線軌跡、軸向截面速度場、橫向截面速度場以及切向速度沿管長的分布分別見圖3～6，以進口速度80 m/s為例進行闡述，溫度場等同。從圖3～6可看出，扭曲片強化傳熱管內(nèi)的流體由活塞流逐漸變成旋轉(zhuǎn)流，且有直徑50倍長度的影響區(qū)域。管內(nèi)速度場在接近扭曲片時產(chǎn)生切向速度，經(jīng)過扭曲片后，切向速度逐漸減弱，從16 m/s逐漸衰減為零。紐帶強化傳熱管內(nèi)的流動跡線顯示了速度變化規(guī)律，速度呈周期性變化，在10 m/s到25 m/s區(qū)域，沿管長周期性減小。管中心區(qū)域速度值相對較大，而管壁區(qū)域速度值相對較小。兩種強化傳熱管的截面速度分布規(guī)律接近，且紐帶強化傳熱管的平均切向速度大于扭曲片強化傳熱管的平均切向速度。

圖3 扭曲片（a）和紐帶（b）強化傳熱管的流線軌跡Fig.3 Pathes of filament lines of the twisted tape tube(a) and bond tube(b).

圖4 扭曲片（a）和紐帶（b）強化傳熱管的軸向截面速度場Fig.4 Axial section velocity fields of the twisted tape tube(a) and bond tube(b).

圖5 扭曲片（a）和紐帶（b）強化傳熱管的橫向截面速度場Fig.5 Transverse section velocity fields of the twisted tape tube(a) and bond tube(b).

圖6 扭曲片（a）和紐帶（b）強化傳熱管的切向速度沿管長的分布規(guī)律Fig.6 Distributions of the tangential velocities of the twisted tape tube(a) and bond tube(b) along axial direction.

兩種強化傳熱管截面內(nèi)的湍動強度見圖7。由圖7可見，扭曲片強化換熱管壁面和中心區(qū)域的湍動強度分布的不均勻性較大，在扭曲片區(qū)域及其端部的湍動強度更大；紐帶強化傳熱管的湍動強度在中心區(qū)域相對較小，而在紐帶區(qū)域略大，分布較均勻。

圖7 扭曲片（a）和紐帶（b）強化傳熱管截面內(nèi)的湍動強度Fig.7 Turbulence intensity distribution of the twisted tape tube(a) and bond tube(b).

2.2 兩種強化傳熱管的壓降場分析

兩種強化傳熱管的總壓降沿管長的分布見圖8。由圖8可見，紐帶強化傳熱管的壓降大于扭曲片強化傳熱管。扭曲片強化傳熱管在扭曲片區(qū)域的壓降迅速下降，在扭曲片管后部的壓降逐漸下降，兩者下降斜率不同，前部較大，后部較??；盡管紐帶強化傳熱管的壓降沿管長有一些周期性變化，但整體趨勢是逐漸下降的。

圖8 扭曲片（a）和紐帶（b）強化傳熱管的總壓降隨管長的變化規(guī)律Fig.8 Changes of the total pressure drops with the lengths of the twisted tape tube(a) and bond tube(b) along axial direction.

2.3 兩種強化傳熱管的溫度場分析

兩種強化傳熱管的軸向截面和橫向截面溫度場的分布分別見圖9～10。由圖9可見，兩種強化傳熱管的溫度場比較接近，管內(nèi)溫度沿軸向逐漸降低，在管出口位置溫度達到最低。由于在裂解爐線性廢鍋中發(fā)生的是裂解氣的冷卻過程，管壁附近溫度較低，而管內(nèi)區(qū)域的溫度相對較高。從圖10可看出，扭曲片強化傳熱管在管中心區(qū)域的溫度相對較高，僅扭曲片兩端部分的溫度相對較低，呈對稱分布；紐帶強化傳熱管橫向截面溫度場的分布與扭曲片強化傳熱管的相似，但溫度分布更加均勻。

圖9 扭曲片（a）和紐帶（b）強化傳熱管軸向截面的溫度場分布Fig.9 Distributions of temperature fields in the twisted tape tube(a) and bond tube(b) in axial section.

圖10 扭曲片（a）和紐帶（b）強化傳熱管橫向截面的溫度場分布Fig.10 Distributions of temperature fields in the twisted tape tube（a） and bond tube（b） in transverse section.

2.4 兩種強化傳熱管的評價

扭曲片和紐帶強化傳熱管的壓降及壓降比見圖11。從圖11可看出，隨雷諾數(shù)的增加，兩種強化傳熱管的壓降逐漸增加，紐帶強化傳熱管的壓降遠大于扭曲片強化傳熱管的壓降；在相同的雷諾數(shù)下，扭曲片強化傳熱管的壓降是紐帶強化傳熱管壓降的0.35～0.40倍。

扭曲片和紐帶強化傳熱管的努塞爾數(shù)與努塞爾數(shù)比見圖12。從圖12可看出，兩種強化傳熱管的努塞爾數(shù)均隨雷諾數(shù)的增大逐漸增加，且數(shù)值相差很小。在相同雷諾數(shù)下，扭曲片強化傳熱管的努塞爾數(shù)是紐帶強化傳熱管的0.95～0.98倍。

圖11 扭曲片和紐帶強化傳熱管的壓降（a）及壓降比（b）Fig.11 Pressure drops in the twisted tape tube and bond tube(a)，and the pressure drop ratio(b).

圖12 扭曲片和紐帶強化傳熱管的努塞爾數(shù)（a）與努塞爾數(shù)比（b）Fig.12 Nu of the twisted tape tube and bond tube(a) and the Nu ratio(b).

扭曲片和紐帶強化傳熱管的強化傳熱因子與強化傳熱因子的比見圖13。從圖13可看出，兩種強化傳熱管的強化傳熱因子均隨雷諾數(shù)的增加而略有減小，且扭曲片強化傳熱管的強化傳熱因子大于紐帶強化傳熱管的強化傳熱因子。這說明扭曲片強化傳熱管的強化傳熱性能優(yōu)于紐帶強化傳熱管，即扭曲片強化傳熱管在較低壓降下能達到較好的傳熱效果。

圖13 扭曲片和紐帶強化傳熱管的強化傳熱因子（a）與強化傳熱因子的比（b）Fig.13 Strengthened heat transfer factors of the twisted tape tube and bond tube(a) and their ratio(b).

3 具有扭曲片強化傳熱管的裂解爐線性廢鍋的工藝計算

根據(jù)兩種強化傳熱管數(shù)值模擬結(jié)果得到的壓降和努塞爾數(shù)比值的關(guān)系，利用商業(yè)傳熱軟件中紐帶強化傳熱管模型代替扭曲片強化傳熱管模型，對裂解原料為加氫尾油、石腦油和乙烷的150 kt/a具有扭曲片強化傳熱管的裂解爐線性廢鍋進行工藝計算。從商業(yè)流程模擬軟件中將裂解產(chǎn)物物性數(shù)據(jù)輸入到商業(yè)傳熱軟件，其中，加氫尾油和石腦油的物性參數(shù)通過虛擬組分的方法進行計算；之后，輸入線性廢鍋的結(jié)構(gòu)參數(shù)和操作參數(shù)，利用商業(yè)軟件中的流動和傳熱模型完成工藝計算。分別以加氫尾油、石腦油和乙烷為裂解原料時，線性廢鍋的運行周期分別為50，60，80 d；在線性廢鍋中加入3個扭曲片強化傳熱管后，線性廢鍋的運行周期延長一倍。加入扭曲片前后在同一周期內(nèi)線性廢鍋的壓降、出口溫度及傳熱量的變化見表1～6。以加氫尾油為例，從表1和表2可看出，在運行初期，加入扭曲片強化傳熱管前后，出口溫度相差5 ℃左右；而在運行中期，加入扭曲片強化傳熱管前后，出口溫度則相差30 ℃左右；而在運行后期，加入扭曲片強化傳熱管前后，出口溫度相差35 ℃。這主要是由于在運行初期，幾乎沒有焦層，扭曲片主要強化管程；在運行中期，焦層厚度開始增加，焦層和管程的熱阻相當(dāng)，扭曲片強化傳熱和抑制結(jié)焦的作用得到充分發(fā)揮；而在運行后期，熱阻主要集中在焦層上，扭曲片主要的作用是抑制結(jié)焦。從表1～6還可看出，在以加氫尾油、石腦油和乙烷為裂解原料時，加入扭曲片強化傳熱管后線性廢鍋的傳熱量分別增加了8.2%，6.1%，2.4%。這是因為原料不同，在線性廢鍋中結(jié)焦的厚度和焦層熱阻也有所不同。

表1 線性廢鍋的壓降、出口溫度及傳熱量（加氫尾油為原料）Table 1 Pressure drop，outlet temperature and heat flux of a linear waste heat boiler(hydrogenated tail oil as feedstock)

表2 加扭曲片強化傳熱管后線性廢鍋的壓降、出口溫度及傳熱量（加氫尾油為原料）Table 2 Pressure drop， outlet temperature and heat flux of linear waste heat boiler with the twisted tape tube(hydrogenated tail oil as feedstock)

表3 線性廢鍋的壓降、出口溫度及傳熱量（石腦油為原料）Table 3 Pressure drop， outlet temperature and heat flux of linear waste heat boiler(naphtha as feedstock)

表4 加扭曲片強化傳熱管后裂解線性廢鍋的壓降、出口溫度及傳熱量（石腦油為原料）Table 4 Pressure drop，outlet temperature and heat flux of linear waste heat boiler with the twisted tape tube(naphtha as feedstock)

表5 線性廢鍋的壓降、出口溫度及傳熱量（乙烷為原料）Table 5 Pressure drop，outlet temperature and heat flux of linear waste heat boiler(ethane as feedstock)

表6 加扭曲片強化傳熱管后裂解線性廢鍋的壓降、出口溫度及傳熱量（乙烷為原料）Table 6 Pressure drop，outlet temperature and heat flux of linear waste heat boiler with the twisted tape tube(ethane as feedstock)

4 結(jié)論

1）扭曲片強化傳熱管的強化傳熱性能優(yōu)于紐帶強化傳熱管。扭曲片強化傳熱管的壓降和努塞爾數(shù)均小于紐帶強化傳熱管，扭曲片強化傳熱管的壓降為紐帶強化傳熱管的0.35～0.40倍，努塞爾數(shù)為紐帶強化管的0.95～0.98倍。扭曲片強化傳熱管的強化傳熱因子大于紐帶強化傳熱管。

2）將扭曲片強化傳熱管的數(shù)值模擬結(jié)果應(yīng)用于裂解爐線性廢鍋的工藝計算中，裂解原料越輕，傳熱量增加越小，裂解原料為分別加氫尾油、石腦油和乙烷時，線性廢鍋的傳熱量分別增加了8.2%，6.1%，2.4%。

符 號 說 明

c1k-ε模型方程參數(shù)

c2k-ε模型方程參數(shù)

Dm分子擴散系數(shù)，m2/s

Dt有效擴散系數(shù)，m2/s

Dturb湍動擴散系數(shù)，m2/s

ek動量方程中的模型參數(shù)

f 強化傳熱管的摩擦系數(shù)

f0光管的摩擦系數(shù)

H 由對流引起的單位焓，J/kg

Hj物質(zhì)j的焓，J/kg

K 由熱傳導(dǎo)引起的單位焓，J/kg

k 單位質(zhì)量的湍動能，m2/s2

Nu 強化傳熱管的努塞爾數(shù)

Nu0光管的努塞爾數(shù)

nfactor強化傳熱因子

p 壓力，Pa

Rε附加生成項， 代表平均應(yīng)變率對ε 的影響

T 溫度，K

t 單位時間，s

ui在i方向上的速度分量，m/s

uj在j方向上的速度分量，m/s

uk在k方向上的速度分量，m/s

xi在i方向上的長度，m

xj在j方向上的長度，m

xk在k方向上的長度，m

yj各組分的摩爾分數(shù)

ε 單位質(zhì)量的湍動能耗散率，m2/s3

η 湍流流動時間和平均湍流剪切時間的比率

λm分子導(dǎo)熱系數(shù)，J/（m·s·K）

λt有效導(dǎo)熱系數(shù)，J/（m·s·K）

λturb湍動導(dǎo)熱系數(shù)，J/（m·s·K）

μ 計算溫度下空氣動力黏度，Pa·s

μi在i方向上的速度分量，m/s

μj在j方向上的速度分量，m/s

μt有效黏度系數(shù)，kg/（m·s）

μm分子黏性系數(shù)，kg/（m·s）

μturb湍動黏性系數(shù)，kg/（m·s）

ρ 流體密度，kg/m3

σεk-ε模型方程參數(shù)

［1］ Masoud R，Sayed R S，Ammar A. Experiment and CFD Studies on Heat Transfer and Friction Factor Characheristics of a Tube Equipped with Modified Twisted Inserts［J］. Chem Eng Process，2009，48（3）：762 - 770.

［2］ Sarada Naga S，Sita Rama Raju AV，Kalyani Radha K，et al. Enhancement of Heat Transfer Using Varying Width Twisted Tape Inserts［J］. Int J Eng，Sci Technol，2010，2（6）：107 - 118.［3］ 洪蒙納，鄧先和. 含旋流片縮放管內(nèi)流動與傳熱數(shù)值模擬［J］. 石油化工設(shè)備，2007，36（6）：42 - 45.

［4］ 王楊君，鄧先和. 內(nèi)插旋流片的管內(nèi)流動與換熱的數(shù)值模擬［J］. 化工學(xué)報，2007，58（10）：2455 - 2461.

［5］ 周水洪，鄧先和，王楊君，等. 旋流片強化換熱器殼程傳熱的數(shù)值模擬與實驗［J］. 華南理工大學(xué)學(xué)報：自然科學(xué)版，2007，35（4）：77 - 81.

［6］ Smith E，Chinaruk T，Pongjet P. Experimental Investigation of Heat Transfer and Flow Friction in a Circular Tube Fitted with Regularly Spaced Twisted Tape Elements［J］. Int Commun Heat Mass Transfer，2006，33（10）：1225 - 1233.

［7］ Smith E，Chinaruk T，Petpices E，et al. Convective Heat Transfer in a Circular Tube with Short Length Twisted Tape Insert［J］. Int Commun Heat Mass Transfer，2009，36（4）：365 - 371.

［8］ 王國清. 強化裂解爐管傳熱研究［D］. 北京：北京化工大學(xué)，2001.

［9］ 王國清，曾清泉. 乙烯裂解爐管強化傳熱［J］. 石油化工，2001，30（7）：528 - 530.

［10］ 中國石油化工集團公司，中國科學(xué)院金屬研究所，中國石油化工集團公司北京化工研究院.一種熱交換管及其制造方法和應(yīng)用：中國，1260469A［P］. 2000-07-19.

［11］ 王國清，張利軍，許士興，等. 扭曲片管強化傳熱技術(shù)在裂解爐中的應(yīng)用［J］. 乙烯工業(yè)，2006，18（2）：19 - 22.

［12］ 袁霞光. 三種強化傳熱爐管的流動和傳熱特性的數(shù)值模擬［J］. 石油化工，2011，40（7）：796 - 799.

［13］ 翟俊霞，涂善東，李大驥. 含靜態(tài)混合元件乙烯裂解爐管內(nèi)流場的數(shù)值模擬［J］. 蘭州理工大學(xué)學(xué)報，2004，30（5）：55 - 58.

［14］ 韓檬，張衛(wèi)江，劉豐合. 含扭曲片乙烯裂解爐管內(nèi)流場的數(shù)值模擬［J］. 石油學(xué)報：石油加工，2007，23（3）：63 - 69.

（編輯 王 萍）

Numerical Simulation and Process Calculation of Linear Waste Heat Boiler of Cracking Furnace with Twisted Tapes

Liu Junjie，Wang Guoqing，Zhang Lijun，Zhou Xianfeng
（SINOPEC Beijing Research Institute of Chemical Industry，Beijing 100013，China）

Both Strengthened heat transfer tube with twisted tape(twisted tape tube) and with bond(bond tube)，which were used in linear waste heat boiler in ethylene unit，were simulated numerically，the results were applied to commercial heat transfer software，and the process calculation of the linear waste heat boiler with the twisted tape tube was carried out. The velocity field，pressure drop field and temperature field of the two kinds of strengthened heat transfer tubes were calculated，and the pressure drop，Nu，strengthened heat transfer factors，pressure drop ratio and Nu ratio at different Re were obtained. The simulation results showed that the heat transferring effect of the twisted tape tube was better than that of the bond tube，the pressure drop and Nu of the former were lower than those of the latter，and the strengthened heat transfer factor of the former was bigger than that of the latter. The process calculation results indicated that when the feedstocks were hydrogenated tail oil，naphtha and ethane，compared to common linear waste heat boiler，the capacity of heat transmission of the linear waste heat boiler with the twisted tape tube increased by 8.2%，6.1% and 2.4%，respectively.

bond；twisted tape；heat transfer enhancement；cracking furnace；linear waste heat boiler；numerical simulation；ethylene unit

1000 - 8144（2015）08 - 0919 - 09

TQ 018

A

2015 - 05 - 15；［修改稿日期］ 2015 - 06 - 19。

劉俊杰（1976—），男，山東省寧津縣人，博士，高級工程師，電話 010 - 59202725，電郵 liujj.bjhy@sinopec.com。