張斌 莊才敖 鄭宇

摘要：質(zhì)子交換膜燃料電池作為燃料電池中重要的一種類型，可應(yīng)用于汽車及小規(guī)模的發(fā)電站與便攜式移動(dòng)能源，是當(dāng)前新能源領(lǐng)域的研究熱點(diǎn)。主要針對(duì)車用質(zhì)子交換膜燃料電池的三種仿真參數(shù)模型就當(dāng)前國(guó)內(nèi)外研究進(jìn)展進(jìn)行論述。

Abstract： Proton exchange membrane fuel cells， as an important type of fuel cells， can be applied to automobiles and small-scale power stations and portable mobile energy sources， which is a current research focus in the field of new energy. This paper mainly discusses the current research progress at home and abroad for three simulation parameter models of proton exchange membrane fuel cells for vehicles.

關(guān)鍵詞：質(zhì)子交換膜燃料電池（PEMFC）;仿真模型;技術(shù)路線

Key words： proton exchange membrane fuel cell （PEMFC）;simulation model;technical route

中圖分類號(hào)：TM911.4? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-957X（2021）06-0210-03

0? 引言

燃料電池作為一項(xiàng)具有廣闊應(yīng)用前景的能源技術(shù)，具有高效率、低排放、低噪音及模塊化等優(yōu)勢(shì)[1]。質(zhì)子交換膜燃料電池作為其中極其重要的一個(gè)種類，有望大規(guī)模應(yīng)用于汽車、發(fā)電、移動(dòng)電源等眾多領(lǐng)域，減輕當(dāng)前面臨的汽車尾氣排放及燃煤污染物排放所帶來(lái)的環(huán)境問題與溫室氣體造成的氣候問題[2]。質(zhì)子交換膜燃料電池的單電池結(jié)構(gòu)主要包括陽(yáng)極流場(chǎng)、膜電極結(jié)構(gòu)（陽(yáng)極擴(kuò)散層、陽(yáng)極催化層、質(zhì)子交換膜、陰極催化層、陰極擴(kuò)散層）、陰極流場(chǎng)[3]。其中陰陽(yáng)極流場(chǎng)結(jié)構(gòu)影響進(jìn)入擴(kuò)散層流體的壓力、流速以及電流傳導(dǎo)與熱傳導(dǎo)等因素，一直是燃料電池研究的熱點(diǎn)[4]。而核心組件即是質(zhì)子交換膜和膜兩側(cè)的多孔電極，質(zhì)子交換膜必須保證足夠高的質(zhì)子導(dǎo)電性以及能夠化學(xué)和機(jī)械穩(wěn)定性[5]。多孔電極即發(fā)生電化學(xué)反應(yīng)的催化劑薄層，保證反應(yīng)氣體、電子、質(zhì)子都能參與反應(yīng)。氣體擴(kuò)散層則為反應(yīng)氣體從流場(chǎng)進(jìn)入反應(yīng)區(qū)域提供路徑并與雙極板相連使電子形成電路以及為膜電極結(jié)構(gòu)提供機(jī)械支撐。雙極板雖然不出現(xiàn)在單電池中，卻是構(gòu)成電池組的因素，將一個(gè)電池的陽(yáng)極與相鄰電池的陰極進(jìn)行電氣連接并隔絕氣體防止?jié)B透又為電池組提供了結(jié)構(gòu)支撐以及起到傳導(dǎo)熱量的作用。燃料電池各組件的材料特性、反應(yīng)機(jī)理等比較復(fù)雜，試驗(yàn)成本較高，數(shù)值仿真技術(shù)可以優(yōu)化燃料電池系統(tǒng)的設(shè)計(jì)與匹配，對(duì)試驗(yàn)研究起到指導(dǎo)與輔助作用[6]。常用的仿真模型包括分布參數(shù)模型、集總參數(shù)模型以及混合參數(shù)模型[7]。

1? 分布參數(shù)模型

分布參數(shù)模型即以燃料電池仿真涉及的各種守恒方程，包括質(zhì)量守恒、動(dòng)量守恒、電荷守恒、能量守恒、組分守恒等守恒方程來(lái)描述燃料電池內(nèi)部參數(shù)[7]，包括電池內(nèi)部流體（氫氧燃料電池主要為H2、O2;直接甲醇燃料電池為甲醇溶液）的壓力、流速、濃度、溫度等參數(shù)以及液態(tài)水與電流密度的分布情況等。分布參數(shù)模型可以詳細(xì)的描繪電池內(nèi)部情況，但一般忽略了燃料電池配套設(shè)施比如燃料供給系統(tǒng)、冷卻系統(tǒng)等的影響。常用的商用仿真軟件包括Fluent、COMSOL Multiphysics等。

1.1 流場(chǎng)結(jié)構(gòu)對(duì)PEMFC影響

在分布參數(shù)模型下，羅鑫等基于COMSOL Multiphysics軟件比較了不同流場(chǎng)形式對(duì)PEMFC性能的影響，通過對(duì)陰陽(yáng)極流體濃度與電流密度的比較發(fā)現(xiàn)交指型流道的排水效果最好，功率密度也最高[8]，見圖1。蛇形流場(chǎng)是目前廣泛應(yīng)用的流場(chǎng)形式，單蛇形流場(chǎng)由于流道長(zhǎng)度較長(zhǎng)，壓降較大[9];不如多通道蛇形流場(chǎng)適用于面積較大的流場(chǎng);相對(duì)于二通道的流場(chǎng)結(jié)構(gòu)，三通道流場(chǎng)結(jié)構(gòu)會(huì)使反應(yīng)氣體及電流密度分布更加均勻[10]。

1.2 PEMFC內(nèi)部傳質(zhì)現(xiàn)象

水的凝結(jié)速度可以表示為：

其中rw表示水凝結(jié)速度，K表示絕對(duì)滲透率，ρl表示液態(tài)水密度，s為水的飽和度，pc表示孔隙介質(zhì)毛細(xì)壓力，ε為孔隙率，μl表示液態(tài)水粘度。

陽(yáng)極流道入口與陰極流道入口的位置會(huì)影響PEMFC內(nèi)部水濃度分度。陽(yáng)極流道入口與陰極流道入口成對(duì)角線布置（逆流）時(shí)，平均水含量更高。因此采用逆流布置可以提高燃料電池的自增濕效果[11]，見圖2。但當(dāng)燃料電池內(nèi)水濃度過高時(shí)產(chǎn)生的“水淹”現(xiàn)象會(huì)對(duì)電池性能產(chǎn)生嚴(yán)重的不利影響，此時(shí)可以通過優(yōu)化流道結(jié)構(gòu)，如采用變寬度流道結(jié)構(gòu)[11];以及改變進(jìn)氣相對(duì)濕度[12]等手段加以控制。

1.3 PEMFC特殊工況分析

對(duì)于燃料電池的一些特殊工況需要進(jìn)行動(dòng)態(tài)的研究分析，比如啟動(dòng)工況與變負(fù)載工況等。啟動(dòng)工況常用電流密度隨時(shí)間變化曲線來(lái)表征，電流密度在短時(shí)間內(nèi)上升并達(dá)到穩(wěn)定狀態(tài)則啟動(dòng)成功[13]，見圖3。流場(chǎng)結(jié)構(gòu)形式[14]、進(jìn)氣參數(shù)、電池堆內(nèi)單電池個(gè)數(shù)[15]、質(zhì)子交換膜特性[16]等都會(huì)對(duì)燃料電池啟動(dòng)性能產(chǎn)生影響。

當(dāng)負(fù)載發(fā)生階躍變化時(shí)，電池電壓容易產(chǎn)生“overshoot”（陡升）與“undershoot”（陡降）現(xiàn)象[17]。如燃料電池汽車在加載瞬間電池電壓會(huì)發(fā)生陡降的現(xiàn)象。這種不利現(xiàn)象會(huì)降低電池堆的性能表現(xiàn)甚至導(dǎo)致電池堆的不可逆損害[7]。這種現(xiàn)象主要由氣體供應(yīng)滯后和電池內(nèi)部液態(tài)水的積累所造成[7]。當(dāng)電池電流密度突變時(shí)還會(huì)導(dǎo)致電池溫度的動(dòng)態(tài)變化[18]。

2? 集總參數(shù)模型

集總參數(shù)模型仿真速度快、操作簡(jiǎn)便，可以用來(lái)分析燃料電池的動(dòng)態(tài)響應(yīng)性與系統(tǒng)控制，但無(wú)法得到燃料電池內(nèi)部參數(shù)的分布情況[7]。最常使用的仿真軟件是Matlab/Simulink仿真平臺(tái)。

對(duì)于PEMFC的系統(tǒng)控制主要包括了機(jī)理模型與經(jīng)驗(yàn)?zāi)Ｐ蚚19]以及以經(jīng)驗(yàn)?zāi)Ｐ蜑榛A(chǔ)，考慮質(zhì)子交換膜燃料電池配套系統(tǒng)的半經(jīng)驗(yàn)半機(jī)理模型[20]?？梢曰谒惴刂迫剂想姵囟训臏囟萚21]、流道內(nèi)流體流量的控制[22];利用集總參數(shù)模型研究PEMFC結(jié)構(gòu)與參數(shù)對(duì)性能的影響如流道布置方式對(duì)傳熱傳質(zhì)效果的影響[23];采用動(dòng)態(tài)集總參數(shù)模型分析PEMFC系統(tǒng)的動(dòng)態(tài)特性[24]以及建立PEMFC輔助系統(tǒng)的模型分析[7]。

3? 混合參數(shù)模型

仿真模型的研究由一維拓展到二維和三維模型，從PEMFC的特定區(qū)域拓展到單電池以及電池堆的整體性能仿真[6]。對(duì)于PEMFC的仿真建模而言，分布參數(shù)與集總參數(shù)模型各有優(yōu)勢(shì)也各有不足。因此，將兩種參數(shù)模型相結(jié)合的混合參數(shù)模型在PEMFC的仿真分析中更具優(yōu)勢(shì)。首先利用分布參數(shù)模型建立燃料電池堆的仿真模型，可以得到電池內(nèi)部參數(shù)的分布情況;利用集總參數(shù)模型建立輔助系統(tǒng)如空壓機(jī)、增濕器、燃料供給系統(tǒng)等的仿真模型[7]，便可以充分發(fā)揮兩種參數(shù)模型的優(yōu)勢(shì)，得到PEMFC工作時(shí)，隨著運(yùn)行參數(shù)的改變，內(nèi)部物理量的動(dòng)態(tài)響應(yīng)。常用的仿真方法即利用Simulink的S-函數(shù)與Fluent軟件日志文件建立兩者的接口[25]，模擬更接近真實(shí)的PEMFC系統(tǒng)運(yùn)行狀態(tài)。

目前有關(guān)燃料電池領(lǐng)域的混合參數(shù)模型研究還很匱乏，但基于混合參數(shù)模型仿真模擬的巨大優(yōu)勢(shì)，可以想見未來(lái)的研究潛力巨大。

4? 質(zhì)子交換膜燃料電池在汽車上的應(yīng)用

當(dāng)下的新能源汽車以純電動(dòng)汽車和插電式混合動(dòng)力汽車為主，儲(chǔ)能裝置為鋰電池，受鋰電池技術(shù)的影響，純電動(dòng)汽車的續(xù)航里程短、冬季續(xù)航斷崖式下跌和充電時(shí)間長(zhǎng)，其便利性一直不如傳統(tǒng)燃油車。鋰電池隨著充放電次數(shù)的增加，電池容量也會(huì)不斷衰減，最直接的表現(xiàn)就是純電動(dòng)汽車使用三年后，續(xù)航里程嚴(yán)重縮短。鋰電池報(bào)廢后，也會(huì)產(chǎn)生較大的環(huán)境污染。

質(zhì)子交換膜燃料電池汽車解決了大部分純電動(dòng)汽車上的問題，燃料電池的是一種高效的能源轉(zhuǎn)換裝置，實(shí)際轉(zhuǎn)換效率在40-60%[26]。主要參與的反應(yīng)物是氧氣和氫氣，氧氣可直接從空氣中獲取，只需另外加注氫氣，反應(yīng)后的主要產(chǎn)物是水，是非常清潔的能源裝置。我國(guó)氫燃料電池汽車?yán)塾?jì)推廣應(yīng)用超過2000輛，投入運(yùn)行加氫站有12 座，且在北上廣等地均開展了示范應(yīng)用[27]。車用質(zhì)子交換膜燃料電池在新能源汽車的應(yīng)用前景十分廣闊，能夠有利解決環(huán)境污染問題，同時(shí)改善我國(guó)的能源結(jié)構(gòu)。

參考文獻(xiàn)：

[1]Frano Barbir. PEM燃料電池：理論與實(shí)踐[M].二版.機(jī)械工業(yè)出版社，2016：130-131.

[2]Bockris JO'M， Veziroglu TN. A Solar-Hydrogen Energy System for Environmental Compatibility. Environment Conservation 1985;12（No. 2）：105-18.

[3]Sasaki K， Li H W， Hayashi A， et al. Hydrogen Energy Engineering[M]. Springer， 2016.

[4]陳士忠，夏忠賢，王藝澄，張旭陽(yáng)，吳玉厚.蛇形流場(chǎng)PEMFC性能影響因素的數(shù)值模擬[J].電源技術(shù)，2017，41（02）：230-233.

[5]Gottesfeld S， Zawodzinski TA. Polymer Electrolyte Fuel Cells. In： Alkire RC， Gerischer H， Kolb DM， Tobias CW， editors. Advances in Electrochemical Science and Engineering， vol. 5. New York： Wiley-VCH; 1997.

[6]翟雙，周蘇，陳鳳祥，等.質(zhì)子交換膜燃料電池分布參數(shù)模型數(shù)值仿真研究進(jìn)展[J].同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，40（6）：932-936，954.

[7]翟雙，周蘇.質(zhì)子交換膜燃料電池仿真方法及若干現(xiàn)象研究[M].同濟(jì)大學(xué)出版社，2018：24-28.

[8]羅鑫，陳士忠，夏忠賢.不同流場(chǎng)的PEMFC性能研究[J].電池，2017，47（04）：208-211.

[9]Watkins DS， Dircks KW，Epp DG. Novel Fuel Cell Fluid Flow Field Plate; 1991. U.S. Patent 4， 988， 583.

[10]趙康，趙志軍，程驍.PEMFC不同流道蛇形流場(chǎng)的三維建模仿真[J].電池，2018，48（05）：318-321.

[11]張卿雷.基于CFD方法的質(zhì)子膜燃料電池內(nèi)部傳質(zhì)分析及優(yōu)化設(shè)計(jì)[D].上海交通大學(xué)，2017.

[12]王忠豪.質(zhì)子交換膜燃料電池水管理模型研究[D].重慶理工大學(xué)，2018.

[13]袁峰.低溫環(huán)境下質(zhì)子交換膜燃料電池冷啟動(dòng)性能研究[D].浙江工業(yè)大學(xué)，2019.

[14]Santamaria AD， Bachman J， Park JW. Cold-start of parallel and interdigitated flow-field polymer electrolyte membrane fuel cell [J]. Electrochimicaa Acta， 2013， 107（10）：327-338.

[15]Zhou Y， Luo Y， Yu S， et al. Modeling of cold start processes and performance optimization forproton exchange membrane fuel cell stacks[J]. Journal of Power Sources， 2014， 247（2）： 738-748.

[16]Lin R， Lin X， Weng Y， et al. Evolution of thermal drifting during and after cold start of proton exchange membrane fuel cell by segmented cell technology[J]. International Journal of Hydrogen Energy， 2015， 40（23）：7370-7381.

[17]Meng H. Numerical investigation of a PEM fuel cell using a two-phase non-intermal mixed-domain model[J]. Journal of Power Sources， 2007， 171： 738.

[18]周蘇，李壯運(yùn)，翟雙，陳鳳祥.PEMFC電堆建模及特殊工況下的動(dòng)態(tài)分析[J].太陽(yáng)能學(xué)報(bào)，2011，32（07）：1123-1128.

[19]王金龍，王登峰，賈迎春，陳靜.質(zhì)子交換膜燃料電池?cái)?shù)學(xué)建模進(jìn)展[J].河北工業(yè)科技，2005，22（6）：315-320.

[20]郟怡穎，趙治國(guó).車用質(zhì)子交換膜燃料電池動(dòng)態(tài)建模與仿真分析[J].機(jī)械與電子，2012（04）：53-57.

[21]Li X， Deng Z H， Wei D， et al. Novel variable structure control for the temperature of pem fuel cell stack based on the dynamic thermal affine model[J]. Energy Conversion and Management， 2011， 52（11）： 3256-3274.

[22]Zhang J Z， Liu G D， Yu W S， et al. Adaptive control of the airflow of a pem fuel cell system[J]. Journal of Power Source， 2008， 179（2）： 649-659.

[23]Yu S S， Im S， Kim S， et al. A parametric study of the performance of a planar membrane humidifier with a heat mass exchange model for design optimization[J]. International Journal of Heat and Mass Transfer， 2011（54）：1344-1351.

[24]楊朵.質(zhì)子交換膜燃料電池動(dòng)態(tài)建模與仿真[C].中國(guó)自動(dòng)化學(xué)會(huì)系統(tǒng)仿真專業(yè)委員會(huì)、中國(guó)系統(tǒng)仿真學(xué)會(huì)仿真技術(shù)應(yīng)用專業(yè)委員會(huì).第18屆中國(guó)系統(tǒng)仿真技術(shù)及其應(yīng)用學(xué)術(shù)年會(huì)論文集（18th CCSSTA 2017）.中國(guó)自動(dòng)化學(xué)會(huì)系統(tǒng)仿真專業(yè)委員會(huì)、中國(guó)系統(tǒng)仿真學(xué)會(huì)仿真技術(shù)應(yīng)用專業(yè)委員會(huì)：中國(guó)自動(dòng)化學(xué)會(huì)系統(tǒng)仿真專業(yè)委員會(huì)，2017：87-91.

[25]賀明艷，周蘇，黃自萍，孫澎濤.基于SIMULINK/FLUENT的PEMFC系統(tǒng)的協(xié)同仿真[J].系統(tǒng)仿真學(xué)報(bào)，2011，23（01）：38-43.

[26]王一.燃料電池動(dòng)力裝置水熱管理系統(tǒng)研究及儲(chǔ)能方案論證[D].哈爾濱工程大學(xué)，2015.

[27]蔣誠(chéng).車用燃料電池?zé)峁芾硐到y(tǒng)設(shè)計(jì)開發(fā)與性能驗(yàn)證[D].吉林大學(xué)，2020.