耦合催化层结块的质子交换膜燃料电池建模方法与流程

1.本发明属于电化学能量转换及燃料电池领域，具体涉及一种耦合催化层结块的质子交换膜燃料电池建模方法。


背景技术：

2.氢能是可再生能源的能量转化中枢，随着各国碳中和政策的陆续颁布，人类社会开始逐步向无碳化氢能社会探索。质子交换膜燃料电池(pemfc)作为氢能网络中的重要终端应用装置，可将氢能直接转化为电能，具有清洁无污染，无噪音，功率密度高，转化效率高，启动速度快以及运行温度低等优点，被广泛认为是下一代能量转化装置，可应用于交通，便携式能源，固定电站以及分布式能源系统等领域。然而，质子交换膜燃料电池还处于商业化初期，如何实现高性能、长寿命以及低成本的电池是目前的研究重点。
3.质子交换膜燃料电池包含有复杂的多相反应传质过程，如气体扩散，电子/质子传递，气液相变，膜吸/放水以及液体渗流等。反应气体、电子和质子分别在不同通道内传输，共同在三相界面处发生电化学反应。催化层中贵金属催化剂的使用是造成燃料电池成本高的主要原因之一，减少贵金属催化剂的用量是目前的研究热点。然而，催化剂的减少造成反应面积降低，增加了反应物的传输阻力，因此不可避免的恶化了电池的性能。理解催化层中复杂的多相反应传质过程可以有效指导催化层的设计优化，降低成本，提高电池的性能和寿命。
4.构建高精度、高鲁棒性质子交换膜燃料电池模型可以有效减少实验成本，获得电池内多相传输反应传质机制，从而指导燃料电池的设计与优化。如前所述，催化层中包含有复杂的多相反应传质过程，对其进行合理的模型描述存在明显挑战。常用的界面模型和均质模型未能考虑催化层组分间的相互耦合关系以及氧气跨膜传输的巨大阻力，常常过高的预估电池性能。而介观和微观分子动力学模型由于时间/空间尺度的限制很难反映电池的宏观性能。
5.在背景技术部分中公开的上述信息仅仅用于增强对本发明背景的理解，因此可能包含不构成在本国中本领域普通技术人员公知的现有技术的信息。


技术实现要素：

6.针对现有技术中存在的问题，本发明提供一种耦合催化层结块的质子交换膜燃料电池一维建模及求解方法。综合考虑电池内部的气液两相流动，扩散以及相变过程，耦合所建立的催化层结块模型，引入氧气跨膜阻力和修正电化学模型，改善燃料电池模型中未考虑局部氧气传输的不足，可考虑催化层的构成对氧气传输及宏观性能的影响。并相应地提出一种基于该模型的高精度、高鲁棒性的迭代求解方法及流程。
7.本发明的目的是通过以下技术方案予以实现，耦合催化层结块的质子交换膜燃料电池建模方法包括以下步骤：
8.第一步骤，构建的质子交换膜燃料电池包括阳极流道和阴极流道和膜电极，所述
膜电极包括阳极和阴极侧的气体扩散层gdl、微孔层mpl，催化层cl以及质子交换膜pem，质子交换膜燃料电池中的各组分包括括氢气、氧气、氮气以及水蒸气，构建各组分的组分扩散方程：
[0009][0010]
其中t为时间；为哈密顿算子；ε和s分别为多孔电极的孔隙率和液态水饱和度；ci为各组分摩尔浓度；为各组分的有效扩散系数；mi为各组分的摩尔质量；si表示由电化学反应或相变造成的组分源项，其中，
[0011][0012][0013][0014][0015]
式中ja和jc分别为阴阳极两侧的电化学反应速率；f为法拉利常数；sv为水蒸气组分扩散方程的源项；s
v-l
为气态水转化为液态水的源项；s
m-v
为膜态水转化为气态水的源项；i为h2、或o2、或n2、或水蒸气(v)；agdl和cgdl分别代表阳极和阴极侧的气体扩散层；ampl和cmpl分别代表阳极和阴极侧的微孔层；acl和ccl分别代表阳极和阴极侧的催化层；gdls和mpls以及cls分别为阴阳极两侧的gdl，mpl和cl；
[0016]
第二步骤，构建质子交换膜燃料电池中多孔电极的液压守恒方程，其中，式中ρ1，p1，μ1，k1分别代表液态水的密度，压力，粘度和相对渗透率；k为多孔介质的固有渗透率；s1为方程的源项，其为：
[0017][0018]
其中，表示h2o的摩尔质量；第三步骤，构建质子交换膜燃料电池中膜态水守恒方程：
[0019]
[0020]
式中ρ
pem
和ew分别表示膜的密度和当量重量；ω，λ和表示离聚物的体积分数，膜态水含量和有效扩散系数；s
mw
为膜态水方程的源项，其为：
[0021][0022]
式中，s
eod
和s
p
分别代表由于电渗拖曳和液压渗透作用产生的源项；
[0023]
第四步骤，构建质子交换膜燃料电池中的能量守恒方程，
[0024][0025]
式中c
p，l
，c
p，g
和c
p，s
分别表示液体，气体和固体的比热容；k
eff
代表有效导热系数；ρg和ρs分别为气体和固体的密度；t为温度；s
t
为能量方程中的源项，分别由不同组件的反应热，欧姆热，相变热以及活化热组成，
[0026][0027]
式中j为垂直于膜方向上的体积电流密度；δsa和δsc分别为阴极和阳极两侧的熵差；η
ohm
和η
act
分别为欧姆和活化过电势；h为相变潜热；bps表示阴阳极两侧的双极板；
[0028]
第五步骤，构建电化学模型，在恒电流计算模式下，电池的输出电压为可逆电压减去活化、欧姆和浓差过电势：v＝e
rev-η
ohm-η
act，con
，
[0029][0030][0031]
其中e
rev
为可逆电压；η
act，con
为活化和浓差过电势；和为进口氢气和氧气的气压；和分别为电子通过极板和多孔电极骨架的面积比电阻；和分别代表质子通过催化层和膜中的面积比电阻；r为理想气体常数；i为运行的电流密度，
[0032]
阳极侧的活化和浓差过电势为：
[0033][0034][0035]
其中j
0，a
和分别为阳极体积交换电流密度和参考体积交换电流密度；为氢气浓度；为参考氢气浓度；αa和αc分别为阳极和阴极的转换系数；ηa为阳极侧的活化和浓差过电势；
[0036]
第六步骤，构建阴极催化层结块模型，氧气从催化层孔隙中到pt催化剂表面间存在氧气阻力，氧气阻力为：
[0037][0038]
其中r
l
为氧气穿过液态水薄膜的扩散阻力；r
ion，int
为氧气进入离聚物表面的界面阻力；为氧气在离聚物中的扩散阻力；为氧气从离聚物中到pt颗粒表面的界面阻力；
[0039]
第七步骤，建立所述催化层结块模型的耦合关系，pt表面参加电化学反应的氧气浓度为：
[0040][0041]
其中和分别表示pt表面和离聚物界面上的氧气浓度；i
local
为离聚物表面的局部电流密度，a
ion
为离聚物的比表面积；
[0042][0043]
式中为亨利系数；为催化层孔隙中的氧气浓度，
[0044][0045][0046]
其中i
0，c
和为交换电流密度和参考交换电流密度；ηc为阴极侧的活化和浓差过电势；w为能量参数；a
pt
为单位体积pt催化剂的表面积；θ
pt/o
为pt氧化的覆盖率；
[0047]
第八步骤：构建组分扩散方程、液压守恒方程、膜态水守恒方程和能量守恒方程的边界条件，其中组分扩散方程和液态守恒方程的两侧边界分别为狄利克雷边界和纽曼边界，分别如下：
[0048][0049][0050][0051][0052]
式中agc/agdl为阳极侧流道与扩散层的交接面；acl/pem为阳极侧催化层与质子交换膜的交接面；cgc/cgdl为阴极侧流道与扩散层的交接面；ccl/pem为阴极侧催化层与质子交换膜的交接面；c
i，in
和c
i，out
为组分的进口和出口浓度；p
l
和pg代表液体和气体压力；ci为前文所述的各组分摩尔浓度；c
i，agc/agdl
为阳极侧流道与扩散层的交接面上各组分摩尔浓度；c
i，cgc/cgdl
为阴极侧流道与扩散层的交接面上各组分摩尔浓度，
[0053]
膜态水守恒方程和能量守恒方程的边界条件各自采用纽曼边界和狄利克雷边界条件：
[0054][0055]
t
agc/agdl
＝t0，t
cgc/cgdl
＝t0，
[0056]
其中acl/ampl代表阳极侧催化层与微孔层的交界面；ccl/cmpl代表阴极侧催化层与微孔层的交界面；t
agc/agdl
为阳极侧流道与扩散层的交接面的温度；t
cgc/cgdl
为阴极侧流道与扩散层的交接面的温度；t0为燃料电池的运行温度；λ为前文的膜态水含量。
[0057]
所述的一种耦合催化层结块的质子交换膜燃料电池建模及求解方法中，采用有限容积法离散所述组分扩散方程、液压守恒方程、膜态水守恒方程和能量守恒方程，其中扩散项基于二阶中心差分进行离散，纽曼边界条件采用附加源项法进行处理。
[0058]
所述的一种耦合催化层结块的质子交换膜燃料电池建模及求解方法中，引入gauss-seidel方法对所述组分扩散方程、液压守恒方程、膜态水守恒方程和能量守恒方程迭代求解，其中，
[0059]
步骤1，设置电流密度，输入运行和结构参数，对各物理量场进行初始化；
[0060]
步骤2，基于物理量更新模型中所有需要的物性和源项；
[0061]
步骤3，依次求解各组分扩散方程；
[0062]
步骤4，求解液压守恒方程并内迭代获得多孔电极内的液态水饱和度；
[0063]
步骤5，求解膜态水守恒方程；
[0064]
步骤6，求解相互耦合的电化学模型和催化层结块模型；
[0065]
步骤7，判断电流密度是否到达极限电流密度，若达到则结束求解，若未达到则继续接下来的流程；
[0066]
步骤8，求解能量守恒方程；
[0067]
步骤9，判断所述组分扩散方程、液压守恒方程、膜态水守恒方程和能量守恒方程
是否达到残差标准，若达到则更新电流密度，以上一电流密度计算结果作为初始值开始下一电流密度的计算；若未满足收敛准则，对物理量进行超松弛或亚松驰更新，返回步骤2继续进行计算。
[0068]
所述的一种耦合催化层结块的质子交换膜燃料电池建模方法中，第六步骤中，氧气穿过液态水薄膜的扩散阻力r为：
[0069][0070]
其中为氧气在液态水中的扩散系数；为液态水的当量厚度：
[0071][0072]aion
＝4π(rc δ
ion
)2n
pt/c
，
[0073]
其中a
ion
为离聚物的比表面积；rc为碳载体的半径；δ
ion
和n
pt/c
分别为离聚物膜的厚度以及单位体积pt/c的数目，分别为：
[0074][0075][0076]
其中ε
ion
和ε
pt/c
分别为离聚物和pt/c的体积分数，
[0077][0078][0079]
式中，m
pt
，ξ
pt/c
和ζ
i/c
代表铂载量，pt/c质量比和i/c质量比；ρ
pt
和ρc分别为铂和碳的密度；εc为碳载体的体积分数；催化层的孔隙率由离聚物和pt/c的体积分数获得：
[0080]
ε
p
＝1-ε
ion-ε
pt/c
，
[0081]
氧气进入离聚物表面的界面阻力r
ion，int
通过以下关系获得：
[0082][0083]
式中k1为界面阻力系数；为氧气在离聚物中扩散的扩散系数。
[0084]
所述的一种耦合催化层结块的质子交换膜燃料电池建模方法中，氧气在离聚物中的扩散阻力为：
[0085]
[0086]
氧气从离聚物中到pt颗粒表面的界面阻力为：
[0087][0088]
式中，θ
pt/o
为pt氧化的覆盖率；k2表示pt表面吸附阻力系数；r
pt
表示pt颗粒的半径，a
ecsa
为电化学反应面积。
[0089]
本发明的有益效果在于：
[0090]
本发明构建了一维多相非等温质子交换膜燃料电池模型，综合考虑质子交换膜燃料电池内部的气液两相流动，扩散以及多种状态水间的相变及守恒过程。耦合所建立的催化层结块模型，引入氧气跨膜阻力和修正电化学模型，可考虑催化层的构成对氧气传输，电荷传输及宏观性能的影响，克服了高电流密度下传统模型不能反映极大浓差损失的问题。并基于该模型，相应地提出了一种高精度、高鲁棒性的迭代求解方法及流程。本发明可为质子交换膜燃料电池的精准预测和设计优化提供指导。
[0091]
上述说明仅是本发明技术方案的概述，为了能够使得本发明的技术手段更加清楚明白，达到本领域技术人员可依照说明书的内容予以实施的程度，并且为了能够让本发明的上述和其它目的、特征和优点能够更明显易懂，下面以本发明的具体实施方式进行举例说明。
附图说明
[0092]
通过阅读下文优选的具体实施方式中的详细描述，本发明各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。说明书附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。显而易见地，下面描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图；而且在整个附图中，用相同的附图标记表示相同的部件。
[0093]
在附图中：
[0094]
图1为质子交换膜燃料电池示意图；
[0095]
图2为催化层结块模型示意图；
[0096]
图3为模型迭代求解方法及流程示意图；
[0097]
图4为本发明一个实施例的模型验证极化曲线；
[0098]
图5为本发明一个实施例的不同pt载量下的极化曲线；
[0099]
图6为本发明一个实施例的氧气传输阻力随pt载量的变化曲线；
[0100]
图7为本发明一个实施例的pt催化剂表面氧气浓度随pt载量的变化曲线。
[0101]
以下结合附图和实施例对本发明作进一步的解释。
具体实施方式
[0102]
下面将参照附图1至图7更详细地描述本发明的具体实施例。虽然附图中显示了本发明的具体实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施
例所限制。相反，提供这些实施例是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
[0103]
需要说明的是，在说明书及权利要求当中使用了某些词汇来指称特定组件。本领域技术人员应可以理解，技术人员可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名词的差异来作为区分组件的方式，而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”或“包括”为一开放式用语，故应解释成“包含但不限定于”。说明书后续描述为实施本发明的较佳实施方式，然所述描述乃以说明书的一般原则为目的，并非用以限定本发明的范围。本发明的保护范围当视所附权利要求所界定者为准。
[0104]
为便于对本发明实施例的理解，下面将结合附图以具体实施例为例做进一步的解释说明，且各个附图并不构成对本发明实施例的限定。
[0105]
耦合催化层结块的质子交换膜燃料电池建模方法包括如下步骤：
[0106]
第一步骤：质子交换膜燃料电池一般由极板，流道和膜电极组成，其中膜电极包含有气体扩散层(gdl)，微孔层(mpl)，催化层(cl)以及质子交换膜(pem)。催化层是发生电化学反应的场所，主要由孔隙、铂(pt)族催化剂、碳载体以及离聚物等部分构成。构建质子交换膜燃料电池中各组分的组分扩散方程，包括氢气，氧气，氮气以及水蒸气。多孔电极中的对流作用较弱，为求解方便忽略方程中的对流项，在阴阳极两侧的gdl，mpl和cl内进行求解：
[0107][0108]
其中t(s)为时间；为哈密顿算子；ε和s分别为多孔电极的孔隙率和液态水饱和度；ci(mol m-3
)为各组分摩尔浓度；d
ieff
(m
2 s-1
)为各组分的有效扩散系数，采用经典bruggeman修正公式，基于孔隙率和液态水饱和度对组分本身的扩散系数进行修正；mi(kg mol-1
)为各组分的摩尔质量；si(kg m-3
s-1
)表示由电化学反应或相变造成的组分源项，相应的具体表达式为：
[0109][0110][0111][0112][0113]
式中ja和jc(am-3
)分别为阴阳极两侧的电化学反应速率；f为法拉利常数，取96487c mol-1
；sv为水蒸气组分扩散方程的源项；s
v-1
为气态水转化为液态水的源项；s
m-v
为膜态水转化为气态水的源项；i为h2、或o2、或n2、或水蒸气(v)；agdl和cgdl分别代表阳极和
阴极侧的气体扩散层；ampl和cmpl分别代表阳极和阴极侧的微孔层；acl和ccl分别代表阳极和阴极侧的催化层；gdls和mpls以及cls分别为阴阳极两侧的gdl，mpl和cl
[0114]
作为优选，s
v-1
和s
m-v
可由下列模型求得：
[0115][0116][0117]
其中，表示h2o的摩尔质量；
[0118]
作为优选，传质系数γ
v-1
和γ
m-v
(s-1
)分别选择为100s-1
和1.3s-1
。水蒸气的饱和浓度c
sat
(mol m-3
)和平衡状态下的膜态水含量λ
eq
可分别由经典经验关联式获得：
[0119][0120][0121][0122]
式中p
sat
(pa)为水蒸气的饱和压力；a为水蒸气的活度；
[0123]
第二步骤：由于孔隙参数的不同，不同多孔电极交接面上液压连续，液态水饱和度并不连续。因此，基于液态水质量守恒和darcy定律，构建质子交换膜燃料电池中多孔电极的液压守恒方程，均在阴阳极两侧的gdl，mpl和cl内进行求解：
[0124][0125]
式中ρ1(kg m-3
)，p1(pa)，μ1(pa s)，k1分别代表液态水的密度，压力，粘度和相对渗透率；k为多孔介质的固有渗透率；s1(kg m-3
s-1
)为方程的源项，如下：
[0126][0127]
作为优选，相对渗透率可通过如下公式进行计算：
[0128]
k＝s
3.0
ꢀꢀ
(2-3)
[0129]
第三步骤：构建质子交换膜燃料电池中膜态水守恒方程。膜态水存在于电解质中，与质子的导电能力息息相关，分别在cl和pem中进行求解：
[0130]
[0131]
式中ρ
pem
(kg m-3
)和ew(kg mol-1
)分别表示膜的密度和当量重量；ω，λ和(m
2 s-1
)表示离聚物的体积分数，膜态水含量和有效扩散系数，与膜态水含量有关，也需要通过离聚物的体积分数进行修正；s
mw
(mol m-3 s-1
)为膜态水方程的源项，可以通过下式计算：
[0132][0133][0134]
式中，s
eod
(mol m-3 s-1
)和s
p
(mol m-3 s-1
)分别代表由于电渗拖曳和液压渗透产生的源项；
[0135]
作为优选，s
eod
和s
p
可由下式进行计算：
[0136][0137][0138]
式中，nd为电渗拖曳系数；(pa)和(pa)为阳极和阴极两侧的液压；k
mem
(m2)，δ
mem
(m)和δ
cl
(m)分别为质子交换膜的固有渗透率，厚度以及催化层的厚度；
[0139]
第四步骤：构建质子交换膜燃料电池中的能量守恒方程
[0140]
分别考虑气态，固态和液态物质的共同作用，由于多孔介质内的对流作用几乎可以不计，忽略该方程中的对流项，并在所有的计算区域内求解：
[0141][0142]
式中c
p，l
，c
p，g
和c
p，s
(kj kg-1 k-1
)分别表示液体，气体和固体的比热容；k
eff
(w m-1 k-1
)代表有效导热系数；ρg(kg m-3
)和ρs(kg m-3
)分别为气体和固体的密度；t(k)为温度；s
t
(w m-3
)为能量方程中的源项，分别由不同组件的反应热，欧姆热，相变热以及活化热组成：
[0143][0144]
式中j(a cm-3
)为垂直于膜方向上的体积电流密度；δsa和δsc(j mol-1 k-1
)分别为阴阳极两侧反应的熵差；η
ohm
和η
act
(v)分别为欧姆损失和活化损失；h(j mol-1
)为相变潜
热；bps表示阴阳极两侧的双极板。
[0145]
第五步骤：构建电化学模型，在恒电流计算模式下，电池的输出电压为可逆电压减去活化，欧姆和浓差过电势：
[0146]
v＝e
rev-η
ohm-η
act，con
ꢀꢀ
(5-1)
[0147][0148][0149]
其中e
rev
(v)为可逆电压；η
act，con
(v)为活化和浓差过电势；和(pa)为进口氢气和氧气的气压；和(ωm2)分别为电子通过极板和多孔电极骨架的面积比电阻；和分别代表质子通过催化层和膜中的面积比电阻；r(j mol-1 k-1
)为理想气体常数；i(a m-2
)为运行的电流密度。
[0150]
阳极侧的活化和浓差过电势由巴特勒-沃尔默方程(bulter-volmer)获得：
[0151][0152][0153]
其中j
0，a
和(a m-3
)分别为阳极体积交换电流密度和参考体积交换电流密度；(mol m-3
)为氢气浓度；(mol m-3
)为参考氢气浓度；αa和αc分别为阳极和阴极的转换系数，通常取0.5；ηa为阳极侧的活化和浓差过电势。
[0154]
第六步骤：构建阴极催化层结块模型。氧气从催化层孔隙中到pt催化剂表面间存在着巨大的局部传质阻力，这是高电流密度工况下引起较大浓差损失的主要原因。如图2所示，本发明中将这一过程中的浓差阻力分为以下四种：
[0155][0156]
其中r1(s m-1
)为氧气穿过液态水薄膜的扩散阻力；r
ion，int
(s m-1
)为氧气进入离聚物表面的界面阻力；(s m-1
)为氧气在离聚物中的扩散阻力；(s m-1
)为氧气从离聚物中到pt颗粒表面的界面阻力。
[0157]
进一步地，获得结块模型中的关键阻力参数。假设液态水均匀的覆盖于结块表面，液态水薄膜的扩散阻力r1计算如下：
[0158][0159]
其中(m
2 s-1
)为氧气在液态水中的扩散系数；为液态水的当量厚度：
[0160][0161]aion
＝4π(rc δ
ion
)2n
pt/c
ꢀꢀ
(6-4)
[0162]
其中a
ion
(m
2 m-3
)为离聚物的比表面积；rc(m)为碳载体的半径；δ
ion
(m)和n
pt/c
(m-3
)分别为离聚物膜的厚度以及单位体积pt/c的数目。可分别由以下模型获得：
[0163][0164][0165]
其中ε
ion
和ε
pt/c
分别为离聚物和pt/c的体积分数。本发明通过催化层中的铂载量和pt/c比来获得pt/c的体积分数。在计算离聚物的体积时考虑其吸/放水引发的体积膨胀与收缩。
[0166][0167][0168]
式中m
pt
(mg cm-2
)，ξ
pt/c
和ζ
i/c
代表铂载量，pt/c质量比和i/c质量比；ρ
pt
和ρc(kg m-3
)分别为铂和碳的密度；εc为碳载体的体积分数；催化层的孔隙率可有离聚物和pt/c的体积分数获得：
[0169]
ε
p
＝1-ε
ion-ε
pt/c
ꢀꢀ
(6-9)
[0170]
另外，离聚物界面的阻力可以通过以下关系获得：
[0171][0172]
式中k1为界面阻力系数；(m
2 s-1
)为氧气在离聚物中扩散的扩散系数，作为优选，可由下式获得：
[0173][0174]
假设pt催化剂均匀离散的分布在碳载体表面，则氧气在离聚物内部的扩散距离大于离聚物的平均厚度，需对其进行几何修正。相似地，氧气在pt表面的吸附阻力通过上述扩散阻力进行修正。
[0175]
[0176][0177]
式中θ
pt/o
为pt氧化的覆盖率；k2表示pt表面吸附阻力系数；r
pt
(m)表示pt颗粒的半径，可通过经验关联式计算获得：
[0178][0179]
式中a
ecsa
(m
2 kg-1
)为电化学反应面积。
[0180]
第七步骤，建立催化层结块模型与pemfc模型间的耦合关系。在稳态条件下，假设进入pt颗粒表面的氧气全部被电化学反应消耗，pt表面参加电化学反应的氧气浓度可通过通量守恒来确定：
[0181][0182]
其中(mol m-3
)和(mol m-3
)分别表示pt表面和离聚物界面上的氧气浓度；i
local
(a m-2
)为离聚物表面的局部电流密度，可通过体积电流密度jc和离聚物的比表面积进行折算，见公式(7-1)。
[0183][0184]
式中(pa m3mol-1
)为亨利系数；(mol m-3
)为催化层孔隙中的氧气浓度。
[0185]
作为优选，氧气在离聚物表面的亨利系数可由离聚物中的水含量确定：
[0186][0187]
由于阳极侧氢气的反应动力学远大于阴极侧的氧还原反应，本发明仅考虑阴极侧催化层内传质阻力的影响。基于此，仅需对阴极侧的bulter-volmer方程进行修正如下：
[0188][0189][0190]
其中i
0，c
和(a m-2
)为交换电流密度和参考交换电流密度；ηc(v)为阴极侧的活化和浓差过电势；w为能量参数；a
pt
(m2m-3
)为单位体积pt催化剂的表面积。
[0191]
作为优选，a
pt
可由电化学反应面积和pt载量计算获得：
[0192][0193]
第八步骤：构建上述控制方程的边界条件。本发明分别建立了组分扩散方程，液压方程，膜态水方程和能量方程的边界条件。其中组分扩散方程和液态方程的两侧边界分别
seidel方法对各个控制方程进行迭代求解。如图3所示，具体流程如下：(1)设置电流密度，输入运行和结构参数，对各物理量场进行初始化；(2)基于物理量更新模型中所有需要的物性和源项；(3)依次求解各组分扩散方程(h2，o2，n2，h2o)；(4)求解液压方程并内迭代获得多孔电极内的液态水饱和度；(5)求解膜态水守恒方程；(6)求解相互耦合的电化学模型和催化层结块模型；(7)判断是否到达极限电流密度，若达到则程序结束求解，若未达到则继续接下来的流程；(8)求解能量方程；(9)判断各控制方程是否达到残差标准，若达到则更新电流密度，以上一电流密度计算结果作为初始值开始下一电流密度的计算；若未满足收敛准则则对物理量进行超松弛或业松驰更新，返回流程(2)继续进行计算。
[0209]
为进一步理解本发明，图4为本发明一个实施例的模型验证极化曲线示意图。本实例中的基本运行工况为铂载量0.4mg cm-2
，阴阳极进口湿度分别为59％和84％。电解质质量分数为0.2，进口气体压力均为大气压。结果表明实验数据与模型计算数据吻合良好，证明了模型的可靠性。
[0210]
为进一步理解本发明的运行和结构参数，本发明的一个具体计算实例中，相关结构参数为：流道长宽高分别为0.1m，0.001m和0.001m；脊的宽度和高度为0.001m和0.0015m；gdl，mpl，cl以及pem的厚度分别为1.9
×
10-4
m，4
×
10-5
m，10-5
m，5.08
×
10-5
m；固有渗透率为10-11
，10-12
，10-13
和10-20
；agdl，ampl和acl的孔隙率分别是0.6，0.4和0.3；cgdl和cmpl的孔隙率为0.6和0.4，ccl的孔隙率根据催化层的具体构成来计算；gdl，mpl，cl的接触角为120
°
，110
°
和100
°
；铂载量为0.1mg cm-2
，pt/c比和i/c比分别为0.3和0.7。
[0211]
相关的运行参数为：运行温度353.15k；阴阳极进口压力均为2.0atm；进口相对湿度均为1.0；蒸发冷凝速率选择100s-1
，膜吸/放水速率选择为1.3s-1
；阴阳极的参考电流密度分别为1.5
×
10-5
a m-2
和20a m-2
；参考氢气和氧气浓度分别为56.4mol m-3
和3.39mol m-3
。
[0212]
基于上述具体实施例中的结构和运行参数，根据本发明所提出的模型和迭代求解方法进行逐步计算。本实施例下各控制方程在固定电流密度条件下对应的残差曲线均可降至1-10
以下，证明了模型和求解方法的合理性。图5为上述实施例的不同pt载量下的极化曲线，图6及图7为上述实施例的氧气传输阻力及pt催化剂表面氧气浓度随pt载量的变化曲线。该模型可以考虑详细的催化层构成对氧气传输和电池性能的影响，可为工程上质子交换膜燃料电池的优化设计提供指导。
[0213]
尽管以上结合附图对本发明的实施方案进行了描述，但本发明并不局限于上述的具体实施方案和应用领域，上述的具体实施方案仅仅是示意性的、指导性的，而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本发明权利要求所保护的范围的情况下，还可以做出很多种的形式，这些均属于本发明保护之列。