质子交换膜燃料电池动态建模及特性分析.pdf
《质子交换膜燃料电池动态建模及特性分析.pdf》由会员分享,可在线阅读,更多相关《质子交换膜燃料电池动态建模及特性分析.pdf(6页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、书书书第 17 卷第 4 期电化学Vol 17No 42011 年 11 月ELECTROCHEMISTRYNov 2011文章编号:1006-3471(2011)04-0438-06质子交换膜燃料电池动态建模及特性分析贾秋红1,2,韩明3*,邓斌1,张财志3(1 西南交通大学新型驱动技术中心,四川 成都 60031;2 重庆理工大学重庆汽车学院,重庆 400050;3 淡马锡理工学院工程学院,新加坡 529757)收稿日期:2011-07-09,修订日期:2011-08-17*通讯作者,Tel:(86-65)68705652,E-mail:minghantp edu sg重庆市科委自然科学基
2、金计划项目(CSTC,2010BB4302)和重庆市教育委员会科学技术研究项目(KJ110810)资助摘要:根据分析质子交换膜燃料电池电压的影响因素,利用理想气体状态方程的微分形式求得氢气和氧气分压,并且结合热力学电动势及燃料电池的各种损耗,建立 PEMFC 的动态模型 合理设定相关参数,得到 PEM-FC 的稳态模型,其仿真程度与实验结果有较好的吻合,表明该模型的正确性、可操作性和有效性 利用建立的动态模型分析在载荷突变的情况下,输出电压和氢气分压的动态响应特性,从而达到准确模拟实际燃料电池系统的工作性能,评价预测燃料电池的承载能力和指导燃料电池堆操作的目的关键词:质子交换膜燃料电池;动态模
3、型;氢气分压;载荷突变;动态响应特性中图分类号:TM911文献标识码:A燃料电池是新世纪最有前途的清洁能源之一,有可能取代热机成为主要的能量转换方式 其中质子交换膜燃料电池(Proton Exchange Mem-brane Fuel Cell,PEMFC)具有功率密度高、能量转换效率高、可低温快速启动和环境友好等优点,成为了目前世界上发展最为迅速的燃料电池1-2 到目前为止,国内外已就 PEMFC 特性作了大量的研究探索,并取得了重要的研究进展,其中大多是依据稳态数学模型 相比之下,有关动态模型的研究尚存在不足3-4 而 PEMFC 动态特性研究对于了解其运行机理、优化电池结构以及系统控制非
4、常重要 因此,近年来,燃料电池动态特性的研究已成为该领域的一个热点课题 Khan4 等发展了一个同时考虑电池电化学瞬态过程、传热传质和反应气体流动过程的 PEMFC 非稳态模型;Hernan-dez5 等应用非线性状态空间模拟方法研究了流道内的液态水对 PEMFC 动态特性的影响;Buch-holz6 等用非线性系统识别方法对 PEM 电池堆的动态特性进行了研究;Francol7 等发展了 PEMFC瞬态分析的多尺度动态机理模型;莫志军8 等利用“黑箱”模型就质子交换膜燃料电池进行了建模与动态仿真研究;胡桂林9 等发展了一个依据计算流体动力学的非稳态、非等温的三维两相流数学模型 以上动态模型的
5、建立主要是采用分布参数模型,仿真精确度高、能够模拟电池内部状态的空间分布,但仿真占用时间多,收敛速度慢 本文在 PEMFC 经验模型的基础上,通过理想气体状态方程推导求得氢气和氧气分压的计算公式,建立了一种 PEMFC 动态特性建模和仿真的简单实用方法,利用 Matlab/Simulink 工具对电池堆在载荷突变情形下电池电压、反应气体压力的动态响应特性模拟实际燃料电池系统的工作性能 研究结果表明,该方法可以应用于燃料电池的动态特性仿真分析、优化设计和燃料电池系统的自动控制1PEMFC 工作原理及电化学模型1 1PEMFC 的工作原理一个典型的 PEMFC 由集流板、流场板、气体扩散层、催化层
6、和质子交换膜等组成 PEMFC 以全氟磺酸型的固体聚合物为电解质,铂/碳或铂-铑/碳为电催化剂,氢或净化重整气为燃料,空气或纯氧为氧化剂,带有气体流动通道的石墨板或表面改性的金属板为双极板 电池电极是一种多孔气体扩散电极,一般由扩散层和催化层组成,扩散层一般由碳纸或碳布制作,其作用是支撑催化层、收集电流,并为电化学反应提供电子通道、气体通道和排水通道 催化层则是发生电化学反应的场所,是电极的核心部分当分别向阳极和阴极供给氢气与氧气时,反应气体经扩散层扩散,进入多孔阳极的氢原子被催化剂吸附并分离为氢离子和电子,氢离子经由质子交换膜转移到阴极(正极),电子在电极内传递至负极(阳极)集流板经外电路负
7、载流向阴极,阴极催化层氢离子和氧原子结合成水分子,生成的水通过电极随反应尾气排出电极反应如下:阳极:H22H+2e阴极:1/2O2+2e+2H+H2O总反应:1/2O2+H2H2O(1)1 2 电池的电化学模型PEM 燃料电池单电池在反应生成液态水的情况下其理论标准电势为1 229 V,但由于不可逆损失,实际电池电势随平衡电势的降低而下降 燃料电池的不可逆损失常被称为极化过电势或者过电压 极化的产生主要有3 种原因:活化极化;欧姆极化;浓差极化,从而导致燃料电池实际电压小于理论电势质子交换膜燃料电池单电池的输出电压Vcell10 可用式(2)表示:Vcell=ENernst act ohm c
8、on(2)式中,ENernst是热力学电动势;act是活化过电压;ohm是欧姆过电压;con是浓差过电压由 N 个燃料电池单体串联在一起的燃料电池堆,其电压 Vst可用式(3)表示:Vst=NVcell(3)1)热力学电动势根据 Nernst 方程由氢气、氧气组成的燃料电池,热力学电动势为3:ENernst=G2F+S2F(T Tref)+RT2F ln(PH2)+12ln(PO2(4)式中,G 表示电池反应 Gibbs 自由能的变化(J/mol);F 是法拉第常数(96487C);S 是电池反应的熵变化(J/mol);R 是气体常数(8 314 J/Kmol);PH2是氢气在阳极催化剂/气体
9、界面的分压(atm);PO2是氧气在阴极催化剂/气体界面的分压(atm);T 是电池温度(K);Tref表示参考温度(K);如以标准大气压和常温下的 G、S 和 Tref代入上式,得到简化的表达式为:ENernst=1 229 8 5 104(T 298 15)+4 308 105 T (ln(PH2)+12ln(PO2)(5)2)活化过电压由电池阴极和阳极产生的活化过电压 act11 表示如下:act=0 9514+0 0032 T 1 87 104T ln(i)+7 4 105 T ln(CO2)(6)式中,CO2为气液界面溶解氧的浓度,由 Henry定律可知:CO2=PO2/5 08 1
10、06exp(498/T)(7)3)浓差极化过电压浓差过电压决定反应物传递到反应区以及产物的排出速率,有限的传递速率导致电化学反应区氧气和氢气浓度低于主流区的浓度 文献 12-13 将浓差过电压 con表示为:con=meni(8)式中,m 是与温度有关的质量传递系数,n 是催化层中电化学反应产物的增长率 m 和 n 在通常情况下约为 3 105V 和 8 103cm2mA114 4)欧姆极化过电压欧姆电阻主要是由于质子交换膜的膜电阻 Rm和阻碍电子通过膜的阻抗 Rc引起的,可由式(9)表示:ohm=iRohm=i(Rm+Rc)(9)由于离子电荷传输比电子电荷传输更难,其在欧姆电阻中占主要作用1
11、4,因此不妨忽略电子通过膜的阻抗作用根据文献 15,质子交换膜的阻抗 Rm可用式(10)表示:Rm=rMA(10)式中,是质子交换膜的厚度(m);A 是膜的有效面积(cm2)就 Nafion 系列,其质子交换膜的电阻率 rM可用以下的经验式(11)表示16:rM=181 6 1+0 03(iA)+0 062(T303)2(iA)2 5 0 634 3(iA)exp(4 18T 303T)(11)式中,表示质子交换膜的含水量,是一个可调整的参数,乃为气体相对湿度和化学计量数的函数,由电池实际工作状态决定,其值一般在 14 934第 4 期贾秋红等:质子交换膜燃料电池动态建模及特性分析23 之间1
12、7,本文取 202PEMFC 的动态模型建立燃料电池是一种在稳态下能够提供可靠功率的能源,但它不能为快速的瞬态电流载荷作出及时响应 其中主要的原因即在燃料电池内部的电化学反应和热动力学与机械附件系统未能达到同步匹配 为了能够准确评估和预测燃料电池在不同瞬态条件下(燃料电池终端的电气故障、燃料电池车的启动,加速和制动等)的状况,建立燃料电池动态模型的仿真是非常必要的18 2 1假设条件根据实际的实验操作条件,本模型作如下假设:1)燃料电池各气体均为理想气体;2)电池气体流道内气体压力处处相等;3)电池内的温度处处相等,水为液态水;4)各单电池的性能参数相同,电池堆性能参数是各单电池性能参数的累加
13、2 2气体摩尔流量1)气体反应的摩尔流量根据法拉第定律和式(1)、(2)和(3),在给定电流 I 的条件下,由 N 个单电池组成的电池堆参加电化学反应消耗掉的氢和氧摩尔流量(qH2,r,qO2,r)可分别用下式表示:qH2,r=I2FNqO2,r=I4FN(12)2)输入气体的摩尔流量qH2,in=qH2,r/Uf(13)式中,Uf是氢的利用率,一般情况下,氢的利用率为 80%90%19,本文取 85%2 3气体分压依照假设条件,各气体均应满足理想气体状态方程,对氢气,有方程:PH2VH2=nH2RT(14)式中 PH2为氢气的气体分压(atm);VH2和 nH2分别表示氢气在电池堆中所占的体
14、积及其摩尔数;R 为通用气体常数(8 314 J/molK)以式(14)对时间求导数得:ddtPH2=RTVH2ddtnH2因ddtnH2=qH2,in qH2,out qH2,r即ddtPH2=RTVH2(qH2,in qH2,out qH2,r)(15)此处,qH2,in和 qH2,out分别表示入口处和出口处氢的摩尔流量按文献 20,出口处氢的摩尔流量可表示为qH2,out=KH2 PH2(16)式中,KH2为氢气的阀摩尔值常数 kmol/(satm)所以ddtPH2=RTVH2(qH2,in KH2 PH2qH2,r)将上式两边作拉普拉斯变换并结合式(12)和(13)得:PH2=RTR
15、TKH2+VH2 s(qH2,in qH2,r)=1/KH21+VH2RTKH2 s(1Uf11+es)12FN(17)同理可得:PO2=1/KO21+VO2RTKO2 s(2 11Uf11+es)14FN(18)2 4PEMFC 的动态模型综合式(2),(5)(11)和(17)(18)可得到PEMFC 的动态模型 本文使用 Matlab/Simulink 进行模型的搭建及子系统封装,其建立的 PEMFC 仿真模块框图如图 1 所示3PEMFC 模型实验验证3 1仿真参数的设定为了能够实现仿真模型的运行,一些必要的物理参数须有确定的值,本模型参数取值均取自实际的实验条件 表(1)列出 1 2
16、kW PEMFC 的动态模型仿真参数3 2稳态电压电流特性分析及实验验证假设 e、H2、O2等为零系统达到稳态,然后在Matlab/Simulink 环境下计算稳态模型 V I 特性曲线,并与实验结果相比较,验证模型的准确性 对1 2 kW PEMFC,实验测试的参数设定为:温度为323 15 K,阳 极 和 阴 极 气 体 的 进 口 压 力 分 别为2 atm(203 kPa)和3 5 atm(355 kPa),相对湿度044电化学2011 年图 1PEMFC 的 Simulink 仿真模块框图Fig 1Simulation module diagram of Simulink for t
17、he PEMFC表 11 2 kW PEMFC 动态模型仿真参数Tab 1Simulation parameters of dynamic model forthe 1 2 kW PEMFCParameterValueParameterValueT323 15KH_O1 145N40A97 75F96485C/molH226 1sKO22 52mol/(satm)O22 91sKH20 843mol/(satm)e0 8 s阳极为100%,阴极为90%由 G200 型燃料电池测试设备(Greenlight Innovation 公司)测试 图 2 给出PEMFC 稳态模型仿真与实验结果比较可以
18、看出,在电流密度较低的情况下,实验结果与仿真结果有一定的偏差,主要是由于模型建立时简化一些公式引起的,但其偏差程度多在误差范围内 其它情况,实验结果与仿真计算结果吻合较好,证明了该模型建立的准确性,继续进一步的完善研究,当可指导燃料电池的操作4动态特性分析4 1动态模型的考虑本文建立的动态模型主要考虑了氢气、氧气的动态响应时间项以及电化学反应时间的动态响应常数项在载荷突变时对燃料电池动态特性的图 2PEMFC 稳态仿真与实验的 V I 特性曲线Fig 2A comparison in the V I characteristic curves ofthe PEMFC between simul
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 质子 交换 燃料电池 动态 建模 特性 分析
限制150内