arrow
Volume 38, Issue 1
Local Pressure Correction for the Stokes System

Malte Braack & Utku Kaya

J. Comp. Math., 38 (2020), pp. 125-141.

Published online: 2020-02

Export citation
  • Abstract

Pressure correction methods constitute the most widely used solvers for the time-dependent Navier-Stokes equations. There are several different pressure correction methods, where each time step usually consists in a predictor step for a non-divergence-free velocity, followed by a Poisson problem for the pressure (or pressure update), and a final velocity correction to obtain a divergence-free vector field. In some situations, the equations for the velocities are solved explicitly, so that the most numerically expensive step is the elliptic pressure problem. We here propose to solve this Poisson problem by a domain decomposition method which does not need any communication between the sub-regions. Hence, this system is perfectly adapted for parallel computation. We show under certain assumptions that this new scheme has the same order of convergence as the original pressure correction scheme (with  global projection). Numerical examples for the Stokes system show the effectivity of this new pressure correction method. The convergence order $\mathcal O(k^{ 2})$ for resulting velocity fields can be observed in the norm $l^2(0,T ;L^2(\Omega))$.

  • AMS Subject Headings

76D07, 65M55

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address

braack@math.uni-kiel.de (Malte Braack)

kaya@math.uni-kiel.de (Utku Kaya)

  • BibTex
  • RIS
  • TXT
@Article{JCM-38-125, author = {Braack , Malte and Kaya , Utku}, title = {Local Pressure Correction for the Stokes System}, journal = {Journal of Computational Mathematics}, year = {2020}, volume = {38}, number = {1}, pages = {125--141}, abstract = {

Pressure correction methods constitute the most widely used solvers for the time-dependent Navier-Stokes equations. There are several different pressure correction methods, where each time step usually consists in a predictor step for a non-divergence-free velocity, followed by a Poisson problem for the pressure (or pressure update), and a final velocity correction to obtain a divergence-free vector field. In some situations, the equations for the velocities are solved explicitly, so that the most numerically expensive step is the elliptic pressure problem. We here propose to solve this Poisson problem by a domain decomposition method which does not need any communication between the sub-regions. Hence, this system is perfectly adapted for parallel computation. We show under certain assumptions that this new scheme has the same order of convergence as the original pressure correction scheme (with  global projection). Numerical examples for the Stokes system show the effectivity of this new pressure correction method. The convergence order $\mathcal O(k^{ 2})$ for resulting velocity fields can be observed in the norm $l^2(0,T ;L^2(\Omega))$.

}, issn = {1991-7139}, doi = {https://doi.org/10.4208/jcm.1906-m2018-0210}, url = {http://global-sci.org/intro/article_detail/jcm/13688.html} }
TY - JOUR T1 - Local Pressure Correction for the Stokes System AU - Braack , Malte AU - Kaya , Utku JO - Journal of Computational Mathematics VL - 1 SP - 125 EP - 141 PY - 2020 DA - 2020/02 SN - 38 DO - http://doi.org/10.4208/jcm.1906-m2018-0210 UR - https://global-sci.org/intro/article_detail/jcm/13688.html KW - Stokes system, Navier-Stokes, Pressure correction, Finite elements. AB -

Pressure correction methods constitute the most widely used solvers for the time-dependent Navier-Stokes equations. There are several different pressure correction methods, where each time step usually consists in a predictor step for a non-divergence-free velocity, followed by a Poisson problem for the pressure (or pressure update), and a final velocity correction to obtain a divergence-free vector field. In some situations, the equations for the velocities are solved explicitly, so that the most numerically expensive step is the elliptic pressure problem. We here propose to solve this Poisson problem by a domain decomposition method which does not need any communication between the sub-regions. Hence, this system is perfectly adapted for parallel computation. We show under certain assumptions that this new scheme has the same order of convergence as the original pressure correction scheme (with  global projection). Numerical examples for the Stokes system show the effectivity of this new pressure correction method. The convergence order $\mathcal O(k^{ 2})$ for resulting velocity fields can be observed in the norm $l^2(0,T ;L^2(\Omega))$.

Malte Braack & Utku Kaya. (2020). Local Pressure Correction for the Stokes System. Journal of Computational Mathematics. 38 (1). 125-141. doi:10.4208/jcm.1906-m2018-0210
Copy to clipboard
The citation has been copied to your clipboard