Abstract
We study the interaction of a planar air shock with a perturbed, monodispersed, particle curtain using point-particle simulations. In this Eulerian-Lagrangian approach, equations of motion are solved to track the position, momentum, and energy of the computational particles while the carrier fluid flow is computed in the Eulerian frame of reference. In contrast with many Shock-Driven Multiphase Instability (SDMI) studies, we investigate a configuration with an initially high particle volume fraction, which produces a strongly two-way coupled flow in the early moments following the shock-solid phase interaction. In the present study, the curtain is about 4 mm in thickness and has a peak volume fraction of about 26%. It is composed of spherical particles of d = 115μm in diameter and a density of 2500 kg.m−3, thus replicating glass particles commonly used in multiphase shock tube experiments or multiphase explosive experiments. We characterize both the evolution of the perturbed particle curtain and the gas initially trapped inside the particle curtain in our planar three-dimensional numerical shock tube. Control parameters such as the shock strength, the particle curtain perturbation wavelength and particle volume fraction peak-to-trough amplitude are varied to quantify their influence on the evolution of the particle cloud and the initially trapped gas. We also analyze the vortical motion in the flow field. Our results indicate that the shock strength is the primary contributor to the cloud particle width. Also, a classic Richtmyer-Meshkov instability mixes the gas initially trapped in the particle curtain and the surrounding gas. Finally, we observe that the particle cloud contribute to the formation of longitudinal vortices in the downstream flow.