气体在微细多孔介质中流动阻力研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号:083362 气体在微细多孔介质中的流动阻力研究1 黄寓理 姜培学 胥蕊娜 （热科学与动力工程教育部重点实验室，清华大学热能工程系，北京 100084 ） Tel ：010 Email: Jiangpx@tsinghua.edu.cn 摘要：本文对空气、氢气和氦气流过微细多孔介质内部的流动阻力特性进行了实验研究，分析不同气 体在不同颗粒直径条件下摩擦因子与等效雷诺数的关系。结果表明：对于空气、氢气和氦气，当颗粒 直径为 225 μm 和 125 μm 时，由于可压缩性的影响，摩擦因子略大于经验公式结果。对于空气，当颗 粒直径为 90 μm 和 40 μm 时，实验得到的摩擦因子小于经验公式结果。对于氢气和氦气，当颗粒直径 为 125 μm 时，经过可压缩性修正的摩擦因子实验值略小于经验公式结果；颗粒直径为 90 μm 和 40 μ m 时，摩擦因子小于经验公式结果。通过实验数据得到的渗透率明显大于经验公式计算结果。因此， 对于空气、氢气和氦气，当颗粒平均直径小于 90 μm 时，微细多孔介质中的流动特性存在稀薄气体效 应。 关键词： 微细多孔介质；实验研究；摩擦因子 1 引言 流体在微细多孔介质中的流动研究具有重要的科学价值和实际意义，在发散冷却、 储氢、燃料电池、催化作用、吸附作用、石油开采、危险品处理及微型推进器等领域有 广泛的应用。对于微细多孔介质中的流动已经有非常多的研究，Scheidegger[1]对早期理 论和研究进行了详细的综述。近年来，很多学者研究并改进了气体在多孔介质内的输运 理论[2-6] ，由于气体在多孔介质内的流动受可压缩性的影响较为明显，因此表现出与液 体不同的流动特性，因此经典 Darcy 定律需要被修正。一些研究发现，在同一多孔介质 中，通过实验获得的气体渗透率要大于液体的渗透率。同时，不同气体通过实验得到的 渗透率也不相同[1, 5, 7] 。Klinkenberg[7]首先研究了这种现象并且提出了对于 Darcy 定律的 经验修正。根据 Klinkenberg 的研究, 渗透率的增大是由于气体流过多孔介质壁面时出现 速度滑移所引起的。定义为气体分子平均自由程与孔隙直径之比的 Knudsen 数 Kn ，将 流动划分成四个区域：连续介质区（Kn0.001 ），温度跳跃与速度滑移区（0.001 Kn 0.1 ）, 过渡区（0.1 Kn 10 ）和自由分子区（Kn 10 ）。其中温度跳跃与速度滑移区中在连续 方程和能量方程的基础上，壁面处应使用适当的速度滑移和温度跳跃边界条件进行修正 [8] 。胥蕊娜等在已有经验关联式基础上使用 Kn 数对微多孔介质内空气滑移流动进行了 修正，修正后的关联式与实验结果及数值计算结果符合较好[9] 。但是，不同气体在微细 多孔介质中流动特性差异的比较，以及气体性质对微尺度效应的影响，在他们的研究中 都没有涉及。 本文实验研究了空气、氢气及氦气在微细多孔介质中的流动阻力特性及气体性质、 1 本文得到国家自然科学基金重点项目和973 项目(2007CB210107)资助 颗粒尺度对于流动的影响。实验段为锡青铜颗粒烧结多孔管，颗粒平均直径分别为 225 μm 、125 μm 、90 μm 和 40 μm 。 2 空气，氢气及氦气流过微细多孔介质实验研究 2.1 实验系统及数据处理 气体流过烧结微细多孔结构流动阻力的实验系统如图 1所示。气体从高压气瓶流出， 经过质量流量控制器控制流量，质量流量计测量流量后，进入实验段。对于氢气的实验， 需要在实验段出口处设置气体收集装置。 图 1. 实验系统图 实验段为烧结多孔管，烧结多孔结构孔隙率，颗粒直径见表 1。烧结多孔介质的水 力当量直径用 d [9] =2d ε/3(1-ε)求出 。表 1 还给出了使用该关系式计算的不同颗粒直径下 h p 孔隙当量直径，以及不同气体在多孔管内流动的平