Abstract
In this paper we study efficient iterative methods for solving the system of linear equations arising from the fully implicit Runge-Kutta discretizations of a class of partial differential-algebraic equations. In each step of the time integration, a block two-by-two linear system is obtained and needed to be solved numerically. A preconditioning strategy based on an alternating Kronecker product splitting of the coefficient matrix is proposed to solve such linear systems. Some spectral properties of the preconditioned matrix are established and numerical examples are presented to demonstrate the effectiveness of this approach.
Similar content being viewed by others
References
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)
Benzi, M., Golub, G.H.: A preconditioner for generalized saddle point problems. SIAM J. Matrix. Anal. Appl. 26, 20–41 (2004)
Berman, A., Plemmons, R.J.: Nonnegative matrices in the Mathematical Sciences. SIAM, Philadelphia (1994)
Briggs, W., Henson, V.E., McCormick, S.F.: A multigrid tutorial. SIAM, Philadelphia (2000)
Chen, H.: A splitting preconditioner for the iterative solution of implicit Runge-Kutta and boundary value methods. BIT 54, 607–621 (2014)
Chen, H.: Generalized Kronecker product splitting iteration for the solution of implicit Runge-Kutta and boundary value methods. Numer. Linear. Algebra. Appl. 22, 357–370 (2015)
Dax, A.: The convergence of linear stationary iterative processes for solving singular unstructured systems of linear equations. SIAM Rev. 32, 611–635 (1990)
Colli Franzone, P., Deuflhard, P., Erdmann, B., Lang, J., Pavarino, L.F.: Adaptivity in space and time for reaction-diffusion systems in electrocardiology. SIAM J. Sci. Comput. 28, 942–962 (2006)
Elman, H.C., Ramage, A., Silvester, D.J.: IFISS: a Matlab toolbox for modelling incompressible flow. ACM Trans. Math. Software 33, Article 14 (2007)
Elman, H.C., Silvester, D.J., Wathen, A.: Finite elements and fast iterative solvers with applications in incompressible fluid dynamics. Numer. Math. Sci. Comput. Oxford University Press, Oxford (2005)
Ethier, M., Bourgault, Y.: Semi-implicit time-discretization schemes for the Bidomain model. SIAM J. Numer. Anal. 46, 2443–2468 (2008)
Gerardo-Giorda, L., Mirabella, L., Nobile, F., Perego, M., Veneziani, A.: A model-based block-triangular preconditioner for the Bidomain system in electrocardiology. J. Comput. Phys. 228, 3625–3639 (2009)
Gerardo-Giorda, L., Mirabella, L.: Spectral analysis of a block-triangular preconditioner for the Bidomain system in electrocardiology. Electron. Trans. Numer. Anal. 39, 186–201 (2012)
Hairer, E., Wanner, G.: Solving ordinary differential equations II. Stiff and differential algebraic problems. Springer, Berlin (1996)
Horn, R.A., Johnson, C.R.: Topics in matrix analysis. Cambridge University Press, Cambridge (1991)
Keener, J., Sneyd, J.: Mathematical physiology, 2nd ed., vol. I-II. Springer, New York (2009)
Mardal, K.A., Nielsen, B.F., Cai, X., Tveito, A.: An order optimal solver for the discretized bidomain equations. Numer. Linear. Algebra. Appl. 14, 83–98 (2007)
Marsh, M.E., Ziaratqahi, S.T., Spiteri, R.J.: The secrets to the success of the Rush-Larsen method and its generalizations. IEEE Trans. Biomed. Eng. 59, 2506–2515 (2012)
Nilssen, T.K., Staff, G.A., Mardal, K.A.: Order-optimal preconditioners for fully implicit Runge-Kutta schemes applied to the bidomain equations. Numer. Meth. Part. Diff. Equ. 27, 1290–1312 (2011)
Pavarino, L.F., Scacchi, S.: Multilevel additive Schwarz preconditioners for the Bidomain reaction-diffusion system. SIAM J. Sci. Comput. 31, 420–443 (2008)
Pavarino, L.F., Scacchi, S.: Parallel multilevel Schwarz and block preconditioners for the Bidomain parabolic-parobolic and parabolic-elliptic formulations. SIAM J. Sci. Comput. 33, 1897–1919 (2011)
Pennacchio, M., Simoncini, V.: Efficient algebraic solution of reaction-diffusion systems for the cardiac excitation process. J. Comput. Appl. Math. 145, 49–70 (2002)
Pennacchio, M., Simoncini, V.: Algebraic multigrid preconditioners for the bidomain reaction-diffusion system. Appl. Numer. Math. 59, 3033–3050 (2009)
Pennacchio, M., Simoncini, V.: Fast structured AMG preconditioning for the bidomain model in electrocadiology. SIAM J. Sci. Comput. 33, 721–745 (2011)
Perego, M., Veneziani, A.: An efficient generalization of the Rush-Larsen method for solving electro-physiology membrane equations. Electron. Trans. Numer. Anal. 35, 234–256 (2009)
Plank, G., Liebmann, M., Weber Dos Santos, R., Vigmond, E.J., Haase, G.: Algebraic multigrid preconditioner for the cardiac Bidomain model. IEEE Trans. Biomed. Engrg. 54, 585–596 (2007)
Ruge, J.W., Stüben, K.: Algebraic multigrid. In: McCormick, S.F. (ed.) Multigrid methods, Frontiers Appl. Math., vol. 3, pp 73–130. SIAM, Philadelphia (1987)
Saad, Y., Schultz, M.H.: GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Statist. Comput. 7, 856–869 (1986)
Simoncini, V., Benzi, M.: Spectral properties of the Hermitian and skew-Hermitian splitting preconditioner for saddle point problems. SIAM J. Matrix Anal. Appl. 26, 377–389 (2004)
Strang, G.: On the construction and comparison of difference schemes. SIAM J. Numer. Anal. 5, 506–517 (1968)
Sundnes, J., Lines, G.T., Mardal, K.A., Tveito, A.: Multigrid block preconditioning for a coupled system of partial differential equations modeling the electrical activity in the heart. Comput. Methods Biomech. Biomed. Engrg. 5, 397–409 (2002)
Sundnes, J., Lines, G.T., Tveito, A.: An operator splitting method for solving the Bidomain equations coupled to a volume conductor model for the torso. Comput. Math. Biosci. 194, 233–248 (2005)
Varga, R.S.: Matrix iterative analysis. Prentice-Hall, Englewood Cliffs (1962)
Vigmond, E.J., Weber dos Santos, R., Prassl, A.J., Deo, M., Plank, G.: Solvers for the cardiac bidomain equations. Progr. Biophys. Molecular Biol. 96, 3–18 (2008)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Chen, H. A splitting preconditioner for implicit Runge-Kutta discretizations of a partial differential-algebraic equation. Numer Algor 73, 1037–1054 (2016). https://doi.org/10.1007/s11075-016-0128-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-016-0128-5
Keywords
- Iterative methods
- Preconditioning
- Implicit Runge-Kutta methods
- Differential-algebraic equation
- Bidomain equations