質(zhì)子交換膜燃料電池大電流密度下運(yùn)行工況優(yōu)化研究*

張洪凱 詹志剛 何曉波 帥 露 隋邦杰 潘 牧

(武漢理工大學(xué)材料復(fù)合新技術(shù)國家重點(diǎn)實(shí)驗(yàn)室1) 武漢 430070)(武漢理工大學(xué)能源與動(dòng)力工程學(xué)院2) 武漢 430063) (武漢理工大學(xué)汽車工程學(xué)院3) 武漢 430070)(加拿大維多利亞大學(xué)機(jī)械系集成能源系統(tǒng)實(shí)驗(yàn)室4) 維多利亞市 V8W 2Y2)

質(zhì)子交換膜燃料電池大電流密度下運(yùn)行工況優(yōu)化研究*

張洪凱1,2)詹志剛1,2)何曉波1,2)帥 露2)隋邦杰3,4)潘 牧1)

(武漢理工大學(xué)材料復(fù)合新技術(shù)國家重點(diǎn)實(shí)驗(yàn)室1)武漢 430070)(武漢理工大學(xué)能源與動(dòng)力工程學(xué)院2)武漢 430063) (武漢理工大學(xué)汽車工程學(xué)院3)武漢 430070)(加拿大維多利亞大學(xué)機(jī)械系集成能源系統(tǒng)實(shí)驗(yàn)室4)維多利亞市 V8W 2Y2)

利用計(jì)算流體動(dòng)力學(xué)軟件，針對(duì)在大電流密度、陰極無加濕、陽極加濕，陰陽極有背壓等操作特點(diǎn)下運(yùn)行的PEM電池進(jìn)行數(shù)值模擬，綜合考慮電池性能和電流密度分布的均勻性，尋求最優(yōu)工況點(diǎn).結(jié)果表明，對(duì)于給定的操作壓力，存在一個(gè)最佳性能溫度，壓力增加，最佳溫度隨之增加，在壓力為100，200，300 kPa時(shí)，最佳性能溫度分別為75，60，50 ℃；由于陰極不加濕，在最佳性能溫度點(diǎn)膜含水量分布及電流密度分布不均勻；在較低溫度下運(yùn)行時(shí)電流密度分布均勻性隨著壓力的增加逐漸提高，但在較高溫度下運(yùn)行時(shí)則相反.

質(zhì)子交換膜燃料電池；電流密度；水傳輸；運(yùn)行工況；優(yōu)化

0 引 言

在化石能源面臨枯竭和環(huán)境污染日益嚴(yán)重的今天，尋找一種高效、無污染的新型能源成為了亟待解決的問題，質(zhì)子交換膜(PEM)燃料電池因其高效率、高功率密度、低污染等優(yōu)點(diǎn)廣泛受到人們的關(guān)注.PEM燃料電池可以應(yīng)用于汽車、航天、移動(dòng)充電設(shè)備等領(lǐng)域，被認(rèn)為市場(chǎng)潛力十分巨大[1].由于市場(chǎng)應(yīng)用的需求，PEM電池體積不斷減小而功率密度不斷增加，日本豐田汽車公司于2014年推出了MIRAI 燃料電池汽車，其電堆的質(zhì)量功率密度和體積功率密度達(dá)到2.0 kW/kg和3.1 kW/L[2].電池堆在大功率下運(yùn)行，電流密度往往也較高，可能高達(dá)1.5～2.5 A/cm2，因此，一般需要增加背壓，以提高反應(yīng)氣體濃度；同時(shí)電池生成水較多，合理利用生成水，可以簡(jiǎn)化或去掉加濕系統(tǒng).因此研究膜電極內(nèi)部水、氣傳輸及分布規(guī)律，尋找最優(yōu)工況點(diǎn)，對(duì)電池高效、穩(wěn)定、持久運(yùn)行，具有重要意義.

Su等[3-4]通過模擬方法探索了溫度、濕度、壓力、氣體過量系數(shù)等操作參數(shù)對(duì)電池性能的影響.Kim等[5]分析了反應(yīng)氣體加濕度對(duì)燃料電池性能的影響，認(rèn)為陽極充分加濕而陰極不加濕的條件下質(zhì)子交換膜仍能被充分潤濕.Amirinejad等[6]研究操作條件對(duì)電池性能的影響發(fā)現(xiàn)，在操作壓力較高的條件下提高溫度會(huì)使得電池性能提高.Santarelli等[7]通過實(shí)驗(yàn)發(fā)現(xiàn)，只有陽極加濕的燃料電池性能會(huì)隨著操作壓力的提高而明顯提高，同時(shí)還能夠提高電池穩(wěn)定性.Jang等[8]通過研究發(fā)現(xiàn)，陽極加濕溫度提高會(huì)使得電池性能提高.

文中針對(duì)在大電流密度、陰極無加濕、陽極加濕，陰陽極有背壓等操作特點(diǎn)下運(yùn)行的PEM電池，利用計(jì)算流體動(dòng)力學(xué)軟件，進(jìn)行數(shù)值模擬，尋找電池最佳工況點(diǎn).獲得的結(jié)論對(duì)于PEM電池的設(shè)計(jì)和操作具有參考意義.

1 計(jì)算模型

1.1 幾何模型

根據(jù)實(shí)際金屬板電池結(jié)構(gòu)，建立三維單流道電池模型，因其結(jié)構(gòu)對(duì)稱，為減少計(jì)算量，選取1/2流道區(qū)域建模，見圖1.外圍尺寸為50 mm×1.1 mm×1.435 mm，活性面積為55 mm2，氣體和冷卻水流道寬度與深度為0.5 mm×0.4 mm，集流板厚度0.1 mm，氣體擴(kuò)散層厚度0.2 mm，催化層厚度0.01 mm，膜厚度0.015 mm.在進(jìn)行網(wǎng)格靈敏性驗(yàn)證后，確定總網(wǎng)格總數(shù)為89 000個(gè).

圖1 幾何模型

1.2 控制方程

質(zhì)子交換膜燃料電池是一個(gè)多相、多尺度的復(fù)雜系統(tǒng).其內(nèi)部流道及多孔介質(zhì)內(nèi)氣體的擴(kuò)散、水的相變及流動(dòng)，催化層內(nèi)部的電化學(xué)反應(yīng)、水在質(zhì)子交換膜內(nèi)的傳遞都同時(shí)發(fā)生并且相互耦合.描述上述現(xiàn)象及過程的主要控制方程包括如下的質(zhì)量守恒方程、動(dòng)量守恒方程、組份守恒方程、能量守恒方程、電荷守恒方程等.

(1)

-εp+·(εμu)+Su

(2)

(3)

(4)

(5)

(6)

1.3 操作條件

文中主要研究陽極和陰極加濕度分別為100%，0%條件下大電流運(yùn)行時(shí)電池最優(yōu)性能.陰陽極過量系數(shù)分別為1.5和2.5，出口背壓50，100，150，200 kPa，操作溫度和反應(yīng)氣體的溫度為50，60，70，80，90 ℃.主要電化學(xué)參數(shù)及反應(yīng)氣體物性見表1[9-10]，電池主要部件材料物性見表2[11-12].

表1 主要物性參數(shù)

表2 材料參數(shù)

2 結(jié)果與討論

圖2為模擬與實(shí)驗(yàn)測(cè)試伏安曲線對(duì)比，由圖2可知,模擬值與實(shí)測(cè)值基本一致，同時(shí)也驗(yàn)證了該模型可靠性.

圖2 模擬與實(shí)驗(yàn)測(cè)試伏安曲線

圖3為1.5 A/cm2時(shí)不同操作條件下電池的電壓.隨著操作壓力的增加，因?yàn)榉磻?yīng)氣體濃度增加，電池性能也逐漸增加.在相同壓力下，存在一個(gè)最佳性能的操作溫度.例如，在壓力為300 kPa時(shí)，隨著溫度從50 ℃開始增加，電池性能漸漸增加，在75 ℃左右到達(dá)最大值，之后又逐漸降低，這是因?yàn)殡姵匦阅苁艿窖鯕鉂舛群湍櫇裥詢蓚€(gè)因素的影響.在溫度較低時(shí)，液態(tài)水容易生成，催化層中液態(tài)水相飽和度較高，使得膜潤濕性較好，內(nèi)阻較低；但同時(shí)較高的液態(tài)水相飽和度阻礙了反應(yīng)氣體的傳輸，不利于電化學(xué)反應(yīng)的進(jìn)行；溫度繼續(xù)增加，兩種作用的趨勢(shì)相反，膜逐漸變得過于干枯，內(nèi)阻增加，電池性能下降；在75 ℃左右時(shí)這種綜合作用得到了最佳的性能.

圖3 1.5 A/cm2時(shí)不同操作條件下電池性能

圖4為300 kPa下不同溫度時(shí)催化層中間平面液態(tài)水相飽和度分布.沿氣體流道從進(jìn)口到出口方向，水汽逐漸積累，相飽和度也逐漸增加；右半邊處于冷卻水流道下而左半邊處于氣體流道下，溫度差異導(dǎo)致右邊液相飽和度總體高于左半邊.在50 ℃時(shí)，液相飽和度(即液態(tài)水占據(jù)的孔隙體積比例)高達(dá)17%，并且范圍較大，這種狀態(tài)必然有利于膜的潤濕而阻礙氧氣的傳輸.隨著溫度的升高，液相飽和度逐漸降低，80 ℃時(shí)僅局部有少量液態(tài)水，90 ℃時(shí)所有反應(yīng)面積上都為干的.

圖4 300 kPa下不同溫度催化層中間平面相飽和度

圖5為300 kPa下膜中中間平面上磺酸基團(tuán)水含量分布.其總體分布和液相飽和度的分布是一致的，隨著溫度的升高，水含量降低，即膜的潤濕性降低，質(zhì)子電阻升高.

圖5 300 kPa下不同溫度膜中間平面水含量

圖6 St=2.5，θ=70 ℃，隨電流密度的變化

圖6是過量系數(shù)為2.5，溫度為70 ℃時(shí)Δq隨著操作壓力的變化關(guān)系，Δq定義為反應(yīng)生產(chǎn)水與陰極出口氣體以飽和狀態(tài)帶走的水的差值.由圖6可知，Δq隨著壓力升高而增加，說明升高操作壓力可以使較多的水保留在MEA中，有益于膜的潤濕；降低操作壓力，被空氣帶走的水增加，MEA剩余的水分減少，因此膜容易變干，為保持膜較佳的潤濕狀態(tài)進(jìn)而保持膜的性能，需要操作降低溫度.在操作壓力為200 kPa和150 kPa時(shí)，電池性能最優(yōu)時(shí)對(duì)應(yīng)的溫度分別為60 ℃和55 ℃，這些點(diǎn)構(gòu)成了電池性能最優(yōu)的溫度控制曲線，見圖3.

圖7是操作壓力為100 kPa下不同溫度膜中間平面水含量.同樣的電流密度與溫度下，與圖5相比較，膜的潤濕狀態(tài)顯著下降，反應(yīng)在圖3中，是電池性能顯著下降.圖3中在100 kPa等操作壓力下電池性能并非單調(diào)變化，溫度增高時(shí)陰極側(cè)膜潤濕性降低，但陽極側(cè)保持100%加濕，實(shí)際帶入的水汽量增加，膜陽極側(cè)以及膜整體的潤濕有所改善；即使膜中水含量一定，隨著溫度的升高，水在膜內(nèi)擴(kuò)散系數(shù)會(huì)增加，水在膜內(nèi)的分布會(huì)更加均勻，因此性能反而有增高趨勢(shì)，但總體上增量不大.

圖7 100 kPa時(shí)不同溫度膜中間面的水含量

人們除了希望PEM電池能有最佳的電輸出性能外，同時(shí)也希望在活性面積上電流密度等各物理場(chǎng)能均勻分布，這樣才可以使得電池長久、穩(wěn)定運(yùn)行.圖7和圖5對(duì)比，同樣的溫度下，壓力較高時(shí)膜中水含量分布較為均勻；同樣的壓力下，在溫度較低時(shí)因催化層液態(tài)水較多，使得膜中水含量較高且分布較為均勻.圖8為50 ℃與75 ℃時(shí)不同壓力膜中間面的電流密度分布.其中橫坐標(biāo)為活性面積上出現(xiàn)的電流密度范圍，從0.9到2.1 A/cm2，以0.1 A/cm2分檔計(jì)算.在50 ℃、 300 kPa下，大部分區(qū)域工作電流密度為1.5～1.8 A/cm2，而隨著操作壓力的降低，工作電流密度逐漸分散在更大區(qū)域內(nèi)，也即分布更加不均勻，這和前面分析的膜的水含量分布規(guī)律是一致的.在75 ℃ 100 kPa下，大部分區(qū)域工作電流密度為1.5～1.8 A/cm2，顯得比較均勻，這是因?yàn)榇藭r(shí)膜偏干，膜干的“比較均勻”，而隨著操作壓力的增加，膜總體的潤濕性提高，但并不均勻，使得工作電流密度逐漸在更大區(qū)域內(nèi)分散，也即分布更加不均勻，這種分布狀態(tài)不利于對(duì)電池長久、溫度運(yùn)行，因此還需要在結(jié)構(gòu)設(shè)計(jì)和運(yùn)行管理方面進(jìn)行改進(jìn).

圖8 50 ℃與75 ℃時(shí)不同壓力膜中間面的電流密度分布

3 結(jié) 論

1) 在1.5 A/cm2、陽極100%加濕、陰極無加濕條件下，隨著操作壓力的增加，電池性能增加.

2) 對(duì)于給定的操作壓力，存在一個(gè)最佳性能溫度，壓力增加，最佳溫度也隨之增加，在壓力為100，200，300 kPa時(shí)，最佳性能溫度分別為75，60，50 ℃.

3) 由于陰極不加濕，在最佳溫度時(shí)膜中的含水量分布及電流密度分布不均勻；在較低溫度下運(yùn)行時(shí)電流密度分布均勻性隨著壓力的增加逐漸提高，但在較高溫度下運(yùn)行時(shí)則相反.

[1]YUN W, KEN S, CHEN Y, et al. A review of polymer electrolyte membrane fuel cells: technology, applications, and needs on fundamental research[J]. Applied Energy,2011,88:981-1007.

[2]DONG H J, KWANG N K, SEUNG M B, et al. The effect of relative humidity of the cathode on the performance and the uniformity of PEM fuel cells[J]. International Journal of Hydrogen Energy,2011,36:12499-12511.

[3]SU A, FERNG Y M, SHIH J C. CFD investigating the effects of different operating conditions on the performance and the characteristics of a high-temperature PEMFC[J]. Energy,2010,35:16-27.

[4]SUKKEE U, WANG C Y. Computational study of water transport in proton exchange membrane fuel cells[J]. Journal of Power Sources,2006,156:211-223.

[5]KIM H Y, KIM K. Numerical study on the effects of gas humidity on proton-exchange membrane fuel cell performance[J]. Hydrogen Energy,2016(1)：1-8.

[6]AMIRINEJAD M, ROWSHANZAMIR S, MOHAMMAD H, et al. Effects of operating parameters on performance of a proton exchange membrane fuel cell[J]. Journal of Power Sources,2006,161:872-875.

[7]SANTARELLI M G, TORCHIO M F. Experimental analysis of the effects of the operating variables on the performance of a single PEMFC[J]. Energy Conversion and Management,2007,28:40-51.

[8]JANG J H, CHIU H C, YAN W M, et al. Effects of operating conditions on the performances of individual cell and stack of PEM fuel cell[J]. Journal of Power Sources,2008,180:476-483.

[9]CORINNA H, FRANK K, ALEXANDER D. Study of the influence of key test parameters on the performance of a PEMFC stack[J]. Solid State Ionics,2015,275:75-79.

[10]HOMAYOON K, MEHRZAD S, MOHAMMADREZA H, et al. Model development and optimization of operating conditions to maximize PEMFC performance by response surface methodology[J]. Energy Conversion and Management,2015,93:9-22.

[11]HUANG C P, JIANG R C, ELBACCOUCH M, et al. On-board removal of CO and other impurities in hydrogen for PEM fuel cell applications[J]. J Power Sources,2006,162:563-71.

[12]MIN C H, HE Y L, LIU X L, et al. Parameter sensitivity examination and discussion of PEM fuel cell simulation model validation part II: results of sensitivity analysis and validation of the model[J]. Journal of Power Sources,2006,160:374-385.

Optimization of Operating Conditions for PEMFC Running at Large Current Density

ZHANG Hongkai1,2)ZHAN Zhigang1,2)HE Xiaobo1,2)SHUAI Lu2)SUI Bangjie3,4)PAN Mu1)

(State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan 430070, China)1)(School of Power and Energy Engineering, Wuhan University of Technology, Wuhan 430063, China)2)(School of Automotive Engineering, Wuhan University of Technology, Wuhan 430070, China)3)(Institute for Integrated Energy Systems and Department of Mechanical Engineering, University of Victoria, Victoria V8W 2Y2, Canada)4)

Considering the performance and uniformity of the current density distribution, CFD software is used in this paper to simulate the PEM cells, which operate at large current density, no cathode humidification but full anode humidification, and backpressure, to seek for the optimized operating point. The main conclusions are as follows: for a given operating pressure, there is an optimal temperature which makes the fuel cell have the best performance, and the optimal temperature increases with the pressure; the best temperatures for 100,200 and 300 kPa are 75,60 and 50 ℃, respectively. Since the cathode is not humidified, the membrane water content distribution and current density distribution are not uniform at the optimum performance temperature point; the current density distribution uniformity increases with pressure at lower temperatures, but it has the opposite trend at higher temperatures.

PEMFC; current density; water transport; operating conditions; optimization

2017-02-23

*國家自然科學(xué)基金項(xiàng)目(21676207)、國家重點(diǎn)研發(fā)計(jì)劃“新能源汽車”試點(diǎn)專項(xiàng)課題(2016YFB0101207)資助

TM911.4

10.3963/j.issn.2095-3844.2017.03.033

張洪凱(1993—)：男，碩士生，主要研究領(lǐng)域?yàn)槿剂想姵貍鬏敩F(xiàn)象和新能源動(dòng)力裝置及系統(tǒng)