Consistency and local conservation in finite-element simulation of flow in porous media

被引:0
作者
Laurent White
Jeremy Brandman
Dimitar Trenev
机构
[1] Corporate Strategic Research,ExxonMobil Research and Engineering
来源
Computational Geosciences | 2021年 / 25卷
关键词
Finite-element method; Flow in porous media; Mass conservation;
D O I
暂无
中图分类号
学科分类号
摘要
Computational models of subsurface flow involve the coupling between fluid velocity and transport equations. While solving these problems independently can be challenging, we expose the additional complication associated with coupling them. We focus on numerical methods that can readily accommodate unstructured meshes, with emphasis on the mixed finite-element method for the velocity-pressure formulation. The transport equation is discretized using continuous and discontinuous finite elements, and also finite volumes. We theoretically identify the discrete requirements, in space and time, to be satisfied for the constant-preserving test to pass when coupling these two discrete entities: a unit tracer concentration should be discretely preserved to machine precision under the right conditions. In our work, the constant-preserving test is referred to as the consistency test. In particular, we show that a discrete formulation that relies entirely on continuous finite elements is consistent. The overarching conclusion is that consistency is achieved when the discrete tracer equation reduces to the discrete pressure equation upon setting the tracer concentration to one. Equivalently, we conclude that consistency is achieved when coupling discrete formulations that rely on identical local-conservation properties. In explaining this result, we also make yet another attempt at debunking the myth that continuous finite elements are not locally conservative and cannot be used for flow simulation. We emphasize that continuous finite elements are locally conservative in a way that is derived from the discrete form, but not in a finite-volume sense, that is, not by computing the integral of the normal velocity on any element boundary.
引用
收藏
页码:1123 / 1138
页数:15
相关论文
共 189 条
[1]  
Arnold D(2013)Hierarchical benchmark case study for history matching, uncertainty quantification and reservoir characterisation Computers & Geosciences 50 4-15
[2]  
Demyanov V(2015)Multiphase modeling of geologic carbon sequestration in saline aquifers Ground Water 53 362-377
[3]  
Tatum D(2002)Discrete fluxes and mass balance in finite elements J. Hydraul. Eng. 128 87-92
[4]  
Christie M(2013)Numerical simulations of solute transport in highly heterogeneous formations: A comparison of alternative numerical schemes Adv. Water Resour. 52 178-189
[5]  
Rojas T(1982)Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations Comput. Methods Appl. Mech. Eng. 32 199-259
[6]  
Geiger S(2013)On the application of the continuous Galerkin finite element method for conservation problems SIAM J. Sci. Comput. 35 A2953-A2975
[7]  
Corbett P(2016)Design and analysis of ADER-type schemes for model advection-diffusion-reaction equations J. Comput. Phys. 327 553-575
[8]  
Bandilla KW(2009)Parctical modeling appraoches for geological storage of carbon dioxide Ground Water 47 627-638
[9]  
Celia MA(1984)Discontinuous upwinding and mixed finite elements for two-phase flows in reservoir simulation Comput. Methods Appl. Mech. Eng. 47 93-118
[10]  
Birkholzer JT(1998)The local discontinuous Galerkin method for time-dependent convection-diffusion systems SIAM J. Numer. Anal. 35 2440-2463