Simulation of astrophysical flows in isentropic approximation
One of the difficulties of numerical simulations of cold supersonic astrophysical flows is a big difference in different types of energy. Gravitational and/or kinetic energy of the gas could be much larger than its internal energy. In such a case, it is possible to get large numerical errors in the simulations. To avoid this difficulty, conservation of entropy equation was used instead of energy conservation equation. The entropy conservation equation does not contain the gravitational and kinetic energy. The application of the isentropic set of equations is correct when the flow does not contain shocks or the amplitude of the shocks (shock wave Mach number) is not large. We estimate the violation of the energy conservation low when the “shock wave” is isentropic.