A robust linearization scheme for finite volume based discretizations for simulation of two-phase flow in porous media

In this work we consider a mathematical model for two-phase flow in porous media. The fluids are assumed immiscible and incompressible and the solid matrix non-deformable. The mathematical model for the two-phase flow is written in terms of the global pressure and a complementary pressure (obtained by using the Kirchhoff transformation) as primary unknowns. For the spatial discretization, finite volumes have been used (more precisely the multi-point flux approximation method) and in time the backward Euler method has been employed. We present here a new linearization scheme for the nonlinear system arising after the temporal and spatial discretization. We show that the scheme is linearly convergent. Numerical experiments are presented that sustain the theoretical results.

[1]  Luca Bergamaschi,et al.  MIXED FINITE ELEMENTS AND NEWTON-TYPE LINEARIZATIONS FOR THE SOLUTION OF RICHARDS' EQUATION , 1999 .

[2]  Jeremy E. Kozdon,et al.  Multidimensional upstream weighting for multiphase transport on general grids , 2012, Computational Geosciences.

[3]  Florin A. Radu,et al.  Convergence of MPFA on triangulations and for Richards' equation , 2008 .

[4]  Michael G. Edwards,et al.  Unstructured, Control-Volume Distributed, Full-Tensor Finite-Volume Schemes with Flow Based Grids , 2002 .

[5]  Peter Knabner,et al.  Fully coupled generalized hybrid-mixed finite element approximation of two-phase two-component flow in porous media. Part I: formulation and properties of the mathematical model , 2013, Computational Geosciences.

[6]  Ragnar Winther,et al.  Robust convergence of multi point flux approximation on rough grids , 2006, Numerische Mathematik.

[7]  P. Knabner,et al.  Numerical Methods for Elliptic and Parabolic Partial Differential Equations , 2003, Texts in Applied Mathematics.

[8]  F. Radu,et al.  Mixed finite elements for the Richards' equation: linearization procedure , 2004 .

[9]  Zhangxin Chen,et al.  Degenerate two-phase incompressible flow III. Sharp error estimates , 2001, Numerische Mathematik.

[10]  M. Wheeler,et al.  Mixed Finite Elements for Elliptic Problems with Tensor Coefficients as Cell-Centered Finite Differences , 1997 .

[11]  J. Nordbotten.  Cell‐centered finite volume discretizations for deformable porous media , 2014 .

[12]  Florin A. Radu,et al.  Newton method for reactive solute transport with equilibrium sorption in porous media , 2010, J. Comput. Appl. Math..

[13]  Zhangxin Chen.  Degenerate Two-Phase Incompressible Flow: I. Existence, Uniqueness and Regularity of a Weak Solution , 2001 .

[14]  Eun-Jae Park,et al.  Mixed finite element methods for nonlinear second-order elliptic problems , 1995 .

[15]  Peter Knabner,et al.  Fully coupled generalised hybrid-mixed finite element approximation of two-phase two-component flow in porous media. Part II: numerical scheme and numerical results , 2012, Computational Geosciences.

[16]  Brett J. Larkin,et al.  An ANSI C program to determine in expected linear time the vertices of the convex hull of a set of planar points , 1991 .

[17]  S. Kräutle,et al.  The semismooth Newton method for multicomponent reactive transport with minerals , 2011 .

[18]  Ivar Aavatsmark,et al.  Monotonicity for MPFA methods on triangular grids , 2011 .

[19]  Florin A. Radu,et al.  Mixed finite element discretization and Newton iteration for a reactive contaminant transport model with nonequilibrium sorption: convergence analysis and error estimates , 2011 .

[20]  G. Chavent.  Mathematical models and finite elements for reservoir simulation , 1986 .

[21]  Jan M. Nordbotten.  Convergence of a Cell-Centered Finite Volume Discretization for Linear Elasticity , 2015, SIAM J. Numer. Anal..

[22]  Roland Masson,et al.  Convergence of Finite Volume MPFA O type Schemes for Heterogeneous Anisotropic Diffusion Problems on General Meshes , 2010 .

[23]  Ivar Aavatsmark,et al.  Sufficient criteria are necessary for monotone control volume methods , 2009, Appl. Math. Lett..

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

[25]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[26]  Todd Arbogast.  The existence of weak solutions to single porosity and simple dual-porosity models of two-phase incompressible flow , 1992 .

[27]  I. Aavatsmark,et al.  An Introduction to Multipoint Flux Approximations for Quadrilateral Grids , 2002 .

[28]  Ivar Aavatsmark,et al.  Monotonicity of control volume methods , 2007, Numerische Mathematik.

[29]  F. Brezzi,et al.  Mathematical models and finite elements for reservoir simulation : single phase, multiphase and multicomponent flows through porous media , 1988 .

[30]  Mauricio Sepúlveda,et al.  Error estimates for the finite volume discretization for the porous medium equation , 2010, J. Comput. Appl. Math..

[31]  Marián Slodicka,et al.  A Robust and Efficient Linearization Scheme for Doubly Nonlinear and Degenerate Parabolic Problems Arising in Flow in Porous Media , 2001, SIAM J. Sci. Comput..

[32]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[33]  Jan M Nordbotten,et al.  Practical Modeling Approaches for Geological Storage of Carbon Dioxide , 2009, Ground water.

[34]  Peter Knabner,et al.  Newton—Type Methods for the Mixed Finite Element Discretization of Some Degenerate Parabolic Equations , 2006 .