A Rigorous Examination of Four-Way Coupled Euler-Lagrange Simulations of Shock-Particle Curtain Interaction

The classical planar shock-curtain problem is computationally investigated in an Euler-Lagrange framework. A matrix of simulations is conducted in which each force component the particles are subjected to along with the Reynolds Stresses closures are systematically varied. The quest is to progressively turn on individual force components enabling us to closely monitor spatio-temporal effects on the particle bed. Classical corrections (such as Sangani's volume fraction corrections [3]) seem to overshoot the observed curtain evolution compared to most recent models such as Osnes's [2] and Briney's [1]. The mean quasi-steady is the main force driver, whereas the fluctuating quasi-steady recently published [2] seemed to greatly diffuse the particle interfaces. The unsteady forces have substantial contributions to the bed expansion, especially in the middle of the bed. The quasi-steady effects decrease when these unsteady effects are included. Finally, collisions are observed to play a role in the downstream propagation.