A new alternating positive semidefinite splitting preconditioner for saddle point problems from time-harmonic eddy current models

Abstract

Based on the special positive semidefinite splittings of the saddle point matrix, we propose a new alternating positive semidefinite splitting (APSS) iteration method for the saddle point problem arising from the finite element discretization of the hybrid formulation of the time-harmonic eddy current problem. We prove that the new APSS iteration method is unconditionally convergent for both cases of the simple topology and the general topology. The new APSS matrix can be used as a preconditioner to accelerate the convergence rate of Krylov subspace methods. Numerical results show that the new APSS preconditioner is superior to the existing preconditioners.

Keywords

Time-harmonic eddy current problem / saddle point problem / alternating positive semidefinite splitting (APSS) / convergence analysis / preconditioner / iteration method

References

[1] Arioli M, Manzini G. A null space algorithm for mixed finite-element approximations of Darcy’s equation. Comm Numer Methods Engrg, 2002, 18: 645–657

[2] Bai Z Z. Optimal parameters in the HSS-like methods for saddle-point problems. Numer Linear Algebra Appl, 2009, 16: 447–479

[3] Bai Z Z. Block alternating splitting implicit iteration methods for saddle-point problems from time-harmonic eddy current models. Numer Linear Algebra Appl, 2012, 19: 914–936

[4] Bai Z Z, Golub G H. Accelerated Hermitian and skew-Hermitian splitting iteration methods for saddle-point problems. IMA J Numer Anal, 2007, 27: 1–23

[5] 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, 2003, 24: 603–626

[6] 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, 2004, 98: 1–32

[7] Bai Z Z, Parlett B N, Wang Z Q. On generalized successive overrelaxation methods for augmented linear systems. Numer Math, 2005, 102: 1–38

[8] Bai Z Z, Wang Z Q. On parameterized inexact Uzawa methods for generalized saddle point problems. Linear Algebra Appl, 2008, 428: 2900–2932

[9] Benzi M. Solution of equality-constrained quadratic programming problems by a projection iterative method. Rend Mat Appl, 1993, 13: 275–296

[10] Benzi M, Golub G H. A preconditioner for generalized saddle point problems. SIAM J Matrix Anal Appl, 2004, 26: 20–41

[11] Benzi M, Golub G H, Liesen J. Numerical solution of saddle point problems. Acta Numer, 2005, 14: 1–137

[12] Brezzi F, Fortin M. Mixed and Hybrid Finite Element Methods. Springer Ser Comput Math, Vol 15. New York: Springer-Verlag, 1991

[13] Cao Y, Dong J L, Yu Y M. A relaxed deteriorated PSS preconditioner for nonsymmetric saddle point problems from the steady Navier-Stokes equation. J Comput Appl Math, 2015, 273: 41–60

[14] Cao Y, Du J, Niu Q. Shift-splitting preconditioners for saddle point problems. J Comput Appl Math, 2014, 272: 239–250

[15] Chen C R, Ma C F. A generalized shift-splitting preconditioner for singular saddle point problems. Appl Math Comput, 2015, 269: 947–955

[16] Gould N, Orban D, Rees T. Projected Krylov methods for saddle-point systems. SIAM J Matrix Anal Appl, 2014, 35: 1329–1343

[17] Ke Y F, Ma C F. Spectrum analysis of a more general augmentation block preconditioner for generalized saddle point matrices. BIT, 2016, 56: 489–500

[18] Kotiuga P. Topological considerations in coupling magnetic scalar potentials to stream functions describing surface currents. IEEE Trans Magn, 1989, 25: 2925–2927

[19] Krukier L A, Krukier B L, Ren Z R. Generalized skew-Hermitian triangular splitting iteration methods for saddle-point linear systems. Numer Linear Algebra Appl, 2014, 21: 152–170

[20] Ren Z R, Cao Y. An alternating positive semidefinite splitting preconditioner for saddle point problems from time-harmonic eddy current models. IMA J Numer Anal, 2016, 36: 922–946

[21] Rodríguez A, Herńandez R. Iterative methods for the saddle-point problem arising from the HC/EI formulation of the eddy current problem. SIAM J Sci Comput, 2009, 31: 3155–3178

[22] Rodríguez A, Hiptmair R, Valli A. A hybrid formulation of eddy current problems. Numer Methods Partial Differential Equations, 2005, 21: 742–763

[23] Rodríguez A, Valli A. Eddy Current Approximation of Maxwell Equations: Theory, Algorithms and Applications. Milan: Springer, 2010

[24] Van der Vorst H A. Iterative Krylov Methods for Large Linear Systems. Cambridge: Cambridge Univ Press, 2003

[25] Zhang G F, Ren Z R, Zhou Y Y. On HSS-based constraint preconditioners for generalized saddle-point problems. Numer Algorithms, 2011, 57: 273–287

RIGHTS & PERMISSIONS

Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature