An Efficient ELLAM Implementation for Modeling Solute Transport in Fractured Porous Media

The goal of this study is to introduce an adaptation of the Eulerian-Lagrangian localized adjoint method (ELLAM) for the simulation of mass transport in fractured porous media, and to evaluate the performance of ELLAM in such a case. The fractures are represented explicitly using the discrete fracture model. The velocity field was calculated using the mixed hybrid finite element method. A sound ELLAM implementation is developed to address numerical artifacts (spurious oscillations and numerical dispersion) arising from the presence of fractures. The efficiency of the developed ELLAM implementation was further improved by taking advantage of the parallel computing on shared memory architecture for the tasks related to particles tracking and linear system resolving. The performance of ELLAM was tested by comparing its results with the Eulerian discontinuous Galerkin method based on several benchmark problems dealing with different fracture configurations. The results highlight the robustness and accuracy of ELLAM, as it allows the use of large time steps, and overcomes the Courant-Friedrichs-Lewy (CFL) restriction. The outcome also shows that ELLAM is more efficient when fracture density is increased.

[1]  M. Lacroix,et al.  Interpolation techniques applied to the Eulerian-Lagrangian solution of the convection-dispersion equation in natural coordinates , 1997 .

[2]  Frederick Delay,et al.  Simulating Solute Transport in Porous or Fractured Formations Using Random Walk Particle Tracking: A Review , 2005 .

[3]  R. Naidu.  Recent Advances in Contaminated Site Remediation , 2013, Water, Air, & Soil Pollution.

[4]  Martin J. Blunt,et al.  Streamline-based simulation of advective–dispersive solute transport , 2004 .

[5]  A. Younes,et al.  On the use of large time steps with ELLAM for transport with kinetic reactions over heterogeneous domains , 2009 .

[6]  Georg Kosakowski,et al.  A class of analytical solutions for multidimensional multispecies diffusive transport coupled with precipitation‐dissolution reactions and porosity changes , 2012 .

[7]  A. Firoozabadi,et al.  An efficient numerical model for multicomponent compressible flow in fractured porous media , 2014 .

[8]  P. Huggenberger,et al.  Variable-density flow in heterogeneous porous media--laboratory experiments and numerical simulations. , 2009, Journal of contaminant hydrology.

[9]  Laurent Stainier,et al.  On the use of large time steps with an energy momentum conserving algorithm for non-linear hypoelastic constitutive models , 2004 .

[10]  B. Berkowitz.  Characterizing flow and transport in fractured geological media: A review , 2002 .

[11]  William G. Gray,et al.  Eulerian-Lagrangian localized adjoint methods with variable coefficients in multiple dimensions. , 1990 .

[12]  Joachim Moortgat,et al.  Implicit finite volume and discontinuous Galerkin methods for multicomponent flow in unstructured 3D fractured porous media , 2016, 1608.06894.

[13]  Bruce D. Jones,et al.  Smoothed particle hydrodynamics and its applications for multiphase flow and reactive transport in porous media , 2015, Computational Geosciences.

[14]  Marwan Fahs,et al.  A new approach to avoid excessive numerical diffusion in Eulerian–Lagrangian methods , 2007 .

[15]  Philippe Ackerer,et al.  Solving the advection-diffusion equation with the Eulerian-Lagrangian localized adjoint method on unstructured meshes and non uniform time stepping , 2005 .

[16]  A. Younes,et al.  An efficient numerical model for hydrodynamic parameterization in 2D fractured dual-porosity media , 2014 .

[17]  H. Langtangen,et al.  Mixed Finite Elements , 2003 .

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

[19]  B. D. Rogers,et al.  SMOOTHED PARTICLE HYDRODYNAMICS (SPH) IN CFD , 2014 .

[20]  Hong Wang,et al.  An Optimal-Order Error Estimate for an ELLAM Scheme for Two-Dimensional Linear Advection-Diffusion Equations , 2000, SIAM J. Numer. Anal..

[21]  Cass T. Miller,et al.  An ELLAM approximation for advective–dispersive transport with nonlinear sorption , 2006 .

[22]  T. Babadagli,et al.  Modeling miscible injection in fractured porous media using random walk simulation , 2012 .

[23]  A. Werner,et al.  Seawater intrusion in fractured coastal aquifers: A preliminary numerical investigation using a fractured Henry problem , 2015 .

[24]  Hussein Hoteit,et al.  Compositional Modeling of Discrete-Fractured Media Without Transfer Functions by the Discontinuous Galerkin and Mixed Methods , 2006 .

[25]  T. F. Russell,et al.  A Finite‐Volume ELLAM for Three‐Dimensional Solute‐Transport Modeling , 2003, Ground water.

[26]  Hong Wang,et al.  A multiscale Eulerian–Lagrangian localized adjoint method for transient advection–diffusion equations with oscillatory coefficients , 2009 .

[27]  P. Ackerer,et al.  A new efficient Eulerian–Lagrangian localized adjoint method for solving the advection–dispersion equation on unstructured meshes , 2006 .

[28]  Clint Dawson,et al.  Discontinuous Galerkin Methods for Modeling Hurricane Storm Surge , 2011 .

[29]  P. Ackerer,et al.  Solving the advection–dispersion equation with discontinuous Galerkin and multipoint flux approximation methods on unstructured meshes , 2008 .

[30]  J. S. Wang,et al.  Flow and Transport in Fractured Rocks , 1991 .

[31]  Roussos Dimitrakopoulos,et al.  Discretizing two‐dimensional complex fractured fields for incompressible two‐phase flow , 2011 .

[32]  Pallav Sarma,et al.  New transfer functions for simulation of naturally fractured reservoirs with dual porosity models , 2006 .

[33]  Abbas Firoozabadi,et al.  An implicit numerical model for multicomponent compressible two-phase flow in porous media , 2015 .

[34]  Hussein Hoteit,et al.  An efficient numerical model for incompressible two-phase flow in fractured media , 2008 .

[35]  Peter Dietrich,et al.  Flow and transport in fractured porous media , 2005 .

[36]  Richard W. Healy,et al.  Solution of the advection-dispersion equation in two dimensions by a finite-volume Eulerian-Lagrangian localized adjoint method , 1998 .

[37]  L. Durlofsky,et al.  An Efficient Discrete-Fracture Model Applicable for General-Purpose Reservoir Simulators , 2004 .

[38]  Mayur Pal,et al.  Control-Volume Distributed Multi-Point Flux Approximation Coupled with a Lower-Dimensional Fracture Model , 2015, ANSS 2015.

[39]  Philip John Binning,et al.  A forward particle tracking Eulerian-Lagrangian localized adjoint method for solution of the contaminant transport equation in three dimensions , 2002 .

[40]  G. I. Barenblatt,et al.  Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata] , 1960 .

[41]  Thierry A. Mara,et al.  A new benchmark semi-analytical solution for density-driven flow in porous media , 2014 .

[42]  Fikret Kačaroğlu,et al.  Review of Groundwater Pollution and Protection in Karst Areas , 1999 .

[43]  Richard E. Ewing,et al.  An optimal-order estimate for Eulerian-Lagrangian localized adjoint methods for variable-coefficient advection-reaction problems , 1996 .

[44]  P. Ackerer,et al.  Mixed finite elements for solving 2‐D diffusion‐type equations , 2010 .

[45]  A. Younes,et al.  Modeling 2D Multispecies Reactive Transport in Saturated/Unsaturated Porous Media with the Eulerian–Lagrangian Localized Adjoint Method , 2012, Water, Air, & Soil Pollution.

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

[47]  Jacques Bodin,et al.  Simulation of solute transport in discrete fracture networks using the time domain random walk method , 2003 .

[48]  Richard W. Healy,et al.  A finite‐volume Eulerian‐Lagrangian Localized Adjoint Method for solution of the advection‐dispersion equation , 1993 .

[49]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[50]  Knut-Andreas Lie,et al.  Discontinuous Galerkin methods for advective transport in single-continuum models of fractured media , 2009 .

[51]  Philip John Binning,et al.  A split operator approach to reactive transport with the forward particle tracking Eulerian Lagrangian localized adjoint method , 2004 .

[52]  H. Suk.  Generalized semi-analytical solutions to multispecies transport equation coupled with sequential first-order reaction network with spatially or temporally variable transport and decay coefficients , 2016 .

[53]  Satish Karra,et al.  Particle tracking approach for transport in three-dimensional discrete fracture networks , 2015, Computational Geosciences.

[54]  Albert J. Valocchi,et al.  A multidimensional streamline-based method to simulate reactive solute transport in heterogeneous porous media , 2010 .

[55]  Frederick Delay,et al.  Time domain random walk method to simulate transport by advection‐dispersion and matrix diffusion in fracture networks , 2001 .

[56]  H. Redkey,et al.  A new approach. , 1967, Rehabilitation record.

[57]  Timothy D. Scheibe,et al.  A smoothed particle hydrodynamics model for reactive transport and mineral precipitation in porous and fractured porous media , 2007 .

[58]  Richard E. Ewing,et al.  A Family of Eulerian-Lagrangian Localized Adjoint Methods for Multi-dimensional Advection-Reaction Equations , 1999 .

[59]  Philip John Binning,et al.  A Finite Volume Eulerian‐Lagrangian Localized Adjoint Method for Solution of the Contaminant Transport Equations in Two‐Dimensional Multiphase flow Systems , 1996 .

[60]  A. Werner,et al.  On the effects of preferential or barrier flow features on solute plumes in permeable porous media , 2016 .

[61]  Hussein Hoteit,et al.  Compositional Modeling by the Combined Discontinuous Galerkin and Mixed Methods , 2004 .

[62]  Shuyu Sun,et al.  Adaptive mixed finite element methods for Darcy flow in fractured porous media , 2016 .

[63]  W. Kinzelbach,et al.  The Random Walk Method in Pollutant Transport Simulation , 1988 .

[64]  H. Kazemi,et al.  NUMERICAL SIMULATION OF WATER-OIL FLOW IN NATURALLY FRACTURED RESERVOIRS , 1976 .

[65]  Michael A. Celia,et al.  A probabilistic collocation Eulerian-Lagrangian localized adjoint method on sparse grids for assessing CO2 leakage through wells in randomly heterogeneous porous media , 2015 .

[66]  Muhammad Sahimi,et al.  Flow and Transport in Porous Media and Fractured Rock: From Classical Methods to Modern Approaches , 1995 .

[67]  German Engineering.  New benchmark. , 2007, Nursing management.

[68]  Joachim Moortgat,et al.  Adaptive implicit finite element methods for multicomponent compressible flow in heterogeneous and fractured porous media , 2017 .

[69]  H. Suk.  Modified Mixed Lagrangian–Eulerian Method Based on Numerical Framework of MT3DMS on Cauchy Boundary , 2016, Ground water.

[70]  Dong Liang,et al.  An ELLAM Approximation for Highly Compressible Multicomponent Flows in Porous Media , 2002 .

[71]  M. Belayneh,et al.  Finite Element - Node-Centered Finite-Volume Two-Phase-Flow Experiments With Fractured Rock Represented by Unstructured Hybrid-Element Meshes , 2007 .

[72]  Chaoshui Xu,et al.  A new computer code for discrete fracture network modelling , 2010, Comput. Geosci..

[73]  S. P. Neuman,et al.  Trends, prospects and challenges in quantifying flow and transport through fractured rocks , 2005 .

[74]  Thomas Graf,et al.  Fracture network optimization for simulating 2D variable-density flow and transport , 2015 .

[75]  Jean-Raynald de Dreuzy,et al.  Random Walk Methods for Modeling Hydrodynamic Transport in Porous and Fractured Media from Pore to Reservoir Scale , 2016, Transport in Porous Media.

[76]  J. Bear,et al.  An adaptive pathline-based particle tracking algorithm for the Eulerian–Lagrangian method , 2000 .

[77]  André Fourno,et al.  A new estimation of equivalent matrix block sizes in fractured media with two-phase flow applications in dual porosity models , 2017 .

[78]  Anis Younes,et al.  Hybrid and multi-point formulations of the lowest-order mixed methods for Darcy's flow on triangles , 2008 .

[79]  A. Younes,et al.  Solving density driven flow problems with efficient spatial discretizations and higher-order time integration methods , 2009 .

[80]  Wolfgang Nowak,et al.  A method for implementing Dirichlet and third‐type boundary conditions in PTRW simulations , 2014 .

[81]  Carl Tim Kelley,et al.  Numerical simulation of water resources problems: Models, methods, and trends , 2013 .

[82]  Frederick Delay,et al.  Flow in double-porosity aquifers: Parameter estimation using an adaptive multiscale method , 2014 .