Numerical simulation of shock wave propagation in 2-D channels with obstacles filled with chemically reacting gas suspensions

Abstract

Modification of the serial Fortran code for solving unsteady 2-D Euler equations for the mixture of compressible gas and polydisperse particles was carried out using OpenMP technology. Modified code was verified and parallel speed-up was measured. Analysis showed that the data on parallel efficiency is in a good agreement with the Amdahls law, which gives the estimate for serial code fraction about 30%. Parallel code was used for the numerical simulation of two test-cases, namely shock wave propagation in 2-D channel with obstacles filled with reactive Al-O2 gas particle mixture and heterogeneous detonation propagation in polydisperse suspensions. For the first test-case the data on particles distribution in the flow was obtained, the existense of particle free zones inside the vortices was demonstrated and the attenuation of a shock wave was studied. In the second test, numerical simulation of detonation shock wave propagation in plain 2-D channel for the three polydisperse mixtures was carried out and data on detonation regimes was also obtained.