侯普輝

摘 要：排放法規(guī)的日趨嚴(yán)格，柴油機(jī)尾氣后處理器SCR作為處理NOX排放物的一種有效手段，吸引了越來(lái)越多的目光的關(guān)注。柴油機(jī)尾氣后處理器SCR內(nèi)部的氣體流動(dòng)狀況及尿素噴射效果與其使用效率以及柴油機(jī)的動(dòng)力性和經(jīng)濟(jì)性都有很大的影響。文章通過(guò)建立三維SCR流動(dòng)和噴霧計(jì)算模型，計(jì)算得出SCR內(nèi)部壓力、流場(chǎng)分布以及尿素噴射霧化分布情況。通過(guò)本次仿真，找到內(nèi)部影響SCR工作效率的關(guān)鍵設(shè)計(jì)因素，從而更好的優(yōu)化產(chǎn)品設(shè)計(jì)，滿足性能及一致性的要求。結(jié)果表明：柴油機(jī)尾氣后處理器SCR內(nèi)部孔道的分布對(duì)于性能有一定的影響，需要通過(guò)生產(chǎn)進(jìn)行控制從而滿足尾氣后處理的效率及產(chǎn)品一致性要求。

關(guān)鍵詞：柴油機(jī);SCR;仿真

中圖分類(lèi)號(hào)：U664.1 ?文獻(xiàn)標(biāo)識(shí)碼：B ?文章編號(hào)：1671-7988（2020）05-113-04

Absrtact： With the increasingly strict emission regulations， diesel exhaust processor SCR， as an effective mean to deal with NOX emissions， has attracted more and more attention. The gas flow condition and urea injection effect inside the diesel exhaust processor SCR have great influence on its efficiency and the dynamic and economic of diesel engine.By establishing a three-dimensional SCR flow and spray calculation model， the internal pressure， flow field distribution and urea spray atomization distribution are calculated. Through this simulation， we can find the key design factors that affect SCR's efficiency， so as to better optimize product design and meet the requirements of performance and consistency. The results show that the distribution of the internal hole of the diesel exhaust processor SCR has some effect on performance， and it needs to be controlled by production to meet the efficiency and product consistency requirements of exhaust gas after treatment.

Keyword： Diesel Engine; SCR; Simulation

緒論

隨著排放法規(guī)的日趨嚴(yán)格，柴油機(jī)尾氣后處理SCR作為有效處理柴油機(jī)NOX的技術(shù)手段，吸引了越來(lái)越多目光的關(guān)注。

針對(duì)某發(fā)動(dòng)機(jī)，為了滿足國(guó)五排放的要求，采用了燃燒優(yōu)化+SCR路線：先通過(guò)機(jī)內(nèi)燃燒改善降低微粒排放，再用SCR來(lái)降低NOx排放。該路線對(duì)SCR的轉(zhuǎn)化效率提出了很高的要求，在歐洲采用較多，能夠滿足歐Ⅳ和歐V排放標(biāo)準(zhǔn)。

本論文通過(guò)建立CFD仿真模型并進(jìn)行不同穿孔對(duì)比分析，優(yōu)化SCR系統(tǒng)的設(shè)計(jì)參數(shù)，提高產(chǎn)品性能與效率。

1 項(xiàng)目背景及目的

1.1 選擇性催化還原（SCR）簡(jiǎn)介

選擇性催化還原技術(shù)（SCR）是針對(duì)柴油車(chē)尾氣排放中NOx的一項(xiàng)處理工藝，即在催化劑的作用下，噴入還原劑氨或尿素，把尾氣中的NOx還原成N2和H2O。從而達(dá)到既節(jié)能、又減排的目的，該項(xiàng)技術(shù)是歐洲主流技術(shù)路線，歐洲長(zhǎng)途載貨車(chē)和大型客車(chē)幾乎全部采用這一技術(shù)。SCR中發(fā)生的化學(xué)反應(yīng)如下：

尿素水解：（NH2）2CO+H2O→2NH3+CO2

NOX還原：NO+NO2+2NH3→2N2+3H2O

NH3氧化：4NH3+3O2→2N2+6H2O

在SCR系統(tǒng)中發(fā)生的復(fù)雜的物理和化學(xué)反應(yīng)包括：尿素水溶液的噴射、霧化、蒸發(fā)、尿素的水解和熱解氣相化學(xué)反應(yīng)以及NOX在催化劑表面與NH3發(fā)生的催化表面化學(xué)反應(yīng)。

1.2 項(xiàng)目背景

針對(duì)某型號(hào)柴油發(fā)動(dòng)機(jī)，為滿足國(guó)五排放法規(guī)的要求，初步設(shè)計(jì)出如下圖1所示的筒式SCR后處理器。

采用某品牌的空氣輔助式尿素噴射系統(tǒng)，同時(shí)采用兩塊載體串聯(lián)封裝的方法來(lái)反應(yīng)降低NOx排放。

樣件試制階段，通過(guò)將排氣后處理器（EGP）切割并對(duì)比觀察內(nèi)部結(jié)構(gòu)，發(fā)現(xiàn)EGP內(nèi)部穿孔的分布呈現(xiàn)不同的狀態(tài)，如下圖所示，這里將不同的兩種分布形式定義為a型式和b型式。

造成穿孔型式不用的原因是因?yàn)樵翔F皮在沖孔是模具方向沒(méi)有進(jìn)行控制，從而導(dǎo)致孔的排布產(chǎn)生差異。

為了驗(yàn)證產(chǎn)品性能，同時(shí)也為了對(duì)比分析不同穿孔排布對(duì)整個(gè)后處理系統(tǒng)的影響，本論文中分別建立了SCR系統(tǒng)的內(nèi)部流動(dòng)及尿素噴射模型，通過(guò)對(duì)不同穿孔的排布型式進(jìn)行對(duì)比分析，研究出最優(yōu)化的穿孔型式設(shè)計(jì)，從而保證產(chǎn)品性能及一致性要求。

2 分析模型建立

2.1 控制方程

2.1.1 質(zhì)量守恒方程

流體力學(xué)中，任何流體的運(yùn)動(dòng)都必須遵循質(zhì)量守恒定律，質(zhì)量守恒定律可以描述為在單位時(shí)間里，某一微元體的質(zhì)量增加量與同一時(shí)段內(nèi)流進(jìn)該微元體內(nèi)的凈質(zhì)量相同，由此可以得出連續(xù)方程這一運(yùn)動(dòng)學(xué)方程，質(zhì)量守恒方程的表達(dá)式如下：

2.1.2 動(dòng)量守恒方程

動(dòng)量守恒定律對(duì)流體運(yùn)動(dòng)規(guī)律的數(shù)學(xué)表述用動(dòng)量方程表示，由牛頓第二定律推導(dǎo)得出，動(dòng)量守恒定律可描述為：微元體內(nèi)流體動(dòng)量關(guān)于時(shí)間的變化率等同于外界環(huán)境作用在該微元體上所有力的總和，依據(jù)此定律，可導(dǎo)出x、y和z三個(gè)方向上的動(dòng)量守恒方程：

式中：p為作用于流體微元體上的壓強(qiáng);τxx、τxy、τxz是因分子黏性作用而產(chǎn)生的作用在微元體表面上的黏性應(yīng)力τ的分量;Fx、Fy、Fz為微元體上的體力，x、y、z軸依據(jù)笛卡爾直角坐標(biāo)系建立，同時(shí)確定z軸的正方向豎直向上，與重力方向相反。

2.2 計(jì)算模型及算法

2.2.1 流動(dòng)模型

氣流在SCR內(nèi)的流動(dòng)采用湍流計(jì)算模型，我們稱湍流又為紊流，當(dāng)流體中的雷諾數(shù)超過(guò)臨界值時(shí)，就會(huì)呈現(xiàn)出流動(dòng)特性隨機(jī)變化的無(wú)序的混亂狀態(tài)。本文采取的是當(dāng)前工程中應(yīng)用廣泛的標(biāo)準(zhǔn)k-ε雙方程模型，k表示湍動(dòng)能，ε表示湍動(dòng)能耗散率，計(jì)算公式如下。

2.2.2 噴霧計(jì)算模型

FLUENT中通過(guò)定義顆粒的初始位置、速度、尺寸以及溫度來(lái)使用該模型。根據(jù)對(duì)物理屬性的定義而確定顆粒初始條件，用來(lái)初始化軌道和傳熱傳質(zhì)計(jì)算。

拉格朗日離散相模型遵循歐拉-拉格朗日方法，流體為連續(xù)相，而離散相是計(jì)算流場(chǎng)中粒子的運(yùn)動(dòng)得到。顆粒的平衡方程（x方向）為：

式中：u為流體相速度;up為顆粒速度;為顆粒的單位質(zhì)量曳力;gx為x方向重力加速度;ρp為顆粒密度;ρ1為流體密度;Fx為x方向的其他作用力;Re為相對(duì)雷諾數(shù)（顆粒雷諾數(shù)）;CD為曳力系數(shù);μ為流體動(dòng)力黏度;dp為顆粒直徑;對(duì)于球形顆粒，在一定的雷諾數(shù)范圍內(nèi)，上式中的a1、a2、a3為常數(shù)。

2.3 CFD模型及邊界條件

2.3.1 CFD模型

本文采用Fluent軟件對(duì)柴油機(jī)后處理器進(jìn)行計(jì)算與分析，運(yùn)用Ansys劃分計(jì)算網(wǎng)格，并控制網(wǎng)格質(zhì)量滿足計(jì)算要求。

運(yùn)用離散相模型模擬尿素噴射過(guò)程中微粒的運(yùn)動(dòng)與分布，讓尿素顆粒以一定的粒徑分布和質(zhì)量流量條件從網(wǎng)格入口處進(jìn)入，跟蹤尿素顆粒的運(yùn)動(dòng)軌跡。

2.3.2 邊界條件

對(duì)于流動(dòng)分析，選取壓降較大的ESC工況點(diǎn)C100進(jìn)行計(jì)算。

尿素噴射仿真則在ESC循環(huán)的B50工況點(diǎn)下進(jìn)行。

采用質(zhì)量流量進(jìn)口邊界條件和壓力出口邊界條件，由于柴油機(jī)后處理器尾端和大氣相通，因此設(shè)置其壓力出口邊界條件為0。

SCR內(nèi)部載體按照多孔介質(zhì)模型進(jìn)行計(jì)算，通過(guò)擬合實(shí)驗(yàn)數(shù)據(jù)測(cè)得載體的粘性和慣性阻力系數(shù)如表3所示：

3 計(jì)算結(jié)果分析

3.1 流動(dòng)結(jié)果分析

整個(gè)排氣后處理器系統(tǒng)的壓降如表所示，總壓降為12.35Kpa滿足設(shè)計(jì)的要求（<15Kpa）。

通過(guò)對(duì)比壓降的記過(guò)可以看出，穿孔分布形式的變化對(duì)整個(gè)排氣系統(tǒng)的壓降沒(méi)有影響。系統(tǒng)壓降最大的過(guò)程在載體出口到整個(gè)后處理器系統(tǒng)的出口，達(dá)到了4.08Kpa，說(shuō)明此部分存在設(shè)計(jì)優(yōu)化空間，可以進(jìn)一步的降低壓降，從而提高柴油機(jī)的燃油經(jīng)濟(jì)性。

3.2 尿素水溶液霧化結(jié)果分析

下圖是尿素噴射過(guò)程中尿素分解率的變化曲線，可以看出針對(duì)不同的穿孔分布類(lèi)型，尿素分解率發(fā)生較為明顯的變化，說(shuō)明穿孔的分布設(shè)計(jì)對(duì)SCR系統(tǒng)中尿素的分解效率有較大的影響。

通過(guò)對(duì)比可以看出，a型式穿孔設(shè)計(jì)對(duì)于的尿素分解率更高，達(dá)到了93.2%，與b型式相比尿素分解率提高了大約4%，更高的尿素分解效率更有利于尿素對(duì)尾氣中NOx的去除。因此需要確認(rèn)采用a型式進(jìn)行生產(chǎn)，從而保證更好的轉(zhuǎn)化效率。

3.3 載體端面尿素分布結(jié)果分析

研究表明載體端面上尿素噴射的分布均勻性對(duì)SCR整體的轉(zhuǎn)化效率有明顯的影響，因此在此引入均勻性系數(shù)（Uniformity Index-UI）。

Vi代表NH3的濃度，A是網(wǎng)格面積，i代表網(wǎng)格編號(hào)，Vavg則是NH3的平均濃度，Atot代表所有網(wǎng)格表面積。

UI是一個(gè)介于0-1之間的數(shù)值，接近于1則表示均勻性越高，反之則均勻性越差，此產(chǎn)品開(kāi)發(fā)設(shè)計(jì)要求需要保證UI>0.9。

下圖是載體端面尿素顆粒的分布云圖，通過(guò)計(jì)算均勻性

系數(shù)可以發(fā)現(xiàn)a型穿孔對(duì)應(yīng)的均勻性系數(shù)較高，達(dá)到0.939。

針對(duì)b型穿孔設(shè)計(jì)，UI計(jì)算值為0.898，沒(méi)能滿足設(shè)計(jì)要求（>0.9）。

4 結(jié)論

本論文中分別建立了SCR系統(tǒng)的內(nèi)部流動(dòng)及尿素噴射模型，通過(guò)對(duì)不同穿孔的排布型式進(jìn)行對(duì)比分析，研究結(jié)果表明：

（1）某柴油機(jī)SCR系統(tǒng)壓降及流場(chǎng)分布滿足系統(tǒng)設(shè)計(jì)要求。

（2）a、b兩種不同類(lèi)型穿孔對(duì)SCR內(nèi)部壓力及流場(chǎng)沒(méi)有影響。

（3）a、b兩種類(lèi)型穿孔對(duì)尿素水溶液噴射分解效率影響有明顯差異，a型穿孔分布優(yōu)于b型的穿孔設(shè)計(jì)，尿素分解率提高約4%。

（4）a型設(shè)計(jì)能使載體端面達(dá)到更均勻的NH3分布，UI可以達(dá)到0.939，滿足產(chǎn)品設(shè)計(jì)開(kāi)發(fā)需求。