Header

UZH-Logo

Maintenance Infos

High order asymptotic preserving deferred correction implicit-explicit schemes for kinetic models


Abgrall, Rémi; Torlo, Davide (2020). High order asymptotic preserving deferred correction implicit-explicit schemes for kinetic models. SIAM Journal on Scientific Computing, 42(3):B816-B845.

Abstract

This work introduces an extension of the residual distribution (RD) framework to stiff relaxation problems. The RD is a class of schemes which is used to solve a hyperbolic system of partial differential equations. To our knowledge, it has been used only for systems with mild source terms, such as gravitation problems or shallow water equations. What we propose is an implicit-explicit (IMEX) version of the RD schemes that can resolve stiff source terms, without refining the discretization up to the stiffness scale. This can be particularly useful in various models, where the stiffness is given by topological or physical quantities, e.g., multiphase flows, kinetic models, or viscoelasticity problems. We will focus on kinetic models that are BGK approximation of hyperbolic conservation laws. The extension to more complicated problems will be carried out in future works. The provided scheme is able to catch different relaxation scales automatically, without losing accuracy; we prove that the scheme is asymptotic preserving and this guarantees that, in the relaxation limit, we recast the expected macroscopic behavior. To get a high order accuracy, we use an IMEX time discretization combined with a deferred correction procedure, while naturally RD provides high order space discretization. Finally, we show some numerical tests in one and two dimensions for stiff systems of equations.

Abstract

This work introduces an extension of the residual distribution (RD) framework to stiff relaxation problems. The RD is a class of schemes which is used to solve a hyperbolic system of partial differential equations. To our knowledge, it has been used only for systems with mild source terms, such as gravitation problems or shallow water equations. What we propose is an implicit-explicit (IMEX) version of the RD schemes that can resolve stiff source terms, without refining the discretization up to the stiffness scale. This can be particularly useful in various models, where the stiffness is given by topological or physical quantities, e.g., multiphase flows, kinetic models, or viscoelasticity problems. We will focus on kinetic models that are BGK approximation of hyperbolic conservation laws. The extension to more complicated problems will be carried out in future works. The provided scheme is able to catch different relaxation scales automatically, without losing accuracy; we prove that the scheme is asymptotic preserving and this guarantees that, in the relaxation limit, we recast the expected macroscopic behavior. To get a high order accuracy, we use an IMEX time discretization combined with a deferred correction procedure, while naturally RD provides high order space discretization. Finally, we show some numerical tests in one and two dimensions for stiff systems of equations.

Statistics

Citations

Altmetrics

Downloads

12 downloads since deposited on 28 Sep 2020
12 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Institute of Mathematics
Dewey Decimal Classification:340 Law
610 Medicine & health
510 Mathematics
Scopus Subject Areas:Physical Sciences > Computational Mathematics
Physical Sciences > Applied Mathematics
Uncontrolled Keywords:Applied Mathematics, Computational Mathematics
Language:English
Date:1 January 2020
Deposited On:28 Sep 2020 14:04
Last Modified:29 Sep 2020 20:00
Publisher:Society for Industrial and Applied Mathematics
ISSN:1064-8275
Additional Information:Copyright © 2020, Society for Industrial and Applied Mathematics
OA Status:Green
Publisher DOI:https://doi.org/10.1137/19m128973x
Project Information:
  • : FunderH2020
  • : Grant ID642768
  • : Project TitleModCompShock - Modelling and Computation of Shocks and Interfaces

Download

Green Open Access

Download PDF  'High order asymptotic preserving deferred correction implicit-explicit schemes for kinetic models'.
Preview
Content: Published Version
Language: English
Filetype: PDF
Size: 1MB
View at publisher