A novel three-dimensional element free Galerkin (EFG) code for simulating two-phase fluid flow in porous materials

Abstract In the past few decades, numerical simulation of multiphase flow systems has received increasing attention because of its importance in various fields of science and engineering. In this paper, a three-dimensional numerical model is developed for the analysis of simultaneous flow of two fluids through porous media. The numerical approach is fairly new based on the element-free Galerkin (EFG) method. The EFG is a type of mesh-less method which has rarely been used in the field of flow in porous media. The weak forms of the governing partial differential equations are derived by applying the weighted residual method and Galerkin technique. The penalty method is utilized for imposition of the essential boundary conditions. To create the discrete equation system, the EFG shape functions are used for spatial discretization of pore fluid pressures and a fully implicit scheme is employed for temporal discretization. The obtained numerical results indicate that the EFG method has the capability to substitute the classical FE and FD approaches from the accuracy point of view, provided that the efficiency of the EFG is improved. The developed EFG code can be used as a robust numerical tool for simulating two-phase flow processes in the subsurface layers in various engineering disciplines.

[1]  Karsten Pruess,et al.  Robust numerical methods for saturated-unsaturated flow with dry initial conditions in heterogeneous media , 1995 .

[2]  Mary F. Wheeler,et al.  Discontinuous galerkin finite element methods applied to two-phase, air-water flow problems , 2005 .

[3]  Akira Murakami,et al.  Mesh-free Method for Soil-water Coupled Problem within Finite Strain and Its Numerical Validity , 2005 .

[4]  Ming-Wan Lu,et al.  Least‐squares collocation meshless method , 2001 .

[5]  Jonathan F. Sykes,et al.  Compositional simulation of groundwater contamination by organic compounds: 1. Model development and verification , 1993 .

[6]  B. Ataie‐Ashtiani,et al.  Comparison of Numerical Formulations for Two-phase Flow in Porous Media , 2010 .

[7]  Guirong Liu.  Mesh Free Methods: Moving Beyond the Finite Element Method , 2002 .

[8]  Michel Vauclin,et al.  Experimental and numerical analysis of two-phase infiltration in a partially saturated soil , 1986 .

[9]  George F. Pinder,et al.  A Multiphase Approach to the Modeling of Porous Media Contamination by Organic Compounds: 1. Equation Development , 1985 .

[10]  Ping Lin,et al.  Numerical analysis of Biot's consolidation process by radial point interpolation method , 2002 .

[11]  B. Rivière,et al.  Adaptive simulations of two-phase flow by discontinuous Galerkin methods , 2006 .

[12]  P. Huyakorn,et al.  A three‐dimensional finite‐element model for simulating water flow in variably saturated porous media , 1986 .

[13]  Emil O. Frind,et al.  Two‐phase flow in heterogeneous porous media: 1. Model development , 1991 .

[14]  H. Shodja,et al.  Analysis of two-phase flow of compressible immiscible fluids through nondeformable porous media using moving finite elements , 1993 .

[15]  Charles R. Faust,et al.  Simulation of three‐dimensional flow of immiscible fluids within and below the unsaturated zone , 1989 .

[16]  Maren Paul,et al.  Simulation of two-phase flow processes in heterogeneous porous media with adaptive methods , 2003 .

[17]  B. Schrefler,et al.  The Finite Element Method in the Static and Dynamic Deformation and Consolidation of Porous Media , 1998 .

[18]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[19]  O. Oung,et al.  Numerical modelling of two-phase flow in a geocentrifuge , 2003, Environmental Modelling & Software.

[20]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[21]  Emil O. Frind,et al.  Two‐phase flow in heterogeneous porous media: 2. Model application , 1991 .

[22]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[23]  S. Atluri,et al.  A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics , 1998 .

[24]  E. Oñate,et al.  A FINITE POINT METHOD IN COMPUTATIONAL MECHANICS. APPLICATIONS TO CONVECTIVE TRANSPORT AND FLUID FLOW , 1996 .

[25]  Philip John Binning,et al.  A mass conservative numerical solution for two‐phase flow in porous media with application to unsaturated flow , 1992 .

[26]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[27]  Philip John Binning,et al.  Practical implementation of the fractional flow approach to multi-phase flow simulation , 1999 .

[28]  George F. Pinder,et al.  On the Simulation of Nonaqueous Phase Organic Compounds in the Subsurface , 1986 .

[29]  François Lehmann,et al.  Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media , 1998 .

[30]  J. Philip.  THE THEORY OF INFILTRATION: 1. THE INFILTRATION EQUATION AND ITS SOLUTION , 1957 .

[31]  K. Soga,et al.  Some numerical issues using element‐free Galerkin mesh‐less method for coupled hydro‐mechanical problems , 2009 .

[32]  Ali Pak,et al.  Three-dimensional simulation of fully coupled hydro-mechanical behavior of saturated porous media using Element Free Galerkin (EFG) method , 2012 .

[33]  Charles R. Faust,et al.  Transport of Immiscible Fluids Within and Below the Unsaturated Zone: A Numerical Model , 1985 .

[34]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[35]  K. Rathfelder,et al.  Mass balance errors in modeling two-phase immiscible flows: causes and remedies , 1993 .

[36]  J. F. Sykes,et al.  Numerical Modeling of Immiscible Organic Transport at the Hyde Park Landfill , 1986 .

[37]  B. Rivière,et al.  Fully implicit discontinuous finite element methods for two-phase flow , 2007 .

[38]  Christoph Grüninger.  Discontinuous Galerkin methods for two-phase flows in porous media , 2010 .

[39]  Hubert J. Morel-Seytoux,et al.  A Two‐Phase Numerical Model for Prediction of Infiltration: Applications to a Semi‐Infinite Soil Column , 1985 .

[40]  Jack C. Parker,et al.  A parametric model for constitutive properties governing multiphase flow in porous media , 1987 .

[41]  J. Parker,et al.  An efficient finite element method for modeling multiphase flow , 1989 .

[42]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[43]  Kang Tai,et al.  RADIAL POINT INTERPOLATION COLLOCATION METHOD (RPICM) FOR THE SOLUTION OF TWO PHASE FLOW THROUGH POROUS MEDIA , 2003 .

[44]  George F. Pinder,et al.  A Multiphase Approach to the Modeling of Porous Media Contamination by Organic Compounds: 2. Numerical Simulation , 1985 .

[45]  J. Philip,et al.  THE THEORY OF INFILTRATION: 2. THE PROFILE OF INFINITY , 1957 .

[46]  R. H. Brooks,et al.  Hydraulic properties of porous media , 1963 .