Stability Analysis of Decoupled Time-stepping Schemes for the Specialized Conduction System/myocardium Coupled Problem in Cardiology

¹ Université Tunis El Manar, FST 2092 Tunis Tunisia
² INRIA Bordeaux Sud-Ouest, Carmen project, 200 rue de la vieille tour 33405 Talence Cedex, France
³ IHU-Liryc, Electrophysiology and Heart Modeling Institute, foundation Bordeaux Université, F-33600 Pessac-Bordeaux, France

^* Corresponding author. E-mail: nejib.zemzemi@inria.fr

Abstract

The Purkinje network is the rapid conduction system in the heart. It ensures the physiological spread of the electrical wave in the ventricles. In this work, we consider a problem that models the coupling between the Purkinje network and the myocardium. We first prove the stability of the space semi-discretized problem. Then we present four different strategies for solving the Purkinje/ myocardium coupling. The strategies are based on different time discretization of the coupling terms. The first scheme is fully coupled, where the coupling terms are considered implicit. The second and the third schemes are based on Gauss-Seidel time-splitting schemes where one coupling term is considered explicit and the other is implicit. The last is a Jacobi-like time-splitting scheme where both coupling terms are considered explicit. Our main result is the proof of the stability of the three considered schemes under the same restriction on the time step. Moreover, we show that the energy of the problem is slightly affected by the time-splitting schemes. We illustrate the theoretical result by different numerical simulations in 2D. We also conduct 3D simulations using physiologically detailed ionic models.