An exponential integrator for advection-dominated reactive transport in heterogeneous porous media

We present an exponential time integrator in conjunction with a finite volume discretisation in space for simulating transport by advection and diffusion including chemical reactions in highly heterogeneous porous media representative of geological reservoirs. These numerical integrators are based on the variation of constants solution and solving the linear system exactly. This is at the expense of computing the exponential of the stiff matrix comprising the finite volume discretisation. Using real Leja points or a Krylov subspace technique compared to standard finite difference-based time integrators. We observe for a variety of example applications that numerical solutions with exponential methods are generally more accurate and require less computational cost. They hence comprise an efficient and accurate method for simulating non-linear advection-dominated transport in geological formations.

[1]  Richard D. Hornung,et al.  Adaptive Mesh Refinement and Multilevel Iteration for Flow in Porous Media , 1997 .

[2]  L. Margolin,et al.  On balanced approximations for time integration of multiple time scale systems , 2003 .

[3]  P. Knabner,et al.  Numerical Methods for Elliptic and Parabolic Partial Differential Equations , 2003, Texts in Applied Mathematics.

[4]  H. Tchelepi,et al.  Multi-scale finite-volume method for elliptic problems in subsurface flow simulation , 2003 .

[5]  M. Belayneh,et al.  Fluid flow partitioning between fractures and a permeable rock matrix , 2004 .

[6]  C. Steefel,et al.  Approaches to modeling of reactive transport in porous media , 1996 .

[7]  Todd A. Ehlers,et al.  REVIEWS IN MINERALOGY AND GEOCHEMISTRY , 2005 .

[8]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[9]  J. Shadid,et al.  Studies on the accuracy of time-integration methods for the radiation-diffusion equations , 2004 .

[10]  Marco Vianello,et al.  A massively parallel exponential integrator for advection-diffusion models , 2009, J. Comput. Appl. Math..

[11]  P. Kitanidis,et al.  Characterization of mixing and dilution in heterogeneous aquifers by means of local temporal moments , 2000 .

[12]  R. Eymard,et al.  Finite Volume Methods , 2019, Computational Methods for Fluid Dynamics.

[13]  Lucia Faravelli,et al.  Stochastic finite elements for crash problems , 1990 .

[14]  A. Ostermann,et al.  Implementation of exponential Rosenbrock-type integrators , 2009 .

[15]  Homa Karimabadi,et al.  Self-adaptive time integration of flux-conservative equations with sources , 2006, J. Comput. Phys..

[16]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[17]  Homa Karimabadi,et al.  A new asynchronous methodology for modeling of physical systems: breaking the curse of courant condition , 2005 .

[18]  Witold F. Krajewski,et al.  Modeling radar-rainfall estimation uncertainties using parametric and non-parametric approaches , 2008 .

[19]  Timothy D. Scheibe,et al.  Mixing‐induced precipitation: Experimental study and multiscale numerical analysis , 2008 .

[20]  J. W. Thomas.  Numerical Partial Differential Equations , 1999 .

[21]  Carl I. Steefel,et al.  Reactive transport in porous media , 1996 .

[22]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

[23]  M. Dentz,et al.  Modeling non‐Fickian transport in geological formations as a continuous time random walk , 2006 .

[24]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[25]  J. Shadid,et al.  Studies of the Accuracy of Time Integration Methods for Reaction-Diffusion Equations ∗ , 2005 .

[26]  Timothy D. Scheibe,et al.  EFFECTS OF INCOMPLETE MIXING ON MULTICOMPONENT REACTIVE TRANSPORT , 2009 .

[27]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[28]  Thomas Y. Hou,et al.  Convergence of a Nonconforming Multiscale Finite Element Method , 2000, SIAM J. Numer. Anal..

[29]  John A. Trangenstein,et al.  Adaptive Mesh Re nement and Multilevel Iteration for Flow inPorous , 1997 .

[30]  ARTICLE IN PRESS Journal of Computational and Applied Mathematics ( ) – , 2005 .

[31]  Håvard Berland,et al.  NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET , 2005 .

[32]  J. Knight,et al.  Analytical solution of a spatially variable coefficient advection–diffusion equation in up to three dimensions , 1999 .

[33]  Marco Caliari,et al.  Accurate evaluation of divided differences for polynomial interpolation of exponential propagators , 2007, Computing.

[34]  B. Minchev,et al.  A review of exponential integrators for first order semi-linear problems , 2005 .

[35]  Marco Vianello,et al.  The ReLPM Exponential Integrator for FE Discretizations of Advection-Diffusion Equations , 2004, International Conference on Computational Science.

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

[37]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[38]  L. Bergamaschi,et al.  Interpolating discrete advection-diffusion propagators at Leja sequences , 2004 .

[39]  Marco Vianello,et al.  Comparing Leja and Krylov Approximations of Large Scale Matrix Exponentials , 2006, International Conference on Computational Science.

[40]  Martin J. Blunt,et al.  Streamline‐based dual‐porosity simulation of reactive transport and flow in fractured reservoirs , 2004 .

[41]  Mike Christie,et al.  Uncertainty quantification for porous media flows , 2006, J. Comput. Phys..

[42]  Gene H. Golub,et al.  Matrix computations , 1983 .

[43]  Lloyd N. Trefethen,et al.  Fourth-Order Time-Stepping for Stiff PDEs , 2005, SIAM J. Sci. Comput..

[44]  K. Stüben.  A review of algebraic multigrid , 2001 .

[45]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[46]  M. Blunt,et al.  Simulating Flow in Heterogeneous Systems Using Streamtubes and Streamlines , 1996 .

[47]  Margot Gerritsen,et al.  MODELING FLUID FLOW IN OIL RESERVOIRS , 2005 .