Manuel Torrilhon

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.

Birte Schmidtmann, Benjamin Seibold, Manuel Torrilhon

Weighted essentially non-oscillatory (WENO) and finite volume (FV) methods employ different philosophies in their way to perform limiting. We show that a generalized view on limiter functions, which considers a two-dimensional, rather than a one-dimensional dependence on the slopes in neighboring cells, allows to write WENO3 and $3^\text{rd}$-order FV schemes in the same fashion. Within this framework, it becomes apparent that the classical approach of FV limiters to only consider ratios of the slopes in neighboring cells, is overly restrictive. The hope of this new perspective is to establish new connections between WENO3 and FV limiter functions, which may give rise to improvements for the limiting behavior in both approaches.

Mohsen Sadr, Manuel Torrilhon, M. Hossein Gorji

Maximum-Entropy Distributions offer an attractive family of probability densities suitable for moment closure problems. Yet finding the Lagrange multipliers which parametrize these distributions, turns out to be a computational bottleneck for practical closure settings. Motivated by recent success of Gaussian processes, we investigate the suitability of Gaussian priors to approximate the Lagrange multipliers as a map of a given set of moments. Examining various kernel functions, the hyperparameters are optimized by maximizing the log-likelihood. The performance of the devised data-driven Maximum-Entropy closure is studied for couple of test cases including relaxation of non-equilibrium distributions governed by Bhatnagar-Gross-Krook and Boltzmann kinetic equations.

Birte Schmidtmann, Mariia Astrakhantceva, Manuel Torrilhon

In this paper we present a new family of approximate Riemann solvers for the numerical approximation of solutions of hyperbolic conservation laws. They are approximate, also referred to as incomplete, in the sense that the solvers avoid computing the characteristic decomposition of the flux Jacobian. Instead, they require only an estimate of the globally fastest wave speeds in both directions. Thus, this family of solvers is particularly efficient for large systems of conservation laws, i.e. with many different propagation speeds, and when no explicit expression for the eigensystem is available. Even though only fastest wave speeds are needed as input values, the new family of Riemann solvers reproduces all waves with less dissipation than HLL, which has the same prerequisites, requiring only one additional flux evaluation.

Peter Lewintan, Lambert Theisen, Manuel Torrilhon

In this paper, we finally prove the well-posedness of the linearized R13 moment model, which describes, e.g., rarefied gas flows. As an extension of the classical fluid equations, moment models are robust and have been frequently used, yet they are challenging to analyze due to their additional equations. By effectively grouping variables, we identify a 2-by-2 block structure, allowing us to analyze well-posedness within the abstract LBB framework for saddle point problems. Due to the unique tensorial structure of the equations, in addition to an interesting combination of tools from Stokes' and linear elasticity theory, we also need new coercivity estimates for tensor fields. These Korn-type inequalities are established by analyzing the symbol map of the symmetric and trace-free part of tensor derivative fields. Together with the corresponding right inverse of the tensorial divergence, we obtain the existence and uniqueness of weak solutions. This result also serves as the basis for future numerical analysis of corresponding discretization schemes.

Daniel Doehring, Jesse Chan, Hendrik Ranocha et al.

We introduce the concept of volume term adaptivity for high-order discontinuous Galerkin (DG) schemes solving time-dependent partial differential equations. Termed v-adaptivity, we present a novel general approach that exchanges the discretization of the volume contribution of the DG scheme at every Runge-Kutta stage based on suitable indicators. Depending on whether robustness or efficiency is the main concern, different adaptation strategies can be chosen. Precisely, the weak form volume term discretization is used instead of the entropy-conserving flux-differencing volume integral whenever the former produces more entropy than the latter, resulting in an entropy-stable scheme. Conversely, if increasing the efficiency is the main objective, the weak form volume integral may be employed as long as it does not increase entropy beyond a certain threshold or cause instabilities. Thus, depending on the choice of the indicator, the v-adaptive DG scheme improves robustness, efficiency and approximation quality compared to schemes with a uniform volume term discretization. We thoroughly verify the accuracy, linear stability, and entropy-admissibility of the v-adaptive DG scheme before applying it to various compressible flow problems in two and three dimensions.

Eda Yilmaz, Georgii Oblapenko, Manuel Torrilhon

In the present work, an approach to the moment closure problem on the basis of orthogonal polynomials derived from Gram matrices is proposed. Its properties are studied in the context of the moment closure problem arising in gas kinetic theory, for which the proposed approach is proven to have multiple attractive mathematical properties. Numerical studies are carried out for model gas particle distributions and the approach is compared to other moment closure methods, such as Grad's closure and the maximum-entropy method. The proposed ``Gramian'' closure is shown to provide very accurate results for a wide range of distribution functions.

Birte Schmidtmann, Manuel Torrilhon

We are interested in the numerical solution of large systems of hyperbolic conservation laws or systems in which the characteristic decomposition is expensive to compute. Solving such equations using finite volumes or Discontinuous Galerkin requires a numerical flux function which solves local Riemann problems at cell interfaces. There are various methods to express the numerical flux function. On the one end, there is the robust but very diffusive Lax-Friedrichs solver; on the other end the upwind Godunov solver which respects all resulting waves. The drawback of the latter method is the costly computation of the eigensystem. This work presents a family of simple first order Riemann solvers, named HLLX$ω$, which avoid solving the eigensystem. The new method reproduces all waves of the system with less dissipation than other solvers with similar input and effort, such as HLL and FORCE. The family of Riemann solvers can be seen as an extension or generalization of the methods introduced by Degond et al. \cite{DegondPeyrardRussoVilledieu1999}. We only require the same number of input values as HLL, namely the globally fastest wave speeds in both directions, or an estimate of the speeds. Thus, the new family of Riemann solvers is particularly efficient for large systems of conservation laws when the spectral decomposition is expensive to compute or no explicit expression for the eigensystem is available.