本文亮點：1、建立了110kW的大功率燃料電池模型，模型經(jīng)試驗驗證，可靠性高 2、探究不同運行參數(shù)對燃料電池一致性的影響研究，為后續(xù)燃料電池熱管理提供指導
摘 要 隨著燃料電池堆朝著大功率發(fā)展，其工作時單體間的不一致性更加明顯，長時間處于惡劣工作條件的單體壽命會明顯短于其他單體，并導致電池堆的壽命大幅縮減。為探究不同運行參數(shù)對大功率燃料電池單體工作性能一致性的影響規(guī)律，首先，建立了包括流體網(wǎng)絡模型、燃料電池電壓模型和燃料電池熱阻模型三個部分的110 kW大功率燃料電池模型。其次，開展燃料電池穩(wěn)態(tài)試驗，對所建立的燃料電池模型進行試驗驗證，仿真與試驗結(jié)果誤差在5%以內(nèi)。最后，基于模型仿真，以電壓最大偏差率為評價指標，分別探究工作電流、冷卻水流量和冷卻水進口溫度三個運行參數(shù)對燃料電池電壓一致性的影響規(guī)律。仿真結(jié)果表明，工作電流對燃料電池單體電壓一致性的影響程度更大，其次是冷卻水進口溫度，最后是冷卻水流量。本研究有助于大功率燃料電池發(fā)動機的結(jié)構(gòu)優(yōu)化設計以及為熱管理控制策略開發(fā)提供指導。
關(guān)鍵詞 燃料電池；電壓一致性；最大偏差率
隨著人類對環(huán)境保護的重視程度提高，一系列政策的制定推動著燃料電池普及與應用。其中質(zhì)子交換膜燃料電池，憑借其功率密度高、體積小以及工作溫度低等特點，廣泛運用于燃料電池車以及儲能場景中。然而燃料電池對工作溫度要求苛刻，長期處于過高或過低溫度工作的電池，其壽命將顯著短于其他工作在正常溫度的電池，這也會導致燃料電池堆的“短板效應”顯現(xiàn)，使得電池堆的壽命將由壽命最短的單體決定。燃料電池堆在工作時由于結(jié)構(gòu)布置、運行狀態(tài)和冷卻方式不同等原因，存在溫度分布不均勻的問題。并且隨著車用燃料電池的應用，燃料電池也將朝著大功率發(fā)展，大功率燃料電池的熱分布特性將變得不同，局部一致性變差的問題更加凸顯。
目前對于燃料電池單體一致性的研究主要采用仿真方式進行。現(xiàn)有常見的燃料電池堆仿真研究方法包括建立燃料電池堆的一維模型、二維模型和三維模型。Baschuk等建立了含50個單體電池的流體網(wǎng)絡模型，仿真結(jié)果表明，通過增加歧管的尺寸，減少每個雙極板的氣體流動通道的數(shù)量，以及合理地改變氣體流動通道中的流阻，有助于燃料電池堆內(nèi)反應物的均勻分布。趙巖等采用燃料電池流道簡化方法和壓力速度耦合算法數(shù)學模型對30~50節(jié)燃料電池進行仿真，仿真結(jié)果表明采用Z形結(jié)構(gòu)的空氣分布均勻性比U形好，并且增加歧管寬度有利于改善氣體分布。Yang等基于Z形歧管，考慮每個燃料電池內(nèi)部的流量分布以及反應、相變和傳輸過程，建立了包含5節(jié)單體的1+1維多相電堆子模型和流量分布綜合電堆模型，仿真結(jié)果表明，增加入口壓力對反應物分布均勻性的影響不大，但它提高了電堆工作性能，并且采用更大的歧管橫截面積能夠使電池之間的反應物分布更均勻。同時Yang等建立了考慮了摩擦和局部壓降損失以及膜水、液態(tài)水、氣體組分和熱量的耦合傳輸過程的二維電堆模型，通過仿真發(fā)現(xiàn)如果歧管的橫截面積相對較小，U形結(jié)構(gòu)可以比Z形結(jié)構(gòu)獲得更好的反應物分布均勻性，這有助于減小整個電堆的幾何尺寸。Asgharian等建立了考慮雙極板流道的歧管10節(jié)燃料電池三維模型，通過CFD仿真，結(jié)果表明歧管的直徑增加能夠使流道中氣體分布更均勻。Huang等建立含100、200、300個單體燃料電池的歧管三維模型，基于CFD仿真的手段，結(jié)果表明，隨著燃料電池個數(shù)增加和工作電流增大，歧管流量分布更加不均勻。
目前針對燃料電池的建模仿真，大部分研究只考慮了燃料電池氣體歧管結(jié)構(gòu)、電壓或溫度其中某一個或者某幾個因素。針對燃料電池堆性能一致性的研究，多數(shù)研究主要聚焦在流量分布上，而電池堆實際工作時，氣體流量、電壓和溫度之間存在耦合關(guān)系，相互影響。并且現(xiàn)有研究的燃料電池存在功率小的問題，而針對大功率電池堆性能一致性的變化規(guī)律尚未明確。
因此針對以上存在的問題，本工作以某款110 kW大功率燃料電池為研究對象，建立考慮流量、電壓和溫度分布的燃料電池多參數(shù)耦合模型，通過燃料電池穩(wěn)態(tài)試驗對模型進行驗證，并且基于該模型仿真結(jié)果探究不同運行參數(shù)對燃料電池單體電壓一致性的影響規(guī)律。本研究旨在提高燃料電池工作時單體性能的一致性并提高燃料電池的使用壽命，為后續(xù)燃料電池堆結(jié)構(gòu)優(yōu)化設計和熱管理策略設計提供指導。
1 燃料電池理論模型建立
本工作所建立的燃料電池模型包括：流體網(wǎng)絡模型、燃料電池電壓模型和燃料電池熱阻模型三部分。首先將燃料電池初始溫度以及總的氣體流量代入流體網(wǎng)絡模型計算并更新單體氣體分布結(jié)果，其次將單體氣體分布結(jié)果代入燃料電池電壓模型計算并更新單體電壓結(jié)果，最后將電壓結(jié)果代入燃料電池熱阻模型計算并更新溫度，而溫度也會影響氣體流量和電壓結(jié)果，因此需要進行迭代計算，并且在模型當中額外考慮了流道結(jié)構(gòu)的局部損失、冷卻水流量分布和燃料電池接觸熱阻等因素的影響。
1.1 流體網(wǎng)絡模型
本工作所研究的110 kW大功率燃料電池堆的陽極歧管、陰極歧管和冷卻水歧管布置均采用“U形”結(jié)構(gòu)布置即流體的進口和出口在同一側(cè)，其中陽極歧管和陰極歧管進出口均在電堆的右側(cè)，冷卻水進出口在電堆的左側(cè)，如圖1所示。
圖1 大功率燃料電池堆歧管布置方式
流體網(wǎng)絡模型主要用于計算陽極氣體、陰極氣體以及冷卻水的流量分布，由于計算流程和內(nèi)容相似，本工作以陽極和陰極氣體為例。如圖2所示的“U形”流體網(wǎng)絡圖，反應氣體從進氣歧管進入，然后分配到每一個電池，最后從排氣歧管流出，在流動過程中反應物的物性參數(shù)會受到單體溫度的影響，同時遵循能量守恒定律和質(zhì)量守恒定律，通過計算得到反應氣體在每塊電池中的質(zhì)量流量，為后續(xù)構(gòu)建燃料電池電壓模型提供輸入。
圖2 “U形”燃料電池堆流體網(wǎng)絡圖
1.2 燃料電池電壓模型
1.3 燃料電池熱阻模型
通過燃料電池電壓模型得到每塊電池的電壓以后，即可計算出每塊電池的產(chǎn)熱，本研究建立了帶有散熱模塊的電池堆熱阻模型。散熱形式是在相鄰兩個單電池中間加入一個冷卻板。
燃料電池主要產(chǎn)熱的區(qū)域是在陰極催化層內(nèi)，H+通過之間的交換膜，到達陰極催化層與O2發(fā)生還原反應，產(chǎn)生由于反應熵變引起的可逆熱，燃料電池堆產(chǎn)生的熱量Q為：
圖3 燃料電池熱阻模型示意圖
燃料電池的接觸熱阻主要有四個部分：氣體擴散層和雙極板之間的接觸熱阻、氣體擴散層和微孔層之間的接觸熱阻、微孔層和催化層之間的接觸熱阻以及催化層和質(zhì)子交換膜之間的接觸熱阻。其中氣體擴散層和微孔層之間的熱阻常用在微孔層和擴散層之間增加一層微孔層-擴散層集成區(qū)來代替。由于催化層與氣體擴散層和質(zhì)子交換膜接觸好、接觸熱阻較小，計算時可以忽略這部分熱阻。
同時根據(jù)流體網(wǎng)絡模型計算得到單體冷卻水流量分布，可以計算單體散熱量。燃料電池主要通過冷卻水與雙極板對流換熱來帶走燃料電池產(chǎn)生的熱量，其帶走的熱量為：
根據(jù)公式(17)~(21)相應的熱阻以及產(chǎn)熱和散熱量對每塊單體建立熱平衡方程，通過求解熱平衡方程組，即可得到每塊電池的內(nèi)部溫度分布結(jié)果。
1.4 計算流程
通過上文可以得出，電池堆的溫度分布將會影響單個電池的性能以及反應物的一些物性參數(shù)，而單體電池之間不同輸出性能下產(chǎn)生的熱量又會影響整個電池堆的溫度分布，所以燃料電池堆模型是需要迭代計算的。將所建立的燃料電池模型在MATLAB中計算，其計算流程如圖4所示。
圖4 燃料電池模型計算流程圖
（1）在初始計算中，輸入氫氣、空氣、冷卻水進口總流量和溫度等運行參數(shù)以及歧管尺寸、單體數(shù)量等結(jié)構(gòu)參數(shù)，并初始化燃料電池單體溫度、氫氣流量、空氣流量和冷卻液流量分布；
（2）根據(jù)單體溫度、氫氣流量、空氣流量和冷卻液流量分布，通過流體網(wǎng)絡模型計算并更新每塊電池的氫氣流量、空氣流量和冷卻液流量分布；
（3）將單體的氫氣流量、空氣流量以及溫度分布結(jié)果代入燃料電池電壓模型，計算并更新每塊電池的電壓分布；
（4）根據(jù)燃料電池電壓和冷卻水流量分布，分別計算單體產(chǎn)熱量和散熱量，通過燃料電池熱阻模型，計算并更新每塊電池的溫度分布；
（5）比較本次溫度計算結(jié)果與上次溫度結(jié)果的差值是否小于誤差要求，如果大于誤差，則將本次計算的溫度分布代入步驟(2)中的流體網(wǎng)絡模型重新進行計算；如果滿足要求，則輸出結(jié)果。
2 燃料電池模型仿真結(jié)果與試驗驗證
本試驗對象選用某款110 kW大功率燃料電池并基于開發(fā)的燃料電池熱管理測試系統(tǒng)進行試驗，燃料電池部分參數(shù)如表1所示。
表1 燃料電池部分參數(shù)
根據(jù)GB/T 24554—2022《燃料電池發(fā)動機性能試驗方法》，首先對燃料電池進行熱機過程預處理，然后進行穩(wěn)態(tài)特性試驗，試驗過程功率和工作電壓隨電流的變化和部分運行參數(shù)隨電流的變化如圖5所示。
圖5 部分運行參數(shù)隨電流的變化
將圖5所示的燃料電池相關(guān)運行參數(shù)代入所建立的燃料電池模型中，通過MATLAB進行計算，可以得到每塊單體電池的氣體流量、電壓、溫度等結(jié)果。由于在試驗中單體電壓和溫度結(jié)果較氣體流量方便獲取，本工作將從試驗中獲得單體電壓和溫度進行仿真模型有效性驗證。
其中不同電流下燃料電池單體平均電壓的仿真與試驗結(jié)果如圖6所示，可以看出隨著工作電流的增大，燃料電池單體平均電壓減小，并且仿真和試驗結(jié)果誤差在3%以內(nèi)。
圖6 燃料電池單體平均電壓的仿真和試驗結(jié)果對比
按照離氣體進口的距離從近到遠對燃料電池單體進行編號。將試驗臺架記錄的電池單體電壓與仿真電壓進行對比，結(jié)果如圖7所示，燃料電池單體電壓仿真與試驗結(jié)果的誤差在5%以內(nèi)，證明了模型可靠性較高。同時隨著工作電流的增大，燃料電池單體的電壓變小，燃料電池單體電壓的不一致性越來越明顯，并且從整體上看，靠近氣體進口側(cè)的電壓比遠離氣體進口側(cè)的電壓高，這符合“U形”進氣結(jié)構(gòu)的特點，并且靠近燃料電池兩端的單體電壓會低于其他單體電壓。
圖7 不同工況下燃料電池單體仿真電壓與試驗電壓對比
在試驗時，一般把冷卻水出口位置的溫度作為燃料電池堆的平均溫度。將仿真計算得到的每塊單體電池溫度進行平均處理得到電池堆的平均溫度，其結(jié)果如圖8所示，可以看出隨著工作電流的增大，電池堆的平均溫度也在增加，并且仿真平均溫度和試驗冷卻水出口溫度誤差在2%以內(nèi)。
圖8 燃料電池仿真平均溫度和冷卻水出口溫度對比
同時通過試驗測試平臺配備的紅外測溫設備可以獲得燃料電池工作時表面的溫度數(shù)據(jù)。通過對獲得的溫度數(shù)據(jù)進行處理和修正，得到如圖9所示的結(jié)果，絕大部分仿真誤差在5%以內(nèi)，證明了模型可靠性較高。同時從圖中可以看出，燃料電池單體溫度隨著工作電流增大而增加，并且溫度分布變得更不均勻，從整體上看靠近氣體進口側(cè)的溫度比遠離氣體進口側(cè)的溫度高，并且靠近兩端的單體溫度要低于平均溫度，同時工作電流較大時，兩端單體溫度會明顯低于平均溫度。這與單體電壓的變化規(guī)律相似。
圖9 不同工況下燃料電池單體仿真溫度與紅外溫度對比
3 運行參數(shù)對單體性能一致性的影響
從上述試驗結(jié)果和仿真結(jié)果可以看出，燃料電池在工作時，單體性能存在不一致性，而在試驗過程中電流變化的同時，冷卻水溫度與流量等參數(shù)也在發(fā)生變化。為更好地探究不同參數(shù)對燃料電池單體性能一致性的影響規(guī)律，本工作以660 A的試驗工況（工作電流為660 A，冷卻水流量為230 kg/min，冷卻水進口溫度為72 ℃）為基準，分別探究電流、冷卻水流量和冷卻水進口溫度三個運行參數(shù)對單體性能一致性的影響。為使仿真更貼近實際使用情況，保證燃料電池工作在正常的溫度區(qū)間（65～85 ℃），本工作探究的燃料電池工作電流為60～780 A，冷卻水流量為120～420 kg/min，冷卻水進口溫度為50～80 ℃。
通過所建立的燃料電池模型進行仿真。由于燃料電池堆壽命取決于性能最差的單體，并且電壓可以反映燃料電池的工作性能，所以有必要對單體電池電壓最大值和最小值與平均值的偏離程度進行監(jiān)測，本工作以電壓最大偏差率圖片來評價燃料電池單體的一致性，其中電壓最大偏差率圖片計算公式為：
3.1 電流對單體電壓一致性的影響
當電流作為唯一變量時，探究電流從60 A到780 A時燃料電池電壓一致性的變化規(guī)律，其仿真結(jié)果如圖10所示。可以看出隨著電流增大，燃料電池的電壓減小，當燃料電池在較小的工作電流下時（300 A以下），單體電壓一致性隨著電流增大變化不明顯，但是處于大工作電流下時（300 A以上），單體電壓一致性隨著電流增大而變差，相較于300 A時，電流每增大10%，最大偏差率平均增大13%。
圖10 電流為唯一變量的仿真結(jié)果
究其原因為當工作電流為唯一變量時，燃料電池單體散熱基本不變，燃料電池堆的溫度主要由產(chǎn)熱決定，而產(chǎn)熱跟工作電流有關(guān)，工作電流增加則會導致燃料電池的歐姆損失和活化損失增大，如圖11所示，可以得出歐姆損失和活化損失以及對應的極差均隨著電流的增大而變大，但是歐姆損失的極差變化幅度隨著電流的增大而增大，而活化損失的極差變化幅度隨著電流的增大而減小，當工作電流從60 A變化到300 A，歐姆損失極差增加了1 mV，活化損失極差增加了0.5 mV；當工作電流從300 A變化到720 A，歐姆損失極差增加了5 mV，活化損失極差增加了0.3 mV。而歐姆損失主要在中等電流密度（0.15~1.2 A/cm2）下起作用，這也是導致當燃料電池在300 A電流以下時單體電壓不一致性的變化幅度隨著工作電流的增大而減小，在300 A電流以上時，單體電壓不一致性的變化幅度隨著工作電流的增大而增加的原因。因此單體電壓的一致性將由歐姆損失和活化損失決定，同時工作電流對單體歐姆損失的影響較單體活化損失更大，因此這也使得燃料電池的工作性能一致性隨工作電流的升高而降低。
圖11 電流為唯一變量的歐姆和活化損失極差
3.2 冷卻水流量對單體電壓一致性的影響
當冷卻水流量作為唯一變量時，探究流量從120 kg/min到420 kg/min時燃料電池電壓一致性的變化規(guī)律，其仿真結(jié)果如圖12所示?？梢钥闯?，隨著冷卻水流量增大，燃料電池的電壓變小，單體電壓一致性變好，相較于420 kg/min，流量每減少10%，最大偏差率平均增大8%。但是當流量增大到一定程度時，此時再增大流量對單體電壓一致性改善不大。
圖12 冷卻水流量為唯一變量的仿真結(jié)果
究其原因為當冷卻水流量為唯一變量時，此時燃料電池產(chǎn)熱一定，燃料電池堆的溫度主要由散熱決定，而散熱主要與冷卻水流量有關(guān)，呈反比關(guān)系。當流量大到一定程度時，對降低燃料電池的溫升作用不明顯，但是會額外增加泵的能耗。而單體溫度的一致性則會影響單體電壓的一致性，由于冷卻水流量增大，燃料電池堆的溫度降低，則會導致燃料電池的歐姆損失增大和活化損失增大，但是會使單體間溫差變小，使得歐姆損失和活化損失的極差變小，如圖13所示，當冷卻水流量從120 kg/min變化到420 kg/min時，歐姆損失極差降低6 mV，活化損失極差降低了1 mV，因此燃料電池工作性能的一致性隨著冷卻水流量增加而提高。因此在實際工作時，冷卻水流量的選擇應該結(jié)合工作電流、系統(tǒng)能耗以及燃料電池的工作溫度區(qū)間進行確定。
圖13 冷卻水流量為唯一變量的歐姆和活化損失極差
3.3 冷卻水進口溫度對單體電壓一致性的影響
當冷卻水進口溫度作為唯一變量時，探究溫度從50 ℃到80 ℃時燃料電池電壓一致性的變化規(guī)律，其仿真結(jié)果如圖14所示。隨著冷卻水進口溫度增大，燃料電池的電壓上升，雖然燃料電池的溫升沒有明顯變化，但是單體電壓一致性變好，相較于80 ℃，進口溫度每減少10%，最大偏差率平均增大20%。
圖14 冷卻水進口溫度為唯一變量的仿真結(jié)果
究其原因為當冷卻水進口溫度為唯一變量時，燃料電池堆的產(chǎn)熱和散熱基本不變，溫升基本不變，因此冷卻水進口溫度將決定燃料電池堆穩(wěn)定后的溫度，這將會影響燃料電池的電壓。隨著冷卻水進口溫度降低，燃料電池堆的溫度降低，燃料電池的歐姆損失和活化損失增大，當冷卻水進口溫度從50 ℃變化到80 ℃，歐姆損失極差降低3 mV，活化損失極差提高0.9 mV。所以單體工作性能不一致性將由歐姆損失和活化損失決定，雖然單體活化損失的一致性隨著冷卻水進口溫度升高而降低，但是由于單體歐姆損失的一致性隨著冷卻水進口溫度升高而提高，并且較單體活化損失，冷卻水進口溫度對歐姆損失的一致性影響更大，因此這也使得燃料電池工作性能一致性隨冷卻水進口溫度升高而提高。雖然冷卻水進口溫度升高，可以改善單體性能以及提高其一致性，但是這會使燃料電池的工作溫度過高，超過正常的工作溫度區(qū)間，反而導致質(zhì)子交換膜的破壞，進而使得燃料電池壽命減短。因此在實際工作時，冷卻水進口溫度應當合理選擇，一般控制在65～70 ℃。
4 結(jié) 論
本工作建立了110 kW大功率燃料電池模型，通過試驗對模型進行驗證，并基于該模型探究了電流、冷卻水流量和冷卻水進口溫度對單體電壓一致性的影響，主要獲得了以下幾點結(jié)論：
圖15 冷卻水進口溫度為唯一變量的歐姆和活化損失極差
（1）建立燃料電池多參數(shù)耦合模型，并進行燃料電池穩(wěn)態(tài)試驗驗證，仿真結(jié)果與試驗結(jié)果對比，平均電壓和溫度誤差在5%以內(nèi)，模型精度較高；
（2）基于試驗與仿真的單體電壓和溫度分布結(jié)果，單體電壓和溫度的分布規(guī)律相似，整體上呈現(xiàn)靠近氣體進口側(cè)高，遠離氣體進口側(cè)低的特點，但是靠近燃料電池堆兩端的單體電壓和溫度會明顯低于其他單體；
（3）工作電流、冷卻水流量和冷卻水進口溫度的變化主要影響了電池堆的溫度，導致歐姆損失和活化損失發(fā)生改變，而其中活化損失受到的影響更大，進而使單體電壓一致性發(fā)生改變；
（4）在探究運行參數(shù)對單體電壓一致性影響的研究中，按影響程度的大小排序依次為工作電流、冷卻水進口溫度、冷卻水流量。
