燃料电池流道内液滴流动建模与分析.docVIP

PAGE 1 　 燃料电池流道内液滴流动建模与分析 摘要：在燃料电池运行过程中，流道内液态水不能迅速排出将导致电堆进出口压差增大，性能下降，甚至失效。基于液滴存在于通道内造成的流速变化，建立了液滴流动的力平衡模型，推导了液态水累积极限半径。针对液滴累计过程的半径变化周期，对流道的进出口压差建立积分平均模型，获得液滴存在导致的通道进出口压差增长预测。采用两相流仿真验证模型合理性，通过模型分析了关键参数的影响。结果表明，由积水带来的压差提升效应在更深的流道中相对明显，随着流速和接触角的增加则呈先增大后减小的趋势。 　　质子交换膜燃料电池(PEMFC)由于具备零排放、低运行温度、低噪声、高功率密度等优点，有望成为未来的理想能源供给装置。近十年来，PEMFC越来越受到世界各国的重视，并在如固定电力系统、汽车和其他发电机组等方面获得了广泛应用[1-2]。由于PEMFC中电化学反应的主要产物为水，液态水的累积会对电堆性能造成覆盖气体进入催化层界面、影响反应平衡导致反应速率下降等不良影响[3-4]。同时，流道内的液态水会提高流道的流阻，从而使压力分布不均或增加进气空压机的能耗，降低电堆输出效率[5]。所以，电堆中的液态水的累积效应及其排出流动规律就成为国内外专家学者关注的重点。 　　Minor等[6]采用聚二甲基硅氧烷(PDMS)、玻璃和碳纸制造矩形气体通道实验模型，观察PEMFC气体分配通道中液滴的流动行为。Ding等[7]通过局部通道的两相流模型模拟了在微通道中从气体扩散层(GDL)表面出现水滴的动态形成过程。Chen等[8]考虑了微观结构对GDL表面的影响，将复杂的GDL结构描述为分散在顶部空白空间中的立方体阻块进行了仿真研究。Hossain等[9]研究了GDL孔径和孔隙分布对液态水流动的影响。 　　上述研究都集中在液态水的流动特性，但是气体通道内的液态水流动行为与通道的几何尺寸具有很强的相关性，液态水的形成会对通道内压力产生重要的影响，进而影响液态水的排除。现有的理论解析模型难以准确地预测流道中液态水影响下整个通道的气体压力变化情况。而直接采用多相流计算机仿真的方法由于受限于庞大的计算资源需求和收敛稳定性问题，对于燃料电池电堆的应用仍有难度。 　　本文基于液滴附着受力平衡条件，建立累积失稳判断条件，分析流道内水滴累积的极限半径，并基于此累积过程对流道存在液态水情况下的流阻变化情况进行了建模分析，并讨论流道高度、吹扫速度和壁面物性等关键参数对流道进出气压差的影响。 　　1理论模型 　　1.1液滴累积平衡条件 　　对于流道内的液态水，其粘附壁面累积的原因为表面张力，而排出的动力来源为反应气的吹扫作用。如图1所示，液滴平衡条件为所受空气阻力与表面张力的平衡： 　　由于液滴一旦接触到更亲水的双极板表面，其横截面内的力平衡将被破坏，表面张力作用会使液滴向极板运动。液滴一旦运动起来，便不会持续阻碍反应气传输且更易于离开气体通道，此时静力平衡描述将失效。故估计特定流速下形成的最大液滴尺寸，应考虑在GDL表面的情况。考虑液滴累积到极限的情况满足如下条件： 　　(1)液滴在气流方向上变形使得前进角达到极大和后退角达到极小； 　　(2)液滴仍未开始运动，位于GDL表面。 　　当液滴在GDL表面达到极限半径，表面张力不能再维持液滴附着的情况。将表面张力的粘附作用放大到极限条件为： 　　1.2附加阻力分析 　　由于液态水的存在，液滴受到的空气阻力体现为对液滴的推动作用。而与此同时，通道内的气体也受到反作用力，即为通道上覆盖液态水表面的摩擦作用。而其他未受液态水影响区域壁面的摩擦作用不变。 　　未考虑液态水情况的局部压差变化可以通过式(12)计算： 　　由于实际情况下液滴往往会在某时刻接触到侧壁，从而不能保持完整形状，此时其迎风面积会小于理论值，对应的阻力计算也应有所衰减。故定义了式(20)中的l为阻力衰减系数，用于反应实际气体通道内液滴的平均阻力变化情况，取值范围为0～1。 　　2模型验证 　　为了验证模型，通过商业CFD软件Fluent进行两相流仿真，将得到的压差结果进行对标。在此采用几何较复杂的蜿蜒流道验证模型的有效性。单个周期蜿蜒流道的几何参数设定和两相流仿真结果如图2所示。仿真采用的空气流入速度为5m/s，液态水通过注水孔模拟进入通道的速度为0.25m/s；在固定流道宽度为1mm情况下，针对图中的几何参数，设计了3因素3水平的正交实验，具体的实验参数与结果对比如表1所示。通过CFD仿真与模型计算的结果对比可以发现，模型可以准确反应几何条件变化时压差的改变趋势，大部分情况的预测偏差低于10%。造成结果偏大的原因为对接触前进/后退角的缩放；而在流道高度较高时，对形成完整液滴的空间充分，阻力衰减系数的引入会导致模型预测结果偏小。