Symmetry-preserving enforcement of low-dissipation method based on boundary variation diminishing principle

A class of high-order shock-capturing schemes, PnTm-BVD (Deng et al., 2019; Deng et al., 2020) schemes, have been devised to solve the Euler equations with substantially reduced numerical dissipation, which enable high-resolution simulations to resolve flow structures of wider range scales. In such simulations with low dissipation, errors of round-off level might grow and contaminate the numerical solutions. A typical example of such problems is the loss of symmetry in the numerical solutions for physical problems of symmetric configurations even if the schemes are mathematically in line with the symmetry rules. In this study, the mechanisms of symmetry-breaking in a finite volume framework with the P4T2-BVD reconstruction scheme are thoroughly examined. Particular attention has been paid to remove the possible causes due to the lack of associativity in floating-point arithmetic which is associated with round-off errors. Modifications and new techniques are proposed to completely remove the possible causes for symmetry breaking in different components of the P4T2-BVD finite volume solver. Benchmark tests that have symmetric solution structures are used to verify the proposed methods. The numerical results demonstrate the perfect symmetric solution structures.

Automated summary

A class of high-order shock-capturing schemes, PnTm-BVD, has been developed to solve Euler equations with reduced numerical dissipation. This study thoroughly examines the mechanisms of symmetry-breaking in a finite volume framework with the P4T2-BVD reconstruction scheme and proposes modifications to completely remove the causes for symmetry breaking. Benchmark tests confirm perfect symmetric solution structures.