質子交換膜燃料電池雙極板流動結構仿生設計探究

摘要：針對質子交換膜燃料電池（PEMFC）內部燃料分布不均勻以及水管理能力不足導致的性能受限問題，基于海棠葉的葉脈結構，設計了一種新型仿生流場，采用數值方法深入研究了其性能。首先，通過觀察海棠葉脈的自然流動形態，設計仿生流場；然后，通過數值模擬對比分析仿生流場、傳統蛇形流場和平行流場的流場性能；最后，探究極限電流密度、峰值功率密度及內部氧氣和電流密度的分布特性。研究結果表明：與傳統的蛇形流場和平行流場相比，仿生流場的壓降更低，催化層氧氣濃度分布更均勻，供氣質量更高；相較于平行流場和蛇形流場，仿生流場的電池極限電流密度分別提升了35.42%、6.81%，峰值功率密度分別增加了29.71%、5.22%，且除水能力較平行流場有所提升。該研究為PEMFC流場設計的優化提供了新思路，對提升燃料電池的整體性能和效率具有重要的實際應用價值。

關鍵詞：質子交換膜燃料電池；仿生流場；海棠葉脈；雙極板

中圖分類號：TM911 文獻標志碼：A

DOI：10.7652/xjtuxb202504007 文章編號：0253-987X（2025）04-0070-11

Biomimetic Design Exploration of Flow Field Structure in Proton Exchange Membrane Fuel Cell Bipolar Plates

SHI Zihang1， YANG Yajing2， WEI Yanju1， ZHANG Chenyang1

（1. School of Energy and Power Engineering， Xi’an Jiaotong University， Xi’an 710049， China;

2. School of Astronautics Engineering， Xi’an Jiaotong University， Xi’an 710049， China）

Abstract：In response to performance limitations caused by uneven fuel distribution and inadequate water management capabilities in proton exchange membrane fuel cells （PEMFCs）， a novel bionic flow field is designed based on the vein structure of begonia leaves， and its performance is systematically investigated through numerical simulation. Firstly， the bionic flow field is developed by mimicking the natural flow patterns of begonia leaf veins. Then， numerical simulation techniques are employed to compare and analyze the flow field performance of the biomimetic flow field， traditional serpentine flow field， and parallel flow field. Finally， the distribution characteristics of limiting current density， peak power density， internal oxygen concentration， and current density are explored. The research results demonstrate that compared to traditional serpentine and parallel flow fields， the biomimetic flow field exhibits lower pressure drop， more uniform oxygen concentration distribution in the catalyst layer， and higher gas supply quality. Furthermore， the biomimetic flow field shows a 35.42% and 6.81% increase in limiting current density， a 29.71% and 5.22% increase in peak power density compared to the parallel and serpentine flow fields， respectively， with improved water retention capabilities compared to the parallel flow field. This research provides a new perspective for optimizing the design of flow fields in PEMFCs， offering significant practical value for enhancing the overall performance and efficiency of fuel cells.

Keywords：proton exchange membrane fuel cell; biomimetic flow field; begonia leaf veins; bipolar plate

近年來，能源消耗量的持續攀升加劇了環境污染^［1］。由于擁有較高能量密度和轉換效率、體積小、重量輕及快速充能等優點^［2-4］，燃料電池已逐漸成為極具潛力的可再生能源裝置。雙極板作為質子交換膜燃料電池（PEMFC）內部的關鍵組件，其結構不僅決定了流場的形狀，還負責將反應氣體高效傳輸至催化層，調控氣體的分布與流動狀態^［5-6］。

燃料電池傳統流場主要包括平行、蛇形、點狀以及交錯4種構型^［7］，其中平行和點狀流場壓降較小但容易發生堵塞，這在一定程度上阻礙了物質的傳遞效率^［8］；與平行流場相比，蛇形與交錯流場除水能力和電化學性能明顯提高，但均存在壓降顯著的缺點^［9-10］。目前，對傳統流場的改進思路主要集中在流場結構的設計和優化上^［11-13］，旨在優化反應氣體的分布狀態及流動路徑。

除了優化傳統流場結構外，還可以運用仿生流場的概念來設計流場。Kloess等^［14］根據自然界的葉片設計了一款流場，結果表明：與蛇形和交叉流動模式相比，仿生流動模式性能更好，壓降最低，物質分布非常均勻，功率密度峰值相較于傳統流場提升了30%。Chen等^［15］受葉脈結構啟發設計了一款仿生流場，相較于傳統平行流場，反應物的分布均勻性更好，消耗比更高，歐姆熱生成量更少，且制造成本可控。Kahraman等^［16］遵循穆雷定律設計了仿生流場，與標準蛇形設計相比，新設計的流場在0.4 V工作電壓下，功率提升了42.1%，電流分布均勻性、溫度優化分布及有效除水這3方面性能也得到了提升。

上述研究表明：仿葉脈流場可有效提升燃料電池電流密度、功率輸出及穩定性，改進水傳輸與排除，減少水淹現象。然而，現有設計多聚焦于形態模仿，缺乏對流動特性與性能提升關聯的深入分析。本文基于海棠葉脈設計新型仿生流場，通過仿真研究其性能，并在相同條件下與平行流場和蛇形流場對比，分析了仿生流場的優勢及流動特性對電化學性能的影響。

1 模型建立

1.1 幾何模型

在漫長的自然選擇和進化中，植物葉脈形成了高效且均勻的物質傳輸網絡^［17］，葉片中的導管類似毛細管，通過毛細現象產生毛細壓力，促進水分的被動吸收，由于導管內水分的內聚力小于與管壁的附著力，因此能夠推動水分沿葉脈上升，確保葉片得到充足水分；而葉脈中的水分流動呈不可壓縮穩定層流狀態，有助于水分高效且均勻地分布到葉片各部分。海棠，因其出色的抗寒、抗旱、抗澇、抗鹽堿和抗病蟲能力，成為重點研究對象，其高效營養運輸能力可為流場結構設計提供重要啟示。本文采用熒光溶液示蹤法，通過電子顯微鏡觀察海棠葉脈網格中水分分布隨時間演變的動態過程，如圖1所示，其中白色為熒光溶液流經區域，黑色為其他區域。

基于觀察到的海棠葉脈流動特點，本文在極板設計過程中，創新性地采用將流入流道與單元流場相結合的設計理念。具體來說，流入流道的設計靈感來源于海棠葉脈結構，而單元流場則借鑒了離散化微脈網格區域的設計思想。

圖2展示了燃料電池仿生流場的極板模型、極板剖面和流體域。圖2（a）為仿葉脈結構的燃料電池雙極板模型。圖2（b）為仿葉脈結構的流場，包括流場入口、流場出口和兩組流道，流場入口設置在葉柄處，出口設置在葉尖處，兩組流道對稱設置在主脈兩側，流道包括主流入流道、多個二級流入流道、多個二級流出流道、主流出流道和若干個單元流場通道。圖2（c）為兩個單元間的流道構造。圖2（d）為仿生流場流體域與擴散層接觸的一面，圖中紅色部分是流道下方流場與擴散層被隔開的區域，藍色部分是流場與氣體擴散層相接觸的區域。在流場設計中，不同的單元流場之間采用并聯連接方式，每一個單元的流場入口都通過流入流道與流場入口相連，每一個單元的流場出口都通過流出流道與流場出口相連。相較于各單元串聯，并聯結構的反應物濃度均勻性更好；同時，采用逆流模式操作電池，使氧氣和氫氣的入口位于電池對角線的兩端，既利用了逆流布置的自加濕功能，也進一步提升了反應物濃度分布的均勻性^［18］。

與PEMFC常用的平行流場和蛇形流場構型相比，上述仿生流場構型內部壓力適中，具有多條流路設計，在一定程度上有效解決了平行流場中由于壓差較小而導致的堵塞問題，提高了物質傳遞效率，同時也緩解了蛇形流場中單一流路導致的壓降過大問題。此外，仿生流場獨特的單元流場并聯設計能夠在理想情況下實現反應物的單元化，從而促進反應物濃度分布更均勻。

通過對流道進行流動仿真，不斷優化，最終對仿生流場的結構特征做出如下設置：①由于入口形狀為矩形，因此流道形狀可近似為一個細長四棱錐；②流場入口尺寸和流場面積之間不存在直接關聯；③為了簡化流場構造，使流場兩側流動對稱，將二級流入流道和主流入流道分支角度設置為45°；④二級流入流道末端寬度為0.2mm，為了盡可能保證每個單元流場中的流量相同，將流道兩側壁面與水平線的夾角設置為1.6°；⑤將氣體擴散層與流體流入通道隔開，為了避免被隔開的部分出現反應氣體饑渴區，留下一條寬為0.2mm的縫隙使氣體與擴散層接觸，縫隙與兩側單元流場的距離保持在0.1mm以上；⑥單元流場輪廓近似于邊長為7mm的菱形，各個流場輪廓形狀根據流入流出流道、流場邊緣和單元流場出入口微調；⑦為了減小流動阻力，對靠近主流出流道的單元流場輪廓適當進行圓角處理。

PEMFC仿真模型的幾何參數如表1所示，參考文獻［19-20］選擇合適的擴散層厚度、質子交換膜厚度及流場高度等參數，且在建模時保持不變。

為了更好地評估葉脈仿生流場燃料電池性能，構建實際工程應用中廣泛使用的平行流場和雙通道蛇形流場燃料電池模型用于對比。圖3（a）和3（b）分別展示了平行流場和雙通道蛇形流場的流體域模型。為使這兩種流場燃料電池中反應物與氣體擴散層的接觸面積與仿生流場燃料電池一致，平行流場通道寬度統一設定為1.00mm，通道間距為1.00mm；蛇形流場通道寬度設為0.80mm，通道間距為0.75mm。 除此之外，其余模型幾何參數均與葉脈仿生流場燃料電池模型保持一致，以確保對比的公正性和準確性。

1.2 數學模型

PEMFC在運行過程中涉及復雜的多物理場耦合過程，為了精準地模擬PEMFC內部流動、傳熱、傳質和電化學反應過程，本節構建三維、穩態的數學模型，并對其進行了嚴格驗證。模型涵蓋如下所示的控制方程。

質量守恒方程^［21-23］可寫為

（ερ）t+ Δ ·（ερu）=S m （1）

式中：ε為孔隙率；ρ為密度；u為速度；S m 為質量源項，為模擬中引入或移除質量的速率^［24］，在氣體擴散層和流道中S m =0。陽極催化層的質量源項S ma 和陰極催化層的質量源項S mc 分別表示如下

S ma =－M（ H 2）2Fj a （2）

S mc =M（ H2O ）2F－M（ O 2）4Fj c （3）

式中：下標 a、c 分別表示陽極和陰極；M（ H 2）、M（ H2O ）、M（ O2 ）分別為氫氣、水蒸氣、氧氣的摩爾質量；j為體積電流密度；F為法拉第常數。

動量守恒方程^［21-23］可寫為

（ερu）t+ Δ ·（ερuu）=－ε Δ p+ Δ ·（εμ Δ ·u）+S u （4）

式中：p為流體壓力；μ為流體動力黏度；S u 為動量源項。氣體流場中S u =0。催化劑層（ CL ）和氣體擴散層（ GDL ）的孔隙中，S u ^［24］可由下式表示

S u =－μKεu（5）

式中：K為GDL和CL的滲透率。

能量守恒方程可寫為

（ερcpT）t+ Δ ·（ερcpuT）= Δ ·（λ Δ T）+S Q （6）

式中：cp為比定壓熱容；T為溫度；λ為熱導率； S Q 為能量源項。S Q 表示PEMFC運行時耗散的能量^［24］，涵蓋了歐姆損失、化學反應損失、相變損失和過電位損失等，可由下式表示

S Q =I2R ohm +ηI change +h phase +h reaction （7）

式中：I為電流密度；R ohm 為電池內阻；η為過電位；I change 為交換電流密度；h phase 為水的相變熱；h reaction 為化學反應焓。

組分守恒方程可寫為

（ερck）t+ Δ ·（ερuck）= Δ ·（ρD ^eff k Δ ck）+Sk（8）

式中：ck為第k種組分的濃度；D ^eff k為自由流質量擴散系數的有效值；Sk為組分源項。

電荷守恒方程可寫為

Δ ·（σ s" Δ φ s ）+S s =0（9）

Δ ·（σ m" Δ φ m ）+S m =0（10）

式中：σ s 、σ m 分別為固相、膜相電導率；φ s 、φ m 分別為固相、膜相電勢；S s 、S m 分別為固相、膜相電流源項。

交換電流密度是對陰、陽極化學反應速率的直觀表達，由于電路中的電子傳導效率以及膜電極組件中的離子傳導效率共同決定了反應速率，因此可以由活化過電勢得到交換電流密度。Butler-Volmer方程為電極上的氧化還原反應過程以及上述化學反應速率與電位之間的關系提供了有效數學描述，常被用于求解交換電流密度^［25］，可由下式表示

R a =（ξ a j ^ref "a ）c（ H 2）c ^ref （ H 2）^{γ a} （ e ^{α a Fη a /RT}－ e ^{－α c Fη a /RT}）（11）

R c =（ξ c j ^ref "c ）c（ O 2）c ^ref （ O 2）^{γ c} （－ e ^{α a Fη c /RT}+ e ^{－α c Fη c /RT}）（12）

式中：ξ為活化比表面積；j ^ref 為參考交換電流密度；c ^ref （ H 2）、c ^ref （ O 2）、c（ H 2）、c（ O 2）分別為 H 2和 O 2的參考濃度及濃度；γ為濃度指數；α為電荷傳輸系數；η為活化過電位；R為理想氣體常數。

1.3 邊界條件與模型假設

在相同邊界條件下，對3種流場的PEMFC進行仿真。仿真采用質量流量入口和壓力出口設定，入口流動為層流，出口為壓力出口，以確保化學反應穩定。離子流僅限于膜電極組件內部，外部不穿越邊界。為簡化計算，膜相電勢設為0，電子僅在陰陽極外部電流中流動；陽極電勢設為0，陰極電勢為電池電壓；其他壁面采用無滑移邊界條件。燃料電池模型的主要參數如表2所示。

為了簡化計算，同時確保化學反應穩定，假設燃料電池不受外力影響，工作環境保持穩態；忽略材料性質的影響和不可逆損失；氫氣和氧氣均設為不可壓縮的理想氣體且在PEMFC流道內層流流動；模擬過程不考慮重力的影響。

1.4 模型驗證

為確保所構建燃料電池模型的科學性和可靠性，基于實驗數據對模型進行了詳盡且嚴格的驗證工作。以湯夢琪等^［26］的實驗數據為基礎，構建了一個單流道燃料電池模型，模型結構的示意圖如圖4所示，具體的幾何參數如表3所示。

在設定模型邊界條件時，嚴格遵循1.2和1.3節中確定的源項配置、模型設定、簡化策略以及網格劃分方案，以確保模型的準確性和一致性。驗證模型所選取的操作參數、主要物理性質參數及仿真環境均與既有實驗研究的數據完全一致。

圖5展示了經過網格獨立性驗證后的數值結果與文獻［26］實驗結果的對比。可以看出，在低電流密度（小于1.2A·cm^-2）下，模擬與實驗結果一致性較好；但隨著電流密度增大，偏差增大，特別是在高電流密度時，由于電滲作用，水的生成增加且難以排除，導致“水淹”現象發生。這是由于當前模型未考慮液態水在氣體擴散層中的占比，造成水的狀態未能被準確模擬，可能導致低壓條件下的性能被高估。總體來看，模擬與實驗結果的誤差不超過3.78%，擬合度較好，驗證了仿真模型的可靠性。

為了驗證本文仿生流場模型的網格獨立性，設置4種網格劃分方案。根據模型中質子交換膜（PEM）和CL部件的形狀規則，采用六面體網格對其進行劃分，以保證計算精度和效率。仿生流場的形狀呈現高度復雜性，故流體域部分采用自由四面體網格劃分以適應復雜的幾何形態。GDL需要和氣體流場共享拓撲網格，因此采取共節點網格策略使用自由四面體網格劃分。對于平行流場，由于形狀規則，不同部件均采用六面體網格劃分。蛇形流場中圓角部分采用自由四面體網格劃分，以適應其幾何變化。流體域直通道和其他部件均采用正六面體網格掃掠，以提高網格質量和計算精度。不同流場結構的PEMFC網格劃分方案見表4。

采用上述4種網格劃分方案對不同流場的燃料電池進行獨立性檢驗，在工作電壓0.5V條件下比較電流密度的變化趨勢，如圖6所示。結果顯示：仿生流場網格數達274萬時，與583萬網格的相對誤差為0.15%；平行流場網格數為187萬時，與385萬網格的誤差為0.99%；蛇形流場網格數為234萬時，與477萬網格的誤差僅為0.03%。為兼顧精度和計算效率，后續模型計算均采用方案2。

2 結果討論

2.1 氧氣流動與氧氣濃度的分布

在燃料電池運行中，陰極GDL與CL交界面的氧氣濃度分布是影響電化學反應和水生成的關鍵因素。為確保反應高效進行，必須將交界面氧氣濃度維持在適當且均勻的水平。

圖7給出了工作電壓為0.6V時，不同流場陰極GDL與CL交界面處的氧氣濃度分布。由圖可見，仿生流場、平行流場和蛇形流場的氧氣濃度峰值分別為4.40、4.03、4.98mol·m^－3，最小值分別為0.46、1.58×10^－3、1.19×10^－2mol·m^－3，仿生流場的氧氣濃度峰值最低，且是唯一未出現氧饑渴區的。

進一步分析仿生流場中的氧氣濃度分布可見，每個單元流場的入口氧氣濃度較高，出口較低，但遠離入口區域濃度仍較高，整體分布較為均勻。流道與反應區隔離使得部分區域氧氣濃度較低，氧氣分布呈現先骨架化流動，再離散區域擴散的特征。相比之下，平行流場在入口及連接流道處氧氣濃度最高，但周圍和中部供氧不足；蛇形流場氧氣濃度差最大，濃度從入口向出口逐漸遞減，梯度明顯。

圖8展示了電壓為0.6V時，氧氣在陰極流道內的流動軌跡和速度分布。由圖可見，在仿生流場中，主流道和各二級流道的流速較高，表明氧氣自進入流場后就迅速充滿整個“骨架”，抵達每個單元流場入口，隨后順利進入各個單元流場。不同單元流場之間的流速相近，為氧氣的均勻輸送奠定了基礎。

由圖8還可以看出，在平行流場中，由于中部豎直流道流速較慢，氧氣通量降低，催化層與擴散層交界處易形成大面積供氧不足區域；而在雙通道蛇形流場中，在相同的質量流量下，由于全流程僅有兩個流道，流場流通截面積較小，因此流道中氧氣的流速非常高，展現出良好的氧氣運輸能力。

為了評估燃料電池反應物及電流密度在陰極側擴散層與催化層交界面中分布的均勻性，引入面加權均勻度指數^［27］，可寫為

ω a =1－∑ni=1（Xi－ a ）Ai2 a ∑ni=1Ai（13）

式中：Xi、 a 分別為區域Ai內物質X的濃度/電流密度及平均濃度/平均電流密度；n為區域數。

圖9展示了工作電壓從1.0V降至0.3V時，不同流場陰極側GDL與CL交界面氧氣濃度分布均勻度指數的變化趨勢。結果顯示：當工作電壓高于0.85V時，3種流場的氧氣濃度均勻度相近且較高；當電壓降至0.85V以下，平行流場的均勻度指數開始下降，速度逐漸減緩，而仿生流場和蛇形流場的均勻度指數呈單調下降趨勢。值得注意的是，在所有工作電壓下，仿生流場的氧氣濃度均勻度指數均優于平行流場和蛇形流場。

在燃料電池的催化層，氧氣濃度受到流場氧氣分配和氧氣消耗速率的影響。高工作電壓下，電流較小，氧氣消耗速率低，各流場的氧氣較為充足，均勻度較高。隨著電壓降低，電流增大，氧氣消耗速率上升。當電壓降至0.4V時，仿生流場和平行流場的電流密度增長放緩，氧氣消耗速率減慢，均勻度指數下降速度減緩。相反，蛇形流場在低電壓下仍保持較快的電流增長，導致氧氣濃度均勻度指數持續快速下降。

綜上所述，3種流場中，仿生流場的氧氣濃度分布最為均勻，供氣質量最佳，這一結果對于優化燃料電池性能、提高氧氣利用效率具有重要意義。

2.2 氣體擴散層與催化層交界面處水的濃度分布

燃料電池的排水性能對反應物傳質、化學反應和離子傳導有重要影響。為了評估新型流場的排水能力，需對比分析陰極GDL與CL交界面處水的濃度。

圖10展示了0.6V工作電壓下陰極GDL與CL交界面處水的濃度分布，結果顯示水的濃度與氧氣流場的流動狀態密切相關。在仿生流場中，由于氧氣流入通道與氣體擴散層的隔離，流道下方的擴散層和催化層缺乏對流現象，導致水的濃度較高。此外，由于水的排出方向與氣體流向一致，單元流場出口及流出流道下方催化層的水的濃度也較高。

在平行流場中，水的濃度較高的區域主要集中在邊緣和中部。邊緣區域因缺乏流道，水難以排出；中部由于氧氣流速較低，缺乏對流，導致水排除不充分。相比之下，由于流速較高，蛇形流場強烈的對流現象能有效排水，水的濃度從入口到出口逐漸增高，除邊緣少量高濃度區域外，其余部分的水均得到了有效去除。經過計算，3種流場中陰極側GDL與CL交界面處水的濃度均值分別為13.206、14.547、13.098mol·m^－3，表明蛇形流場的除水能力最佳，仿生流場的除水能力僅次于蛇形流場。

2.3 氣體擴散層與催化層交界面的電流密度分布

電流密度分布的均勻性對燃料電池的穩定運行至關重要。過高或過低的電流密度會導致局部反應過快或過慢，不僅降低整體性能，還會縮短質子交換膜的壽命，影響燃料電池的長期使用。

圖11展示了0.6V工作電壓下，不同流場陰極側GDL與CL交界面的電流密度分布。由圖可見，仿生流場和蛇形流場的電流密度分布較為均勻，幾乎沒有過低區域，而平行流場中部由于氧氣流速較低，導致電流密度過低，反應較慢。進一步分析顯示，仿生流場的電流密度峰值與最低值之差為0.2320A·cm^－2，平行流場為0.8364A·cm^－2，蛇形流場為0.7310A·cm^－2，表明仿生流場沒有局部反應過快或過慢的區域。與平行流場和蛇形流場相比，仿生流場的峰值電流密度分別降低了31.35% 和46.31%，這進一步證明其電流密度分布優越，有助于提升燃料電池性能。

工作電壓從1.0V降低至0.3V的過程中，不同流場陰極側氣體GDL與CL交界面處電流密度分布均勻度指數隨電壓的變化趨勢如圖12所示。由圖可見，隨著工作電壓降低，3種流場結構的電流密度均呈現上升趨勢，化學反應速率加快，導致氧氣消耗速率顯著增強。這一現象使得催化層中不同區域的氧氣濃度梯度逐漸增大，進而影響了電流密度分布，使其均勻性逐漸降低。

在所有模擬工況下，仿生流場均表現出較好的電流密度分布均勻性，這是由于仿生流場的設計能夠使氧氣分布更為均勻，有效降低了濃差極化損失，為電流密度的均勻分布奠定了堅實基礎。這一特性使得仿生流場在燃料電池運行過程中的性能更為穩定，從而具有潛在的應用優勢。

2.4 流場壓降

圖13展示了不同流場陰極側氣體通道的壓力分布。由圖可見，仿生流場的壓降為444.36Pa，平行流場為864.80Pa，蛇形流場為5 021.40Pa，顯示出仿生流場在壓降控制上的明顯優勢。這一優勢主要歸結于仿生流場獨特的結構特點，由于其內部具有眾多分叉，使得氣體在流場內的流速相對較低，且流動時的拐角基本小于45°，因此有效減少了局部壓力損失。蛇形流場僅包含兩條獨立的流通路徑，氣體流速較高，且存在多個90°拐角，因而導致大量的局部和沿程壓力損失。平行流場雖然流速不高，但其拐角同樣為90°，流道拐角處產生較多旋渦，因此也造成了較大的局部壓力損失。綜合比較后可知：平行流場的壓降最小，仿生流場次之，蛇形流場壓降最大，這一結果驗證了仿生流場在燃料電池設計中的優越性，特別是在承受更高入口流速方面展現出巨大潛力，為燃料電池流場結構的設計和優化提供了有益參考。

2.5 電化學性能

極化曲線與功率密度曲線作為燃料電池電化學性能的重要表征方式，能夠體現出不同流場設計方案對燃料電池性能的影響。圖14顯示了不同流場的極化與功率密度曲線隨電流密度的變化情況。

由圖14可見，在低電流密度區域，3種流場的功率密度與極化曲線差異不大，這一現象主要歸因于數值模擬過程中，與膜的特性緊密相關的活化損失占主導地位，而在高電壓工況下，膜的特性變化通常被忽略。隨著電壓降低，濃差極化效應逐漸占據主導地位，3種流場間的性能差異逐漸擴大，特別是在低工作電壓條件下，平行流場的輸出性能相較于蛇形和仿生流場呈現顯著下降趨勢；而仿生流場則展現出最為優越的輸出性能，相較于平行流場和蛇形流場，仿生流場的極限電流密度分別提高了35.42%、6.81%，峰值功率密度分別提高了29.71%、5.22%。

3 結 論

為深入探究不同流場結構對質子交換膜燃料電池性能的影響，本文基于海棠葉脈仿生結構設計了一種新型仿生流場，采用數值模擬軟件構建了仿生流場、平行流場及蛇形流場3種燃料電池仿真模型，并進行對比分析，得到如下結論。

（1）與傳統的蛇形流場和平行流場相比，葉脈仿生流場壓降更低，催化層氧氣濃度分布更均勻，供氣質量更高，有利于降低泵送功率和提升燃料電池電流密度的分布均勻性。

（2）相較于傳統的蛇形流場和平行流場，仿生流場有著較強的除水能力，雖略低于蛇形流場，但強于平行流場，在低電壓工況下運行時，較強的除水能力能有效減少水淹現象。

（3）仿生流場的電化學輸出性能最為優越，相較于平行流場和蛇形流場，仿生流場的極限電流密度分別提高了35.42%、6.81%，峰值功率密度分別提高了29.71%、5.22%。

參考文獻：

［1］FICHMAN B T. Annual energy review 2011 ［EB/OL］. （2012-09-01）［2024-05-12］. https：//www.osti.gov/biblio/1212312/.

［2］葉好， 韓丁波， 張劍波， 等. 車用燃料電池質子交換膜研究進展 ［J］. 時代汽車， 2024（6）： 126-128.

YE Hao， HAN Dingbo， ZHANG Jianbo， et al. Advances in proton exchange membrane research for automotive fuel cells ［J］. Auto Time， 2024（6）： 126-128.

［3］WANG Yun， RUIZ DIAZ D F， CHEN K S， et al. Materials， technological status， and fundamentals of PEM fuel cells： a review ［J］. Materials Today， 2020， 32： 178-203.

［4］XIAO Fei， WANG Yucheng， WU Zhipeng， et al. Recent advances in electrocatalysts for proton exchange membrane fuel cells and alkaline membrane fuel cells ［J］. Advanced Materials， 2021， 33（50）： 2006292.

［5］胡翀. 質子交換膜燃料電池多物理場耦合仿真研究與流場結構優化設計 ［D］. 鎮江： 江蘇科技大學， 2023.

［6］HEIDARY H， KERMANI M J， DABIR B. Influences of bipolar plate channel blockages on PEM fuel cell performances ［J］. Energy Conversion and Management， 2016， 124： 51-60.

［7］趙永志， 蒙波， 陳霖新， 等. 氫能源的利用現狀分析 ［J］. 化工進展， 2015， 34（9）： 3248-3255.

ZHAO Yongzhi， MENG Bo， CHEN Linxin， et al. Utilization status of hydrogen energy ［J］. Chemical Industry and Engineering Progress， 2015， 34（9）： 3248-3255.

［8］趙強， 郭航， 葉芳， 等. 質子交換膜燃料電池流場板研究進展 ［J］. 化工學報， 2020， 71（5）： 1943-1963.

ZHAO Qiang， GUO Hang， YE Fang， et al. State of the art of flow field plates of proton exchange membrane fuel cells ［J］. CIESC Journal， 2020， 71（5）： 1943-1963.

［9］LI W Z， YANG W W， ZHANG W Y， et al. Three-dimensional modeling of a PEMFC with serpentine flow field incorporating the impacts of electrode inhomogeneous compression deformation ［J］. International Journal of Hydrogen Energy， 2019， 44（39）： 22194-22209.

［10］MANSO A P， MARZO F F， BARRANCO J， et al. Influence of geometric parameters of the flow fields on the performance of a PEM fuel cell： a review ［J］. International Journal of Hydrogen Energy， 2012， 37（20）： 15256-15287.

［11］JIAO Kui， BACHMAN J， ZHOU Yibo， et al. Effect of induced cross flow on flow pattern and performance of proton exchange membrane fuel cell ［J］. Applied Energy， 2014， 115： 75-82.

［12］LIU Haichao， YANG Weimin， TAN Jing， et al. Numerical analysis of parallel flow fields improved by micro-distributor in proton exchange membrane fuel cells ［J］. Energy Conversion and Management， 2018， 176： 99-109.

［13］ABDULLA S， PATNAIKUNI V S. Detailed analysis of polymer electrolyte membrane fuel cell with enhanced cross-flow split serpentine flow field design ［J］. International Journal of Energy Research， 2019， 43（7）： 2806-2820.

［14］KLOESS J P， WANG Xia， LIU J， et al. Investigation of bio-inspired flow channel designs for bipolar plates in proton exchange membrane fuel cells ［J］. Journal of Power Sources， 2009， 188（1）： 132-140.

［15］CHEN Tao， XIAO Yong， CHEN Tiezhu. The impact on PEMFC of bionic flow field with a different branch ［J］. Energy Procedia， 2012， 28： 134-139.

［16］KAHRAMAN H， HA石自航1，楊亞晶2，魏衍舉1，張晨楊1IMOGLU C， EVIK C， et al. A different flow field design approach for performance improvement of a PEMFC ［J］. Acta Physica Polonica： A， 2017， 131（3）： 484-486.

［17］COCHARD H， NARDINI A， COLL L. Hydraulic architecture of leaf blades： where is the main resistance？ ［J］. Plant， Cell amp; Environment， 2004， 27（10）： 1257-1267.

［18］KHAZAEE I， SABADBAFAN H. Effect of humidity content and direction of the flow of reactant gases on water management in the 4-serpentine and 1-serpentine flow channel in a PEM （proton exchange membrane） fuel cell ［J］. Energy， 2016， 101： 252-265.

［19］WANG Yun， CHEN K S， MISHLER J， et al. A review of polymer electrolyte membrane fuel cells： technology， applications， and needs on fundamental research ［J］. Applied Energy， 2011， 88（4）： 981-1007.

［20］陶慧芳. 基于多物理場耦合的氫燃料電池極板流道分析及優化 ［D］. 濟南： 山東大學， 2022.

［21］GONG Fan， YANG Xiaolong， ZHANG Xun， et al. The study of Tesla valve flow field on the net power of proton exchange membrane fuel cell ［J］. Applied Energy， 2023， 329： 120276.

［22］YUAN Wei， TANG Yong， PAN Minqiang， et al. Model prediction of effects of operating parameters on proton exchange membrane fuel cell performance ［J］. Renewable Energy， 2010， 35（3）： 656-666.

［23］FAN Linhao， ZHANG Guobin， JIAO Kui. Characteristics of PEMFC operating at high current density with low external humidification ［J］. Energy Conversion and Management， 2017， 150： 763-774.

［24］弗朗諾·巴爾伯. PEM燃料電池： 理論與實踐 ［M］. 李東紅， 連曉峰， 譯. 北京： 機械工業出版社， 2016.

［25］LIM B H， MAJLAN E H， DAUD W R W， et al. Numerical analysis of modified parallel flow field designs for fuel cells ［J］. International Journal of Hydrogen Energy， 2017， 42（14）： 9210-9218.

［26］湯夢琪. 質子交換膜燃料電池雙極板熱導率的仿真研究 ［D］. 北京： 清華大學， 2019.

［27］TAO Hongge， CHEN Huanxin， XIE Junlong， et al. An alternative approach to quantifying fluid flow uniformity based on area-weighted average velocity and mass-weighted average velocity ［J］. Energy and Buildings， 2012， 45： 116-123.

（編輯 李慧敏）