On direct inverse of Stokes, Helmholtz and Laplacian operators in view of time-stepper-based Newton and Arnoldi solvers in incompressible CFD

Abstract

Factorization of the incompressible Stokes operator linking pressure and velocity is revisited. The main purpose is to use the inverse of the Stokes operator with a large time step as a preconditioner for Newton and Arnoldi iterations applied to computation of steady three-dimensional flows and to study of their stability. It is shown that the Stokes operator can be inversed within an acceptable computational effort. This inverse includes fast direct inverses of several Helmholtz operators and iterative inverse of the pressure matrix. It is shown, additionally, that fast direct solvers can be attractive for the inverse of the Helmholtz and Laplace operators on fine grids and at large Reynolds numbers, as well as for other problems where convergence of iterative methods slows down. Implementation of the Stokes operator inverse to time-stepping-based formulation of the Newton and Arnoldi iterations is discussed.

Cited by in corpus (4)

• Linear instability of the lid-driven flow in a cubic cavity
• Semi-implicit direct forcing immersed boundary method for incompressible viscous thermal flow problems: a Schur complement approach
• Computational modeling of magnetoconvection: effects of discretization method, grid refinement and grid stretching
• Non-modal disturbances growth in a viscous mixing layer flow