A hybrid SIAC -- data-driven post-processing filter for discontinuities in solutions to numerical PDEs
<|reference_start|>A hybrid SIAC -- data-driven post-processing filter for discontinuities in solutions to numerical PDEs: We introduce a hybrid filter that incorporates a mathematically accurate moment-based filter with a data driven filter for discontinuous Galerkin approximations to PDE solutions that contain discontinuities. Numerical solutions of PDEs suffer from an $\mathcal{O}(1)$ error in the neighborhood about discontinuities, especially for shock waves that arise in inviscid compressible flow problems. While post-processing filters, such as the Smoothness-Increasing Accuracy-Conserving (SIAC) filter, can improve the order of error in smooth regions, the $\mathcal{O}(1)$ error in the vicinity of a discontinuity remains. To alleviate this, we combine the SIAC filter with a data-driven filter based on Convolutional Neural Networks. The data-driven filter is specifically focused on improving the errors in discontinuous regions and therefore {\it only includes top-hat functions in the training dataset.} For both filters, a consistency constraint is enforced, while the SIAC filter additionally satisfies $r-$moments. This hybrid filter approach allows for maintaining the accuracy guaranteed by the theory in smooth regions while the hybrid SIAC-data-driven approach reduces the $\ell_2$ and $\ell_\infty$ errors about discontinuities. Thus, overall the global errors are reduced. We examine the performance of the hybrid filter about discontinuities for the one-dimensional Euler equations for the Lax, Sod, and sine-entropy (Shu-Osher) shock-tube problems.<|reference_end|>
@article{terrab2024a,
    title={A hybrid SIAC -- data-driven post-processing filter for discontinuities
  in solutions to numerical PDEs},
    author={Soraya Terrab, Samy Wu Fung, Jennifer K. Ryan},
    journal={arXiv preprint arXiv:2408.05193},
    year={2024},
    archivePrefix={arXiv},
    eprint={2408.05193},
    primaryClass={math.NA cs.NA}
}