第38卷第5期化學(xué)工程Vol. 38 No.52010年5月CHEMICAL ENGINEERING( CHINA)May 2010Texaco氣化爐冷態(tài)流場(chǎng)和湍流混合的數(shù)值模擬曹月叢，陳彩霞,高晉生(華東理工大學(xué)資源與環(huán)境工程學(xué)院，上海200237)摘要:在直徑為1000 mm,高4000 mm的雙通道射流噴嘴的氣化爐上,以Texaco氣化爐冷模試驗(yàn)為基準(zhǔn)對(duì)象,將空氣經(jīng)環(huán)隙和中心射人氣化爐,以氫氣為示蹤劑，預(yù)測(cè)環(huán)隙和中心射流的混合程度。結(jié)合流體質(zhì)量與動(dòng)量守恒方程和k-e湍流模型,用SIMPLER算法計(jì)算,對(duì)氣化爐內(nèi)的冷態(tài)流場(chǎng)和湍流混合進(jìn)行了模擬,模擬了爐內(nèi)速度分布量綱一濃度分布、混合分?jǐn)?shù)分布和軸向衰減的情況。結(jié)果顯示:氣化爐內(nèi)濃度分布極不均勻;爐內(nèi)存在富氧和貧氧區(qū);環(huán)院和中心通道射流動(dòng)最比加大,混合分?jǐn)?shù)沿軸向衰減加快,達(dá)到充分混合的時(shí)間縮短。模擬結(jié)果與冷模試驗(yàn)結(jié)果的比較表明計(jì)算值與試驗(yàn)值吻合良好。關(guān)鍵詞:Texaco;氣化爐;數(shù)值模擬;冷態(tài)流場(chǎng);湍流中圖分類號(hào):TQ 021. I文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1005-9954(2010)05 0030.04Numerical simulation of cold-flow and turbulent mixing in Texaco gasifierCAO Yue-cong, CHEN Cai-xia, GAO Jin-sheng(School of Resource & Environmental Engineering, East China University ofScience and Technology , Shanghai 200237, China)Abstract:A dual-channel jet nozle cold model for Texaco gasifier with diameter of 1 000 mm and height of4 000 mm was used as the benchmark. The air was injected into the gasifier through the annular and concentricchannel of the nozzle ，and the hydrogen was injected as a tracer for tracking the extent of mixing. The SIMPLERalgorithm was used to solve the fluid mass and momentum conservation equations which were coupled with theturbulent model. The cold-low field and turbulent mixing in the gasifier were simulated. The simulation calculatedthe distributions of gas velocity, dimension one concentration and mixture fraction, and the decay along the axialdirection. The results show that the distribution of gas concentration is non-uniform, and a rich and a poor oxygenregion exist in the gasifier. The mixture fraction decreases quickly along the axial direction, and the time requiredfor perfeet mixing shortens with increasing jet momentum ratio of annular to concentric jet. The simulation resulsagree with the experiment of cold flow.Key words :Texaco; gasifier; numerical simulation; cold-flow; turbulent flowTexaco煤氣化是由美國(guó)Texaco公司在重油氣擬。 上述數(shù)值模擬不能很好地模擬氣化爐內(nèi)的湍流化基礎(chǔ)_上開(kāi)發(fā)成功的第二代煤氣化技術(shù),自1993年混合情況和組分濃度分布,模擬也沒(méi)有得到有效驗(yàn)在我國(guó)實(shí)現(xiàn)工業(yè)化運(yùn)行。Texaco氣化過(guò)程的試驗(yàn)研證,需要開(kāi)發(fā)氣流床整體數(shù)值模型。究主要集中在冷態(tài)過(guò)程,有冷態(tài)流場(chǎng)的速度分Texaco氣流床氣化爐內(nèi)湍流流動(dòng)與混合過(guò)程起布"、濃度分布2]等。近年來(lái),用商業(yè)軟件如Fluent決定作用。作為開(kāi)發(fā)整體數(shù)值模型的第1步,本文對(duì)Texaco氣化爐的數(shù)值模擬引人注目。賀阿特對(duì)氣化爐內(nèi)的冷態(tài)流動(dòng)和湍流混合進(jìn)行了模擬,結(jié)等[3]用簡(jiǎn)化PDF模型對(duì)渣油復(fù)雜氣化過(guò)程進(jìn)行了合流體流動(dòng)質(zhì)量與動(dòng)量方程和k-8湍流模型,計(jì)算二維模擬;劉向軍等(引人渦量-流函數(shù)、k-e模型模出流場(chǎng)速度分布、濃度分布和混合分?jǐn)?shù)的變化,并對(duì)擬Texaco氣化爐內(nèi)煤氣化過(guò)程;吳玉新等(S)用簡(jiǎn)化揩蟻結(jié)果與試驗(yàn)結(jié)果PDF模型對(duì)Texaco水煤漿氣化爐進(jìn)行三維數(shù)值?；局袊?guó)煤化工床整體數(shù)值模型HHCNMHG基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)日(20876049)作者簡(jiǎn)介:曹月叢(1983- -),女碩士研究生,研究方向?yàn)槊簹饣^(guò)程的數(shù)值模擬;陳彩霞,通訊聯(lián)系人,電話: (021 )64251092, E-mil:exchen@ ecust. edu. cn。曹月叢等Texaco 氣化爐冷態(tài)流場(chǎng)和湍流混合的數(shù)值模擬●31●方面邁開(kāi)了堅(jiān)實(shí)的一步。.H,先后加入中心和環(huán)隙流體中后試驗(yàn)測(cè)定的同一位置處H2質(zhì)量分?jǐn)?shù),.和qm.o。分別為噴嘴中心1試驗(yàn)分析通道和環(huán)隙通道流體的質(zhì)量流量。1.1試驗(yàn)裝置及試驗(yàn)條件充分混合時(shí)的混合分?jǐn)?shù)試驗(yàn)的主要裝置為直徑1000mm,高4000mm.f. =nmo/(qms.n +9m.on)(5)的有機(jī)玻璃模型爐。噴嘴中心通道直徑為49.5 mm ,外環(huán)直徑為112.5 mm。H2為示蹤劑,空2數(shù)值模擬氣為氣化介質(zhì),測(cè)定氣化爐內(nèi)的濃度分布。根據(jù)工2.1控制方程業(yè)射流速度和環(huán)隙中心射流動(dòng)量比,不同體積流量控制方程包括質(zhì)量、動(dòng)量和能最守恒方程,這里的空氣分別通過(guò)噴嘴的中心和環(huán)隙,射人氣化爐。為冷態(tài)模擬,且是均-氣相,所以可以不考慮能量守氫氣先通過(guò)中心進(jìn)人氣化爐,再通過(guò)環(huán)隙進(jìn)入氣化恒和顆粒源項(xiàng)。由瞬時(shí)方程得出的Favre平均方程爐,保證二次所進(jìn)的氫氣體積流量相等。試驗(yàn)條件如下門(mén):見(jiàn)表1[6] ,qv,;為中心射流體積流量, u;為中心射流連續(xù)性方程d(pn) =0.(6)流速,qr,。為環(huán)隙射流體積流量,u。為環(huán)隙射流流x速。動(dòng)量方程表1試驗(yàn)條件武(西可)=-啞+(新+到)一點(diǎn)()ax; 'Table 1 Experimental conditions(7)中心射流環(huán)隙射流氣相組分方程9v;/9v./u。[/動(dòng)量比(m'.h1) (m.s~')_ (m'.h") (m.s")2(pJ)=-2(puY)(8)1 382.560.030.01.0444湍流動(dòng)能方程2574.090.0799.030.0 0. 46403765.0120.00.261 1然=點(diǎn)(告路)-,啊-e (9)“.xsxσx，x-4382.5985.037.01.5880湍動(dòng)能耗散率5574.0985. 00.705 56765.0120. 037.0 0. 397 0王(出de)+(影)畫(huà)吧+C(周)(10)1.2氣化爐內(nèi)濃度分 布的計(jì)算采用與試驗(yàn)一致的定 量方法[6] ,定義了量綱一式中:σk,σ,為交換系數(shù),分別為0.9和1.22;C,Cz為模型常數(shù),分別為1.44和1.92;p為氣體密度,濃度η:n:=(y;-y.)/y.(1)kg/m3 ;u為氣體速度, m/s;xr,x,xp為三維空間的坐式中:y為爐內(nèi)任意位置中心流體的摩爾分?jǐn)?shù),標(biāo),m;p為靜壓,Pa;上標(biāo)”為動(dòng)量的變動(dòng),“一”為均為以噴嘴中心通道流體為基準(zhǔn)計(jì)算的充分混合時(shí)的值;Y,為不同氣體組分的產(chǎn)率;μ為湍流黏度，摩爾分?jǐn)?shù)。.kg/(m.s);k為湍流動(dòng)能,m2/s2 ;&為湍流動(dòng)能耗y. =qv,/(qv.s +qr,o) .(2)散率,m-2.s-3 ;p u"u"和p u"Y,為湍流擴(kuò)散項(xiàng)。y: =C;9v,/(cq9v. +c.9v.)(3)2.2氣相湍流模型 .式中:c;和c。分別為H2先后加入中心和環(huán)隙流體湍流是通:過(guò)速度分量和一些標(biāo)量的時(shí)空波動(dòng)來(lái)后試驗(yàn)測(cè)定的同- -位置H,濃度。表征的。在控制方程中,雷諾時(shí)均應(yīng)力和雷諾通量1.3混合分?jǐn)?shù)的計(jì)算通過(guò)k-8模型計(jì)算”]。模型中的參數(shù)使用標(biāo)準(zhǔn)混合分?jǐn)?shù)定義為系數(shù)。9ni+9m. w9m.s +wo9m.onf=- 9m,wq9m.,.n(4)2.3中國(guó)煤化工上HCNMH G系表示"”。微分式中q9m.,為氣化爐內(nèi)某點(diǎn)中心流體的質(zhì)量流量,9m,o方程被轉(zhuǎn)化為有限差分萬(wàn)程計(jì)算,區(qū)域網(wǎng)格劃分后，為相應(yīng)點(diǎn)上環(huán)隙流體的質(zhì)量流量,w;和w。分別為用SIMPLER算法進(jìn)行網(wǎng)格內(nèi)的迭代計(jì)算。●32.化學(xué)工程2010年第38卷第5期.3網(wǎng)格劃分及模擬條件4.2濃度分布的模擬結(jié)果3.1網(wǎng)格劃分對(duì)試驗(yàn)2和試驗(yàn)6進(jìn)行模擬,計(jì)算出w0;和w。,模擬直徑1 000 mm,高4000 mm的氣化爐,拱然后轉(zhuǎn)化成c,和c。,利用式(1)-(3),計(jì)算出量綱頂為直徑1 000 mm的半球體,直段高度3 500 mm，- -濃度n,試驗(yàn)結(jié)果參考文獻(xiàn)[6]。在氣化爐的不出口在直簡(jiǎn)底部。若以31 x31的網(wǎng)格,計(jì)算結(jié)果不同軸向高度處,研究了η沿徑向的分布。試驗(yàn)2和夠精確。所以,網(wǎng)格劃分以51 x51 x 130,局部進(jìn)6的模擬與試驗(yàn)結(jié)果的比較見(jiàn)圖2和圖3。其中,r行網(wǎng)格細(xì)化,分別在出口高度0- -20 mm處,軸向?yàn)闅饣癄t內(nèi)離中心軸線的徑向距離。序號(hào)1- 4對(duì)3個(gè)網(wǎng)格,在高度3000-4000mm處，軸向29個(gè)應(yīng)的軸向高度分別為離噴嘴距離17. 2,41.1,65.5,網(wǎng)格。105.5 cm。3.2射流條件2.0[ ■試驗(yàn)值模擬值模擬時(shí)H2體積流量是中心通道空氣體積流量的1%,由表1計(jì)算出中心和環(huán)隙通道空氣的質(zhì)量+pr n流量和H2的質(zhì)量分?jǐn)?shù),列于表2中。模擬條件為-1.06 102030400 10203040 0 102040 0 1020304050r/cm298 K, 100 kPa。迭代次數(shù)為3000次。圖2試驗(yàn)2的模擬與試驗(yàn)結(jié)果比較Fig.2 Contrast of simulation to experimental reult of Tet 2表2射流條件Table 2 Jet conditions.[■試驗(yàn)值1[3[中心射流環(huán)隙射流一模擬值9m.s.n/w/10-*min/w./10-*動(dòng)量比(kg.s")(kg.s)-0.1 0.1376.9710. 2873.3281.044-0 102030400 102030400 10203040 0 102030405020. 2066. 9560.2874. 9930. 464030. 2756. 9456. 6550.261 1圈3試驗(yàn)6的模擬與試驗(yàn)結(jié)果比較40.1376. 9710.3542. 6981.588 0Fig.3 Contrast of simulation to experimental result of Test 650.2064.0480. 705560.2750.397 0試驗(yàn)2和6的模擬與試驗(yàn)結(jié)果基本吻合,軸向和徑向濃度梯度顯著,隨著軸向距離的增加,混合越充分,當(dāng)η;>0,為富氧區(qū);當(dāng)η; <0,為貧氧區(qū)。4模擬結(jié)果及討論為便于與試驗(yàn)結(jié)果比較,只對(duì)試驗(yàn)2和6進(jìn)行了模4.1流場(chǎng)分布圖1所示分別為氣化爐內(nèi)的流線圖和速度分布擬。圖。從圖中可以看出氣化爐內(nèi)的射流區(qū),回流區(qū)和4.3混合分?jǐn)?shù)變化的模擬結(jié)果管流區(qū)。噴嘴附近速度最大,中心軸線附近速度大,對(duì)試驗(yàn)3和試驗(yàn)5進(jìn)行模擬,計(jì)算出w;和wo,根據(jù)表2,利用式(4) ,計(jì)算出混合分?jǐn)?shù),以便與試驗(yàn)距離噴嘴遠(yuǎn)的地方,速度漸小。結(jié)果[')對(duì)比。在氣化爐的不同軸向高度處,研究了混合分?jǐn)?shù)f沿徑向的分布。試驗(yàn)3和5的模擬結(jié)果與試驗(yàn)結(jié)果的比較如圖4和圖5所示。序號(hào)1-4209對(duì)應(yīng)的軸向高度分別為離噴嘴17.2,41.1,65. 5,23.316.7 :1.04 ■試驗(yàn)值13.3 昌08一模擬值10.0.66.67圖4r廿3.33中國(guó)煤化工304001020 30 4050(a)流線(b)速度^fYHCNMHG圈1流線與速度分布田4 試驗(yàn)3的模擬與試驗(yàn)結(jié)果比較Fig. 1 Streamline and velocity distributionFig. 4 Contast of simulation to experimental result of Test 3曹月叢等Texaco 氣化爐冷態(tài)流場(chǎng)和湍流混合的數(shù)值模擬●33●1.09●試驗(yàn)值|[3[由圖7得出,環(huán)隙與中心通道射流動(dòng)量比加" 0.0.8 --模擬值大,混合分?jǐn)?shù)沿軸向衰減加快,達(dá)到充分混合點(diǎn)的距離距噴嘴越近,混合情況就好。所以,適度01020304001020304001020304001020304050增加環(huán)隙射流的動(dòng)量,有利于改善爐內(nèi)的混合過(guò)程。圈5試驗(yàn)5的模擬與試驗(yàn)結(jié)果比較Fig.5 Contrast of simulation to experimental result of Test 55結(jié)論試驗(yàn)3與試驗(yàn)5的模擬結(jié)果基本與試驗(yàn)結(jié)果相(1)模擬了Texaco氣化爐內(nèi)的冷態(tài)濃度分布和符合,在距離噴嘴65.5 cm附近,即在拱頂以下一定混合情況。結(jié)果表明,計(jì)算值與試驗(yàn)值吻合較好。距離附近,有回流區(qū),用k-E模型各向同性的假定不證明了本文的模型及計(jì)算程序是可靠的。能精確模擬回流區(qū)各向異性的湍流結(jié)構(gòu),致使模擬(2)冷態(tài)氣化爐內(nèi)的軸向和徑向濃度分布極不結(jié)果與試驗(yàn)結(jié)果-致性不太好,還有待于進(jìn)一步的均勻,混合分?jǐn)?shù)的變化也表明,氣化爐內(nèi)混合情況不研究。當(dāng)f>f.時(shí),為富氧區(qū),f