2021/06/03 by Hiro Wakimura, Wakimura, Hiro, Shinichi Takagi +3 · 1 citation
Computer Science · Mathematics · Physics and Astronomy · #Computational Physics (physics.comp-ph) #FOS: Mathematics #FOS: Physical sciences #Numerical Analysis (math.NA) #cs.NA #math.NA #physics.comp-ph
paper · pdf · doi:10.48550/arxiv.2106.01872
32 pages, 8 figures
arxiv created 2021/06/03 · arxiv updated 2021/06/04
A class of high-order shock-capturing schemes, PnTm-BVD (Deng et al., J. Comp. Phys., 386:323-349, 2019; Comput. & Fluids, 200:104433, 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.