An Assessment of the Drag Models in the Case of a Shock Interacting With a Fixed Bed of Point Particles

In this work, three-dimensional Euler–Lagrange (EL) point-particle simulations of a shock wave interacting with a fixed bed of particles are carried out. The results from the particle-resolved (PR) simulations are used to assess the performance of the point-particle drag models during short time scales. We demonstrate that in a one-way coupled regime, the point-particle simulations recover the dominant gas dynamic features of the flow and are in a good agreement with the exact Riemann solution of a shock traveling through a sudden area contraction. Although the PR simulations are inviscid, we show that a dissipative drag is necessary to predict the mean behavior of the gas. As a model for the inviscid shock-induced (SI) drag two different models are presented in lieu of the quasi-steady drag. Finally, two-way coupled simulations are performed at four different particle volume fractions {0.10, 0.15, 0.20, 0.25} and three different incident shock Mach numbers {1.22, 1.66, 3.0} and compared against the data from PR inviscid simulations. At a lower Mach number (1.22), averaged flow quantities from the two-way coupled simulations agree well with the PR simulations. As the Mach number increases, we observe that the discrepancies between the point-particle and the PR simulations grow. A sensitivity analysis of the drag models involved reveals a strong influence of the inviscid-unsteady force on the gas quantities especially in the case of a strong shock interacting with a dense bed of particles. The use of Mach correlation beyond the subcritical regime coupled with the model for volume fraction correction is identified as a probable cause for the additional drag.