Recently, most research on gas-liquid flow through packed particulate beds has focused on macro scale parameters, such as the bed pressure drop, liquid holdup, and soon. However, very few efforts have focused on the local flow behavior, which is essential for understanding the concept of flow mal-distribution. Therefore, in this study, the CFD (Computational Fluid Dynamics) - VOF (Volume-of-Fluid) method, which is able to track the position of the interface between multiphase fluids and the effects of surface tension, was employed to investigate the gas liquid two-phase flow across a spherical particle. By observing the whole process of the liquid flow pattern as well as vortex production, we found that the gas vortex was strongly related to the liquid flow pattern, whose shape and scale vary with the liquid flow stages. Moreover, parametric studies were conducted to determine how the vortices beneath and on the surface of the particle change. The fully revealed mechanism of vortex formation on the bottom of a particle could significantly benefit the understanding of gas-liquid flow hysteresis inside packed beds and thus be applied to improve the reactor performance. (C) 2016 Elsevier B.V. All rights reserved.