包凌志，孫柏剛，汪 熙

（北京理工大學(xué) 機(jī)械與車輛學(xué)院，北京 100081，中國）

當(dāng)前，化石燃料的能源危機(jī)和排放法規(guī)的日益嚴(yán)苛是內(nèi)燃機(jī)行業(yè)所面臨的巨大挑戰(zhàn)，因此內(nèi)燃機(jī)必須向高效率、低排放，甚至零排放的方向發(fā)展[1]。H2燃燒的產(chǎn)物為水，沒有碳排放，所以也被認(rèn)為是一種未來的、清潔的能源載體[2]。

當(dāng)前氫能的應(yīng)用主要有2種方式：氫燃料電池和氫燃料內(nèi)燃機(jī)。氫燃料電池具有高效、無噪音的特點(diǎn)，且生成物只為水。對(duì)于另1種方式，H2具有寬廣的可燃范圍、極快的燃燒速度和極低的點(diǎn)火能量，非常適合作為傳統(tǒng)內(nèi)燃機(jī)的替代燃料[3]。相比于燃料電池，氫內(nèi)燃機(jī)成本低、工作范圍廣、可以利用工業(yè)廢氫。寶馬、豐田等開發(fā)的氫內(nèi)燃機(jī)汽車進(jìn)一步證明了氫內(nèi)燃機(jī)具有良好的動(dòng)力性、經(jīng)濟(jì)性和實(shí)用性[4]。

氫內(nèi)燃機(jī)的進(jìn)氣組織方式主要分為缸內(nèi)直噴 （direct injection， DI）和 缸 外 進(jìn) 氣 道 噴 射（port fuel injection，PFI）[3]2種。PFI 氫內(nèi)燃機(jī)是實(shí)現(xiàn)氫能源應(yīng)用最簡潔、經(jīng)濟(jì)和可靠的途徑[5]。然而采用進(jìn)氣道噴射時(shí)，H2與空氣在進(jìn)氣道混合時(shí)會(huì)占用氣缸容積，在化學(xué)當(dāng)量比下可以達(dá)到30%，導(dǎo)致功率下降，同時(shí)易發(fā)生回火和早燃的問題[6]。DI氫內(nèi)燃機(jī)則能解決以上問題，有效避免異常燃燒，大幅提升動(dòng)力性（比汽油機(jī)高17%），是當(dāng)前研究的熱點(diǎn)和未來技術(shù)發(fā)展方向[7]。

由于H2絕熱火焰溫度高，依據(jù)澤爾多維奇機(jī)理在燃燒的過程中，氫內(nèi)燃機(jī)會(huì)產(chǎn)生NOx排放，隨著溫度的升高和氧濃度的增加，NOx生成量將呈指數(shù)型增長。在直噴氫內(nèi)燃機(jī)中，由于噴射時(shí)刻晚，H2與空氣的混合時(shí)間短，混合氣濃度分布不均勻， NOx最高可達(dá)8 000 ppm[8]。因此如何控制氮氧排放，使NOx近乎降低為0，是直噴氫內(nèi)燃機(jī)發(fā)展亟需解決的問題。

控制NOx排放的方法主要分為缸內(nèi)凈化和缸外處理2種方法。對(duì)直噴氫內(nèi)燃機(jī)，通過調(diào)整噴射相位、優(yōu)化噴霧結(jié)構(gòu)、采用多次噴射[9]的方法，可以有效組織缸內(nèi)混合氣分布，避免局部過濃區(qū)域的出現(xiàn)。同時(shí)，采用噴水[10]或者廢氣再循環(huán)[11]（exhaust gas recirculation，EGR），可以降低缸內(nèi)燃燒溫度，從而降低排放。缸外處理方面，在大負(fù)荷工況下，利用傳統(tǒng)的三元催化器[12]，或者采用氮氧化物儲(chǔ)存還原系統(tǒng)（NOxstorage-reduction，NSR），可以達(dá)到98%的NOx轉(zhuǎn)化率。

本文基于一臺(tái)自然吸氣的2.0 L直噴氫內(nèi)燃機(jī)，以探索其原排（無后處理）近零排放（NOx≤10×10-6，即20 ppm）特性為目標(biāo)，重點(diǎn)研究了過量空氣系數(shù)、轉(zhuǎn)速和點(diǎn)火角對(duì)排放的影響，并通過多參數(shù)的耦合調(diào)節(jié)，得到了樣機(jī)的最大動(dòng)力性下的零排放工作邊界。

1 試驗(yàn)方法

試驗(yàn)用直噴氫內(nèi)燃機(jī)是基于一臺(tái)2.0 L自然吸氣直列四缸直噴汽油機(jī)改造而來，具體參數(shù)如表 1所示。改造工作主要包括：對(duì)原機(jī)缸蓋擴(kuò)孔、安裝高壓H2噴嘴、搭建高壓供氫系統(tǒng)、更換冷型火花塞、設(shè)計(jì)加工金屬進(jìn)氣歧管等。

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

直噴氫內(nèi)燃機(jī)通常需要在壓縮沖程中將H2噴入缸內(nèi)，因此噴射壓力選擇為6 ～ 10 MPa，由于H2體驅(qū)動(dòng)方式，文獻(xiàn)[13]對(duì)流量和噴霧特性開展了定容試驗(yàn)分析。

試驗(yàn)采用的噴嘴是外開軸針式噴嘴，具有90°的軸針錐角和外開式環(huán)形噴孔，噴嘴結(jié)構(gòu)示意圖如圖 1所示。

圖1 噴嘴內(nèi)部結(jié)構(gòu)示意圖

H2從壓力為35 MPa高壓儲(chǔ)氫系統(tǒng)出來后，經(jīng)過兩級(jí)減壓閥，減壓到噴射所需壓力。為解決大流量噴嘴工作過程中，供氫壓力波動(dòng)的問題，試驗(yàn)設(shè)計(jì)了一個(gè)容量為250 ml的氫軌，能使壓力波動(dòng)率小于0.3%，整體試驗(yàn)臺(tái)架布置如圖 2所示。

圖 2 直噴氫內(nèi)燃機(jī)臺(tái)架布置

試驗(yàn)臺(tái)架上安裝有CW250電渦流測(cè)功機(jī)，可測(cè)得輸出扭矩和功率。氫內(nèi)燃機(jī)上安裝了曲軸轉(zhuǎn)角傳感器和Kistler 6118火花塞式缸內(nèi)壓力傳感器，并利用Kibox燃燒分析儀實(shí)時(shí)輸出缸壓、放熱率等燃燒參數(shù)。

試驗(yàn)采用CMF010、CMF025科氏流量計(jì)（精度 ± 0.1% FSa）和熱線式空氣流量計(jì)（精度± 0.1% FSa）測(cè)量試驗(yàn)H2和空氣流量。采用高精度的排氣分析儀Horiba 7110和AVL DiGas 4000分別測(cè)量三元催化器前后的NOx排放。試驗(yàn)采用Motohawk快速原型控制系統(tǒng)，對(duì)發(fā)動(dòng)機(jī)的節(jié)氣門、噴射、點(diǎn)火等進(jìn)行電子控制。

試驗(yàn)過程中，保證冷卻水溫度85℃，機(jī)油溫度95℃。在不同噴射壓力下，通過調(diào)整噴氫脈寬，調(diào)節(jié)過量空氣系數(shù)λ。待發(fā)動(dòng)機(jī)在當(dāng)前工況穩(wěn)定運(yùn)轉(zhuǎn)2 min后，測(cè)量并計(jì)算200個(gè)循環(huán)的試驗(yàn)數(shù)據(jù)平均值。

2 試驗(yàn)結(jié)果及分析

2.1 過量空氣系數(shù)對(duì)NO x排放的影響

試驗(yàn)條件為發(fā)動(dòng)機(jī)轉(zhuǎn)速n= 2 000 r/min，平均有效壓力 ( brake mean effective pressure, BMEP） = 0.3 MPa，此時(shí)發(fā)動(dòng)機(jī)的扭矩為48.2 Nm，功率穩(wěn)定在10 kW。噴氫相位固定在壓縮上止點(diǎn)前曲柄轉(zhuǎn)角（CA）220°，保證H2-空氣的充分預(yù)混，噴射壓力設(shè)定為10 MPa，過量空氣系數(shù)λ從1掃略到3.1，催化器前NOx隨濃度的變化如圖 3所示。氮氧排放隨著過量空氣系數(shù)增加，先增大到3 800×10-6（即3 800 ppm）后減少，在λ＞2.5后逐漸趨于零排放。直噴氫內(nèi)燃機(jī)的排放對(duì)過量空氣系數(shù)十分敏感，在λ為1.0 ～ 2.2的區(qū)間內(nèi)，NOx排放都很高，因此近零排放的控制策略應(yīng)避開這一區(qū)間，選擇λ≥ 2.5的稀燃區(qū)域。

圖3 NO x排放隨過量空氣系數(shù)變化

指示熱效率變化如圖4所示，隨著過量空氣系數(shù)的增加，均先增加后趨于穩(wěn)定。這主要因?yàn)橄∪己蠊?jié)氣門開度增加，泵氣損失減少。另一方面，排氣溫度隨濃度降低均勻減少，從當(dāng)量比時(shí)的462 ℃減少至λ= 3.1時(shí)的296 ℃。這時(shí)傳熱損失和排氣能量損失也相對(duì)降低，熱效率有所提升，最高指示熱效率可以達(dá)到40%。

圖5顯示了50% 燃燒點(diǎn)AI50和燃燒持續(xù)期的變化趨勢(shì)，可以看出AI50在所有工況中均保持在CA = 8°～11°之間，證明此時(shí)點(diǎn)火提前角為最佳扭矩點(diǎn)火角（maximum brake torque timing，MBT）。燃燒持續(xù)期在λ＜ 1.6時(shí)，基本穩(wěn)定在曲柄轉(zhuǎn)角為 20°；當(dāng)λ= 1.82時(shí)，增加至曲柄轉(zhuǎn)角為25°，最后變?yōu)榍D(zhuǎn)角為 32°。 這個(gè)趨勢(shì)主要與H2火焰?zhèn)鞑ニ俣扔嘘P(guān)，稀燃后火焰?zhèn)鞑ニ俣葴p慢，燃燒持續(xù)期增加。

缸壓曲線的變化如圖 6所示，隨著過量空氣系數(shù)的減少，節(jié)氣門開度變小，進(jìn)氣量減少，爆發(fā)壓力逐漸降低，但最大壓力升高率卻相對(duì)增加。λ= 1時(shí)，缸壓曲線呈現(xiàn)火花點(diǎn)火發(fā)動(dòng)機(jī)明顯的定容燃燒的趨勢(shì)，此時(shí)火焰?zhèn)鞑ニ俣瓤?，燃燒放熱率高。如圖 7所示，在化學(xué)當(dāng)量比處的瞬時(shí)放熱率（ηex）峰值可以達(dá)到60 J每度曲柄轉(zhuǎn)角（CA），遠(yuǎn)高于其他混合氣濃度下的放熱率，由于燃燒持續(xù)期短，化學(xué)當(dāng)量比工況下的點(diǎn)火提前角應(yīng)接近上止點(diǎn)點(diǎn)火。

2.2 λ = 2.5時(shí)的稀燃特性

由前面的結(jié)論，在進(jìn)氣量最大的條件下，混合氣濃度選定為λ= 2.5，就能在保證發(fā)動(dòng)機(jī)動(dòng)力性的同時(shí)實(shí)現(xiàn)近零排放。為驗(yàn)證這一結(jié)果，稀燃特性試驗(yàn)的噴氫相位為上止點(diǎn)前曲柄轉(zhuǎn)角160°，即進(jìn)氣門關(guān)閉時(shí)刻，防止進(jìn)氣回流; 此時(shí)H2與空氣混合時(shí)間短，在高轉(zhuǎn)速工況下，混合氣為分層燃燒。保持節(jié)氣門開度100%，噴射壓力為6 MPa，進(jìn)氣可變氣門正時(shí)（variable valve timing，VVT）為初始位置，排氣VVT提前曲柄轉(zhuǎn)角60°，調(diào)節(jié)噴氫脈寬使得過量空氣系數(shù)穩(wěn)定在2.5，點(diǎn)火提前角取最佳扭矩點(diǎn)火角 （MBT），試驗(yàn)結(jié)果如圖8所示。隨著轉(zhuǎn)速增加，功率先增加，在4 500 r/min時(shí)達(dá)到28.2 kW后減少，而扭矩則逐漸下降，扭矩最高為n= 2 000 r/min時(shí)的74.3 Nm。

圖4 指示熱效率（ηshow）、排氣溫度（θex） 隨過量空氣系數(shù)（λ）變化

圖 6 不同過量空氣系數(shù)缸壓隨曲軸轉(zhuǎn)角變化

圖 7 不同過量空氣系數(shù)瞬時(shí)放熱率隨曲軸轉(zhuǎn)角變化

平均有效壓力（BMEP）和NOx排放的變化如圖 9所示，隨著轉(zhuǎn)速增加，BMEP從0.46 MPa下降到0.18 MPa，而NOx排放在1 500 ～ 4 000 r/min都穩(wěn)定在100×10-6（即100 ppm）左右，在4 500 r/min急劇增長到793×10-6（即793 ppm），并在之后一直保持較高水平。這里主要是因?yàn)椴缓线m的進(jìn)排氣VVT角度造成的，在高轉(zhuǎn)速下進(jìn)氣門關(guān)閉時(shí)刻過晚，排氣開啟時(shí)刻過早，導(dǎo)致進(jìn)氣量不足，燃燒惡化，排放增加，排氣能量增加，熱效率下降。

平均指示壓力變動(dòng)系數(shù)COVIMEP是氫發(fā)動(dòng)機(jī)用于表示燃燒循環(huán)變動(dòng)最常用的參數(shù)，計(jì)算公式為：

圖8 稀燃工況下功率扭矩隨轉(zhuǎn)速變化曲線

圖9 稀燃工況下平均有效壓力（BEMP）和NO x隨轉(zhuǎn)速變化

其中：σIMEP是平均指示壓力的標(biāo)準(zhǔn)偏差，PIMEP是平均指示壓力的平均值， COVIMEP一般以3%為界限，低于3% 時(shí)發(fā)動(dòng)機(jī)平穩(wěn)工作。

圖 10顯示了循環(huán)變動(dòng)系數(shù)在稀燃工況下的變化曲線。在4 000 r/min 之前，循環(huán)變動(dòng)系數(shù)小于3%，發(fā)動(dòng)機(jī)工作在穩(wěn)定狀態(tài)；隨著轉(zhuǎn)速提高，燃燒情況逐漸惡化，導(dǎo)致循環(huán)變動(dòng)系數(shù)很高。因此近零排放在保證發(fā)動(dòng)機(jī)平穩(wěn)運(yùn)轉(zhuǎn)的前提下，主要探索1 500 ～ 4 000 r/min范圍內(nèi)的工況。

采用稀燃可以大幅降低NOx排放，當(dāng)?shù)捎趪姎湎辔坏耐七t，混合時(shí)間短，在缸內(nèi)不能形成均質(zhì)混合氣，并會(huì)出現(xiàn)局部濃區(qū)，這點(diǎn)在高轉(zhuǎn)速工況更加明顯，導(dǎo)致λ= 2.5時(shí)， NOx大部分在100×10-6（ 即100 ppm）左右，仍然不能達(dá)到近零排放。因此需要對(duì)其他參數(shù)進(jìn)行調(diào)整，進(jìn)一步控制排放。

2.3 點(diǎn)火提前角對(duì)NO x排放的影響

控制點(diǎn)火提前角是調(diào)整燃燒相位、控制排放的有效手段。試驗(yàn)保持λ= 2.5，節(jié)氣門全開，噴射相位為上止點(diǎn)前曲柄轉(zhuǎn)角160°，從圖11中可以看出，在各個(gè)轉(zhuǎn)速工作下，點(diǎn)火角曲柄轉(zhuǎn)角往后推遲約5°后，NOx排放從原先的90×10-6（即90 ppm）左右，迅速下降至12×10-6（即12 ppm），達(dá)到近零排放的目標(biāo)。指示熱效率的變化如圖12所示，推遲點(diǎn)火角后，雖然排放降低，但犧牲了約1%的熱效率，這主要是因?yàn)橥七t點(diǎn)火角后，放熱率曲線整體后移，爆發(fā)壓力降低，缸內(nèi)溫度降低。

圖 10 稀燃工況下平均有效壓力變動(dòng)系數(shù)（COVBMEP） 隨轉(zhuǎn)速變化

圖11 NO x排放推遲點(diǎn)火角前后隨轉(zhuǎn)速變化

如圖13所示，推遲點(diǎn)火提前角還會(huì)帶來循環(huán)變動(dòng)系數(shù)的增加，燃燒不穩(wěn)定性提高，但平均指示壓力變動(dòng)系數(shù) （COVIMEP）都滿足小于3%穩(wěn)定運(yùn)轉(zhuǎn)的要求。因此，對(duì)于直噴氫內(nèi)燃機(jī)，適當(dāng)調(diào)整點(diǎn)火角，可以在滿足工況要求的條件下，進(jìn)一步降低排放。

2.4 直噴氫內(nèi)燃機(jī)近零排放邊界探索

結(jié)合之前的結(jié)論，近零排放試驗(yàn)固定噴射相位為上止點(diǎn)前曲柄轉(zhuǎn)角160°，從而給予H2-空氣充足的混合時(shí)間，保證發(fā)動(dòng)機(jī)工作穩(wěn)定性。節(jié)氣門開度為100%，轉(zhuǎn)速范圍為1 500 ～ 4 000 r/min，以NOx近零排放（NOx≤20×10-6）為目標(biāo)條件，噴射壓力為6 MPa和8 MPa，設(shè)置初始λ= 2.5，并通過不斷調(diào)整噴氫脈寬和點(diǎn)火提起角，得到最大動(dòng)力性的工作邊界，試驗(yàn)結(jié)果如圖14所示。在起始段，平均有效壓力穩(wěn)定在0.47 MPa左右，但隨著轉(zhuǎn)速增加，BMEP逐漸降低，在n= 4 000 r/min時(shí)達(dá)到0.32 MPa。增大噴射壓力，噴霧流速和貫穿距增加，缸內(nèi)氣體流速和湍流強(qiáng)度增強(qiáng)，更容易形成均質(zhì)稀燃混合氣，從而有效降低排放，提升動(dòng)力。從圖14中可以看出，提高噴射壓力，BMEP在各個(gè)工況點(diǎn)下，大約都能提升20 kPa。因此, 不同噴射壓力下，近零排放的控制策略都不相同。后續(xù)的試驗(yàn)還將對(duì)更高噴射壓力進(jìn)行測(cè)試，從而拓展零排放的工作范圍。

功率變化如圖15所示，隨轉(zhuǎn)速逐漸增加，在轉(zhuǎn)速為2 500 r/min后趨于穩(wěn)定，噴射壓力為6 MPa時(shí)，達(dá)到20 kW左右，而噴射壓力為8 MPa時(shí)，最大功率為21.5 kW。

圖12 指示熱效率推遲點(diǎn)火角前后隨轉(zhuǎn)速變化

圖13 推遲點(diǎn)火角前后平均指示壓力變動(dòng)系數(shù)（COVIMEP）隨轉(zhuǎn)速變化

圖14 不同噴射壓力下近零排放工作區(qū)域

圖15 不同噴射壓力下功率隨轉(zhuǎn)速變化曲線

圖14、圖15曲線下方近零排放區(qū)域可以通過控制節(jié)氣門開度、控制濃度、調(diào)整負(fù)荷的方式達(dá)到相應(yīng)的工況點(diǎn)，在此區(qū)域內(nèi)都為近零排放。

過量空氣系數(shù)的變化如圖16所示，在6 MPa的噴射壓力下，3 000 r/min之前，λ均控制在2.5左右，而在3 500 r/min時(shí)，調(diào)整點(diǎn)火角已經(jīng)無法控制NOx，為保證排放低于20×10-6（即20 ppm），將λ調(diào)整至2.84。而噴射壓力為8 MPa時(shí)，λ= 2.73就可以滿足排放需求。

如圖17所示，指示熱效率ηshow先增加后減少，在2 000 ～3 000 r/min工作區(qū)域，熱效率都在39%左右。但是轉(zhuǎn)速增加后，在高轉(zhuǎn)速工況時(shí)熱效率急劇降級(jí)，這主要是因?yàn)檫@時(shí)的混合氣過稀，燃燒不穩(wěn)定，此時(shí)為了實(shí)現(xiàn)低排放的目標(biāo)而損失了大量的燃料經(jīng)濟(jì)性。

圖16 不同噴射壓力下過量空氣系數(shù)隨轉(zhuǎn)速變化

圖17 不同噴射壓力下指示熱效率隨轉(zhuǎn)速變化曲線

近零排放的平均有效壓力的循環(huán)變動(dòng)系數(shù)變化如圖18所示，隨著轉(zhuǎn)速增加，高轉(zhuǎn)速下進(jìn)氣流速加快，同時(shí)過量空氣系數(shù)也逐漸增大，在稀燃工況下發(fā)動(dòng)機(jī)工作不穩(wěn)定性增加，高噴射壓力有利于混合氣的形成，使用8 MPa噴射壓力后，循環(huán)變動(dòng)系數(shù)大約能降低0.5%，在近零排放區(qū)域內(nèi)總體工作穩(wěn)定。

圖18 不同噴射壓力下平均有效壓力（COVIMEP）隨轉(zhuǎn)速變化

3 結(jié) 論

本文基于一臺(tái)自然吸氣的2.0 L直噴氫內(nèi)燃機(jī)，探索了其近零排放的工作邊界，得到如下結(jié)論：

1） 噴氫內(nèi)燃機(jī)的NOx排放對(duì)過量空氣系數(shù)十分敏感，在均質(zhì)混合氣的前提下，當(dāng)λ≥2.5時(shí)， NOx排放為0。

2） 噴氫時(shí)刻設(shè)置在進(jìn)氣門關(guān)閉前，會(huì)造成進(jìn)氣回流，內(nèi)燃機(jī)動(dòng)力性下降。推遲噴H2后，由于混合不均勻短，當(dāng)λ= 2.5時(shí)，NOx排放仍為90×10-6，且轉(zhuǎn)速增加后排放劇烈惡化。

3） 適當(dāng)推遲點(diǎn)火角，可以在犧牲少量熱效率和發(fā)動(dòng)機(jī)平穩(wěn)性的同時(shí)，實(shí)現(xiàn)近零排放（NOx≤20×10-6）。

4） 通過稀燃和推遲點(diǎn)火角，得到了直噴氫內(nèi)燃機(jī)近零排放邊界，最大功率可達(dá)21.5 kW，最高指示熱效率為39%，后續(xù)采用增壓系統(tǒng)、增大噴射壓力、優(yōu)化稀燃燃燒參數(shù)還可以繼續(xù)擴(kuò)展此區(qū)域的工作范圍。

參考文獻(xiàn)(References)

[1] Kalghatgi G. Is it really the end of internal combustion engines and petroleum in transport [J].Appl Energ, 2018, 225: 965-974.

[2] Abe J O, Popoola A P I, Ajenifuja E, et al. Hydrogen energy, economy and storage: Review and recommendation [J].Int'l J Hydro Energ, 2019, 44(29): 15072-15086.

[3] Verhelst S, Wallner T. Hydrogen-fueled internal combustion engines [J].Prog Energ Combust, 2009, 35(6): 490-527.

[4] Verhelst S. Recent progress in the use of hydrogen as a fuel for internal combustion engines [J].Int'l J Hydro Energ, 2014, 39(2): 1071-1085.

[5] LUO Qinghe, HU Jibin, SUN Baigang, et al. Experimental investigation of combustion characteristics and NOxemission of a turbocharged hydrogen internal combustion engine [J].Int'l J Hydro Energ, 2019, 44(11): 5573-5584.

[6] Luo Qinghe, Hu Jibin, Sun Baigang, et al. Effect of equivalence ratios on the power, combustion stability and NOx controlling strategy for the turbocharged hydrogen engine at low engine speeds [J].Int'l J Hydro Energ, 2019. 44(31): 17095-17102.

[7] Wimmer A, Wallner T, Ringler J, et al. H2-direct injection – A highly promising combustion concept [C] //SAE 2005 World Congress & Exhibition 2005, Michigan USA, 2005: 2001-2005.

[8] Welch A, Mumford D, Munshi S, et al. Challenges in developing hydrogen direct injection technology for internal combustion engines [C] //Powertrains, Fuels and Lubricants Meeting2008, Illinois USA, 2008: 2001-2008.

[9] Wallner T, Scarcelli R, Nande A M, et al. Assessment of multiple injection strategies in a direct-injection hydrogen research engine [J].SAE Int'l J Engines, 2009, 2(1): 1701-1709.

[10] Younkins M, Wooldridge M S, Boyer B A. Port injection of water into a di hydrogen engine [C]// SAE 2015 World Congress & Exhibition 2015, Michigan USA, 2015: 2011-2015.

[11] Tsujimura T, Suzuki Y. Development of a large-sized direct injection hydrogen engine for a stationary power generator [J].Int'l J Hydro Energ, 2019, 44(22): 11355-11369.

[12] BAO Lingzhi, SUN Baigang, LUO Qinghe, et al. Simulation and experimental study of the NOxreduction by unburned H2in TWC for a hydrogen engine [J].Int'l J Hydro Energ, 2020, 45(39): 20491-20500.

[13] WANG Xi, SUN Baigang, LUO Qinghe, et al. Visualization research on hydrogen jet characteristics of an outward-opening injector for direct injection hydrogen engines [J].Fuel, 2020, 280: 118710.