Modeling subsurface water resource systems involving heterogeneous porous media using the variational multiscale formulation

Summary This work compares two popular mixed finite element formulations used to model subsurface flow and transport in heterogeneous porous media, namely, the lowest order Raviart–Thomas and the variational multiscale stabilized formulations. Comparison is made based on performance for several problems of engineering relevance that involve highly heterogeneous material properties (permeability ratios of up to 1 × 10 5 ), open flow boundary conditions (pressure driven flows), and large scale domains in two dimensions. Numerical experiments are performed to show the degree to which mass conservation is violated when a flow field computed using either element is used as the advection velocity in a transport model. The results reveal that the equal-order interpolation under the variational multiscale formulation shows considerable mass production or loss for problems that involve flow tangential to layers of differing permeability. But the violation of local mass balance is marginal for problems in which flow is orthogonal to the layers of differing permeability. For problems involving pressure dependent viscosity, we show that models with a high degree of pressure dependence exhibit improved performance for the variational multiscale method. The results are useful in establishing rudimentary estimates of the error produced by using the variational mutliscale formulation for several different types of problems related to subsurface water resource systems.

[1]  W. Wakeham,et al.  The Transport Properties of Carbon Dioxide , 1990 .

[2]  E. N. Da C. Andrade.  Viscosity of liquids , 1952, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[3]  I. Babuska.  Error-bounds for finite element method , 1971 .

[4]  W. D. Evans,et al.  PARTIAL DIFFERENTIAL EQUATIONS , 1941 .

[5]  Thomas J. R. Hughes,et al.  Weak imposition of Dirichlet boundary conditions in fluid mechanics , 2007 .

[6]  Mary F. Wheeler,et al.  Compatible algorithms for coupled flow and transport , 2004 .

[7]  Rainer Helmig,et al.  CO2 leakage through an abandoned well: problem-oriented benchmarks , 2007 .

[8]  Carl Barus,et al.  Isothermals, isopiestics and isometrics relative to viscosity , 1893, American Journal of Science.

[9]  K. R. Rajagopal,et al.  A numerical study of fluids with pressure‐dependent viscosity flowing through a rigid porous medium , 2009, ArXiv.

[10]  D. Z. Turner,et al.  A stabilized formulation for the advection–diffusion equation using the Generalized Finite Element Method , 2011 .

[11]  Michael A. Heroux,et al.  Amesos: A Set of General Interfaces to Sparse Direct Solver Libraries , 2006, PARA.

[12]  T. Hughes,et al.  Two classes of mixed finite element methods , 1988 .

[13]  Todd Arbogast,et al.  A Locally Conservative Eulerian-Lagrangian Method for a Model Two-Phase Flow Problem in a One-Dimensional Porous Medium , 2012, SIAM J. Sci. Comput..

[14]  R. Sani,et al.  Résumé and remarks on the open boundary condition minisymposium , 1994 .

[15]  Albert J. Valocchi,et al.  Non-negative mixed finite element formulations for a tensorial diffusion equation , 2008, J. Comput. Phys..

[16]  Mario J. Martinez,et al.  Comparison of Galerkin and control volume finite element for advection–diffusion problems , 2006 .

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

[18]  B. Strazisar,et al.  Modeling carbon dioxide sequestration in saline aquifers: Significance of elevated pressures and salinities , 2005 .

[19]  Olivier Pironneau,et al.  Pressure boundary condition for the time‐dependent incompressible Navier–Stokes equations , 2006 .

[20]  Evan H Abramson,et al.  Viscosity of carbon dioxide measured to a pressure of 8 GPa and temperature of 673 K. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  L. D. Marini,et al.  Two families of mixed finite elements for second order elliptic problems , 1985 .

[22]  Harsha Nagarajan,et al.  Enforcing the non‐negativity constraint and maximum principles for diffusion with decay on general computational grids , 2010, ArXiv.

[24]  L. Durlofsky.  Accuracy of mixed and control volume finite element approximations to Darcy velocity and related quantities , 1994 .

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

[26]  Arif Masud,et al.  A multiscale/stabilized finite element method for the advection–diffusion equation , 2004 .

[27]  Frédéric Valentin,et al.  A Petrov-Galerkin enriched method : A mass conservative finite element method for the Darcy equation , 2007 .

[28]  Stefan Bachu,et al.  Semianalytical solution for CO2 leakage through an abandoned well. , 2005, Environmental science & technology.

[29]  O. Ladyzhenskaya.  The Boundary Value Problems of Mathematical Physics , 1985 .

[30]  Mark A. Herkommer,et al.  Data-volume reduction of data gathered along lines using the correlation coefficient to determine breakpoints , 1985 .

[31]  P. W. Bridgman,et al.  The Viscosity of Liquids under Pressure. , 1925, Proceedings of the National Academy of Sciences of the United States of America.

[32]  Thomas J. R. Hughes,et al.  A stabilized mixed finite element method for Darcy flow , 2002 .

[33]  P. Raviart,et al.  A mixed finite element method for 2-nd order elliptic problems , 1977 .

[34]  Annamaria Mazzia,et al.  Comparison of 3D Flow Fields Arising in Mixed and Standard Unstructured Finite Elements , 2003, International Conference on Computational Science.

[35]  Thomas J. R. Hughes,et al.  A stabilized mixed discontinuous Galerkin method for Darcy flow , 2006 .

[36]  Philip M. Gresho,et al.  The implementation of normal and/or tangential boundary conditions in finite element codes for incompressible fluid flow , 1982 .

[37]  D. Z. Turner,et al.  A stabilized mixed finite element method for Darcy flow based on a multiscale decomposition of the solution , 2006 .

[38]  J. L. Gupta.  Fluid Motion Past a Porous Circular Cylinder With Initial Pressure Gradient , 1980 .