楊智遠(yuǎn)　江國(guó)和　魏海軍　張旭升　王濤

摘要：

由于船舶機(jī)艙空間狹小，船用二沖程柴油機(jī)SCR反應(yīng)器及其管路布置較為困難，利用CFD對(duì)柴油機(jī)100%負(fù)荷工況下SCR反應(yīng)器前混合段內(nèi)霧化后的尿素溶液與排煙的混合過程進(jìn)行模擬。通過模擬得到6S35MEB9船用柴油機(jī)SCR系統(tǒng)在混合器葉片角度為15°、導(dǎo)流器擴(kuò)展角度為75°時(shí)可以實(shí)現(xiàn)短距離內(nèi)還原劑與排煙的良好混合，且整套系統(tǒng)壓降控制在1 400 Pa以下。對(duì)6S35MEB9船用柴油機(jī)全尺寸SCR系統(tǒng)進(jìn)行100 h臺(tái)架試驗(yàn)，在中國(guó)船級(jí)社NOx排放指南要求的E3推進(jìn)工況下，實(shí)現(xiàn)原機(jī)NOx比排放量由18.15 g/（kW·h）減少至3.17 g/ （kW·h） 。本設(shè)計(jì)為實(shí)船SCR系統(tǒng)的配備提供了理論和應(yīng)用基礎(chǔ)。

關(guān)鍵詞：

船舶低速柴油機(jī)； SCR系統(tǒng)； 混合器； 導(dǎo)流器； 優(yōu)化； 試驗(yàn)驗(yàn)證

中圖分類號(hào)： U664.121.1； TK421.5

文獻(xiàn)標(biāo)志碼： A

Abstract：

In view of the fact that the arrangement of the twostroke diesel engine SCR reactor and its pipe is difficult due to the narrow space in a ship engine room， CFD is used to simulate the mixing process of atomized urea solution and exhaust gas in the mixing section before SCR reactor under 100% load condition of diesel engine. The simulation results show that， when the blade angle of the mixer is 15° and the expansion angle of the diffuser is 75°， the SCR system of 6S35MEB9 marine diesel engine can achieve a good mixing of reducing agent and exhaust smoke in short distance， and the pressure drop of the whole system is controlled less than 1 400 Pa. Through 100 h bench scale test for the fullscale SCR system of 6S35MEB9 marine diesel engine under the E3 propulsion condition required by the NOx emission guidelines of China Classification Society， the NOx brake specific emission of the original engine is reduced from 18.15 g / （kW·h） to 3.17 g / （kW·h）. The design provides a theoretical and practical basis for the real ship SCR system.

Key words：

marine slow speed diesel engine； SCR system； mixer； diffuser； optimization； test verification

0引言

為使船用二沖程低速柴油機(jī)NOx排放量達(dá)到《MARPOL公約》附件VI規(guī)定的排放要求（比排放量限值為3.4 g/（kW·h）），目前最有效的機(jī)外控制技術(shù)之一是尿素選擇性催化還原（urea selective catalytic reduction， UreaSCR）技術(shù)[1]。然而，在實(shí)際船舶應(yīng)用中，船舶柴油機(jī)SCR系統(tǒng)的初次安裝成本和使用成本遠(yuǎn)遠(yuǎn)高于其他廢氣處理技術(shù)[23]。為節(jié)約船舶所有人成本并確保SCR系統(tǒng)在全壽命期內(nèi)高效運(yùn)行，國(guó)內(nèi)外學(xué)者[47]對(duì)SCR反應(yīng)器中催化劑性能、反應(yīng)器結(jié)構(gòu)及其控制策略等進(jìn)行了大量的研究。文獻(xiàn)[810]利用數(shù)值模擬結(jié)合試驗(yàn)驗(yàn)證的方法，研究了SCR系統(tǒng)中尿素溶液的霧化混合性能。文獻(xiàn)[11]在SCR反應(yīng)器前加裝混合器并對(duì)其進(jìn)行了測(cè)試，證明加裝混合器能夠提高SCR系統(tǒng)的NOx轉(zhuǎn)化效率并減少氨氣泄漏。文獻(xiàn)[3]通過建立SCR系統(tǒng)三維動(dòng)態(tài)數(shù)學(xué)模型，對(duì)SCR系統(tǒng)尿素溶液霧化過程和混合過程進(jìn)行了模擬，結(jié)果表明，高壓噴射加靜態(tài)混合器可顯著提高SCR系統(tǒng)的氣液混合均勻度，提高SCR系統(tǒng)脫硝的效率。

本文以6S35MEB9船用二沖程低速柴油機(jī)為研究對(duì)象，利用CFD對(duì)柴油機(jī)100%負(fù)荷工況下SCR反應(yīng)器前混合段內(nèi)霧化后的尿素溶液與排煙的混合過程進(jìn)行模擬，對(duì)混合器、導(dǎo)流器參數(shù)進(jìn)行優(yōu)化后，搭建SCR反應(yīng)系統(tǒng)，隨后進(jìn)行臺(tái)架試驗(yàn)來(lái)驗(yàn)證本設(shè)計(jì)是否滿足tierⅢ的要求。

1數(shù)學(xué)模型和評(píng)價(jià)指標(biāo)

1.1SCR系統(tǒng)

SCR系統(tǒng)主要包括混合段、擴(kuò)張段、SCR反應(yīng)器和收縮段。反應(yīng)器尺寸為4 770 mm×2 048 mm×2 086 mm。催化劑雙排布置，每排催化劑前布置吹灰管路。仿真模型采用六面體和四面體混合網(wǎng)格，網(wǎng)格節(jié)點(diǎn)數(shù)為4 046 614，網(wǎng)格單元數(shù)為14 757 491。SCR反應(yīng)器的網(wǎng)格劃分見圖1。

1.2控制方程

UreaSCR系統(tǒng)的物理化學(xué)反應(yīng)包括系統(tǒng)內(nèi)部的傳熱傳質(zhì)、尿素水溶液蒸發(fā)與分解、催化還原化學(xué)反應(yīng)等，其反應(yīng)機(jī)理十分復(fù)雜[12]。該系統(tǒng)數(shù)學(xué)模型主要包括連續(xù)性方程、動(dòng)量守恒方程、能量守恒方SCR反應(yīng)器網(wǎng)格劃分示意圖

程和湍流方程。這里，湍流方程采

用kεf湍流模型[13]，其k方程、ε方程、ζ方程和f方程分別為

dkdt=Pk-ε+xk（v+Cμζkts）kxk

dεdt=Cε1Pk-Cε2εT+xkv+Cμζktsσεεxk

dζdt=f-ζkPk+xkv+Cμζktsσζζxk

f-L22f2xk2=C1-1+C′1Pkε2/3-ζts+C2Pkk

式中： k為湍動(dòng)能；Pk為應(yīng)力項(xiàng)；ε為湍動(dòng)能耗散率；v為運(yùn)動(dòng)黏度； ζ為黏度尺度比；ts為湍流時(shí)間尺度；T為流體溫度；L為湍流長(zhǎng)度尺度；f為松弛因子；σε和σζ分別為ε方程和ζ方程所對(duì)應(yīng)的普朗特?cái)?shù)；C1、C′1、C2、Cε1、Cε2、Cμ為經(jīng)驗(yàn)常數(shù)。

1.3催化反應(yīng)器多孔介質(zhì)模型

SCR催化劑主要有蜂窩式催化劑和板式催化劑2種。本系統(tǒng)選用蜂窩式催化劑，它屬于多孔介質(zhì)。對(duì)多孔介質(zhì)的仿真主要設(shè)置兩個(gè)阻力系數(shù)，即黏性阻力系數(shù)和慣性阻力系數(shù)。這里，組分傳輸方程、連續(xù)性方程、在慣性坐標(biāo)

系i方向上的動(dòng)量守恒方程[13]分別為

φgρgωm，gt=

φgρgωm，gμgz+φgρgzρgDeffωm，gz+Mm，gl（vl，kl（cm，Ts））

ρt+xi（ρui）=S

t（ρui）+xj（ρuiuj）=-pxi+τijxj+ρgi+Fi

式中：φg為催化劑孔道體積分?jǐn)?shù)；ρg為氣體密度；ωm，g為氣體中組分m的質(zhì)量分?jǐn)?shù)；μg為氣體平均速度；Deff為氣體的平均擴(kuò)散速度；Mm，g為氣體中組分m的摩爾質(zhì)量；l為反應(yīng)l的反應(yīng)速率；vl，k為反應(yīng)l中組分m的當(dāng)量系數(shù)；cm為組分m的濃度；Ts為催化劑溫度；ρ為流體密度；ui和uj分別為流體在i、j方向上的速度分量；S為流體常數(shù)；ρ為流體密度；p為流體在i方向上的壓力；τij為應(yīng)力張量分量；gi為流體在i方向上的重力加速度；Fi為流體在i方向上受的力。

1.4評(píng)價(jià)標(biāo)準(zhǔn)及目標(biāo)

尿素溶液與排煙的混合均勻性對(duì)SCR系統(tǒng)工作性能的影響較為明顯。為對(duì)比各種幾何參數(shù)，如混合器葉片角度、導(dǎo)流器擴(kuò)張角度等因素對(duì)流動(dòng)的影響，須給出流動(dòng)均勻性的評(píng)判準(zhǔn)則，以便對(duì)不同參數(shù)方案的流動(dòng)均勻性進(jìn)行比較。

均勻性指數(shù)[14]是用于評(píng)價(jià)流體在各截面上的濃度和速度均勻分布情況的，用γ表示：

γ=1-Ni=1αi-Xi2X

式中：αi為網(wǎng)格i的流體濃度或速度；為截面上的流體平均濃度或平均速度；Xi為網(wǎng)格i的面積；X為截面面積；N為截面上網(wǎng)格總數(shù)。γ的取值范圍為0～1，γ值越大，截面處的濃度或速度越均勻。

本設(shè)計(jì)確保催化劑入口截面處還原劑質(zhì)量濃度均勻性指數(shù)大于0.85，SCR反應(yīng)器入口截面處尿素溶液與排煙混合組分的速度大小偏差在-10%～10%范圍內(nèi)。

壓降計(jì)算公式為

ΔP=Pin-Pout

式中：ΔP為壓強(qiáng)損失；Pin為進(jìn)口壓強(qiáng)；Pout為出口壓強(qiáng)。

為降低SCR系統(tǒng)本身對(duì)主機(jī)工作性能的影響，本設(shè)計(jì)確保整個(gè)SCR系統(tǒng)壓降低于1 500 Pa。

2混合器、導(dǎo)流器設(shè)計(jì)與優(yōu)化

2.1邊界條件

對(duì)6S35MEB9船用二沖程低速柴油機(jī)在100%負(fù)荷工況下SCR反應(yīng)器進(jìn)行模擬，該工況下的功率為3 250 kW，NOx的比排放量為18.15 g/（kW·h）。排煙的入口流速為21.942 m/s，溫度為225 ℃；質(zhì)量分?jǐn)?shù)為40%的尿素溶液的質(zhì)量流量為69 kg/h；尿素噴嘴有6孔，孔徑為1.4 mm；霧化空氣壓力為2.4 bar（1 bar≈1×105 Pa）；催化劑黏性阻力系數(shù)為2 000 001 m-2，慣性阻力系數(shù)為861 m-1，孔隙率為0.88。

2.2混合器葉片角度優(yōu)化

混合器安裝在SCR反應(yīng)器與尿素噴嘴之間的管路內(nèi)，使混合氣體產(chǎn)生強(qiáng)烈湍流，提高尿素溶液與排煙的混合均勻性?；旌掀饔芍虚g圓轂和8個(gè)帶有角度的長(zhǎng)方形葉片組成，見圖2。

混合器葉片角度α是葉片截面與軸線的夾角。隨著α增大，由混合器導(dǎo)致的SCR系統(tǒng)壓降會(huì)急劇升高，故選擇α為10°、15°、20°、25°進(jìn)行比較研究。

圖3為計(jì)算得到的還原劑質(zhì)量濃度均勻性指數(shù)、尿素溶液與排煙的混合組分速度均勻性指數(shù)和壓降隨α的變化曲線。當(dāng)α為15°時(shí)，這3個(gè)指標(biāo)值分別為0.66、0.92和96 Pa。通過對(duì)圖3的分析可知，α的增加：一方面可以增大排煙湍流強(qiáng)度，有利于混合；另一方面會(huì)引起排煙旋流程度逐漸加劇，徑向速度和切向速度增加，導(dǎo)致速度均勻性降低。因此，綜合考慮均勻性指數(shù)和合理的壓降，選擇α為15°的混合器進(jìn)行下一步優(yōu)化。

2.3導(dǎo)流器設(shè)計(jì)及優(yōu)化

導(dǎo)流器布置在混合器之后，引導(dǎo)混合組分均勻地流向催化劑表面，進(jìn)一步促進(jìn)尿素溶液與排煙的混合。導(dǎo)流器三維結(jié)構(gòu)見圖4。導(dǎo)流器擴(kuò)張角度β（見圖5）是影響混合的關(guān)鍵參數(shù)之一。由于受臺(tái)架空間限制，選擇擴(kuò)張角度β為70°、75°、80°進(jìn)行比較研究。

圖6為還原劑質(zhì)量濃度均勻性指數(shù)、尿素溶液與排煙的混合組分速度均勻性指數(shù)和壓降隨β變化的曲線。當(dāng)β為70°或75°時(shí)，還原劑質(zhì)量濃度均勻性

指數(shù)和尿素溶液與排煙的混合組分速度均勻性指數(shù)相差較小；當(dāng)β為75°時(shí)這2個(gè)均勻性指數(shù)均比β為80°時(shí)的大，同時(shí)引起的壓降較β為80°時(shí)的低36 Pa。因此，在減少反應(yīng)器整體尺寸的前提下，為保證尿素溶液與排煙的混合更加均勻并減小擴(kuò)張段和導(dǎo)流器結(jié)構(gòu)壓降，選擇β為75°。

2.4混合器、導(dǎo)流器和反應(yīng)器綜合分析

根據(jù)流場(chǎng)分析結(jié)論和實(shí)際安裝條件，選定混合器葉片角度為15°，導(dǎo)流器擴(kuò)張角度為75°，混合器位于噴嘴前400 mm，尿素溶液噴嘴與反應(yīng)器入口法蘭之間的距離為3 000 mm，進(jìn)行系統(tǒng)整體混合均勻性分析。

圖7和8分別為在混合器尿素噴嘴上游催化劑入口截面處還原劑質(zhì)量濃度分布、尿素溶液與排煙混合組分速度分布的云圖。計(jì)算結(jié)果表明，在催化劑入口截面處，還原劑質(zhì)量濃度均勻性指數(shù)為0.86，混合組分速度均勻性指數(shù)為0.98，混合器和導(dǎo)流器壓降為195 Pa，滿足設(shè)計(jì)要求。

2.5整體壓降模擬分析

壓降是衡量SCR反應(yīng)器功率損失的重要指標(biāo)。由圖9可知，反應(yīng)器的入口壓強(qiáng)為103 605 Pa，出口壓強(qiáng)為102 160 Pa。含有混合器、導(dǎo)流器、壓縮空氣管路附件的計(jì)算模型的進(jìn)出口壓差為1 445 Pa，小于要求的1 500 Pa的標(biāo)準(zhǔn)，因此整體壓降滿足設(shè)計(jì)要求。

3試驗(yàn)驗(yàn)證

本SCR系統(tǒng)是針對(duì)MAN 6S35MEB9船用二沖程低速發(fā)動(dòng)機(jī)在100%負(fù)荷工況下全尺寸設(shè)計(jì)的，為進(jìn)一步驗(yàn)證SCR系統(tǒng)的反應(yīng)效果，進(jìn)行100 h的耐久性試驗(yàn)，臺(tái)架試驗(yàn)見圖10。

本次試驗(yàn)選用零號(hào)柴油，委托上海滬江柴油機(jī)排放檢測(cè)科技有限公司進(jìn)行排放測(cè)試，采用美國(guó)CAI600測(cè)試設(shè)備，在E3循環(huán)下對(duì)MAN 6S35MEB9發(fā)動(dòng)機(jī)在25%、50%、75%、100%負(fù)荷工況下的SCR系統(tǒng)工作性能分別測(cè)試60 min。測(cè)試涵蓋了柴油機(jī)本身所有的熱工參數(shù)，SCR反應(yīng)器壓降，尿素噴射量，NOx、CO、CO2、CH的排放量，測(cè)試結(jié)果見圖11和12。結(jié)果表明，本SCR系統(tǒng)將原機(jī)的比排放量由18.15 g/（kW·h）降至3.17 g/（kW·h），滿足tier Ⅲ的標(biāo)準(zhǔn)。

4結(jié)論

（1）隨著混合器葉片角度的增加，催化劑入口截面處還原劑質(zhì)量濃度均勻性指數(shù)逐漸增大，但排煙旋流程度逐漸加劇，使混合器入口截面處尿素溶液與排煙混合組分速度均勻性指數(shù)降低，同時(shí)造成較大的壓強(qiáng)損失（即壓降）。為合理控制SCR反應(yīng)器整體壓降，并提高還原劑質(zhì)量濃度和混合組分速度在催化劑入口截面處的均勻性，選擇混合器葉片角度為15°即可獲得較好效果。

（2）在減少反應(yīng)器整體尺寸的前提下，為保證尿素溶液與排煙之間的混合更加均勻并減小擴(kuò)張段和導(dǎo)流器結(jié)構(gòu)壓降，選擇導(dǎo)流器擴(kuò)張角度為75°。

（3）對(duì)MAN 6S35MEB9船用低速發(fā)動(dòng)機(jī)全尺寸SCR系統(tǒng)的臺(tái)架測(cè)試表明，該反應(yīng)器設(shè)計(jì)較為合理，反應(yīng)器的工作效率為75%～85%。在氨逃逸體積分?jǐn)?shù)低于1×10-5的前提下，該設(shè)計(jì)使原機(jī)排放滿足tier Ⅲ的要求，對(duì)其他船用柴油機(jī)的脫硝研發(fā)具有指導(dǎo)意義。

參考文獻(xiàn)：

[1]

蔡永祥. 國(guó)Ⅳ柴油機(jī)OBD系統(tǒng)及其監(jiān)控單元電磁兼容性設(shè)計(jì)[D]. 武漢： 武漢理工大學(xué)， 2010.

[2]BAIK J H， YIM S D， NAM InSik， et al. Control of NOx emissions from diesel engine by selective catalytic reduction （SCR） with urea[J]. Topics in Catalysis， 2004， 30/31： 3741.

[3]朱元清， 周松， 王金玉， 等. 船舶柴油機(jī)UreaSCR系統(tǒng)工作過程計(jì)算研究[J]. 中南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2013， 44（5）： 20932100.

[4]ZHENG Guanyu， PALMER G， SALANTA G， et al. Mixer development for urea SCR applications[R]. New York： SAE， 2009. DOI： 10.4271/2009012879.

[5]RAJADURAI S. Improved NOx reduction using wire mesh thermolysis mixer in urea SCR system[R]. New York： SAE， 2008. DOI： 10.4271/2008012636.

[6]ZHANG Xiaogang， ROMZEK M， MORGAN C. 3D numerical study of mixing characteristics of NH3 in front of SCR[R]. New York： SAE， 2006. DOI： 10.4271/2006013444.

[7]陳悅， 呂林. 船用柴油機(jī)SCR系統(tǒng)混合器仿真與試驗(yàn)研究[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2016， 37（1）： 2429.

[8]BIRKHOLD F， MEINGAST U， WASSERMANN P， et al. Modeling and simulation of the injection of ureawatersolution for automotive SCR DeNOxsystems[J].Applied Catalysis B： Environmental， 2007， 70（1/2/3/4）： 119127.

[9]STRM H， LUNDSTRM A， ANDERSSON B. Choice of ureaspray models in CFD simulations of ureaSCR systems[J]. Chemical Engineering Journal， 2009， 150（1）： 6982.

[10]張文娟， 帥石金， 董紅義， 等. 尿素SCRNOx催化器流動(dòng)、還原劑噴霧及表面化學(xué)反應(yīng)三維數(shù)值模擬[J]. 內(nèi)燃機(jī)學(xué)報(bào)， 2007， 25（5）： 433438.

[11]ZHAN R， LI Wei， EAKLE S T， et al. Development of a novel device to improve urea evaporation， mixing and distribution to enhance SCR performance[R]. New York： SAE， 2010. DOI： 10.4271/2010011185.

[12]MADIA G S. Measures to enhance the NOx conversion in ureaSCR systems for automotive applications[D]. Zurich， Switzerland： Swiss Federal Institute of Technology， 2002.

[13]AVL List GmbH. FIRE v2008Aftertreatment Manual[M]. Graz， Austria： AVL List GmbH， 2008： 6379.

[14]朱天宇， 李德波， 方慶艷， 等. 燃煤鍋爐SCR煙氣脫硝系統(tǒng)流場(chǎng)優(yōu)化的數(shù)值模擬[J]. 動(dòng)力工程學(xué)報(bào)， 2015， 35（6）： 481488.

（編輯趙勉）