Smoothness-Increasing Accuracy-Conserving (SIAC) Postprocessing for Discontinuous Galerkin Solutions over Structured Triangular Meshes

Smoothness-Increasing Accuracy-Conserving (SIAC) Postprocessing for Discontinuous Galerkin Solutions over Structured Triangular Meshes

Mirzaee, H.
Li, L.
Ryan, J.K.
Kirby, R.M.

Electrical Engineering, Mathematics and Computer Science

Delft Institute of Applied Mathematics

2011-09-22

Theoretically and computationally, it is possible to demonstrate that the order of accuracy of a discontinuous Galerkin (DG) solution for linear hyperbolic equations can be improved from order $k$+1 to 2$k$+1 through the use of smoothness-increasing accuracy-conserving (SIAC) filtering. However, it is a computationally complex task to perform this in an efficient manner, which becomes an even greater issue considering nonquadrilateral mesh structures. In this paper, we present an extension of this SIAC filter to structured triangular meshes. The basic theoretical assumption in the previous implementations of the postprocessor limits the use to numerical solutions solved over a quadrilateral mesh. However, this assumption is restrictive, which in turn complicates the application of this postprocessing technique to general tessellations. Additionally, moving from quadrilateral meshes to triangulated ones introduces more complexity in the calculations as the number of integrations required increases. In this paper, we extend the current theoretical results to variable coefficient hyperbolic equations over structured triangular meshes and demonstrate the effectiveness of the application of this postprocessor to structured triangular meshes as well as exploring the effect of using inexact quadrature. We show that there is a direct theoretical extension to structured triangular meshes for hyperbolic equations with bounded variable coefficients. This is a challenging first step toward implementing SIAC filters for unstructured tessellations. We show that by using the usual B-spline implementation, we are able to improve on the order of accuracy as well as decrease the magnitude of the errors. These results are valid regardless of whether exact or inexact integration is used. The results here demonstrate that it is still possible, both theoretically and computationally, to improve to 2$k$+1 over the DG solution itself for structured triangular meshes.

high-order methods
discontinuous Galerkin
smoothness-increasing accuracy-conserving filtering
accuracy enhancement

http://resolver.tudelft.nl/uuid:735eafc0-aaca-453b-8baf-7e5e5e95a0f0

https://doi.org/10.1137/110830678

Society for Industrial and Applied Mathematics

0036-1429

SIAM Journal on Numerical Analysis, 49 (5), 2011

Institutional Repository

journal article

(c) 2011 Society for Industrial and Applied Mathematics