李石彪，葛明生，沈黎明，羅福強(qiáng)

（1.江蘇大學(xué)汽車與交通工程學(xué)院，江蘇 鎮(zhèn)江 212013；2.江蘇四達(dá)動力集團(tuán)，江蘇 無錫 214000）

柴油機(jī)具有較好的可靠性、經(jīng)濟(jì)性和較高的熱效率，因而成為農(nóng)業(yè)機(jī)械和工程機(jī)械中較普遍的動力源之一［1］。日益嚴(yán)格的排放法規(guī)讓人們更加關(guān)注降低主要排放污染物的途徑［2－3］，排放物中最難平衡和處理的是NOx和炭煙，通過優(yōu)化燃燒室的結(jié)構(gòu)參數(shù)降低NOx和炭煙是機(jī)內(nèi)凈化的一種良好途徑。本研究設(shè)計(jì)3種不同結(jié)構(gòu)參數(shù)的燃燒室，設(shè)計(jì)時控制縮口率和徑深比不同，通過對3種燃燒室的數(shù)值模擬，了解缸內(nèi)油氣混合物的氣流運(yùn)動，分析排放物NOx和soot的生成歷程，為燃燒室結(jié)構(gòu)的進(jìn)一步設(shè)計(jì)和優(yōu)化提供了依據(jù)。

1 發(fā)動機(jī)參數(shù)和燃燒室結(jié)構(gòu)

本研究選用4缸增壓直噴柴油機(jī)，排量2.43L，壓縮比17.5∶1，標(biāo)定轉(zhuǎn)速3 000r/min，標(biāo)定功率55kW。

綜合考慮降低NOx和炭煙排放，設(shè)計(jì)了3種燃燒室（見圖1）。綜合壓縮比、縮口率（（d－D）/D）、口徑比（D/B）、徑深比（D/H）等結(jié)構(gòu)參數(shù)，設(shè)計(jì)時保持3種燃燒室有相同的壓縮比和燃燒室容積，控制3種燃燒室的口徑比在0.55±0.005以內(nèi)，3種燃燒室有不同縮口率（1號、2號、3號燃燒室分別為4.3，10.2，－11.8）和徑深比（1號，2號，3號燃燒室分別為3.0，2.6，3.6）。

縮口率會影響燃燒室內(nèi)的氣流活躍程度、渦流持續(xù)期、噴霧貫穿距離、噴霧錐角、破裂長度和霧化液滴的尺寸分布，進(jìn)而影響NOx和炭煙的生成。徑深比影響燃燒室中擠流和逆擠流的大小，設(shè)計(jì)的3種燃燒室縮口率和徑深比有較大差異，且有較高的凸臺，可以充分利用凸臺的導(dǎo)流作用形成較強(qiáng)的擠流和逆擠流，促進(jìn)油氣的混合［4］，促進(jìn)擴(kuò)散燃燒。

2 計(jì)算模型的選取及邊界條件的確定

2.1 模型的建立

本計(jì)算用Pro/e建立燃燒室模型，壓縮余隙為1mm，用Hypermesh對3個燃燒室進(jìn)行六面體網(wǎng)格劃分，1號燃燒室上止點(diǎn)網(wǎng)格數(shù)為111 120個，2號燃燒室上止點(diǎn)網(wǎng)格數(shù)為112 400個，3號燃燒室上止點(diǎn)網(wǎng)格數(shù)為114 720個，導(dǎo)入CFD軟件中生成180°～540°移動網(wǎng)格。

選取修正的κ－ε雙方程湍流模型［5］來模擬內(nèi)燃機(jī)缸內(nèi)流場，噴霧計(jì)算中用Walljet1模型模擬噴霧撞壁過程［6］，采用 Wave離散模型［7］，假設(shè)噴射的油滴與噴嘴出口直徑尺寸相同，其他分散的小油滴都是由于空氣流動促使液氣互相作用而形成的。液滴蒸發(fā)則采用 Dukowicz進(jìn)行模擬［8－9］，本計(jì)算選用Shell著火模型計(jì)算自燃時間［10］，選擇的燃燒模型為湍流渦破碎 Eddy－Break－Up模型［11］。排放模型分別選擇Extended－Zeldovich氮氧排放模型和Kennedy－Hiroyaso－Magnussen炭煙排放模型。

2.2 計(jì)算條件

計(jì)算工況為標(biāo)定轉(zhuǎn)速3 000r/min；計(jì)算范圍從218°（進(jìn)氣門關(guān)閉）到441°（排氣門打開），初始溫度383K，初始壓力185kPa；噴油器7孔，噴孔直徑0.145mm，安裝傾斜角為13°，位置相對于燃燒室中心，沿X軸正方向偏移2mm，沿Y軸正方向偏移1mm；油嘴伸出量為3.4mm，單缸循環(huán)噴油量34.5mg。噴油提前角7°，噴油持續(xù)期353°～369°。氣缸蓋底平面溫度551K，氣缸壁面溫度450K，活塞頂凹坑溫度593K。

3 計(jì)算結(jié)果及分析

3.1 模型的驗(yàn)證

2號燃燒室缸內(nèi)壓力實(shí)測值與計(jì)算值的對比見圖2。從圖2可以看出，試驗(yàn)值和計(jì)算值較吻合。

3.2 缸內(nèi)流場

從圖3可以看出，在2°BTDC時，2號燃燒室由于有較大的縮口率，且有較大的凹坑和較高的凸臺，油束的兩側(cè)已經(jīng)有一大一小旋向相反的擠壓渦流，大的渦流在燃燒室的凹坑壁面以及凸臺的導(dǎo)流作用下，沿著凸臺邊緣以較快的速度流動，能夠較好地促進(jìn)油氣的混合。1號燃燒室除了旋向相反的擠壓渦流之外，靠近燃燒室邊緣處仍有小的渦流，小渦流能促進(jìn)邊緣的油氣混合。與1號和2號燃燒室相比，3號燃燒室僅在油束的周圍形成兩個旋向相反的較大的擠壓渦流，在敞口拐角處有小渦流，因?yàn)?號燃燒室口徑較大，本身并不組織進(jìn)氣渦流，受進(jìn)氣等因素的影響及靠近上止點(diǎn)時受到擠壓作用產(chǎn)生少量渦流。在5°ATDC時，隨著活塞向下止點(diǎn)運(yùn)動，1號和2號燃燒室之前的大擠壓渦流開始變小、破碎，在油束附近形成多個小渦流，渦流中心慢慢向油束靠近，凸臺附近氣流運(yùn)動加強(qiáng)，可以有效地將附著在凸臺及凸臺附近的油霧帶到燃燒室頂隙；但2號燃燒室油束附近的渦流與燃燒室凹坑處的渦流相互干擾，從而不利于凹坑處空氣的有效利用和油氣的充分燃燒，主要原因是徑深比過小即凹坑較深，且在活塞向下運(yùn)動時，燃燒室內(nèi)的高壓氣體進(jìn)入擠流區(qū)導(dǎo)致油束相互干擾。3號燃燒室在膨脹過程中擠壓作用減弱，導(dǎo)致燃燒室內(nèi)的渦流進(jìn)一步減弱。

3.3 缸內(nèi)湍動能

湍動能是氣流活躍程度評價的重要指標(biāo)，湍動能的大小在很大程度上反映油氣混合的程度，從而影響著后期的擴(kuò)散燃燒。從圖4可以看出3種燃燒室的湍動能不盡相同，但皆在上止點(diǎn)后約10°達(dá)到峰值，且從上止點(diǎn)到上止點(diǎn)后10°保持著較高的湍動能，即在這段時間內(nèi)燃燒室保持較活躍的氣流運(yùn)動，能較好地促進(jìn)油氣混合，促進(jìn)擴(kuò)散燃燒的進(jìn)行。從圖中可以看出1號和2號燃燒室比3號燃燒室湍動能大，缸內(nèi)的氣流運(yùn)動更活躍，主要原因是1號和2號燃燒室有較大的縮口率，能很好地組織氣流運(yùn)動。

3.4 缸內(nèi)壓力、溫度和放熱規(guī)律比較

由圖5可見，2號燃燒室放熱最早最快，在366.5°時達(dá)到放熱峰值，1號燃燒室次之，3號燃燒室放熱峰值出現(xiàn)最晚，比2號燃燒室推遲3.5°。主要原因是2號燃燒室縮口率最大，進(jìn)氣渦流較強(qiáng)，油氣混合較快且充分，霧化最好，放熱最早。2號燃燒室缸內(nèi)最高壓力高于1號和3號燃燒室；最高缸內(nèi)溫度略低于1號燃燒室，但溫度升高早，當(dāng)其達(dá)到最高溫度時，1號燃燒室缸內(nèi)溫度仍快速升高；3號燃燒室缸內(nèi)溫度最高。原因主要是2號燃燒室內(nèi)初期燃燒的可燃混合氣數(shù)量相對較多，放熱更早且集中，壓力升高更早且缸內(nèi)壓力峰值最高，1號燃燒室次之。

3.5 NOx和soot的生成

1號燃燒室和2號燃燒室缸內(nèi)平均溫度相差不大，但2號燃燒室初期氣流速度最小，燃燒室中心和頂隙處的氣流運(yùn)動較弱且無較大渦流，油氣混合較差，燃燒不夠充分，生成的NOx最少（見圖6）。而3號燃燒室缸內(nèi)平均溫度要高于1號、2號燃燒室，且高溫持續(xù)期較長，其次，3號燃燒室燃燒初期氣流速度最快，油氣混合充分，促進(jìn)前期燃燒，因此，3號燃燒室的NOx生成量較多。2號燃燒室的soot生成量要多于1號燃燒室（見圖7），原因是2號燃燒室有較大的縮口率和較小徑深比，從而有較長的渦流持續(xù)期，因此，2號燃燒室后期氣流速度最大，氣流最活躍，但由于燃燒室內(nèi)油束相互干擾，影響油束的霧化質(zhì)量，影響擴(kuò)散燃燒，soot大量生成。3號燃燒室的soot生成量最少，原因是3號燃燒室有較大的徑深比，能充分地利用新鮮空氣加速soot的氧化。

3.6 不同進(jìn)氣渦流比下排放特性

從圖8至圖10可以看出，3種燃燒室在不同渦流比下，soot和NOx排放的變化趨勢不同。進(jìn)氣渦流比在2.0～2.6范圍內(nèi)，隨著進(jìn)氣渦流比的增大，1號和2號燃燒室NOx排放量減少，且縮口率越大減少的幅度越大。而3號燃燒室隨著渦流比的增大，NOx排放大幅增加，其中渦流比從2.0變化到2.3時，NOx排放增加約11%。主要原因是1號與2號燃燒室的縮口率比3號燃燒室大，相對于3號燃燒室增強(qiáng)了縱向的進(jìn)氣渦流和壓縮渦流，并不需要太大的進(jìn)氣渦流燃燒室的油氣已經(jīng)較好地混合，較大的進(jìn)氣渦流導(dǎo)致油束相互干擾，從而影響前期的預(yù)混合和燃燒，降低NOx排放。因此對于1號和2號燃燒室而言，較大的進(jìn)氣渦流比能降低NOx排放。3號燃燒室結(jié)構(gòu)本身并不組織進(jìn)氣渦流，由氣缸蓋和進(jìn)氣道形狀產(chǎn)生的進(jìn)氣渦流能很好地促進(jìn)油氣的混合，從而改善前期燃燒，更大的進(jìn)氣渦流能使油氣燃燒充分，缸內(nèi)壓力和溫度升高，NOx排放增加。從圖中可以看出，渦流比對soot前期生成有較大的影響，但是對后期soot的氧化過程影響小，可見進(jìn)氣渦流對后期的擴(kuò)散燃燒影響較小。

4 試驗(yàn)研究

4.1 試驗(yàn)條件和試驗(yàn)設(shè)備

對3種燃燒室進(jìn)行臺架試驗(yàn)，試驗(yàn)采用WE42N水力測功機(jī)、EIM0310D發(fā)動機(jī)測控系統(tǒng)及FCM－13瞬態(tài)油耗測量儀。進(jìn)水溫度為60～75℃，出水溫度為（88±5）℃，柴油溫度控制在（38±5）℃，機(jī)油溫度（95±5）℃，采用輔助冷卻手段。試驗(yàn)工況與計(jì)算工況相同（3 000r/min，75%負(fù)荷），測量3種燃燒室NOx和炭煙排放特性。

4.2 試驗(yàn)結(jié)果

從圖11可以看出，3號燃燒室的NOx排放最高，1號燃燒室次之，2號燃燒室最低，而3號燃燒室炭煙排放最低。試驗(yàn)結(jié)果與模擬結(jié)果變化趨勢相同，1號燃燒室有較為折中的排放特性。

5 結(jié)論

a）當(dāng)燃燒室采用較小徑深比和較大縮口率時，前期稍小的氣流速度和較弱的氣流運(yùn)動使缸內(nèi)溫度不至太高，氧氣濃度較低，有利于減少NOx的生成；當(dāng)采用較小的縮口率和較大的徑深比時，后期稍大的氣流速度有利于油氣混合，促進(jìn)后期擴(kuò)散燃燒，加速soot的氧化；

b）不同進(jìn)氣渦流比下，縮口率存在較大差異的3種燃燒室NOx排放呈現(xiàn)不同的變化趨勢，主要取決于3種燃燒室本身的結(jié)構(gòu)特征，負(fù)縮口率的3號燃燒室在較大的進(jìn)氣渦流比下油氣混合更充分，燃燒也更充分，NOx排放增加；進(jìn)氣渦流主要影響前期燃燒，對NOx的生成影響較大，對soot的氧化過程影響較小。

［1］ Jaichandar S，Annamalai K.Effects of open combustion chamber geometries on the performance of pongamia biodiesel in a DI diesel engine［J］.Fuel，2012，98：272－279.

［2］ 陳 利.燃燒室形狀對170柴油機(jī)性能影響的分析［D］.天津：天津大學(xué)，2008.

［3］ 許正偉.直噴柴油發(fā)動機(jī)燃燒室形狀和噴射規(guī)律研究［D］.成都：西華大學(xué)，2009.

［4］ 宮長明，劉巽俊，譚滿志，等.低排放的直噴式柴油機(jī)燃燒室研究［J］.農(nóng)業(yè)機(jī)械學(xué)報(bào)，1999，30（1）：18－22.

［5］ 周 磊，趙長祿，張付軍，等.利用燃燒模擬對柴機(jī)燃燒室的優(yōu)化設(shè)計(jì)［J］.燃燒科學(xué)與技術(shù)，2004，10（5）：466－470.

［6］ 陶文銓.數(shù)值傳熱學(xué)［M］.2版.西安：西安交通大學(xué)出版社，2001.

［7］ Sachina E M，Sazhin S S，Heikal M R，et al.The shell autoignition model：applications to gasoline and diesel fuels［J］.Fuel，1999，78：389－401.

［8］ Kaminage Takash，Kusaka Jin.Improvement of combustion and exhaust gas emission in a passenger car diesel engine by modification of combustion chamber design［C］.SAE Paper 2006－01－3436.

［9］ 焦運(yùn)景，張惠明，田 遠(yuǎn)，等.直噴式柴油機(jī)燃燒室?guī)缀涡螤顚ε欧庞绊懙亩嗑S數(shù)值模擬研究［J］.內(nèi)燃機(jī)工程，2007，28（4）：11－15.

［10］ Halstead M P，Kirsch L J，Prothero A，et al.A mathematical Model for Hydrocarbon Autoignition at High Pressures［C］//Proc R Soc.London：［s.n.］，1975（A346）：515－538.

［11］ Magnussen B F，Hjertager B H.On mathematical modeling of turbulent combustion with special emphasis on soot formation and combustion［C］//Proceedings of the Sixteenth International Symposium on Combustion.Pittsburgh：The Combustion Institute，1977.