Analysis of time-stepping methods for the monodomain model
From MaRDI portal
Reaction-diffusion equations (35K57) PDEs in connection with biology, chemistry and other natural sciences (35Q92) Numerical methods for stiff equations (65L04) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Biomedical imaging and signal processing (92C55)
Abstract: To a large extent, the stiffness of the bidomain and monodomain models depends on the choice of the ionic model, which varies in terms of complexity and realism. In this paper, we compare and analyze a variety of time-stepping methods: explicit or semi-implicit, operator splitting, exponential, and deferred correction methods. We compare these methods for solving the bidomain model coupled with three ionic models of varying complexity and stiffness: the phenomenological Mitchell-Schaeffer model, the more realistic Beeler-Reuter model, and the stiff and very complex ten Tuscher-Noble-Noble-Panfilov (TNNP) model. For each method, we derive absolute stability criteria of the spatially discretized monodomain model and verify that the theoretical critical time-steps obtained closely match the ones in numerical experiments. We also verify that the numerical methods achieve an optimal order of convergence on the model variables and derived quantities (such as speed of the wave, depolarization time), and this in spite of the local non-differentiability of some of the ionic models. The efficiency of the different methods is also considered by comparing computational times for similar accuracy. Conclusions are drawn on the methods to be used to solve the monodomain model based on the model stiffness and complexity, measured respectively by the eigenvalues of the model's Jacobian and the number of variables, and based on strict stability and accuracy criteria.
Recommendations
- Analysis of a Monte Carlo boundary propagation method
- scientific article; zbMATH DE number 4068217
- Asymptotic behaviour of time stepping methods for phase field models
- Stability and convergence of time-stepping methods for a nonlocal model for diffusion
- scientific article; zbMATH DE number 2187914
- Accurate monotonicity-preserving schemes with Runge-Kutta time stepping
- Analysis of a variable time-step discretization for the Penrose-Fife phase relaxation problem
- Convergence of the time-discretized monotonic schemes
Cites work
- A numerical method for the solution of the bidomain equations in cardiac tissue
- A two-current model for the dynamics of cardiac membrane
- An efficient generalization of the Rush--Larsen method for solving electro-physiology membrane equations
- Computing the Electrical Activity in the Heart
- Deferred correction methods for initial boundary value problems
- High-order time stepping for the incompressible Navier-Stokes equations
- scientific article; zbMATH DE number 940566 (Why is no real title available?)
- scientific article; zbMATH DE number 6154734 (Why is no real title available?)
- Implicit-Explicit Methods for Time-Dependent Partial Differential Equations
- Mathematical physiology
- On the accuracy of operator splitting for the monodomain model of electrophysiology
- On the Construction and Comparison of Difference Schemes
- Operator splitting and adaptive mesh refinement for the Luo-Rudy I model
- Semi-Implicit Time-Discretization Schemes for the Bidomain Model
- Solving Ordinary Differential Equations I
- Taking math to heart: mathematical challenges in cardiac electrophysiology
Cited in
(11)- Moving-habitat models: a numerical approach
- A curvilinear isogeometric framework for the electromechanical activation of thin muscular tissues
- Gating-enhanced IMEX splitting methods for cardiac monodomain simulation
- On efficiency and accuracy in cardioelectric simulation
- Semi-Implicit Time-Discretization Schemes for the Bidomain Model
- High-order operator splitting for the bidomain and monodomain models
- scientific article; zbMATH DE number 7524472 (Why is no real title available?)
- An explicit local space-time adaptive framework for monodomain models in cardiac electrophysiology
- Arbitrary high-order unconditionally stable methods for reaction-diffusion equations with inhomogeneous boundary condition via deferred correction
- Explicit stabilized multirate methods for the monodomain model in cardiac electrophysiology
- A splitting, discontinuous Galerkin solver for the cell-by-cell electroneutral Nernst-Planck framework
This page was built for publication: Analysis of time-stepping methods for the monodomain model
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2204153)
