氫燃料電池極板S 型流道結構的仿真優化

劉利娜，饒國燃，張冬潔

（北京理工大學珠海學院工業自動化學院，珠海 519088）

引言

在當前能源和環境問題日益嚴峻的局面下，氫能作為一種清潔高效的能源材料，逐步受到人們的關注。作為氫能利用的手段之一，質子交換膜燃料電池因具有比功率高、能量轉換效率高、操作溫度低、耐久性較好等優點而備受青睞[1,2]。在現行商業化的燃料電池設計思路中，膜電極和雙極板是組成燃料電池的兩大核心組件，其中雙極板承擔著輸送反應氣體、排出反應產物、集流、支撐的作用[3]。雙極板的流場設計十分重要，合理的流場設計不僅要保證燃料和氧化劑供應充足且分配均勻，還要把反應生成的熱和水及時排出，避免發生超溫或水淹[4]，維持燃料電池正常運行。

有關雙極板流場的研究使用較多的有直型流道流場、蛇形流場和交指形流場，研究者一般采用數值模擬方法研究岸寬、流道深度和寬度等形貌尺寸的影響。例如，沈俊等[5]采用Fluent 軟件分別研究了在流道中加凸臺和擴散層開孔兩種方法對強化電池內部傳質的效果，得出兩種方法均能強化電池內部傳質，提高電池性能的結論。Ge[6]等提出了一種二維、穩態燃料電池模型，研究了操作條件、氣體流向、膜的厚度等對燃料電池性能的影響，并認為提高陰極壓力可改善陽極加濕的問題。Atyabi[7]等新提出了一種新型流場，并采用數值模擬方法將新型流場與傳統平直流場進行對比分析，得出了新型流場性能優于傳統流場的結論。

本文基于以上研究內容，以S 型流道為研究對象，通過解構S 型流道結構，采用數值模擬方法從全流場層面建立了單流道燃料電池的三維模型，分析了S 型流道主要設計參數對電池性能的影響，同時考慮到陰極是水生成的主要區域，而單電池組裝時陰陽極不同對應方式會造成陰極流場結構的不同，故進一步研究了S 型流道對應方式對電池性能的影響，從而為雙極板流場結構改進與優化提供參考。

1 模型和方法

1.1 基于FLUENT 的燃料電池物理化學過程建模

質子交換膜燃料電池內部的反應是發生在多體之間的并由多物理場參與的復雜的電化學反應。一個完整的三維全尺寸燃料電池多體模型包括了陽極集流板（AC）、陽極流道（AH）、陽極擴散層（GDL）、陽極催化層（CAT）、質子交換膜（MEM）、陰極催化層（CAT）、陰極擴散層（GDL）、陰極流道（CH）和陰極集流板（CC）等九個部分，其結構示意圖如圖1 所示。Fluent 軟件中的PEMFC 模塊可以針對全尺寸燃料電池模型來進行多物理場耦合求解，其中涉及到的控制方程有質量守恒方程，動量守恒方程，能量守恒方程，組分擴散方程，電荷守恒方程和氣體擴散方程等[8-10]，具體方程如下（1）～（7）所示。通過設置合理的邊界條件并耦合求解上述方程就能得到燃料電池內部各組分的分布情況和各反應氣體的流動情況。在反應氣體流量分配均勻的情況下，極板上每條流道中的反應氣體流動、傳質與化學反應強度大體相同，為了節約計算成本，本文只選取其中一條流道進行分析。

圖1 質子交換膜燃料電池的結構示意圖

式中：

ρ、ε、P、μ、cp、keff、T、σ 和φ—分別代表密度、孔隙率、壓力、粘度、定壓比熱容、有效熱導率、溫度、電導率和電勢；

Sm、Su、SQ、Sk、Rs和Rm—分別為質量源項、動量源項、能量源項、組分源項、電子電流源項和質子電流源項；

ck、Dkeff、Dk—分別為組分k 的濃度、有效擴散系數和在多孔介質中的擴散系數。

1.2 S 型流道結構

在質子交換膜燃料電池極板流場設計中，直型流道、蛇形流道是兩類常見的基本流道型式。前者工質流動阻力小，結構簡單，但傳質效果相對較差；后者工質流動路徑復雜，阻力較大，且在大型極板流場設計中較難應用，但傳質和排水效果相對較好。S 型流道結構可看成是直型流道與蛇形流道的一種結合，整體上呈現直通型，比直型流道具有更高的湍流度（一定程度上強化了傳質），也比蛇形流道具有較小的流動阻力。此外，在金屬雙極板流場設計中，S 型流道還易于布置雙極板中間的冷卻液流動通道，例如本田在其專利中即披露了一款采用S型流道的金屬雙極板燃料電池結構。本文的研究針對S型流道展開。直型、蛇形和S 型流道的結構示意如圖2所示。

圖2 直型、蛇形和S 型極板流道結構

本文將S 型流道結構進行解構，其型線可由三個具有相切關系的圓弧描述，如圖3 所示。其中圓A、B 在同一直線上，相距L，且均與圓C 相切（圓A 和B 并不要求相切），三個圓的半徑并不要求相同，設其半徑分別為R1、R3、R2。將三段相切的弧線平移w 可得到流道寬度為w 的S 型流道。

圖3 S 型極板流道型線描述

定義r=d/L 為S 型流道型線的相對彎曲度，其中d 為彎曲度，即為S 型線底部到頂部的距離；L 為拉伸度，即為兩個圓心的距離。當R1≠R3時，d 不唯一，故本文僅考慮R1=R3的情形（事實上，R1與R3不同時的S 型流道結構效果處于兩種本文所考慮的情形之間，故可忽略），一般要求R1+R2≥L/2，此時，

雖然流道的其余參數即流道寬度和深度對流體流動特性和燃料電池性能都存在影響，但已有不少學者進行了研究，本文主要討論S 型流道拉伸度和彎曲度對工質流動、傳質特性和燃料電池性能的影響。在保證流道寬度和深度不變的前提下，通過改變彎曲度和拉伸度的值搭建不同的單流道燃料電池模型來對上述問題進行研究，為了便于模型搭建，假定R1=R2。

除了S 型流道型線之外，本文進一步討論了燃料電池膜電極兩側陰陽極板上S 型流道的對應方式對燃料電池性能的影響。在保證流道深度、寬度等條件相同的情況下，只改變陰陽極S 型流道的對應方式，建立了兩種不同的對應模型，根據陰陽極流道結構相對于膜電極是否對稱，將兩種模型分別稱為正對應模型和反對應模型，如圖4 所示。從圖中可以看出，兩種模型的陽極流道完全相同，而陰極流道彎曲方向相反。本文將通過研究兩種模型的單流道燃料電池極化曲線對上述問題進行研究。

圖4 單流道燃料電池幾何模型

正反對應模型單流道燃料電池的整體幾何尺寸為50 mm（長）×4 mm（寬）×2.66 mm（高），其中膜電極的有效活性面積為50 mm×4 mm，流道截面尺寸為1 mm×1 mm，其余幾何尺寸如表1 所示。而進行S 型流道設計參數影響研究時，電池組裝方式采用正對應方式建模，且保證流道截面尺寸相同。

表1 單流道燃料電池的幾何尺寸

邊界條件的設置對燃料電池模型的求解結果非常重要。將陰、陽極進口邊界全部設置為質量流量進口，通過改變工作電壓來得到不同工況下的電流密度，從而得到燃料電池的極化曲線。出口邊界采用壓力出口，氫氧進氣方式采用逆流，即陰陽極流道內氣體流動方向相反。各計算區域材料物性參數保持不變，主要材料和邊界條件的設置情況如表2 所示。

表2 主要材料參數與邊界條件

燃料電池的多體模型具有跨尺度的特點，具體表現在其流道長度方向的尺寸和寬深度方向的尺寸一般相差一個數量級，這樣給網格的劃分和模型的求解帶來巨大的困難。為排除網格數量對計算結果的影響，將膜電極中陰極氣體擴散層、陰極催化層、質子交換膜、陽極催化層和陽極氣體擴散層等五個區域的網格進行細分，并進行獨立性驗證，使網格數量滿足計算要求。求解方法采用SIMPLE 算法，全部方程的求解采用BCGSTAB 方法和F-cycle 循環，并調節松弛因子來加速結果收斂。

2 結果與討論

2.1 模型驗證

為驗證模型的準確性和可靠性，將模擬得到的極化曲線與Wei Yuan 的實驗數據進行對比[10]，如圖5 所示。在進行模擬計算時所用各種材料的物性參數、邊界條件等均參照文獻進行選取和設置，電池的運行工況即操作條件也與實驗保持一致。從圖中可以看出，模擬結果與實驗數據基本吻合，說明模型具有準確性和可靠性，可滿足計算要求。

圖5 模擬結果與實驗數據對比

2.2 S 型流道型線結構的影響

為研究S 型流道拉伸度對燃料電池性能的影響，采用控制變量法保證流道深度、寬度和彎曲度等參數不變，只改變拉伸度可得到六種不同的S 型流道設計方案，如表3 所示。全部采用正向組裝方式建立相應的單流道燃料電池模型，在相同工作電壓下比較六種方案的電池性能，并與相同條件下直型流道的結果進行對比。

表3 不同拉伸度的設計方案

圖6為不同拉伸度對燃料電池性能的對比，從圖中可以看出，氧氣的流動阻力和單電池的電流隨著拉伸度的增加均在減小，阻力的減小程度明顯大于電流的變化。同時計算得到了相同條件下直型流道的流動阻力為138 Pa，電流為0.594 A，而方案6 中S 型流道的流動阻力為143 Pa，稍大于直型流道，電流為0.591 A，與直型流道幾乎相同，這是因為拉伸度越大相同空間尺寸下可布置的S 型流道單元數目越少，對氧氣的流動擾動越小，那么S 型流道越來越接近直型流道，所以其性能也越接近直型流道。

圖6 不同拉伸度的對比

類似的，為研究S 型流道彎曲度對燃料電池性能的影響，表4 列出了不同彎曲度的設計方案，其對氧氣流動阻力和電池性能的對比如圖7 所示，從圖中可以看出，彎曲度越大對氧氣的流動擾動效果越強，傳質效果也越好，氧氣的流動阻力也越大，反之，彎曲度越小，流動阻力也越小，那么S 型流道也越接近直型流道。而單電池的電流并沒有隨著彎曲度的增加一直升高，這是由于擾動效果增強的同時加快了化學反應速率，使得反應生成的水增多，流道越彎曲排水越困難，更容易發生“水淹”，所以電流有所降低。方案10 中S 型流道的電流最大，最大電流為0.992 A，與直型流道相比，電流提升了的0.67倍。

表4 不同彎曲度的設計方案

圖7 不同彎曲度的對比

綜合比較彎曲度和拉伸度對電池性能的影響作用，可以發現彎曲度對氧氣的流動擾動作用較大，傳質效果較好，所以單電池的電流較大。而相對彎曲度的大小是由兩者共同決定的，單從相對彎曲度分析，流動阻力隨著相對彎曲度的增加

單調遞增，而單電池的電流隨著相對彎曲度的增加先升高后降低，存在一個最優值。特別的對于方案1 和方案8，方案1 的電流為0.624 A，方案8 的電流為0.615 A，雖然兩種方案中彎曲度和拉伸度均不相同，可相對彎曲度的大小卻是相同的，所以兩種方案的電流相差不大。

2.3 陰陽極S 型流道對應方式的影響

正反對應模型單電池的極化曲線如圖8 所示。從圖中可以看出，電池的電壓隨著電流密度的升高而逐漸降低，電池的功率密度隨著電流密度的升高先有所升高然后降低。在低電流密度下，兩者的極化曲線幾乎重合，說明兩者的性能相差不大，而在高電流密度下，兩者的性能有了明顯不同，尤其在電流密度大于1 A·cm-2時，反對應模型的電壓和功率密度均高于正對應模型，說明在高電流密度下，反對應模型的電池性能優于正對應模型，因此優先采用反對應的方式來組裝單電池。這是由于正對應模型陰陽極流道結構關于膜電極相互對稱，氫氣和氧氣與膜電極的公共接觸面積相同，同時陰陽極板與膜電極的公共接觸面積也完全相同，導致膜電極的總暴露面積較小，而反對應模型，氫氧流道相互交叉，氫氣和氧氣與膜電極公共接觸區域較小，但陰陽極板與膜電極的總暴露面積較大，提高了膜電極的利用率，從而使電池的性能有了提升。

圖8 兩種模型單電池的極化曲線

3 結論

本文基于S 型流道建立了多種不同的單流道燃料電池三維幾何模型，并采用數值模擬方法分別研究了S 型流道主要設計參數和S 型流道對應方式對燃料電池性能的影響，研究表明：

1）S 型流道的彎曲度和拉伸度都會影響燃料電池的性能，與拉伸度的影響作用相比，彎曲度的影響作用較大，兩者共同決定了相對彎曲度的大小，一般相對彎曲度越大，流動擾動越大，傳質效果越好，流動阻力也越大，同時也會造成排水比較困難電池性能下降，故存在最優的相對彎曲度值，使得流道在流動阻力和傳質效果上達到平衡。

2）通過對比正反對應模型的極化曲線可以發現，在高電流密度下，反對應模型膜電極的利用率有所提升，使得其單電池的性能優于正對應模型。因此，在單電池組裝時，優先選用反對應方式。