張應(yīng)兵，姚煒

(安徽江淮汽車股份有限公司技術(shù)中心，安徽合肥 230022)

?

基于缸內(nèi)過程模擬的發(fā)動機優(yōu)化設(shè)計

張應(yīng)兵，姚煒

(安徽江淮汽車股份有限公司技術(shù)中心，安徽合肥 230022)

針對某GDI發(fā)動機的多種設(shè)計方案，利用CFD手段模擬了缸內(nèi)的流動、噴霧、點火燃燒等過程，對流動和燃燒情況進行評估，給出了設(shè)計方案的優(yōu)化選擇建議。

GDI發(fā)動機；缸內(nèi)過程模擬；優(yōu)化設(shè)計

0 引言

汽油直接噴射(GDI)技術(shù)是汽油機發(fā)展的主流方向之一，相對于傳統(tǒng)的氣道噴射發(fā)動機，它更加環(huán)保和高效。但缸內(nèi)直接噴射也要求更加匹配的缸內(nèi)氣體流動和燃燒過程，對發(fā)動機的設(shè)計工作提出了更高的要求。

在文中提及的GDI汽油機設(shè)計過程中，工程師提出了4個燃燒室設(shè)計方案，如何選定最優(yōu)的方案成為困擾項目組的難題。最終作者利用缸內(nèi)過程的CFD模擬直觀了解氣體流動、噴霧和點火燃燒的過程，評價各設(shè)計方案的特點，從而達到最優(yōu)化設(shè)計的效果。

1 設(shè)計方案介紹

圖1 氣道方案示意圖

設(shè)計工程師提供了如圖1所示的兩種氣道設(shè)計方案A和B，其中方案B向氣門桿部的突起更大，希望帶來較大的滾流比。

另外預(yù)計的壓縮比變動區(qū)間為9.7～10.2，因此對表1所示的4個設(shè)計方案進行了建模分析。

表1 設(shè)計方案列表

2 CAE建模

2.1 三維模型與網(wǎng)格

為了建立網(wǎng)格模型，作者準備了完整的燃燒室表面文件，包含進排氣道、進排氣門燃燒室頭部、火花塞、缸套和活塞頭部等位置。此次分析采用了AVL的CFD軟件FIRE，網(wǎng)格主尺寸為1 mm，火花塞部位進行了適當(dāng)加密。

圖2和表2分別是此次分析的氣門升程曲線和氣門正時的相關(guān)數(shù)據(jù)。因為網(wǎng)格劃分的需要，此次分析中最小氣門升程被定義為0.2 mm。已有的經(jīng)驗表明當(dāng)氣門升程小于0.2 mm時，缸內(nèi)氣體的流動狀態(tài)受氣道影響極小。最終的計算網(wǎng)格是以六面體為主的混合網(wǎng)格，數(shù)目在(70～150)萬之間。

圖2 氣門升程曲線圖

表2 氣門正時

2.2 邊界條件和計算設(shè)置

為了使分析結(jié)果具有可比性，此次模擬進出口均使用了靜壓邊界，且設(shè)定排氣道出口為0.36 MPa，溫度為830 ℃；進氣道入口的壓力為0.21 MPa，溫度為40 ℃；循環(huán)噴油量均為71.5 mg。

計算中將流動介質(zhì)設(shè)置為可壓縮模式，并選用了k-ε湍流方程以及標準壁面方程。動量的差分格式為MINMOD，湍流、能量和組分的差分格式則選用了一階Upwind，連續(xù)性方程選用中心差分格式。

2.3 噴霧及燃燒模型

噴霧和燃燒子模型的選取是缸內(nèi)燃燒模擬最重要的部分。

此次計算的噴霧模擬中選用了Dukowicz霧化模型以及Wave破碎模型；壁面油膜是模擬的重點關(guān)注項，選用了FIRE自帶的Bai Gosman/Expert模型。

燃燒模型則選用了ECFM模型。這是擬序火焰模型的一種，由標準CFM模型和MCFM模型擴展得到，使用兩步化學(xué)反應(yīng)機制，可以用于模擬火花點火發(fā)動機的預(yù)混燃燒、敲缸和排放物生成。

3 計算結(jié)果分析

3.1 排氣流動

圖3顯示了不同曲軸轉(zhuǎn)角時排氣道的速度場分布，曲軸轉(zhuǎn)角170°、250°和340°分別對應(yīng)排氣初始、最大氣門升程和排氣末期3個不同時刻。

可以看到排氣初期由于氣缸和排氣道壓力差較大，排氣流速較大；隨著排氣門打開，壓力差減小，氣道內(nèi)的流速逐漸降低，并出現(xiàn)明顯的流動分離情況；而在排氣末期，由于喉口面積減小，缸內(nèi)氣體流出速度逐漸增大，氣道內(nèi)的流動分離也逐漸減小。

圖3 排氣道流速分布圖

3.2 進氣流動

對于缸內(nèi)直噴的發(fā)動機，進氣滾流是設(shè)計開發(fā)中需要重點關(guān)注的指標，作者通過滾流的旋轉(zhuǎn)角速度來評價缸內(nèi)滾流，將其除以曲軸的旋轉(zhuǎn)角速度以得到一個量綱為一的滾流比。圖4和圖5分別顯示了4個方案整個計算循環(huán)中滾流比和ω滾流比的變化曲線。由于每個方案的ω滾流比呈對稱分布，因此圖5的ω滾流比只取了其中一條曲線。

從圖4—5可以看到滾流比受氣道影響較大，優(yōu)化后的氣道滾流比相對原氣道有了較大的提升，最大的瞬態(tài)滾流比從2.6提高到3.4左右；但滾流比曲線的走勢沒有改變，滾流比極值時刻基本相同。壓縮比會對瞬態(tài)滾流比造成一定影響，但影響不大。而ω滾流比沒有明顯的變化規(guī)律，顯然它受氣道和壓縮比共同影響。

需要注意的是CFD模擬得到的滾流數(shù)值不能直接與穩(wěn)態(tài)試驗臺的數(shù)值進行比較。

圖4 瞬態(tài)滾流比曲線圖

圖5 瞬態(tài)ω滾流比曲線圖

3.3 進氣量

除了滾流比外，汽油機的進氣量也是作者考慮的重點，圖6顯示了不同方案氣缸內(nèi)物質(zhì)質(zhì)量隨曲軸轉(zhuǎn)角的變化趨勢。由于此次分析中不同方案的噴油量相同，因此該曲線可以看作進氣量的變化曲線。

圖6 缸內(nèi)質(zhì)量變化曲線圖

可以非常明顯地看到：氣道對進氣量的影響很大，氣道A的進氣量明顯大于氣道B；而壓縮比也會對進氣量造成一定的影響，增大壓縮比后缸內(nèi)滾流增大，使得氣體回流減少，因此缸內(nèi)質(zhì)量會略有增加。

3.4 點火預(yù)測

汽油機點火與火花塞能量和點火時刻的局部流動密切相關(guān)。為了保證點火，火花塞必須具有一定的氣流速度和湍動能。圖7和圖8分別顯示了不同設(shè)計方案點火時刻火花塞附近的速度和湍動能分布，可以明顯看到方案一和方案三相似，而方案二和方案四相似，這說明氣道對流速和湍動能的影響要大于壓縮比。方案三的流速略大于方案一，說明增大壓縮比能夠在一定程度上增大火花塞附件的氣流速度。

增大湍動能有利于點火，因此相對方案一、三，方案二、四會更加容易著火。

圖7 點火時刻流速分布圖

圖8 點火時刻湍動能分布圖

3.5 缸內(nèi)壓力和溫度

缸內(nèi)的壓力和溫度是兩個非常重要和直觀的參數(shù)，它們直接影響了發(fā)動機的動力性、經(jīng)濟型和可靠性，圖9和圖10分別顯示了4個設(shè)計方案的缸內(nèi)壓力和溫度。

圖9 缸內(nèi)壓力變化曲線圖

圖10 缸內(nèi)溫度變化曲線圖

增大壓縮比后缸內(nèi)的壓力和溫度反而大幅降低，與作者的開發(fā)經(jīng)驗并不相符。為了深入研究大壓縮比時缸內(nèi)壓力降低的原因，額外考察了缸內(nèi)的氣體燃料質(zhì)量和循環(huán)放熱率，如圖11和12所示。

圖11 缸內(nèi)氣體燃料質(zhì)量變化曲線圖

圖12 燃燒放熱率變化曲線圖

針對圖12的積分顯示，4個設(shè)計方案燃燒過程的放熱量是基本相同的，這說明各方案消耗的燃料量相同。圖11中顯示壓縮沖程中方案三、四的氣體燃料質(zhì)量明顯大于方案一、二，考慮到氣缸內(nèi)的燃料總質(zhì)量相同，因此認為增大壓縮比時噴霧的碰壁現(xiàn)象得到了明顯改善。

氣體燃料較少使得空燃比降低，稀燃的現(xiàn)象更加明顯，稀燃是有利于燃燒火焰?zhèn)鞑サ?，因此方案一、二的燃燒放熱率更大，氣缸?nèi)的壓力和溫度明顯提高。

碰壁的液體燃料由于燃燒時的高溫得到了較快的蒸發(fā)，在做功行程的末期已經(jīng)完全燃燒，故4個方案的總放熱量相同。

3.6 NO排放預(yù)測

NO是汽油機的重要污染物，因此此次計算采用了Zeldovich模型對NO生成進行了評估，圖13顯示了各方案的NO質(zhì)量曲線。

圖13 NO質(zhì)量變化曲線圖

N2含量和溫度是NO生成時最重要的影響因素，由于各設(shè)計方案中空氣的含量基本相同，即N2質(zhì)量基本相同，因此此次計算中溫度成為最重要的影響因子。由于方案一、二的缸內(nèi)溫度明顯大于方案三、四，其生成的NO也明顯增多。

200 K的溫差帶來了近10倍的NO質(zhì)量差異，說明NO生成效率對2 000～3 000 K之間的溫度非常敏感，因此降低缸內(nèi)溫度將是減小NO裸機排放的重要手段。

4 總結(jié)

綜上所述：就此次分析的4個方案來看，增大壓縮比能顯著降低NO排放，而大滾流比的氣道則更有利于發(fā)動點火，因此設(shè)計工程師選定方案四作為最終方案，而方案二則被列為備選。

目前已經(jīng)根據(jù)此例中的設(shè)計方案四完成了樣機組建，初步的試驗表明它能夠滿足各項設(shè)計指標。

【1】周龍保.內(nèi)燃機學(xué)[M].北京：機械工業(yè)出版社，1999.

【2】解茂昭.內(nèi)燃機計算燃燒學(xué)[M].大連：大連理工大學(xué)出版社，1995.

【3】楊杰，蔡少娌，齊運亮，等.滾流強度對LPG發(fā)動機混合氣形成的影響[J].山東建筑大學(xué)學(xué)報,2009,24(5):396-401.

YANG J,CAI S L,QI Y L,et al.The Effects of the Tumble on the Formation of Stratified Mixture Structure in a Liquid Phase LPG Injection Engine[J].Journal of Shandong Jianzhu University,2009,24(5):396-401.

【4】藍志寶，秦際宏，楊曉，等.進氣道對增壓汽油機流動特性及性能的影響[J].車用發(fā)動機,2016(1):83-87.

LAN Z B,QIN J H,YANG X，et al.Influence of Intake Port on Air Flow Characteristics and Performance of Turbocharged Gasoline Engine[J].Vehicle Engine,2016(1):83-87.

【5】徐玉梁，傅光琦，孫吉龍，等.車用四氣門汽油機工作過程三維動態(tài)模擬[J].農(nóng)業(yè)機械學(xué)報,2007,38(7):53-57.

XU Y L,FU G Q,SUN J L,et al.Dynamic Simulation on the Working Process of Four-valve Gasoline Engines for Vehicles[J].Transactions of the Chinese Society for Agricultural Machinery,2007,38(7):53-57.

Optimization of Engine Design Based on In-cylinder Process Simulation

ZHANG Yingbing，YAO Wei

(Research & Development Center,Anhui Jianghuai Automobile Co.,Ltd.,Hefei Anhui 230022,China)

To choose the best one from several design schemes, CFD method was used to simulate the flow, injection, combustion process of a GDI engine that was in development. Evaluations about in-cylinder tumble and combustion were done, some advises were given for port design, injector chosen and timing optimization.

GDI engine; In-cylinder process simulation; Optimization design

2016-07-13

張應(yīng)兵(1979—)，男，工程師，研究方向為發(fā)動機產(chǎn)品的設(shè)計開發(fā)。E-mail:13349110312@qq.com。

10.19466/j.cnki.1674-1986.2016.10.009

U464

A

1674-1986(2016)10-041-04