268驻极体过滤器气相流场及初始压降的数值计算.doc

驻极体过滤器气相流场及初始压降的数值计算 安徽工业大学 建筑工程学院 钱付平 韩云龙 (陈光 黄志甲 摘要：基于多孔介质模型对二维驻极体过滤器内气相流场进行了数值计算，得到了不同结构（填充密度和过滤器深度）驻极体过滤器内的静压分布。为了观察多孔介质模型在研究过滤器流场的适应性，本文还计算了不包括多孔介质模型时，相同填充密度不同纤维排列驻极体过滤器内的气相流场。结果表明：基于多孔介质模型时，静压沿着过滤器深度方向线性变化，填充密度越大，初始压降越大。因此，如果多孔介质的粘性阻力系数选择合理，该模型可以应用于计算过滤器的气相流场。但是，该模型并没有考虑过滤器内部纤维排列方式对其内部流场的影响，而实际上，过滤器内部纤维的排列行为将显著影响过滤器的流场和初始压降。因此如果要进一步研究过滤器内复杂的气固两相流，还需对该模型进行深入研究。 关键词：驻极体过滤器; 多孔介质模型; 气相流场; 初始压降; 数值模拟 1 前言 据有关部门监测显示，中国48.1%的城市空气质量处于中度或重度污染，颗粒物仍是影响中国空气质量的首要污染物。目前中国大多数城市人口长期生活在可吸入颗粒物超标的空气环境中，人口超过百万的大型城市，空气中二氧化硫和颗粒物超标比例高，空气质量达标比例低。现代社会中，人们几乎90%的时间是在室内度过的[为了提高室内空气品质，减少室内污染物水平，目前普遍采用的一种方式就是引入更多的室外新鲜空气。颗粒物是影响室内的重要因素之一。因此，国外早在二十多年前就开始了对颗粒物的研究，室内颗粒物的浓度及其影响因素也就成了一个重要的研究方向及课题。中华人民共和国国家标准洁净厂房设计规范GBJ73-2001[4]中规定，高效空气过滤器的压降达到初始压降的两倍时则应予以更换。所以从这个意义上来看，准确地预测驻极体过滤器的压降对于过滤器的使用是非常重要的。而且，准确地预测驻极体过滤器内气相流场对于进一步研究其内部微细颗粒的运动也具有积极意义。但是通过实验来研究驻极体过滤器结构尺寸对压降的影响，费时费力。另一方面，随着计算机技术的飞速发展，国外已有部分学者利用计算流体力学（CFD）技术研究驻极体过滤器内部气固两相流场[5-10]，CFD计算和实验比较，明显的优点在于可以在相对便宜的情况下改变流动和结构参数同时可以获得比较准确的结果。但是，绝大部分研究者在计算时都将过滤器简化成一种规则排列的纤维结构，其结果和实际情况会存在一定的差距。实际上，驻极体过滤器内部是一种多孔介质（Porous medium），而且这些纤维大都为一种随机的排列结构。 本文利用数值计算作为工具并基于多孔介质模型来研究结构尺寸及操作参数和气相流场的关系。在此基础上得出一个能研究驻极体过滤器内压力场和速度场的新方法，该方法的探索将有助于优化驻极体过滤器的结构设计。 2 CFD模型 2.1多孔介质模型 多孔介质模型采用经验公式定义多孔介质上的流动阻力。从本质上说，多孔介质模型就是在动量方程中增加了一个代表动量消耗的源项，。该源项可以模拟多孔介质的作用。由两部分组成：一个是粘性损失项，一个是惯性损失项： （1） 式中是第个（、或方向）动量方程的源项，和为给定矩阵。对于本文的多孔介质模型而言，上述源项可以简化成如下形式： （2） 式中为多孔介质渗透性（permeability）。渗透性和捕集效率一样是反映过滤器性能的重要参数。 在计算多孔介质模型是，有两个参数需要在计算前设定，即粘性阻力系数（，）和惯性阻力系数（）。针对流体流过纤维过滤器的情况，可以应用Blake-Kozeny方程[11]计算过滤器内的有关多孔介质参数。Blake-Kozeny方程实际上是半经验Ergun方程的一种简化形式，半经验Ergun方程如下所示： （3） 上式中，为纤维直径，为过滤器的孔隙率，且，其中为填充密度或称固体体积份额（solid volume fraction，SVF）。 对于高效过滤器绝大部分都是工作在层流条件下，因此上式可以简化为如下的形式： （4） 由（4）式可知，各方向上的粘性阻力系数为： （5） Gebart[12]则通过二维的数值模拟认为纤维过滤器中的这种粘性阻力系数实际上包括平行于纤维的粘性阻力系数和垂直于纤维的粘性阻力系数，即： ， （6） Davies[13]通过实验得出了有关初始压降的关联式，已有研究者证明该关联式对于SVF在0.6%-30%范围内是准确的： （7） 式中为过滤介