PIROCK: a swiss-knife partitioned implicit-explicit orthogonal Runge-Kutta Chebyshev integrator for stiff diffusion-advection-reaction problems with or without noise

Assyr Abdulle (1) , Gilles Vilmart (2, 3)

1 EPFL - Ecole Polytechnique Fédérale de Lausanne
2 IPSO - Invariant Preserving SOlvers
3 IRMAR - Institut de Recherche Mathématique de Rennes

A partitioned implicit-explicit orthogonal Runge-Kutta method (PIROCK) is proposed for the time integration of diffusion-advection-reaction problems with possibly severely stiff reaction terms and stiff stochastic terms. The diffusion terms are solved by the explicit second order orthogonal Chebyshev method (ROCK2), while the stiff reaction terms (solved implicitly) and the advection and noise terms (solved explicitly) are integrated in the algorithm as finishing procedures. It is shown that the various coupling (between diffusion, reaction, advection and noise) can be stabilized in the PIROCK method. The method, implemented in a single black-box code that is fully adaptive, provides error estimators for the various terms present in the problem, and requires from the user solely the right-hand side of the differential equation. Numerical experiments and comparisons with existing Chebyshev methods, IMEX methods and partitioned methods show the efficiency and flexibility of our new algorithm.

Mots clés: orthogonal Runge-Kutta Chebyshev methods (ROCK); stabilized second-order integration method; partitioned Runge-Kutta methods.; partitioned Runge-Kutta methods

Domaines: Analyse numérique [math.NA]

Assyr Abdulle, Gilles Vilmart. PIROCK: a swiss-knife partitioned implicit-explicit orthogonal Runge-Kutta Chebyshev integrator for stiff diffusion-advection-reaction problems with or without noise. Journal of Computational Physics, 2013, 242, pp.869-888. ⟨10.1016/j.jcp.2013.02.009⟩. ⟨hal-00739757v2⟩