何 勇，黃雅蘭

（邵陽(yáng)學(xué)院機(jī)械與能源工程學(xué)院，湖南 邵陽(yáng) 422000）

0 前言

柴油機(jī)微粒捕集器簡(jiǎn)稱(chēng)DPF（Diesel Particulate Filter,DPF）。從20世紀(jì)70年代開(kāi)始，國(guó)外學(xué)者針對(duì)柴油機(jī)排氣中的大量微粒對(duì)DPF 的研制開(kāi)發(fā)展開(kāi)研究，到如今DPF是世界上公認(rèn)降低柴油機(jī)PM排放最有效、技術(shù)相對(duì)成熟的后處理裝置之一，性能良好的DPF對(duì)排氣中PM的過(guò)濾效率可達(dá)95%[1]，現(xiàn)在已逐漸成為滿(mǎn)足排放法規(guī)要求的必要技術(shù)配置。在一系列再生技術(shù)中，噴油助燃再生這一主動(dòng)再生技術(shù)，再生燃燒器通過(guò)燃料供給系統(tǒng)和噴油系統(tǒng)適時(shí)噴射一定量的燃油，伴隨供氣系統(tǒng)提供充足的空氣以及點(diǎn)火系統(tǒng)將其點(diǎn)燃，燃燒器便會(huì)均勻產(chǎn)生相應(yīng)的含氧燃?xì)?，引發(fā)過(guò)濾體內(nèi)的微粒燃燒[2]。噴油助燃再生技術(shù)采用與柴油機(jī)一樣的燃料，取用便利，燃料供給系統(tǒng)和噴油系統(tǒng)結(jié)構(gòu)簡(jiǎn)單，再生窗口較寬，燃料燃燒后能量利用率高、響應(yīng)速度快，再生時(shí)間極短且再生效率高[3]。但噴油助燃再生的燃燒過(guò)程與組織具有一定的難度，燃燒器在冷啟動(dòng)時(shí)會(huì)由于燃燒不理想引起二次污染，不僅要讓過(guò)濾體內(nèi)積聚的微?？焖偾彝耆紵乙ＷC過(guò)濾體不因過(guò)熱而損壞，這就需要在油氣比、燃油噴射的霧化質(zhì)量、火焰穩(wěn)定性等方面進(jìn)行優(yōu)化與匹配[4]。

因此，在本課題組前期所設(shè)計(jì)的DPF 噴油助燃再生旋流燃燒器的基礎(chǔ)上，采用數(shù)值模擬技術(shù)，研究四種不同的直葉片安裝角度對(duì)DPF 噴油助燃再生旋流燃燒器流場(chǎng)特性的影響規(guī)律，為進(jìn)一步優(yōu)化設(shè)計(jì)旋流燃燒器提供了依據(jù)。

1 計(jì)算模型介紹

1.1 旋流燃燒器的物理模型

本文中所研究的DPF 噴油助燃再生的旋流燃燒器，其主要由油氣混合室、燃燒室、捕集體段以及總排氣出口段四個(gè)部分組成。所設(shè)計(jì)的旋流燃燒器采用的結(jié)構(gòu)形式為軸流式，在其噴油嘴前端設(shè)計(jì)安裝旋流直葉片來(lái)進(jìn)行供風(fēng)，可以確保旋流燃燒器產(chǎn)生高質(zhì)量的高速旋轉(zhuǎn)的流場(chǎng)，有利于促進(jìn)柴油的霧化與蒸發(fā)以及油氣的充分混合，從而進(jìn)一步保證其在中心軸方向形成相應(yīng)長(zhǎng)度的回流區(qū)來(lái)增強(qiáng)火焰的穩(wěn)定性。

本文所設(shè)計(jì)的軸流式旋流直葉片旋流燃燒器的簡(jiǎn)化模型如圖1 所示，圖中噴油嘴安裝孔半徑R 為12.5 mm，內(nèi)環(huán)半徑R1為31 mm，外環(huán)半徑R2為43 mm，葉片厚度δ 為2 mm，葉片寬度a 為19 mm，由于本文主要對(duì)四種不同的葉片安裝角度θ 這一結(jié)構(gòu)參數(shù)進(jìn)行相關(guān)的研究，因此葉片安裝角度θ 將分別為30°、37°、45°以及50°，直葉片旋流器的長(zhǎng)度t 為20 mm，旋流直葉片數(shù)量n為14個(gè)。

圖1 軸流式旋流直葉片旋流燃燒器的簡(jiǎn)化模型

1.2 流體計(jì)算模型

DPF 在實(shí)際進(jìn)行噴油助燃再生時(shí)，其整體結(jié)構(gòu)、系統(tǒng)具有一定的復(fù)雜性，若完全參照實(shí)際物理模型、參數(shù)來(lái)進(jìn)行數(shù)值仿真計(jì)算，不論是運(yùn)行參數(shù)還是結(jié)構(gòu)參數(shù)在投入實(shí)際工作時(shí)不具有固定性，為了滿(mǎn)足較為寬廣的再生窗口，就需要在相應(yīng)的區(qū)間內(nèi)來(lái)取值，所帶來(lái)的成本和困難性將大幅度提高。因此，有必要對(duì)所研究的模型作進(jìn)一步的假設(shè)及簡(jiǎn)化處理，以滿(mǎn)足工程仿真切實(shí)可行的要求。

對(duì)于假設(shè)及簡(jiǎn)化后的流體計(jì)算模型采用網(wǎng)格劃分軟件作離散處理，因?yàn)楸灸Ｐ椭械牟都w簡(jiǎn)化處理為多孔介質(zhì)模型，在數(shù)值計(jì)算時(shí)會(huì)對(duì)其作參數(shù)設(shè)置，所以需對(duì)捕集體部分與燃燒室、總排氣出口部分作分割處理。出于計(jì)算成本、精度以及流體域的結(jié)構(gòu)等方面的考慮，捕集體部分用六面體結(jié)構(gòu)化網(wǎng)格劃分，其他部分用非結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，總網(wǎng)格數(shù)量約33萬(wàn)，網(wǎng)格質(zhì)量都大于0.3。

1.3 計(jì)算模型與邊界條件

1.3.1 湍流模型

FLUENT 提供的湍流模型常用的有：標(biāo)準(zhǔn)k-ε模型、RNGk-ε模型、Realizablek-ε模型、雷諾應(yīng)力模型（RSM）、k-ε模型、壓力修正k-ε模型以及大漩渦模擬模型[5]。本文中所研究的旋流燃燒器等結(jié)構(gòu)擁有復(fù)雜的流場(chǎng)變化與近壁面流動(dòng)細(xì)節(jié)，這些特點(diǎn)屬于雙方程模型的典型應(yīng)用場(chǎng)景。Realizablek-ε模型滿(mǎn)足對(duì)雷諾應(yīng)力的約束條件，在雷諾應(yīng)力上可保持和實(shí)際湍流情況的相同，這是標(biāo)準(zhǔn)k-ε模型和RNGk-ε模型二者都無(wú)法實(shí)現(xiàn)的，這一優(yōu)勢(shì)在FLUENT 的計(jì)算中能夠更加精確地模擬平面或圓形射流的擴(kuò)散速度，與此同時(shí)在旋轉(zhuǎn)流、帶有方向的壓強(qiáng)梯度的邊界層以及分離流等問(wèn)題計(jì)算中其結(jié)果更加符合實(shí)際情況[6]。因此，本文中的研究將選擇采用更加適合的Realizablek-ε湍流模型。

1.3.2 入口邊界條件

本研究在設(shè)置邊界條件時(shí)選擇速度入口，對(duì)上述選擇的湍流模型需要在設(shè)置中填入湍流強(qiáng)度I與湍流尺度l的數(shù)值，由于當(dāng)前缺少求解湍流動(dòng)能k、湍流耗散率ε的準(zhǔn)確公式，所以需要依據(jù)近似公式（1）和（2）來(lái)進(jìn)行估算。

以上兩式中：uavg代表平均流動(dòng)速度；I代表湍流強(qiáng)度，其可以依據(jù)公式（3）進(jìn)行求解；Cμ代表湍流模型內(nèi)設(shè)定的經(jīng)驗(yàn)常數(shù)，一般取值為0.09；l代表湍流尺度，其可以依據(jù)公式（4）進(jìn)行求解。

以上兩式中：ρ代表流體密度；v代表流體流動(dòng)速率；dH代表入口管道的特征長(zhǎng)度；μ代表空氣的動(dòng)力黏度系數(shù)；L代表流動(dòng)的特征長(zhǎng)度。當(dāng)管道是標(biāo)準(zhǔn)圓形截面，且內(nèi)部流動(dòng)充分發(fā)展，若dH與L相同時(shí)使用其相應(yīng)的水力直徑代替，求解公式如公式（5）所示。

式中：A代表通過(guò)流體的截面面積；lc代表過(guò)流截面上流體潤(rùn)濕周長(zhǎng)，對(duì)應(yīng)到本研究則為氣體入口截面上全部過(guò)流截面的周長(zhǎng)。

本文主要對(duì)旋流燃燒器的流場(chǎng)開(kāi)展研究，考慮到柴油機(jī)的排氣成分具有一定的復(fù)雜性且氣體波動(dòng)較大，故將尾氣簡(jiǎn)化處理為新鮮空氣。選取柴油機(jī)某一工況，并采集這一工況下的尾氣數(shù)據(jù)作為尾氣入口的邊界條件，所采集的數(shù)據(jù)和根據(jù)上述公式求解兩個(gè)入口的湍流強(qiáng)度I與湍流尺度l詳見(jiàn)表1。

表1 入口邊界條件

2 仿真結(jié)果與分析

針對(duì)四種不同葉片安裝角度的再生用軸流式旋流燃燒器模型，選擇合適的湍流模型，根據(jù)對(duì)應(yīng)的入口邊界條件，利用CFD 技術(shù)仿真計(jì)算了其流場(chǎng)，并主要分析了四種不同葉片安裝角度對(duì)中心回流區(qū)、重附著區(qū)的影響。

2.1 四種不同葉片安裝角度對(duì)中心回流區(qū)的影響

由文獻(xiàn)可知，對(duì)于旋流燃燒器性能指標(biāo)的評(píng)價(jià)，一般常用回流區(qū)的大小和旋流強(qiáng)度這兩個(gè)參數(shù)來(lái)進(jìn)行考察，其中，回流區(qū)的大小定義為軸向速度是負(fù)數(shù)的區(qū)域，旋流強(qiáng)度則定義為軸向負(fù)速率的大小[7]。圖2 為四種不同葉片安裝角度下旋流燃燒器Z＝0 平面上速度分布云圖，由圖2 可分析出，四種不同葉片安裝角度下，四種旋流燃燒器的油氣混合室在X 方向中心軸線(xiàn)區(qū)域都能形成一定長(zhǎng)度的中心回流區(qū)，并且該中心回流區(qū)一直延伸至燃燒室前部，與柴油機(jī)尾氣入口處的氣流相遇。此外，從圖2 可明顯地看出，葉片安裝角度分別為45°、50°的兩種旋流燃燒器的中心回流區(qū)明顯大于葉片安裝角度分別為30°、37°的兩種旋流燃燒器的中心回流區(qū)，中心回流區(qū)可有效增強(qiáng)旋流燃燒器內(nèi)油氣的卷吸、混合，對(duì)于提高燃燒的火焰穩(wěn)定性、燃燒的徹底程度及減少氮氧化物的生成量起到積極的作用[8]。

圖2 四種不同葉片安裝角度下旋流燃燒器Z＝0平面上速度云圖

進(jìn)一步分析圖2 可知，隨著葉片安裝角度的增大，旋流燃燒器中心回流區(qū)的回流速度整體上也在逐漸變大；葉片安裝角度不斷增大，其在中心軸線(xiàn)上的回流速度大體上也不斷變大。當(dāng)葉片安裝角度為30°時(shí)，最大回流速度出現(xiàn)在大約138 mm 處，為2.5 m/s；當(dāng)葉片安裝角度為37°時(shí)，最大回流速度出現(xiàn)在大約175 mm 處，為4.2 m/s；當(dāng)葉片安裝角度為45°時(shí)，最大回流速度出現(xiàn)在大約81 mm 處，為3.8 m/s；當(dāng)葉片安裝角度為50°時(shí)，最大回流速度出現(xiàn)在大約113 mm處，為4.0 m/s。葉片安裝角度由30°增大至37°時(shí)，最大回流速度也隨之變大，但葉片安裝角度由37°增大至50°時(shí)，最大回流速度呈現(xiàn)小幅度先下降再上升的趨勢(shì)，整體波動(dòng)較小。

此外，伴隨著葉片安裝角度的增大，旋流燃燒器的中心回流區(qū)長(zhǎng)度也隨之增大，整體上升趨勢(shì)逐漸平緩，對(duì)應(yīng)于四種葉片安裝角度，中心回流區(qū)長(zhǎng)度分別為177.9 mm、194.2 mm、200.7 mm、202.9 mm。合適長(zhǎng)度的中心回流區(qū)長(zhǎng)度不但可以增強(qiáng)旋流燃燒器內(nèi)燃油與新鮮空氣在質(zhì)量、能量及動(dòng)量之間的傳遞[9]，而且也可以讓燃燒室擁有適當(dāng)?shù)膶?duì)流換熱，以進(jìn)一步提高旋流燃燒器的效率。

2.2 四種不同葉片安裝角度對(duì)重附著區(qū)的影響

四種不同葉片安裝角度下，四種旋流燃燒器在燃燒室突擴(kuò)結(jié)構(gòu)處均能形成順時(shí)針旋轉(zhuǎn)型局部回流，這一局部漩渦回流處于燃燒器的突擴(kuò)位置及貼近燃燒器壁面，通常被稱(chēng)作重附著區(qū)。重附著區(qū)的存在可顯著加強(qiáng)燃燒室末端冷熱氣流的充分混合及穩(wěn)定性[10]，從而有利于提高再生時(shí)火焰的穩(wěn)定性。其中，重附著區(qū)長(zhǎng)度則決定這一局部漩渦回流區(qū)的大小，重附著區(qū)長(zhǎng)度隨著葉片安裝角度的增大而逐漸減小。

3 結(jié)束語(yǔ)

四種不同葉片安裝角度下，四種旋流燃燒器的油氣混合室在X 方向中心軸線(xiàn)區(qū)域都能形成一定長(zhǎng)度的中心回流區(qū)，并且該中心回流區(qū)一直延伸至燃燒室前部；葉片安裝角度分別為45°、50°的兩種旋流燃燒器的中心回流區(qū)明顯大于葉片安裝角度分別為30°、37°的兩種旋流燃燒器的中心回流區(qū)；隨著葉片安裝角度的增大，旋流燃燒器的中心回流區(qū)長(zhǎng)度也隨之增大，整體上升趨勢(shì)逐漸平緩，中心回流區(qū)的回流速度整體上也在逐漸變大；葉片安裝角度由30°增大至37°時(shí)，最大回流速度也隨之變大，但葉片安裝角度由37°增大至50°時(shí)，最大回流速度呈現(xiàn)小幅度先下降再上升的趨勢(shì)，整體波動(dòng)較小。四種不同葉片安裝角度下，四種旋流燃燒器在燃燒室突擴(kuò)結(jié)構(gòu)處均能形成重附著區(qū)，重附著區(qū)長(zhǎng)度隨著葉片安裝角度的增大而逐漸減小。