秦 巖，劉 宇，李 君，王慶黨，張衛(wèi)國(guó)，張建川

（1.吉林大學(xué)汽車仿真與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室，吉林 長(zhǎng) 春 130022；2.長(zhǎng)城汽車股份有限公司，河北 保 定 071000）

近年來(lái)，我國(guó)的環(huán)境污染程度日益嚴(yán)重，柴油機(jī)微粒物（PM）排放成為城市環(huán)境的主要污染源之一。隨著排放法規(guī)的日益嚴(yán)格，微粒捕集器（DPF）得到廣泛應(yīng)用［1］。對(duì)于DPF，除了捕集效率，載體的壓降水平也是衡量其性能的一個(gè)重要指標(biāo)。本研究從原理上闡明了壁流式蜂窩陶瓷結(jié)構(gòu)的DPF載體壓降產(chǎn)生原理，并利用GT－Power軟件進(jìn)行了模擬計(jì)算，進(jìn)一步對(duì)壓降的產(chǎn)生機(jī)理以及壓降分配進(jìn)行了研究分析。

1 DPF壓降的理論分析

從圖1中可以看出，DPF壓降損失主要分為3大類：氣流收縮和擴(kuò)張損失、入口和出口通道摩擦損失、壁面和炭煙層損失。由于氣流的收縮和擴(kuò)張帶來(lái)的損失較小，而載體入口和出口通道的摩擦損失以及壁面和炭煙層帶來(lái)的壓降占據(jù)了總壓降的絕大部分，所以，在一般的研究中忽略氣體膨脹和收縮帶來(lái)的壓降損失。

對(duì)于載體壁面和炭煙層，根據(jù)達(dá)西定律可以得出壁面的壓降損失Δpwall和炭煙層的壓降損失Δpsoot：

式中：a為入口孔道邊長(zhǎng)；b為出口通道邊長(zhǎng)；w為壁面厚度；ws為炭煙層厚度；μ為廢氣的動(dòng)力黏度；Q為氣體體積流率；L為載體出入口孔道長(zhǎng)度；Df為載體直徑；k0和ksoot分別為壁面和炭煙層的滲透率。

另外，由流體力學(xué)可以得出由入口和出口通道的摩擦損失帶來(lái)的壓降損失Δpinlet和Δpoutlet：

式中：F＝28.454；Vtrap為DPF載體體積。所以，DPF的壓降損失Δp為

結(jié)合式（1）至（4）可知，對(duì)于新鮮的DPF或再生完成的DPF載體，在炭煙初始加載的過(guò)程中，在入口通道壁面上還沒(méi)有出現(xiàn)炭煙層，此時(shí)炭煙層厚度ws為0。對(duì)于物理參數(shù)固定的DPF載體，在發(fā)動(dòng)機(jī)穩(wěn)態(tài)工況下，式（1）Δpwall僅與壁面滲透率k0有關(guān)，式（2）中 Δpsoot為0，式（3）和式（4）中，Δpinlet和Δpoutlet都為定值。所以對(duì)于式（5），在炭煙層形成之前，DPF壓降Δp僅與壁面滲透率k0相關(guān)。隨著炭煙的加載，在載體入口通道壁面逐漸形成炭煙層，此時(shí)DPF壓降還與炭煙層的厚度ws和滲透率ksoot相關(guān)。

Shigeki Daido等對(duì)堇青石和碳化硅材質(zhì)的DPF進(jìn)行了實(shí)驗(yàn)研究［3］，準(zhǔn)確地描述了壁流式DPF的炭煙加載過(guò)程和壓降特性。圖2示出壁面的剖視圖，從圖中可以看出，炭煙的捕集過(guò)程分為3個(gè)階段：①深床過(guò)濾階段（炭煙在壁面積累）；②深床過(guò)濾與餅狀層過(guò)濾階段（炭煙在壁面和炭煙層中積累）；③餅狀層過(guò)濾階段（炭煙在炭煙層中積累）。

捕集過(guò)程中的壓降曲線見(jiàn)圖3，圖中①，②，③依次代表上文中炭煙捕集的3個(gè)階段。在深床過(guò)濾階段，壓降急劇上升，這是因?yàn)榧词购苌俚奶繜熞矔?huì)造成壁中微孔通道變窄，壁面滲透性減小，增大了流動(dòng)阻力；在餅狀層過(guò)濾階段，壓降上升很緩慢，壓降上升速率幾乎為一常數(shù)；而在由深床過(guò)濾向餅狀層過(guò)濾的過(guò)渡階段，壓降上升速率由大變小，直到炭煙層獨(dú)立進(jìn)行炭煙加載。

2 DPF壓降特性仿真分析

2.1 模型的建立與校正

表1列出DPF載體的相關(guān)參數(shù)，圖4示出建立的微粒捕集器GT－Power仿真模型，模型包括DPF、DPF前廢氣進(jìn)入、排氣管道、出入口擴(kuò)張管、DPF后廢氣排出、壓力損失輸出、微粒物沉積量輸出、DPF捕集效率輸出等內(nèi)容。

表1 DPF載體參數(shù)

由式（5）可知，新鮮DPF的壓降僅與載體壁面滲透率有關(guān)。為了提高仿真計(jì)算的準(zhǔn)確性，本研究根據(jù)試驗(yàn)數(shù)據(jù)對(duì)DPF載體的壁面滲透率k0進(jìn)行了校核。選取了3個(gè)工況對(duì)壓降進(jìn)行校核，具體試驗(yàn)數(shù)據(jù)見(jiàn)表2，其中排溫是指DPF入口排氣溫度。

表2 DPF壓降模型校正所用試驗(yàn)數(shù)據(jù)

從壓降公式可以看出，對(duì)于物理參數(shù)確定的新鮮DPF，壓降只與壁面滲透率相關(guān)，而此參數(shù)不能通過(guò)公式運(yùn)算得出。本研究通過(guò)壓降校核，標(biāo)定壁面滲透率數(shù)值，使壓降數(shù)值與試驗(yàn)值相吻合。圖5示出壁面滲透率k0＝9.4×10－14時(shí)模擬值與試驗(yàn)值的對(duì)比。從圖中可知，當(dāng)k0＝9.4×10－14時(shí)，3個(gè)工況DPF壓降的模擬值與試驗(yàn)值的誤差均在5%以內(nèi)，所以k0取9.4×10－14是合理的。

2.2 壓降特性仿真分析

本研究入口排氣流量設(shè)為108kg／h，排氣溫度為563K。為了簡(jiǎn)化計(jì)算，同時(shí)排除其他的干擾因素，排氣成分設(shè)為標(biāo)準(zhǔn)空氣（air）和炭煙（soot）的混合物。其中air的質(zhì)量分?jǐn)?shù)為0.999 8，soot質(zhì)量分?jǐn)?shù)為0.000 2。計(jì)算時(shí)間為3 600s。

圖6示出計(jì)算過(guò)程中DPF壓降隨時(shí)間的變化。從圖中可知，壓降的變化歷程與圖3中的試驗(yàn)壓降曲線是相符的。0～400s為深床過(guò)濾階段；400～900s為深床過(guò)濾與餅狀層過(guò)濾階段；900s以后為餅狀層過(guò)濾階段。在深床過(guò)濾階段，載體壁面上尚未有炭煙層形成，DPF的壓降只與壁面滲透率k0相關(guān)。載體壁面的滲透率隨時(shí)間的變化見(jiàn)圖7，隨著炭煙的捕集，在0～400s壁面的滲透率急劇減小，所以在此階段內(nèi)，DPF的壓降急劇上升。圖8示出炭煙加載過(guò)程中，炭煙層的滲透率隨時(shí)間的變化。從圖中可得，炭煙層的滲透率不隨時(shí)間變化，基本為一定值，所以在400～900s時(shí)，壓降趨緩。

在900s以后，壁面對(duì)炭煙的捕集達(dá)到飽和，此時(shí)炭煙層單獨(dú)起到捕集炭煙的作用，壁面的滲透率趨于一個(gè)穩(wěn)定值。在900s以后，炭煙層的厚度基本呈線性增長(zhǎng)（見(jiàn)圖9）。所以綜合來(lái)看，此時(shí)DPF壓降基本緩慢線性上升［4］。

圖10示出DPF壓降在載體壁面、炭煙層、入口通道和出口通道各部分的分配情況。從圖中可以看出，由DPF載體壁面帶來(lái)的壓降占DPF總壓降的很大一部分。在3 360s左右，炭煙加載為20g左右時(shí)，DPF的總壓降為9.09kPa，壁面壓降為6.88kPa，炭煙層壓降為1.34kPa，入口通道壓降為0.53kPa，出口通道壓降為0.34kPa。

從壁面壓降曲線可以看出，其與圖7的壁面滲透率是相符的，900s以后，由于壁面滲透率保持不變，由式（1）可得，壁面的壓降趨于穩(wěn)定。由于炭煙層的逐漸形成，所以炭煙層厚度ws逐漸增大（見(jiàn)圖9）；另一方面，炭煙層的滲透率保持不變，所以炭煙層壓降趨于線性緩慢增大。由式（3）可得，入口通道壓降的緩慢增大是由于壁面炭煙層厚度逐漸增大所致。由式（4）可得，對(duì)于發(fā)動(dòng)機(jī)穩(wěn)態(tài)工況，出口通道壓降基本為一定值（見(jiàn)圖10）。

3 結(jié)論

a）通過(guò)對(duì)DPF壓降特性進(jìn)行理論分析，闡明了壓降產(chǎn)生的機(jī)理：DPF的壓降主要由可視為多孔介質(zhì)的載體壁面和炭煙層以及載體的入口通道和出口通道的摩擦造成；

b）DPF的炭煙加載過(guò)程分為3個(gè)階段，在深床過(guò)濾階段壓降急劇上升，后逐漸趨緩，直到餅狀層過(guò)濾階段，壓降接近線性增長(zhǎng)，且增長(zhǎng)緩慢；

c）載體壁面滲透率隨炭煙的積累迅速減小，當(dāng)壁面達(dá)到飽和時(shí)，趨于一個(gè)穩(wěn)定值；炭煙層的滲透率保持不變［5］；

d）在DPF的壓降分配方面，炭煙加載到一定程度時(shí)，壁面的壓降占了絕大部分，為80%左右，炭煙層占10%左右，入口通道和出口通道共占10%左右。

［1］ 張 輝.轎車柴油機(jī)微粒捕集器工作過(guò)程數(shù)值模擬及再生控制策略研究［D］.長(zhǎng)春：吉林大學(xué)，2011.

［2］ Yoon C S，Song S H，Chun K M.Measurement of Soot Mass and Pressure Drop Using a Single Channel DPF to Determine Soot Permeability and Density in the Wall Flow Filter［C］.SAE Paper 2007－01－0311.

［3］ Shigeki Daido，Nobuyuki Takagi.Visualization of the PM Deposition and Oxidation Behavior Inside the DPF Wall［C］.SAE Paper 2009－01－1473.

［4］ Xiaogang Zhang，Paul Tennison，William Ruona.3D Numerical Study of Pressure Loss Characteristics and Filtration Efficiency through a Frontal Unplugged DPF［C］.SAE Paper 2010－01－0538.