Exploiting Superconvergence Through Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering

There has been much work in the area of superconvergent error analysis for finite element and discontinuous Galerkin (DG) methods. The property of superconvergence leads to the question of how to exploit this information in a useful manner, mainly through superconvergence extraction. There are many methods used for superconvergence extraction such as projection, interpolation, patch recovery and B-spline convolution filters. This last method falls under the class of Smoothness-Increasing Accuracy-Conserving (SIAC) filters. It has the advantage of improving both smoothness and accuracy of the approximation. Specifically, for linear hyperbolic equations it can improve the order of accuracy of a DG approximation from k + 1 to 2k + 1, where k is the highest degree polynomial used in the approximation, and can increase the smoothness to k − 1. In this article, we discuss the importance of overcoming the mathematical barriers in making superconvergence extraction techniques useful for applications, specifically focusing on SIAC filtering

SIAC Filtering for Nonlinear Hyperbolic Equations

We present the results of the symmetric and one-sided Smoothness-Increasing Accuracy-Conserving (SIAC) filter applied to a discontinuous Galerkin (DG) approximation for two examples of nonlinear hyperbolic conservation laws. The traditional symmetric SIAC filter relies on having a translation invariant mesh, periodic boundary conditions and linear equations. However, for practical applications that are modelled by nonlinear hyperbolic equations, this is not feasible. Instead we must concentrate on a filter that allows error reduction for nonuniform/unstructured meshes and non-periodic boundary conditions for nonlinear hyperbolic equations. This proceedings is an introductory exploration into the feasibility of these requirements for efficient filtering of nonlinear equations

A modified particle method for semilinear hyperbolic systems with oscillatory solutions

We introduce a modified particle method for semi-linear hyperbolic systems with highly oscillatory solutions. The main feature of this modified particle method is that we do not require different families of characteristics to meet at one point. In the modified particle method, we update the ith component of the solution along its own characteristics, and interpolate the other components of the solution from their own characteristic points to the ith characteristic point. We prove the convergence of the modified particle method essentially independent of the small scale for the variable coefficient Carleman model. The same result also applies to the non-resonant Broadwell model. Numerical evidence suggests that the modified particle method also converges essentially independent of the small scale for the original Broadwell model if a cubic spline interpolation is used

A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations

Recent work by McClarren & Hauck [29] suggests that the filtered spherical harmonics method represents an efficient, robust, and accurate method for radiation transport, at least in the two-dimensional (2D) case. We extend their work to the three-dimensional (3D) case and find that all of the advantages of the filtering approach identified in 2D are present also in the 3D case. We reformulate the filter operation in a way that is independent of the timestep and of the spatial discretization. We also explore different second- and fourth-order filters and find that the second-order ones yield significantly better results. Overall, our findings suggest that the filtered spherical harmonics approach represents a very promising method for 3D radiation transport calculations.Comment: 29 pages, 13 figures. Version matching the one in Journal of Computational Physic