Abstract:The locality of solution features in cardiac electrophysiology simulations calls for adaptive methods. Due to the overhead incurred by established mesh refinement and coarsening, however, such approaches failed in accelerating the computations. Here we investigate a different route to spatial adaptivity that is based on nested subset selection for algebraic degrees of freedom in spectral deferred correction methods. This combination of algebraic adaptivity and iterative solvers for higher order collocation time stepping realizes a multirate integration with minimal overhead. This leads to moderate but significant speedups in both monodomain and cell-by-cell models of cardiac excitation, as demonstrated at four numerical examples.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to implement an efficient adaptive algorithm through the use of the Spectral Deferred Correction (SDC) method in cardiac electrophysiology simulations, so as to reduce the computational cost and improve the simulation efficiency. ### Specific problem description: 1. **Challenges in cardiac electrophysiology simulations**: - In cardiac electrophysiology simulations, the characteristics of the solution are local, such as the thin depolarization and repolarization wave fronts of cardiomyocytes. - Traditional time - stepping schemes and spatial discretization methods are inefficient in dealing with these local features, resulting in high computational overhead. 2. **Limitations of existing methods**: - Although the classical space - time grid adaptive method can reduce the number of degrees of freedom, due to frequent grid refinement and coarsening, error estimation, and reassembly of the mass matrix and stiffness matrix, the overall computational efficiency has not been significantly improved. - The Method of Lines with fixed spatial discretization is simple, but it wastes too much fine discretization in most of the domain. 3. **The proposed new method**: - The paper proposes a space - time adaptive method based on algebraic degree - of - freedom selection, combined with the SDC method of high - order time integration. - This method gradually restricts the computational region during the iteration process through nested sub - domain selection, thereby reducing the computational amount of later iterations. - This method achieves multi - rate integration, that is, different convergence orders are used in different regions instead of different time steps. 4. **Objectives**: - Through the above methods, to achieve a moderate but significant speed increase in the simulation of cardiac excitation propagation in single - domain models and cell - cell models. - Verify the effectiveness and efficiency of this method through four numerical examples. ### Summary: The paper aims to solve the problem of low computational efficiency in cardiac electrophysiology simulations by introducing a new space - time adaptive method combined with the SDC method. This method can significantly reduce the computational cost while maintaining accuracy, and is especially suitable for the simulation of cardiac excitation propagation.

Efficient adaptivity for simulating cardiac electrophysiology with spectral deferred correction methods

Parallel Solution In Simulation Of Cardiac Excitation Anisotropic Propagation

Application of Efficient Numerical Methods in Solution of Ordinary Differential Equations for Modeling Electrical Activity in Cardiac Cells

High-order parallel-in-time method for the monodomain equation in cardiac electrophysiology

Simulation of cardiac electrophysiology on next-generation high-performance computers

An explicit local space-time adaptive framework for monodomain models in cardiac electrophysiology

An Explicit Local Space-Time Adaptive Framework for Monodomain Models

Efficient solution of bidomain equations in simulation of cardiac excitation anisotropic propagation

Simulation of Active Cardiac Electromechanical Dynamics

A staggered-in-time and non-conforming-in-space numerical framework for realistic cardiac electrophysiology outputs

Efficient Numerical Methods for Solving the Membrane Equations of Human Ventricular Cell Model

Fast and Reliable Reduced-Order Models for Cardiac Electrophysiology

Integration of Different Cardiac Electrophysiological Models into a Single Simulation Pipeline

Space--time Isogeometric Analysis of cardiac electrophysiology

A non-conforming-in-space numerical framework for realistic cardiac electrophysiological outputs

Anisotropic space-time adaptation for reaction-diffusion problems

Efficient fully implicit time integration methods for modeling cardiac dynamics

A highly parallel algorithm for simulating the elastodynamics of a patient-specific human heart with four chambers using a heterogeneous hyperelastic model

A Cyclical Fast Iterative Method for Simulating Reentries in Cardiac Electrophysiology Using an Eikonal-Based Model

A fast computational model for the electrophysiology of the whole human heart

Adaptive time step selection for spectral deferred correction