“…An improper treatment of the gravitational force and thus the steady states may lead to a solution, which either oscillates around the equilibrium or deviates from the equilibrium in time. This problem is well known in the numerical community and received a great attention in the literature; see, e.g., [1,7,9,17,36,37,39,55,56,60,61]. For the compressible Euler equations that describe the gas phase, we use a well-balanced version of the second-order semi-discrete CU scheme introduced in [9].…”