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...

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Medical & biological engineering & computing Ročník 57; číslo 6; s. 1367 - 1379
Hlavní autori: Chen, Xing-Ji, Luo, Ching-Hsing, Chen, Min-Hung, Zhou, Xiang
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Berlin/Heidelberg Springer Berlin Heidelberg 01.06.2019
Springer Nature B.V
Predmet:
ISSN:0140-0118, 1741-0444, 1741-0444
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí: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.
Bibliografia: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