Stable time integration suppresses unphysical oscillations in the bidomain model

Frontiers in Physics

Spiteri, Raymond J.¹  Torabi Ziaratgahi, Saeed²  Sundnes, Joakim³  Marsh, Megan E.⁴ 

[1] Numerical Simulation Laboratory, Department of Computer Science, University of Saskatchewan, Saskatoon, SK, Canada;Numerical Simulation Laboratory, Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, SK, Canada;Simula Research Laboratory, Fornebu, Norway;Solido Design Automation, Saskatoon, SK, Canada

关键词: bidomain model;    operator splitting;    unphysical oscillations;    SDIRK method;    L-stability;   

DOI  :  10.3389/fphy.2014.00040

学科分类：物理（综合）

来源: Frontiers

【 摘 要 】

The bidomain model is a popular model for simulating electrical activity in cardiac tissue. It is a continuum-based model consisting of non-linear ordinary differential equations (ODEs) describing spatially averaged cellular reactions and a system of partial differential equations (PDEs) describing electrodiffusion on tissue level. Because of this multi-scale, ODE/PDE structure of the model, operator-splitting methods that treat the ODEs and PDEs in separate steps are natural candidates as numerical solution methods. Second-order methods can generally be expected to be more effective than first-order methods under normal accuracy requirements. However, the simplest and the most commonly applied second-order method for the PDE step, the Crank--Nicolson (CN) method, may generate unphysical oscillations. In this paper, we investigate the performance of a two-stage, L-stable singly diagonally implicit Runge--Kutta method for solving the PDEs of the bidomain model. Numerical experiments show that the enhanced stability property of this method leads to more physically realistic numerical simulations compared to both the CN and backward Euler methods.

【 授权许可】

CC BY