Locally conservative Galerkin and finite volume methods for two-phase flow in porous media

In this paper we propose a locally conservative Galerkin (LCG) finite element method for two-phase flow simulations in heterogeneous porous media. The main idea of the method is to use the property of local conservation at steady state conditions in order to define a numerical flux at element boundaries. It provides a way to apply standard Galerkin finite element method in two-phase flow simulations in porous media. The LCG method has all the advantages of the standard finite element method while explicitly conserving fluxes over each element. All the examples presented show that the formulation employed is accurate and robust, while using less CPU time than finite volume method.

[1]  Michael Andrew Christie,et al.  Tenth SPE Comparative Solution Project: a comparison of upscaling techniques , 2001 .

[2]  Haiying Wang,et al.  Locally Conservative Fluxes for the Continuous Galerkin Method , 2007, SIAM J. Numer. Anal..

[3]  Perumal Nithiarasu,et al.  The locally conservative Galerkin (LCG) method for solving the incompressible Navier–Stokes equations , 2008 .

[4]  Chenchen Wang,et al.  Permeability analysis of fractured vuggy porous media based on homogenization theory , 2010 .

[5]  Long Chen.  FINITE ELEMENT METHOD , 2013 .

[6]  Matthew W. Farthing,et al.  Locally conservative, stabilized finite element methods for variably saturated flow , 2008 .

[7]  J. Z. Zhu,et al.  The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique , 1992 .

[8]  Robert L. Lee,et al.  THE CONSISTENT GALERKIN FEM FOR COMPUTING DERIVED BOUNDARY QUANTITIES IN THERMAL AND/OR FLUIDS PROBLEMS , 1987 .

[9]  P. Nithiarasu.  A SIMPLE LOCALLY CONSERVATIVE GALERKIN (LCG) FINITE-ELEMENT METHOD FOR TRANSIENT CONSERVATION EQUATIONS , 2004 .

[10]  Thomas J. R. Hughes,et al.  The Continuous Galerkin Method Is Locally Conservative , 2000 .

[11]  A. Firoozabadi,et al.  Control‐volume method for numerical simulation of two‐phase immiscible flow in two‐ and three‐dimensional discrete‐fractured media , 2004 .

[12]  P. Nithiarasu,et al.  An element‐wise, locally conservative Galerkin (LCG) method for solving diffusion and convection–diffusion problems , 2008 .

[13]  Thomas J. R. Hughes,et al.  Conservation properties for the Galerkin and stabilised forms of the advection–diffusion and incompressible Navier–Stokes equations , 2005 .

[14]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[15]  John W. Barrett,et al.  Total Flux Estimates for a Finite-Element Approximation of Elliptic Equations , 1986 .

[16]  Chenchen Wang,et al.  Numerical Calculation of Equivalent Permeability Tensor for Fractured Vuggy Porous Media Based on Homogenization Theory , 2011 .