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ženo v:
Podrobná bibliografie
Vydáno v:Medical & biological engineering & computing Ročník 57; číslo 6; s. 1367 - 1379
Hlavní autoři: Chen, Xing-Ji, Luo, Ching-Hsing, Chen, Min-Hung, Zhou, Xiang
Médium: Journal Article
Jazyk:angličtina
Vydáno: Berlin/Heidelberg Springer Berlin Heidelberg 01.06.2019
Springer Nature B.V
Témata:
ISSN:0140-0118, 1741-0444, 1741-0444
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
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.
Bibliografie: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