PFI氫內(nèi)燃機(jī)的氫-空氣混合特性

段俊法， 劉福水， 孫柏剛

(1. 北京理工大學(xué)機(jī)械與車輛學(xué)院， 北京 100081； 2. 華北水利水電大學(xué)機(jī)械學(xué)院， 河南 鄭州 450045)

PFI氫內(nèi)燃機(jī)的氫-空氣混合特性

段俊法1,2， 劉福水1， 孫柏剛1

(1. 北京理工大學(xué)機(jī)械與車輛學(xué)院， 北京 100081； 2. 華北水利水電大學(xué)機(jī)械學(xué)院， 河南 鄭州 450045)

以1臺4缸進(jìn)氣道燃料噴射內(nèi)燃機(jī)為基礎(chǔ)，建立了包含進(jìn)氣管道和氣缸的氫內(nèi)燃機(jī)進(jìn)氣系統(tǒng)三維仿真模型。提出用混合氣不均勻性系數(shù)來評價(jià)氫-空氣混合狀況，仿真研究了氫內(nèi)燃機(jī)邊界條件下的進(jìn)氣氣流運(yùn)動規(guī)律和混合過程，發(fā)現(xiàn)氫-空氣混合狀況主要取決于氣流的搬運(yùn)作用和擴(kuò)散面積。探討了氫氣噴射相位和噴射壓力對氫-空氣混合速度和混合均勻性的影響，并以不均勻性系數(shù)為指標(biāo)優(yōu)化了PFI氫內(nèi)燃機(jī)不同轉(zhuǎn)速和當(dāng)量比下的噴射相位與噴射壓力。

氫內(nèi)燃機(jī)； 進(jìn)氣道燃料噴射； 混合氣形成； 均勻性； 噴射相位

隨著石油資源日益緊缺和環(huán)境污染日益加劇，研究者更加關(guān)注內(nèi)燃機(jī)替代燃料。氫氣可以從水分解獲得,燃燒的產(chǎn)物也只有水,因而可以作為循環(huán)利用的清潔能源。氫氣易于點(diǎn)燃，燃燒速度快，可以直接應(yīng)用于傳統(tǒng)內(nèi)燃機(jī)而不需要對內(nèi)燃機(jī)作較大的改動，因而是較為理想的內(nèi)燃機(jī)代用燃料[1-2]。氫氣作為燃料應(yīng)用于內(nèi)燃機(jī)可以采用進(jìn)氣管單點(diǎn)噴射、進(jìn)氣道燃料噴射及缸內(nèi)直噴等幾種方式。進(jìn)氣道燃料噴射(PFI)相對于缸內(nèi)直噴,技術(shù)難度較低，可精確控制氫氣噴射相位，獲得良好的燃燒和排放特性，是目前最常用的氫內(nèi)燃機(jī)燃料供給方式[3-5]。

燃料的混合狀況影響氫內(nèi)燃機(jī)的燃燒特性。火花塞處的氫-空氣混合氣濃度越高，越容易被點(diǎn)燃，燃燒速度越快[6-7]；氫氣和空氣的混合越充分，燃燒室內(nèi)的混合氣濃度越均勻，混合氣局部過濃的可能性越小，氫內(nèi)燃機(jī)的NOx排放越低。

本研究基于一款PFI氫內(nèi)燃機(jī)，建立了包含進(jìn)氣管道和氣缸的三維CFD仿真模型，仿真分析了氫氣射流運(yùn)動規(guī)律，研究了氫-空氣混合特性和氫氣噴射相位對火花塞局部濃度和燃燒室混合氣均勻性的影響，并以混合氣均勻性為目標(biāo)優(yōu)化了PFI氫內(nèi)燃機(jī)的氫氣噴射相位。

1 方法和模型

氫氣和空氣都是無色的，很難采用試驗(yàn)手段觀測到氫氣和空氣的混合過程。因而采用仿真方法分析氫氣噴射、擴(kuò)散、流動和混合的全過程，得到不同氫-空氣混合氣的時(shí)空分布特性。

1.1 仿真模型的建立

以1臺4缸PFI氫內(nèi)燃機(jī)為研究對象建立三維CFD仿真模型。內(nèi)燃機(jī)主要參數(shù)見表1。

表1 內(nèi)燃機(jī)參數(shù)

PFI氫內(nèi)燃機(jī)的氫氣以一定壓力噴射進(jìn)入進(jìn)氣道，被進(jìn)氣氣流運(yùn)送到氣缸內(nèi)，為了分析氫-空氣混合的全過程，需要建立整個(gè)氣流運(yùn)動路徑，含進(jìn)氣管、進(jìn)氣歧管、氣門和氣缸的仿真模型，為了獲得滿意的仿真效果，網(wǎng)格的尺寸很小，因而網(wǎng)格數(shù)很大。

為了縮短計(jì)算時(shí)間，建立仿真模型時(shí)忽略氫內(nèi)燃機(jī)各缸進(jìn)氣不一致，將進(jìn)氣管簡化成壓力穩(wěn)定的集氣腔，建立從集氣腔到氣缸的單缸幾何模型?？紤]氣缸和進(jìn)氣歧管的對稱性，取實(shí)體幾何模型的一半建立了仿真模型(見圖1)。經(jīng)簡化后仿真模型的網(wǎng)格數(shù)仍達(dá)17萬個(gè)。

仿真模型的邊界條件和初始條件主要通過氫內(nèi)燃機(jī)試驗(yàn)獲得，包括進(jìn)氣溫度、壓力、流量等。氫氣噴射壓力設(shè)定為0.3～0.6 MPa，傳熱系數(shù)、摩擦因數(shù)等參數(shù)和邊界條件根據(jù)經(jīng)驗(yàn)公式得到。

1.2 模型校核與驗(yàn)證

建立三維仿真模型后，校核了最大和最小網(wǎng)格尺寸及計(jì)算步長。隨著計(jì)算步長的減小，計(jì)算結(jié)果波動逐步減小，當(dāng)計(jì)算結(jié)果的波動小于2%時(shí)，認(rèn)為計(jì)算結(jié)果可以接受。

為了驗(yàn)證模型，采用仿真模型計(jì)算了不同工況的內(nèi)燃機(jī)進(jìn)氣量，進(jìn)行了試驗(yàn)驗(yàn)證。由表2可見，仿真結(jié)果和試驗(yàn)數(shù)據(jù)[8]吻合，最大誤差不超過5%。

表2 仿真和試驗(yàn)的內(nèi)燃機(jī)進(jìn)氣量

2 結(jié)果分析

2.1 氫氣在進(jìn)氣道內(nèi)的噴射和流動

圖2示出了氫氣在進(jìn)氣道內(nèi)噴射和流動的過程。在噴射初期，氫氣在一定壓力下噴入氣道時(shí)，形成具有一定貫穿距離和錐角的射流，射流碰到進(jìn)氣管壁后沿著管壁散開，在流束的中部，當(dāng)量燃空比很高，幾乎沒有空氣存在，而在流束的邊緣，氫氣向空氣中擴(kuò)散，混合氣濃度逐漸降低。

隨著噴射的持續(xù)，氫氣沿著流束表面逐步向空氣中擴(kuò)散，燃空比沿著射流表面的法向逐步降低，擴(kuò)散的范圍逐步增大。在噴射后期進(jìn)氣門打開，流動的空氣將較高濃度的混合氣向氣缸內(nèi)輸送，高濃度混合氣和更多的空氣接觸并向其中擴(kuò)散，形成濃度較低的混合氣。

氫氣噴流和汽油射流不同，汽油射流在流束的外側(cè)形成很多微小的液滴，這些液滴和流束通過表面蒸發(fā)向空氣中擴(kuò)散。而氫氣射流沒有出現(xiàn)打散的顆粒，因而射流的貫穿和空氣的輸送作用是增大擴(kuò)散面積、實(shí)現(xiàn)良好混合的基礎(chǔ)。

2.2 氫-空氣的缸內(nèi)流動和混合

在氣門開啟后，氫氣和空氣的混合氣進(jìn)入氣缸(見圖3)。在進(jìn)氣初期，氣流受到氣門和燃燒室壁面的導(dǎo)向作用而形成了兩個(gè)旋轉(zhuǎn)方向相反的渦團(tuán)，氣門下方旋向?yàn)檎臏u團(tuán)在前期占據(jù)主導(dǎo)，滾流比為正。隨著氣門升程的增長，由氣門上方進(jìn)入的氣流形成的渦團(tuán)能量越來越大，滾流比逐漸變?yōu)樨?fù)值，并在進(jìn)氣下止點(diǎn)形成了穩(wěn)定的滾流。

由圖4可見，進(jìn)氣門關(guān)閉后，由于滾流的影響，原進(jìn)氣門下方較濃的混合氣沿著滾流的旋轉(zhuǎn)方向向排氣門側(cè)運(yùn)動，而排氣門側(cè)缸壁附近較稀薄的混合氣則運(yùn)動到了進(jìn)氣門側(cè)。滾流越強(qiáng)，氫氣的流動范圍越大，越容易擴(kuò)散與空氣混合。

在壓縮沖程的末期，滾流逐步耗散，形成眾多微小的湍流，此時(shí)氫氣和空氣的氣團(tuán)較小，接觸面積進(jìn)一步增大，從而促進(jìn)兩者的混合。

2.3 噴射時(shí)刻對混合的影響

2.3.1 混合狀況的評價(jià)方法

氫氣與空氣的混合越充分，分子的分布越均勻，在一定溫度下的反應(yīng)時(shí)間越短，越有利于燃燒的組織。因此可以用混合氣的均勻性評價(jià)混合狀況。

將氣缸分成眾多的有限體積網(wǎng)格，認(rèn)為各網(wǎng)格內(nèi)氫氣與空氣的濃度相同，而網(wǎng)格之間的濃度不同，定義缸內(nèi)混合氣不均勻性系數(shù)nindex為

(1)

(2)

式中：φi為單元i的當(dāng)量燃空比；Vi為單元i的體積；V表示整個(gè)氣缸的體積；φavg表示缸內(nèi)平均當(dāng)量燃空比。

由式(1)可知，不均勻性系數(shù)nindex越小，各網(wǎng)格的燃空比越接近，氫氣與空氣的混合越充分。

2.3.2 噴氫相位對混合狀況的影響

噴氫相位影響混合氣當(dāng)量燃空比的空間分布。圖5示出了轉(zhuǎn)速3 000 r/min，當(dāng)量燃空比為1，噴射壓力為0.3 MPa時(shí)不同噴氫相位的缸內(nèi)混合氣分布云圖。

由圖5可知，噴氫相位不同，混合氣的混合狀況有較大差異。噴氫結(jié)束角為388°和611°時(shí)，所有體積微元的當(dāng)量燃空比在0.6～1.4之間；而558°時(shí)體積微元的當(dāng)量燃空比在0.76～1.24之間，且大部分體積微元的當(dāng)量燃空比在0.92～1.16之間；噴氫結(jié)束角為458°時(shí)混合氣濃度分布也比較均勻，混合較為充分。

以不均勻性系數(shù)表征各種噴氫相位時(shí)的混合狀況更為準(zhǔn)確。圖6示出了轉(zhuǎn)速 3 000 r/min，當(dāng)量燃空比為1，噴射壓力為0.3 MPa時(shí)不均勻性系數(shù)和噴氫相位的關(guān)系。由圖可知，噴氫結(jié)束角為558°時(shí)不均勻性系數(shù)為0.04，僅相當(dāng)于噴氫結(jié)束角388°時(shí)的36.2%,說明其混合更為充分。

混合氣均勻性隨噴氫結(jié)束角先減小后增大，原因是噴氫結(jié)束過早或過晚時(shí)的進(jìn)氣流速較低，氫氣團(tuán)被搬運(yùn)的能量小，氫氣與空氣的接觸面積小，擴(kuò)散速度降低，因而混合不夠充分，不均勻性系數(shù)增大。

通過仿真可得不同工況下混合最為充分的噴射相位(見圖7)。隨著轉(zhuǎn)速的提高和當(dāng)量燃空比的增大，使混合氣最均勻的噴射角增大，這主要是由氫氣量的增加和擴(kuò)散時(shí)間的縮短所造成的。

2.3.3 不同工況的最優(yōu)噴射壓力

隨著內(nèi)燃機(jī)轉(zhuǎn)速和負(fù)荷的提高，氫氣供給需求增加，需要增加噴氫持續(xù)期或提高噴射壓力。圖8示出了仿真得到的不同工況下使混合最為均勻的噴射壓力。

在轉(zhuǎn)速為1 000 r/min時(shí)，當(dāng)量燃空比增加最優(yōu)噴射壓力保持不變。此時(shí)保持最優(yōu)噴氫結(jié)束角增大，則噴氫開始角隨當(dāng)量燃空比增加而提前。這是由于在低轉(zhuǎn)速下，以曲軸轉(zhuǎn)角計(jì)量的噴射持續(xù)期較短，噴射壓力升高會使其進(jìn)一步縮短，造成短時(shí)間內(nèi)高濃度的混合氣，不利于形成均勻混合氣，故應(yīng)使用較低的噴射壓力。

由于噴射壓力不變，為獲得較大的當(dāng)量燃空比，噴氫持續(xù)期增大，因而最優(yōu)的噴射結(jié)束角相應(yīng)滯后以提高殘留在進(jìn)氣道內(nèi)混合氣的當(dāng)量燃空比，噴射開始角也由于噴射持續(xù)期的增加而提前。

在轉(zhuǎn)速升高為3 000 r/min 時(shí)，隨當(dāng)量燃空比的增加，最佳噴射壓力逐漸升高。這是由于噴射壓力低而轉(zhuǎn)速較高時(shí)，以曲軸轉(zhuǎn)角計(jì)的噴射持續(xù)期太長而無法避免在氣道內(nèi)形成高濃度的混合氣，因此必須提高噴射壓力。當(dāng)轉(zhuǎn)速提高后，進(jìn)氣回流角滯后，因此噴射結(jié)束角也隨轉(zhuǎn)速升高而滯后。

轉(zhuǎn)速進(jìn)一步升高到5 000 r/min，如果進(jìn)一步提前噴射，會造成進(jìn)氣道內(nèi)混合氣濃度過高而使得進(jìn)氣的當(dāng)量燃空比不均勻度增大。只能考慮進(jìn)一步增大噴氫壓力，此時(shí)因噴氫壓力的增大，噴氫閥前后的壓差超過臨界壓差，流動不穩(wěn)定，因而混合不夠充分。

3 結(jié)論

a) 氫氣與空氣的混合和汽油不同，不存在射流的破碎，主要通過空氣的流動搬運(yùn)氫氣，并通過湍流

實(shí)現(xiàn)更大的接觸面，增大氫氣的擴(kuò)散速度，獲得更好的混合性能；

b) 在各種工況下，都存在使混合氣不均勻性系數(shù)最小的噴氫相位和壓力，隨著轉(zhuǎn)速和燃空比的增大，使混合更為充分的噴氫開始角不斷提前，噴射壓力不斷增大。

[1] Delorme A,Rousseau A,Sharer P,et al.Evolution of hydrogen fueled vehicles compared to conventional vehicles from 2010 to 2045[C].SAE Paper 2009-01-1008.

[2] 楊振中.氫燃料內(nèi)燃機(jī)燃燒與優(yōu)化控制[D].杭州:浙江大學(xué),2001.

[3] Gerrit Kiesgen,Manfred Klüting,Christian Bock.The New 12-Cylinder Hydrogen Engine in the 7 Series:The H2 ICE Age Has Begun[C].SAE Paper 2006-01-0431.

[4] William F Stockhausen, Robert J Natkin,Daniel M Kabat,et al.Ford Hydrogen Engine Powered P2000 Vehicle[C] .SAE Paper 2002-01-0243.

[5] Swain M R,Schade G J,Swain M N.Design and Testing of a Dedicated Hydrogen Fueled Engine[C].SAE Paper 961077,1996.

[6] 李向榮，魏 嶸，孫柏剛，等.燃燒科學(xué)與技術(shù)[M].北京：北京理工大學(xué)出版社，2011.

[7] Nakamura Y,Yamamoto K,Nakajima N,et al.Effects of Hydrogen Enhancement on Efficiency and NOxEmissions of Lean and EGR-Diluted Mixtures in a SI Engine [C] .SAE Paper 2007-24-0122.

[8] 段俊法,劉福水,孫柏剛.PFI氫內(nèi)燃機(jī)回火機(jī)理和控制[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(3):1-5,37.

[編輯： 姜曉博]

Hydrogen-air Mixture Characteristics of PFI Hydrogen Internal Combustion Engine

DUAN Jun-fa1,2, LIU Fu-shui1, SUN Bai-gang1

(1. School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China; 2. School of Mechanical Engineering, North China University of Water Resources and Power, Zhengzhou 450045, China)

A 3D simulation model of air intake system including intake pipe and cylinder was built based on a 4-cylinder and intake port fuel injection hydrogen internal combustion engine. The non-uniformity coefficient was put forward to evaluate the mixing condition of hydrogen-air mixture. The intake air flow movement and mixing process of hydrogen internal combustion engine in the boundary conditions was simulated and researched with the model. It was found that the key factors for mixing condition were airflow carrying function and diffusion area. In addition, the influence of hydrogen injection pressure and timing on hydrogen-air mixing velocity and uniformity was investigated and the injection timing and pressure of different engine speed and equivalent ratio were optimized according to the non-uniformity coefficients.

hydrogen internal combustion engine； intake port fuel injection； mixture formation； uniformity； injection timing

2014-01-20；

2014-03-15

國家自然科學(xué)基金資助項(xiàng)目(51276019)；河南省教育廳科學(xué)技術(shù)研究重點(diǎn)支持項(xiàng)目(14B470016)

段俊法(1972—)，男，副教授，博士，主要研究方向?yàn)閮?nèi)燃機(jī)燃燒、排放和環(huán)境效應(yīng)；duanjunfa@126.com。

孫柏剛(1969—)，男，副教授，博士生導(dǎo)師，主要研究方向?yàn)閮?nèi)燃機(jī)燃燒、排放和環(huán)境效應(yīng)；sunbg@bit.edu.cn。

10.3969/j.issn.1001-2222.2014.03.007

TK431

B

1001-2222(2014)03-0029-04