Error estimates for a mixed finite element discretization of a two-phase porous media flow model with dynamic capillarity

Abstract We analyze a fully discrete numerical scheme for the model describing two-phase immiscible flow in porous media with dynamic effects in the capillary pressure. We employ the Euler implicit method for the time discretization. The spatial discretization is based on the mixed finite element method (MFEM). Specifically, the lowest order Raviart–Thomas elements are applied. In this paper, the error estimates for the saturation, fluxes and phase pressures in L ∞ ( 0 , T ; L 2 ( Ω ) ) are derived for the temporal and spatial discretization to show the convergence of the scheme. Finally, we present some numerical results to support the theoretical findings.

[1]  Peter Knabner,et al.  Error estimates for a mixed finite element discretization of some degenerate parabolic equations , 2008, Numerische Mathematik.

[2]  William G. Gray,et al.  Thermodynamic basis of capillary pressure in porous media , 1993 .

[3]  Iuliu Sorin Pop,et al.  Uniqueness of weak solutions for a pseudo-parabolic equation modeling two phase flow in porous media , 2015, Appl. Math. Lett..

[4]  Iuliu Sorin Pop,et al.  A modified L-scheme to solve nonlinear diffusion problems , 2019, Comput. Math. Appl..

[5]  Timo J. Heimovaara,et al.  Nonequilibrium capillarity effects in two‐phase flow through porous media at different scales , 2011 .

[6]  I. Pop,et al.  Convergence of an MPFA finite volume scheme for a two‐phase porous media flow model with dynamic capillarity , 2018 .

[7]  Stephan Luckhaus,et al.  Flow of Oil and Water in a Porous Medium , 1984 .

[8]  I. Pop,et al.  Degenerate two-phase porous media flow model with dynamic capillarity , 2016 .

[9]  Louis J. Durlofsky,et al.  Coarse scale models of two phase flow in heterogeneous reservoirs: volume averaged equations and their relationship to existing upscaling techniques , 1998 .

[10]  Clément Cancès,et al.  An a posteriori error estimate for vertex-centered finite volume discretizations of immiscible incompressible two-phase flow , 2013, Math. Comput..

[11]  R. Hilfer,et al.  Dimensional analysis of two-phase flow including a rate-dependent capillary pressure-saturation relationship , 2008 .

[12]  Richard E. Ewing,et al.  The approximation of the pressure by a mixed method in the simulation of miscible displacement , 1983 .

[13]  Zhangxin Chen,et al.  Fully Discrete Finite Element Analysis of Multiphase Flow in Groundwater Hydrology , 1997 .

[14]  I. S. Pop,et al.  Two-Phase Flow in Porous Media: Dynamic Capillarity and Heterogeneous Media , 2016, Transport in Porous Media.

[15]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[16]  Stefan Karpinski,et al.  Analysis of a linearization scheme for an interior penalty discontinuous Galerkin method for two‐phase flow in porous media with dynamic capillarity effects , 2017 .

[17]  Jan M. Nordbotten,et al.  A robust, mass conservative scheme for two-phase flow in porous media including Hölder continuous nonlinearities , 2018 .

[18]  C. J. Duijn,et al.  Travelling wave solutions for the Richards equation incorporating non-equilibrium effects in the capillarity pressure , 2018, Nonlinear Analysis: Real World Applications.

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

[20]  Andro Mikelić,et al.  A global existence result for the equations describing unsaturated flow in porous med , 2010 .

[21]  BEN SCHWEIZER,et al.  Two-phase flow equations with a dynamic capillary pressure , 2012, European Journal of Applied Mathematics.

[22]  Stefan Karpinski,et al.  Analysis of an interior penalty discontinuous Galerkin scheme for two phase flow in porous media with dynamic capillary effects , 2017, Numerische Mathematik.

[23]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[24]  Yun Fan,et al.  Equivalent formulations and numerical schemes for a class of pseudo-parabolic equations , 2013, J. Comput. Appl. Math..

[25]  Peter Knabner,et al.  Analysis of an Upwind-Mixed Hybrid Finite Element Method for Transport Problems , 2014, SIAM J. Numer. Anal..

[26]  Alexander A. Weiss,et al.  Dynamic capillary effects in heterogeneous porous media , 2007 .

[27]  Significance of Dynamic Effect in Capillarity during Drainage Experiments in Layered Porous Media , 2010 .

[28]  Iuliu Sorin Pop,et al.  Two-phase porous media flows with dynamic capillary effects and hysteresis: Uniqueness of weak solutions , 2015, Comput. Math. Appl..

[29]  Iuliu Sorin Pop,et al.  Travelling wave solutions for degenerate pseudo-parabolic equation modelling two-phase flow in porous media , 2013 .

[30]  Jan M. Nordbotten,et al.  A robust linearization scheme for finite volume based discretizations for simulation of two-phase flow in porous media , 2015, J. Comput. Appl. Math..

[31]  Anthony Michel,et al.  A Finite Volume Scheme for Two-Phase Immiscible Flow in Porous Media , 2003, SIAM J. Numer. Anal..

[32]  F. Radu,et al.  A study on iterative methods for solving Richards’ equation , 2015, Computational Geosciences.

[33]  R. Helmig.  Multiphase Flow and Transport Processes in the Subsurface: A Contribution to the Modeling of Hydrosystems , 2011 .

[34]  Rainer Helmig,et al.  Variational inequalities for modeling flow in heterogeneous porous media with entry pressure , 2009 .

[35]  Todd Arbogast,et al.  Numerical methods for the simulation of flow in root-soil systems , 1993 .

[36]  Raphaèle Herbin,et al.  Mathematical study of a petroleum-engineering scheme , 2003 .

[37]  Richard E. Ewing,et al.  Galerkin Methods for Miscible Displacement Problems with Point Sources and Sinks - Unit Mobility Ratio Case, , 1984 .

[38]  Peter Knabner,et al.  Order of Convergence Estimates for an Euler Implicit, Mixed Finite Element Discretization of Richards' Equation , 2004, SIAM J. Numer. Anal..

[39]  Ben Schweizer,et al.  A well-posed hysteresis model for flows in porous media and applications to fingering effects , 2010 .

[40]  Florin A. Radu,et al.  Convergence analysis for a mixed finite element scheme for flow in strictly unsaturated porous media , 2014 .

[41]  David A. DiCarlo,et al.  Experimental measurements of saturation overshoot on infiltration , 2004 .

[42]  Mazen Saad,et al.  Analysis of a finite volume-finite element method for Darcy-Brinkman two-phase flows in porous media , 2018, J. Comput. Appl. Math..

[43]  Jan M. Nordbotten,et al.  Geological Storage of CO2: Modeling Approaches for Large-Scale Simulation , 2011 .

[44]  Iuliu Sorin Pop,et al.  A class of pseudo‐parabolic equations: existence, uniqueness of weak solutions, and error estimates for the Euler‐implicit discretization , 2011 .

[45]  Malgorzata Peszynska,et al.  Numerical methods for unsaturated flow with dynamic capillary pressure in heterogeneous porous media , 2008 .