A combination of Residual Distribution and the Active Flux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1D Euler equations

Rémi Abgrall

DOI: https://doi.org/10.48550/arXiv.2011.12572

2021-10-26

Abstract:We show how to combine in a natural way (i.e. without any test nor switch) the conservative and non conservative formulations of an hyperbolic system that has a conservative form. This is inspired from two different class of schemes: the Residual Distribution one \cite{MR4090481}, and the Active Flux formulations \cite{AF1, AF3, AF4,AF5,RoeAF}. The solution is globally continuous, and as in the active flux method, described by a combination of point values and average values. Unlike the "classical" active flux methods, the meaning of the pointwise and cella averaged degrees of freedom is different, and hence follow different form of PDEs: it is a conservative version of the cell average, and a possibly non conservative one for the points. This new class of scheme is proved to satisfy a Lax-Wendroff like theorem. We also develop a method to perform non linear stability. We illustrate the behaviour on several benchmarks, some quite challenging.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to naturally combine the conservative and non - conservative forms of hyperbolic systems to solve the Euler equations in fluid mechanics. Specifically, the author proposes a new numerical method that can handle partial differential equations (PDEs) in both conservative and non - conservative forms simultaneously without any testing or switching. This method combines the advantages of the Residual Distribution (RD) framework and the Active Flux method.

### Main problems and goals

1. **Combining conservative and non - conservative forms**:
- The author aims to show how to naturally combine the conservative and non - conservative forms of hyperbolic systems without any testing or switching.
- This method can maintain the continuity of the solution globally and describe the solution by a combination of point values and cell - averaged values.

2. **Satisfying the Lax - Wendroff theorem**:
- The new method is proven to satisfy a result similar to the Lax - Wendroff theorem, ensuring convergence to the correct weak solution in the limit of mesh refinement.

3. **Non - linear stability**:
- The author develops a method for non - linear stability analysis to ensure the stability of the numerical solution when dealing with discontinuities and high - gradient regions.

4. **Application and verification**:
- Verify the effectiveness of the new method through multiple benchmark tests, including some challenging problems such as the Sod shock - tube problem, the Shu - Osher problem, and the Le Blanc problem, etc.

### Mathematical models

The following forms of the Euler equations are considered in the paper:

- **Conservative form**:
\[
\frac{\partial}{\partial t} \begin{pmatrix} \rho \\ \rho u \\ E \end{pmatrix} + \frac{\partial}{\partial x} \begin{pmatrix} \rho u \\ \rho u^2 + p \\ u(E + p) \end{pmatrix} = 0
\]

- **Primitive variable form**:
\[
\frac{\partial}{\partial t} \begin{pmatrix} \rho \\ u \\ p \end{pmatrix} + \begin{pmatrix} \frac{\partial (\rho u)}{\partial x} \\ u \frac{\partial u}{\partial x} + \frac{1}{\rho} \frac{\partial p}{\partial x} \\ u \frac{\partial p}{\partial x} + (e + p) \frac{\partial u}{\partial x} \end{pmatrix} = 0
\]

- **Entropy form**:
\[
\frac{\partial}{\partial t} \begin{pmatrix} p \\ u \\ s \end{pmatrix} + \begin{pmatrix} u \frac{\partial p}{\partial x} + (e + p) \frac{\partial u}{\partial x} \\ u \frac{\partial u}{\partial x} + \frac{1}{\rho} \frac{\partial p}{\partial x} \\ u \frac{\partial s}{\partial x} \end{pmatrix} = 0
\]

where \(\rho\) is the density, \(u\) is the velocity, \(p\) is the pressure, and \(E = e+\frac{1}{2}\rho u^{2}\)