SIAC Accuracy Enhancement of Stochastic Galerkin Solutions for Wave Equations with Uncertain Coefficients

Abstract

This article establishes the usefulness of the Smoothness-Increasing Accuracy-Increasing (SIAC) filter for reducing the errors in the mean and variance for a wave equation with uncertain coefficients solved via generalized polynomial chaos (gPC) whose coefficients are approximated using discontinuous Galerkin (DG-gPC). Theoretical error estimates that utilize information in the negative-order norm are established. While the gPC approximation leads to order of accuracy of $m-1/2$ for a sufficiently smooth solution (smoothness of $m$ in random space), the approximated coefficients solved via DG improves from order $k+1$ to $2k+1$ for a solution of smoothness $2k+2$ in physical space. Our numerical examples verify the performance of the filter for improving the quality of the approximation and reducing the numerical error and significantly eliminating the noise from the spatial approximation of the mean and variance. Further, we illustrate how the errors are effected by both the choice of smoothness of the kernel and number of function translates in the kernel. Hence, this article opens the applicability of SIAC filters to other hyperbolic problems with uncertainty, and other stochastic equations.

19 pages, 3 figures