Multi-dimensional filtering: Reducing the dimension through rotation

Julia Docampo, Jennifer K. Ryan, Mahsa Mirzargar, Robert M. Kirby

    Research output: Contribution to journalArticlepeer-review

    15 Downloads (Pure)

    Abstract

    Over the past few decades there has been a strong effort towards the development of Smoothness-Increasing Accuracy-Conserving (SIAC) filters for Discontinuous Galerkin (DG) methods, designed to increase the smoothness and improve the convergence rate of the DG solution through this post-processor. These advantages can be exploited during flow visualization, for example by applying the SIAC filter to the DG data before streamline computations [Steffan et al., IEEE-TVCG 14(3): 680-692]. However, introducing these filters in engineering applications can be challenging since a tensor product filter grows in support size as the field dimension increases, becoming computationally expensive. As an alternative, [Walfisch et al., JOMP 38(2);164-184] proposed a univariate filter implemented along the streamline curves. Until now, this technique remained a numerical experiment. In this paper we introduce the line SIAC filter and explore how the orientation, structure and filter size affect the order of accuracy and global errors. We present theoretical error estimates showing how line filtering preserves the properties of traditional tensor product filtering, including smoothness and improvement in the convergence rate. Furthermore, numerical experiments are included, exhibiting how these filters achieve the same accuracy at significantly lower computational costs, becoming an attractive tool for the scientific visualization
    community.
    Original languageEnglish
    Pages (from-to)A2179–A2200
    JournalSIAM Journal on Scientific Computing
    Volume39
    Issue number5
    DOIs
    Publication statusPublished - 27 Sept 2017

    Keywords

    • Discontinuous Galerkin
    • post-processing
    • SIAC filtering
    • accuracy enhancement
    • error reduction

    Cite this