Modeling particle collisions in moderately dense curtain impacted by an incident shock wave

Abstract

The interactions between an incident shock and a moderately dense particle curtain are simulated with the Eulerian–Lagrangian method. A customized solver based on OpenFOAM is extended with an improved drag model and collision model and then validated against two benchmark experiments. The results show that the collision model has a limited impact on curtain morphology compared with the improved drag model. In this work, parametric studies are performed considering different particle sizes, volume fractions, and curtain thicknesses. Smaller particle sizes and larger volume fractions lead to stronger reflected shock and weaker transmitted shock. Attention is paid to the particle collision effects on the curtain evolution behaviors. According to our results, for the mono-dispersed particle curtain, the collision effects on curtain front behaviors are small, even when the initial particle volume fraction is as high as 20%. This is due to the positive velocity gradient across the curtain after the shock wave passage, leading to the faster motion of downstream particles than the upstream ones, and hence, no collision occurs. For the bi-dispersed particle curtain, the collision effects become important in the mixing region of different-size particles. Collisions decelerate small particles while accelerating large ones and cause velocity scattering. Moreover, increasing the bi-dispersed curtain thickness leads to multiple collision force peaks, which is the result of the delayed separation of different particle groups. Our results indicate that the collision model may be unnecessary to predict curtain fronts in mono-dispersed particles, but in bi-dispersed particles, the collision effects are important and, therefore, must be modeled.