Optimized Semi-implicit Methods for Modeling Cardiac Propagation

Authors

  • Riasat Khan Department of Electrical and Computer Engineering North South University, Dhaka, Bangladesh
  • Kwong T. Ng Department of Electrical and Computer Engineering New Mexico State University, Las Cruces, NM 88003, USA

Keywords:

Derivative-free optimization method, monodomain model, operator split method, pattern search algorithm, semi-implicit scheme

Abstract

Computer simulation of cardiac electrophysiology is now considered a powerful tool for exploring the causes of cardiac arrhythmias. Cardiac electric propagation has been studied using the monodomain model to describe wave propagation of action potential in the heart. The governing nonlinear reaction-diffusion partial differential equation is solved with the semi-implicit (implicit-explicit) method that does not have the stability limit of the explicit timestepping scheme. Both first order and second order semi-implicit techniques for temporal discretization are considered in this paper. Second order finite difference technique is used to discretize the spatial derivatives. An explicit finite difference scheme with 512×512 nodes and 0.1 μs time step is used as the benchmark for error calculation. APPSPACK, a parallel pattern search optimization software, is used to obtain the optimal semi-implicit parameters that give the lowest root-meansquare error. Results are presented for the semi-implicit techniques with or without the operator split or protective zone method. They demonstrate that the optimized second order semi-implicit method gives the best overall performance.

Downloads

Download data is not yet available.

References

C. S. Henriquez and W. Ying, The Bidomain Model of Cardiac Tissue: From Microscale to Macroscale. Springer, Boston, 2001.

S. L. Cloherty, N. H. Lovell, S. Dokos, and B. G. Celler, “A 2D monodomain model of rabbit sinoatrial node,” 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Istanbul, pp. 44-47, 2001.

M. Potse, B. Dube, J. Richer, A. Vinet, and R. M. Gulrajani, “A comparison of monodomain and bidomain reaction-diffusion models for action potential propagation in the human heart,” IEEE Tran. on Biomed. Eng., vol. 53, pp. 2425-2435, 2006.

A. R. Conn, K. Scheinberg, and L. N. Vicente, Introduction to Derivative-Free Optimization. SIAM, Philadelphia, 2009.

C. H. Luo and Y. Rudy, “A model of the ventricular cardiac action potential, depolarization, repolarizetion, and their interaction,” Circulation Research, vol. 68, pp. 1501-1526, 1991.

U. M. Ascher, S. J. Ruuth, and B. T. R. Wetton, “Implicit-Explicit methods for time-dependent partial differential equations,” SIAM J. Numeric. Anal., vol. 32, pp. 797-823, 1995.

A. Onose and B. Dumitrescu, “Adaptive cyclic and randomized coordinate descent for the sparse total least squares problem,” 23rd European Signal Processing Conference (EUSIPCO), Nice, pp. 1696-1700, 2015.

G. A. Gray and T. G. Kolda, “Algorithm 856: APPSPACK 4.0: Asynchronous parallel pattern search for derivative-free optimization,” ACM Trans. Math. Softw., vol. 32, pp. 485-507, 2006.

S. Krishnamoorthi, M. Sarkar, and W. S. Klug, “Numerical quadrature and operator splitting in finite element methods for cardiac electrophysiology,” International Journal for Numerical Methods in Biomedical Engineering, vol. 29, pp. 1243-1266, 2013.

R. Khan and K. T. Ng, “Higher order finite difference modeling of cardiac propagation,” IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Kansas City, MO, pp. 1945- 1951, 2017.

Downloads

Published

2019-10-01

How to Cite

[1]
Riasat Khan and Kwong T. Ng, “Optimized Semi-implicit Methods for Modeling Cardiac Propagation”, ACES Journal, vol. 34, no. 10, pp. 1616–1319, Oct. 2019.

Issue

Section

Articles