Related papers: Explicit stabilized multirate methods for the mono…
Explicit stabilized methods are highly efficient time integrators for large and stiff systems of ordinary differential equations especially when applied to semi-discrete parabolic problems. However, when local spatial mesh refinement is…
Stabilized Runge-Kutta methods are especially efficient for the numerical solution of large systems of stiff nonlinear differential equations because they are fully explicit. For semi-discrete parabolic problems, for instance, stabilized…
Stabilized explicit methods are particularly efficient for large systems of stiff stochastic differential equations (SDEs) due to their extended stability domain. However, they loose their efficiency when a severe stiffness is induced by…
This paper considers master equations for Markovian kinetic schemes that possess the detailed balance property. Chemical kinetics, as a prime example, often yields large-scale, highly stiff equations. Based on chemical intuitions, Sumiya et…
Singularly perturbed systems (SPSs) are prevalent in engineering applications, where numerically solving their initial value problems (IVPs) is challenging due to stiffness arising from multiple time scales. Classical explicit methods…
This work considers multirate generalized-structure additively partitioned Runge-Kutta (MrGARK) methods for solving stiff systems of ordinary differential equations (ODEs) with multiple time scales. These methods treat different partitions…
We present the meshfree Mixed Collocation Method (MCM) to solve the monodomain model for numerical simulation of cardiac electrophysiology. We apply MCM to simulate cardiac electrical propagation in 2D tissue sheets and 3D tissue slabs as…
The monodomain model is widely used in in-silico cardiology to describe excitation propagation in the myocardium. Frequently, operator splitting is used to decouple the stiff reaction term and the diffusion term in the monodomain model so…
In [1] is proposed a simplified DeC method, that, when combined with the residual distribution (RD) framework, allows to construct a high order, explicit FE scheme with continuous approximation avoiding the inversion of the mass matrix for…
In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong…
In this paper we propose a novel class of methods for high order accurate integration of multirate systems of ordinary differential equation initial-value problems. The proposed methods construct multirate schemes by approximating the…
This work proposes a stabilized space--time method for the monodomain equation coupled with the Rogers--McCulloch ionic model, which is widely used to simulate electrophysiological wave propagation in the cardiac tissue. By extending the…
This work aims at presenting a Discontinuous Galerkin (DG) formulation employing a spectral basis for two important models employed in cardiac electrophysiology, namely the monodomain and bidomain models. The use of DG methods is motivated…
We present novel entropy-conservative and entropy-stable multirate Runge-Kutta methods based on Paired Explicit Runge-Kutta (P-ERK) schemes with relaxation for conservation laws and related systems of partial differential equations.…
We present a novel approach which aims at high-performance uncertainty quantification for cardiac electrophysiology simulations. Employing the monodomain equation to model the transmembrane potential inside the cardiac cells, we evaluate…
Implicit-Explicit (IMEX) methods are flexible numerical time integration methods which solve an initial-value problem (IVP) that is partitioned into stiff and nonstiff processes with the goal of lower computational costs than a purely…
To address the issues of stability and accuracy for reaction-diffusion equations, the development of high order and stable time-stepping methods is necessary. This is particularly true in the context of cardiac electrophysiology, where…
The numerical tools to simulate the bidomain model in cardiac electrophysiology are constantly developing due to the great clinical interest and scientific advances in mathematical models and computational power. The bidomain model consists…
In this work we analyze the resort to high order exponential solvers for stiff ODEs in the context of cardiac electrophysiology modeling. The exponential Adams-Bashforth and the Rush-Larsen schemes will be considered up to order 4. These…
We propose entropy-preserving and entropy-stable partitioned Runge--Kutta (RK) methods. In particular, we extend the explicit relaxation Runge--Kutta methods to IMEX--RK methods and a class of explicit second-order multirate methods for…