A comparison of coupled and uncoupled solvers for the cardiac Bidomain model

P. Colli Franzone; L. F. Pavarino; S. Scacchi

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2013)

  • Volume: 47, Issue: 4, page 1017-1035
  • ISSN: 0764-583X

Abstract

top
The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is coupled through the non-linear reaction term with a stiff system of ordinary differential equations (ODEs), the so-called membrane model, describing the ionic currents through the cellular membrane. A novel uncoupled solver for the Bidomain system is here introduced, based on solving twice the parabolic PDE and once the elliptic PDE at each time step, and it is compared with a usual coupled solver. Three-dimensional numerical tests have been performed in order to show that the proposed uncoupled method has the same accuracy of the coupled strategy. Parallel numerical tests on structured meshes have also shown that the uncoupled technique is as scalable as the coupled one. Moreover, the conjugate gradient method preconditioned by Multilevel Hybrid Schwarz preconditioners converges faster for the linear systems deriving from the uncoupled method than from the coupled one. Finally, in all parallel numerical tests considered, the uncoupled technique proposed is always about two or three times faster than the coupled approach.

How to cite

top

Colli Franzone, P., Pavarino, L. F., and Scacchi, S.. "A comparison of coupled and uncoupled solvers for the cardiac Bidomain model." ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique 47.4 (2013): 1017-1035. <http://eudml.org/doc/273115>.

@article{ColliFranzone2013,
abstract = {The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is coupled through the non-linear reaction term with a stiff system of ordinary differential equations (ODEs), the so-called membrane model, describing the ionic currents through the cellular membrane. A novel uncoupled solver for the Bidomain system is here introduced, based on solving twice the parabolic PDE and once the elliptic PDE at each time step, and it is compared with a usual coupled solver. Three-dimensional numerical tests have been performed in order to show that the proposed uncoupled method has the same accuracy of the coupled strategy. Parallel numerical tests on structured meshes have also shown that the uncoupled technique is as scalable as the coupled one. Moreover, the conjugate gradient method preconditioned by Multilevel Hybrid Schwarz preconditioners converges faster for the linear systems deriving from the uncoupled method than from the coupled one. Finally, in all parallel numerical tests considered, the uncoupled technique proposed is always about two or three times faster than the coupled approach.},
author = {Colli Franzone, P., Pavarino, L. F., Scacchi, S.},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {operator splitting; multilevel preconditioners; parallel computing},
language = {eng},
number = {4},
pages = {1017-1035},
publisher = {EDP-Sciences},
title = {A comparison of coupled and uncoupled solvers for the cardiac Bidomain model},
url = {http://eudml.org/doc/273115},
volume = {47},
year = {2013},
}

TY - JOUR
AU - Colli Franzone, P.
AU - Pavarino, L. F.
AU - Scacchi, S.
TI - A comparison of coupled and uncoupled solvers for the cardiac Bidomain model
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2013
PB - EDP-Sciences
VL - 47
IS - 4
SP - 1017
EP - 1035
AB - The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is coupled through the non-linear reaction term with a stiff system of ordinary differential equations (ODEs), the so-called membrane model, describing the ionic currents through the cellular membrane. A novel uncoupled solver for the Bidomain system is here introduced, based on solving twice the parabolic PDE and once the elliptic PDE at each time step, and it is compared with a usual coupled solver. Three-dimensional numerical tests have been performed in order to show that the proposed uncoupled method has the same accuracy of the coupled strategy. Parallel numerical tests on structured meshes have also shown that the uncoupled technique is as scalable as the coupled one. Moreover, the conjugate gradient method preconditioned by Multilevel Hybrid Schwarz preconditioners converges faster for the linear systems deriving from the uncoupled method than from the coupled one. Finally, in all parallel numerical tests considered, the uncoupled technique proposed is always about two or three times faster than the coupled approach.
LA - eng
KW - operator splitting; multilevel preconditioners; parallel computing
UR - http://eudml.org/doc/273115
ER -

References

top
  1. [1] T.M. Austin, M.L. Trew and A.J. Pullan, Solving the cardiac Bidomain equations for discontinuous conductivities. IEEE Trans. Biomed. Eng.53 (2006) 1265–1272. 
  2. [2] S. Balay, K. Buschelman, W.D. Gropp, D. Kaushik, M. Knepley, L. Curfman McInnes, B.F. Smith and H. Zhang, PETSc Users Manual.Tech. Rep. ANL-95/11 - Revision 2.1.5, Argonne National Laboratory (2002). 
  3. [3] S. Balay, K. Buschelman, W.D. Gropp, D. Kaushik, M. Knepley, L. Curfman McInnes, B.F. Smith and H. Zhang, PETSc home page. http://www.mcs.anl.gov/petsc (2001). 
  4. [4] M. Boulakia, S. Cazeau, M.A. Fernandez, J.-F. Gerbeau and N. Zemzemi, Mathematical modeling of electrocardiograms: a numerical study. Ann. Biomed. Eng.38 (2010) 1071–1097. 
  5. [5] R.H. Clayton, O. Bernus, E.M. Cherry, H. Dierckx, F.H. Fenton, L. Mirabella, A.V. Panfilov, F.B. Sachse, G. Seemann and H. Zhang, Models of cardiac tissue electrophysiology: Progress, challenges and open questions. Progr. Biophys. Molec. Biol.104 (2011) 22–48. 
  6. [6] P. Colli Franzone and L.F. Pavarino, A parallel solver for reaction-diffusion systems in computational electrocardiology. Math. Mod. Meth. Appl. Sci.14 (2004) 883–911. Zbl1068.92024MR2069498
  7. [7] P. Colli Franzone, L.F. Pavarino and S. Scacchi, Mathematical and numerical methods for reaction–diffusion models in electrocardiology, in Modeling of Physiological flows, edited by D. Ambrosi, A. Quarteroni and G. Rozza. Springer (2011) 107–142. Zbl1211.92008
  8. [8] P. Colli Franzone, L.F. Pavarino and B. Taccardi, Simulating patterns of excitation, repolarization and action potential duration with cardiac bidomain and monodomain models. Math. Biosci.197 (2005) 33–66. Zbl1074.92004MR2167484
  9. [9] P. Colli Franzone, P. Deuflhard, B. Erdmann, J. Lang and L.F. Pavarino, Adaptivity in space and time for reaction-diffusion systems in Electrocardiology. SIAM J. Sci. Comput.28 (2006) 942–962. Zbl1114.65110MR2240798
  10. [10] P. Deuflhard, B. Erdmann, R. Roitzsch and G.T. Lines, Adaptive finite element simulation of ventricular fibrillation dynamics. Comput. Visual. Sci.12 (2009) 201–205. MR2507225
  11. [11] M. Dryja, M.V. Sarkis and O.B. Widlund, Multilevel Schwarz methods for elliptic problems with discontinuous coefficients in three dimensions. Numer. Math.72 (1996) 313–348. Zbl0857.65131MR1367653
  12. [12] M. Dryja and O.B. Widlund, Multilevel additive methods for elliptic finite element problems. Parallel algorithms for partial differential equations (Kiel 1990) Notes Numer. Fluid Mech. 31 (1991) 58–69. Zbl0783.65057MR1167868
  13. [13] M. Dryja and O.B. Widlund, Domain decomposition algorithms with small overlap. SIAM J. Sci. Comput.15 (1994) 604–620. Zbl0802.65119MR1273155
  14. [14] M. Ethier and Y. Bourgault, Semi-implicit time-discretization schemes for the Bidomain model. SIAM J. Numer. Anal.46 (2008) 2443–2468. Zbl1182.92009MR2421042
  15. [15] M.A. Fernandez and N. Zemzemi, Decoupled time–marching schemes in computational cardiac electrophysiology and ECG numerical simulation. Math. Biosci.226 (2010) 58–75. Zbl1193.92024MR2675617
  16. [16] M. Fink, S.A. Niederer, E.M. Cherry, F.H. Fenton, J.T. Koivumaki, G. Seemann, T. Rudiger, H. Zhang, F.B. Sachse, D. Beard, E.J. Crampin and N.P. Smith, Cardiac cell modelling: observations from the heart of the cardiac physiome project. Prog. Biophys. Mol. Biol.104 (2011) 2–21. 
  17. [17] L.G. Giorda, L. Mirabella, F. Nobile, M. Perego and A. Veneziani, A model-based block-triangular preconditioner for the Bidomain system in electrocardiology. J. Comput. Phys.228 (2009) 3625–3639. Zbl1187.92053MR2511070
  18. [18] L. Gerardo Giorda, M. Perego and A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology. Math. Model. Numer. Anal.45 (2011) 309–334. Zbl1274.92022MR2804641
  19. [19] I.J. LeGrice, B.H. Smaill, L.Z. Chai, S.G. Edgar, J.B. Gavin and P.J. Hunter, Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. Amer. J. Physiol. Heart Circ. Physiol. 269 (1995) H571–H582. 
  20. [20] S. Linge, J. Sundnes, M. Hanslien, G.T. Lines and A. Tveito, Numerical solution of the bidomain equations. Philos. Trans. R. Soc. A367 (2009) 1931–1950. Zbl1185.65169MR2512073
  21. [21] C. Luo and Y. Rudy, A model of the ventricular cardiac action potential: depolarization, repolarization, and their interaction. Circ. Res.68 (1991) 1501–1526. 
  22. [22] K.-A. Mardal, B.F. Nielsen, X. Cai and A. Tveito, An order optimal solver for the discretized bidomain equations. Numer. Linear Algebra Appl.14 (2007) 83–98. Zbl1199.65111MR2292297
  23. [23] G. Karypis and V. Kumar, MeTis: Unstructured Graph Partitioning and Sparse Matrix Ordering System, Version 4.0. http://www.cs.umn.edu/~metis/. University of Minnesota, Minneapolis, MN (2009). 
  24. [24] M. Munteanu and L.F. Pavarino, Decoupled Schwarz algorithms for implicit discretization of nonlinear Monodomain and Bidomain systems. Math. Mod. Meth. Appl. Sci.19 (2009) 1065–1097. Zbl1178.65116MR2553178
  25. [25] M. Munteanu, L.F. Pavarino and S. Scacchi. A scalable Newton-Krylov-Schwarz method for the Bidomain reaction-diffusion system. SIAM J. Sci. Comput.31 (2009) 3861–3883. Zbl1205.65261MR2556566
  26. [26] M. Murillo and X.-C. Cai, A fully implicit parallel algorithm for simulating the non-linear electrical activity of the heart. Numer. Linear Algebra Appl.11 (2004) 261–277. Zbl1114.65112MR2065816
  27. [27] J.S. Neu and W. Krassowska, Homogenization of syncytial tissues. Crit. Rev. Biomed. Eng.21 (1993) 137–199. 
  28. [28] P. Pathmanathan, M.O. Bernabeu, R. Bordas, J. Cooper, A. Garny, J.M. Pitt-Francis, J.P. Whiteley and D.J. Gavaghan, A numerical guide to the solution of the bidomain equations of cardiac electrophysiology. Progr. Biophys. Molec. Biol.102 (2010) 136–155. 
  29. [29] L.F. Pavarino and S. Scacchi, Multilevel additive Schwarz preconditioners for the Bidomain reaction-diffusion system. SIAM J. Sci. Comput.31 (2008) 420–443. Zbl1185.65179MR2460784
  30. [30] L.F. Pavarino and S. Scacchi, Parallel Multilevel Schwarz and Block Preconditioners for the Bidomain Parabolic-Parabolic and Parabolic-Elliptic Formulations. SIAM J. Sci. Comput.33 (2011) 1897–1919. Zbl1233.65070MR2831039
  31. [31] M. Pennacchio, G. Savaré and P.C. Franzone. Multiscale modeling for the bioelectric activity of the heart. SIAM J. Math. Anal.37 (2006) 1333–1370. Zbl1113.35019MR2192297
  32. [32] M. Pennacchio and V. Simoncini, Efficient algebraic solution of reaction-diffusion systems for the cardiac excitation process. J. Comput. Appl. Math.145 (2002) 49–70. Zbl1006.65102MR1914350
  33. [33] M. Pennacchio and V. Simoncini, Algebraic multigrid preconditioners for the bidomain reaction-diffusion system. Appl. Numer. Math.59 (2009) 3033–3050. Zbl1171.92017MR2560833
  34. [34] M. Pennacchio and V. Simoncini, Fast structured AMG preconditioning for the bidomain model in electrocardiology. SIAM J. Sci. Comput.33 (2011) 721–745. Zbl1227.92034MR2785969
  35. [35] G. Plank, M. Liebmann, R. Weber dos Santos, E.J. Vigmond and G. Haase, Algebraic Multigrid Preconditioner for the Cardiac Bidomain Model. IEEE Trans. Biomed. Eng. 54 (2007) 585–596. 
  36. [36] M. Potse, B. Dubè, J. Richer, A. Vinet and R. Gulrajani, A comparison of Monodomain and Bidomain reaction–diffusion models for action potential propagation in the human heart. IEEE Trans. Biomed. Eng.53 (2006) 2425–2434. 
  37. [37] P.-A. Raviart, The use of numerical integration in finite element methods for solving parabolic equations. In Topics in Numerical Analysis, edited by J.J.H. Miller. Academic Press (1973) 233–264. Zbl0293.65086MR345428
  38. [38] Z. Qu and A. Garfinkel, An advanced algorithm for solving partial differential equation in cardiac conduction. IEEE Trans. Biomed. Eng.46 (1999) 1166–1168. 
  39. [39] A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations. Springer (1997). Zbl0803.65088MR1299729
  40. [40] S. Scacchi, A hybrid multilevel Schwarz method for the bidomain model. Comput. Methods Appl. Mech. Eng.197 (2008) 4051–4061. Zbl1194.78048MR2458128
  41. [41] S. Scacchi, A multilevel hybrid Newton-Krylov-Schwarz method for the Bidomain model of electrocardiology. Comput. Methods Appl. Mech. Eng.200 (2011) 717–725. Zbl1225.92011MR2749030
  42. [42] S. Scacchi, P. Colli Franzone, L.F. Pavarino and B. Taccardi, Computing cardiac recovery maps from electrograms and monophasic action potentials under heterogeneous and ischemic conditions. Math. Mod. Methods Appl. Sci.20 (2010) 1089–1127. Zbl1194.92050MR2673412
  43. [43] K.B. Skouibine, N. Trayanova and P. Moore, A numerically efficient model for the simulation of defibrillation in an active bidomain sheet of myocardium. Math. Biosci.166 (2000) 85–100. Zbl0963.92019
  44. [44] B.F. Smith, P. Bjørstad and W.D. Gropp, Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations, Cambridge University Press (1996). Zbl0857.65126MR1410757
  45. [45] J.A. Southern, G. Plank, E.J. Vigmond and J.P. Whiteley, Solving the coupled system improves computational efficiency of the Bidomain equations. IEEE Trans. Biomed. Eng.56 (2009) 2404–2412. 
  46. [46] J. Sundnes, G.T. Lines, K.A. Mardal and A. Tveito, Multigrid block preconditioning for a coupled system of partial differential equations modeling the electrical activity in the heart. Comput. Methods Biomech. Biomed. Eng.5 (2002) 397–409. 
  47. [47] J. Sundnes, G.T. Lines and A. Tveito, An operator splitting method for solving the bidomain equations coupled to a volume conductor model for the torso. Math. Biosci.194 (2005) 233–248. Zbl1063.92018MR2142490
  48. [48] H. Si, http://tetgen.berlios.de/. Weierstrass Institute for Applied Analysis and Stochastics, Berlin, Germany. 
  49. [49] V. Thomée, Galerkin Finite Element Methods for Parabolic Problems. Springer (1997). Zbl1105.65102MR1479170
  50. [50] A. Toselli and O.B. Widlund, Domain Decomposition Methods: Algorithms and Theory. Comput. Math. Springer-Verlag, Berlin 34 (2004). Zbl1069.65138
  51. [51] J.A. Trangenstein and C. Kim, Operator splitting and adaptive mesh refinement for the Luo-Rudy I model. J. Comput. Phys.196 (2004) 645–679. Zbl1056.92014MR2054352
  52. [52] E.J. Vigmond, F. Aguel and N.A. Trayanova, Computational techniques for solving the bidomain equations in three dimensions. IEEE Trans. Biomed. Eng.49 (2002) 1260–1269. 
  53. [53] E.J. Vigmond, R. Weber dos Santos, A.J. Prassl, M. Deo and G. Plank, Solvers for the cardiac bidomain equations. Progr. Biophys. Molec. Biol. 96 (2008) 3–18. 
  54. [54] R. Weber dos Santos, G. Plank, S. Bauer and E.J. Vigmond, Parallel multigrid preconditioner for the cardiac bidomain model. IEEE Trans. Biomed. Eng. 51 (2004) 1960–1968. 
  55. [55] J.P. Whiteley, An efficient numerical technique for the solution of the monodomain and bidomain equations. IEEE Trans. Biomed. Eng.53 (2006) 2139–2147. 
  56. [56] M. Zaniboni, 3D current-voltage-time surfaces unveil critical repolarization differences underlying similar cardiac action potentials: A model study. Math. Biosci.233 (2011) 98–110. Zbl1226.92013MR2856847
  57. [57] X. Zhang, Multilevel Schwarz methods. Numer. Math.63 (1992) 521–539. Zbl0796.65129MR1189535

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.