Preconditioners for two-phase incompressible Navier-Stokes flow

Niall Bootland, Alistair Bentley, Christopher Kees, Andrew Wathen

Research output: Contribution to journalArticle

Abstract

We consider iterative methods for solving the linearized Navier–Stokes equations arising from two-phase flow problems and the efficient preconditioning of such systems when using mixed finite element methods. Our target application is simulation within the Proteus toolkit; in particular, we will give results for a dynamic dam-break problem in two dimensions. We focus on a preconditioner motivated by approximate commutators which has proved effective, displaying mesh-independent convergence for the constant coefficient single-phase Navier–Stokes equations. This approach is known as the “pressure convection–diffusion” (PCD) preconditioner [H. C. Elman, D. J. Silvester, and A. J. Wathen, Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics, 2nd ed., Oxford University Press, Oxford, 2014]. However, the original technique fails to give comparable performance in its given form when applied to variable coefficient Navier–Stokes systems such as those arising in two-phase flow models. Here we develop a generalization of this preconditioner appropriate for two-phase flow, requiring a new form for PCD. We omit considerations of boundary conditions to focus on the key features of two-phase flow. Before considering our target application, we present numerical results within the controlled setting of a simplified problem using a variety of different mixed elements. We compare these results with those for a straightforward extension to another commutator-based method known as the “least-squares commutator” (LSC) preconditioner, a technique also discussed in the aforementioned reference. We demonstrate that favorable properties of the original PCD and LSC preconditioners (without boundary adjustments) are retained with the new preconditioners in the two-phase situation.

LanguageEnglish
PagesB843–B869
Number of pages27
JournalSIAM Journal on Scientific Computing
Volume41
Issue number4
DOIs
Publication statusPublished - 22 Aug 2019

Fingerprint

Incompressible Navier-Stokes
Electric commutators
Stokes Flow
Two phase flow
Preconditioner
Two-phase Flow
Commutator
Convection-diffusion
Navier Stokes equations
Least Squares
Navier-Stokes Equations
Fluid dynamics
Iterative methods
Mixed Element
Dams
Dam Break
Iterative Solvers
Target
Navier-Stokes System
Mixed Finite Element Method

Keywords

  • preconditioner
  • two-phase flow
  • Navier-Stokes
  • Schur complement
  • finite elements

Cite this

Bootland, Niall ; Bentley, Alistair ; Kees, Christopher ; Wathen, Andrew. / Preconditioners for two-phase incompressible Navier-Stokes flow. In: SIAM Journal on Scientific Computing. 2019 ; Vol. 41, No. 4. pp. B843–B869.
@article{481d02b40fc74cac93dec2a35f80e770,
title = "Preconditioners for two-phase incompressible Navier-Stokes flow",
abstract = "We consider iterative methods for solving the linearized Navier–Stokes equations arising from two-phase flow problems and the efficient preconditioning of such systems when using mixed finite element methods. Our target application is simulation within the Proteus toolkit; in particular, we will give results for a dynamic dam-break problem in two dimensions. We focus on a preconditioner motivated by approximate commutators which has proved effective, displaying mesh-independent convergence for the constant coefficient single-phase Navier–Stokes equations. This approach is known as the “pressure convection–diffusion” (PCD) preconditioner [H. C. Elman, D. J. Silvester, and A. J. Wathen, Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics, 2nd ed., Oxford University Press, Oxford, 2014]. However, the original technique fails to give comparable performance in its given form when applied to variable coefficient Navier–Stokes systems such as those arising in two-phase flow models. Here we develop a generalization of this preconditioner appropriate for two-phase flow, requiring a new form for PCD. We omit considerations of boundary conditions to focus on the key features of two-phase flow. Before considering our target application, we present numerical results within the controlled setting of a simplified problem using a variety of different mixed elements. We compare these results with those for a straightforward extension to another commutator-based method known as the “least-squares commutator” (LSC) preconditioner, a technique also discussed in the aforementioned reference. We demonstrate that favorable properties of the original PCD and LSC preconditioners (without boundary adjustments) are retained with the new preconditioners in the two-phase situation.",
keywords = "preconditioner, two-phase flow, Navier-Stokes, Schur complement, finite elements",
author = "Niall Bootland and Alistair Bentley and Christopher Kees and Andrew Wathen",
year = "2019",
month = "8",
day = "22",
doi = "10.1137/17M1153674",
language = "English",
volume = "41",
pages = "B843–B869",
journal = "SIAM Journal on Scientific Computing",
issn = "1064-8275",
number = "4",

}

Preconditioners for two-phase incompressible Navier-Stokes flow. / Bootland, Niall; Bentley, Alistair; Kees, Christopher; Wathen, Andrew.

In: SIAM Journal on Scientific Computing, Vol. 41, No. 4, 22.08.2019, p. B843–B869.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Preconditioners for two-phase incompressible Navier-Stokes flow

AU - Bootland, Niall

AU - Bentley, Alistair

AU - Kees, Christopher

AU - Wathen, Andrew

PY - 2019/8/22

Y1 - 2019/8/22

N2 - We consider iterative methods for solving the linearized Navier–Stokes equations arising from two-phase flow problems and the efficient preconditioning of such systems when using mixed finite element methods. Our target application is simulation within the Proteus toolkit; in particular, we will give results for a dynamic dam-break problem in two dimensions. We focus on a preconditioner motivated by approximate commutators which has proved effective, displaying mesh-independent convergence for the constant coefficient single-phase Navier–Stokes equations. This approach is known as the “pressure convection–diffusion” (PCD) preconditioner [H. C. Elman, D. J. Silvester, and A. J. Wathen, Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics, 2nd ed., Oxford University Press, Oxford, 2014]. However, the original technique fails to give comparable performance in its given form when applied to variable coefficient Navier–Stokes systems such as those arising in two-phase flow models. Here we develop a generalization of this preconditioner appropriate for two-phase flow, requiring a new form for PCD. We omit considerations of boundary conditions to focus on the key features of two-phase flow. Before considering our target application, we present numerical results within the controlled setting of a simplified problem using a variety of different mixed elements. We compare these results with those for a straightforward extension to another commutator-based method known as the “least-squares commutator” (LSC) preconditioner, a technique also discussed in the aforementioned reference. We demonstrate that favorable properties of the original PCD and LSC preconditioners (without boundary adjustments) are retained with the new preconditioners in the two-phase situation.

AB - We consider iterative methods for solving the linearized Navier–Stokes equations arising from two-phase flow problems and the efficient preconditioning of such systems when using mixed finite element methods. Our target application is simulation within the Proteus toolkit; in particular, we will give results for a dynamic dam-break problem in two dimensions. We focus on a preconditioner motivated by approximate commutators which has proved effective, displaying mesh-independent convergence for the constant coefficient single-phase Navier–Stokes equations. This approach is known as the “pressure convection–diffusion” (PCD) preconditioner [H. C. Elman, D. J. Silvester, and A. J. Wathen, Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics, 2nd ed., Oxford University Press, Oxford, 2014]. However, the original technique fails to give comparable performance in its given form when applied to variable coefficient Navier–Stokes systems such as those arising in two-phase flow models. Here we develop a generalization of this preconditioner appropriate for two-phase flow, requiring a new form for PCD. We omit considerations of boundary conditions to focus on the key features of two-phase flow. Before considering our target application, we present numerical results within the controlled setting of a simplified problem using a variety of different mixed elements. We compare these results with those for a straightforward extension to another commutator-based method known as the “least-squares commutator” (LSC) preconditioner, a technique also discussed in the aforementioned reference. We demonstrate that favorable properties of the original PCD and LSC preconditioners (without boundary adjustments) are retained with the new preconditioners in the two-phase situation.

KW - preconditioner

KW - two-phase flow

KW - Navier-Stokes

KW - Schur complement

KW - finite elements

UR - https://epubs.siam.org/loi/sisc

U2 - 10.1137/17M1153674

DO - 10.1137/17M1153674

M3 - Article

VL - 41

SP - B843–B869

JO - SIAM Journal on Scientific Computing

T2 - SIAM Journal on Scientific Computing

JF - SIAM Journal on Scientific Computing

SN - 1064-8275

IS - 4

ER -