Novel Approach for Solving the Discrete Stokes Problems Based on Augmented Lagrangian and Global Techniques with Applications for Stokes Problems

doi:10.1007/s42967-025-00488-0

Abstract

Abstract: In this paper, a novel augmented Lagrangian preconditioner based on the global Arnoldi for accelerating the convergence of Krylov subspace methods is applied to linear systems of equations with a block three-by-three structure, and these systems typically arise from discretizing the Stokes equations using mixed finite-element methods. Spectral analyses are established for the exact versions of these preconditioners. Finally, the obtained numerical results claim that our novel approach is more efficient and robust for solving the discrete Stokes problems. The efficiency of our new approach is evaluated by measuring the computational time.

Key words: Stokes equation, Saddle point problem, Krylov subspace method, Global Krylov subspace method, Augmented Lagrangian-based preconditioning

CLC Number:

65F10
65F08

A. Badahmane, A. Ratnani, H. Sadok. Novel Approach for Solving the Discrete Stokes Problems Based on Augmented Lagrangian and Global Techniques with Applications for Stokes Problems[J]. Communications on Applied Mathematics and Computation, 2026, 8(3): 1076-1094.

TrendMD

References

[1] Arnoldi, W.E.: The principle of minimized iterations in the solution of the matrix eigenvalue problem. Quart. Appl. Math. 9, 17-29 (1951)
[2] Badahmane, A., Bentbib, A.H., Sadok, H.: Preconditioned global Krylov subspace methods for solving saddle point problems with multiple right-hand sides. Electron. Trans. Numer. Anal. 51, 495-511 (2019)
[3] Bai, Z.-Z., Golub, G.H.: Accelerated Hermitian and skew-Hermitian splitting iteration methods for saddle-point problems. IMA J. Numer. Anal. 27, 1-23 (2007)
[4] Bai, Z.-Z., Golub, G.H., Ng, M.K.: Hermitian and skew-Hermitian splitting methods for non-Hermitian positive definite linear systems. SIAM J. Matrix Anal. Appl. 24, 603-626 (2003)
[5] Bai, Z.-Z., Golub, G.H., Pan, J.-Y.: Preconditioned Hermitian and skew-Hermitian splitting methods for non-Hermitian positive semidefinite linear systems. Numer. Math. 98, 1-32 (2004)
[6] Bai, Z.-Z., Parlett, B.N., Wang, Z.-Q.: On generalized successive overrelaxation methods for augmented linear systems. Numer. Math. 102, 1-38 (2005)
[7] Bai, Z.-Z., Wang, Z.-Q.: On parameterized inexact Uzawa methods for generalized saddle point problems. Linear Algebra Appl. 428, 2900-2932 (2008)
[8] Benzi, M., Golub, G.H., Liesen, J.: Numerical solution of saddle point problems. Acta Numer. 14, 1-137 (2005)
[9] Bramble, J.H., Pasciak, J.E., Vassilev, A.T.: Analysis of the inexact Uzawa algorithm for saddle point problems. SIAM J. Numer. Anal. 34, 1072-1092 (1997)
[10] Ebadi, G., Alipour, N., Vuik, C.: Deflated and augmented global Krylov subspace methods for the matrix equations. Appl. Numer. Math. 99, 137-150 (2016)
[11] Elman, H.C., Silvester, D.J., Wathen, A.J.: Finite Elements and Fast Iterative Solvers with Applications in Incompressible Fluid Dynamics. Oxford University Press, New York (2005)
[12] Golub, G.H., Wu, X., Yuan, J.-Y.: SOR-like methods for augmented systems. BIT 41, 71-85 (2001)
[13] Guo, P., Li, C.-X., Wu, S.-L.: A modified SOR-like method for the augmented systems. J. Comput. Appl. Math. 274, 58-69 (2015)
[14] Jbilou, K., Sadok, H., Tinzefte, A.: Oblique projection methods for linear systems with multiple right-hand sides. Electron. Trans. Numer. Anal. 20, 119-138 (2005)
[15] Saad, Y.: A flexible inner-outer preconditioned GMRES algorithm. SIAM J. Sci. Comput. 14, 461-469 (1993)
[16] Saad, Y.: Iterative Methods for Sparse Linear Systems. SIAM, Philadelphia, PA (2003)
[17] Saad, Y., Schultz, M.H.: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comp. 7, 856-869 (1986)
[18] Wang, N.-N., Li, J.-C.: A class of new extended shift-splitting preconditioners for saddle point problems. J. Comput. Appl. Math. 357, 123-145 (2019)
[19] Zhang, J.-J., Shang, J.-J.: A class of Uzawa-SOR methods for saddle point problems. Appl. Math. Comput. 216, 2163-2168 (2010)
[20] Zheng, Q.-Q., Ma, C.-F.: Preconditioned AHSS-PU alternating splitting iterative methods for saddle point problems. Comput. Appl. Math. 273, 217-225 (2016)