A numerical model based on the kinetic theory of granular media has been tested to simulate the unsteady two-phase flow in solid–liquid fluidized beds. We focus here on the local unsteady flow features as predicted by the model and comparison with experimental data of Zenit and Hunt (2000) are discussed. In most of the cases, a series of two-dimensional rotational structures develops along the bed height, which controls the instantaneous solid fraction field as well as the level of fluctuating kinetic energy of the particles. Interestingly, with high inertia particles and at high solid concentration, these structures disappear to the benefit of a one-dimensional oscillating flow. The shape of these concentration waves is very similar to tho...