Simulation of the layer inversion phenomenon in binary liquid--fluidized beds by DEM–CFD with a drag law for polydisperse systems

Liquid-fluidized beds of particular binary mixtures exhibit the layer inversion phenomenon, a peculiar result of the mechanical equilibrium developing in such multiphase systems. Because of the crucial role of the hydrodynamic interaction, these are ideal test cases for assessing fluid-particle models in multi-particle CFD simulations. In the present work the layer inversion phenomenon is reproduced via Discrete Element Method (DEM) simulations, in which the drag force model for polydisperse systems described in Cello et al. (2010) is used. The simulations serve primarily to assess the suitability of the DEM–CFD approach and particularly the validity of the above drag model, although the results available may also prove very useful to investigate the prevailing mechanisms. To analyze a sufficiently broad variety of cases, involving different solid species, size ratios and operating conditions, three systems selected amongst the published literature are considered. The comparison of simulation and experimental observation is carried out in terms of overall bed and mixed layer interface heights, voidage and component distributions along bed height and time evolution of the species centers of mass. Simulations with liquid velocities below, above and at the critical inversion conditions are carried out. Additionally, simulations using a drag force model traditionally used for monodisperse beds are reported for comparison. The results not only demonstrate the importance of correctly accounting for the local size distribution in the bed, but also prove the validity of the overall computational approach. The predictions of the simulations are in good to excellent agreement with experiments, depending on the system considered, both in terms of critical velocity and, most notably, expansion of the individual components in the bed. The analysis of the hydrodynamics in the bed allows to investigate the local particle flow field, highlighting the presence of a steady irregular motion of the solids in apparently chaotic vortices continuously forming and disappearing, which is thought to be the mechanism responsible for mixing. The fluid–particle interaction forces exhibit a constant profile along bed height, even in the presence of a strongly non-uniform concentration profile.