Key result
Adaptive fourth-order compact schemes allow ~45% coarser meshes than second-order methods in cardiac electrophysiology models.
Why the study?
Anisotropic reaction–diffusion equations with rapidly varying reaction terms require fine meshes and small time steps to avoid spurious oscillations, motivating development of more efficient numerical schemes.
A novel fourth-order compact scheme with adaptive time steps allows for 45% coarser meshes without compromising accuracy in simulating anisotropic cardiac electrophysiology.
May enable efficient cardiac electrophysiology simulations; leaves open clinical translation and validation.
Many problems in biology and engineering are governed by anisotropic reaction–diffusion equations with a very rapidly varying reaction term. These characteristics of the system imply the use of very fine meshes and small time steps in order to accurately capture the propagating wave avoiding the appearance of spurious oscillations in the wave front. This work develops a fourth‐order compact scheme for anisotropic reaction–diffusion equations with stiff reactive terms. As mentioned, the scheme accounts for the anisotropy of the media and incorporates an adaptive time step for handling the stiff reactive term. The high‐order scheme allows working with coarser meshes without compromising numerical accuracy rendering a more efficient numerical algorithm by reducing the total computation time and memory requirements. The order of convergence of the method has been demonstrated on an analytical solution with Neumann boundary conditions. The scheme has also been implemented for the solution of anisotropic electrophysiology problems. Anisotropic square samples of normal and ischemic cardiac tissue have been simulated by means of the monodomain model with the reactive term defined by Luo–Rudy II dynamics. The simulations proved the effectiveness of the method in handling anisotropic heterogeneous non‐linear reaction–diffusion problems. Bidimensional tests also indicate that the fourth‐order scheme requires meshes about 45% coarser than the standard second‐order method in order to achieve the same accuracy of the results.
No takes yet. Share an insight, caveat, or question.
Heidenreich et al. (2009) studied Anisotropic electrophysiology problems (normal and ischemic cardiac tissue). Fourth-order compact scheme with adaptive time step vs. Standard second-order method was evaluated on Mesh coarseness required to achieve the same accuracy. A fourth-order compact scheme with adaptive time step required meshes about 45% coarser than the standard second-order method to achieve the same accuracy in cardiac electrophysiology models.
Synapse has enriched 5 closely related papers on similar clinical questions. Consider them for comparative context: