Truncation Error Estimation in the p-Anisotropic Discontinuous Galerkin Spectral Element Method
Provides a more efficient and accurate error estimation method for p-adaptation in high-order CFD simulations, benefiting computational scientists.
The paper presents a novel anisotropic truncation error estimator for p-adaptation in DGSEM, derived from τ-estimation. The method reduces computational cost and improves accuracy for higher-order extrapolations, validated on compressible Navier-Stokes equations.
In the context of Discontinuous Galerkin Spectral Element Methods (DGSEM), $τ$-estimation has been successfully used for p-adaptation algorithms. This method estimates the truncation error of representations with different polynomial orders using the solution on a reference mesh of relatively high order. In this paper, we present a novel anisotropic truncation error estimator derived from the $τ$-estimation procedure for DGSEM. We exploit the tensor product basis properties of the numerical solution to design a method where the total truncation error is calculated as a sum of its directional components. We show that the new error estimator is cheaper to evaluate than previous implementations of the $τ$-estimation procedure and that it obtains more accurate extrapolations of the truncation error for representations of a higher order than the reference mesh. The robustness of the method allows performing the p-adaptation strategy with coarser reference solutions, thus further reducing the computational cost. The proposed estimator is validated using the method of manufactured solutions in a test case for the compressible Navier-Stokes equations.