Exponential Adams Bashforth ODE solver for stiff problems

Yves Coudière (1) , Charlie Douanla Lontsi (1) , Charles Pierre (2)

1 CARMEN - Modélisation et calculs pour l'électrophysiologie cardiaque
2 LMAP - Laboratoire de Mathématiques et de leurs Applications [Pau]

We analyze in this paper a class of explicit exponential methods for the time integration of stiff differential problems. Precisely, we considered Adams exponential integrators with general varying stabilizers. General stabilization brings flexibility and computational facilities for ODE systems and for semilinear evolution PDEs coupled with ODE systems. Stability and convergence are proven, by introducing a new framework that extends multistep linear methods. Dahlquist stability is numerically investigated. A(α)-stability is observed under a condition on the stabilizer, which is a singular property for explicit schemes. The methods are numerically studied for two stiff models in electrophysiology. Its performances are compared with several classical methods. We conclude that for stiff ODE systems, it provides a cheaper way to compute accurate solutions at large time steps than implicit solvers.

Mots clés: stiff equations; explicit high-order multistep methods; exponential integrators of Adams type; stability and convergence; Dahlquist stability .

Domaines: Analyse numérique [math.NA]

Yves Coudière, Charlie Douanla Lontsi, Charles Pierre. Exponential Adams Bashforth ODE solver for stiff problems. 2016. ⟨hal-01394036v1⟩