李曙陽，伍永福

（1. 呼和浩特中燃城市燃?xì)獍l(fā)展有限公司，內(nèi)蒙古 呼和浩特 010050； 2. 內(nèi)蒙古科技大學(xué)能源與環(huán)境學(xué)院，內(nèi)蒙古 包頭 014010）

氧氣噴射角度對(duì)煤制甲醇轉(zhuǎn)化爐內(nèi)速度場(chǎng)的影響

李曙陽1，伍永福2

（1. 呼和浩特中燃城市燃?xì)獍l(fā)展有限公司，內(nèi)蒙古 呼和浩特 010050； 2. 內(nèi)蒙古科技大學(xué)能源與環(huán)境學(xué)院，內(nèi)蒙古 包頭 014010）

為了研究氧氣噴射角度對(duì)煤制甲醇轉(zhuǎn)化爐中速度場(chǎng)的影響，利用ANSYS和Fluent軟件構(gòu)建煤制甲醇轉(zhuǎn)化爐模型，研究了氧氣噴射角度在0°、8°、11°和15°時(shí)對(duì)轉(zhuǎn)化爐內(nèi)速度分布變化，發(fā)現(xiàn)氧氣噴射角度在0°時(shí)，轉(zhuǎn)化爐內(nèi)速度梯度最小，高速區(qū)域最小且沿軸向延展，火焰最穩(wěn)定，不易出現(xiàn)火焰淬滅現(xiàn)象。

煤制甲醇；速度場(chǎng)；氧氣噴射角度；數(shù)值模擬

煤制甲醇是煤化工的重要分支，目前普遍采用的煤制甲醇工藝是以焦?fàn)t煤氣為主要原料，通過在轉(zhuǎn)化爐中的轉(zhuǎn)化，生成甲醇。為了解各組分在轉(zhuǎn)化爐中的轉(zhuǎn)化情況，利用數(shù)值模擬方法研究了全氧轉(zhuǎn)化爐在不同氧氣噴射角度時(shí)的速度分布情況。

1 轉(zhuǎn)化爐模型的建立

1.1 物理模型的建立及其網(wǎng)格化

本文數(shù)學(xué)模擬是針對(duì)全氧轉(zhuǎn)化爐，爐型上部為蓋帽、上端兩個(gè)進(jìn)口為進(jìn)甲烷氣體（斷面直徑為2.5 mm），下面兩個(gè)斷面進(jìn)口為進(jìn)氧氣（斷面直徑為1.5 mm）。燒嘴將轉(zhuǎn)化爐隔開分為上下兩個(gè)爐膛，分別為上部的氣體混合室和下半部的燃燒室。氣體進(jìn)口方向垂直于混合氣體流向[1]，使用ANSYS14.0構(gòu)建的轉(zhuǎn)化爐物理模型如圖1所示，計(jì)算區(qū)域網(wǎng)格化后如圖2所示。以爐內(nèi)傳熱模型為核心，由燃燒模型、爐內(nèi)爐氣流動(dòng)與爐內(nèi)傳熱模型、鋼坯加熱模型耦合而成，這些模型與它們的定解條件構(gòu)成全部的數(shù)學(xué)模型及其求解問題。燃燒過程模型是爐內(nèi)傳熱模型的前提，爐內(nèi)爐氣流動(dòng)與爐內(nèi)傳熱模型互為關(guān)聯(lián)，鋼坯加熱模型是爐內(nèi)傳熱模型的結(jié)果?？紤]到數(shù)學(xué)模型的計(jì)算機(jī)數(shù)值求解的可行性和數(shù)學(xué)模型的可靠性，在數(shù)學(xué)模型的建立過程中，進(jìn)行如下假設(shè)：1）以轉(zhuǎn)化爐的中心截面為對(duì)稱面，忽略爐肋管與鋼坯之間的傳熱。2）爐頭、爐尾無逸氣和吸氣。3）將焦?fàn)t煤氣視為甲烷單一組分。

由于在轉(zhuǎn)化爐中焦?fàn)t煤氣和氧氣混合時(shí)，焦?fàn)t煤氣與氧氣都是由各自的通道分別射入混合室，初始射流之間有一定的距離，兩者相遇后混合但不燃燒，可建立湍流模型。常見的湍流模型有Spakart-Allmaras模型、k-ε模型、大渦模擬模型。本研究使用Fluent軟件建立的湍流模型為k-ε模型，選用的算法是經(jīng)典的SIMPLE算法，并用二階迎風(fēng)格式計(jì)算各組分精度。本模型的邊壁條件為：5個(gè)管口：分別為2個(gè)焦?fàn)t煤氣入口、2個(gè)氧氣入口和一個(gè)混合氣體出口；其余為壁面。各入口選擇速度入口作為邊界條件（Velocity-inlet）；由于對(duì)轉(zhuǎn)化爐出口的煙氣沒有要求，常將其邊界條件選為自由流體出口邊界條件，本研究中出口采用自由出口邊界條件（out-flow）。進(jìn)出口處的湍流參數(shù)主要有湍流強(qiáng)度參數(shù)I(Turbulent Intensity)和水力直徑Φ(Hydraulic Diameter)[2]，具體見公式（1）和公式2)。由于轉(zhuǎn)化爐為鋼坯制成，本研究中將轉(zhuǎn)化爐爐體假定為絕緣材料，壁面條件設(shè)定為Wall，各參數(shù)采用FLUENT中的默認(rèn)條件。

式中：

Re —雷諾數(shù)；

Φ —水力直徑，m。

圖1 轉(zhuǎn)化爐燃燒器物理模型Fig.1 The physical model of reforming furnace burner

圖2 轉(zhuǎn)化爐燃燒器模型網(wǎng)格劃分Fig.2 The grid division of reforming furnace burner model

1.2 湍流流動(dòng)數(shù)學(xué)模型的建立

煤制甲醇轉(zhuǎn)化爐熱過程中的爐膛內(nèi)氣體流動(dòng)、對(duì)流傳熱、傳質(zhì)及燃燒過程均為湍流過程，因此，為了能夠定量描述爐膛內(nèi)的湍流均流場(chǎng)，必須建立均流場(chǎng)的控制微分方程組。目前以半經(jīng)驗(yàn)理論為基礎(chǔ)，廣泛研究和應(yīng)用的湍流輸運(yùn)模型基本上是圍繞輸運(yùn)通量進(jìn)行了的，對(duì)轉(zhuǎn)化爐內(nèi)氣體流動(dòng)的研究通常是引入湍流粘性系數(shù)的概念，以普朗特混合長(zhǎng)度理論為依據(jù)而建立kε-雙方程模型。湍流流動(dòng)kε-雙方程模型在實(shí)際工程應(yīng)用中取得了巨大的成就，給出了大批與實(shí)驗(yàn)吻合的數(shù)值解，使越來越多的科技人員認(rèn)識(shí)到湍流kε-雙方程模型的可靠性和實(shí)用性，是目前應(yīng)用最廣泛的工程模型。用kε-雙方程湍流模型來描述煤制甲醇轉(zhuǎn)化爐爐膛內(nèi)爐氣流動(dòng)過程是可行的[5-6]。模型的相關(guān)方程如公式(3-8)所示，

連續(xù)方程：

動(dòng)量方程：

能量方程：

質(zhì)量方程：

k方程：

ε方程：

式中： μ1=Cμρk2/ε；μeff= μ0+μ1；Cμ=0.09；CD=0.08～0.38；C1=1.44；C2=1.92。ρ為密度（kg/m3）；t為時(shí)間（s）；xj為直角坐標(biāo)j方向的坐標(biāo)；uj為直角坐標(biāo)系j方向上的速度（m/s）；P為壓力（Pa）；μ0為分子粘性系數(shù)（Pa·S）；μ1為湍流動(dòng)力粘性系數(shù)（Pa·S）；μeff為湍流有效粘性系數(shù)（Pa·S）；k為湍流脈動(dòng)動(dòng)能；Cμ為無量綱系數(shù)；CD為無量綱系數(shù)；l為湍流的脈動(dòng)普照朗特混合長(zhǎng)度；ε為湍流脈動(dòng)動(dòng)能耗散系數(shù)；C1、C2為無量綱系數(shù)；H為熱焓（kJ/(kg·K)）；Гh為熱擴(kuò)散系數(shù)（W/(m·K)）；是熱焓源項(xiàng)（W/m3）；ms為組分（%）；Гs為組分?jǐn)U散系數(shù)（m2/s）；為反應(yīng)源項(xiàng)（1/s）。

2 數(shù)值模擬結(jié)果及分析

為了考察模型的可靠性，對(duì)模型的收斂情況進(jìn)行考察，經(jīng)初始化運(yùn)算，經(jīng)過對(duì)模型欠松弛因子的設(shè)定，最終得到所建立模型的殘差曲線，如圖3所示。由圖3可以看出，所建立模型的殘差曲線逐漸收攏逼近于X軸，可預(yù)見的是，隨著運(yùn)算的進(jìn)行，模型誤差逐漸趨于0，模型具有收斂性。

在煤制甲醇轉(zhuǎn)化爐中，焦?fàn)t煤氣與氧氣提前預(yù)混可以使焦?fàn)t煤氣得到充分的燃燒，可預(yù)見的是，在預(yù)混環(huán)節(jié)混合越均勻焦?fàn)t煤氣燃燒越徹底，煤制甲醇的轉(zhuǎn)化率越高。

圖3 轉(zhuǎn)化爐燃燒器模型的殘差曲線Fig.3 The residual curve of reforming furnace burner model

當(dāng)氧氣以不同角度射入混合倉時(shí)，氣體混合程度會(huì)發(fā)生變化，氧氣進(jìn)口的角度改變會(huì)改變混合氣體的速度，濃度，壓力的分布。本研究主要討論氧氣進(jìn)口角度不同時(shí)轉(zhuǎn)化爐中壓力場(chǎng)的分布情況。為了考察氧氣進(jìn)口角度對(duì)轉(zhuǎn)化爐壓力場(chǎng)的影響，模擬了氧氣噴射角度在0°、8°、11°和15°時(shí)轉(zhuǎn)化爐內(nèi)壓力場(chǎng)分布情況。模擬結(jié)果如圖4-7所示，對(duì)應(yīng)的各入射角度在Y軸上壓力分布曲線如圖8-11所示。由圖4、6、8、10可以看出，氧氣噴射角度不同時(shí)，在轉(zhuǎn)化爐內(nèi)，速度場(chǎng)的分布差異較大，尤其是0°入射角時(shí)，低速區(qū)占比明顯高于8°、11°、15°時(shí)的低速區(qū)；就轉(zhuǎn)化爐內(nèi)的流速梯度而言，0°入射角時(shí)速度梯度也明顯低于8°、11°、15°時(shí)的流速梯度。

圖4 氧氣0°入射時(shí)轉(zhuǎn)化爐燃燒器內(nèi)的速度場(chǎng)分布Fig.4 The velocity field of reforming furnace burner as the injection angle of oxygen is 0°

圖5 氧氣0°入射時(shí)轉(zhuǎn)化爐燃燒器YZ平面上速度場(chǎng)分布Fig.5 The velocity field of reforming furnace burner on YZ plane as the injection angle of oxygen is 0°

圖6 氧氣8°入射時(shí)轉(zhuǎn)化爐燃燒器內(nèi)的速度場(chǎng)分布Fig.6 The velocity field of reforming furnace burner as the injection angle of oxygen is 8°

圖7 氧氣8°入射時(shí)轉(zhuǎn)化爐燃燒器YZ平面上速度場(chǎng)分布Fig.7 The velocity field of reforming furnace burner on YZ lane as the injection angle of oxygen is 8°

圖8 氧氣11°入射時(shí)轉(zhuǎn)化爐燃燒器內(nèi)的速度場(chǎng)分布Fig.8 The velocity field of reforming furnace burner as the injection angle of oxygen is 11°

圖9 氧氣11°入射時(shí)轉(zhuǎn)化爐燃燒器YZ平面上速度場(chǎng)分布Fig.9 The velocity field of reforming furnace burner on YZ plane as the injection angle of oxygen is 11°

圖10 氧氣15°入射時(shí)轉(zhuǎn)化爐燃燒器內(nèi)的速度場(chǎng)分布Fig.10 The velocity field of reforming furnace burner as the injection angle of oxygen is 15°

圖11 氧氣15°入射時(shí)轉(zhuǎn)化爐燃燒器YZ平面上速度場(chǎng)分布Fig.11 The velocity field of reforming furnace burner on YZ plane as the injection angle of oxygen is 15°

由圖5、7、9、11、12、13、14、15可以看出，氧氣噴射角度在0°時(shí)，轉(zhuǎn)化爐內(nèi)的大部分區(qū)域處于低速區(qū)，低速區(qū)的流速范圍是0～6.691 m/s，在轉(zhuǎn)化爐中央?yún)^(qū)域存在一個(gè)中流速區(qū)域，區(qū)域流速范圍是6.691～13.38 m/s，在上下爐膛分界處存在一個(gè)中心高流速區(qū)域，該區(qū)域流速范圍是13.38～20.07 m/s，主要沿軸向延伸，徑向分布范圍有限，在軸向上，該高速區(qū)的最高流速在12.5 m/s左右；

圖12 氧氣0°入射時(shí)轉(zhuǎn)化爐燃燒器Y軸上速度分布Fig.12 The velocity field of reforming furnace burner on Y axis as the injection angle of oxygen is 0°

圖13 氧氣8°入射時(shí)轉(zhuǎn)化爐燃燒器Y軸上速度分布Fig.13 The velocity field of reforming furnace burner on Y axis as the injection angle of oxygen is 8°

圖14 氧氣11°入射時(shí)轉(zhuǎn)化爐燃燒器Y軸上速度分布Fig.14 The velocity field of reforming furnace burner on Y axis as the injection angle of oxygen is 11°

圖15 氧氣15°入射時(shí)轉(zhuǎn)化爐燃燒器Y軸上速度分布Fig.15 The velocity field of reforming furnace burner on Y axis as the injection angle of oxygen is 15°

氧氣噴射角度在8°時(shí)，轉(zhuǎn)化爐內(nèi)的大部分區(qū)域處于低速區(qū)，低速區(qū)的流速范圍是0～7.505 m/s，高于0°噴射角時(shí)的低速速度范圍，所占區(qū)域也明顯小于0°噴射角時(shí)的分布區(qū)域，在轉(zhuǎn)化爐中央?yún)^(qū)域存在一個(gè)中流速區(qū)域，區(qū)域流速范圍是7.505～15.01 m/s，這一速度也高于0°噴射角時(shí)速度范圍，而所占區(qū)域則高于0°噴射角時(shí)的中流速區(qū)域，在上下爐膛分界處存在一個(gè)中心高流速區(qū)域，該區(qū)域流速范圍是15.01～22.51 m/s，呈圓球狀，在軸向上，該高速區(qū)的最高流速在14.3 m/s左右。

當(dāng)氧氣噴射角度在11°和15°時(shí)，轉(zhuǎn)化爐內(nèi)的速度場(chǎng)分布與8°時(shí)的分布規(guī)律相類似，且呈現(xiàn)爐內(nèi)速度梯度進(jìn)一步加劇、中高速區(qū)域進(jìn)一步增加、高速區(qū)徑向延展、中高速區(qū)軸向延展且上移的基本規(guī)律。而爐內(nèi)速度梯度的增加雖然能在一定程度上加快燃燒速度，但是也必將加大轉(zhuǎn)化爐內(nèi)的燃燒和火焰不穩(wěn)定性，再加上中高速流速區(qū)域的向上爐膛的移動(dòng)，對(duì)燃燒過程起到極大的抑制作用[7]，另一方面，在上下爐膛分界處所存在的高速區(qū)的徑向延展，將可能導(dǎo)致燃燒的淬滅，使轉(zhuǎn)化過程終止[8]，這對(duì)應(yīng)焦?fàn)t煤氣制甲醇過程顯然是不利的。

3 結(jié) 論

利用k-ε模型構(gòu)建的煤制甲醇轉(zhuǎn)化爐數(shù)學(xué)模型模擬結(jié)果顯示，在焦?fàn)t煤氣進(jìn)氣口方向?yàn)?°時(shí)，不同的氧氣噴射角度會(huì)對(duì)轉(zhuǎn)化爐爐膛內(nèi)的速度場(chǎng)分布產(chǎn)生較大的影響；在氧氣噴射角度為0°、8°、11°、15°時(shí)，在上下爐膛分界處均存在一個(gè)高速流速區(qū)；相比較于8°、11°、15°的氧氣入射角，0°氧氣入射角的高速流速區(qū)域更小，主要沿軸向延展，高速區(qū)域的速度最低、范圍最小；加大氧氣入射角度會(huì)加大轉(zhuǎn)化爐內(nèi)的中高速區(qū)域占比、中高速區(qū)域上移，高速區(qū)域范圍擴(kuò)大和徑向延展，這將可能導(dǎo)致轉(zhuǎn)化爐中的燃燒淬滅和轉(zhuǎn)化過程的終止，氧氣的噴射角以0°為宜。

［1］張固. 甲醇新鮮合成氣氫碳比的優(yōu)化方法探討[J]. 化工設(shè)計(jì), 2004, 05: 6-10.

［2］龔曙光, 黃云清. 有限元分析與ANSYS APDL編程及高級(jí)應(yīng)用[M].北京：機(jī)械工業(yè)出版社, 2009.

［3］趙博寧, 李健. 天然氣預(yù)熱溫度對(duì)高溫蓄熱式加熱爐內(nèi)燃?xì)饬鲃?dòng)影響的數(shù)值模擬研究[J]. 化學(xué)工程與裝備, 2010, 08: 27-29.

［4］李彥鵬, 王金枝, 許晉源. 切向燃燒鍋爐燃燒器區(qū)流場(chǎng)及壁面負(fù)壓的試驗(yàn)與數(shù)值研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2000, 01: 66-70.

［5］嚴(yán)陣. 低熱值高溫煤氣燃燒器流場(chǎng)冷態(tài)試驗(yàn)及模擬研究[D]. 華中科技大學(xué), 2011.

［6］蘇俊林, 徐曉英, 王震坤. 切向鼓風(fēng)式燃燒器三維流場(chǎng)的數(shù)值模擬[J]. 冶金能源, 2008, 03: 30-32.

［7］龍鈺, 楊向平, 尹恩杰, 等. 制氫轉(zhuǎn)化爐內(nèi)綜合數(shù)值模擬研究[J].煉油技術(shù)與工程, 2010, 08: 22-25.

［8］周桂娟, 蔡建光. 不同CO2抽出比的PSA尾氣在制氫轉(zhuǎn)化爐內(nèi)燃燒的數(shù)值模擬研究[J]. 石油化工設(shè)備技術(shù), 2012, 03: 35-38；42.

Influence of Injection Angle of Oxygen on Velocity Field in Conversion Furnace of Coal-methanol

LI Shu-yang1，WU Yong-fu2
（1. Hohhot Zhongran City Gas Development Co., Ltd.，Inner Mongolia Hohhot 010050, China； 2. School of Energy and Environmental Engineering，Inner Mongolia University of Science and Technology, Inner Mongolia Baotou 014010, China）

In order to study the influence of injection angle of oxygen on velocity field in conversion furnace of coal-methanol, the mathematical model of conversion furnace of coal-methanol was built with ANSYS and Fluent. The pressure field of conversion furnace of coal-methanol was simulated as the injection angle was 0°，8°，11° or 15°.The results show that the minimum velocity gradient can appear as the injection angle is 0°,the high speed area is the smallest and the region extends along the axial direction, the combustion of coke oven gas is the stablest,the quenching of flame is not easy to happen.

Coal-methanol; Velocity field; Injection angle of oxygen; Numerical simulation

TQ 545

A

： 1671-0460（2015）01-0174-04

國(guó)家自然科學(xué)基金項(xiàng)目，項(xiàng)目號(hào)：51464041。

2014-12-10

李曙陽（1979-），男，內(nèi)蒙古自治區(qū)呼和浩特人，工程師，碩士，1999年畢業(yè)于內(nèi)蒙古工業(yè)大學(xué)城市燃?xì)夤こ虒I(yè)，研究方向：從事城市燃?xì)馔顿Y建設(shè)與運(yùn)營(yíng)管理技術(shù)工作。E-mail：13948910707@163.com。

伍永福（1974-），男，副教授，博士，熱工設(shè)備工藝過程仿真模擬及優(yōu)化。E-mail：wyf07@imust.cn。