潘鐘鍵，周凱，程立志

(1.長沙學院 機電工程學院，長沙 410003)(2.山河智能裝備股份有限公司 航空研究院，長沙 410100)(3.中南大學 機電工程學院，長沙 410083)

0 引 言

活塞發(fā)動機是目前地面交通工具、通用航空裝備的主要動力[1]，特別是燃油經(jīng)濟性好的重油(航空煤油或者柴油)發(fā)動機，其排放問題越來越受到人們的關(guān)注[2]，尤其是針對即將開放的通用航空市場。通用航空近幾年增長迅猛，截止到2017年，美國的通用航空飛機數(shù)量超過了300 000架，一些輕型飛機從航空汽油發(fā)動機更換成航空重油發(fā)動機，但排放污染問題仍一直存在。國內(nèi)外在發(fā)動機的排放上進行了大量研究，提出了多種控制排放的方法，但絕大部分都是從發(fā)動機后處理、噴油控制策略等方面進行研究。國外，M.Capobianco等[3]對帶有旁通閥的渦輪增壓器進行了非定常流動性能試驗，研究排放物的生成規(guī)律；T.Tang等[4]研究了柴油氧化催化劑和顆粒物過濾之間的關(guān)系以改善排放；N.Stamoudis等[5]對重油蒸發(fā)燃料液滴模型的對流傳熱進行了適當處理，精確模擬了燃燒排放物的生成。國內(nèi)，楊帥等[6]設(shè)計了一套電子EGR系統(tǒng)，確保不同工況下的EGR率，來降低NOx的排放；樓狄明等[7]、孫萬臣等[8]從氧化傳感器和微粒捕集器方面對排放進行了大量研究，以改善柴油排放顆粒物的生成；倪計民等[9]采用可變渦輪增壓和廢氣循環(huán)的方法對排放進行改善；鄒憲等[10]對不同噴油策略下的排放和油耗進行了研究。

本文從燃燒室的基本參數(shù)出發(fā)，通過改變結(jié)構(gòu)參數(shù)，對發(fā)動機某一瞬態(tài)下的排放進行研究，總結(jié)不同燃燒室結(jié)構(gòu)參數(shù)對排放的影響，以期為發(fā)動機的排放研究提供參考。

1 噴霧數(shù)學模型建立

在直噴發(fā)動機缸內(nèi)，高壓噴入的重油會發(fā)生一系列物理變化，例如重油液滴的蒸發(fā)、破碎、碰撞、耗散等。為了更好地研究燃燒過程中的重油特性，針對上述過程對其建立數(shù)學模型。

1.1 重油液滴蒸發(fā)模型

采用Dukowicz模型建模，并做如下假設(shè)：①整個重油液滴內(nèi)部的溫度相等；②球形液滴；③液滴周圍具有穩(wěn)態(tài)氣相流場；④液滴表面準靜態(tài)；⑤氣液相交處熱力學平衡。

根據(jù)以上假設(shè)，重油液滴被噴入氣缸內(nèi)，重油液滴吸收的能量主要用于自身加熱和蒸發(fā)，方程為

(1)

式中：md為液滴的質(zhì)量；cpd為重油定壓比熱容；L為液滴的蒸發(fā)潛熱；α為對流傳熱系數(shù)；As為重油液滴表面積。

根據(jù)熱力學平衡條件，質(zhì)量流的表達方程為

(2)

式中：fvs為蒸發(fā)質(zhì)量流量；qs為重油液滴表面熱流密度；Q為氣流傳給液滴的熱量。

(3)

設(shè)Le=1且質(zhì)量傳遞系數(shù)：

(4)

由單滴關(guān)聯(lián)性理論得到，Nusselt數(shù)值由下式所得

Nμ=2+0.6Re1/2Pr1/3

(5)

式中：Re和Pr可通過實驗獲得。

綜合式(1)～式(5)，用Nusselt數(shù)替代對流傳熱系數(shù)α，重油液滴從周圍吸收的熱量為

Q=DdπλNμ(T∞-Ts)

(6)

式中：Dd為液滴的直徑；λ為波長。

1.2 重油破碎模型

重油噴入氣缸后發(fā)生破碎，初級破碎發(fā)生在距高壓噴嘴較近的區(qū)域，次級破碎發(fā)生區(qū)域較遠。選用Wave模型估算液滴半徑。

r=2πσB/ρ2U2

(7)

式中：B為量級常數(shù)1；U為氣液兩相運動速差；σ為液滴表面張力。

實際計算公式為

(8)

式中：Λ為對應波長；a為振幅臨界值；Ω為最大擾動波增長速率。

當不穩(wěn)定波振幅大于臨界值，重油液滴發(fā)生破碎分裂。

1.3 碰壁模型

由于活塞式發(fā)動機體積較小，其工作的氣缸直徑亦較小，重油噴射過程中，不可避免地與氣缸壁面發(fā)生碰撞，碰撞后產(chǎn)生多種運動軌跡，例如粘附、反彈、貼壁運動等。

采用walljet液滴碰壁模型，并假設(shè)液滴碰壁后，壁面液膜對液滴運動不產(chǎn)生影響，碰撞后的運動形式根據(jù)參數(shù)We來判斷，液滴撞擊壁面后粒徑大小與韋伯數(shù)相關(guān)，其方程為

(9)

式中：d1為撞壁前液滴直徑；d0為撞壁后液滴直徑。

當We值低于50時，用反射碰壁模型；當We值大于50時，采用射流碰壁模型。

1.4 湍流擴散模型

(10)

式中：Rni為速度分量隨機數(shù)(0～1)；k為液滴的湍動能；erf-1高斯函數(shù)反函數(shù)。

重油液滴與氣體湍流作用的時間tturb為

(11)

重油液滴在湍流團中的作用時間t>tturb時，液滴進入下一湍流渦團[11]。

2 燃燒室模型的建立與分析

將重油燃燒的噴霧燃燒模型應用到燃燒分析軟件中，F(xiàn)ire軟件是發(fā)動機燃燒模擬方面的專業(yè)軟件，由美國AVL公司推出，目前大部分發(fā)動機研發(fā)公司都在使用該軟件對發(fā)動機燃燒進行模擬[12]。根據(jù)發(fā)動機公司提供的發(fā)動機參數(shù)，輸入到Fire相應的模塊中，例如噴油孔直徑、上止點間隙等，但燃燒室的形狀需要發(fā)動機公司提供具體圖紙獲得。在AVL Fire軟件中，模擬燃燒過程是從進氣門關(guān)閉到排氣門開啟為止，該發(fā)動機提供的原始數(shù)據(jù)是進氣門關(guān)閉角度為下止點后36°，排氣門開啟角度為下止點前48°，按通用設(shè)計原則，定義上止點為720 °CA，下止點為540 °CA，燃燒模擬過程是從576 °CA到848 °CA，設(shè)定發(fā)動機轉(zhuǎn)速為2 300 r/min。

在進行精確的燃燒模擬之前，需要對燃燒室進行網(wǎng)格劃分，燃燒過程中，某曲軸轉(zhuǎn)角下的動態(tài)網(wǎng)格，由Fire軟件對其進行網(wǎng)格劃分，如圖1所示。

圖1 740 °CA時的燃燒室模型網(wǎng)格

為了對燃燒室進行優(yōu)化設(shè)計，在原有方案的基礎(chǔ)上，提供新的優(yōu)化方案，具體數(shù)據(jù)如表1所示。

表1 燃燒室基本尺寸

燃燒室基本尺寸示意圖如圖2所示。為了研究啞鈴型燃燒室基本尺寸對排放的影響，在優(yōu)化時，選定兩組方案，一組改變?nèi)紵业淖畲笾睆胶涂s口直徑，另一組在前一組的基礎(chǔ)上繼續(xù)改變中間深度和最大深度。對兩組設(shè)計方案進行燃燒數(shù)值模擬分析。

圖2 燃燒室基本尺寸

對三種方案進行燃燒分析，并進行對比。由于燃燒主要在活塞達到上止點附近進行，在活塞下行10 °CA時出現(xiàn)最大壓力和最大溫度，由于NOx和Soot的生成與溫度有直接關(guān)系，選定730 °CA作為分析的瞬時狀態(tài)，該狀態(tài)下NOx和Soot排放量的多少直接反映最終排放物的多少，因此對該角度進行分析，缸內(nèi)溫度場如圖3所示。

(a) 初始方案

(b) 方案一

(c) 方案二

從圖3可以看出：燃燒室的結(jié)構(gòu)尺寸發(fā)生變化后，缸內(nèi)溫度值也發(fā)生了變化，火焰?zhèn)鞑ヒ?guī)律基本沒有發(fā)生較大的變化，最大溫度出現(xiàn)位置基本和火焰?zhèn)鞑ヂ窂较嗤环桨敢坏臏囟容^原方案降低了71.4 K，方案二較原方案溫度降低了101.2 K，溫度的降低對NOx和Soot的生成起到了一定的抑制效果。

730 °CA下的NOx濃度分布情況如圖4所示，用質(zhì)量分數(shù)來表示其生成速率和濃度分布，分布位置與缸內(nèi)溫度場基本相似，溫度是NOx生成的主要原因。

(a) 初始方案

(b) 方案一

(c) 方案二

從圖4可以看出：NOx化合物最大生成量有所下降，方案一在原方案基礎(chǔ)上下降了10.7%，方案二在原方案基礎(chǔ)上下降23%，總體下降趨勢明顯。

對730 °CA下修改前后的方案進行Soot生成進行分析，如圖5所示，用質(zhì)量分數(shù)表示Soot的生成速率和濃度分布情況。Soot生成主要發(fā)生在擴散燃燒階段，因此其分布與溫度場分布略有不同，從圖5可以看出：方案一在原方案基礎(chǔ)上Soot的最大濃度上升了5.6%，方案二在原方案基礎(chǔ)上Soot的最大濃度下降了2.8%。

(a) 初始方案

(c) 方案二

綜上所述，方案二帶來的NOx排放下降明顯，Soot生成也略有降低，從總的結(jié)果來看，排放趨于降低趨勢。

3 實 驗

按方案二的燃燒室尺寸進行制造，完成裝配后進行臺架測試，接入AVL分析儀，實驗如圖6所示。

對發(fā)動機進行測試，記錄發(fā)動機基本參數(shù)和排放數(shù)據(jù)，記錄數(shù)據(jù)如表2所示。

圖6 發(fā)動機排放測試

轉(zhuǎn)速/rpm功率/kW扭矩/(N·m)不透光煙度/m-1NOx/ppm原機方案二原機方案二2 700712500.250.248636852 500622350.170.166715282 300431800.070.074783752 100301360.020.02298232

從表2可以看出：發(fā)動機功率扭矩基本沒有發(fā)生變化，方案二所對應的排放與原機測得數(shù)據(jù)確有下降趨勢，下降數(shù)據(jù)與仿真計算的數(shù)據(jù)基本相符，表明仿真計算方法真實可信。

4 結(jié) 論

(1) 縮小發(fā)動機燃燒室的最大直徑尺寸和縮口尺寸能降低發(fā)動機缸內(nèi)燃燒溫度，能降低NOx的生成，但Soot生成量有所增加。最大直徑尺寸降低4.4%，溫度降低71.4 K，NOx生成降低10.7%，Soot的最大生成量上升了5.6%。

(2) 燃燒室最大深度和中間深度也能對缸內(nèi)溫度以及排放產(chǎn)生影響，中間深度增加17.2%，最大深度增加14.7%，NOx生成降低23%，Soot最大生成降低2.8%。

(3) 對修改后的燃燒室進行臺架實驗，發(fā)現(xiàn)發(fā)動機功率和扭矩基本無變化，但排放物生成確有減少，證明燃燒室的優(yōu)化是可行。