期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 466, 期 -, 页码 -出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2022.111396
关键词
Discrete fracture matrix; Computational geometry; Mimetic finite difference; High performance computing; Flow and transport in fractured media; Delaunay triangulation
资金
- LANL LDRD program [20180621ECR, 20220019DR]
- Center for Nonlinear Studies
- Department of Energy (DOE) Basic Energy Sciences program [LANLE3W1]
- National Nuclear Security Administration of U.S. Department of Energy [89233218CNA000001]
- LAUR [LA-UR-21-31458]
This paper presents a comprehensive workflow for simulating single-phase flow and transport in fractured porous media. The workflow includes mesh generation, discretization of governing equations, and implementation of numerical methods for high-performance computing. The paper also provides methods to improve mesh quality and verification tests based on analytic solutions.
We present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据