天然氣重整器內(nèi)部傳遞過程的數(shù)值模擬

《天然氣重整器內(nèi)部傳遞過程的數(shù)值模擬》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、維普資訊http://www.cqvip.com弟28巷弟3別大然氣工業(yè)天然氣重整器內(nèi)部傳遞過程的數(shù)值模擬*楊國剛岳丹婷袁金良。(1．大連海事大學(xué)輪機(jī)學(xué)院2．瑞典蘭德大學(xué)熱能工程系)楊國剛等．天然氣重整器內(nèi)部傳遞過程的數(shù)值模擬．天然氣工業(yè)，2008，28(3)：121—123．摘要天然氣通過重整反應(yīng)可以生成氫氣，成為燃料電池氫原料的重要來源。目前緊湊型天然氣重整器的開發(fā)設(shè)計，主要集中在材料方面的研究，而內(nèi)部機(jī)理方面的研究較少。對天然氣重整器內(nèi)部流動與傳熱現(xiàn)象進(jìn)行了模擬與分析。分析中考慮了化學(xué)反應(yīng)的影響，采用了耦合的邊界條件以及可變的熱物性參數(shù)

2、，建立了描述流通管道和多孔催化劑層內(nèi)部流動與傳熱現(xiàn)象的三維數(shù)學(xué)模型。采用SIMPLE算法對模型進(jìn)行求解，得到了化學(xué)反應(yīng)速率、混合氣濃度以及溫度等的空間分布。結(jié)果表明，化學(xué)反應(yīng)限制在厚度約為1mm的催化劑薄層內(nèi)進(jìn)行，沿主流動方向，CH濃度由19．5降到7左右，H濃度由4．1上升到13左右，混合氣溫度上升約1O℃。研究結(jié)果對天然氣重整器的開發(fā)、結(jié)構(gòu)設(shè)計具有參考意義。主題詞燃料電池天然氣重整器傳熱傳質(zhì)數(shù)值模擬天然氣(甲烷占97以上)通過重整反應(yīng)生成的氫氣，可成為燃料電池氫原料的重要來源。近年出現(xiàn)的緊湊型天然氣重整器由燃燒與重整兩個平行管道組成，具有

3、熱傳導(dǎo)的距離以及氣體擴(kuò)散路徑短的優(yōu)點，降低了重整器中的傳熱、傳質(zhì)阻力。因此，該重整器不僅能量利用率較高，也使得體積大大縮小l_1]。目前緊湊型天然氣重整器的內(nèi)部機(jī)理方面的研究較少口]。筆者將對單個重整管道的不同傳遞過程進(jìn)行三維模擬，不涉及燃燒管道，燃燒管道對嘲是供的鞋瞳(q)重整管道的影響將采用熱邊界條件來描述(相鄰燃圖1緊湊型重整器結(jié)構(gòu)示意圖燒管道固壁有熱量流入，其他為絕熱條件)。一式中：△。。。表示反應(yīng)溫度為1000K時的反應(yīng)熱，、過程描述與數(shù)學(xué)模型方程(1)稱為水蒸氣重整反應(yīng)，為吸熱反應(yīng)；方程(2)重整器結(jié)構(gòu)見圖1。取一半寬度管道進(jìn)行研

4、究，稱為水一氣變換反應(yīng)，為放熱反應(yīng)?？偡磻?yīng)吸熱，熱坐標(biāo)原點設(shè)定在復(fù)合管道的左下角處，坐標(biāo)、量由燃燒管道內(nèi)的催化燃燒反應(yīng)提供。、分別表示管道的長、高、寬，U、、w分別為對應(yīng)3個控制方程包括連續(xù)方程、動量方程、能量方程等。方向的速度分量。連續(xù)方程：V(peff)一0(3)重整的化學(xué)過程用2個獨(dú)立的平衡式來表示：動量方程：(p“w)一一VP+V(／zffV)+Sd．CH4+H2O—CO+3H2(Ah1。。。K一226kJ／mo1)(4)(1)式中：為有效密度，V為混合物的體積平均速度張Co+H2O——CO2+H2(Ah1oo0K一一35kJ／mo1

5、)量；P為壓強(qiáng)；源項Sd．的計算公式l5為：(2)Sdi一(ffV／／~eff)一peffBVlVl(5)*本文受到國家自然科學(xué)基金資助項目(編號：50706004)、遼寧省2007年博士啟動基金的資助。作者簡介：楊國剛，1972年生，副教授，工學(xué)博士；畢業(yè)于大連理工大學(xué)，主要從事氣體流動與傳熱現(xiàn)象的數(shù)值模擬研究；發(fā)表文章1O余篇。地址：(116026)遼寧省大連市凌海路2號。電話：(0411)81702713，13050561150。E-mail：yanggg168@SO—hu．com維普資訊http://www.cqvip.com天然氣工

6、業(yè)2008年3月式中：、floB分別為有效黏度、多孔層的有效滲表1重整反應(yīng)器管道尺寸表cm透率、公式常數(shù)。能量方程表示為：(pfffffvT)一(keffT)+ST(6)ST—RAh+RAh(7)式中：S為與重整反應(yīng)有關(guān)的熱源項；R為水蒸氣重整反應(yīng)與變換反應(yīng)的摩爾速率；Ah為反應(yīng)焓。三、結(jié)果與討論組分守恒方程采用以下的統(tǒng)一形式：圖2為管道中心處CH質(zhì)量濃度沿主流動方(I。。ff)一(I?！癉“)+S≠(8)向與管道高度方向的分布。Y軸單位為相對高度這里，為燃料氣體的質(zhì)量分?jǐn)?shù)。(即實際高度Y與復(fù)合管道總高h(yuǎn)的比值)?？梢钥催吔鐥l件為：出，CH濃

7、度沿著主流動與高度兩個方向逐漸降低。(1)在同壁底部(一0)處沿主流動方向，CH質(zhì)量濃度由人口處的0．195，降U—V—W一0，kf—01一q，J一0至出口處的0．056～0．084；沿高度方向，在人口很短Y的一段距離(約1ram)，反應(yīng)可以發(fā)生在催化劑層很(為CH、H2、CO、H2O、C()2)深的區(qū)域；而在催化劑層的其他區(qū)域，反應(yīng)主要集中(2)在固壁頂部與側(cè)面處在靠近燃料管道的薄層內(nèi)(約1ram)進(jìn)行。U——W一0，q一0，l／一01O(3)在平板中心位置(一0)處8U一一一w一一一0azazazaz采用控制容積積分法對數(shù)學(xué)模型進(jìn)行離散，

8、對流項采用二次迎風(fēng)插值格式，擴(kuò)散項采用中心差分格式。采用SIMPIE算法解決速度與壓力耦合問題。oo10．2二、算例研究xfm、圖2CH的質(zhì)量濃度沿主流動方向分布圖