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

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:400

Article

Kusch, Jonas¹  Schlachter, Louisa² 

[1] Karlsruhe Inst Technol, Computat Sci & Math Methods, Englerstr 2, D-76128 Karlsruhe, Germany

[2] TU Kaiserslautern, Fachbereich Math, Erwin Schrodinger Str, D-67663 Kaiserslautern, Germany

关键词: Uncertainty quantification;    Polynomial Chaos;    Stochastic Galerkin;    Filter;    Hyperbolicity;    Multi-element;   

DOI  :  10.1016/j.cam.2021.113714

来源: Elsevier

【 摘 要 】

In this article we study intrusive uncertainty quantification schemes for systems of conservation laws with uncertainty. While intrusive methods inherit certain advantages such as adaptivity and an improved accuracy, they suffer from two key issues. First, intrusive methods tend to show oscillations, especially at shock structures and second, standard intrusive methods can lose hyperbolicity. The aim of this work is to tackle these challenges with the help of two different strategies. First, we combine filters with the multi-element approach for the hyperbolicity-preserving stochastic Galerkin (hSG) scheme. While the limiter used in the hSG scheme ensures hyperbolicity, the filter as well as the multi-element ansatz mitigate oscillations. Second, we derive a multi-element approach for the intrusive polynomial moment (IPM) method. Even though the IPM method is inherently hyperbolic, it suffers from oscillations while requiring the solution of an optimization problem in every spatial cell and every time step. The proposed multi-element IPM method leads to a decoupling of the optimization problem in every multi-element. Thus, we are able to significantly decrease computational costs while improving parallelizability. Both proposed strategies are extended to adaptivity, allowing to adapt the number of basis functions in each multi-element to the smoothness of the solution. We finally evaluate and compare both approaches on various numerical examples such as a NACA airfoil and a nozzle test case for the two-dimensional Euler equations. In our numerical experiments, we observe the mitigation of spurious artifacts. Furthermore, using the multi-element ansatz for IPM significantly reduces computational costs. (C) 2021 Elsevier B.V. All rights reserved.

【 授权许可】

Free