Combination of “quadratic adaptive algorithm” and “hybrid operator splitting” or uniformization algorithms for stability against acceleration in the Markov model of sodium ion channels in the ventricular cell model

The Markovian model has generally been used for cardiac electrophysiological simulations. However, the Markovian model is so stiff that speeding up the computation of the algorithms with variable time-steps always results in simulation instability. In particular, the unstable simulations always occu...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Medical & biological engineering & computing Jg. 57; H. 6; S. 1367 - 1379
Hauptverfasser: Chen, Xing-Ji, Luo, Ching-Hsing, Chen, Min-Hung, Zhou, Xiang
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Berlin/Heidelberg Springer Berlin Heidelberg 01.06.2019
Springer Nature B.V
Schlagworte:
ISSN:0140-0118, 1741-0444, 1741-0444
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:The Markovian model has generally been used for cardiac electrophysiological simulations. However, the Markovian model is so stiff that speeding up the computation of the algorithms with variable time-steps always results in simulation instability. In particular, the unstable simulations always occur at a low voltage rate or current change, while transition rates in the Markovian model are changing markedly. The uniformization (UNI) method allows for a Markovian model simulation with high stability but also a high computation cost. To save computation costs with variable time-steps, we propose a speed increasing idea that is a compromise to the trade-off between stability and acceleration by combining Chen-Chen-Luo’s “quadratic adaptive algorithm” (CCL) method with “hybrid operator splitting” (HOS) into the solver (CCL + HOS solver). The computation cost of this CCL + HOS solver is approximately 24 times lower than the CCL + UNI solver, and the CCL + HOS solver can function 295 times faster in comparison to the HOS solver with a fixed time-step (DT). The suggested optimal solver should be CCL + HOS solver with a maximum time-step at 0.1 ms due to its high speed with low error. Additionally, the CCL method has much better performance and stability than the hybrid method in this single-cell model simulation.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
content type line 23
ISSN:0140-0118
1741-0444
1741-0444
DOI:10.1007/s11517-019-01956-5