纳微米级孔隙多孔介质内气体流动特征.docxVIP

纳微米级孔隙多孔介质内气体流动特征 多孔介质内液体流动的研究已成为许多应用科学和技术领域的基础，如土壤力学、地下水水文学、石油工程、水处理、工业过滤、陶瓷工程、冶金和防毒剂（气体）。因此，对多孔介质中液体流动的详细特征具有重要的实际意义。一般来说，为了模拟流量流动，假定采用连续假设或分子假设。连续假设适用于大多数流量变化。然而，随着系统长度的减小，例如纳米级孔，间隙通道的几何尺寸非常小，几乎可以与气体分子的平均自由方程相比较。此时，kin数kn（分子自由方程和特征长度的定义）相对较大。当气体通过多孔介质流动时，分子和孔壁之间的碰撞非常严重，壁面发生了滑移，对气体的传输有很大影响。因此，有必要考虑这种纳米规模的进给速度，并不固定连续性假设。 物理学中对流体的描述分为三个层次: (1) 分子层次, 立足分子动力学, 研究分子碰撞的微观机理; (2) 动力学层次 (介观层次) , 用非平衡态统计物理的方法研究流体, 经典的动力学方程是玻尔兹曼方程; (3) 流体力学层次, 是宏观层次, 所用的基本方程是Navier-Stokes方程.目前, 多孔介质内的流体流动描述一般采用第3个层次.然而对内部孔隙结构复杂的纳微米级孔隙多孔介质, 其孔隙尺度主要分布在纳微米级, 很难用以上三种流动层次来描述流动特征. 2001年, Karniadakis和Beskok得出了在连续介质、滑移、对流和不同分子类型下的渗透率变化表达式, 得到了普遍适用于连续流区、滑移流区、过渡流区和自由分子流区的理想气体流动方程.2007年, Florence等对Beskok-Karniadakis模型进行了扩展, 提出纳微米级孔隙多孔介质渗透率变化只与克努森数有关, 并给出了渗透率受克努森数影响的关系式, 但此关系式具有很强的非线性特征.本文基于Beskok-Karniadakis模型, 对纳微米级孔隙内的气体流动方程进行了改进, 提出多项式修正系数, 对Beskok-Karniadakis模型进行简化和推广, 建立了适用于不同流态的纳微米级孔隙内气体流动数学模型, 并应用此模型推导出了纳微米级孔隙多孔介质内的气体流量公式. 1 纳米级多孔气体流动数学模型 1.1 纳米级孔隙介质流态 2002年, Civan研究指出气体在微孔介质中的流动状态取决于多孔介质本身的物理性质和气体分子平均自由程, 并归纳总结Liepmann、Stahl、Kaviany等的研究成果, 提出利用克努森数划分气体流动区域, 把多孔介质内的气体流动分为三个区域: (1) 克努森流区域; (2) 过渡或滑移流区域; (3) 黏滞流区域. 克努森数Kn是一个被广泛认可的量纲为一的参数, 用来判断流体是否适合连续性假设, 其定义为 式中:λ为气体分子平均自由程, m;r为孔喉半径, m. 式中:KB为玻尔兹曼常数, 1.3805×10-23J·K-1;δ为气体分子的碰撞直径;T为温度;p为压力. 表1为不同气体组分的分子碰撞有效直径, 由表1及式 (2) 可以得到气体分子平均自由程与压力和温度的关系, 见图1. 图2为不同尺度下的克努森数随压力变化的关系.由图2可以看出, 随着压力增大, 在双对数坐标系, 克努森数呈线性减小.根据克努森数数值的大小, 可将多孔介质内的流动划分为不同的流态, 即为连续流 (黏滞流) 、滑移流、过渡流及自由分子流 (克努森流) .当Kn≤0.001时, 流动为连续流, 无滑移边界, 处于连续流区, 可用传统的达西定律来描述;当0.001Kn≤0.1时, 流动为滑移流, 存在滑移边界, 处于滑移区;当0.1Kn≤10时, 流动为过渡流, 存在滑移效应, 处于过渡区;当Kn10时, 流动为自由分子流, 处于自由分子区, 用菲克定律描述.对纳微米级孔隙多孔介质, 在压力为0.1~100 MPa范围内, 克努森数小于10, 因此本文所研究流态不涉及自由分子流.微米级孔隙内的流动为连续流和滑移流, 当孔隙直径大于50μm时, 流体流动均为连续流动, 压力增高使得部分滑移流转换为连续流;纳米级孔隙内的流动多以过渡流、滑移流为主, 压力增高使得部分转换为滑移流或连续流. 1.2 利用besk溶渗-karniadakis模型原理介绍 流体在多孔介质内的流动问题通常可以用达西定律进行描述.然而, 达西定律是一个基于宏观观测的实验定律.对于纳微米级孔隙多孔介质, 其内部孔隙结构复杂, 气体流动过程中与纳微米孔道表面发生剧烈碰撞、扩散作用, 同时由于滑移等现象的存在, 不能用简单的达西方程所描述. Beskok-Karniadakis模型得出了在连续介质、滑移、对流和不同分子类型下的渗透率的变化, 从而得到渗流速度为 式中:K0为多孔介质渗透率;μ为气体黏度;x为两个渗流截面间的距离;α为稀疏