A dual adaptive explicit time integration algorithm for efficiently solving the cardiac monodomain equation. (12th April 2021)

• Record Type: Journal Article
• Authors: Mountris, Konstantinos A.; Pueyo, Esther
• Is Part Of: International journal for numerical methods in biomedical engineering. Volume 37:Number 7(2021)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2021-04-12
• Subjects: adaptive explicit integration -- cardiac electrophysiology -- operator splitting; Biomedical engineering -- Periodicals; Imaging systems in medicine -- Periodicals; Numerical analysis -- Periodicals; Engineering mathematics -- Periodicals; 610.28
• DOI: 10.1002/cnm.3461
• Languages: English
• ISSNs: 2040-7939

Abstract

Abstract: The monodomain model is widely used in in‐silico cardiology to describe excitation propagation in the myocardium. Frequently, operator splitting is used to decouple the stiff reaction term and the diffusion term in the monodomain model so that they can be solved separately. Commonly, the diffusion term is solved implicitly with a large time step while the reaction term is solved by using an explicit method with adaptive time stepping. In this work, we propose a fully explicit method for the solution of the decoupled monodomain model. In contrast to semi‐implicit methods, fully explicit methods present lower memory footprint and higher scalability. However, such methods are only conditionally stable. We overcome the conditional stability limitation by proposing a dual adaptive explicit method in which adaptive time integration is applied for the solution of both the reaction and diffusion terms. We perform a set of numerical examples where cardiac propagation is simulated under physiological and pathophysiological conditions. Results show that the proposed method presents preserved accuracy and improved computational efficiency as compared to standard operator splitting‐based methods. Abstract : We propose a fully explicit method for solving the monodomain model for in‐silico cardiology by applying adaptive time integration for both the reaction and diffusion terms. Benchmark simulations show the ability to preserve accuracy and improve computational efficiency compared to standard operator‐splitting.