English

A hybrid SIAC -- data-driven post-processing filter for discontinuities in solutions to numerical PDEs

Numerical Analysis 2025-03-17 v2 Numerical Analysis

Abstract

We present a hybrid filter that is only applied to the approximation at the final time and allows for reducing errors away from a shock as well as near a shock. It is designed for discontinuous Galerkin approximations to PDEs and combines a rigorous moment-based Smoothness-Increasing Accuracy-Conserving (SIAC) filter with a data-driven CNN filter. While SIAC improves accuracy in smooth regions, it fails to reduce the O(1)\mathcal{O}(1) errors near discontinuities, particularly in inviscid compressible flows with shocks. Our hybrid SIAC-CNN filter, trained exclusively on top-hat functions, enforces consistency constraints globally and higher-order moment conditions in smooth regions, reducing both 2\ell_2 and \ell_\infty errors near discontinuities and preserving theoretical accuracy in smooth regions. We demonstrate its effectiveness on the Euler equations for the Lax, Sod, and Shu-Osher shock-tube problems.

Cite

@article{arxiv.2408.05193,
  title  = {A hybrid SIAC -- data-driven post-processing filter for discontinuities in solutions to numerical PDEs},
  author = {Soraya Terrab and Samy Wu Fung and Jennifer K. Ryan},
  journal= {arXiv preprint arXiv:2408.05193},
  year   = {2025}
}
R2 v1 2026-06-28T18:08:50.766Z