Saddle point problems frequently appear in many mathematical and engineering applications. Most systems of partial differential equations with constraints give rise to saddle point linear systems. Typical examples include mixed finite element formulations to solve fluid flows and/or elasticity problems under full incompressibility. The inversion of saddle point problems is challenging due to inherent numerical instability in the direct inversion methods. Many direct and iterative methods have been proposed to overcome this challenges, such as the Schur complement and the Uzawa’s method. In the context of mixed finite element for incompressible flows using stable H(div)-L2 spaces for velocity and pressure, we propose an iterative method that can effectively solve a saddle point problem iteratively by summing a small compressibility to the original matrix. The preconditioning matrix is symmetric positive definite, which allows the usage of Cholesky decomposition and/or CG-like iterative solvers to compute the incremental solution for the velocities unknowns. A procedure to compute the average pressure of each element of the incompressible problem is developed using the unbalanced fluxes caused by the compressibility perturbation. The average is updated during the iterative process as a function of the velocity increment at each iteration.

Back to Top

Published on 01/07/24

Accepted on 01/07/24

Submitted on 01/07/24

Volume Numerical Methods and Algorithms in Science and Engineering, 2024

DOI: 10.23967/wccm.2024.063

Licence: CC BY-NC-SA license

Are you one of the authors of this document?