A node-centered finite volume method for a fracture model on triangulations

Abstract In this paper, a node-centered finite volume method based on triangulations for a fracture model is presented, in which we restrict the pressure to the linear finite element space while the velocity can be approximated by constant vectors element by element. The numerical scheme is established just associated with the pressure to avoid the saddle-point problem. Error estimates of O(h) accuracy for the discrete H1 semi-norm and the discrete L2 norm of pressure p and the (L2)2 norm of velocity u are developed on general triangulations. Under an additional assumption about essentially symmetric control volumes, the error estimates for the pressure p can be improved to O(h3/2). Finally, numerical experiments are carried out to verify the accuracy and convergence rates for the proposed finite volume scheme.

[1]  Zhiqiang Cai,et al.  On the finite volume element method , 1990 .

[2]  A. Poynter.  Thematic mapping using microcomputers: a commercial map producer's viewpoint , 1985 .

[3]  Jean E. Roberts,et al.  Space-time Domain Decomposition and Mixed Formulation for solving reduced fracture models , 2015 .

[4]  Vincent Martin,et al.  Modeling fractures as interfaces with nonmatching grids , 2012, Computational Geosciences.

[5]  Michael G. Edwards,et al.  Finite volume discretization with imposed flux continuity for the general tensor pressure equation , 1998 .

[6]  Zhiqiang Cai,et al.  The finite volume element method for diffusion equations on general triangulations , 1991 .

[7]  Fernando A. Morales,et al.  The narrow fracture approximation by channeled flow , 2010 .

[8]  Zhiqiang Cai,et al.  On the accuracy of the finite volume element method for diffusion equations on composite grids , 1990 .

[9]  R. Eymard,et al.  Finite Volume Methods , 2019, Computational Methods for Fluid Dynamics.

[10]  Zhengjia Sun,et al.  A block-centered finite difference method for reduced fracture model in Karst aquifer system , 2017, Comput. Math. Appl..

[11]  Jerome Droniou,et al.  FINITE VOLUME SCHEMES FOR DIFFUSION EQUATIONS: INTRODUCTION TO AND REVIEW OF MODERN METHODS , 2014, 1407.1567.

[12]  Mary F. Wheeler,et al.  Convergence of a symmetric MPFA method on quadrilateral grids , 2007 .

[13]  R. Herbin,et al.  An Error Estimate for a Nite Volume Scheme for a Diiusion Convection Problem on a Triangular Mesh , 1995 .

[14]  I. Faille,et al.  A control volume method to solve an elliptic equation on a two-dimensional irregular mesh , 1992 .

[15]  D. Rose,et al.  Some errors estimates for the box method , 1987 .

[16]  Vincent Martin,et al.  Modeling Fractures and Barriers as Interfaces for Flow in Porous Media , 2005, SIAM J. Sci. Comput..

[17]  R. Eymard,et al.  Finite volume approximation of elliptic problems and convergence of an approximate gradient , 2001 .

[18]  Philippe Angot,et al.  ASYMPTOTIC AND NUMERICAL MODELLING OF FLOWS IN FRACTURED POROUS MEDIA , 2009 .

[19]  Alfio Quarteroni,et al.  A multiscale Darcy–Brinkman model for fluid flow in fractured porous media , 2011, Numerische Mathematik.

[20]  Ivar Aavatsmark,et al.  Discretization on Non-Orthogonal, Quadrilateral Grids for Inhomogeneous, Anisotropic Media , 1996 .

[21]  Alessio Fumagalli,et al.  Numerical modelling of multiphase subsurface flow in the presence of fractures , 2011 .