Local linear stability of dual-pairing summation-by-parts methods for nonlinear conservation laws
This addresses a critical issue in computational fluid dynamics for researchers and engineers, providing a novel strategy to improve simulation reliability, though it builds incrementally on prior entropy-stable methods.
The paper tackles the problem of designing high-order numerical methods for nonlinear conservation laws that are both entropy-stable and locally energy-stable, showing that entropy-stable volume upwind filters in dual-pairing summation-by-parts methods can ensure local energy stability, with numerical experiments validating this in 1D and 2D simulations including turbulent scales.
A recent study by Gassner et al. [J. Sci. Comput. 90:79 (2022)] demonstrates that local energy stability--that is, ensuring the asymptotic numerical growth rate does not exceed the continuous growth rate--is crucial for achieving accurate numerical simulations of nonlinear conservation laws. While nonlinear entropy stability is necessary for numerical stability (i.e., ensuring the boundedness of nonlinear numerical solutions), local energy stability is essential to prevent unresolved high-frequency wave modes from dominating the simulation. Currently, it remains an open question whether high-order numerical methods for nonlinear conservation laws can be simultaneously entropy-stable and locally energy-stable. In this work, we examine the local energy-stability properties of recently developed entropy-stable, high-order accurate dual-pairing (DP) SBP methods, as introduced by Duru et al. [arXiv: 2411.06629], for nonlinear conservation laws. Our analysis indicates that the entropy-stable volume upwind filter inherent in these methods can ensure local energy stability. This approach offers a novel numerical strategy for designing reliable high-order methods for nonlinear conservation laws that are provably entropy-stable and locally energy-stable. The theoretical findings are supported by numerical experiments involving the inviscid Burgers equation and nonlinear shallow water equations, in 1D and 2D. Furthermore, we present accurate numerical simulations of 2D barotropic shear instability, with fully developed turbulence, demonstrating the efficiency of the DP SBP method in resolving turbulent scales.