质子交换膜燃料电池流场设计（节选自课程设计作业）

3.2 仿真建模

3.2.1 几何

PEMFC电池组由多个单体电池堆叠而成，单体电池结构如图3-1右图所示，燃料气和氧气分别从阳极、阴极流道进入，穿过多孔扩散层和多孔电极到达质子交换膜附近，在催化剂的作业下发生电极反应，电子则通过流板导出，经过负载到达正极。在此，PEMFC电池组首先被抽象为如图3-1右图的单体电池，单体电池平行流道则可抽象为平行的单通道，如图3-2所示，

图3-1 燃料电池组示意图

图3-2流道模型示意图
（1）流道。最下方的长方体表示阳极流道，通燃料气，流道两侧的空间表示导电流场板，因其电阻相比与多孔电极可忽略，所以流场板与多孔电极接触界面的电位可视为电池负极电位。
（2）气体扩散层。阳极流道上方紧邻的长方体表示气体扩散层（GDL），气体扩散层内不发生电化学反应，仅仅存在传质与动量传递，气体扩散层一部分面积与流道接触，作为质量交换的界面，另一部分面积与流场板接触，起到传导电子的作用。
（3）多孔电极。与气体扩散层相邻的长方体表示多孔电极，在多孔电极中既发生电化学反应也存在质量传递与动量传递。
（4）质子交换膜。在两个多孔电极之间的部分表示质子交换膜，膜内存在质子传质，而气体无法通过。

3.2.2 物理场设置

（1）二次电流分布
二次电流分布用于求解电流分布，因此所有能够导电的部件都应考虑。二次电流分布由气体扩散层、多孔电极、质子交换膜三部分构成。气体扩散层设置为具有一定电导率为的电极，多孔电极设置为包含多孔电极反应的电极，其中电解质电导率为，电极电导率为，电解质体积分数为。质子交换膜设置为电导率为的电解质。

（2）气体动量传递
动量传递作用于流道、气体扩散层和多孔电极，对于流道，通过Navier-Stokes方程求解流速u和压力p，对于气体扩散层和多孔电极，必须要考虑渗流作用，因此使用Brinkman方程进行求解。流道的一端设置为“入口”，另一端设置为“出口”，多孔电极侧面边界设置为“对称”，其余边界设置为“壁