基于格子玻尔兹曼方法的梯度化扩散层降解研究.pdf

第 12 卷 第 2 期 新 能 源 进 展 Vol.12 No.2 2024 年 4 月 ADVANCES IN NEW AND RENEWABLE ENERGY Apr.2024 *收稿日期：2023-07-24 修订日期：2023-11-07 基金项目：国家自然科学基金项目（52176084）；国家重点研发计划项目（2022YFE0207600）；天津市杰出青年基金项目（23JCJQJC00290）；清华大学智能绿色车辆与交通全国重点实验室开放课题项目（KFY2403）通信作者：王誉霖，E-mail： 引用本文：齐丽霞,许浩锴,董晓燕,等.基于格子玻尔兹曼方法的梯度化扩散层降解研究J.新能源进展,2024,12(2):175-181.Citation：QI Lixia,XU Haokai,DONG Xiaoyan,et al.Gradient diffusion layer degradation based on lattice Boltzmann methodJ.Advances in new and renewable energy,2024,12(2):175-181.文章编号：2095-560X（2024）02-0175-07 基于格子玻尔兹曼方法的梯度化扩散层降解研究*齐丽霞，许浩锴，董晓燕，秦世位，王誉霖（天津商业大学，天津市制冷技术重点实验室，天津 300134）摘 要：采用随机算法重建了不同老化程度的梯度化气体扩散层二维微观结构。随后，采用多相伪势格子玻尔兹曼方法，研究不同降解程度的梯度设计扩散层内液态水的传输行为。结果表明，随着扩散层降解的加剧，扩散层排水性能显著降低，大量液态水残留在扩散层内。此外，梯度化设计的扩散层在进口区域初始分布较多聚四氟乙烯有利于扩散层的水管理，同时在老化后还具有较高的毛细压力进而控制水驱面积，展现出较高的排水性能及抗老化性能。研究结果可为制作高抗老化性能和高排水性能的扩散层提供研究基础和设计思路。关键词：气体扩散层；格子玻尔兹曼方法；老化；排水性能 中图分类号：TK91 文献标志码：A DOI：10.3969/j.issn.2095-560X.2024.02.008 Gradient Diffusion Layer Degradation Based on Lattice Boltzmann Method QI Lixia,XU Haokai,DONG Xiaoyan,QIN Shiwei,WANG Yulin(Tianjin Key Lab of Refrigeration Technology,Tianjin University of Commerce,Tianjin 300134,China)Abstract:This study employed a random algorithm to reconstruct the two-dimensional microstructure of gradient gas diffusion layers with varying degrees of aging.Subsequently,the multiphase pseudo-potential lattice Boltzmann method was utilized to investigate the transport behavior of liquid water in gradient-designed diffusion layers with varying degrees of degradation.The results indicate that as the degradation of the diffusion layer intensifies,the drainage performance significantly decreases,and a considerable amount of liquid water remains in the diffusion layer.In addition,the gradient-designed diffusion layer initially distributes more polytetrafluoroethylene in the inlet area,enhancing water management within the diffusion layer.At the same time,after aging,the diffusion layer exhibits higher capillary pressure to control the water-driving area and shows higher drainage and anti-aging performance.This study provides a research basis and design for making diffusion layers with high aging resistance and drainage performance.Keywords:gas diffusion layer;lattice Boltzmann method;aging;drainage performance 0 引 言 质子交换膜燃料电池是一种通过电化学反应将反应物中的化学能直接转化为电能的装置，其具有环保、高效、响应快等优点，因此被认为是最有前景的能源技术之一1。然而，质子交换膜燃料电池大规模商业化仍面临巨大挑战2-3，其耐久性是主要障碍之一4。近年来，人们对燃料电池催化层和双极板的老化进行了深入的研究。然而，气体扩散层（gas diffusion layer,GDL）的老化问题却很少受到重视。最近的研究表明，扩散层的降解不容忽视。扩散层主要有以下两个功能：为反应物输送至催化层提供气体通道，并将反应产生的液态水排出燃料电池5。扩散层是一种复杂的微孔结构，通常由碳纸或碳布制成，其内部是由多个碳纤维交错叠加而成的多孔微观结构。为了获得疏水特性和改善排176 新 能 源 进 展 第 12 卷 水性能，通常在制备扩散层时加入一定质量分数的聚四氟乙烯（polytetrafluoroethylene,PTFE）。聚四氟乙烯侵入扩散层并附着在碳纤维表面，其中交错排列的碳纤维形成的小孔更容易被聚四氟乙烯侵占。同时，真实的聚四氟乙烯结构会占据扩散层的孔隙，降低孔隙率。扩散层的排水性能主要取决于其多孔结构和疏水特性6。MOLAEIMANESH 等7通过格子玻尔兹曼方法（lattice Boltzmann method,LBM）研究了有无聚四氟乙烯的扩散层。结果表明，具有疏水特性的扩散层能极大地促进液体的排出。然而，在燃料电池运行过程中，扩散层的结构会因老化而损坏。老化的两个原因为碳纤维的化学降解和聚四氟乙烯的机械降解。在燃料电池的运行过程中，碳纤维与液体发生反应并被氧化导致降解，也称为碳腐蚀。降解导致扩散层中的孔隙形态和表面特征发生变化，出现裂缝和大孔隙，进而改变液态水的流动行为。聚四氟乙烯降解表现为从碳纤维表面脱落，导致碳纤维疏水性能丧失，液体在亲水性碳纤维上的流动行为改变。YU 等8通过扫描电子显微镜、X 射线光电子能谱和热重分析对扩散层进行了分析。结果表明，扩散层老化会破坏碳纤维，导致扩散层的疏水性明显降低，但没有深入研究碳纤维结构的变化对排水性能的影响机理。此外，通过数值模拟研究扩散层在不同老化程度下排水性能的相关文献较少。基于以上问题，本研究采用随机算法建立聚四氟乙烯浓度梯度化设计的扩散层二维结构。同时，对扩散层进行一定程度的降解，得到扩散层降解前后的二维结构。随后，利用 SHAN 等9提出的多相伪势格子玻尔兹曼方法研究梯度化扩散层多孔结构中的气液两相流动行为。通过对扩散层进行聚四氟乙烯双梯度化设计和三梯度化设计，对比常规扩散层获得其排水性能和抗老化性能的提升，进而对不同老化程度的梯度化扩散层的排水性能进行较为准确的评价，本研究旨在通过梯度化结构提升扩散层的抗老化性能以及气体扩散能力。1 扩散层重建方法 为了获得与激光扫描共聚焦显微镜（confocal laser scanning microscope,CLSM）和 3D-X 射线断层 扫 描 观 察 到 的 相 似 的 GDL 以 及 微 孔 层（microporous layer,MPL）微观结构，依据 DAINO等10采用的随机算法构建扩散层的二维微观结构，如图 1（a、b）所示。扩散层中的碳纤维被视为直径为 7 m、长度假定为无限长的圆柱体，是随机生成和交错的，相对于扩散层的厚度方向可被忽略。扩散层的初始孔隙率设定为 65%，聚四氟乙烯的质量分数设定为 10%。扩散层计算域的尺寸大小为 200 m 100 m，每个体素为 1 m 1 m，晶格域为 200 100 体素。图 1 CLSM（a）和 3D-X 射线断层扫描（b）下的 GDL 微观结构；梯度化扩散层及老化最终扩散层二维微观结构（c f）Fig.1 GDL microstructure under CLSM(a)and 3D-X-ray tomography(b);two-dimensional microstructure of gradient diffusion layer and aging final diffusion layer(c-f)第 2 期 齐丽霞等：基于格子玻尔兹曼方法的梯度化扩散层降解研究 177 采用随机算法构建了沿厚度方向呈梯度分布的聚四氟乙烯扩散层，包括双梯度化分布和三梯度化分布。双梯度化扩散层的进口（靠近催化层侧）和出口区域（靠近流场侧），三梯度化扩散层的入口、中间和出口区域分别具有不同的聚四氟乙烯含量。随后，将常规扩散层和梯度化扩散层进行 20 次老化处理，为方便观察，研究了具有代表性的初始状态、老化前期状态、老化中期状态、老化后期状态、最终老化状态等 5 个阶段下的液态水传输行为，相应编号如表 1 所示。同时，假设气体扩散层中碳纤维的最终降解量为 30%，聚四氟乙烯的降解量为 15%10-11。图 1（c f）显示了聚四氟乙烯总含量为 10%的梯度化气体扩散层以及老化后扩散层的代表性重建微观结构。表 1 聚四氟乙烯含量为 10%的常规扩散层及梯度化扩散层设计 Table 1 Design of conventional diffusion layer and gradient diffusion layer with PTFE content of 10%设计 方案编号 聚四氟乙烯分布含量/%常规扩散层 1 5 10（均匀分布）双梯度扩散层 6 10 12（进口）；8（出口）11 15 14（进口）；6（出口）16 20 16（进口）；4（出口）三梯度扩散层 21 25 12（进口）；10（中间）；8（出口）26 30 14（进口）；10（中间）；6（出口）31 35 16（进口）；10（中间）；4（出口）2 数值模型 2.1 控制方程 与计算流体动力学的传统数值方法相比，格子玻尔兹曼方法具有更好的数值稳定性以及结构多样性12，在研究具有复杂多孔结构的气体扩散层内的两相流方面更具优势。因此采用多相伪势格子玻尔兹曼方法9研究含聚四氟乙烯的扩散层多孔微观结构中的两相气液流动问题。使用单弛豫时间碰撞算子和 D2Q9 晶格模型。分布函数的演化方程为：()()()()eq,1,iiikkiifxet ttfx tfx tfx t+=(1)式中：t 为时刻；t为时间步长；ie为晶格速度；x为粒子位置；为松弛参数。平衡态分布函数(),kifx teq,为：()()()2eqeq24ss2eq2s111212iikkiie ue uccfuc+=eq,(2)式中：cs为音速；i为权重因子；为流体的宏观密度；u 为宏观速度，对于 D2Q9 模型，s3cc=，其中 c=x/t，参考文献13。为了显示流体对固体的不同润湿性，引入外力项kF。kF为组分 k 上的外力之和，包括表面张力F1k，流固作用力 F2k和外力 F3k。kkkkFFFF=+123 (3)流体表面张力 F1k表示为：()()()()(),kkkkkkkxkFxxGx xxxx=1(4)式 中：()kkx为 有 效 密 度，定 义 为()()001 expkkkx=；kkG为流体之间的相互作用强度，仅考虑最近的邻格粒子的作用力，其控制着两相流体的不溶性和表面张力。流固作用力表示为：()()()()()s,kkksxFxxGx xnxxx=2(5)式中：ns为指示函数，ns=1 时表示固相，ns=0 时表示液相。通过调整流固相互作用强度ksG可以获得不同的接触角。2.2 格子玻尔兹曼方法验证 2.2.1 气泡试验 表面张力是多孔介质中两相流的作用力，反映了液体分子在气液界面上的黏附力，其与多孔介质中的毛细性和润湿性密切相关。在模拟中，采用伪势多相 LBM 通过气泡试验来验证模型的准确性，178 新 能 源 进 展 第 12 卷 将液体设定在气体中心，同时将所有边界设置为周期性边界，液体在计算域的中心进行动态演化，并记录液滴内外压差和液滴达到平衡时的半径。对于悬浮在气体中的液滴，应符合下述拉普拉斯定律，即液滴内外压差与表面张力成正比，与液滴半径成反比：PR=(6)式中：为表面张力；P 为液滴的内外压差；R 为液滴半径。图 2 结果表明在不同的粒子间相互作用力下，液气界面的压差 P 和最终平衡状态下的气泡半径的倒数 1/R 呈线性关系，符合拉普拉斯定律。0.000.020.040.060.080.100.000.010.020.030.04LBM 液滴内外压差结果线性拟合P/(lmlu1ts2)(1/R)/lu1(a)0.000.050.100.15-0.50.00.51.01.52.02.5密度差/(lmlu3)/a.u.kkG(b)图 2 拉普拉斯验证：（a）液滴内外压差与液滴半径倒数之间的关系；（b）不同kkG值的气体最大密度和最小密度之间的差异变化 Fig.2 Laplaces test:(a)the relationship between the pressure difference inside and outside the droplet and the reciprocal radius of the droplet;(b)the difference between the maximum and minimum density of gas with different kkGvalues 2.2.2 静态接触角验证 接触角（）可以表征固体表面的润湿性。当液体在亲水性固体上时，接触角小于 90，少量液体在可湿性固体表面形成薄膜扩散。当液体处于疏水性固体中时，接触角大于 90，液体容易从非湿润固体表面脱落，通过改变 Gks获得流体和固体表面之间不同的湿润性。为了进一步验证模型的可靠性，对模型进行了静态接触角验证。设定液滴位于水平固体板上，并将计算区域的左侧和右侧设置为周期性边界条件，上部和底部设置为无滑移边界条件，观察液滴在不同 Gks下的动态演化过程，记录液固界面的接触角，如图 3 所示，可以观察到通过改变 Gks，接触角在逐渐变化，并且呈现近似线性关系。图 3 显示的结果可用于指定不同的 Gks改变固体表面的湿润特性。上述两个验证表明了本文模型的准确性。-0.4-0.20.00.20.420406080100120140160/()Gks/a.u.图 3 静态接触角随 Gks的变化曲线 Fig.3 Change curve of static contact angle with Gks 3 结果与分析 3.1 聚四氟乙烯双梯度的扩散层 对聚四氟乙烯含量为 10%的双梯度扩散层及其5 次代表性老化过程中液态水流动特性进行研究，同时与常规扩散层的相同老化过程中液态水流动特性进行了对比。图 4 为常规扩散层老化过程中液态水流动过程的剖面图14。图 5 展示了双梯度扩散层老化过程中液态水入侵过程的剖面图。可以看出无论是常规还是双梯度扩散层，液态水在扩散层内都是在毛细压力的作用下进行毛细指进运动，并倾向于入侵较大的和毛细阻力较低的孔隙，随后逐渐形第 2 期 齐丽霞等：基于格子玻尔兹曼方法的梯度化扩散层降解研究 179 成固定的流动通道，最终突破并排出扩散层。后续的液态水也会优先通过该流动通道突破。此外，随着扩散层老化程度的加剧，碳纤维和聚四氟乙烯降解导致扩散层结构破坏，造成局部毛管压力降低，进而减小了液态水入侵阻力，导致液态水饱和度含量随老化的加剧而增加，严重影响扩散层的排水性能。但相对于常规扩散层，双梯度扩散层在老化后，进口段（靠近催化层侧）仍具有较高的聚四氟乙烯含量，因此该处毛管压力依然相对较高，同时液态水突破速度依然相对较快，因此相对于同等老化程度的常规扩散层，双梯度扩散层具有较低的液态水饱和度。图 4 聚四氟乙烯均匀分布的扩散层不同老化程度的两相流14 Fig.4 Two-phase flow of polytetrafluoroethylene in uniformly distributed diffusion layers with different aging degrees14 图 5 聚四氟乙烯双梯度分布的扩散层不同老化程度的气液两相流形态 Fig.5 Two-phase flow of polytetrafluoroethylene with double gradient distribution of diffusion layers with different aging degrees 不同双梯度化扩散层在老化前后液态水饱和度和有效孔隙率变化情况如图 6 所示。当双梯度扩散层靠近进口区域的聚四氟乙烯含量为 14%且出口区域含量为 6%时，如图 5（b）中方案 11 15 所示，相对于常规扩散层和其他双梯度扩散层，其在老化过程中始终具有最低的液态水饱和度含量以及最大的有效孔隙率。在最终老化状态下，其液态水饱和度为 28.1%，与老化前相比，液态水饱和度增加了8.9%，而常规扩散层（图 4）在老化后液态水饱和度增加达 10.2%14。因此，双梯度化扩散层表现出更好的抗老化性能。180 新 能 源 进 展 第 12 卷 初始前期中期后期最终24283236 12%/8%14%/6%16%/4%老化阶段液态水饱和度/%414243444546有效孔隙率/%图 6 双梯度扩散层不同老化程度下液态水饱和度和有效孔隙度曲线 Fig.6 Liquid saturation and effective porosity curves of double gradient diffusion layers with different aging degrees 3.2 聚四氟乙烯三梯度的扩散层 对三梯度扩散层的抗老化性能开展研究，如图 7所示，发现同常规扩散层和双梯度扩散一样，随着老化的加剧，碳纤维和聚四氟乙烯降解导致扩散层结构破坏，造成局部毛管压力降低，进而减小了液态水入侵阻力，液态水入侵区域增多，导致液态水饱和度含量随老化的加剧而增加，严重影响扩散层的排水性能。进一步对比图 7（a）和图 7（c）发现当扩散层进口区域分布较多的聚四氟乙烯时，在最终老化状态下，进口区域液态水突破路径较少，改善了液态水的滞留情况，同时相较于双梯度而言，聚四氟乙烯三梯度分布拥有更平缓的毛管压力分布，更有利于液态水的排出。图 7 聚四氟乙烯三梯度分布的扩散层不同老化程度的气液两相流形态 Fig.7 Two-phase flow of polytetrafluoroethylene diffusion layers with different aging degrees in three gradient distribution图 8 显示了三梯度扩散层在老化前后液态水饱和度含量和有效孔隙率的变化情况。图中可见在相同老化状态下，随着进口区域的聚四氟乙烯含量的增加，扩散层内的液态水含量先增加后减小，当扩散层进口区域聚四氟乙烯含量为 14%、中间区域为10%和出口区域为 6%时图 7（b），该扩散层在老化前后具有最低的液态水饱和度含量和最大的有效孔隙率。其最终老化状态下液态水饱和度含量为27.9%，与老化前的扩散层相比，仅仅增加了 8.6%，因此该三梯度化扩散层具有较好的抗老化性能及气体扩散能力。第 2 期 齐丽霞等：基于格子玻尔兹曼方法的梯度化扩散层降解研究 181 初初前前中前后前最最24283236 12%/10%/8%14%/10%/6%16%/10%/4%老化阶段液态水饱和度/%40414243444546有效孔隙率/%图 8 三梯度扩散层不同老化程度下的液态水饱和度和有效孔隙度曲线 Fig.8 Liquid saturation and effective porosity curves of three gradient diffusion layers with different aging degrees 4 结 论 采用随机算法重建了聚四氟乙烯总含量为 10%的梯度化气体扩散层的多孔结构，其结构与扫描电镜下观察到的结构一致。随后，利用多相伪势格子玻尔兹曼方法研究了梯度化扩散层在老化过程中液态水的输运特性。模拟结果表明，扩散层的老化对液态水的传输行为有很大影响。梯度化设计的扩散层在进口区域分布较高的聚四氟乙烯能有效提供较高的毛管压力，从而减少液态水突破路径，进而提高突破速度。因此，梯度化扩散层具有较好的排水性能和气体扩散层能力。虽然老化后的扩散层因结构破坏导致毛管压力的降低，致使液态水占据更多孔隙而堵塞气体通道，但梯度化扩散层排水性能下降幅度相较于传统常规聚四氟乙烯分布的扩散层老化更低，表现出更高的抗老化性能。参考文献：1 HE W,ZHANG J F,LI H L,et al.Optimal thermal management of server cooling system based cooling tower under different ambient temperaturesJ.Applied thermal engineering,2022,207:118176.DOI:10.1016/j.applthermaleng.2022.118176.2 WANG Y L,XU H K,WANG X D,et al.Multi-sub-inlets at cathode flow-field plate for current density homogenization and enhancement of PEM fuel cells in low relative humidityJ.Energy conversion and management,2022,252:115069.DOI:10.1016/j.enconman.2021.115069.3 WANG Y L,WANG X A,FAN Y Z,et al.Numerical investigation of tapered flow field configurations for enhanced polymer electrolyte membrane fuel cell performanceJ.Applied energy,2022,306:118021.DOI:10.1016/j.apenergy.2021.118021.4 叶跃坤,池滨,江世杰,等.质子交换膜燃料电池膜电极耐久性的提升J.化学进展,2019,31(12):1637-1652.DOI:10.7536/PC190105.5 程植源,周荣良,李嘉颀,等.气体扩散层孔隙率梯度对质子交换膜燃料电池水管理的影响J.内燃机与动力装置,2022,39(3):41-47.DOI:10.19471/ki.1673-6397.2022.03.007.6 WANG Y L,XU H K,ZHANG Z,et al.Lattice Boltzmann simulation of a gas diffusion layer with a gradient polytetrafluoroethylene distribution for a proton exchange membrane fuel cellJ.Applied energy,2022,320:119248.DOI:10.1016/j.apenergy.2022.119248.7 MOLAEIMANESH G,AKBARI M H.Water droplet dynamic behavior during removal from a proton exchange membrane fuel cell gas diffusion layer by lattice-Boltzmann methodJ.Korean journal of chemical engineering,2014,31(4):598-610.DOI:10.1007/s11814-013-0282-6.8 YU S C,LI X J,LI J,et al.Study on hydrophobicity degradation of gas diffusion layer in proton exchange membrane fuel cellsJ.Energy conversion and management,2013,76:301-306.DOI:10.1016/j.enconman.2013.07.034.9 SHAN X W,CHEN H D.Lattice Boltzmann model for simulating flows with multiple phases and componentsJ.Physical review E,1993,47(3):1815-1819.DOI:10.1103/PhysRevE.47.1815.10 DAINO M M,KANDLIKAR S G.3D phase-differentiated GDL microstructure generation with binder and PTFE distributionsJ.International journal of hydrogen energy,2012,37(6):5180-5189.DOI:10.1016/j.ijhydene.2011.12.050.11 SCHMITTINGER W,VAHIDI A.A review of the main parameters influencing long-term performance and durability of PEM fuel cellsJ.Journal of power sources,2008,180(1):1-14.DOI:10.1016/j.jpowsour.2008.01.070.12 刘高洁,郭照立,施保昌.多孔介质中流体流动及扩散的耦合格子 Boltzmann 模型J.物理学报,2016,65(1):014702.DOI:10.7498/aps.65.014702.13 CHEN W,JIANG F M.Impact of PTFE content and distribution on liquidgas flow in PEMFC carbon paper gas distribution layer:3D lattice Boltzmann simulationsJ.International journal of hydrogen energy,2016,41(20):8550-8562.DOI:10.1016/j.ijhydene.2016.02.159.14 WANG Y L,XU H K,HE W,et al.Lattice Boltzmann simulation of the structural degradation of a gas diffusion layer for a proton exchange membrane fuel cellJ.Journal of power sources,2023,556:232452.DOI:10.1016/j.jpowsour.2022.232452.作者简介：齐丽霞（2000-），女，硕士研究生，主要从事质子交换膜燃料电池扩散层及催化层性能研究。王誉霖（1985-），男，博士，副教授，硕士生导师，主要从事新能源汽车燃料电池技术、多孔介质传热传质及两相流、动力电池热管理研究。