潘建欣， 劉康儒， 寧 楠， 楊 軍， 吳 桐， 謝曉峰

(1. 武漢船用電力推進(jìn)裝置研究所， 湖北 武漢 430064;2. 中國船舶集團(tuán)有限公司第七一八研究所， 河北 邯鄲 056027;3. 清華大學(xué) 核能與新能源技術(shù)研究院， 北京 100084)

1 前 言

在未來社會對高效、清潔、安全的能源體系的追求下，質(zhì)子交換膜燃料電池(PEMFC)成為當(dāng)今備受關(guān)注的綠色電源系統(tǒng)。在PEMFC 系統(tǒng)中，常利用引射器或氫氣循環(huán)泵對電堆出口的氫氣引射回流，在提高氫氣利用率的同時保持電堆性能。其中，氫氣引射器具有體積小、無寄生功率損耗、無噪聲等特點，是實現(xiàn)大功率PEMFC 氫氣循環(huán)的主流方向。然而，對于傳統(tǒng)的單噴嘴引射器而言，在引射器結(jié)構(gòu)固定后，難以滿足燃料電池系統(tǒng)的全功率強(qiáng)變載的需求。

近年來，國內(nèi)外學(xué)者開展了旨在提升引射器性能的結(jié)構(gòu)參數(shù)與仿真優(yōu)化研究。Zhu 等[1]提出了一種用于分析引射器流動特性的速度函數(shù)，利用計算流體力學(xué)(CFD)結(jié)合數(shù)據(jù)回歸和參數(shù)辨識方法，確定速度函數(shù)的指數(shù)，并研究了不同電流下氫氣子系統(tǒng)的再循環(huán)性能。Maghsoodi 等[2]討論了引射器結(jié)構(gòu)對燃料電池系統(tǒng)性能的影響，并利用神經(jīng)網(wǎng)絡(luò)和遺傳算法尋優(yōu)。Yin 等[3]建立引射器三維數(shù)值模型，研究二次流操作條件對引射器性能的影響，結(jié)果表明，二次流的濕度大、溫度高會降低二次流中的氫氣含量。Kou 等[4]建立了80 kW PEMFC 系統(tǒng)的Matlab/Simulink 模型，研究操作條件和幾何結(jié)構(gòu)對引射器性能的影響，但文中所建立的燃料電池模型未考慮水汽的影響。為擴(kuò)大引射器在燃料電池系統(tǒng)中的適用功率范圍，一些學(xué)者提出不同的引射器結(jié)構(gòu)設(shè)計與氫氣子系統(tǒng)設(shè)計優(yōu)化。Le Tri 等[5]將引射器的一維模型與車用PEMFC 集成，研究發(fā)現(xiàn)相較于單引射器狹窄的工作范圍，雙引射器可以實現(xiàn)燃料電池工作范圍的全覆蓋。Brunner等[6]提出一種可變流量引射器，在一次流噴嘴處設(shè)置一針閥，改變噴嘴的開口大小，該引射器在Ballard Mark 9 SSL 電堆中性能表現(xiàn)優(yōu)異。Xue 等[7]提出一種多噴嘴噴射器設(shè)計思路，通過實驗和CFD 仿真，認(rèn)為該設(shè)計可以滿足燃料電池寬范圍的功率輸出，同時實現(xiàn)相對穩(wěn)定的陽極入口壓力。Jenssen 等[8]在Brunner 可變流量引射器的基礎(chǔ)上，采用多級串聯(lián)燃料電池堆設(shè)計，提升引射器低功率下的性能。Chen等[9]提出了一種嵌套式引射器，并設(shè)置一條與引射器平行的旁路以進(jìn)一步擴(kuò)大其工作范圍，該引射器的引射性能可以滿足150 kW 的燃料電池堆在9%～100% 輸出功率范圍內(nèi)的供氫需求。目前對于引射器陽極氫氣回收性能的研究主要集中在利用計算流體力學(xué)軟件對其結(jié)構(gòu)參數(shù)進(jìn)行仿真優(yōu)化，或?qū)嵨镌谒罱ǖ囊淦鳒y試臺架中進(jìn)行氣體流動實驗。雖然CFD 方法有助于深入了解引射器內(nèi)部的流動行為，但其模型僅簡單地模擬了陽極的特性，沒有考慮動態(tài)響應(yīng)過程中電堆膜間發(fā)生的氫氮滲透與水生成等因素對氫氣循環(huán)過程及引射特性的影響。

文中利用gFUELCELL 仿真建模平臺，研究了不同噴嘴開關(guān)策略下四噴嘴引射器的引射能力，得到最優(yōu)噴嘴控制策略，分析該策略下電堆內(nèi)部物理場變化特性。研究得到燃料電池系統(tǒng)全運行功率下，多噴嘴引射器的最優(yōu)控制策略，實現(xiàn)了搭載四噴嘴引射器的150 kW 流量從5%～100% 的近似線性調(diào)控，解決了燃料電池系統(tǒng)低功率區(qū)間引射器失效的問題。研究結(jié)果為多噴嘴引射器在燃料電池系統(tǒng)中的應(yīng)用與選型提供了有價值的參考。

2 四噴嘴引射器氫氣子系統(tǒng)

2.1 引射器工作原理

引射器是一種進(jìn)行質(zhì)量、能量交換的機(jī)械裝置，其結(jié)構(gòu)由一次進(jìn)流噴嘴、二次回流管、吸入腔、混合管和擴(kuò)散管組成。傳統(tǒng)單噴嘴氫氣引射器結(jié)構(gòu)如圖1 所示。在燃料電池系統(tǒng)工作過程中，高壓氫氣經(jīng)減壓閥調(diào)節(jié)后，進(jìn)入一次進(jìn)流噴嘴。在通過一次進(jìn)流噴嘴收斂處，氣體流速增大，壓力減小，在吸入腔內(nèi)形成一個低壓區(qū)域。當(dāng)二次流體回流進(jìn)口壓力高于吸入腔壓力時，二次流體會在壓差作用下被吸入引射器中，與一次進(jìn)流一并流入混合管內(nèi)混合?；旌虾蟮牧黧w在擴(kuò)散管內(nèi)速度逐漸降低，得到一股滿足燃料電池工作壓力要求的引射流體。

圖1 氫氣引射器結(jié)構(gòu)圖Fig.1 Schematic diagram of hydrogen ejector structure

2.2 氫氣子系統(tǒng)

燃料電池氫氣循環(huán)子系統(tǒng)如圖2 所示，四噴嘴氫氣引射器的控制技術(shù)路線與單噴嘴引射器基本相同，主要由壓力調(diào)節(jié)閥、引射器、燃料電池堆、氣水分離器和排氫閥組成。從氫罐供應(yīng)的氫氣經(jīng)過減壓閥和壓力調(diào)節(jié)閥，隨后進(jìn)入氫氣引射器，在引射器中一次流氫氣與二次回流氫氣發(fā)生混合，進(jìn)入燃料電池陽極入口。陽極出口的回流氫氣需經(jīng)過氣水分離器，盡可能消除液態(tài)水對引射器與陽極入口氫氣的影響。四噴嘴氫氣引射器通過設(shè)計4 個不同直徑的一次進(jìn)流噴嘴，調(diào)節(jié)一次進(jìn)流噴嘴處電磁閥門V1-V4 的開閉，實現(xiàn)不同的引射效果。

圖2 燃料電池氫氣循環(huán)子系統(tǒng)Fig.2 Schematic diagram of hydrogen recycle subsystems of fuel cells

2.3 引射器結(jié)構(gòu)設(shè)計

本工作基于150 kW 燃料電池電堆對引射器進(jìn)行結(jié)構(gòu)設(shè)計，電堆的基本參數(shù)如表1 所示。

表1 燃料電池工作參數(shù)Table 1 Operation parameters of the fuel cell

單噴嘴引射器的初始結(jié)構(gòu)設(shè)計采用索科洛夫引射器的設(shè)計方法[10]，對于大功率燃料電池，3～5 噴嘴是目前研究者采用較多的選擇。為降低噴嘴復(fù)雜度，同時實現(xiàn)引射器出口流量的近線性調(diào)控，經(jīng)過前期初步研究，結(jié)合實際制造工藝，本工作選擇最主流的四噴嘴引射器開展研究。四噴嘴引射器的喉部噴嘴直徑按其面積比1:2:4:8 將單噴嘴引射器的直徑進(jìn)行分割，引射器的主要結(jié)構(gòu)參數(shù)如表2所示。

表2 引射器結(jié)構(gòu)參數(shù)Table 2 Structural parameters of the ejector

3 四噴嘴引射器數(shù)學(xué)模型

基于圖2(a)，利用燃料電池系統(tǒng)仿真軟件gFUELCELL 平臺建立了PEMFC 系統(tǒng)模型。模型中包含了電堆模型、氫氣子系統(tǒng)模型、空氣子系統(tǒng)模型以及冷卻子系統(tǒng)模型。在引射器的4個噴嘴前端均設(shè)置了電磁閥。通過對4 個電磁閥開關(guān)的控制，完成各個噴嘴的開閉操作，實現(xiàn)引射器噴嘴工作模式的變化，以適應(yīng)電堆加載過程中的工作需求。圖3 為四噴嘴引射器模型的燃料電池系統(tǒng)模型流程圖。

圖3 包含四噴嘴引射器模型的燃料電池系統(tǒng)模型Fig.3 Schematic diagram of the fuel cell system model with a four nozzle ejector

引射器分為4 個區(qū)域：(1) 噴嘴區(qū)域；(2) 引射區(qū)域；(3) 混合段區(qū)域和(4) 擴(kuò)散區(qū)域。同時，對引射器模型做出了假設(shè)：(1) 模型遵循理想氣體狀態(tài)方程；(2) 進(jìn)入引射器的水均為氣態(tài)水；(3) 不考慮徑向上的溫度和壓力變化；(4) 引射器整體絕熱。計算時，對部分影響因素進(jìn)行簡化，相關(guān)公式采用一定的經(jīng)驗公式進(jìn)行描述[1]。

3.1 噴嘴段區(qū)域計算方程

當(dāng)流體為亞音速流動時，噴嘴入口的質(zhì)量流量和噴嘴喉部馬赫數(shù)計算如下：

離開噴嘴喉部的一次流流體，會在吸入室內(nèi)發(fā)生膨脹，該膨脹的極限壓力等于二次流壓力。通過等熵過程進(jìn)行計算離開噴嘴后射流體的馬赫數(shù)、溫度、速度和射流流體離開噴嘴后的直徑。

3.2 引射段區(qū)域計算方程

進(jìn)入混合段之前，引射流體質(zhì)量流量、速度、馬赫數(shù)和溫度計算如下：

其中，混合段速度分布系數(shù)，利用經(jīng)驗方程進(jìn)行計算：

3.3 混合段區(qū)域計算方程

射流流體和引射流體進(jìn)入混合段后發(fā)生混合。根據(jù)動量守恒方程：

混合段的能量守恒方程如下：

其中，混合段的混合流體速度vmix計算方程如下：

3.4 擴(kuò)散段區(qū)域計算方程

擴(kuò)散段的作用主要是增大流通截面積、提升引射器出口壓力。在擴(kuò)散段中，由等熵過程相關(guān)方程得出出口壓力，其計算式為：

對于引射器本身，其能量守恒方程如下：

為了評價引射器的性能，引入引射比、氫氣化學(xué)計量比2 個重要參數(shù)：

4 結(jié)果與討論

4.1 四噴嘴引射器控制模式對引射性能的影響

根據(jù)4 個噴嘴的閥門開關(guān)策略不同，共有15 種控制模式，根據(jù)噴嘴段一次流量公式，可以看出一次質(zhì)量流量與噴嘴喉部面積呈正比。因此對于四噴嘴引射器，以四噴嘴全開的最大供氫流量占比為100%來計算，可以得到每種控制模式下的最大供氫流量占比。每個控制模式下的噴嘴開關(guān)策略及供氫流量占比如表3 所示。

表3 四噴嘴引射器控制模式表Table 3 Control modes of the four nozzle ejector

從圖4(a)可以看出，單一噴嘴模式下引射器能起到引射能力的工作區(qū)間有一定的局限性，其引射比在工作區(qū)間內(nèi)呈先上升后下降的趨勢。在相對較低的電流密度下，引射比幾乎為0，這是因為當(dāng)運行電流密度相對較小時，陽極側(cè)氫氣消耗量也較小，氫氣源出口的壓力調(diào)節(jié)閥會通過減小閥門開度，降低引射器一次流進(jìn)口壓力，從而保持電堆陽極進(jìn)口處氫氣壓力的穩(wěn)定，與此同時，過低的一次流進(jìn)口壓力導(dǎo)致射流流體速度降低，在吸入室產(chǎn)生的負(fù)壓較小，最終導(dǎo)致在低功率情況下，引射器失效。在其工作區(qū)間內(nèi)，隨著電堆電流密度的增大，陽極側(cè)氫氣消耗量隨之上升，氫氣出口壓力調(diào)節(jié)閥逐漸增大閥門開度，一次流質(zhì)量流量與壓力不斷增加，此時吸入室產(chǎn)生的負(fù)壓不斷增大，引射比不斷上升到達(dá)峰值。

圖4 不同控制模式對引射性能的影響Fig.4 Effects of different control modes on ejector performance

圖4(b)為不同控制模式下陽極側(cè)進(jìn)堆氫氣壓力隨電流密度變化的曲線，對于單一噴嘴控制模式而言，當(dāng)電流密度進(jìn)一步增大，氫氣需求量逐漸到達(dá)該引射器控制模式的上限，導(dǎo)致其無法維持電堆所需的氫氣流量及壓力，引射性能逐漸下降。分析原因是，當(dāng)噴嘴直徑變小，其內(nèi)部流阻變大，若要引射同樣流量的氣體，需要更高的射流入口壓力，而射流入口壓力取決于比例閥特性以及其氫氣前端壓力值。因此，當(dāng)系統(tǒng)中比例閥特性與前端壓力值一定時，就必然決定了固定結(jié)構(gòu)式引射器在小噴嘴直徑下無法滿足高功率的需求。

對比圖4(a)與(c)，可以看出在單一噴嘴模式下，氫氣計量比和引射比的變化趨勢基本相同。引射比增加，更多的氫氣通過引射器再循環(huán)到陽極，氫氣計量比隨之增加。

4.2 四噴嘴引射器在大功率PEMFC系統(tǒng)控制策略

單一噴嘴模式的工作范圍有限，雖然小噴嘴直徑模式可以滿足低電流運行的要求，但在高電流區(qū)間會失效。而大噴嘴直徑模式則相反，其在低電流區(qū)間幾乎不能產(chǎn)生引射。因此，根據(jù)上述特性，結(jié)合15 種控制模式下的進(jìn)堆壓力值與氫氣計量比變化趨勢，得到四噴嘴引射器在燃料電池堆全功率區(qū)間變流量穩(wěn)壓供氫控制策略，如表4 所示。

表4 不同電流區(qū)間對應(yīng)的四噴嘴引射器控制策略Table 4 Control Strategies for the four nozzle ejectors at different current ranges

4.2.1 優(yōu)化控制策略下全功率區(qū)間的四噴嘴引射器與燃料電池性能

根據(jù)表4，對全功率區(qū)間多噴嘴引射器系統(tǒng)控制進(jìn)行仿真。從圖5(a)可以看到，在對四噴嘴引射器進(jìn)行控制策略優(yōu)化后，可以實現(xiàn)0.11～2.11 A·cm-2電流密度范圍內(nèi)的引射，與單噴嘴引射器1.30～1.86 A·cm-2的作用區(qū)間相比，其運行工況覆蓋范圍擴(kuò)大至5%～100%。

圖5 四噴嘴引射器控制策略下的性能變化曲線Fig.5 Performance profiles under the control strategy of four nozzle ejector

圖5(b)為優(yōu)化控制策略下電流拉載過程中氫氣計量比變化曲線，可以看到隨著拉載電流的增加，氫氣計量比逐漸下降，在電流密度為0.1～0.5 A·cm-2時下降較快，分析原因是在該區(qū)間內(nèi)引射比呈上升趨勢，增加的引射比意味著更多的二次回流氫氣，而隨著拉載電流的增大，陽極的氮氣含量累積量也在快速增加，從而導(dǎo)致經(jīng)引射器送入陽極的氫氣中氮氣含量增加，氫氣計量比快速下降。其中，陽極側(cè)氮氣的產(chǎn)生是由于質(zhì)子交換膜間發(fā)生了氣體滲透現(xiàn)象。原理上，質(zhì)子交換膜中反應(yīng)氣體是不會發(fā)生滲透的，防止反應(yīng)物組分在電化學(xué)反應(yīng)發(fā)生前混合。然而，由于膜的多孔結(jié)構(gòu)，水含量以及氫和氧在水中的可溶解性等原因，PEM 之間會發(fā)生微量氣體滲透過程。圖5(c)為優(yōu)化控制策略下電流拉載過程中陽極入口氫氣壓力變化曲線，可以看到在電流拉載的過程中，陽極入口壓力控制在250 kPa 的目標(biāo)值，其中壓力的較大波動是由噴嘴控制模式切換導(dǎo)致，小幅波動受排氫閥間歇排氫的影響。圖5(b)為優(yōu)化策略下電流拉載過程中電堆性能曲線，可以看到隨著電流密度的增加，電堆效率與單片電壓逐漸降低，電堆功率逐漸上升。

圖6 是電流拉載時，陰極/陽極氮氣累積量與氮氣摩爾通量變化曲線。在電流密度小于1.0 A·cm-2，PEM 氮氣摩爾通量為負(fù)值，氮氣從陰極向陽極擴(kuò)散，在陰極入口流量及排氫策略不變的情況下，陰極氮氣含量不斷下降，陽極氮氣含量不斷上升逐漸大于陰極，當(dāng)電流密度大于1.0 A·cm-2時，PEM 氮氣摩爾通量為正值，氮氣開始從陽極向陰極擴(kuò)散。

圖6 電流拉載過程中，陰極/陽極氮氣累積量與PEM 氮氣摩爾通量變化曲線Fig.6 Cathode/anode nitrogen molar accumulation as a function of PEM nitrogen molar flux diffusion during the current pulling process

4.2.2 優(yōu)化控制策略下操作條件對四噴嘴引射器引射性能的影響

4.2.2.1 出口壓力的影響

圖7 為電流拉載過程中，引射器出口壓力變化對引射性能的影響?？梢钥吹?，引射器出口氣體壓力變化對噴嘴模式切換波動幾乎無影響，隨著引射器出口壓力的增大，引射比、氫氣計量比和出口質(zhì)量流量均增大，且輸出電流越大，引射器出口壓力對出口質(zhì)量流量的影響越明顯。

圖7 引射器出口壓力變化對引射性能的影響Fig.7 Effects of outlet pressure variation on ejector performance

4.2.2.2 排氫頻率的影響

圖8 為排氫頻率對引射性能的影響，可以看到，減小排氫頻率可以提升引射器引射比和出口質(zhì)量流量，但同時增加了氫氣回路內(nèi)的氮氣累積量，氫氣計量比降低。而增大排氫頻率會使噴嘴模式切換帶來的出口質(zhì)量流量波動更加明顯，出口壓力波動也較為劇烈。因此，在實際系統(tǒng)運行過程中，需要針對不同電流工況設(shè)置合適的排氫頻率。

圖8 氫氣路尾排頻率對引射性能的影響Fig.8 Effects of exhaust frequency on the eject performance

5 結(jié) 論

對四噴嘴氫氣引射器在大功率PEMFC 氫氣子系統(tǒng)的性能開展建模與研究，采用集總參數(shù)建模方法對四噴嘴引射器進(jìn)行建模，研究了四噴嘴引射器控制模式對引射比的影響，同時優(yōu)化噴嘴控制策略，進(jìn)行全功率區(qū)間氫氣子系統(tǒng)仿真，得到以下結(jié)論：

(1) 每一種噴嘴控制模式下，引射器具有引射能力時，電堆的運行電流區(qū)間有一定的局限性，在引射作用區(qū)間，隨著電堆運行電流的增加，其引射比先升高后下降，存在一個最大值。

(2) 混合段直徑不變，隨著噴嘴喉部面積減小，引射器的作用區(qū)間逐漸向低電流工況移動，同時整體的引射比升高，在作用區(qū)間內(nèi)的引射性能變化也更加劇烈。

(3) 提出了四噴嘴引射器在不同電流區(qū)間的閥門控制策略，使其能夠滿足功率為8～153 kW 的電堆穩(wěn)壓供氫需求。

(4) 引射器出口壓力增大，引射比和出口質(zhì)量流量均增大，引射性能提升。

符號說明：

D — 直徑，mmw — 質(zhì)量分?jǐn)?shù)，%F — 法拉第常數(shù)，C·mol-1 μ — 引射比/引射系數(shù)h — 比焓，J·kg-1，h=CpT(Cp 為比定壓熱容) ρ — 密度，kg·m-3 I — 電流，A 上標(biāo)k — 絕熱指數(shù)，理想氣體k=Cp/Cv(Cv 為比定容熱容) in — 引射器噴嘴進(jìn)口M — 摩爾質(zhì)量，kg·mol-1 out — 引射器出口Mr — 相對分子質(zhì)量，kg·mol-1 下標(biāo)Ma — 馬赫數(shù) i — 噴嘴數(shù)量nv — 混合段速度分布系數(shù) jet — 離開噴嘴后的射流流體p — 壓力，kPa mix — 混合段qm — 質(zhì)量流量，kg·s-1 nozzle — 噴嘴喉部R — 摩爾氣體常數(shù)，J·mol-1·K-1 p — 一次流流體/射流流體S — 化學(xué)計量比 R，C — 臨界比值T — 溫度，K s — 二次流流體/引射流體v — 速度，m·s-1 s，2 — 進(jìn)入混合段的引射流體