{"title":"Global accelerated Hermitian and skew–Hermitian splitting preconditioner for the solution of discrete Stokes problems","authors":"A. Badahmane , A. Ratnani , H. Sadok","doi":"10.1016/j.cam.2025.116620","DOIUrl":null,"url":null,"abstract":"<div><div>In fluid mechanics, numerous applications necessitate solving a sequence of linear systems. These systems typically arise from discretizing the Stokes equations using mixed-finite element methods. The matrices that result from this process often exhibit a saddle point structure, which makes the iterative solution of preconditioned linear systems challenging. To effectively solve the large scale and ill-conditioned linear systems, it is necessary to implement efficient linear solvers. We present a global approach, specifically the preconditioned global conjugate gradient (PGCG), aimed at enhancing the performance of preconditioned Hermitian and skew–Hermitian splitting preconditioner <span><math><mrow><mo>(</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>P</mi><mi>H</mi><mi>S</mi><mi>S</mi></mrow></msub><mo>)</mo></mrow></math></span> and accelerated Hermitian and skew–Hermitian splitting preconditioner <span><math><mrow><mo>(</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>A</mi><mi>H</mi><mi>S</mi><mi>S</mi></mrow></msub><mo>)</mo></mrow></math></span>. The new preconditioners can be utilized to expedite the convergence of the generalized minimal residual (GMRES) method. We evaluate the effectiveness of the preconditioned iterative methods by considering the Central Processing Unit (<span><math><mi>CPU</mi></math></span>) times and numbers of the <span><math><mi>P</mi></math></span>GMRES iterations. The numerical results indicate that incorporating <span><math><mi>P</mi></math></span>GCG method improves the performance of the preconditioners <span><math><msub><mrow><mi>P</mi></mrow><mrow><mi>P</mi><mi>H</mi><mi>S</mi><mi>S</mi></mrow></msub></math></span> and <span><math><msub><mrow><mi>P</mi></mrow><mrow><mi>A</mi><mi>H</mi><mi>S</mi><mi>S</mi></mrow></msub></math></span>.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"467 ","pages":"Article 116620"},"PeriodicalIF":2.1000,"publicationDate":"2025-03-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational and Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0377042725001359","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
In fluid mechanics, numerous applications necessitate solving a sequence of linear systems. These systems typically arise from discretizing the Stokes equations using mixed-finite element methods. The matrices that result from this process often exhibit a saddle point structure, which makes the iterative solution of preconditioned linear systems challenging. To effectively solve the large scale and ill-conditioned linear systems, it is necessary to implement efficient linear solvers. We present a global approach, specifically the preconditioned global conjugate gradient (PGCG), aimed at enhancing the performance of preconditioned Hermitian and skew–Hermitian splitting preconditioner and accelerated Hermitian and skew–Hermitian splitting preconditioner . The new preconditioners can be utilized to expedite the convergence of the generalized minimal residual (GMRES) method. We evaluate the effectiveness of the preconditioned iterative methods by considering the Central Processing Unit () times and numbers of the GMRES iterations. The numerical results indicate that incorporating GCG method improves the performance of the preconditioners and .
期刊介绍:
The Journal of Computational and Applied Mathematics publishes original papers of high scientific value in all areas of computational and applied mathematics. The main interest of the Journal is in papers that describe and analyze new computational techniques for solving scientific or engineering problems. Also the improved analysis, including the effectiveness and applicability, of existing methods and algorithms is of importance. The computational efficiency (e.g. the convergence, stability, accuracy, ...) should be proved and illustrated by nontrivial numerical examples. Papers describing only variants of existing methods, without adding significant new computational properties are not of interest.
The audience consists of: applied mathematicians, numerical analysts, computational scientists and engineers.