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

Research output: Contribution to journalArticle

27 Citations (Scopus)

Abstract

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 Léja 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. © 2010 Elsevier Inc. All rights reserved.

Original languageEnglish
Pages (from-to)3957-3969
Number of pages13
JournalJournal of Computational Physics
Volume229
Issue number10
DOIs
Publication statusPublished - 20 Apr 2010

Keywords

  • Advection-diffusion equation
  • Exponential integration
  • Fast time integrators
  • Krylov subspace
  • Léja points
  • Porous media

Fingerprint Dive into the research topics of 'An exponential integrator for advection-dominated reactive transport in heterogeneous porous media'. Together they form a unique fingerprint.

Cite this