Summary
A novel approach for selecting appropriate reconstructions is implemented to the hyperbolic conservation laws in the high‐order local polynomial‐based framework, for example, the discontinuous Galerkin (DG) and flux reconstruction (FR) schemes. The high‐order polynomial approximation generally fails to correctly capture a strong discontinuity inside a cell due to the Gibbs phenomenon, which is replaced by more stable approximation on the basis of a troubled‐cell indicator such as that used with the total variation bounded (TVB) limiter. This paper examines the applicability of a new algorithm, so‐called boundary variation diminishing (BVD) reconstruction, to the weighted essentially nonoscillatory methodology in the FR framework including the nodal type DG method. The BVD reconstruction adaptively chooses a proper approximation for the solution function so as to minimize the jump between values at the left and right side of cell boundaries. The results of the BVD algorithm are comparable to those with the conventional TVB limiter in terms of oscillation suppression and numerical dissipation in one‐dimensional linear advection and nonlinear system equations, while the TVB limiter performs better in the case with strong discontinuities (the blast wave problem in the Euler equations). Overall, since the present BVD algorithm does not need any ad hoc constant such as the TVB parameter, it could be more reliable than the conventional TVB limiter that is often used in the DG and FR communities for shocks and other discontinuities. The proposed method would lead to a parameter‐free robust algorithm for the local polynomial‐based high‐order schemes.