In this paper we write, analyze and experimentally compare three different numerical schemes dedicated to the one dimensional barotropic Navier-Stokes equations:
a staggered scheme based on the Rusanov one for the inviscid (Euler) system,a staggered pseudo-Lagrangian scheme in which the mesh “follows” the fluid,the Eulerian projection (on a fixed mesh) of the preceding scheme.
All these schemes only involve the resolution of linear systems (all the nonlinear terms are solved in an explicit way). We propose numerical illustrations of their behaviors on particular solutions in which the density has discontinuities (hereafter called Hoff solutions). We show that the three schemes seem to converge to the same solutions, and we compare the evolution of the amplitude of the discontinuity of the numerical solution (with the pseudo-Lagrangian scheme) with the one predicted by Hoff and observe a good agreement.
Accuracy of the conventional Lagrangian scheme for one-dimensional hydrodynamics