PEMFC阴极的LBM模拟第一部分：模型建立与计算区域.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093487 PEMFC 阴极的LBM 模拟 第一部分：模型建立和计算区域 陈 黎 栾辉宝 陶文铨 （西安交通大学能源与动力工程学院，动力工程多相流国家重点实验室，西安，710049） （联系电话：029，E-mail ：wqtao@mail.xjtu.edu.cn ） 摘要：在质子交换膜燃料电池（PEMFC ）阴极中，由于电化学反应速率较低，且反应气体传递缓慢， 阴极成为限制PEMFC 性能提高的主要因素之一。在阴极中发生着各种物理、化学、电学、电化学过 程。从孔隙尺度研究有助于深刻理解这些过程。本文利用格子－Boltzmann 方法（LBM ）研究了PEMFC 阴极中的流动、反应物和生成物的传质、质子传递和电化学反应过程。本文第一部分提出了综合模拟 阴极中发生的过程的VCCP －LBM 方案。第二部分讨论利用VCCP －LBM 模拟得到的结果以及阴极气 体通道中安装方块对流动传质和电池性能的影响。 关键词：PEMFC ；阴极；LBM ；电池性能；流场改善 0 引言 质子交换膜燃料电池（PEMFC ）的阴极氧气还原反应的速率远低于阳极氢气氧化反 应的速率，且阴极中反应气体传输过程缓慢， 因而阴极成为提高PEMFC 性能的主要限 [1] [2-8] 制因素之一 ，已成为众多学者的研究重点并提出了许多模型 。这些模型多为基于守 恒方程的宏观模型，其中 PEMFC 的多孔组件譬如GDL 常被处理为各向同性的多孔介 质，与实际情况相差很大[9] 。此外，在这些模型中应用了大量经验公式，这些公式本身 存在着很多不确定性。 因此，很有必要采用新的数值模拟方法基于孔隙尺度直接模拟阴极中详细的流动和 传质过程。20 世纪80 年代兴起的格子－Boltzmann 方法（LBM ）对研究多孔介质内的 流动传质问题非常有效[12] 。近来，许多学者采用LBM 方法研究了PEMFC 中的流动传 [10, 13-17] [9] 热传质问题 。Lang 等人 基于重建的三维GDL 结构模拟得到了GDL 厚度方向和 平面方向的渗透率。Wang 等人[13]采用LBM 研究了憎水GDL 中的两相流动问题。Park 等人[14]利用LBM 模拟了碳布GDL 中的单相流动问题并计算了GDL 的有效渗透率，他 们发现碳纤维的取向对GDL 的渗透率有很大影响。在其另一篇文献中[15] ，采用多相LBM 研究了液滴在碳纸GDL 中运动过程。Koido 等人[16]利用单相和多相LBM 模拟得到了不 同饱和度下GDL 的渗透率。Lang 等人[17]采用多相LBM 模拟了液态水在憎水气体通道 （GC ）中的传递过程。这些研究工作有助于深刻了解PEMFC 中发生的流动传质过程。 然而，上述文献紧紧局限于PEMFC 某一组件，并没有综合考虑PEMFC 中发生的 流动、传质和电化学反应过程。基于PEMFC 整体的模拟无疑会提供更完整和可靠的信 息。因此，本文采用LBM 综合模拟PEMFC 阴极中的流动，传质和电化学反应过程。 项目基金：国家自然科学重点项目（NO资助 1 LBM 模型介绍 1.1. 传质LBM 模型 对具有化学反应的传质来说，粒子分布函数演化方程为[18] ： 1 f x +e ∆t t −f x t =− f x t −f eq x t +S (1) i ( i , ) i ( , ) ( i ( , ) i ( , )) τ 其中S 和发生在CL 中的电化学反应有关，其取值在下文讨论。平衡态分布函数为： f k ,eq X k J =+K e ⋅u .