论文部分内容阅读
在页岩气开采过程中，由于水的存在，流动通道内会出现气水两相同时流动的情况，而对于微纳尺度的流动而言，表面张力会在流动过程中占据主导地位，水的存在会堵塞流动通道，造成“水堵”现象。因此，气水两相界面的流动机理和力学性质对于分析气体在页岩中的流动状态而言至关重要。
基于SPH的页岩微尺度流动模拟
【摘 要】
：
在页岩气开采过程中，由于水的存在，流动通道内会出现气水两相同时流动的情况，而对于微纳尺度的流动而言，表面张力会在流动过程中占据主导地位，水的存在会堵塞流动通道，造成“水堵”现象。因此，气水两相界面的流动机理和力学性质对于分析气体在页岩中的流动状态而言至关重要。
【机 构】
：
中国石油大学石油工程学院油藏工程系,青岛,266580
【出 处】
：
第十届全国流体力学学术会议
【发表日期】
：
2018年11期
其他文献
基于分子动理论的格子Boltzmann方法(Lattice Boltzmann Method,LBM)是通过微观的分子运动力学推导演化而来的,具有清晰的物理背景、并行特性、易于处理边界条件和实现程序实施,以及计算稳定的优点。浸入边界法(Immersed Boundary Method,IBM)对模拟大变形柔性结构与流体相互作用问题具有极大优越性。由于LBM与IBM使用的都是笛卡尔网格系统,本文有效
炸药爆轰反应流动的数值模拟以拉格朗日方法应用最为广泛，特别是交错网格型拉格朗日方法(运动量定义在网格节点、状态量定义在网格单元中心)，但是，交错网格型拉格朗日方法存在如下几个主要问题：1)流动的内能方程便于求解而非总能量方程被求解，因此系统的总能量难以保持守恒；2)人为粘性用于捕捉爆轰波时爆轰波间断往往被抹平或者产生非物理振荡，且人为粘性表达形式及系数靠经验确定；3)时间方向和空间方向的离散均仅为
浸润边界-格子波尔兹曼通量求解方法(IB-LBFS)是一种在笛卡尔网格中计算复杂外形非定常流动问题的新方法。其应用有限体积方法，通过在单元界面进行基于格子波尔兹曼方法(LBM)的通量重构进行流场求解，相对于传统的LBM方法提高了计算效率和灵活性。通过应用隐式边界条件，实现了流场边界条件的精确满足。
本文介绍了作者在多区域流动问题上的一些进展，包括物理建模、数值算法和计算机模拟。首先我们提出了一系列的数值方法来捕捉多区域系统中存在的界面网络的时间演化。这些方法包括一种新颖的界面网络输运追踪格式，一种高阶精度重新初始化算法，和一种快速界面结构重构方法。这些方法具有高效、高精度、鲁棒和数值一致的优点，并且提供了一种与其他物理耦合的框架。
对于越来越稀薄的流动,Navier-Stokes方程逐步失效部分归因于其单温度平衡近似.为了克服这个问题,基于Bhatnagar-Gross-Krook(BGK)方程和Landau-Teller-Jeans弛豫模型的连续多温度模型已经被提出(K.Xu et al,2008,并成功应用于二维非平衡流计算.在最近的研究中,已经成功发展出应用于平衡流动的两步四阶气体GKS(L.Pan et al,201
提出一种基于速度修正型浸入边界法(IBM)的PISO求解器，并采用开源CFD软件OpenFOAM对其进行验证.IBM通过欧拉网格中一系列的拉格朗日点来标识物体边界，而物体边界与流场通过回复力的方式相互影响，即在流体运动方程中加体积力项来实现.传统IBM在速度场的计算中事先求解回复力的值，因此其边界点难以完全满足无滑移边界条件.
本文以Hartmann的人工粘性激波捕捉方法为基础，提出了一种基于通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法，该方法形式简单，可向非均匀网格进行推广，编程易于实现。数值试验表明，对于等熵涡算例和高阶曲网格的圆柱算例，在大于二阶的情况下，该方法能有效保证解的精度，对光滑区的影响很小。
本文提出了一种模拟大密度比多相流的界面格子玻尔兹曼通量求解器(ILBFS)。不同于传统的有限体积多相流求解器,是利用宏观物理量来计算单元边界通量,ILBFS利用Chapman-Enskog展开重构宏观控制方程的对流项和扩散项,然后利用格子玻尔兹曼方程的局部解来计算单元边界通量。
多尺度气体流动在航空航天工程中广泛存在,例如,在微喷管或再入飞行器周围,流场尺度的分布跨越稀薄和连续流动区域。基于气体动理学理论的随机粒子方法,如DSMC,可以准确捕捉稀薄气体效应。然而,在分子碰撞占主导的连续流动区域,DSMC方法的计算效率变得非常低。为了克服DSMC方法的这一困难,可以在近连续区域采用模型方法简化分子碰撞过程,如BGK模型。