施盛耀,湯 東,李天祥

(江蘇大學(xué) 汽車與交通工程系, 江蘇 鎮(zhèn)江 212013)

0 引言

近幾十年來(lái),我國(guó)的經(jīng)濟(jì)迅猛發(fā)展,柴油機(jī)也獲得了越來(lái)越廣泛的應(yīng)用[1]。但柴油機(jī)排放的顆粒物(PM)作為空氣污染源之一,致使空氣污染問(wèn)題愈發(fā)嚴(yán)峻[2-3]。為此世界各國(guó)制定了嚴(yán)格的排放法規(guī),并強(qiáng)制要求柴油機(jī)安裝顆粒物捕集器(DPF),以減少顆粒物的排放[4-5]。柴油機(jī)排氣環(huán)境極其嚴(yán)峻復(fù)雜,DPF在使用過(guò)程中容易出現(xiàn)破損和老化等失效情況,所以需要使用傳感器對(duì)DPF進(jìn)行監(jiān)測(cè),以保證其正常工作。隨著國(guó)家排放法規(guī)的不斷加嚴(yán),傳統(tǒng)的壓差傳感器響應(yīng)時(shí)間長(zhǎng)、精度不高,無(wú)法滿足診斷要求,因此,需要進(jìn)一步開(kāi)發(fā)車用柴油機(jī)顆粒物傳感器實(shí)現(xiàn)實(shí)時(shí)、快速、精確診斷[6-8]。

顆粒物傳感器根據(jù)測(cè)量原理主要分為光學(xué)、化學(xué)和電學(xué)傳感器。Soldo等[9]基于粒子的尺寸大于散射輻射的波長(zhǎng),粒子的表面反射入射輻射,折射和衍射輻射被檢測(cè)和計(jì)數(shù)的原理開(kāi)發(fā)出了M-DUST顆粒物傳感器,實(shí)現(xiàn)了對(duì)顆粒物的實(shí)時(shí)監(jiān)測(cè)。Shen等[10]對(duì)電化學(xué)技術(shù)進(jìn)行了深入的研究,開(kāi)發(fā)了電流式與電位式傳感器。Hagen等[11]通過(guò)模擬研究發(fā)現(xiàn)顆粒物濃度的增加會(huì)導(dǎo)致傳感器電容值的增大,且顆粒物濃度增加得越快,傳感器電容值上升的幅度也越大,對(duì)電容式顆粒物傳感器進(jìn)行了改良。隨著國(guó)六排放法規(guī)(表1)的發(fā)布,國(guó)內(nèi)專家也開(kāi)展了對(duì)顆粒物傳感器的研究。楊永杰等[12]以光敏區(qū)光亮為媒介,建立顆粒物濃度與電信號(hào)的關(guān)系,發(fā)明了一種光電傳感器。孫吉勇等[13]基于光散射法提出了一種主要結(jié)構(gòu)為凸透鏡的顆粒物傳感器。

表1 輕型車排放標(biāo)準(zhǔn)中顆粒物的排放限值

漏電流式顆粒物傳感器應(yīng)運(yùn)而生,其生產(chǎn)成本低且能對(duì)廢氣中的顆粒物濃度進(jìn)行瞬態(tài)監(jiān)測(cè)[14-15]。如圖1所示,傳感器頭部豎直安裝在DPF后方的排氣管中,基于文丘里管原理使尾氣流入傳感器中。由于中央電極電壓高達(dá)800～1 500 V,在進(jìn)入柱狀中央電極與筒狀接地外殼間形成的濃度測(cè)試區(qū)時(shí),尾氣中的顆粒物被電離成導(dǎo)體,產(chǎn)生漏電流。建立漏電流值和顆粒物濃度之間的關(guān)系即可實(shí)時(shí)測(cè)量尾氣中的顆粒物濃度。

圖1 漏電流式顆粒物傳感器結(jié)構(gòu)示意圖

1 建模方法

1.1 物理場(chǎng)模型

1.1.1流場(chǎng)模型

顆粒物的運(yùn)動(dòng)規(guī)律主要受到其周圍尾氣流動(dòng)規(guī)律的影響,故將其簡(jiǎn)化為湍流運(yùn)動(dòng)。由此可以采用質(zhì)量與動(dòng)能守恒方程來(lái)表示傳感器內(nèi)氣體的運(yùn)動(dòng)規(guī)律。

質(zhì)量守恒方程：

(1)

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

(2)

式中：ρ為流體密度,g/m3;ui=(u1,u2,u3)為平均速度分量;xi=(x1,x2,x3)為笛卡爾坐標(biāo);p為平均壓力,Pa。

用來(lái)計(jì)算湍流運(yùn)動(dòng)的方法有很多種,考慮到計(jì)算效率和求解精度的要求,選擇k-ε湍流模型對(duì)排氣管以及傳感器中的流場(chǎng)進(jìn)行計(jì)算,用來(lái)研究流場(chǎng)分布情況,計(jì)算方程如下：

(3)

(4)

方程(4)中Gk是平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng),通過(guò)式(5)計(jì)算得到[14]：

(5)

式中：k為湍動(dòng)能,J;μt為湍流黏性系數(shù),Pa·s;σk、σε為k方程與ε方程的湍流普朗特?cái)?shù)。

1.1.2穩(wěn)態(tài)靜電場(chǎng)模型

在傳感器的內(nèi)部,通過(guò)帶有1 000 V中央電極與接地外殼的共同作用,在兩者之間的空隙產(chǎn)生一個(gè)帶有穩(wěn)定電場(chǎng)的濃度測(cè)試區(qū)。故選擇穩(wěn)態(tài)靜電場(chǎng)模型施加對(duì)傳感器內(nèi)部氣流運(yùn)動(dòng)的影響。電場(chǎng)強(qiáng)度公式為：

(6)

1.1.3溫度場(chǎng)模型

傳感器內(nèi)溫度場(chǎng)主要是由材料的導(dǎo)熱性能決定的,因此溫度場(chǎng)模型的建立必須考慮到外界環(huán)境溫度、顆粒物進(jìn)氣溫度及各材料的導(dǎo)熱性能等主要方面。傳感器結(jié)構(gòu)如圖2所示,各零件的材料屬性參數(shù)如表2。在直角坐標(biāo)下各向異性的三維導(dǎo)熱微分方程如下：

圖2 顆粒物傳感器結(jié)構(gòu)示意圖

表2 顆粒物傳感器零件材料屬性參數(shù)

(7)

式中：T為瞬態(tài)溫度,K;c為比熱容,J/(kg·K)[15];ρ為質(zhì)量密度,kg/m3;k為導(dǎo)熱系數(shù),W/(m·K);U為系統(tǒng)內(nèi)能,J。

1.2 幾何模型與邊界條件

考慮到顆粒物通過(guò)文丘里管頭部進(jìn)入傳感器內(nèi)部的過(guò)程中,排氣管的形狀結(jié)構(gòu)和傳感器的安裝位置都會(huì)對(duì)流場(chǎng)的分布情況產(chǎn)生影響,故建立傳感器與排氣管的連體模型,如圖3所示。同時(shí)因?yàn)榕艢夤芘c顆粒物傳感器都可以視為軸對(duì)稱物體,為了減少仿真模擬的計(jì)算量可以選擇只建立一半的模型。在模擬流場(chǎng)時(shí),由于排氣管氣體流速較快、停留時(shí)間較短,故可以認(rèn)為排氣管出入口與傳感器入口溫度一致。為進(jìn)一步縮短仿真時(shí)間,可建立如圖4所示的傳感器二維旋轉(zhuǎn)軸模型。

圖3 顆粒物傳感器三維模型示意圖

圖4 顆粒物傳感器二維旋轉(zhuǎn)軸模型示意圖

在模擬軟件COMSOL中,設(shè)置流場(chǎng)時(shí)入口采用速度入口邊界,出口采用壓力出口邊界,壁條件選擇無(wú)滑移。設(shè)置溫度場(chǎng)時(shí),排氣管入口處與傳感器沉沒(méi)在排氣管內(nèi)的外殼設(shè)置為廢氣溫度。裸露在外的傳感器外殼則設(shè)置為大氣溫度,出口設(shè)置為流出邊界。模擬所設(shè)置的排氣流速與排氣溫度參考493柴油機(jī)的怠速至全負(fù)荷工況,其中排氣流速取9.1～32.7 m/s,排氣溫度取100～550 ℃。

2 模擬分析結(jié)果

2.1 廢氣溫度對(duì)顆粒物沉積速度的影響分析

由于模擬時(shí)考慮符合普遍情況,故使顆粒物統(tǒng)一默認(rèn)為100 nm粒徑,顆粒物濃度設(shè)定為1 mg/m3,顆粒物所帶電荷量考慮為攜帶3個(gè)單位的正電荷。為覆蓋多數(shù)應(yīng)用場(chǎng)景,選定排氣溫度為373～823 K,選擇步長(zhǎng)50 K取工況點(diǎn),模擬不同廢氣溫度下傳感器濃度測(cè)試區(qū)內(nèi)的顆粒物沉積速度變化關(guān)系。

圖5—圖7展示了顆粒物沉積速度隨廢氣溫度的變化關(guān)系。從圖中可以看出,濃度測(cè)試區(qū)內(nèi)顆粒物的3種沉積速度均與排氣溫度呈正相關(guān)。因?yàn)樵跉怏w壓強(qiáng)不變的情況下,廢氣分子的平均自由程λ會(huì)隨著排氣溫度的升高而增加,坎寧安校正因子CC也隨之明顯增大。經(jīng)過(guò)研究可知,顆粒物沉積速度主要是由熱泳系數(shù)、傳感器入口處的溫度梯度決定的。同時(shí)排氣溫度導(dǎo)致的廢氣分子的平均自由程λ增大,致使克努森數(shù)Kn呈明顯增加趨勢(shì),熱泳系數(shù)Kth也隨之緩慢上升。因此廢氣溫度升高直接決定顆粒物熱泳沉積速度的增大。顆粒物的電泳沉積速度主要由電遷移率與電場(chǎng)強(qiáng)度2個(gè)部分組成,而傳感器內(nèi)部可視為穩(wěn)態(tài)電場(chǎng),故電場(chǎng)強(qiáng)度保持不變。當(dāng)排氣溫度增大時(shí),坎寧安校正因子的增幅要高于氣體動(dòng)態(tài)黏度,故顆粒物的電遷移率和電泳沉積速度也隨之增大。顆粒物的慣性碰撞沉積速度主要由坎寧安校正因子和氣體動(dòng)態(tài)黏度的比值決定,由于坎寧安校正因子隨溫度上升的增加幅度略高于氣體動(dòng)態(tài)黏度,故顆粒物的慣性碰撞沉積速度隨廢氣溫度的升高呈現(xiàn)緩慢上升的趨勢(shì)。

圖6 電泳沉積速度隨溫度的變化關(guān)系曲線

圖7 慣性碰撞沉積速度隨溫度的變化關(guān)系曲線

2.2 廢氣流速對(duì)顆粒物沉積速度的影響分析

顆粒物條件設(shè)置如2.1小節(jié),廢氣的流速分別設(shè)置為2.7、4、5、6、7 m/s,模擬不同廢氣流速下傳感器濃度測(cè)試區(qū)內(nèi)的顆粒物沉積速度變化關(guān)系。

圖8—圖10展示了顆粒物沉積速度隨廢氣流速的變化關(guān)系。

圖8 熱泳沉積速度隨流速的變化關(guān)系曲線

圖9 電泳沉積速度隨流速的變化關(guān)系曲線

圖10 慣性碰撞沉積速度隨流速的變化關(guān)系曲線

從圖中可以看出,當(dāng)流入傳感器的廢氣流速不斷增大時(shí),顆粒物熱泳沉積速度呈下降趨勢(shì),而電泳沉積速度和慣性碰撞沉積速度呈現(xiàn)增加的趨勢(shì)。對(duì)于熱泳沉積來(lái)說(shuō),當(dāng)傳感器入口廢氣流速升高時(shí),廢氣在濃度測(cè)試區(qū)內(nèi)溫度損失降低,濃度測(cè)試區(qū)內(nèi)的溫度梯度就會(huì)變小,而廢氣本身的溫度與熱泳系數(shù)變化不大,故當(dāng)傳感器入口廢氣流速增大時(shí),顆粒物熱泳沉積速度會(huì)減小。對(duì)于電泳沉積來(lái)說(shuō),當(dāng)傳感器內(nèi)廢氣流速較低時(shí),廢氣從傳感器入口運(yùn)動(dòng)到濃度測(cè)試區(qū)的過(guò)程中會(huì)與傳感器壁面進(jìn)行熱交換,流速較高時(shí)廢氣運(yùn)動(dòng)到濃度測(cè)試區(qū)導(dǎo)致流速較高時(shí)濃度測(cè)試區(qū)的廢氣溫度也較高。氣體流速的增加會(huì)導(dǎo)致氣體溫度增加,故當(dāng)流速增大時(shí),電泳沉積速度呈緩慢增加趨勢(shì)。對(duì)于慣性碰撞沉積來(lái)說(shuō),廢氣在傳感器壁面垂直方向上的速度分量與坎寧安校正因子兩者與顆粒物的慣性碰撞沉積速度呈正相關(guān)。流速的增大會(huì)導(dǎo)致廢氣溫度緩慢增加,從而導(dǎo)致坎寧安校正因子CC緩慢增大。同時(shí),廢氣流速的增大會(huì)使廢氣與傳感器壁面垂直的速度分量不斷增大,故當(dāng)廢氣流速增大時(shí),顆粒物的慣性碰撞沉積速度迅速增大,且流速越大顆粒物的慣性碰撞沉積速度增加幅度越大。

3 傳感器顆粒物沉積試驗(yàn)

為進(jìn)一步研究柴油機(jī)工作過(guò)程中顆粒物在漏電流式顆粒物傳感器內(nèi)的沉積速度變化情況及其對(duì)傳感器信號(hào)輸出值的影響,搭建柴油機(jī)實(shí)驗(yàn)臺(tái)架,對(duì)不同工況下的傳感器信號(hào)變化規(guī)律進(jìn)行分析,結(jié)合沉積機(jī)理與傳感器信號(hào)對(duì)傳感器響應(yīng)時(shí)間進(jìn)行預(yù)測(cè),并對(duì)模擬結(jié)果進(jìn)行驗(yàn)證。

3.1 試驗(yàn)方案

建立柴油機(jī)試驗(yàn)臺(tái)架,通過(guò)測(cè)功機(jī)對(duì)柴油機(jī)輸出功率進(jìn)行測(cè)量,在排氣管上安裝漏電流式顆粒物傳感器,測(cè)量柴油機(jī)不同工況下的排氣濃度并記錄傳感器的輸出信號(hào)值,利用溫度傳感器對(duì)傳感器附近的廢氣溫度進(jìn)行測(cè)量,通過(guò)流量計(jì)測(cè)量排氣管內(nèi)排氣流量,采用AVL415煙度計(jì)對(duì)傳感器輸出信號(hào)進(jìn)行對(duì)比驗(yàn)證。選取不同負(fù)荷不同轉(zhuǎn)速的多個(gè)工況點(diǎn),測(cè)量傳感器對(duì)應(yīng)工況點(diǎn)下的輸出信號(hào)失效時(shí)間。對(duì)每一個(gè)傳感器在每一工況下的信號(hào)輸出進(jìn)行4次重復(fù)采集,計(jì)算4次的平均值進(jìn)行數(shù)據(jù)分析,避免偶然誤差。試驗(yàn)選用的工況點(diǎn)及相應(yīng)轉(zhuǎn)速、負(fù)荷情況如表3所示。

表3 顆粒物傳感器信號(hào)試驗(yàn)的工況點(diǎn)

3.2 試驗(yàn)設(shè)備

選取一臺(tái)四缸四沖程、高壓共軌、增壓中冷493柴油機(jī)搭建試驗(yàn)臺(tái)架,柴油機(jī)主要參數(shù)如表4所示,測(cè)功機(jī)選用杭州中成測(cè)試設(shè)備有限公司的電渦流測(cè)功機(jī),在排氣管尾端安裝AVL415煙度計(jì),保證傳感器的信號(hào)輸出規(guī)律與排氣煙度的變化規(guī)律一致。圖11展示了發(fā)動(dòng)機(jī)的試驗(yàn)臺(tái)架,圖12展示了試驗(yàn)的信號(hào)采集系統(tǒng)。

表4 493柴油機(jī)的主要參數(shù)

3.3 試驗(yàn)結(jié)果分析

圖13展示了不同工況下傳感器失效時(shí)間分布。從圖中可以看出,對(duì)于相同轉(zhuǎn)速下的傳感器,其失效時(shí)間隨著負(fù)荷的增大呈現(xiàn)出急速減小的趨勢(shì),這是由于傳感器失效時(shí)間主要受到顆粒物濃度以及顆粒物沉積速度的影響,當(dāng)柴油機(jī)轉(zhuǎn)速不變,其充量系數(shù)基本保持不變,此時(shí)若增加負(fù)荷,循環(huán)供油量就要增加,導(dǎo)致過(guò)量空氣系數(shù)減小,燃燒室內(nèi)缺氧的區(qū)域擴(kuò)大,混合氣的形成與燃燒就會(huì)變得惡化,從而使顆粒物的濃度增加,而顆粒物濃度的增大則會(huì)導(dǎo)致失效時(shí)間不斷減少。當(dāng)顆粒物濃度相差不大時(shí),溫度較高的工況傳感器失效時(shí)間也較短,這是由于溫度的增加會(huì)導(dǎo)致顆粒物3種沉積機(jī)制所導(dǎo)致的沉積速度均呈增加趨勢(shì),而失效時(shí)間與顆粒物濃度和沉積速度呈負(fù)相關(guān),故溫度較高時(shí)傳感器失效時(shí)間較短。前文對(duì)不同排氣狀態(tài)下顆粒物沉積速度進(jìn)行了模擬,為對(duì)模擬結(jié)果進(jìn)行驗(yàn)證,基于傳感器失效時(shí)間與顆粒物濃度和顆粒物沉積速度之間的關(guān)系,將對(duì)應(yīng)工況下試驗(yàn)所獲得的傳感器失效時(shí)間轉(zhuǎn)換成顆粒物沉積速度,并與相同排氣條件下模擬獲得的顆粒物沉積速度進(jìn)行對(duì)比分析。

圖13 不同工況下傳感器失效時(shí)間分布圖

圖14展示了轉(zhuǎn)速為2 000 r/min時(shí)不同負(fù)荷下模擬與試驗(yàn)顆粒物沉積速度分布規(guī)律。從圖中可知,模擬與試驗(yàn)顆粒物沉積速度規(guī)律相似,隨著發(fā)動(dòng)機(jī)負(fù)荷的增大,模擬和試驗(yàn)顆粒物沉積速度均呈現(xiàn)出不斷增大的趨勢(shì),這是由于隨著負(fù)荷的增大,柴油機(jī)循環(huán)噴油量增加,燃燒更加劇烈,放出的熱量也更多,導(dǎo)致排出廢氣的溫度也不斷升高。由前文可知,廢氣溫度的升高會(huì)導(dǎo)致顆粒物沉積速度增大,故隨著負(fù)荷的增大,模擬和試驗(yàn)顆粒物沉積速度也呈現(xiàn)增大的趨勢(shì)。同時(shí),在低負(fù)荷時(shí),試驗(yàn)所得的顆粒物沉積速度相較于模擬獲得的顆粒物沉積速度有一個(gè)較為明顯的降低,這是由于隨著失效時(shí)間的增加,粒子重返氣流的概率增加。

圖14 模擬與試驗(yàn)顆粒物沉積速度與負(fù)荷率的關(guān)系曲線

4 結(jié)論

1) 氣體溫度升高時(shí),濃度測(cè)試區(qū)內(nèi)溫度梯度增大,顆粒物熱泳沉積速度不斷增大。隨著廢氣流速的上升,廢氣在濃度測(cè)試區(qū)內(nèi)停留的時(shí)間縮短,溫度損失降低,濃度測(cè)試區(qū)內(nèi)的溫度梯度變小,顆粒物熱泳沉積速度不斷減小。

2) 隨著溫度的增加,電遷移率Z不斷增大,致使顆粒物電泳沉積速度加快。流速較高時(shí)廢氣運(yùn)動(dòng)到濃度測(cè)試區(qū)的時(shí)間減少,溫度下降也會(huì)減少,故顆粒物電泳沉積速度增加。

3) 氣體溫度增加時(shí),坎寧安校正因子增幅大于氣體動(dòng)態(tài)黏度增大,顆粒物的慣性碰撞沉積速度隨廢氣溫度的升高呈現(xiàn)緩慢上升的趨勢(shì)。廢氣流速的增大使廢氣與傳感器壁面垂直的速度分量不斷增大,導(dǎo)致顆粒物的慣性碰撞沉積速度迅速增大。

4) 相同轉(zhuǎn)速下,顆粒物濃度隨著負(fù)荷的增加而增大,傳感器的失效時(shí)間隨著負(fù)荷的增大而減小。隨著發(fā)動(dòng)機(jī)負(fù)荷的增大,模擬和試驗(yàn)顆粒物沉積速度均呈現(xiàn)不斷增大的趨勢(shì),在低負(fù)荷時(shí),試驗(yàn)所得的顆粒物沉積速度較模擬獲得的顆粒物沉積速度明顯降低。