Well-Balanced Schemes for the Euler Equations with Gravitation: Conservative Formulation Using Global Fluxes
For computational fluid dynamics researchers, this provides a new well-balanced scheme for Euler equations with gravity, though it is an incremental improvement over existing well-balanced methods.
The authors develop a second-order well-balanced central-upwind scheme for the compressible Euler equations with gravity, using a conservative reformulation with global fluxes. The scheme exactly preserves steady-state solutions and reduces numerical viscosity near steady states, demonstrated in 1D and 2D tests.
We develop a second-order well-balanced central-upwind scheme for the compressible Euler equations with gravitational source term. Here, we advocate a new paradigm based on a purely conservative reformulation of the equations using global fluxes. The proposed scheme is capable of exactly preserving steady-state solutions expressed in terms of a nonlocal equilibrium variable. A crucial step in the construction of the second-order scheme is a well-balanced piecewise linear reconstruction of equilibrium variables combined with a well-balanced central-upwind evolution in time, which is adapted to reduce the amount of numerical viscosity when the flow is at (near) steady-state regime. We show the performance of our newly developed central-upwind scheme and demonstrate importance of perfect balance between the fluxes and gravitational forces in a series of one- and two-dimensional examples.