Rémi Abgrall

Birte Schmidtmann, Rémi Abgrall, Manuel Torrilhon

In this article, we propose a finite volume limiter function for a reconstruction on the three-point stencil. Compared to classical limiter functions in the MUSCL framework, which yield $2^{\text{nd}}$-order accuracy, the new limiter is $3^\text{rd}$-order accurate for smooth solutions. In an earlier work, such a $3^\text{rd}$-order limiter function was proposed and showed successful results [2]. However, it came with unspecified parameters. We close this gap by giving information on these parameters.

Rémi Abgrall, Davide Torlo

This work aims to extend the residual distribution (RD) framework to stiff relaxation problems. The RD is a class of schemes which is used to solve hyperbolic system of partial differential equations. Up to our knowledge, it was used only for systems with mild source terms, such as gravitation problems or shallow water equations. What we propose is an IMEX (implicit--explicit) version of the residual distribution 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, viscoelasticity problems. Moreover, the provided scheme is able to catch different relaxation scales automatically, without losing accuracy. The scheme is asymptotic preserving and this guarantees that in the relaxation limit, we recast the expected macroscopic behaviour. To get a high order accuracy, we use an IMEX time discretization combined with a Deferred Correction (DeC) procedure, while naturally RD provides high order space discretization. Finally, we show some numerical tests in 1D and 2D for stiff systems of equations.

Rémi Abgrall, Philipp Öffner, Yongle Liu

In this paper, we provide a few new properties of Active Flux (AF)/Point-Average-Moment PolynomiAl-interpreted (\pampa) schemes. First, we show, in full generality, that the AF/pampa schemes can be interpreted in such a way that the discontinuous Galerkin (dG) scheme is one of their building blocks. Secondly we provide intrinsic bound preserving properties of the current variant of pampa. This is also illustrated numerically. Last, we show, at least in one dimension, that the pampa scheme has the summation by part (SBP) property.

Rémi Abgrall, Michael Dumbser, Pierre-Henri Maire et al.

We present a new class of structure-preserving semi-discrete continuous-discontinuous Galerkin (CG-DG) finite element schemes for linear and nonlinear hyperbolic systems of partial differential equations on unstructured simplex meshes that automatically satisfy the following properties: i) the new schemes are not only cellwise conservative, but also locally pointwise conservative everywhere, hence they satisfy the integral form of the conservation law on arbitrary control volumes that do not have to coincide with the mesh at all; ii) the new methods naturally satisfy the two basic vector calculus identities $\nabla \cdot \nabla \times \mathbf{A}$ and $\nabla \times \nabla Z$ exactly pointwise locally and globally everywhere on the discrete level; iii) for linear symmetric hyperbolic systems the schemes are naturally energy conservative for the square energy, i.e. nonlinearly stable in the $L^2$ norm. The key ingredient of the new CG-DG schemes is the use of two different but compatible approximation spaces: the classical DG space $\mathcal{U}_h^N$ of discontinuous piecewise polynomials of degree up to $N$ and a classical finite element space $\mathcal{W}_h^{N+1}$ of globally continuous piecewise polynomials of degree $N+1$. In the new CG-DG schemes, the discrete solution $\mathbf{u}_h$ is sought in $\mathcal{U}_h^N$, while a suitable discrete flux field $\tilde{\mathbf{f}}_h$ is computed in $\mathcal{W}_h^{N+1}$. For $N=0$ our new schemes are directly related to cell-centered finite volume schemes with suitable vertex-based fluxes. All claimed properties of the schemes are first mathematically proven and are then also verified via suitable numerical tests. We show applications of our approach to three linear and nonlinear hyperbolic systems.

Jianfang Lin, Rémi Abgrall, Jianxian Qiu

In this paper, we propose a high order residual distribution conservative finite difference scheme for solving steady state conservation laws. A new type of WENO (weighted essentially non-oscillatory) termed as WENO-ZQ integration is used to compute the numerical fluxes and source term based on the point values of the solution, and the principles of residual distribution schemes are adapted to obtain steady state solutions. Extensive numerical examples in both scalar and system test problems in one and two dimensions demonstrate the efficiency, high order accuracy and the capability of resolving shocks of the proposed methods.