Stability of Correction Procedure via Reconstruction With Summation-by-Parts Operators for Burgers' Equation Using a Polynomial Chaos Approach
For researchers in uncertainty quantification and numerical methods, this work provides a stability proof and entropy stable fluxes for CPR applied to PC systems, but the findings are incremental as they extend existing CPR stability theory to a specific equation.
The paper proves stability of the correction procedure via reconstruction (CPR) with summation-by-parts operators for the polynomial chaos (PC) system derived from Burgers' equation with uncertain boundary/initial conditions, and constructs entropy stable numerical fluxes for all PC systems for the first time. Numerical tests show that different methods (finite volume, finite difference, CPR) yield significantly different solutions for an initial shock, attributed to high sensitivity to dissipation.
In this paper, we consider Burgers' equation with uncertain boundary and initial conditions. The polynomial chaos (PC) approach yields a hyperbolic system of deterministic equations, which can be solved by several numerical methods. Here, we apply the correction procedure via reconstruction (CPR) using summation-by-parts operators. We focus especially on stability, which is proven for CPR methods and the systems arising from the PC approach. Due to the usage of split-forms, the major challenge is to construct entropy stable numerical fluxes. For the first time, such numerical fluxes are constructed for all systems resulting from the PC approach for Burgers' equation. In numerical tests, we verify our results and show also the advantage of the given ansatz using CPR methods. Moreover, one of the simulations, i.e. Burgers' equation equipped with an initial shock, demonstrates quite fascinating observations. The behaviour of the numerical solutions from several methods (finite volume, finite difference, CPR) differ significantly from each other. Through careful investigations, we conclude that the reason for this is the high sensitivity of the system to varying dissipation. Furthermore, it should be stressed that the system is not strictly hyperbolic with genuinely nonlinear or linearly degenerate fields.