Bojan Popov

Jean-Luc Guermond, Murtazo Nazarov, Bojan Popov et al.

A new second-order method for approximating the compressible Euler equations is introduced. The method preserves all the known invariant domains of the Euler system: positivity of the density, positivity of the internal energy and the local minimum principle on the specific entropy. The technique combines a first-order, invariant domain preserving, Guaranteed Maximum Speed method using a Graph Viscosity (GMS-GV1) with an invariant domain violating, but entropy consistent, high-order method. Invariant domain preserving auxiliary states, naturally produced by the GMS-GV1 method, are used to define local bounds for the high-order method which is then made invariant domain preserving via a convex limiting process. Numerical tests confirm the second-order accuracy of the new GMS-GV2 method in the maximum norm, where 2 stands for second-order. The proposed convex limiting is generic and can be applied to other approximation techniques and other hyperbolic systems.

Jean-Luc Guermond, Bojan Popov, Ignacio Tomas

We introduce an approximation technique for nonlinear hyperbolic systems with sources that is invariant domain preserving. The method is discretization-independent provided elementary symmetry and skew-symmetry properties are satisfied by the scheme. The method is formally first-order accurate in space. A series of higher-order methods is also introduced. When these methods violate the invariant domain properties, they are corrected by a limiting technique that we call convex limiting. After limiting, the resulting methods satisfy all the invariant domain properties that are imposed by the user (see Theorem~7.24). A key novelty is that the bounds that are enforced on the solution at each time step are necessarily satisfied by the low-order approximation.

Jean-Luc Guermond, Bojan Popov

This paper is concerned with the construction of a fast algorithm for computing the maximum speed of propagation in the Riemann solution for the Euler system of gas dynamics with the co-volume equation of state. The novelty in the algorithm is that it stops when a guaranteed upper bound for the maximum speed is reached with a prescribed accuracy. The convergence rate of the algorithm is cubic and the bound is guaranteed for gasses with the co-volume equation of state and the heat capacity ratio $γ$ in the range $(1,5/3]$

Jean-Luc Guermond, Bojan Popov, Jean Ragusa

We introduce a (linear) positive and asymptotic preserving method or solving the one-group radiation transport equation. The approximation in space is discretization agnostic: the space approximation can be done with continuous or discontinuous finite elements (or finite volumes, or finite differences). The method is first-order accurate in space. This type of accuracy is coherent with Godunov's theorem since the method is linear. The two key theoretical results of the paper are Theorem~4.4 and Theorem~4.8. The method is illustrated with continuous finite elements. It is observed to converge with the rate $\calO(h)$ in the $L^2$-norm on manufactured solutions, and it is $\calO(h^2)$ in the diffusion regime. Unlike other standard techniques, the proposed method does not suffer from overshoots at the interfaces of optically thin and optically thick regions.

Jean-Luc Guermond, Bojan Popov

We propose a numerical method to solve general hyperbolic systems in any space dimension using forward Euler time stepping and continuous finite elements on non-uniform grids. The properties of the method are based on the introduction of an artificial dissipation that is defined so that any convex invariant sets containing the initial data is an invariant domain for the method. The invariant domain property is proved for any hyperbolic system provided a CFL condition holds. The solution is also shown to satisfy a discrete entropy inequality for every admissible entropy of the system. The method is formally first-order accurate in space and can be made high-order in time by using Strong Stability Preserving algorithms. This technique extends to continuous finite elements the work of \cite{Hoff_1979,Hoff_1985}, and \cite{Frid_2001}.