Oscillation mitigation of hyperbolicity-preserving intrusive uncertainty quantification methods for systems of conservation laws (Q2229920)

Uncertainties play a key role in modeling hyperbolic systems of conservation laws if crucial data or parameters might not be available exactly due to measurement errors. Thus, to achieve an adequate description of reality, one needs to include non-deterministic effects in the approximation of deterministic systems. Many approaches, such as methods based on Bayesian inversion, Monte Carlo algorithms or stochastic Galerkin schemes, have been proposed to quantify the uncertainties in order to account for them in predictions and simulation. In the context of these Uncertainty Quantification (UQ) methods, there are so-called non-intrusive and intrusive approaches. In this article, the authors focus on intrusive UQ methods, which aim to increase the overall efficiency but require a modification of the underlying solver for the deterministic problem. The biggest challenge of UQ methods for hyperbolic conservation laws lies in the fact that discontinuities in the physical space propagate into the solution manifold and cause oscillations which translate into nonphysical step-like approximations of the expected value around shocks. The authors propose two strategies to ensure both the mitigation of oscillations and the preservation of hyperbolicity by employing the multi-element ansatz together with the hyperbolicity limiter and filtering to stochastic Galerkin, and by deriving an extension of the IPM method to multi-element basis functions. Both strategies are extended to adaptivity, allowing to adapt the number of basis functions in each multi-element to the smoothness of the solution. They evaluate and compare both approaches on numerical examples such as a NACA airfoil and a nozzle test case for the two-dimensional Euler equations and observe the mitigation of spurious artifacts. Using the multi-element ansatz for IPM also significantly reduces computational costs.