Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
129 tokens/sec
GPT-4o
28 tokens/sec
Gemini 2.5 Pro Pro
42 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

Explicit stabilized multirate methods for the monodomain model in cardiac electrophysiology (2401.01745v3)

Published 3 Jan 2024 in math.NA and cs.NA

Abstract: Fully explicit stabilized multirate (mRKC) methods are well-suited for the numerical solution of large multiscale systems of stiff ordinary differential equations thanks to their improved stability properties. To demonstrate their efficiency for the numerical solution of stiff, multiscale, nonlinear parabolic PDE's, we apply mRKC methods to the monodomain equation from cardiac electrophysiology. In doing so, we propose an improved version, specifically tailored to the monodomain model, which leads to the explicit exponential multirate stabilized (emRKC) method. Several numerical experiments are conducted to evaluate the efficiency of both mRKC and emRKC, while taking into account different finite element meshes (structured and unstructured) and realistic ionic models. The new emRKC method typically outperforms a standard implicit-explicit baseline method for cardiac electrophysiology. Code profiling and strong scalability results further demonstrate that emRKC is faster and inherently parallel without sacrificing accuracy.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (66)
  1. A. Abdulle. Fourth order Chebyshev methods with recurrence relation. SIAM Journal on Scientific Computing, 23(6):2041–2054, 2002.
  2. A. Abdulle. Explicit Stabilized Runge–Kutta Methods. In Encyclopedia of Applied and Computational Mathematics, pages 460–468. Springer, 2015.
  3. Explicit stabilized multirate method for stiff differential equations. Mathematics of Computation, 91(338):2681–2714, 2022.
  4. A. Abdulle and A. A. Medovikov. Second order Chebyshev methods based on orthogonal polynomials. Numerische Mathematik, 18:1–18, 2001.
  5. A. Abdulle and G. Vilmart. PIROCK: A swiss-knife partitioned implicit-explicit orthogonal Runge-Kutta Chebyshev integrator for stiff diffusion-advection-reaction problems with or without noise. Journal of Computational Physics, 242:869–888, 2013.
  6. A matrix–free high–order solver for the numerical solution of cardiac electrophysiology. Journal of Computational Physics, 478:111984, 2023.
  7. Lifex-fiber: An open tool for myofibers generation in cardiac computational models. BMC Bioinformatics, 24(1):143, Apr. 2023.
  8. Unified Form Language: A domain-specific language for weak formulations of partial differential equations. ACM Transactions on Mathematical Software, 40, 2014.
  9. PETSc Web page, 2023.
  10. Models of cardiac tissue electrophysiology: Progress, challenges and open questions. Progress in Biophysics and Molecular Biology, 104(1-3):22–48, 2011.
  11. Myokit: A simple interface to cardiac cellular electrophysiology. Progress in Biophysics and Molecular Biology, 120(1-3):100–114, Jan. 2016.
  12. Higher order triangular finite elements with mass lumping for the wave equation. SIAM J. Numer. Analysis, 38:2047–2078, 2001.
  13. Mathematical cardiac electrophysiology, volume 13. Springer, 2014.
  14. Rush-Larsen time-stepping methods of high order for stiff problems in cardiac electrophysiology. Electronic Transactions on Numerical Analysis, 52:342–357, 2020.
  15. Ionic mechanisms underlying human atrial action potential properties: Insights from a mathematical model. American Journal of Physiology-Heart and Circulatory Physiology, 275(1):H301–H321, 1998.
  16. Mixed-precision explicit stabilized Runge-Kutta methods for single- and multi-scale differential equations. Journal of Computational Physics, 464, 2022.
  17. Parallel distributed computing using Python. Advances in Water Resources, 34(9):1124–1139, 2011.
  18. Simulation of human ischemic stroke in realistic 3D geometry. Communications in Nonlinear Science and Numerical Simulation, 18(6):1539–1557, 2013.
  19. Hypre: A library of high performance preconditioners. In International Conference on Computational Science, pages 632–641. Springer, 2002.
  20. Personalized computational modeling of left atrial geometry and transmural myofiber architecture. Medical Image Analysis, 47:180–190, July 2018.
  21. Personalized computational finite element models for left atrial electromechanics, Sept. 2021.
  22. Impact of mechanical deformation on pseudo-ECG: a simulation study. EP Europace, 18(suppl_4):iv77–iv84, 2016.
  23. Fast Characterization of Inducible Regions of Atrial Fibrillation Models With Multi-Fidelity Gaussian Process Classification. Frontiers in Physiology, 13:757159, Mar. 2022.
  24. New higher-order mass-lumped tetrahedral elements for wave propagation modelling. SIAM J. Sci. Comput., 40:A2830–A2857, 2018.
  25. Epicardial Fibrosis Explains Increased Endo–Epicardial Dissociation and Epicardial Breakthroughs in Human Atrial Fibrillation. Frontiers in Physiology, 11:68, 2020.
  26. Left atrial appendage electrical isolation reduces atrial fibrillation recurrences: simulation study. Circulation: Arrhythmia and Electrophysiology, 14(1):e009230, 2021.
  27. M. J. Grote and J. Diaz. Energy Conserving Explicit Local Time Stepping for Second-Order Wave Equations. SIAM Journal on Scientific Computing, 31(3):1985–2014, 2009.
  28. Runge–Kutta-based explicit local time-stepping methods for wave propagation. SIAM Journal on Scientific Computing, 37(2):A747–A775, 2015.
  29. M. J. Grote and T. Mitkova. Explicit local time-stepping methods for Maxwell’s equations. Journal of Computational and Applied Mathematics, 234(12):3283–3302, Oct. 2010.
  30. M. J. Grote and T. Mitkova. High-order explicit local time-stepping methods for damped wave equations. J. Comput. Appl. Math., 239:270–289, 2013.
  31. Solving Ordinary Differential Equations I, volume 8. Springer-Verlag, Berlin, 2 edition, 2008.
  32. E. Hairer and G. Wanner. Solving Ordinary Differential Equations II, volume 14. Springer-Verlag, Berlin, 2002.
  33. M. Hochbruck and A. Ostermann. Exponential integrators. Acta Numerica, 19:209–286, 2010.
  34. A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology, 117(4):500–544, Aug. 1952.
  35. Convergence analysis of BDDC preconditioners for composite DG discretizations of the cardiac cell-by-cell model. SIAM Journal on Scientific Computing, 45(6):A2836–A2857, 2023.
  36. Hybrid Parallelization of a Large-Scale Heart Model. In Facing the Multicore-Challenge II, pages 120–132. Springer, 2012.
  37. B. Lindberg. IMPEX: A program package for solution of systems of stiff differential equations. Technical report, Dept. of Information Processing, Royal Inst. of Tech., Stockholm, 1972.
  38. B. Lindberg. IMPEX2: A procedure for solution of systems of stiff differential equations. Technical Report TRITA-NA-7303, Department of Information Processing, Royal Institute of Technology, Stockholm, Sweden, 1973.
  39. Efficient time splitting schemes for the monodomain equation in cardiac electrophysiology. International Journal for Numerical Methods in Biomedical Engineering, 39(2), 2023.
  40. The CellML model repository. Bioinformatics (Oxford, England), 24(18):2122–2123, 2008.
  41. A stabilized Runge-Kutta-Legendre method for explicit super-time-stepping of parabolic and mixed equations. Journal of Computational Physics, 257(PA):594–626, 2014.
  42. Verification of cardiac tissue electrophysiology simulators using an N-version benchmark. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 369:4331–4351, 2011.
  43. An explicit local space-time adaptive framework for monodomain models in cardiac electrophysiology. Computer Methods in Applied Mechanics and Engineering, 422:116806, 2024.
  44. The significant effect of the choice of ionic current integration method in cardiac electro-physiological simulations. International Journal for Numerical Methods in Biomedical Engineering, 27(11):1751–1770, Nov. 2011.
  45. Space-discretization error analysis and stabilization schemes for conduction velocity in cardiac electrophysiology. International Journal for Numerical Methods in Biomedical Engineering, 32(10):e02762, 2016.
  46. On Sampling Spatially-Correlated Random Fields for Complex Geometries. In Y. Coudière, V. Ozenne, E. Vigmond, and N. Zemzemi, editors, Functional Imaging and Modeling of the Heart, volume 11504, pages 103–111, Cham, 2019. Springer International Publishing.
  47. Modeling cardiac muscle fibers in ventricular and atrial electrophysiology simulations. Computer Methods in Applied Mechanics and Engineering, 373:113468, 2021.
  48. Modeling the cardiac electromechanical function: A mathematical journey. Bulletin of the American Mathematical Society, 59(3):371–403, 2022.
  49. G. Rosilho De Souza. Numerical Methods for Deterministic and Stochastic Differential Equations with Multiple Scales and High Contrasts. PhD thesis, EPFL, 2020.
  50. G. Rosilho De Souza. mRKC: A multirate Runge–Kutta–Chebyshev code, 2022.
  51. Boundary integral formulation of the cell-by-cell model of cardiac electrophysiology. Engineering Analysis with Boundary Elements, 158:239–251, 2024.
  52. Effect of gap junction distribution, size, and shape on the conduction velocity in a cell-by-cell model for electrophysiology. In O. Bernard, P. Clarysse, N. Duchateau, J. Ohayon, and M. Viallon, editors, Functional Imaging and Modeling of the Heart, volume 13958 of Lecture Notes in Computer Science, pages 117–126, Cham, 2023. Springer.
  53. S. Rush and H. Larsen. A Practical Algorithm for Solving Dynamic Membrane Equations. IEEE Transactions on Biomedical Engineering, BME-25(4):389–392, 1978.
  54. The fibrotic kernel signature: Simulation-free prediction of atrial fibrillation. In O. Bernard, P. Clarysse, N. Duchateau, J. Ohayon, and M. Viallon, editors, Functional Imaging and Modeling of the Heart, volume 13958 of Lecture Notes in Computer Science, pages 87–96, Cham, 2023. Springer.
  55. Basix: A runtime finite element basis evaluation library. Journal of Open Source Software, 7(73):3982, 2022.
  56. Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes. ACM Transactions on Mathematical Software, 48(2):18:1–18:23, 2022.
  57. L. F. Shampine. Lipschitz constants and robust ode codes’, Computational Methods in Nonlinear Mechanics. In Proceedings of the TICOM Second International Conference, JT Oden, Ed., North-Holland Publishing Company, New York, pages 427–449, 1980.
  58. L. F. Shampine. Diagnosing Stiffness for Runge–Kutta Methods. SIAM Journal on Scientific and Statistical Computing, 12(2):260–272, Mar. 1991.
  59. RKC: An explicit solver for parabolic PDEs. Journal of Computational and Applied Mathematics, 88(2):315–326, 1998.
  60. Stiffness analysis of cardiac electrophysiological models. Annals of Biomedical Engineering, 38(12):3592–3604, 2010.
  61. Erratum: Stiffness analysis of cardiac electrophysiological models. Annals of Biomedical Engineering, 40(7):1622–1625, 2012.
  62. K. H. W. J. ten Tusscher and A. V. Panfilov. Alternans and spiral breakup in a human ventricular tissue model. American Journal of Physiology. Heart and Circulatory Physiology, 291(3):H1088–1100, Sept. 2006.
  63. P. J. Van der Houwen and B. P. Sommeijer. On the internal stability of explicit, $m$-stage Runge–Kutta methods for large m𝑚mitalic_m-values. Zeitschrift für Angewandte Mathematik und Mechanik, 60(10):479–485, 1980.
  64. J. G. Verwer. An implementation of a class of stabilized explicit methods for the time integration of parabolic equations. ACM Transactions on Mathematical Software (TOMS), 6(2):188–205, 1980.
  65. J. G. Verwer. Explicit Runge–Kutta methods for parabolic partial differential equations. Applied Numerical Mathematics, 22(1-3):359–379, 1996.
  66. Convergence properties of the Runge–Kutta–Chebyshev method. Numerische Mathematik, 57(1):157–178, 1990.
Citations (1)
View on Semantic Scholar

Summary

We haven't generated a summary for this paper yet.

Ai Sparkles Streamline Icon: https://streamlinehq.com Summarize Now