A generalized Mimetic Finite Difference method and Two-Point Flux schemes over Voronoi diagrams

Al-Hinai, Omar; Wheeler, Mary F.; Yotov, Ivan

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 51 (2017) no. 2, pp. 679-706.

We develop a generalization of the mimetic finite difference (MFD) method for second order elliptic problems that extends the family of convergent schemes to include two-point flux approximation (TPFA) methods over general Voronoi meshes, which are known to satisfy the discrete maximum principle. The method satisfies a modified consistency condition, which utilizes element and face weighting functions. This results in shifting the points on the elements and faces where the pressure and the flux are most accurately approximated. The flux bilinear form is non-symmetric in general, although it reduces to a symmetric form in the case of TPFA. It can be defined as the $L^2$-inner product of vectors in two $H(\Omega ;div)$ discrete spaces, which are constructed via suitable lifting operators. A specific construction of such lifting operators is presented on rectangles. We note that a different choice is made for test and trial spaces, therefore the method can be viewed as a $H(\Omega ;div)$-conforming Petrov–Galerkin Mixed Finite Element method. We prove first-order convergence in pressure and flux, and superconvergence of the pressure under further restrictions. We present numerical results that support the theory.

Classification : 65M60, 65N08, 76S05
Mots clés : Mimetic finite difference, finite volume methods, discrete maximum principle, polyhedral meshes, Voronoi diagrams

I. Aavatsmark, An introduction to multipoint flux approximations for quadrilateral grids. Comput. Geosci. 6 (2002) 405–432. | DOI | MR | Zbl

J. Aghili, S. Boyaval and D. Di Pietro, Hybridization of mixed high-order methods on general meshes and application to the stokes equations. Comput. Methods Appl. Math. 15 (2015) 111–134. | DOI | MR | Zbl

S. Agmon, Lectures on elliptic boundary value problems. Vol. 2. American Mathematical Soc. (2010). | MR

T. Arbogast, M. Wheeler and I. Yotov, Mixed finite elements for elliptic problems with tensor coefficients as cell-centered finite differences. SIAM J. Numer. Anal. 34 (1997) 828–852. | DOI | MR | Zbl

T. Arobogast, C. Dawson, P. Keenan, M. Wheeler and I. Yotov, Enhanced cell-centered finite differences for elliptic equations on general geometry. SIAM J. Scientific Comput. 19 (1998) 404–425. | DOI | MR | Zbl

K. Aziz, Reservoir simulation grids: Opportunities and problems. J. Petroleum Technol. (1993) 658–663

S. Balay, M. Adams, J. Brown, P. Brune, K. Buschelman, V. Eijkhout, W.D. Gropp, D. Kaushik, M. Knepley, L. McInnes, K. Rupp, B. Smith and H. Zhang, PETSc Web page, http://www.mcs.anl.gov/petsc (2014).

J. Baranger, J.-F. Maitre and F. Oudin, Connection between finite volume and mixed finite element methods. RAIRO-Modél. Math. Anal. Numér. 30 (1996) 445–465. | DOI | Numdam | MR | Zbl

S. Behnel, R. Bradshaw, C. Citro, L. Dalcin, D. Seljebotn and K. Smith, Cython: The best of both worlds. Comput. Sci. Eng. 13 (2011) 31–39. | DOI

P. Berry, S Bonduá, V. Bortolotti, C. Cormio and E.M. Vasini, A gis-based open source pre-processor for georesources numerical modeling. Environ. Model. Softw. 62 (2014) 52–64. | DOI

S. Brenner and L. Scott, The Mathematical Theory of Finite Element Methods. Springer, 2nd edition (2002). | MR | Zbl

F. Brezzi and M. Fortin, Mixed and hybrid finite element methods. Springer-Verlag, New York (1991). | MR | Zbl

F. Brezzi, J. Douglas Jr and L. Marini, Two families of mixed finite elements for second order elliptic problems. Numer. Math. 47 (1985) 217–235. | DOI | MR | Zbl

F. Brezzi, K. Lipnikov and M. Shashkov, Convergence of the mimetic finite difference method for diffusion problems on polyhedral meshes. SIAM J. Numer. Anal. 43 (2005) 1872–1896. | DOI | MR | Zbl

F. Brezzi, K. Lipnikov and V. Simoncini, A family of mimetic finite difference methods on polygonal and polyhedral meshes. Math. Models Methods Appl. Sci. 15 (2005) 1533–1551. | DOI | MR | Zbl

F. Brezzi, R. S, Falk and L. D, Marini. Basic principles of mixed virtual element methods. ESAIM: M2AN 48 (2014) 1227–1240. | DOI | Numdam | MR | Zbl

A. Cangiani and G. Manzini, Flux reconstruction and solution post-processing in mimetic finite difference methods. Comput. Methods Appl. Mech. Engrg. 197 (2008) 933–945. | DOI | MR | Zbl

B. Cockburn, J. Gopalakrishnan and R. Lazarov, Unified hybridization of discontinuous galerkin, mixed, and continuous galerkin methods for second order elliptic problems. SIAM J. Numer. Anal. 47 (2009) 1319–1365. | DOI | MR | Zbl

L. Da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. Marini and A. Russo, Basic principles of virtual element methods. Math. Models Methods Appl. Sci. 23 (2013) 199–214. | DOI | MR | Zbl

L. Dalcin, petsc4py. https://code.google.com/p/petsc4py/ (2013).

D.A. Di Pietro, A. Ern and S. Lemaire, An arbitrary-order and compact-stencil discretization of diffusion on general meshes based on local reconstruction operators. Comput. Methods Appl. Math. 14 (2014) 461–472. | DOI | MR | Zbl

J. Douglas Jr and J. Wang, Superconvergence of mixed finite element methods on rectangular domains. Calcolo 26 (1989) 121–133. | DOI | MR | Zbl

J. Droniou, R. Eymard, T. Gallouët and R. Herbin, A unified approach to mimetic finite difference, hybrid finite volume and mixed finite volume methods. Math. Models Methods Appl Sci. 20 (2010) 265–295. | DOI | MR | Zbl

J. Droniou, R. Eymard, T. Gallouët and R. Herbin, Gradient schemes: a generic framework for the discretisation of linear, nonlinear and nonlocal elliptic and parabolic equations. Math. Models Methods Appl. Sci. 23 (2013) 2395–2432. | DOI | MR | Zbl

M. Edwards, Unstructured, control-volume distributed, full-tensor finite-volume schemes with flow based grids. Comput. Geosci. 6 (2002) 433–452. | DOI | MR | Zbl

M. Edwards and C. Rogers, Finite volume discretization with imposed flux continuity for the general tensor pressure equation. Comput. Geosci. 2 (1998) 259–290. | DOI | MR | Zbl

R. Ewing, R. Lazarov and J. Wang, Superconvergence of the velocity along the gauss lines in mixed finite element methods. SIAM J. Numer. Anal. 28 (1991) 1015–1029. | DOI | MR | Zbl

R. Eymard, T. Gallouët and R. Herbin, Finite volume methods. Vol. VII of Handbook of Numerical Analysis. North-Holland, Amsterdam (2000). | MR | Zbl

C.M. Freeman, K.L. Boyle, M. Reagan, J. Johnson, C. Rycroft and G.J. Moridis, Meshvoro: A three-dimensional voronoi mesh building tool for the tough family of codes. Comput. Geosci. 70 (2014) 26–34. | DOI

L.S.K. Fung et al., Parallel unstructured-solver methods for simulation of complex giant reservoirs. SPE J. 13 (2008) 440–446. | DOI

V. Girault and P.-A. Raviart, Finite Element Methods for Navier-Stokes Equations. Number 5 in Springer series in Computational Mathematics. Springer-Verlag (1986). | MR | Zbl

P. Grisvard, Elliptic problems in nonsmooth domains. Vol. 69. SIAM (2011). | MR | Zbl

V. Gyrya, K. Lipnikov, G. Manzini and D. Svyatskiy, M-adaptation in the mimetic finite difference method. Math. Models Methods Appl. Sci. 24 (2014) 1621–1663. | DOI | MR | Zbl

A. Henderson, The paraview guide: A parallel visualization application. kitware inc. Technical report (2007).

R. Ingram, M. Wheeler and I. Yotov, A multipoint flux mixed finite element method on hexahedra. SIAM J. Numer. Anal. 48 (2010) 1281–1312. | DOI | MR | Zbl

E. Jones et al., SciPy: Open source scientific tools for Python (2001).

R. Klausen and A. Stephansen, Convergence of multi-point flux approximation on general grids and media. Int. J. Numer. Anal. Model. 9 (2012) 584–606. | MR | Zbl

R. Klausen and R. Winther, Convergence of multipoint flux approximations on quadrilateral grids. Numer. Methods Partial Differ. Eq. 22 (2006) 1438–1454. | DOI | MR | Zbl

R. Klausen and R. Winther, Robust convergence of multi point flux approximation on rough grids. Numer. Math. 104 (2006) 317–337. | DOI | MR | Zbl

Y. Kuznetsov and S. Repin, New mixed finite element method on polygonal and polyhedral meshes. Russian J. Numer. Anal. Math. Model. 18 (2003) 261–278. | DOI | MR | Zbl

K. Lipnikov, M. Shashkov and I. Yotov, Local flux mimetic finite difference methods. Numer. Math. 112 (2009) 115–152. | DOI | MR | Zbl

K. Lipnikov, G. Manzini and D. Svyatskiy, Analysis of the monotonicity conditions in the mimetic finite difference method for elliptic problems. J. Comput. Phys. 230 (2011) 2620–2642. | DOI | MR | Zbl

K. Lipnikov, G. Manzini and M. Shashkov, Mimetic finite difference method. Physics-compatible numerical methods. J. Comput. Phys. 257 (2014) 1163–1227. | DOI | MR | Zbl

B. Mirtich, Fast and accurate computation of polyhedral mass properties. J. Graphics Tools 1 (1996) 31–50. | DOI

L. Mu, J. Wang and X. Ye, Weak galerkin finite element methods on polytopal meshes. Int. J. Numer. Anal. Model. 12 (2015) 536–566. | MR

M. Nakata, A. Weiser and M. Wheeler, Some superconvergence results for mixed finite element methods for elliptic problems on rectangular domains. The Mathematics of Finite Elements and Applications V, edited by J.R. Whiteman. Academic Press, London (1985) 367–389. | MR | Zbl

C. Palagi and K. Aziz, Use of Voronoi grid in reservoir simluation. SPE Adv. Technol. Ser. 2 (1994) 69–77. | DOI

S. Panday, C.D. Langevin, R.G. Niswonger, M. Ibaraki and J.D. Hughes, Modflow-usg version 1: An unstructured grid version of MODFLOW for simulating groundwater flow and tightly coupled processes using a control volume finite-difference formulation. Technical report, US Geological Survey (2013).

P. Raviart and J. Thomas, A mixed finite element method for 2-nd order elliptic problems. In Mathematical Aspects of Finite Element Methods. Springer (1977) 292–315. | MR | Zbl

T. Russell and M. Wheeler, Finite element and finite difference methods for continuous flows in porous media. In Vol. 1 of Mathematics of Reservoir Simulation, Frontiers in Applied Mathematics, edited by R.E. Ewing. SIAM (1983). | MR | Zbl

C. Rycroft, Voro++: A three-dimensional voronoi cell library in c++ (2009).

Schlumberger Geoquest, Eclipse 100 reference manual. Schlumberger Geoquest (2001).

S. Van Der Walt, S. Colbert and G. Varoquaux, The numpy array: a structure for efficient numerical computation. Comput. Sci. Eng. 13 (2011) 22–30. | DOI

M. Vohralík and B. Wohlmuth, Mixed finite element methods: implementation with one unknown per element, local flux expressions, positivity, polygonal meshes, and relations to other methods. Math. Models Methods Appl. Sci. 23 (2013) 803–838. | DOI | MR | Zbl

A. Weiser and M. Wheeler, On convergence of block-centered finite differences for elliptic problems. SIAM J. Numer. Anal. 25 (1988) 351–375. | DOI | MR | Zbl

M. Wheeler, G. Xue and I. Yotov, A multipoint flux mixed finite element method on distorted quadrilaterals and hexahedra. Numer. Math. 121 (2012) 165–204. | DOI | MR | Zbl

M. Wheeler and I. Yotov, A multipoint flux mixed finite element method. SIAM J. Numer. Anal. 44 (2006) 2082–2106. | DOI | MR | Zbl