Related papers: Rush-Larsen time-stepping methods of high order fo…
The modern Markov chain models of ionic channels in excitable membranes are numerically stiff. The popular numerical methods for these models require very small time steps to ensure stability. Our objective is to formulate and test two…
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…
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…
The work deals with two major topics concerning the numerical analysis of Runge-Kutta-like (RK-like) methods, namely their stability and order of convergence. RK-like methods differ from additive RK methods in that their coefficients are…
A wide range of implicit time integration methods, including multi-step, implicit Runge-Kutta, and Galerkin finite-time element schemes, is evaluated in the context of chaotic dynamical systems. The schemes are applied to solve the Lorenz…
A novel second order family of explicit stabilized Runge-Kutta-Chebyshev methods for advection-diffusion-reaction equations is introduced. The new methods outperform existing schemes for relatively high Peclet number due to their favorable…
We consider quadrature formulas of high order in time based on Radau-type, L-stable implicit Runge-Kutta schemes to solve time dependent stiff PDEs. Instead of solving a large nonlinear system of equations, we develop a method that performs…
Unconditionally stable implicit time-marching methods are powerful in solving stiff differential equations efficiently. In this work, a novel framework to handle stiff physical terms implicitly is proposed. Both physical and numerical…
Models in cardiac electrophysiology are coupled systems of reaction diffusion PDE and of ODE. The ODE system displays a very stiff behavior. It is non linear and its upgrade at each time step is a preponderant load in the computational…
A stability analysis is performed on high-order schemes formulated using the Flux Reconstruction (FR) approach. The one-dimensional advection model equation is used for the assessment of the stability region of these schemes when coupled…
The explicit two-stage fourth-order (TSFO) temporal-spatial coupling method is efficient and compact but suffers severe time-step restrictions for stiff problems with multiple scales. To address Professor Jiequan Li's call for an implicit…
High-order adaptive time-stepping algorithms are of significant practical value and theoretical interest for accelerating long-time fluid-flow simulations and resolving complex dynamical behaviors. While several high-order implicit-explicit…
In this paper, Runge-Kutta-Gegenbauer (RKG) stability polynomials of arbitrarily high order of accuracy are introduced in closed form. The stability domain of RKG polynomials extends in the the real direction with the square of polynomial…
A combination of a steady-state preserving operator splitting method and a semi-implicit integration scheme is proposed for efficient time stepping in simulations of unsteady reacting flows, such as turbulent flames, using detailed chemical…
This paper develops a two-level fourth-order scheme for solving time-fractional convection-diffusion-reaction equation with variable coefficients subjected to suitable initial and boundary conditions. The basis properties of the new…
In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta…
Many time-dependent partial differential equations (PDEs) can be transformed into an ordinary differential equations (ODEs) containing moderately stiff and non-stiff terms after spatial semi-discretization. In the present paper, we…
Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the system, according to…
Most high order computational fluid dynamics (CFD) methods for compressible flows are based on Riemann solver for the flux evaluation and Runge-Kutta (RK) time stepping technique for temporal accuracy. The main advantage of this kind of…
In this paper, we investigate the stability and time-step constraints for solving advection-diffusion equations using exponential time differencing (ETD) Runge-Kutta (RK) methods in time and discontinuous Galerkin (DG) methods in space. We…