A flux reconstruction stochastic Galerkin scheme for hyperbolic conservation laws

Abstract

The study of uncertainty propagation poses a great challenge to design numerical solvers with high fidelity. Based on the stochastic Galerkin formulation, this paper addresses the idea and implementation of the first flux reconstruction scheme for hyperbolic conservation laws with random inputs. Unlike the finite volume method, the treatments in physical and random space are consistent, e.g., the modal representation of solutions based on an orthogonal polynomial basis and the nodal representation based on solution collocation points. Therefore, the numerical behaviors of the scheme in the phase space can be designed and understood uniformly. A family of filters is extended to multi-dimensional cases to mitigate the well-known Gibbs phenomenon arising from discontinuities in both physical and random space. The filter function is switched on and off by the dynamic detection of discontinuous solutions, and a slope limiter is employed to preserve the positivity of physically realizable solutions. As a result, the proposed method is able to capture stochastic cross-scale flow evolution where resolved and unresolved regions coexist. Numerical experiments including wave propagation, Burgers' shock, one-dimensional Riemann problem, and two-dimensional shock-vortex interaction problem are presented to validate the scheme. The order of convergence of the current scheme is identified. The capability of the scheme for simulating smooth and discontinuous stochastic flow dynamics is demonstrated. The open-source codes to reproduce the numerical results are available under the MIT license [1].