Robust Discontinuity Indicators for High-Order Reconstruction of Piecewise Smooth Functions

by   Yipeng Li, et al.

In many applications, piecewise continuous functions are commonly interpolated over meshes. However, accurate high-order manipulations of such functions can be challenging due to potential spurious oscillations known as the Gibbs phenomena. To address this challenge, we propose a novel approach, Robust Discontinuity Indicators (RDI), which can efficiently and reliably detect both C^0 and C^1 discontinuities for node-based and cell-averaged values. We present a detailed analysis focusing on its derivation and the dual-thresholding strategy. A key advantage of RDI is its ability to handle potential inaccuracies associated with detecting discontinuities on non-uniform meshes, thanks to its innovative discontinuity indicators. We also extend the applicability of RDI to handle general surfaces with boundaries, features, and ridge points, thereby enhancing its versatility and usefulness in various scenarios. To demonstrate the robustness of RDI, we conduct a series of experiments on non-uniform meshes and general surfaces, and compare its performance with some alternative methods. By addressing the challenges posed by the Gibbs phenomena and providing reliable detection of discontinuities, RDI opens up possibilities for improved approximation and analysis of piecewise continuous functions, such as in data remap.


page 22

page 24

page 25

page 27

page 37


On explicit form of the FEM stiffness matrix for the integral fractional Laplacian on non-uniform meshes

We derive exact form of the piecewise-linear finite element stiffness ma...

Understand Slope Limiter – Graphically

In this article, we illustrate how the concept of slope limiter can be i...

WLS-ENO Remap: Superconvergent and Non-Oscillatory Weighted Least Squares Data Transfer on Surfaces

Data remap between non-matching meshes is a critical step in multiphysic...

A robust, high-order implicit shock tracking method for simulation of complex, high-speed flows

High-order implicit shock tracking is a new class of numerical methods t...

Piecewise Padé-Chebyshev Reconstruction of Bivariate Piecewise Smooth Functions

We extend the idea of approximating piecewise smooth univariate function...

Corrected subdivision approximation of piecewise smooth functions

Subdivision schemes are useful mathematical tools for the generation of ...

Please sign up or login with your details

Forgot password? Click here to reset