论多孔介质中流体流动问题数值模拟方法.pdfVIP

( ) 2000 年 　第 24 卷 　　　　　　　　　石油大学学报 自然科学版 　　　　　　　　　　Vol . 24 　No . 5 　第 5 期 　　　　　　　　Journal of t he University of Petroleum , China 　　　　　　　　Oct . 2000 　　文章编号 (2000) 050 11106 论多孔介质中流体流动问题的 ·综述 · 数值模拟方法 李 　亨 , 　张锡文 , 　何 　枫 (清华大学工程力学系 ,北京 100084) 　　摘要 :简要介绍了多孔介质中流动问题的数值模拟方法及其发展现状 , 对各种方法的优劣及适用范围进行了评 价 。对多孔介质 纯流体耦合流动区域的数值模拟方法及发展现状进行了综合对比。多孔介质模型方法可以用于多 孔介质中流体流动问题的数值模拟 ,并且非常适于进行多孔介质 纯流体耦合流动区域的数值模拟 。 关键词 :多孔介质 ;流体流动 ;多孔介质 纯流体耦合流动 ;数值模拟 ;发展状况 中图分类号 :O 357 . 3 　　　文献标识码 :A 引　言 对多孔介质中流体流动的数值模拟研究一般包 括宏观流动的数值模拟和微观流动的数值模拟 。在 通常的工程应用中 ,人们只对流体的整体流动特性 ( ) 如速度分布 、压力分布等 感兴趣 ,所以从 1856 年 Darcy 开创了这方面的研究以来 ,大量的学者进行 了有关多孔介质中流体宏观流动的实验研究和数值 模拟研究[ 1 ] 。基于工程应用的目的 ,这里只对多孔 图 1 　多孔材料典型的表征单元体 介质中宏观流动的数值模拟方法的发展情况做一概 物理量的平均可以定义在体积单元 V m 上, 也 述 。 可以定义在组成体积单元的流体部分 Vf 上 。例如 流动速度可以定义在 V m 上 ,称为 Darcy 速度 V , 也 1 　理论基础 ( 可以定义在 V f 上 ,称为固有平均速度 int rinsic av 　　虽然多孔介质中基于孔尺度的流体流动是非常 erage velocity ) v 。这两种速度之 间存在 如 下 的 不规则的 ,但包含足够多孔的空间的平均物理量却 Dup uit ─Forchheimer 关系 : φ ( ) 以一定的规律随空间和时间变化 。在对宏观流动的 　V = v , 1 ( ) 其中, φ为多孔介质的空隙率 。 研究中 ,物理量被定义为表征单元体积 R EV 中物 理量的平均值 。R EV 的选取一是应使其尺度远大 一般情况下, 研究多孔介质的宏观流动时, 都采 于孔的尺度 , 以保证平均结果与 R EV 的大小无关 ; 用 Darcy 速度 V 。 二是应远小于宏观流动区域的尺度 , 以保证能够反 2 　基于传统多孔介质流动控制方程的 映出宏观流动的变化规律 。典型的表征单元体积如