ATTENUATION OF SPURIOUS OSCILLATIONS IN THE NUMERICAL CAPTURE OF SHOCK WAVES

Abstract

Artificial viscosity is often expressed as a superposition of linear and quadratic terms in the first derivative of the velocity field. In trying to find a continuous solution for the hydrodynamic equations, we propose an alternative one-term artificial viscosity which is a linear form of the derivative of the specific volume. It is shown that this artificial viscosity is able to capture the profile of the steady plane shock wave, largely removing the non-physical oscillations originated by the artificial viscosity of von Neumann and Richtmyer. Analytical and numerical calculations for one-dimensional shock using both artificial viscosities are compared.