Improved well-balanced AWENO schemes with hydrostatic reconstruction for the Euler equations under gravitational fields

Abstract

The Euler equations under gravitational fields allow the hydrostatic equilibrium states, which requires that the numerical scheme of the system should also have this characteristic. In our previous work, a well-balanced finite difference conservative AWENO scheme has been constructed to preserve the isothermal equilibrium state accurately (Fu et al., Appl. Numer. Math., 180:1-15, 2022). However, the scheme fails to maintain more complicated hydrostatic equilibrium states, such as the isentropic equilibrium state. In this study, the improved well-balanced AWENO schemes with hydrostatic reconstruction are designed to maintain the general equilibrium states numerically. The key technologies to achieve the well-balanced property are rewriting the source term into an appropriate form based on the equilibrium equation, acting on the conserved variables with an affine-invariant WENO or linearized WENO interpolation operator, employing hydrostatic reconstruction methods to create auxiliary conserved variables, and modifying the numerical flux with auxiliary conserved variables and higher-order conserved terms based on pressure. The proposed well-balanced schemes maintain the equilibrium solution near the machine precision and capture small perturbations of the equilibrium state without spurious oscillations even at the coarse mesh. Furthermore, they can be applicable to the general hydrostatic equilibrium more than the isothermal equilibrium state. Several numerical experiments are tested to show that the developed well-balanced schemes perform well in maintaining equilibrium states and capturing small perturbations.