Numerical simulation of a normally incident shock wave – dense particles layer interaction using the Godunov solver for the Baer–Nunziato equations

Abstract

The work is devoted to the numerical simulation of the known problem of a normally incident shock wave – dense layer of particles interaction and the phenomenon of the pressure rise on the wall under the layer. The novelty of the work is in the numerical approach which is based on the Godunov solver for the Baer-Nunziato equations and the pressure relaxation procedure which takes into account intergranular stresses in the solid phase. The algorithm based on the exact solution of the Riemann problem provides a low numerical dissipation of the solid contacts and is robust at the explicit interfacial boundaries. The algorithm was described in detail; the source code of the Godunov solver for the Baer-Nunziato equations was provided. The full scale experiment of a shock wave – particles layer interaction was simulated. The shape of the pressure curve, obtained on the wall under the particles layer, was explained from the point of view of ongoing wave processes in the layer. A quantitative comparison of the experimental and simulated pressure curves was carried out. Studies of the influence of parameters in the intergranular stresses model on simulation results as well as reversible or irreversible character of loading-unloading process were conducted. Obtained results were compared to the published simulation results by the other authors based on the R.I. Nigmatulin models.