多孔介质孔隙尺度下不可压缩流体流动特性SPH模拟.pdfVIP

第53 卷第2期 大连理工大学学报 Vol. 53 , :So. 2 2013 年 3 月 Journal of Dalian University of Technology Mar. 2 0 1 3 文章编号: 1000-8608(2013)02-0189-05 多孔介质子L隙尺度下不可压缩流体流动特性SPH 模拟 李维仲\赵月帅，宋永臣 (大连理工大学海洋能源利用与节能教育部重点实验室，辽宁大连 116024 ) 摘要:基于孔隙尺度下的多孔介质模型，使用光滑粒子动力学(SPH数值模拟方法研究多 孔介质中不可压缩流体流动现象，分析流体在孔隙尺度下的流动特性.首先通过模拟经典的 Poiseuille 流和Couette 流验证方法的正确性.随后，针对流体在分别由圆柱和方柱为骨架构 造的规则多孔介质内的流动情况，同时使用 SPH 方法和有限无法(FEM) 进行模拟，所得的 结果十分吻合.最后利用 SPH 方法对流体在基于贝雷岩心获得的多孔介质模型和以同尺度 的圆柱为基元构造的一种人工随机多孔介质模型内的流动进行模拟，证明平均流体速度和单 位质量体积力具有很好的线性关系. 关键词:多孔介质;孔隙尺度;光滑粒子动力学{SPH) 中图分类号: 035 1. 2 文献标志码:A 的微观特性，又能再现流体的宏观流动特性.在处 。引言 理复杂边界问题时.LBM 具有极大的优势，因此 在自然界和很多工业应用中，都包含了多孔 它是一种可以用来模拟具有任意几何边界系统 介质内流体流动问题，如石油开采、地下水净化和 (例如多孔介质〉的非常有效的模拟方法.与有限 流化床等领域.因此，深入了解在不同工况下，流 差分方法等相比.LBM 算法简单、易于程序化处 体在多孔介质中的流动特性有着重大的实际意 理、计算效率高和具有天然的并行性，可以在计算 义.在实际工程中，流体在多孔介质内流动问题可 机或者集群上运行，减少了计算代价和计算时间. 以用达西定律进行描述，然而达西定律是一个基 但是.LBM 仍然是一个基于网格的方法，在多尺 于宏观观测的实验定律，流体和多孔介质之间复 度问题中，如果计算区域复杂，尤其是在解决具有 杂的相互作用只能体现在一个宏观的统计参数 h 大畸变、运动物质交界面、变形边界以及自由表面 之中.因而，要详细研究这个参数，了解多孔介质 问题时，容易产生网格扭曲与网格重构问题. 具体构造对h 的影响，发展一种新的能够在孔隙 一个可以替代 LBM 的方法是光滑粒子动力 尺度下直接高效地模拟流体在复杂多孔介质中流 学(SPH)方法[IJ. SPH 方法是一种元网格粒子方 动细节的数值模拟方法是非常必要的. 法，最初主要用于解决天文学问题口-3J