A MCC finite element approximation of incompressible miscible displacement in porous media

In this paper, we analyze an efficient numerical scheme for the coupled system of incompressible miscible displacement in porous media. Mass-conservative characteristic finite element is used for concentration equation, which continuous velocities are needed. Therefore we solve flow equation by mixed finite element method with continuous flux approximation. Then we show this scheme preserves the mass balance identity and optimal l ∞ ( L 2 ) error estimate of the concentration. Numerical experiments are given to confirm the convergence and the property of mass balance.

[1]  D. W. Peaceman.  Fundamentals of numerical reservoir simulation , 1977 .

[2]  H. Rui,et al.  An approximation of incompressible miscible displacement in porous media by mixed finite elements and symmetric finite volume element method of characteristics , 2013 .

[3]  Richard E. Ewing,et al.  The Mathematics of Reservoir Simulation , 2016 .

[4]  Dong Liang,et al.  An ELLAM-MFEM Solution Technique for Compressible Fluid Flows in Porous Media with Point Sources and Sinks , 2000 .

[5]  M. Wheeler,et al.  A characteristics-mixed finite element method for advection-dominated transport problems , 1995 .

[6]  Richard E. Ewing,et al.  A time-discretization procedure for a mixed finite element approximation of miscible displacement in porous media , 1983 .

[7]  Jean E. Roberts,et al.  Numerical methods for a model for compressible miscible displacement in porous media , 1983 .

[8]  Dong Liang,et al.  An Approximation to Miscible Fluid Flows in Porous Media with Point Sources and Sinks by an Eulerian-Lagrangian Localized Adjoint Method and Mixed Finite Element Methods , 2000, SIAM J. Sci. Comput..

[9]  Hongxing Rui,et al.  A second order characteristic finite element scheme for convection-diffusion problems , 2002, Numerische Mathematik.

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

[11]  Todd Arbogast,et al.  A Family of Rectangular Mixed Elements with a Continuous Flux for Second Order Elliptic Problems , 2004, SIAM Journal on Numerical Analysis.

[12]  T. F. Russell,et al.  Convergence analysis of an approximation of miscible displacement in porous media by mixed finite elements and a modified method of characteristics , 1984 .

[13]  Yirang Yuan,et al.  An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method , 2009 .

[14]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[15]  G. R. Shubin,et al.  An unsplit, higher order Godunov method for scalar conservation laws in multiple dimensions , 1988 .

[16]  J. Bear.  Hydraulics of Groundwater , 1979 .

[17]  T. F. Russell,et al.  Finite element and finite difference methods for continuous flows in porous media. , 1800 .

[18]  Hongxing Rui,et al.  A Mass-Conservative Characteristic Finite Element Scheme for Convection-Diffusion Problems , 2008, J. Sci. Comput..

[19]  T. F. Russell,et al.  Time Stepping Along Characteristics with Incomplete Iteration for a Galerkin Approximation of Miscible Displacement in Porous Media , 1985 .

[20]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

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