中圖分類號(hào)：TM911.4 文獻(xiàn)標(biāo)志碼：A

目前，傳統(tǒng)化學(xué)能源儲(chǔ)量日益減少，其燃燒帶來的環(huán)境問題也愈發(fā)嚴(yán)重，推動(dòng)全球向清潔能源轉(zhuǎn)型。在此背景下，低噪聲、零污染排放且具備快速響應(yīng)特性的質(zhì)子交換膜燃料電池（protonexchangemembranefuelcell，PEMFC）展現(xiàn)出獨(dú)特的技術(shù)優(yōu)勢[1-]。該裝置通過氫氧電化學(xué)反應(yīng)直接輸出電能，工作過程中僅生成水與熱能，具有顯著的環(huán)境友好優(yōu)勢[3]。盡管PEMFC的應(yīng)用潛力廣泛，但其商業(yè)化進(jìn)程仍面臨耐久性不足和制造成本偏高的問題，而熱管理就是解決其壽命問題的重要手段[4]

PEMFC在電化學(xué)反應(yīng)中產(chǎn)生的大量廢熱對運(yùn)行穩(wěn)定性構(gòu)成潛在威脅，熱量的持續(xù)積聚會(huì)引發(fā)膜電極組件脫水，最終加速電解質(zhì)膜的結(jié)構(gòu)劣化[5]這種熱失控效應(yīng)不僅會(huì)造成輸出功率減小，還會(huì)縮短電池的使用壽命。由此可見，高效的熱管理對維持PEMFC的性能與穩(wěn)定性具有關(guān)鍵意義。

現(xiàn)階段，PEMFC多在雙極板內(nèi)設(shè)置冷卻流道，或配置專門的冷卻板以實(shí)現(xiàn)散熱，對于高功率密度的電堆裝置，液體冷卻因具備較強(qiáng)的冷卻能力，成為優(yōu)先選擇[6]。其冷卻機(jī)制是冷卻液在冷卻流道內(nèi)循環(huán)流動(dòng)，帶走PEMFC在運(yùn)行過程中產(chǎn)生的熱量。目前，PEMFC的冷卻流道設(shè)計(jì)有平行流道、蛇形流道、點(diǎn)陣形流道、交指形流道4種基本形式。不同的流道配置會(huì)導(dǎo)致不同的流動(dòng)方向，進(jìn)而產(chǎn)生不同的冷卻效果，因此流道的結(jié)構(gòu)和布局是決定冷卻效率的關(guān)鍵因素。

在PEMFC的冷卻流道優(yōu)化研究中，眾多研究者通過改變流道的幾何形狀實(shí)現(xiàn)了熱管理性能的顯著提升。研究表明，三維混沌流道[、特斯拉閥式[8]鋸齒形[9]及波浪形[10]等流道結(jié)構(gòu)，均通過增強(qiáng)流體擾動(dòng)、形成二次渦旋，實(shí)現(xiàn)了傳熱效率提升，其中波浪形流道還可通過調(diào)節(jié)波峰數(shù)量實(shí)現(xiàn)熱傳導(dǎo)與流阻的平衡[1I-12]。此外，擴(kuò)散器型、噴嘴狀[13]及梯形流道[14]則在加速廢熱、液態(tài)水排出方面具有優(yōu)勢。

PEMFC的熱管理效果在很大程度上也取決于冷卻流道的排列方式。相關(guān)研究表明，多螺旋、蜂窩狀等非傳統(tǒng)構(gòu)型的冷卻流道能有效改善溫度分布均勻性，但會(huì)導(dǎo)致壓降的上升；而多通道布局流道可降低局部流動(dòng)阻力并提升溫度均勻性[15-6] C部分新型流道設(shè)計(jì)則實(shí)現(xiàn)了低壓降條件下的高效散熱性能，為流道布局提供了兼顧效率與能耗的多元化解決方案[17-18] ○

研究發(fā)現(xiàn)，在冷卻流道內(nèi)設(shè)置阻塊、凹槽等擾流結(jié)構(gòu)，均能顯著提升PEMFC的傳熱性能。其中，阻塊結(jié)構(gòu)通過增強(qiáng)流體擾動(dòng)提升傳熱效率，但可能導(dǎo)致壓降增加[19]；而凹槽結(jié)構(gòu)在強(qiáng)化傳熱的同時(shí)，可降低壓力損失，為兼顧高效散熱與低壓降的流道設(shè)計(jì)提供了新方向[20]。還有部分研究者運(yùn)用仿生設(shè)計(jì)優(yōu)化PEMFC的冷卻流道結(jié)構(gòu)，有效提升了熱管理性能。仿生葉脈、魚骨及樹狀分形等流道設(shè)計(jì)[2I-23]，不僅優(yōu)化了溫度均勻性，降低了電池溫度與系統(tǒng)壓降，還改善了冷卻液的流動(dòng)分布，實(shí)現(xiàn)了低能耗下的高效散熱。

綜上所述，良好的冷卻流道設(shè)計(jì)能夠?qū)崿F(xiàn)更高效的熱傳遞。然而，部分復(fù)雜流道在實(shí)際加工中面臨工藝難度大、制造成本高的問題，而傳統(tǒng)平行流道仍存在散熱效率低、溫度分布不均的固有缺陷。針對平行流道的熱管理問題，本研究聚焦于多通道分流結(jié)構(gòu)的創(chuàng)新設(shè)計(jì)，提出3種新型冷卻流道，通過數(shù)值模擬探究其在不同工況下的傳熱性能與流動(dòng)特性，以期為解決高功率密度PEMFC的熱管理難題提供兼具工程可行性與高效散熱能力的新方案。

1 模型與參數(shù)

1.1 幾何參數(shù)

質(zhì)子交換膜燃料電池的基本結(jié)構(gòu)通常由陰陽極雙極板（bipolarplate，BP）、氣體擴(kuò)散層（gasdif-fusionlayer，GDL）、催化層（catalystlayer，CL）及質(zhì)子交換膜（protonexchangemembrane，PEM）組成[24]。其熱管理機(jī)制主要依賴于雙極板內(nèi)設(shè)計(jì)的冷卻流道或獨(dú)立冷卻板，通過冷卻液在流道內(nèi)的循環(huán)流通，將電堆運(yùn)行中產(chǎn)生的熱量及時(shí)帶走。這種基于流體動(dòng)力學(xué)的散熱機(jī)制已成為保障PEMFC穩(wěn)定運(yùn)行的關(guān)鍵熱管理措施。圖1展示了設(shè)置冷卻板的燃料電池典型結(jié)構(gòu)。

本研究采用如圖2所示的熱通量模型開展熱力學(xué)分析，該模型包括冷卻板與冷卻流道。在冷卻板底面施加熱通量邊界條件，用以模擬熱量經(jīng)PEMFC雙極板傳遞至冷卻板的過程。該建模方法通過關(guān)鍵參數(shù)的等效替代，實(shí)現(xiàn)了復(fù)雜熱傳導(dǎo)過程的可控模擬。運(yùn)用數(shù)值模擬方法探究平行流道和3種新型流道對PEMFC傳熱性能的影響，流道結(jié)構(gòu)如圖3所示。

1.2 物性參數(shù)與邊界條件

本研究中冷卻板為石墨，冷卻液為液態(tài)水，同時(shí)忽略重力對系統(tǒng)的影響。冷卻液的初始溫度設(shè)定為313K，其入口設(shè)定為速度入口，出口為壓力出口，并且速度方向垂直于進(jìn)出口平面。在計(jì)算過程中，采用SIMPLE算法求解連續(xù)性方程，運(yùn)用Navier-Stokes方程進(jìn)行流動(dòng)迭代計(jì)算，從而獲取數(shù)值模擬結(jié)果。當(dāng)?shù)鷼埐顪p小到 1×10^-8 ，此時(shí)溫度、壓降等各種參數(shù)趨于穩(wěn)定，可認(rèn)為達(dá)到收斂狀態(tài)。具體的邊界條件見表1。

2 模型假設(shè)

本研究針對冷卻流道內(nèi)冷卻液的流動(dòng)傳熱特性，在數(shù)值模擬中作出以下理想化假設(shè)：

1）冷卻液被視為不可壓縮牛頓流體，流動(dòng)狀態(tài)為穩(wěn)態(tài)層流；

2）冷卻液的入口流速和溫度設(shè)定為恒定值；

3）物理參數(shù)假定為不隨溫度變化的常數(shù)；

4）冷卻板表面施加均勻熱通量。

3 數(shù)學(xué)模型方程

3.1 質(zhì)量守恒方程

對于不可壓縮流體，方程可簡化為 ，即

式中： ρ 為流體密度； 是速度矢量； u，v，w 分別是 x，y，z 方向的速度分量

3.2 動(dòng)量守恒方程

式中： ρ 為流體密度； 為速度矢量； p 為流體壓力; μ 為動(dòng)態(tài)黏度; ，F(xiàn)_z ）為體積力矢量。

3.3 能量守恒方程

式中： c_p 是流體的定壓比熱容； T 是溫度; k 是熱導(dǎo)率； 是單位體積內(nèi)的熱量生成率。

3.4對流換熱方程

3.5溫度均勻性指數(shù)（indexofuniformtemperature ，IUT）

式中： I_UT 為溫度均勻性指數(shù)； 表示第 ∣m∣ 個(gè)區(qū)域的溫度； 表示整個(gè)系統(tǒng)的面積加權(quán)平均溫度； A_m 為第 m 個(gè)區(qū)域的面積； n 為區(qū)域的總數(shù)。當(dāng)溫度分布完全均勻時(shí)， I_UT= 1 ；溫度分布越不均勻， I_UT 的值越小。

3.6 努塞爾數(shù) （Nu）

式中： h 是平均傳熱系數(shù)； D_c 為流道水力直徑； k 是流體的導(dǎo)熱系數(shù)。

4網(wǎng)格無關(guān)性驗(yàn)證和模型驗(yàn)證

數(shù)值模擬精度與計(jì)算效率的平衡是仿真模擬的關(guān)鍵問題。針對平行流道模型，在熱通量為5000W/m² 、入口流速為 0.2m/s 的工況下，設(shè)置了6種網(wǎng)格數(shù)量進(jìn)行對比試驗(yàn)，網(wǎng)格1、網(wǎng)格2和網(wǎng)格3、網(wǎng)格4、網(wǎng)格5、網(wǎng)格6分別對應(yīng)的網(wǎng)格單元數(shù)為202145、363865、525585、687305、849025、1 010 745 。如圖4所示，當(dāng)網(wǎng)格數(shù)達(dá)到687305（網(wǎng)格4）時(shí)，繼續(xù)增加至849025（網(wǎng)格5）僅導(dǎo)致最大溫度降低 0.03% 以下、壓降升高 2.05% ，表明此時(shí)網(wǎng)格已經(jīng)滿足仿真需要。為了平衡精度和計(jì)算效率，選取網(wǎng)格4進(jìn)行后續(xù)仿真。

式中： q 是對流換熱量； h 是平均傳熱系數(shù)； A 是換熱面積； T_w 是壁面溫度； T_f 是冷卻液溫度，通常取冷卻液入口溫度。

為檢驗(yàn)數(shù)值模擬結(jié)果的可靠性，在一致工況下將仿真數(shù)據(jù)與BAEK等[的模型結(jié)果進(jìn)行對比。如圖5所示，在5組不同冷卻液質(zhì)量流量條件下，數(shù)值模擬結(jié)果與理論值呈現(xiàn)良好一致性，在冷卻液質(zhì)量流量 2.0E-06m³/s 時(shí)，相對誤差達(dá)到峰值6.39% 。該驗(yàn)證結(jié)果證實(shí)了數(shù)值模擬方法在質(zhì)子交換膜燃料電池?zé)峁芾矸治鲋械挠行耘c可信度。

5 結(jié)果與分析

5.1 速度云圖分析

如圖6所示的速度分布云圖，直觀呈現(xiàn)了不同冷卻流道內(nèi)冷卻液的流動(dòng)特性差異，而這種差異直接影響電池的溫度控制能力與傳熱效率。傳統(tǒng)平行流道受限于直線型結(jié)構(gòu)設(shè)計(jì)，流體流動(dòng)路徑單一且流道結(jié)構(gòu)缺乏變化，導(dǎo)致整體流速偏低且分布不均。這種平緩的流動(dòng)狀態(tài)導(dǎo)致流體在流道內(nèi)停留時(shí)間較長，尤其在中心區(qū)域形成了明顯的低速流動(dòng)死區(qū)，熱量難以及時(shí)被帶走，導(dǎo)致平行流道的傳熱性能不佳。

相比之下，3種新型流道通過結(jié)構(gòu)創(chuàng)新顯著改善了流體流動(dòng)特征。新型流道A和B通過引入轉(zhuǎn)彎結(jié)構(gòu)，顯著改變了流體的流動(dòng)路徑。在彎折區(qū)域，流體速度迅速增加，形成多個(gè)高速區(qū)域。這些高速流動(dòng)區(qū)域加劇了流體湍動(dòng)，破壞了層流邊界層，從而強(qiáng)化了流體與壁面之間的傳熱效率。此外，其分支結(jié)構(gòu)促進(jìn)了不同區(qū)域間的流體混合與熱交換，減少了流動(dòng)死區(qū)，進(jìn)一步提升了冷卻效果。而新型流道C通過延長流動(dòng)路徑并增加彎折次數(shù)，進(jìn)一步提升了流道整體的平均流速，使其成為4種方案中高流速區(qū)域占比最大的流道結(jié)構(gòu)。

5.2 溫度云圖分析

4種流道在流道中心橫截面上的溫度分布云圖如圖7所示，它清晰呈現(xiàn)了不同流道結(jié)構(gòu)對冷卻板溫度場的影響差異。平行流道的溫度分布呈現(xiàn)明顯的中心高溫聚集特征，當(dāng)入口流速為 0.2m/s 時(shí)，最高溫度達(dá)327.36K，低溫區(qū)域僅局限于流道兩側(cè)壁面附近。這種溫度分布是由于平行流道內(nèi)流體流動(dòng)平緩，層流邊界層較厚，熱量主要依賴壁面導(dǎo)熱傳遞，導(dǎo)致中心區(qū)域熱量積聚，溫度均勻性較差。

與平行流道相比，3種新型流道均表現(xiàn)出更加均勻的溫度分布。新型流道A的高溫區(qū)域面積較平行流道有所縮減，其低溫區(qū)域向流道中心擴(kuò)展，最高溫度降至 324.69K 。新型流道B的溫度均勻性進(jìn)一步提升，最高溫度為 324.31K 。其結(jié)構(gòu)在多個(gè)彎折處引發(fā)的流動(dòng)分離與再附著，形成周期性的高低速混合區(qū)域，這種流動(dòng)狀態(tài)促使熱量更均勻地分布，減少了局部熱堆積。

新型流道C的高溫區(qū)域最少且溫度也最低，降至 324.14K 。這主要?dú)w功于其多通道蛇形設(shè)計(jì)。這種設(shè)計(jì)通過多次彎折延長了流動(dòng)距離，增加了流體與壁面的接觸時(shí)間，使得熱量傳遞更加充分。同時(shí)，多通道布局使得流體分布更加均勻，減少了流動(dòng)死區(qū)，進(jìn)一步提高了傳熱效率，使其傳熱性能最佳。

5.3 壓降云圖分析

平行流道及3種新型流道在流道中心橫截面上的壓降分布云圖如圖8所示。從整體來看，不同流道結(jié)構(gòu)內(nèi)的壓降分布差異顯著，平行流道的壓降最低，3種新型流道的壓降均高于平行流道。平行流道的壓降分布呈現(xiàn)出沿流動(dòng)方向逐漸降低的特征，壓力梯度相對平緩，最大壓降為 155.13Pa 。這主要?dú)w因于平行流道內(nèi)流體流動(dòng)順暢，流體速度較低，阻力損失較小。

新型流道A的壓降分布呈現(xiàn)出較小局部高壓區(qū)域，最大壓降達(dá)到 172.77Pa 。這是由于流體在彎折處需要改變流動(dòng)方向，產(chǎn)生額外的流動(dòng)阻力，同時(shí)流體交匯處的流體沖擊與混合效應(yīng)也會(huì)導(dǎo)致壓力損失增加。

新型流道B由于分支通道的延長，壓降相對新型流道A進(jìn)一步升高，最高達(dá)到 187.64Pa 。新型流道C的壓降分布呈現(xiàn)出較為復(fù)雜的變化趨勢，最大壓降為 203.98Pa ，是4種流道中壓降最高的。這是由于其多通道蛇形設(shè)計(jì)，流體需要經(jīng)過多次方向改變，且其分支通道較長，流體在流動(dòng)過程中持續(xù)克服摩擦阻力，能量不斷損失，導(dǎo)致流動(dòng)阻力大幅增加。然而，這種設(shè)計(jì)通過增強(qiáng)流體擾動(dòng)，增強(qiáng)了傳熱效果，雖然壓降增大，但在一定程度上提高了冷卻性能。

5.4熱通量對冷卻性能的影響

如圖9和圖10所示，隨著熱通量從 4000W/m² 增加到 8 000W/m² ，平行流道的最大表面溫度從324.47K 上升到 336.04K ，溫差從 10.29K 增加到 20.68K ，均為4種流道中最高。這表明在熱負(fù)荷增加時(shí)，平行流道的溫度控制能力不足，溫度均勻性變差。相比之下，3種新型流道的最大溫度和溫差均小于平行流道。其中，新型流道C隨著熱通量的升高，最大溫度從321.91K上升到 330.83K 溫差從7.03K增加到 14.07K 。其最大溫度和溫差均為4種流道中最低，這表明新型流道在熱通量增加時(shí)，能更有效地控制最高溫度和溫差，維持更好的溫度均勻性。

如圖11和圖12所示，從平均溫度和IUT隨熱通量升高的變化趨勢來看，各流道的平均溫度隨熱通量增加而顯著上升，IUT則呈下降趨勢。這說明各流道隨熱通量的增加，溫度分布的均勻性逐漸降低。而新型流道的平均溫度始終低于平行流道，IUT始終高于平行流道，表明新型流道在提高溫度均勻性方面優(yōu)于平行流道。而在新型流道中，隨著熱通量的增加，新型流道C的平均溫度最低，新型流道A的平均溫度最高;新型流道B的IUT最高，新型流道A的IUT最小。

如圖13和圖14所示，各流道的平均傳熱系數(shù)（h） 和努塞爾數(shù) （Nu） 隨熱通量增加而增大。但平行流道的 Nu 和 h 均最小。這表明平行流道的傳熱能力提升有限，難以滿足高熱負(fù)荷需求。而新型流道的 h 和 Nu 增幅更為顯著。例如，新型流道C的 h 隨熱通量增加從 93.68W/（m²?K） 增加到161.49W/（m²?K），Nu 從156.14增加到269.15，為4種流道中表現(xiàn)最佳的案例，這表明新型流道能更有效地增強(qiáng)傳熱，適應(yīng)熱負(fù)荷的增加。

5.5冷卻液流速對冷卻性能的影響

如圖15和圖16所示，隨著冷卻液流速從0.2m/s 增加到 1.0m/s ，各流道的最大表面溫度和溫差均逐漸減小，這是由于冷卻液流速增加時(shí)，流體與壁面的相對運(yùn)動(dòng)速度加快，對流換熱能力得以增強(qiáng)。隨著流速的增加，平行流道的最大表面溫度從 327.36K 下降到319.12K，溫差從 12.88K 減少到 5.11K 。這表明增加流速能有效提升平行流道的冷卻性能，提高溫度均勻性。相比之下，3種新型流道的最大溫度和溫差降低幅度更大，溫度控制能力更優(yōu)。其中，新型流道C在 0.2m/s 的最大溫度僅為 324.14K ，溫差為 8.79K ，分別較平行流道的最大溫度和溫差降低了3.22K和 4.09K 為4種流道中表現(xiàn)最佳的案例。這說明新型流道對流速變化的響應(yīng)更為敏感，冷卻性能提升更顯著。

如圖17和圖18所示，從平均溫度和IUT的變化趨勢來看，各流道的平均溫度隨流速增加而降低，IUT則呈上升趨勢。各新型流道的平均溫度較平行流場更小，IUT也更高。新型流道C的平均溫度隨著流速增長從320.24K下降到315.02K，IUT從0.99661上升到0.99906，為表現(xiàn)最佳的案例。相比之下，平行流道的IUT從0.99531上升到0.99813，為表現(xiàn)最差的案例，表明新型流道在提升溫度均勻性方面優(yōu)于平行流道。

如圖19和圖20所示，各流道的 Nu 和 h 均隨流速增加而增大。平行流道的 Nu 和 h 在各流速下均為最小，這表明平行流道的傳熱能力最差。而新型流道的 h 和 Nu 隨著流速的變大增幅更為顯著。其中，新型流道C的 h 隨流速增大從 112.54W/（m²?K） 增加到 130.78W/（m²?K），Nu 從187.57增加到217.97，為表現(xiàn)最好的案例。這是由于新型流道的多通道蛇形設(shè)計(jì)使流體在流動(dòng)中多次改變方向，產(chǎn)生二次流動(dòng)和渦旋，擴(kuò)大了流體與壁面的接觸面積，能更有效地增強(qiáng)傳熱。

如圖21所示，從壓降的變化趨勢來看，各流道的壓降隨流速增加而變大，而各新型流道的壓降在各流速下始終大于平行流道。其中，新型流道C的壓降最大，在 0.2m/s 時(shí)為 203.98Pa ，新型流道B次之，而平行流道最小，壓降為 156.13Pa 。由此可見，新型流道在冷卻性能有較大提升的同時(shí)，壓降也僅較平行流道略有增加，仍然以較低的壓降實(shí)現(xiàn)了較好的傳熱性能，

6 結(jié)論

本研究提出了3種多通道分流結(jié)構(gòu)的新型流道設(shè)計(jì)，旨在提升PEMFC的熱管理性能。探究了不同流道設(shè)計(jì)、熱通量和冷卻液流速對熱管理性能的影響，得出了以下結(jié)論：

1）基于多通道分流結(jié)構(gòu)的新型流道（A、B、C）在傳熱效率與溫度均勻性方面均優(yōu)于傳統(tǒng)平行流道，其中新型流道C的綜合性能最為突出。在0.2m/s 的入口流速工況下，相比平行流道，新型流道C的最大表面溫度降低了3.22K，溫差減小了4.09K ，而壓降僅增加了 48.85Pa ，體現(xiàn)了其在較小流動(dòng)阻力增幅下實(shí)現(xiàn)高效散熱的優(yōu)勢。

2）隨著熱通量從 4000W/m² 增加到8 000W/m² 各流道的最高表面溫度、溫差及平均溫度均呈現(xiàn)增長趨勢，同時(shí)努塞爾數(shù)與表面?zhèn)鳠嵯禂?shù)也同步升高。值得注意的是，溫度均勻性指數(shù)隨熱負(fù)荷增加而持續(xù)降低，反映出高熱通量工況下流道溫度場分布不均的問題加劇。這一現(xiàn)象說明，在應(yīng)對高熱負(fù)荷場景時(shí)，冷卻流道設(shè)計(jì)需優(yōu)先考慮溫度場分布均勻性的優(yōu)化策略，以抑制局部高溫積聚并提升整體熱管理能力。

3）冷卻液流速的增加對冷卻性能有顯著增強(qiáng)作用。隨著入口速度的增加，最大表面溫度與溫差均呈現(xiàn)下降趨勢，溫度均勻性指數(shù)及表面?zhèn)鳠嵯禂?shù)逐步增加，實(shí)現(xiàn)了更均勻的溫度分布和更好的冷卻效果，然而，壓降隨流速升高呈非線性增長。因此，在工程應(yīng)用中需綜合考量冷卻效率與能耗成本，確定適宜的流速條件實(shí)現(xiàn)熱管理性能與能耗的平衡，以達(dá)最佳綜合效益。

參考文獻(xiàn)：

[1]WANGZH，HUKH，ZHANGJ，etal.Gas-liquidmass transfer characteristics of a novel three-dimensional flow fieldbipolarplate forlaser additive manufacturing ofproton exchange membrane fuel cell[J].Renewable Energy， 2023，212：308-309.

[2]LIUMX，F(xiàn)ANWX，LUGL，etal.Study on mass transfer enhancement of locally improved structures and theapplication in serpentine and parallel flow fieldsof

PEM fuel cells[J]．International Journal of Hydrogen Energy，2023，48（50）：19248-19261.

[3]UZUNDURUKAN A，BILGILI M，DEVRIM Y. Examination of compression effects on PEMFC performance by numerical and experimental analyses[J]. International Journal of Hydrogen Energy，2020，45（60）：35085-35096.

[4]HUANGYC，XIAO XL，KANGHF，et al.Thermal management of polymer electrolyte membrane fuel cells ： a critical reviewof heattransfer mechanisms，coling approaches，and advanced cooling techniques analysis[J]. EnergyConversion and Management，2022，254： 115221. 1-115221. 35.

[5] ZHAO J， JIAN Q F， HUANG Z P. Experimental study on heat transfer performance of vapor chambers with potential applications in thermal management of proton exchange membrane fuel cels[J].Applied Thermal Engineering， 2020，180：115847.1-115847.12.

[6]ZHANG G S，KANDLIKAR S G. A critical review of cooling techniques in proton exchange membrane fuel cell stacks[J].International Journal of Hydrogen Energy， 2012，37（3）：2412-2429.

[7]LASBETY，AUVITYB，CASTELAINC，etal.Thermal and hydrodynamic performances of chaotic mini-channel： application to the fuel cellcooling[J]．Heat Transfer Engineering，2007，28（8/9）： 795-803.

[8]SOLTANI A，POURFALLAH M，SEYED SABOUR S M J.Thermal enhancement of PEM fuel cell cooling with novel configurations of Tesla valve and hybrid nanofluids ： a numerical study[J]. International Journal of Hydrogen Energy，2024，69：1263-1275.

[9]AFSHARI E， ZIAEI-RAD M， DEHKORDI M M. Numerical investigation on a novel zigzag-shaped flow channel design for cooling plates of PEM fuel cells[J]. Journal of the Energy Institute，2017，90（5）：752-763.

[10]張，劉良旭.基于燃料電池汽車吸收式熱泵的PEMFC冷卻流道強(qiáng)化換熱性能研究［J].節(jié)能，2024，43 （1） ：22-25.

[11]柳奐宇．質(zhì)子交換膜燃料電池流場設(shè)計(jì)及水熱管理 [D].合肥：合肥工業(yè)大學(xué)，2021.

[12]王琦，徐曉明，司紅磊，等.波形結(jié)構(gòu)冷卻流道對燃料 電池?zé)峁芾硐到y(tǒng)性能影響研究［J].中國工程機(jī)械學(xué) 報(bào)，2022，20（2）：95-100.

[13]IFFI CINDIMH IFFNt al

flow-field design for passive air-cooled polymer electrolyte membrane（PEM） fuel cell stacks[J]. International Journal of Hydrogen Energy，2020，45（20） ： 11704-11713.

[14]梅楠.梯形流道質(zhì)子交換膜燃料電池氣體流道及冷卻 流道研究[D].鎮(zhèn)江：江蘇大學(xué)，2022.

[15]宋江南.質(zhì)子交換膜燃料電池冷卻流場設(shè)計(jì)及優(yōu)化 [D].貴陽：貴州大學(xué)，2023.

[16]BAEK S M，YU SH，NAMJH，et al. A numerical study on uniform cooling of large-scale PEMFCs with different coolant flow field designs[J]．Applied Thermal Engineering，2011，31（8/9）：1427-1434.

[17]王珊．質(zhì)子交換膜燃料電池冷卻流場設(shè)計(jì)及熱模擬 [D].天津：天津大學(xué)，2019.

[18]徐鑫，王珊，韓冬林，等．PEMFC 冷卻流道數(shù)值模擬 及實(shí)驗(yàn)研究[J].電源技術(shù)，2022，46（5）：509-513.

[19]HUANG T M， XIAO YF，YI D X， et al. Enhancing heat transfer in proton exchange membrane fuel cells through obstructed cooling channels：a comprehensive study[J].Applied Thermal Engineering，2024，250： 123326.1-123326. 18.

[20] PENG C，GU H L， ZHANG G X， et al. Numerical study on heat transfer enhancement of a proton exchange membrane fuel cell with the dimpled cooling channel[J]. International Journal of Hydrogen Energy，2023，48（8） ： 3122-3134.

[21]黃浩，鞠洪玲.PEMFC 新型冷卻流道的數(shù)值模擬[J]. 電源技術(shù)，2024，48（1） ：150-155.

[22]FAN X，MENG C，YANG YW，et al. Numerical optimization of the cooling effect of a bionic fishbone channel liquid cooling plat for alarge prismatic lithium-ion battery pack with high discharge rate[J]. Journal of Energy Storage，2023，72（A）：108239.1-108239.13.

[23]CHEN X，CHAI F S，HU S L，et al. Design of PEMFC bipolar plate cooling flow field based on fractal theory [J].Energy Conversion and Management： X，2023， 20：100445.1-100445.12.

[24]BEHROU R，PIZZOLATO A，F(xiàn)ORNER-CUENCA A. Topology optimization as a powerful tool to design advanced PEMFCs flow fields[J]. International Journal of Heat and Mass Transfer，2019，135：72-92.

（責(zé)任編輯：曾 晶）

Design and Performance Study of Novel Cooling Channels for Proton Exchange Membrane Fuel Cells

QI Meiling，HUANG Ying （College of Mechanical Engineering，Guizhou University，Guiyang 55oo25，China）

Abstract：The design ofcooling channels has a significant impact on the thermal management performance of the proton exchange membrane fuel cell（PEMFC）．This study addresses the low heat dissipation performance of traditional parallel channels and proposes three new channel designs to enhance the thermal management performance of PEMFC. The results show that the new channels outperform the paralel channels in terms of cooling effciency and temperature uniformity.Among them，the new channel C has the best performance.At a coolant flow rate of 0.2m/s ， its maximum surface temperature is 3.22K lower than that of the parallel channel， and the temperature difference isreduced by 4.09K ，demonstrating the advantage of achieving efficient heat dissipationwitharelativelysmall increase in pressuredrop.Inaddition，this study further explores theeffects of heat flux and coolant flow rate on thermal management performance.Itis found that high heat flux intensifies the non-uniformity of temperature distribution，highlighting theimportance ofoptimizing thermal distribution uniformity under high heat load.While increasing theflow rate can significantly enhance thecooling effect，it leads to a sharp increase in pressre drop. Therefore，in practical applications，a balance must be struck between cooling eficiency and energy consumption. The new channel designs proposed in this study can effectively improve the thermal management performance of PEMFC，providing theoretical basisand technical support for promoting the engineering application of high-power-density PEMFC.

KeyWords： proton exchange membrane fuel cell（PEMFC）；cooling channel; thermal management; numerica simulation；temperature uniformity