We consider implicit-explicit (IMEX) Runge--Kutta (R-K) schemes for hyperbolic systems with stiff relaxation in the so-called diffusion limit. In such a regime the system relaxes towards a convection-diffusion equation. The first objective of this paper is to show that traditional partitioned IMEX R-K schemes will relax to an explicit scheme for the limit equation with no need of modification of the original system. Of course the explicit scheme obtained in the limit suffers from the classical parabolic stability restriction on the time step. The main goal of this paper is to present an approach, based on IMEX R-K schemes, that in the diffusion limit relaxes to an IMEX R-K scheme for the convection-diffusion equation, in which the diffusion is treated implicitly. This is achieved by a novel reformulation of the problem, and subsequent application of IMEX R-K schemes to it. An analysis of such schemes to the reformulated problem shows that the schemes reduce to IMEX R-K schemes for the limit equation, under the same conditions derived for hyperbolic relaxation [S. Boscarino and G. Russo, SIAM J. Sci. Comput., 31 (2009), pp. 1926--1945]. Several numerical examples including neutron transport equations confirm the theoretical analysis.

Implicit-Explicit Runge--Kutta Schemes for Hyperbolic Systems and Kinetic Equations in the Diffusion Limit

BOSCARINO, SEBASTIANO;RUSSO, Giovanni
2013-01-01

Abstract

We consider implicit-explicit (IMEX) Runge--Kutta (R-K) schemes for hyperbolic systems with stiff relaxation in the so-called diffusion limit. In such a regime the system relaxes towards a convection-diffusion equation. The first objective of this paper is to show that traditional partitioned IMEX R-K schemes will relax to an explicit scheme for the limit equation with no need of modification of the original system. Of course the explicit scheme obtained in the limit suffers from the classical parabolic stability restriction on the time step. The main goal of this paper is to present an approach, based on IMEX R-K schemes, that in the diffusion limit relaxes to an IMEX R-K scheme for the convection-diffusion equation, in which the diffusion is treated implicitly. This is achieved by a novel reformulation of the problem, and subsequent application of IMEX R-K schemes to it. An analysis of such schemes to the reformulated problem shows that the schemes reduce to IMEX R-K schemes for the limit equation, under the same conditions derived for hyperbolic relaxation [S. Boscarino and G. Russo, SIAM J. Sci. Comput., 31 (2009), pp. 1926--1945]. Several numerical examples including neutron transport equations confirm the theoretical analysis.
2013
IMEX Runge--Kutta methods; hyperbolic conservation laws with sources; diffusion equations, stiff systems
File in questo prodotto:
File Dimensione Formato  
Bosc_Par_Rus2013.pdf

solo gestori archivio

Tipologia: Versione Editoriale (PDF)
Dimensione 366.26 kB
Formato Adobe PDF
366.26 kB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11769/39680
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 123
  • ???jsp.display-item.citation.isi??? 108
social impact