Unified Convergence Analysis of Numerical Schemes for a Miscible Displacement Problem

Jérôme Droniou, Robert Eymard, Alain Prignet, Kyle S. Talbot

Research output: Contribution to journalArticleResearchpeer-review

Abstract

This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in (Formula presented.) of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes is compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion.

Original languageEnglish
Pages (from-to)333-374
Number of pages42
JournalFoundations of Computational Mathematics
Volume19
Issue number2
DOIs
Publication statusPublished - 2019

Keywords

  • Convergence analysis
  • Coupled elliptic–parabolic problem
  • Finite differences
  • Gradient discretisation method
  • Mass-lumped finite elements
  • Miscible fluid flow
  • Uniform-in-time convergence

Cite this

Droniou, Jérôme ; Eymard, Robert ; Prignet, Alain ; Talbot, Kyle S. / Unified Convergence Analysis of Numerical Schemes for a Miscible Displacement Problem. In: Foundations of Computational Mathematics. 2019 ; Vol. 19, No. 2. pp. 333-374.
@article{13430f1ff7374f478fd567a58a9696eb,
title = "Unified Convergence Analysis of Numerical Schemes for a Miscible Displacement Problem",
abstract = "This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in (Formula presented.) of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes is compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion.",
keywords = "Convergence analysis, Coupled elliptic–parabolic problem, Finite differences, Gradient discretisation method, Mass-lumped finite elements, Miscible fluid flow, Uniform-in-time convergence",
author = "J{\'e}r{\^o}me Droniou and Robert Eymard and Alain Prignet and Talbot, {Kyle S.}",
year = "2019",
doi = "10.1007/s10208-018-9387-y",
language = "English",
volume = "19",
pages = "333--374",
journal = "Foundations of Computational Mathematics",
issn = "1615-3375",
publisher = "Springer-Verlag London Ltd.",
number = "2",

}

Unified Convergence Analysis of Numerical Schemes for a Miscible Displacement Problem. / Droniou, Jérôme; Eymard, Robert; Prignet, Alain; Talbot, Kyle S.

In: Foundations of Computational Mathematics, Vol. 19, No. 2, 2019, p. 333-374.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Unified Convergence Analysis of Numerical Schemes for a Miscible Displacement Problem

AU - Droniou, Jérôme

AU - Eymard, Robert

AU - Prignet, Alain

AU - Talbot, Kyle S.

PY - 2019

Y1 - 2019

N2 - This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in (Formula presented.) of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes is compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion.

AB - This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in (Formula presented.) of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes is compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion.

KW - Convergence analysis

KW - Coupled elliptic–parabolic problem

KW - Finite differences

KW - Gradient discretisation method

KW - Mass-lumped finite elements

KW - Miscible fluid flow

KW - Uniform-in-time convergence

UR - http://www.scopus.com/inward/record.url?scp=85044362838&partnerID=8YFLogxK

U2 - 10.1007/s10208-018-9387-y

DO - 10.1007/s10208-018-9387-y

M3 - Article

VL - 19

SP - 333

EP - 374

JO - Foundations of Computational Mathematics

JF - Foundations of Computational Mathematics

SN - 1615-3375

IS - 2

ER -