 Advances in Porous Flow 渗流力学进展, 2012, 2, 1-8 http://dx.doi.org/10.12677/apf.2012.21001 Published Online March 2012 (http://www.hanspub.org/journal/apf) Advance in Lattice Boltzmann Method Applied in Direct Methanol Fuel Cell* Weiwei Su, Yousheng Xu# Department of Physics, Zhejiang Normal University, Jinhua Email: #ysyou@zjnu.cn Received: Nov. 25th, 2011; revised: Dec. 20th, 2011; accepted: Dec. 21st, 2011 Abstract: The development of lattice Boltzmann method numerical simulation was reviewed in the article for the direct methanol fuel cell (DMFC). The up-to-the-minute progress of the lattice Boltzmann method applied in multiphase flow, mass and heat transfer and electrochemical reactions in the fuel cells was described in detail. As a computational fluid dynamic method, the lattice Boltzmann method is an effective and efficient numerical scheme to analyze the perform- ance of the fuel cell. The numerical simulations performed better understanding of the complex phenomena of fuel cells. These processes were needed for optimum call design. The latest research development was also reviewed in the article. At the same time, the deficiency of the present model was pointed out. Keywords: Lattice Boltzmann Method; Direct Methanol Fuel Cell; Polymer Electrolyte Membrane Fuel Cell 直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展* 苏卫卫，许友生# 浙江师范大学物理系，金华 Email: #ysyou@zjnu.cn 收稿日期：2011年11 月25 日；修回日期：2011 年12 月20 日；录用日期：2011 年12 月21 日 摘 要：综述了用格子 Boltzmann 方法模拟直接甲醇质子交换膜燃料电池的研究进展，具体介绍了利用格子 Boltzmann 方法模拟燃料电池中的多相流，传热传质和电化学反应过程等的最新进展，认识到格子Boltzmann 方 法作为一种计算流体力学的介观模拟方法，是一种有效的分析燃料电池性能的数值方法，通过数值模拟可以更 好地了解电池的复杂物理化学过程，为优化电池性能提供更充分的理论依据。文章还对研究最新进展予以简单 评述，同时指出了目前研究过程中所建模型存在的不足。 关键词：格子 Boltzmann 方法；直接甲醇燃料电池；质子交换膜燃料电池 1. 引言 随着社会的发展和人口的增长，能源消耗和环境 污染成为一个严峻的社会问题，如何实现可持续发 展，成为人类面临的重大问题。要从根本上解决能源 危机，这就需要实现从传统能源到新型能源的转型， 人类对未来能源的要求应该是清洁无污染，效率高， 经济安全和价格低廉。 燃料电池是一种可以直接将燃料中的化学能转 化为电能，能量转化效率高，环境友好型的电化学装 置。直接甲醇燃料电池(DMFC)属于燃料电池的一种， 使用甲醇，水和氧气作为燃料。DMFC 具有一般的燃 料电池特性外，还具有很多优势，如可在室温快速启 *基金项目：国家自然科学基金 10932010，11072220资助。 #通讯作者。 Copyright © 2012 Hanspub 1  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 动，能量密度高，系统简洁，无电解液流失等。 在燃料电池内发生很多物理化学现象，同时涉及 到传质、传热、多孔介质、多组分、多相、动量和电 化学反应。为了了解这些耦合的物理化学过程，分析 电池材料参数和特定电池设计对电池性能的影响，更 好地优化燃料电池系统，国内外许多研究工作组进行 了实验和数值模拟研究。但是由于燃料电池中的环境 是高效反应的，在电池运行时对其直接进行测量难度 非常大，所以数值模拟成为最常用的研究手段。Inoue 等对质子交换膜燃料电池气体扩散层中的两相流进 行数值研究[1]，林鸿等对质子交换膜燃料电池进行了 三维数值模拟[2]，此外像 Kulikov sky[3]，Ge 和Liu[4]， Xu[5]等对燃料电池中的各个结构进行了数值模拟。 格子 Boltzmann 方法作为计算流体的一种新工 具，可以应用到多孔介质流[6,7]、化学反应流[8]、磁流 体[9]、非牛顿流体[10-12]、多相多质流[13-19]等许多领域。 格子 Boltzmann 方法与其他传统数值方法相比，具有 并行性能好，物理图像清晰，程序易于实施和边界处 理简单等优点，尤其在解决复杂几何结构的问题时有 很大优势。因此，受到许多国内外学者的关注[20]。Joshi 等利用格子 Boltzmann方法模拟固体氧化物燃料电池 中的多组分气体输运现象[21] ；Tabe 等利用格子 Boltzmann 方法研究燃料电池中气体流管道中的两相 流[22]，还有许多数值研究者利用格子 Boltzmann 方法 模拟不同的运行参数对燃料电池的影响，例如 Fei 等 [23]，Hao 和Che ng[24]，Park 等[25]和Niu 等[26]。 2. 直接甲醇燃料电池的基本工作原理和结构 直接甲醇质子交换膜燃料电池主要包括极板、多 孔扩散层和膜电极(MEA)，其中质子交换膜是整个电 池系统的核心部分。直接甲醇燃料电池的基本工作原 理如图 1所示。 甲醇水溶液在阳极流道流动的过程中，穿过阳极 扩散层到达阳极催化层，同时，空气(氧气)在阴极流 道逐渐向阴极扩散层渗透到达阴极催化层。在阳极催 化层甲醇和水反应生成CO2，氢离子并释放出电子。 质子通过电解质膜到达阴极催化层，电子由外电路到 达阴极催化层。在阴极催化层质子和电子与氧气反应 生成水。反应式为： 32 CH OHHOCO 2 6H 6e (1) CH 3 OH + H 2 O CO 2 O 2 +H 2 O(蒸汽) 阳极扩散层 阳极催化层 质子交换膜 阴极催化层 阴极扩散层 H 2 O H + (xH 2 O) H + (yCH 3 OH) 电负 荷 e - e - Figure 1. Schematic view of the DMFC[27] 图1. DMFC工作原理示意图[27] 22 3O6H6e 3HO 2 (2) 3222 3 CH OHOCO2HO 2 0V (3) 由以上反应式可知，直接甲醇燃料电池中的阴极 反应跟质子交换膜燃料电池完全一样，在工作过程中 电池中的游动离子也仅为质子，因此，阴极催化剂均 为Pt/C，电解质膜也都为 Nafion 膜。 3. 直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法数值模拟 3.1. 格子 Boltzmann 方法方程体系 3.1.1. 宏观控制方程 在计算流体力学模拟中直接甲醇燃料电池阳极 和阴极的流场和组分浓度是一样的，下面会分别给出 模型的控制方程。 质量守恒方程： (4) 动量守恒方程： 2 i VVP rVS 2eff llll VCD CS (5) rε表示多孔介质的修正系数。 物料守恒方程： (6) l表示阳极的甲醇，水，CO2和阴极的 O2，水蒸气。 Copyright © 2012 Hanspub 2  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 Cl为组分 l的摩尔浓度，Sl为化学组分在电化学反应 过程中质量产生或消耗的速率。 l S , , 0, a c sj nF sj nF 阳极催化层 阴极催化层 其他区域 (7) 式中，n为组分l反应的得失电子数，s为反应的化学 计量系数，j为体积反应率。 电化学反应率： MeOH MeOH , a ref L aa Lref C jai C exp a F RT (8) 2 2 O O , c ref G cc Gref C jai C exp c F RT MeOH L 2 O G MeOH , (9) 式中，下标 a和c表示阳极和阴极，airef 为参考交换 电流密度， C为甲醇摩尔浓度， C为氧气摩尔 浓度， L ref 2 O ,Gref C为参考甲醇摩尔浓度，C为参考氧气 摩尔浓度，γ和α分别为浓度参数和传递系数，η为过 电位，T为温度，R为气体常数。 3.1.2. 格子 Boltzmann 方程 考虑了外力的流场的格子Boltzmann 方程为： ,i ,, 1 , ii i eq ii fxcett tfxt f xt f xt tF (10) , i fxt 是t时刻处于 x i t 位置，速度为 e的粒子的体平 均密度分布函数，τ是无量纲弛豫时间， 是时间的 增量。体平均平衡态分布函数 , eq i f xt 定义为 2 4 2 ii s ss e ecI cc 2 : 1 eq i ii uu eu f (11) i 为权重， 049 ，当 时，1, 2,3,4i19 i ，当 时， 5, 6, 7i,8 136 i ，cx t为格子速度， 3 s cc， 为多孔介质的孔隙度。 方程(10)中的外力项 i F 定义为[28,29]： 24 1 12 ii i ii eu eF Fw cc 2 ss eF uF c i (12) 采用 D2Q9 格子 Boltzmann 模型，离散的速度 e 为， i ec 0 0 1π1π cos,sin 1~ 4 22 29π29π 2cos,sin 5~8 44 i ii i ii i i i (13) 体平均密度和速度分别定义为， f , 2 001 uvcc cv , 2 i i i et vf G (14) 其中， 0 11 22 t c K ， F 12 t c K (15) 与流场分布函数不同，浓度场的格子 Boltzmann 方程 分布函数为 1 ,,, eq li ililili l il xctgxtgxtg xt tS g (16) 0,0 0 1,0,0,1 1 ~4 1,1 5~8 i i ci i (17) 相应的平衡态分布函数为[30-32]： 2 ,1 eq k lii l s cu gxt Cc (18) 方程(10) 和(16) 中的弛豫时间 和l 分别定义为 20.5 s ct 20.5 lls Dct l D 1, ,;1, , 和，为组分l 的扩散系数。 3.1.3. 边界处理 电池内流场和浓度场采用相同的边界处理格式， 我们以流场为例来说明边界处理方式。如图 2所示， 流体节点以实心圆圈表示，增加的虚拟流体节点以空 心圆圈表示，考虑均匀网格( x y iNjN 1, ,;0 x iNj ) 覆盖流场，水平方向采用非平衡外推格式，竖直方向 采用周期性边界处理格式。由此，需要在竖直方向上 增加两层网格(1, ,;)和( x iN 1 y jN )，代表虚拟流体节点。周期 性边界 处理格 式表达为 Copyright © 2012 Hanspub 3  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 阳 极 阴 极 计算区域 1 5 2 8 4 7 3 6 0 D2Q9 Figure 2. Schematic view of the periodic boundary 图2. 周期性边界处理示意图 ,7,8 , 1 y fiN 2,5,6 ,,0 4,7,8 4 ,1fi (19) 2,5,6 f iN fi ,7,8 , 1 y fiN 4 (20) 式中， 表示虚拟流体节点 上的分布函数 4,7,8 4 ,1fi ,1 y iN f 、7 f 和8 f 与流体节点 ,1i 上的分布函数 4 f 、7 f 和8 f 相等。 水平方向采用非平衡外推格式，考虑碰撞过程， 边界节点 O点的分布函数可以表示为 1 ,,1 eq ii i ,, eq i f OtfOtf Bt f Bt , xj 1, x BN j i (21) 式中，节点ON 为边界节点，节点 与 O点为相邻流体节点。 , f Ot 为碰撞后的分布函数。 3.1.4. 直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法数值模拟研究现状 燃料电池运行时涉及到多孔介质流、热质传递、 电化学反应和多组分多相流等复杂物理过程。采用格 子Boltzmann 方法模拟多相流，可以方便地描述出相 与相之间的相互作用；模拟多孔介质流，可以直接描 述孔隙内的流体流动。格子Boltzmann 方法作为一项 重要的数值计算方法，一直都在发展，在许多领域取 得了开拓性进展[33-38]。 在直接甲醇燃料电池中，甲醇水溶液在阳极催化 层反应生成CO2，如果 CO2气体大量聚集会堵塞流道， 降低燃料电池性能，所以如何去除 CO2气泡是人们较 为关注的问题。Fei 等对 DMFC 阳极流道中的甲醇溶 液和 CO2气体两相流动建立了物理模型如图3，并用 热格子 Boltzmann 模型(TLBM)对阳极微管道中的相 流进行了数值模拟，然后用商业软件对结果进行了验 证，主要目的是得到一个高效的计算方法来解决这一 CO 2 气泡 16.0 μm 1.5 μm x y 上边界（无滑移边界条件） 下边界（无滑移边界条件） 出 口 速 度 边 界 条 件 入 口 速 度 边 界 条 件 流动方向 （ ） （ ） Figure 3. Schematic diagram of the simulation domain of the bub- bly thermal flow in a microchannel in which x-axis is the flow direction and y-axis is opposite to the gravity. The microchannel is 16.0 um long and 1.5 um high. No-slip boundary conditions are applied at the upper and bottom walls, while velocity boundary conditions ar e imposed at the inlet and exit of the microchannel 图3. 微管道中气泡热流的计算区域示意图，x轴是流动方向，y 轴是重力的反方向。微管道长 16.0 um，高1.5 um。上表面和下表 面利用无滑移边界条件，管道的进口和出口施加速度边界条件[35] 问题[23]。研究表明：气泡的尺寸随着壁面温度的增加 而减小，在微管道中气泡尺寸的改变对气泡的流动性 有很明显的影响，小气泡相比较大气泡流动性比较强， 也就是说在壁面温度高的微管道中小气泡移动较快。 关于气泡输运速度和壁面温度之间的关系，TLBM 和 VOF(volume of fluid)两种数值方法得到了一致的结 论，就是即使气泡在入口速度为 0，但只要沿着管道 有温度梯度，气泡可以在水平微管道中自动地移动， 气泡的移动速度随着温度梯度的增加而增加。模拟结 果对 DMFC 中CO2去除技术的研究可以起到积极的作 用，同时证明了热格子 Boltzmann 方法模拟微管道中 气泡热流的可行性和高效性。但在这一工作中，他们 把DMFC 阳极流道看作是微管道进行建模，所建的物 理模型只是对 DMFC阳极流道的一个简化，而且所建 模型没有能够与扩散层及膜当中的传热传质过程较好 地结合起来，没有全面地描述流道中的流场形式。 Tabe 等用格子 Boltzmann方法模拟了质子交换膜 燃料电池阴极中反应物和产物的两相流，研究了水和 气体混合流的动力学特征[39]。在处理Poisson 方程的 离散，节点处的压强边界和入口处的粒子分布函数 时，改进了格子 Boltzmann 方法计算过程使之能够扩 展应用到气体扩散层这一多孔介质中。结果表明：1) 在质子交换膜燃料电池两相流中，考虑毛细管力，粘 力和惯性力的作用，LBM 可以模拟毛细管数 L Ca u 对液态水流动的影响。2) 通过对有无气 体流道两种气体扩散层中的两相流模拟，得到气体流 道和气体扩散层的润湿性对两相流的影响较大，进一 Copyright © 2012 Hanspub 4  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 步说明可以通过控制气体流道和多孔扩散层的润湿 性来得到液态水和空气的最佳分布，从而使催化层的 氧气供应充分。可以看出，文章在研究液态水和空气 之间的相互作用对燃料电池性能影响的现象中，体现 了考虑整个系统物理性质的格子Boltzmann 方法和适 当地选择边界条件的优越性。但是他们的建模没有包 含阳极，催化层和膜，而且模拟中设定温度为恒温， 是一个严重的模拟局限。温度分布的不同会对水的相 变数量有影响，因此不能考虑带有相变的物理现象。 Hao 和Cheng 利用多相自由能格子 Boltzmann 方 法模拟了质子交换膜燃料电池中穿过疏水性气体扩散 层微孔的水滴的形成，并且研究了水滴在气体剪切力 作用下在气体扩散层表面的运动情况[24]。呈现了水滴 的出现，生长，脱落和运动过程，研究了在气体流动 速度和气体扩散层表面润湿性影响下的脱落水滴的尺 寸和去除管道中水滴的时间。结果表明在疏水性强的 气体扩散层表面上，气体流速越大，水滴越容易去除。 此外，利用基于力平衡的解析模型求解了水滴脱落尺 寸，结果跟 LBM 模拟的结果相一致。说明 LBM 是一 种有效的工具来研究 PEMFC 气体流道的水输运现象。 质子交换膜燃料电池的电极通常是由碳布或碳 纸构成的多孔介质。Park 等用格子 Boltzmann 方法模 拟了流体在多孔电极中的流动[40]，计算了多孔介质的 有效渗透率，并与文献[41]中的解析解作对比，结果 一致如图 4。他们还建立了三维纤维多孔介质模型来 研究纤维屑的排列方向对有效渗透率的影响，结果表 明平行与平均流动方向的有效渗透率比垂直于平 1E-4 1E-30.01 1E-3 0.01 0.1 1 0.1 有效渗透率 纤维屑渗透率 Lubrication Model LBM 1 0.1 0.01 1E-3 1E-3 1E-4 0.01 0.1 Figure 4. Comparison of the effective permeability between the LB simulation (circles) and the lubrication model(square), as a func- tion of the dimensionless tow permeability; the nominal poros- ity(εnom) of porous me dium is 0.4 图4. LB模拟和 lubrication 模型的有效渗透率比较[40] 均流动方向的有效渗透率大很多，接近于两倍。 Niu 等用多相，多弛豫时间格子 Boltzmann 模型 模拟研究了质子交换膜燃料电池系统气体扩散层中 的水气输运过程[26]。这一模型基于平均场扩散界面理 论，能够处理大密度比率和不同粘性下的多相流。无 滑移边界和润湿边界分别采用标准的反弹边界条件 和近似平均方法进行处理，能够详细地描述气体扩散 层中的液气输运。不像大多数的格子 Boltzmann 方法 是基于 Bhatnagar-Gross-Krook 碰撞项，文中模型采用 与粘性无关的速度场，这对模拟不同粘性共存的多相 流是非常重要的。通过模拟在润湿壁面上的静电液 滴，并且与理论计算相比较，验证了这一模型的正确 性。然后，文章模拟了质子交换膜燃料电池系统气体 扩散层中的水气流动，研究了在不同条件下与饱和度 有关的输运性质，模拟结果跟文献中的数值模拟和实 验得到的结果相吻合[42,43]。 质子交换膜燃料电池的气体扩散层是由碳布或碳 纸构成的多孔介质。在有关质子交换膜燃料电池的文 献中一般把气体扩散层作为渗透率为常量的各向同性 多孔介质来建模。实际上，纤维的排列方向造成的碳 纸的各向异性影响多孔介质的渗透率。 Park 和Li 通过 格子 Boltzmann 方法解决了质子交换膜燃料电池各向 异性气体扩散层中的多相流现象[25]。把碳布和碳纸的 各向异性孔隙结构分为空隙空间和孔隙空间进行建 模，利用 Stokes/Brinkman 公式求解，从各种显微图像 得到空隙空间和不可渗透的纤维分布。计算得到的多 孔介质的渗透性和文献[44]中的实验数据相一致。各 种纤维分布的模拟结果说明纤维的排列方向影响介质 的渗透率。模拟多相流的格子 Boltzmann 模型是基于 多组分格子 Boltzmann 模型和粒子间相互作用力的结 合，该模型成功地模拟了在多孔介质中液滴运动的复 杂不稳定性，为质子交换膜燃料电池多孔电极中液态 水聚集和去除机理的研究提供了有力的工具。 在质子交换膜燃料电池阴极中，发生着各种化 学、物理、电学、电化学等过程，反应物传递速度慢， 电化学反应速率低是影响PEMFC 性能的主要因素， 需要从孔隙尺度来深入地研究这些过程。陈黎等采用 格子 Boltzmann方法直接从介观层次模拟了质子交换 膜燃料电池阴极催化层的传质和电化学反应过程[45]。 建立了理想二维规则模型，为了验证结果的正确性， 模拟所得的极化曲线和文献[46]的极化曲线吻合较 Copyright © 2012 Hanspub 5  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 好。用电镜扫描 PEMFC 阴极催化层，处理图象得到 了一个实际二维计算模型。利用格子 Boltzmann 方法 进行模拟，得到了二维的氧气浓度和电势分布。表明 利用格子 Boltzmann 方法计算分析PEMFC 阴极催化 层中传递和电化学反应是行之有效的。传统的宏观模 型考虑象曲折度、孔隙率、催化层组分和催化层厚度 等宏观参数对燃料电池催化层性能的影响，这样建立 的模型跟实际结构相比较有一定的差别，因此不能充 分说明催化层结构对燃料电池性能的影响。而 LBM 直接基于实际结构对传质和电化学反应过程进行研 究，不依赖于曲折度、孔隙率这些宏观参数。 0 00.1 0.15 0.2 0.05 –0.1 –0.2 0 –0.1 –0.2 y/H x/H Anode Diffusion Layer Anode Catalyst Layer Membrane 上面介绍了用格子 Boltzmann 方法模拟微管道， 气体扩散层，阴极催化层，但仅仅局限于燃料电池的 某一组件，没有综合考虑电池中发生的传质和流动， 基于燃料电池整体的模拟无疑会得到更完整、可靠的 信息。Delavar 等用格子 Boltzmann 方法，采用 D2Q9 模型对 DMFC全电池进行数值模拟，计算区域包括阳 极流道，阳极催化层，阳极扩散层，质子交换膜，阴 极流道，阴极催化层和阴极扩散层[47]。这一模型通过 模拟反应物和生成物分别在空流道和多孔流道中的流 动模式和浓度场分布来研究电池的性能。从图 5可以 看出，在多孔介质区域流场改变并且速度减小，流体 在阴极流道的速度比在阳极流道快。在模拟流场时， 由于这种速度的不同，存在两个流道交界面的收敛问 题，所以阳极和阴极采用不同的分布函数。从图 6可 以看出，在反应催化层周围，反应物的浓度减小，生 成物的浓度增大。计算结果得到了合理的流场分布和 CO2及甲醇浓度分布，文章还得到了在一定操作条件 1 0 0105 –1 1 0 –1 y/H x/H C a 15 Anode Channel Anode Diffusion and atalyst Layers Cathode Diffusion nd Catalyst Layers Cathode Channel Membrane Figure 5. Velocity contours (U/Umax-anode) for all seven parts of DMFC[47] 图5. DMFC的速度分布(U/U max-a node)[47] (a) 0 00.1 0.15 0.2 0.05 –0.1 –0.2 0 –0.1 –0.2 y/H x/H Anode Diffusion Layer Anode Catalyst Layer Membrane (b) Figure 6. Ins tantaneous carbon dioxide and methanol concentra- tion contours (kgm–3) for anode channel, diffusion and catalyst layers, membrane 图6. 阳极流道，扩散层，催化层和质子交换膜上的瞬时 CO2(a) 和甲醇(b)浓度分布[47] 下的极化曲线，我们可以在此基础上进一步讨论燃料 电池的结构参数和操作条件对电池性能的影响。结果 表明，这一模型能够比较精确地研究 DMFC 的流动模 式和性能，为进一步的研究打下良好基础。 4. 用格子 Boltzmann 方法模拟直接甲醇 质子交换膜燃料电池存在的问题 1) 研究人员用格子 Boltzmann 方法研究了两相 流现象，但是在阴极伴随有水的蒸发和冷凝过程，人 们均忽略了水的相变过程，没有很好地考虑相变与传 热之间的相互影响，没有综合考虑到催化层反应动力 学对传质的影响以及流道中水蒸发和凝结过程中吸 收和释放的潜热的影响。阴极淹没和膜的干涸等现象 都是伴随水的相变过程而产生的一些非正常现象，如 果考虑了相变过程进行深入研究将会是十分有意义 的。 2) 由于进行三维模拟时计算复杂性等因素，进行 Copyright © 2012 Hanspub 6  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 三维数值模拟的工作相对较少，建立的模型大多数都 是二维的，而且进行了简化，忽略了电池内的一些特 殊情况，不能详细地描述电池内的传质传热和电化学 反应过程。 3) 人们把格子 Boltzmann 方法数值模拟的研究 重点放在了两相流上面，忽略了工作温度，空气压力， 甲醇水溶液浓度，流量，电池结构，运行工况等参数 对电池性能的影响。例如，当电池作为交通工具等电 源时，经常处于变工况的工作状态下，了解电池的动 态特性对优化电池的性能也是十分重要的。 5. 结束语 格子 Boltzmann方法作为一种计算流体力学的方 法，应用在许多传统模拟方法不能胜任的领域，比如 多相流，复杂几何结构和伴随有化学反应的流动过 程。人们对直接甲醇质子交换膜燃料电池建立不同的 模型，利用格子 Boltzmann 方法对其进行模拟计算， 涉及多孔介质，微管道，热流，多相流，还有直接对 DMFC 整体进行建模。通过数值模拟，详细揭示了电 池内复杂流动结构、传热传质过程、电化学反应动力 学及其对电池性能的影响，可以为优化燃料电池的性 能提供有力的理论依据，也更具有实际应用价值。但 同时为了更进一步深层次的了解燃料电池，提高其性 能，要综合考虑相变和传热、传质间的相互影响，催 化层内的电化学反应动力学和传质、电荷传递间的耦 合作用。因此，我们还需要结合相变，两相流和电化 学反应等建立先进的燃料电池综合模型，利用格子 Boltzmann 方法的优势对其进行数值模拟。 6. 致谢 感谢给予转载和引用权的资料、图片、文献、研 究思想和设想的所有者。 参考文献 (References) [1] G. Inoue, T. Yoshimoto, Y. Matsukuma, et al. Development of simulated gas diffusion layer of polymer electrolyte fuel cells and evaluation of its structure. Journal of Power Sources, 2008, 175(1): 145-158. [2] 林鸿, 陶文铨. 质子交换膜燃 料电池的三维数值 模拟[J]. 西 安交通大学学报. 2008, 42(1): 41-45. [3] A. A. Kulikovsky, J. Divisek and A. A. Kornyshev. Two-dimen- sional simulation of direct methanol fuel cell—A new (embed- ded) type of current collector. Journal of the Electrochemical Society, 2000, 147(3): 953-959. [4] J. Ge, H. Liu. A three-dimensional mathematical model for liq- uid-fed direct methanol fuel cells. Journal of Power Sources, 2006, 160(1): 413-421. [5] C. Xu, A. Faghri, X. Li, et al. Methanol and water crossover in a passive liquid-feed direct methanol fuel cell. International Jour- nal of Hydrogen Energy, 2010, 35(4): 1469-1777. [6] G. Tang, W. Tao and Y. He. Gas slippage effect on microscale porous flow using the lattice Boltzmann method. Physical Re- view E, 2005, 72(5): Article ID 056301. [7] Z. L. Guo, T. S. Zhao. Lattice Boltzmann model for incom- pressible flows through porous media. Physical Review E, 2002, 66(32B): 036301-036304. [8] S. Chen, S. P. Dawson and G. D. Doolen, et al. Lattice methods and their applications to reacting systems. Computers and Chemi- cal Engineering, 1995, 19(6): 617-646. [9] S. Chen, H. D. Chen, D. Martinez, et al. Lattice Boltzmann model for simulation of magnetohydrodynamics. Physical Re- view Letters, 1991, 67(27): 3776-3779. [10] N. Rakotomalala, D. Salin and P. Watzky. Simulations of viscous flows of complex fluids with a Bhatnagar, Gross, and Krook lat- tice gas. Physics of Fluids, 1996, 8(11): 3200-3202. [11] E. S. Boek, J. Chin and P. V. Coveney. Lattice Boltzmann simu- lation of the flow of non-Newtonian fluids in porous media. In- ternational Journal of Modern Physics B, 2003, 17(1-2): 99-102. [12] S. Gabbanelli, G. Drazer and J. Koplik. Lattice Boltzmann mot- hod for non-Newtonian (powerlaw) fluids. Physical Review E, 2005, 72: Article ID 046312. [13] D. Grunau, S. Chen and K. Eggert. A Lattice Boltzmann model for multiphase fluid-flows. Physics of Fluids, 1993, 5: 2557- 2562. [14] A. Gunstensen, D. Rothman, S. Zaleski, et al. Lattice Boltzmann model of immiscible fluids. Physical Review E, 1991, 43(8): 4320-4327. [15] 李维仲, 李爽. 用格子 Boltzmann 方法模拟液滴合并过程[J]. 热科学与技术, 2007, 6(3): 379-393. [16] M. Swift, E. Orlandini, W. Osborn, et al. Lattice Boltzmann simulations of liquid-gas and binary fluid systems. Physical Re- view E, 1996, 54: 5041-5052. [17] T. Inamuro, T. Ogata, S. Tajima, et al. A lattice Boltzmann mot- hod for incompressible two-phase flows with large density dif- ferences. Journal of Computational Physics, 2004, 198(2): 628- 644. [18] L. Luo. Theory of the lattice Boltzmann mothod: Lattice Boltz- mann models for nonideal gases. Physical Review E, 2000, 62 (4): 4982. [19] X. He, S. Chen and R. Zhang. A lattice Boltzmann scheme for incompressible multiphase flow and its application in simulation of Rayleigh-Taylor instability. Journal of Computational Physics, 1999, 152(2): 642-663. [20] C. L. Tien, A. Majumdar, V. P. Carey, et al. Molecular and mi- croscal transport phnomena: A report on the 2nd US Japan joint seminar, Santa Barbara, California, 1996. Microscale Thermo- physical Engineering, 1997, 1(1): 71-84. [21] A. S. Joshi, K. N. Grew, A. A. Peracchio, et al. Lattice Boltz- mann modeling of 2D gas transport in a solid oxide fuel cell an- ode. Journal of Power Sources, 2007, 164(2): 631-638. [22] Y. Tabe, T. Ochi, K. Kikuta, et al. Numerical simulation of liquid water behavior in separator-channels in PEMFC using LBM. Proceedings of 3rd International Conference on Fuel Cell Sci- ence Engineering and Technology, Ypsilanti, 2005: 1-6. [23] K. Fei, T. S. Chen and C. W. Hong. Direct methanol fuel cell bubble transport simulations via thermal lattice Boltzmann and volume of fluid methods. Journal of Power Sources, 2010, 195(7): 1940-1945. [24] L. Hao, P. Cheng. Lattice Boltzmann simulations of liquid drop- let dynamic behavior on a hydrophobic surface of a gas flow channel. Journal of Power Sources, 2009, 190(2): 435-446. [25] J. Park, X. Li. Multi-phase micro-scale flow simulation in the electrodes of a PEM fuel cell by lattice Boltzmann method. Copyright © 2012 Hanspub 7  直接甲醇质子交换膜燃料电池的格子 Boltzmann 方法研究进展 Copyright © 2012 Hanspub 8 Journal of Power Sources, 2008, 178: 248-257. [26] X. D. Niu, T. Munekata, S. A. Hyodoa, et al. An investigation of water-gas transport processes in the gas-diffusion-layer of a PEM fuel cell by a multiphase multiple-relaxation-time lattice Boltzmann model. Journal of Power Sources, 2007, 172(2): 542- 552. [27] 陈志. 直接甲醇质子交换膜燃料电池中阳极 CO2去除及水、 热的综合管理[D]. 东南大学, 2006. [28] A. Cancelliere, C. Chang, E. Foti, et al. The permeability of a random medium: Comparison of simulation with theory. Physics of Fluids A, 1990, 2(12): 2085-2088. [29] K. J. Malory, L. Furuberg, J. Feder, et al. Dynamics of slow drainage in porous media. Physical Review Letters, 68(14): 2161-2164. [30] M. A. Delavar, M. Farhadi and K. Sedighi. Effect of the heater location on heat transfer and entropy generation in the cavity using the lattice Boltzmann Methode. Heat Transfer Research, 2009, 40(6): 521-536. [31] A. Mezrhab, M. Jami, C. Abid, et al. Lattice-Boltzmann model- ling of natural convection in an inclined square enclosure with partitions attached to its cold wall. International Journal of Heat Fluid Flow, 2006, 27(3): 456-465. [32] A. A. Mohammad. Applied lattice Boltzmann method for trans- port phenomena momentum heat and mass transfer. Calgary: The University of Calgary Press, 2007. [33] S. Chen, G. D. Doolen. Lattice Boltzmann method for fluid flows. Annual Review of Fluid Mechanics, 1998, 30(1): 329- 364. [34] B. Chopard, P. O. Luthi. Lattice Boltzmann computations and applications to physics. Theory Computational Physics, 1999, 217: 115-130. [35] J. R. Ho, C. P. Kuo and W. S. Jiaung. Study of heat transfer in multilayered structure within the framework of dualphaselag heat conduction model using lattice Boltzmann method. International Journal of Heat Mass Transfer, 2003, 46(1): 55-69. [36] J. R. Ho, C. P. Kuo, W. S. Jiaung, et al. Lattice Boltzmann scheme for hyperbolic heat conduction equation. Numerical Heat Trans- fer B, 2002, 41(6): 591-607. [37] R. R. Nourgaliev, T. N. Dinh, T. G. Theofanous, et al. The Lattice Boltzmann equation method: Theoretical interpretation, numerics and implications. International Journal of Multiphase Flow, 2003, 29(1): 117-169. [38] R. G. M. VonderSman, M. H. Ernst and A. C. Berkerbosch. Lattice Boltzmann scheme for cooling of packed cut flowers. International Journal of Heat Mass Transfer, 2000, 43(4): 577- 587. [39] Y. Tabe, Y. Lee, T. Chikahisa, et al. Numerical simulation of liquid water and gas flow in a channel and a simplifed gas diffu- sion layer model of polymer electrolyte membrane fuel cells us- ing the lattice Boltzmann method. Journal of Power Sources, 2009, 193(1): 24-31. [40] J. Park, M. Matsubara and X. Li. Application of lattice Boltz- mann method to a micro-scale flow simulation in the porous electrode of a PEM fuel cell. Journal of Power Sources, 2007, 173: 404-414. [41] F. R. Phelan, Jr., G. Wise. Analysis of transverse flow in aligned fibrous porous media. Composites, 1995, 27A: 25. [42] H. Li, C. Pan and C. T. Miller. Pore-scale investigation of vis- cous coupling effects for two-phase flow in porous media. Physi- cal Review E, 2005, 72(4): Article ID 026705. [43] H. Chen, S. Chen and W. H. Matthaeus. Recovery of the Navier- Stokes equations using a lattice-gas Boltzmann method. Physical Review A, 1992, A45: 5539. [44] J. P. Feser, A. K. Prasad and S. G. Advani. Experimental charac- terization of in-plane permeability of gas diffusion layers. Jour- nal of Power Sources, 2006, 162(2): 1226-1231. [45] 陈黎, 栾辉宝, 陶文铨. PEMFC阴极催化层的 LBM模拟[J]. 电源技术, 2009, 33(9): 794-797. [46] G. Wang, P. P. Mukherjee and C. Y. Wang. Direct numerical simulation (DNS) modeling of PEFC electrodes Part I. Random micro-structure. Electrochimica Acta, 2005, 51(15): 3151-3160. [47] M. A. Delavar, M. Farhadi and K. Sedighi. Numerical simula- tion of direct methanol fuel cells using lattice Boltzmann method. International Journal of Hydrogen Energy, 2010, 35(17): 9306- 9317. |