沈現(xiàn)青，黃昭明

(1.合肥職業(yè)技術(shù)學(xué)院，安徽 合肥238000;2.河海大學(xué)文天學(xué)院，安徽 馬鞍山243031)

0 引 言

近年來，隨著計算機輔助工程技術(shù)(Computer Aided Engineering，CAE)在內(nèi)燃機領(lǐng)域內(nèi)應(yīng)用的日益深入，計算流體力學(xué)(Compute Fluid Dynamics，CFD)的發(fā)展越來越成熟，對內(nèi)燃機缸內(nèi)三維燃燒過程的仿真模擬計算已成為現(xiàn)實.數(shù)值模擬不僅可以節(jié)約大量的發(fā)動機臺架試驗時間和成本，也能獲取試驗中無法測試的發(fā)動機工作過程信息［1，2］.天然氣是一種清潔、高效的新型發(fā)動機能源，在發(fā)動機上燃用能夠取得良好的經(jīng)濟和社會效益，因此對天然氣發(fā)動機的燃燒過程進行仿真計算，獲取燃燒的相關(guān)信息將具有實際意義［3，4］.

基于上述背景，本文采用三維CFD 軟件CONVERGE 建立了天然氣發(fā)動機工作過程的三維計算模型，對燃燒過程進行了數(shù)值模擬，并進行了NOx排放物的計算研究，以為天然氣發(fā)動機的仿真模擬計算和實際應(yīng)用提供技術(shù)參考.

1 天然氣發(fā)動機概述

天然氣是一種無色、無味的氣體，主要成分是甲烷，密度比空氣小，燃燒時火焰呈淡藍色，火焰溫度較高，天然氣的理化特性如表1 所示.

由于天然氣的能量密度較低，目前應(yīng)用于汽車上時采用壓縮鋼瓶的存儲方式，壓縮天然氣(Compressed Natural Gas，CNG)通過一套氣體燃料供給系統(tǒng)，進入發(fā)動機進氣歧管與空氣混合，以形成均質(zhì)混合氣，進入缸內(nèi)燃用.典型的天然氣發(fā)動機氣體供給系統(tǒng)如圖1 所示.

表1 天然氣與汽油的理化特性對比

圖1 天然氣發(fā)動機的供給系統(tǒng)

2 燃燒模型的建立及試驗

目前已有多款CFD 軟件能用來計算內(nèi)燃機的缸內(nèi)過程和進行排放預(yù)測.本文選用了CONVERGE 軟件來進行仿真模擬計算.該款軟件是專門針對發(fā)動機工作過程的計算模擬來進行開發(fā)的，先進的網(wǎng)格劃分理念和網(wǎng)格自適應(yīng)技術(shù)能夠保證發(fā)動機動網(wǎng)格的精確生成.

2.1 試驗用發(fā)動機基本參數(shù)

試驗用發(fā)動機的基本參數(shù)如表2 所示.主要進行了發(fā)動機功率點的缸壓曲線測量和NOx排放量測量的試驗.

表2 試驗用汽油機主要參數(shù)

2.2 幾何模型和網(wǎng)格劃分

從實體模型中提取汽油機建模所需要的幾何元素，包括進氣管、氣道、氣門、缸蓋、燃燒室、活塞、火花塞型面和噴射器型面，對幾何體進行修正處理后如圖2 所示.

圖2 燃燒建模的幾何模型

發(fā)動機工作過程的三維數(shù)值模擬需要合適的網(wǎng)格數(shù)量和形狀.網(wǎng)格數(shù)量過少，雖然計算時間縮短，但不能滿足湍流模型精度要求;網(wǎng)格數(shù)量過多，則計算時間太長，也不利于模擬計算工作的進行，因此需要選擇合適的網(wǎng)格總數(shù)和網(wǎng)格尺寸.同時，對于邊界層和其它需要加密網(wǎng)格的部位(如噴射器部位和點火部位)，應(yīng)確保網(wǎng)格的密度.

采用CONVERGE 的自適應(yīng)技術(shù)劃分汽油機的動網(wǎng)格，基本網(wǎng)格尺寸為0.004m，網(wǎng)格總數(shù)不超過2.5×106個，能夠保證網(wǎng)格良好的收斂性和計算性，網(wǎng)格截面如圖3 所示.

圖3 汽油機動網(wǎng)格的網(wǎng)格截面

2.3 工作過程模型建立

在所完成幾何模型修正和動網(wǎng)格的劃分后，進行邊界條件的設(shè)置，以完成汽油機燃燒過程計算模型的建立.邊界條件包括缸內(nèi)計算邊界條件和氣體燃料噴射的邊界條件等，合理設(shè)置邊界條件以確保燃燒過程仿真計算的準確進行.

發(fā)動機工作過程的仿真模擬計算中，燃燒模型的選定至關(guān)重要.采用CHEMKIN 化學(xué)反應(yīng)動力學(xué)軟件模擬天然氣燃燒的化學(xué)反應(yīng)過程和反應(yīng)機理，生成天然氣反應(yīng)機理的數(shù)據(jù)文件供CONVERGE軟件使用.天然氣的著火方式與汽油一樣，采用點燃方式，需要設(shè)置點火源項參數(shù)，燃燒模型選擇SAGE 模型，NOx排放物的生成模型選用擴展的Zeldovich 氮氧化物模型.所建立的汽油機燃燒過程計算模型如圖4 所示.

表3 燃燒模型和點火源項的參數(shù)輸入

圖4 汽油機燃燒過程的計算模型

圖5 不同曲軸轉(zhuǎn)角下的缸內(nèi)溫度場分布

3 仿真計算結(jié)果分析

3.1 燃燒仿真計算結(jié)果

應(yīng)用所建立的發(fā)動機燃燒過程計算模型對汽油機運行于功率點(85kW/5200r/min)的燃燒過程進行了計算.

圖5 給出了上止點(After Top Dead Center，ATDC)后從15°～95°CA 缸內(nèi)溫度場的分布.圖6為缸壓曲線計算結(jié)果與試驗數(shù)據(jù)的對比圖.

從圖5 和圖6 可以知道，應(yīng)用CONVERGE 軟件建立的計算模型能夠進行燃燒過程的仿真計算，天然氣最高燃燒溫度超過2800K，有利于擴大溫度循環(huán)梯度，提高熱效率;且火焰在傳播過程中呈球形，上止點后95°CA 的缸內(nèi)溫度場分布(圖5(d))顯示出進氣側(cè)部位有終燃混合氣自燃情況發(fā)生，但屬于輕微爆震，有利于燃燒熱效率的提高;缸內(nèi)爆發(fā)壓力曲線的仿真計算數(shù)據(jù)和試驗數(shù)據(jù)吻合較好，說明計算過程準確、有效.

圖6 缸內(nèi)壓力計算數(shù)據(jù)和試驗數(shù)據(jù)對比

圖7 缸內(nèi)NOx 排放物濃度場分布

3.2 NOx 排放量的計算結(jié)果

采用所建立的發(fā)動機燃燒過程計算模型對汽油機運行于功率點(85kW/5200r/min)的NOx排放量進行了計算.排氣門開啟時刻的缸內(nèi)NOx濃度分布如圖7 所示.

汽油機功率點的NOx濃度仿真模擬計算結(jié)果為2.153×10－3kg/kg.

汽油機臺架試驗的NOx排放量單位為ppm(10－6)，需要換算成比排放量，NOx體積濃度排放量與質(zhì)量排放量的關(guān)系［5］，可表示為公式(1):

式中，MNOx為NOx質(zhì)量流量，單位為kg/kg;φNOx為NOx濃度，單位為ppm;kHD為大氣溫度和濕度校正系數(shù)，一般取0.95 ～1.05.

通過公示(1)換算后，汽油機功率點NOx試驗測量值為1.938×10－3kg/kg，與計算結(jié)果差距不大，說明模型能用來預(yù)測NOx排放量.

4 結(jié) 論

(1)CONVERGER 軟件建立的天然氣發(fā)動機燃燒模型的計算結(jié)果表明，天然氣在發(fā)動機上的最高燃燒溫度超過2800K，有利于熱效率的提高.

(2)天然氣發(fā)動機燃燒過程仿真計算出的缸內(nèi)爆發(fā)壓力數(shù)據(jù)與臺架試驗數(shù)據(jù)較為吻合，說明模型的建立及計算過程準確.

(3)天然氣發(fā)動機氮氧化物質(zhì)量分數(shù)計算結(jié)果為2.153×10－3kg/kg，而臺架試驗結(jié)果為1.938×10－3kg/kg，計算精度較高，說明模型能夠預(yù)測天然氣發(fā)動機NOx排放.

［1］ 陳泓，倪計民，葉年業(yè)，等.計算機輔助工程技術(shù)在廢氣再循環(huán)系統(tǒng)開發(fā)中的應(yīng)用研究［J］.機械工程學(xué)報，2012，48(22):155－161.

［2］ 陸小明，張世鷹，葛蘊珊，等.生物柴油發(fā)動機燃燒與排放數(shù)值模擬及試驗研究［J］.內(nèi)燃機，2008，6:30－33.

［3］ 阿里·古巴提.天然氣汽車應(yīng)用技術(shù)研究［D］.西安:長安大學(xué)，2014，5.

［4］ 歐翔飛，羅東曉.國內(nèi)壓縮天然氣汽車產(chǎn)業(yè)發(fā)展分析［J］.天然氣工業(yè)，2007，27(4):129－131.

［5］ 殷勇，項旭昇，陽松林.重型柴油機達國ⅣSCR 試驗研究［J］.汽車科技，2008，2:61－65.

［6］ JIAO，Yun－jing，DONG Hong，ZHANG Hui－ming，et al.Effect of Combustion Chamber Geometry on Combustion Process in Spark Ignition Natural Gas Engine［J］.Chinese Internal Combustion Engine Engineering，2009，30(4):28－33.

［7］ ZHENG，Zun－qing，YAO Ming－fa，LIN Zhi－qiang，et al.Development of Combustion System of an Electric Controlled Lean Burn CNG Engine for Meeting China Stage ⅣEmission Regulations［J］.Chinese Internal Combustion Engine Engineering，2009，30(5):7－11.

［8］ Tomita E，Nobuyuki K，Piao Z Y.et al.Effects of EGR and Early Injection of Diesel Fuel on Combustion Characteristics and Exhaust Emissions in a Methane Dual－Fuel Engine［C］.SAE Paper 2002012723，2002.

［9］ Simon K，Chen.Gas engine combustion principle and applications［C］.SAE Paper 2001012489，2001.

［10］ 尤秋菊，張來斌，王朝暉，等.天然氣發(fā)動機供氣系統(tǒng)的研究現(xiàn)狀［J］.車用發(fā)動機，2006，5(10):1－3.

［11］ 吳聿東.CNG 發(fā)動機燃氣供給系統(tǒng)仿真及噴射控制研究［D］.成都:西華大學(xué)，2013，6.