skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Abstract

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.

Authors:
 [1];  [2]
  1. Mathematics Section, École Polytechnique Fédérale de Lausanne, Station 8, 1015 Lausanne (Switzerland)
  2. École Normale Supérieure de Cachan, Antenne de Bretagne, INRIA Rennes, IRMAR, CNRS, UEB, Av. Robert Schuman, F-35170 Bruz (France)
Publication Date:
OSTI Identifier:
22233591
Resource Type:
Journal Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 242; Other Information: Copyright (c) 2013 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0021-9991
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICAL METHODS AND COMPUTING; ADVECTION; ALGORITHMS; COMPARATIVE EVALUATIONS; DIFFERENTIAL EQUATIONS; DIFFUSION; EFFICIENCY; ERRORS; FLEXIBILITY; NOISE; PARTITION; STOCHASTIC PROCESSES

Citation Formats

Abdulle, Assyr, and Vilmart, Gilles. PIROCK: A swiss-knife partitioned implicit–explicit orthogonal Runge–Kutta Chebyshev integrator for stiff diffusion–advection–reaction problems with or without noise. United States: N. p., 2013. Web. doi:10.1016/J.JCP.2013.02.009.
Abdulle, Assyr, & Vilmart, Gilles. PIROCK: A swiss-knife partitioned implicit–explicit orthogonal Runge–Kutta Chebyshev integrator for stiff diffusion–advection–reaction problems with or without noise. United States. https://doi.org/10.1016/J.JCP.2013.02.009
Abdulle, Assyr, and Vilmart, Gilles. Sat . "PIROCK: A swiss-knife partitioned implicit–explicit orthogonal Runge–Kutta Chebyshev integrator for stiff diffusion–advection–reaction problems with or without noise". United States. https://doi.org/10.1016/J.JCP.2013.02.009.
@article{osti_22233591,
title = {PIROCK: A swiss-knife partitioned implicit–explicit orthogonal Runge–Kutta Chebyshev integrator for stiff diffusion–advection–reaction problems with or without noise},
author = {Abdulle, Assyr and Vilmart, Gilles},
abstractNote = {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.},
doi = {10.1016/J.JCP.2013.02.009},
url = {https://www.osti.gov/biblio/22233591}, journal = {Journal of Computational Physics},
issn = {0021-9991},
number = ,
volume = 242,
place = {United States},
year = {2013},
month = {6}
}