Abstract:
A linearized version of the classical Godunov scheme as applied to nonlinear discontinuity decays is described. It is experimentally shown that this version guarantees an entropy nondecrease, which makes it possible to simulate entropy growth on shock waves. The structure of shock waves after the discontinuity decays is studied. It is shown that the width of the shock waves and the time required for their formation depend on the choice of the Courant number. The accuracy of the discontinuous solutions is tested numerically.