Volume 22, Issue 2

Commun. Comput. Phys., 22 (2017), pp. 494-516.

Published online: 2018-04

Preview Purchase PDF 104 2071
Export citation

Cited by

• Abstract

This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I−∆tL which maps the identity (no preconditioner) for ∆t ≪ 1 to Laplacian preconditioning for ∆t ≫ 1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for ∆t = $\mathcal{O}$(1), away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

• Keywords

• BibTex
• RIS
• TXT
@Article{CiCP-22-494, author = {Cédric Beaume , }, title = {Adaptive Stokes Preconditioning for Steady Incompressible Flows}, journal = {Communications in Computational Physics}, year = {2018}, volume = {22}, number = {2}, pages = {494--516}, abstract = {

This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I−∆tL which maps the identity (no preconditioner) for ∆t ≪ 1 to Laplacian preconditioning for ∆t ≫ 1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for ∆t = $\mathcal{O}$(1), away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2016-0201}, url = {http://global-sci.org/intro/article_detail/cicp/11308.html} }
TY - JOUR T1 - Adaptive Stokes Preconditioning for Steady Incompressible Flows AU - Cédric Beaume , JO - Communications in Computational Physics VL - 2 SP - 494 EP - 516 PY - 2018 DA - 2018/04 SN - 22 DO - http://doi.org/10.4208/cicp.OA-2016-0201 UR - https://global-sci.org/intro/article_detail/cicp/11308.html KW - AB -

This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I−∆tL which maps the identity (no preconditioner) for ∆t ≪ 1 to Laplacian preconditioning for ∆t ≫ 1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for ∆t = $\mathcal{O}$(1), away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

Cédric Beaume. (2020). Adaptive Stokes Preconditioning for Steady Incompressible Flows. Communications in Computational Physics. 22 (2). 494-516. doi:10.4208/cicp.OA-2016-0201
Copy to clipboard
The citation has been copied to your clipboard