PEMEC压缩阴极气体扩散层内两相传输特性的孔隙网络模拟.pdf

江苏大學学报（自然科学版）JOURNAL OF JIANGSU UNIVERSITY(Natural Science Edition)D0I:10.3969/j.issn.1671 7775.2023.06.0052023年11月第44卷第6 期Nov.2023Vol.44No.6开放科学（资源服务）标识码（OSID）：品PEMEC压缩阴极气体扩散层内两相传输特性的孔隙网络模拟陈黎12，罗马吉12，刘（1.武汉理工大学现代汽车零部件技术湖北省重点实验室，湖北武汉430 0 7 0；2.武汉理工大学燃料电池湖北省重点实验室，湖北武汉430 0 7 0)摘要：采用三维孔隙网络模拟方法对Toray090型碳纸随机重构，分析了不同压缩率下质子交换膜电解池阴极气体扩散层（gasdiffusionlayer，G D L）的氢气突破压力、毛细压力分布、相对渗透率和相对扩散率随水饱和度的变化，并对氢气在阴极CDL中传输的传质关系式进行了研究.结果表明：根据孔隙网络模拟计算得到的数据重新拟合氢气传质关系式更适合描述氢气在阴极CDL中的传输；不同压缩率下，氢气在CDL中的传质性能参数（相对渗透率和相对扩散率）几乎不变；随着压缩率的增大，氢气在CDL中传输的阻力也随之增大，导致产生的氢气在质子交换膜界面的压力增大，氢气更容易从阴极跨膜渗漏到阳极.关键词：质子交换膜电解池；孔隙网络模型；气体扩散层；压缩；两相流；氢气传输中图分类号：U473.4；T Q 116.2引文格式：陈黎,罗马吉，刘成PEMEC压缩阴极气体扩散层内两相传输特性的孔隙网络模拟 J.江苏大学学报（自然科学版),2 0 2 3,44(6):6 51-6 56.Pore network simulation of two-phase transport characteristics in(1.Hubei Key Laboratory of Advanced Technology for Automotive Components,Wuhan University of Technology,Wuhan,Hubei 430070,China;2.Hubei Key Laboratory of Fuel Cell,Wuhan University of Technology,Wuhan,Hubei 430070,China)Abstract:Based on Toray 090 carbon paper,the gas diffusion layers in the proton exchange membraneelectrolysis cell under different compression levels were randomly reconstructed,and the breakthroughpressure,capillary pressure curve,relative permeability and relative diffusivity of hydrogen in gasdiffusion layer(GDL)were calculated by the three-dimensional pore network model.The mass transferrelationship of hydrogen in the cathode GDL was also investigated.The results show that the hydrogenmass transfer relation reconstructed from the pore network simulation data is more suitable for describingthe hydrogen transport in the cathode GDL.The mass transfer performance parameters of relativepermeability and relative diffusivity of hydrogen in GDL are almost unchanged under different compressionlevels.With the increasing of compression level,the resistance of hydrogen transmission in GDL isincreased,which leads to the increase of produced hydrogen pressure at the membrane interface,andhydrogen is more likely to leak from the cathode transmembrane to the anode.收稿日期：2 0 2 1-0 8-19基金项目：国家留学基金资助项目（2 0 18 0 6 9 550 7 8）作者简介：陈黎（19 9 7 一），男，江苏南京人，硕士研究生（9 9 548 0 8 6 7 ），主要从事PEM电解池建模仿真研究.罗马吉（19 7 4一），男，江西吉安人，教授（），主要从事氢能与燃料电池研究。刘成1：2文献标志码：A文章编号：16 7 1-7 7 7 5(2 0 2 3)0 6-0 6 51-0 6compressed gas diffusion layer of PEMECCHEN Lil-2,LUO Majil-2,LIU Chengo.1,2652江苏大学学报（自然科学版）第44卷Key words:proton exchange membrane electrolysis cell;pore network model;gas dffusion layer;compression;two-phase flow;hydrogen transport质子交换膜电解池（proton exchange membraneelectrolysiscell，PE M E C)效率高、产物气体纯度高、绿色环保、能耗低、体积小、响应快，是一种将电能转化为化学能的能量转化装置.气体扩散层（gasdiffusion layer,GDL)是PEMEC中的核心组件,可将反应物水输送到反应位点，同时有利于从反应位点除去气体 .为了优化电解池性能,通常电解池在一定的夹持压力下运行.S.TOGHYANI等 2 通过对PEMEC建立三维计算流体动力学（CFD)模型,研究了操作条件和设计参数对PEMEC性能的影响,研究发现在1.6 5V的电压下,GDL的厚度从0.5mm降低到0.2mm,电流密度降低约3%，电解池性能得到提升.S.AL SHAKHSHIR等 3 通过试验研究了活性面积为50 cm的电解池在不同夹持压力下的性能变化,结果表明,PEMEC 的性能随夹持压力的增加而增加电化学阻抗谱（electrochemical impedancespectroscopy,EIS)测量表明,主要是由于高夹持压力下欧姆电阻和活化电阻的降低.K.KISHOR等 4 研究了夹持压力对催化剂活性/孔隙率和GDL/流道之间接触电阻的影响,结果表明,随着夹持压力的增加,接触电阻降低，电化学活性和孔隙率也随之降低.在阴极侧的受压过程中，试验研究发现阴极扩散层的电导率 5 以及阴极侧催化层中离聚物含量的改变导致了PEMEC性能的变化.目前的研究 5 表明，在高夹持压力下，气体扩散层厚度减小确实可以提升电解池的性能，但是其中的作用机制尚未明确，多数研究将其归因于气体扩散层以及催化层电导率的增大，而压缩过程中两相传输特性是否发生变化值得探讨.因此，笔者基于OpenPNM软件包建立PEMEC阴极气体扩散层的三维孔隙网络模型,探究压缩对阴极GDL内水气传输特性的影响，并从介观尺度分析压缩对阴极GDL内传输性能参数的作用机制.1模拟方法孔隙网络模型（porenetworkmodeling，PNM）作为一种模拟3D多孔介质内部多相流流动的高效工具，已经被大量应用于地质科学、燃料电池等领域的研究中 6 .在孔隙网络建模中，将多孔介质视为由喉道连接的孔隙网络，与同样适用于多孔介质建模的格子玻尔兹曼方法（lattice Boltzmann method，LBM)相比,在保证相对精度的情况下,完成相同规模的计算需要的时间更短.孔隙网络模型还可以解析实际的饱和度剖面曲线，在与中子射线的可视化试验对比上也具有优势 6 .1.1几何建模在开源孔隙网络建模软件包OpenPNM中，参照T.G.TRANTER等 7 的方法构建了7 50 m750m375 m的孔隙网络模型,并与试验中测得的Toray090GDL的孔隙率和单项渗透率数据 8 相匹配.孔隙网络几何模型如图1所示，孔径分布以及喉径分布如图2 所示.Z图1Toray090型气体扩散层的孔隙网络几何模型2520%/聘505152535.455556575孔隙直径/um(a)孔径分布14%/聘420图2 Toray090型气体扩散层的孔隙网络模型孔径及喉径分布1.2数学模型在孔隙网络模型中,对于通过两个孔及二者之间喉道的单相流，体积通量q;与ij两孔之间的压20406080100喉道直径/um(6)喉径分布120653第6 期陈黎等：PEMEC压缩阴极气体扩散层内两相传输特性的孔隙网络模拟差成正比,即qij=g;(pi-p,),式中：g,为ij两孔之间的流体流动的传导率;PPi分别为讠i两孔中的压力.在孔隙网络模型中，复杂的喉道横截面被简化为规则形状,喉道中的单相传导率为(2)uL式中：C为无因次常数，其大小取决于喉道横截面形状，当横截面为圆形、四边形和三角形时，分别取0.5000、0.56 32 以及0.6 0 0 0;A为喉道横截面面积;G为形状因子;为流体的黏度;Lo为毛细管长度.在稳态状态下，流入和流出孔喉单元的质量守恒，即4g=Zg(p:-P,)=0.j=1由于孔隙网络模型类似于电阻网络模型，因此，在运用斯托克斯流求解多孔介质的渗透率以及菲克扩散定律求解其扩散率也类似于在电阻网络中运用欧姆定律求解电导率.g，可表示为1111+式中：gp:vgm,为ij两孔的流动传导率;g.为连接ij两孔之间喉道的流动传导率.由式(2）（4)可以迭代求解整个孔隙网络的压力场,从而计算出穿过孔隙网络的体积通量和相应的渗透率。在PEMEC阴极GDL的情形下，氢气相侵人充满液态水的GDL,定义毛细压力P。为两相压力之差，即P.=Pw-PH,式中：PwPH分别为水和氢气的压力.当水侵入充满气相的GDL时需要正毛细压力，而PEMEC运行情况下，则是由氢气侵人液相的GDL,此时毛细压力为负值.1.3多相传输算法在多孔介质中运用孔隙网络模型求解其渗透率和扩散率一般用达西定律和菲克扩散定律来描述，公式如下：Kefr=A,Ap/(QuL),Defr=NAL/(A,c),式中：Kr为有效渗透率;A为CDL人口的截面积；p 为多孔介质两端压差；Q为通过多孔介质的流量;L为多孔介质厚度;Defr为有效扩散率；NA为通(1)过多孔介质的扩散流量；c为两端浓度差。在多孔介质的多相流模拟中,PNM通过求解每个相中的方程来模拟多相流，因此,其中一个相的相对渗透率或扩散率在另一相所占据的孔和喉中会大大降低，通过这种方式，可以在运用渗透算法的同时运行传输算法来确定某一相的相对传输特性.Ker=KKr.,Defr=DDr.p,式中:K为单相渗透率;Kr,为相对渗透率;D为单相扩散率;Dr，,为相对扩散率.Kr,取决于饱和度的大小和饱和度变化的历史，在连续性模型中通常使用经验模型来描述.K.,=Ss,式中：S,为p相的饱和度;=3.0.(3)同样在连续性模型中通常使用经验模型来描述Dr，p，计算公式如下：Dr.,=S,式中:a=2.0 或1.5.实际GDL中相对渗透率或相对扩散率与饱和(4)度的关系往往比经验模型更复杂，为了更好地描述PEMEC中GDL的两相流行为,使用PNM来模拟并拟合Kr,或Dr,与饱和度 S,的关系.2模型验证为验证孔隙网络模型的准确性,对此模型进行注水和排水的仿真研究,将仿真得到的毛细压力和饱和度关系曲线与J.T.GOSTICK等 9 的试验结果进行比较.毛细压力曲线测量系统 9 中水从底部进(5)人，气体则通过顶部进出，由于液体压力保持恒定，需要通过控制气体压力来控制毛细压力，气体压力是通过注射泵来调节,GDL中的饱和度通过分析天平上水的质量变化来测定.当毛细压力达到设定值后，监测天平上的质量变化，当其变化速度足够慢时,则认为此时GDL中的饱和度稳定，记录数据并进行下一个毛细压力点的测量.孔隙网络模型模拟结果与试验结果 9 对比如图3所示，右侧曲线为液态水侵入多孔介质时（进(6)水过程）,计算域中水饱和度随毛细压力的变化;左(7)侧曲线代表空气侵入充满水的多孔介质时也就是水从多孔介质中排出时（排水过程），相应的饱和度变化，(8)(9)(10)(11)654江苏大学学报（自然科学版）第44卷1.0r0.80.60.40.20-30-200-10毛细压力/kPa图3孔隙网络模型模拟结果与试验结果的对比从图3可以看出：第1个区域是毛细管压力突然增加,直到接近水饱和度为0.2,该区域表示由于多孔介质的表观粗糙度，使得表面附近的孔分布与整体孔分布有差异，而试验测量的GDL往往是已经添加了PTFE涂层进行了疏水处理，试验值在低压区域会出现一个肩峰,PNM模拟因为无法施加这一影响而比较平缓;第2 个区域是一个平稳段,在该平稳段中,水饱和度增加到0.8 以上,而毛细管压力几乎保持恒定,这一区域表示多孔材料中平均孔径的孔被水填充;第3区域是毛细管压力的急剧增加,这归因于吸水可用孔空间的减少.由PNM模拟结果与试验结果对比可知:在较大范围的毛细压力内，模拟结果与试验结果具有较好的一致性.表明该孔隙网络模型可靠，3计算结果分析3.1未压缩GDL内两相传输特性规定、y 方向为平行碳纸平面方向即面内方向，z方向为垂直碳纸平面方向即厚度方向.在实际PEMEC运行中，氢气作为介质将在z方向上通过湿润的阴极气体扩散层传递到另一侧,其传输能力的强弱与阴极气体扩散层中水饱和度密切相关.氢气相对渗透率与水饱和度之间的关系曲线如图4所示.1.00.80.2F0图4氢气相对渗透率与水饱和度之间的关系曲线从图4可以看出：在PEMEC的运行过程中，氢气不断产生，氢气的相对渗透率随着多孔介质中水饱和度的降低而不断提高；氢气的相对渗透率在水饱和度值较低时，与经验公式Kr，=S，吻合较好，当水饱和度为0.2 0.7 时,z方向的氢气相对渗透率进水过程文献 9 的试验进水过程PNM模拟一排水过程文献 9 的试验一排水过程PNM模拟10203040一x方向一方向一方向一经验模型式(10),a=3.00.20.4水饱和度值明显高于其经验公式的计算值，且偏差较大，表明在PEMEC运行的大多数时间，实际的氢气渗透量将高于一般经验模型的计算值.将PNM计算得到的z方向氢气相对渗透率与水饱和度关系进行拟合得到新关系式：Kr.H=0.835S%+2.629S2.794S.+0.996,式中：KH为氢气相对渗透率;S为水饱和度，其拟合标准差为0.0 12 0 6,确定系数为0.9 9 9 3，拟合程度更高，这意味着方程（12)更适用于PEMEC中阴极气体扩散层的氢气传输计算.氢气相对扩散率与水饱和度之间的关系曲线如图5所示，在PEMEC运行中氢气在催化层和气体扩散层界面浓度较高处向z方向浓度较低处传递,其传输能力与阴极气体扩散层中水饱和度密切相关.1.00.8F二方向一方向方向0.2一经验模型式(11),a-2.0一经验模型式(11),a=1.500.2图5氢气相对扩散率与水饱和度之间的关系曲线从图5可以看出：在PEMEC刚启动时，水饱和度较高，氢气的相对扩散率较低，氢气难以传输，随着PEMEC的运行，氢气相对扩散率随着水饱和度的降低而大大提高.在z方向，由于气体扩散层的厚度远小于其长宽，气体更易扩散，当水饱和度同比降低时，方向的氢气相对扩散率提高得更快，无论经验公式Kr,=S中a=2.0还是=1.5,都不能体现这一变化.因此PEMEC的氢气相对扩散率与水饱和度的经验公式需改进.通过将PNM所模拟的数据进行拟合，得到z方向氢气的相对扩散率与水饱和度之间的关系式：0.60.8(12)0.40.6水饱和度1.0Dr.n=6.005St 9.912S%+3.893S%-0.964 8.Sw+1.047,式中：D,H为氢气相对扩散率.其拟合标准差为0.0 2 449,确定系数为0.9 9 7 2，表明式(13)能更好描述由浓差驱动的氢气扩散.0.81.0(13)655第6 期陈黎等：PEMEC压缩阴极气体扩散层内两相传输特性的孔隙网络模拟3.2GDL压缩对氢气相对渗透率的影响为研究不同压缩率对阴极GDL两相传输特性的影响,设置了4种不同压缩率的 GDL进行模拟计算,在压缩过程中作如下假设：压缩过程理想化认为是均匀压缩,不考虑压缩过程中纤维的变形.定义压缩率(compress ratio,CR)为GDL压缩厚度变化值与初始厚度值之比,其相关参数如表1所示.表1不同压缩厚度下的GDL参数厚度/um压缩率/%3750338103002026330厚度方向不同压缩率下氢气相对渗透率的变化曲线如图6 所示.1.0元0.80.2F0图6 厚度方向不同压缩率下氢气相对渗透率的变化曲线从图6 可以看出：在低饱和度（0 0.2)以及高饱和度(0.8 1.0)区域,压缩率的变化不会改变厚度方向的氢气相对渗透率；在中等饱和度的区域即在PEMEC运行过程中的大部分情况下，随着厚度方向逐渐进行压缩，氢气的相对渗透率逐步增大.这是因为厚度方向距离的减小作用远大于孔隙率减小的作用,而压缩到一定程度（CR为2 0%以上），压缩对相对渗透率的影响微乎其微.面内方向氢气相对渗透率在不同压缩率下的变化曲线如图7 所示，氢气在不同水饱和度下相对渗透率不会随着压缩率的变化而改变.1.00.80.20图7 面内方向不同压缩率下氢气相对渗透率的变化曲线3.3压缩对氢气相对扩散率的影响不同压缩率下厚度方向和面内方向GDL中氢气相对扩散率随着水饱和度的变化曲线分别如图8、9 所示.1.00.80.6一CR-0-CR=10%CR=20%迁曲度-.CR=30%孔隙率/%0.2F73.552.4071.972.2170.722.1568.431.81CR=0-CR=10%.-CR20%-CR=30%0.20.4水饱和度-CR-0-CR=10%:CR-20%-:-:CR=30%0.20.4水饱和度0图8 厚度方向不同压缩率下氢气相对扩散率的变化曲线1.00.80.2F0图9 面面内方向不同压缩率下氢气相对扩散率的变化曲线0.60.80.60.80.20.21.0从图8、9 可以看出：随着水饱和度的增加，氢气无论是在厚度方向还是在面内方向都难以进行传输,而增大压缩率对氢气的扩散影响甚微.这是因为在压缩过程中,虽然厚度的减小使氢气更容易传输到流道一侧，但是随着压缩率的增大，整体的孔隙率也随之降低,这使得气体传输更加困难,这两者共同作用使得在均匀压缩下氢气在GDL中的扩散不会随着压缩率变化发生太大改变，3.4压缩对GDL内氢气压力的影响不同压缩率下氢气在水中传输时,水饱和度与毛细压力的关系曲线如图10 所示.1.00.80.60.40.20-12-10-8-64-202毛细压力/kPa图10 不同压缩率下，水饱和度与毛细压力的关系曲线1.0从图10 可以看出：随着压缩率的增大，达到某一饱和度所需要的毛细压力绝对值也随之增大，这0.4水饱和度CR-0-CR=10%.CR=20%-CR=30%0.40.6水饱和度CR-0CR-10%CR-20%CR-30%0.60.80.81.01.0656江苏大学学报（自然科学版）第44卷表示在压缩过程中,压缩率增大,GDL 的微结构将变得更加复杂，传质阻力进一步增加，这就需要更大的气相压力排出氢气.随着阴极侧氢气压力的增大，一方面排出扩散层的氢气压力将增大，另一方面会导致催化层与CDL交界面处的压力增大,从而导致界面处更大的氢气浓度，氢气跨膜渗透量也随之增大，从而带来安全隐患.不同孔隙率下阴极GDL中的气相突破压力如图11 所示,压缩率增大,降低了 GDL的孔隙率,GDL结构的复杂程度也随之增大,氢气流出 GDL到达流道需要更大的突破压力，当气体扩散层压缩率从0 增加到30%，孔隙率从0.7 35减少为0.6 8 5突破压力增加了38.3%,其增大代表着阴极侧GDL内整体气相压力的增大，这不仅会导致氢气传质阻力的增大，还会促进氢气的跨膜渗漏过程.8.58.07.57.06.56.05.50.68 0.690.700.710.720.730.74孔隙率图11不同孔隙率下的突破压力曲线4结 论1）孔隙网络模拟计算结果表明：现有的经验模型对于氢气在 GDL内的传输描述并不准确,可以根据孔隙网络模拟计算得到的数据重新拟合获得更准确的氢气传输公式。2）增大压缩率,GDL厚度方向距离减少与结构复杂程度增大的共同作用,使得氢气在GDL内的传质性能参数（相对渗透率和相对扩散率）不会有太大改变.3）压缩率从0 增大到30%，氢气的突破压力增加了38.3%，GDL内的氢气压力显著提升.4）PE M E C在高夹持压力下的运行过程中，GDL内的氢气压力的提升,致使GDL与 CL界面处的氢气压力更高，从而导致界面处更大的氢气浓度，氢气跨膜渗漏量增大等问题，参考文献(References)1 BA R E IB K,D E L A R U A C,M O CK L M,e t a l.L i f e c y-cle assessment of hydrogen from proton exchange mem-brane water electrolysis in future energy systems J.Applied Energy,2019,237:862-872.2 TOGHYANI S,FAKHRADINI S,AFSHARI E,et al.Optimization of operating parameters of a polymer ex-change membrane electrolyzer J.International Journalof Hydrogen Energy,2019,44(13):6403-6414.3AL SHAKHSHIR S,FRENSCH S,KAER S K.On theexperimental investigation of the clamping pressureeffects on the proton exchange membrane water electro-lyser cell performance J.ECS Transactions,2017,77(11):1409-1421.4KISHOR K,PARASHTEKAR A,SAHA S,et al.Ra-tional design strategy for optimization of clamping pres-sure to minimize contact resistance between electrodeand current collector while preserving porosity of elec-trodes in water electrolyzers J.The Canadian Journalof Chemical Engineering,2018,96(4):881-885.5OMRANI R,SHABANI B.Hydrogen crossover in pro-ton exchange membrane electrolysers:the effect of cur-rent density,pressure,temperature,and compressionJ.Electrochimica Acta,D0I:10.1016/j.electacta.2021.138085.6黎方菊，吴伟，汪双凤.PEMFC带沟槽气体扩散层内传输特性孔隙网络模拟 J化工学报，2 0 2 0,7 1(5):1976-1985.LI F J,WU W,WANG S F.Pore network simulation oftransport properties in grooved gas diffusion layer ofPEMFCJ.CIESC Journal,2020,71(5):1976-1985.(in Chinese)7TRANTER T G,GOSTICK J T,BURNS A D,et al.Pore network modeling of compressed fuel cell compo-nents with OpenPNMJ.Fuel Cells,2016,16(4):504-515.8 GOSTICK J T,FOWLER M W,PRITZKER M D,et al.In-plane and through-plane gas permeability of carbon fi-ber electrode backing layersJ.Journal of PowerSources,2006,162(1):228-238.9(GOSTICK J T,IOANNIDIS M A,FOWLER M W,etal.Wettability and capillary behavior of fibrous gas dif-fusion media for polymer electrolyte membrane fuel cellsJ.Journal of Power Sources,2009,194(1):433-444.(责任编辑贾国方）