Knowledge on regime transition from conventional to circulating fluidization is important to scale-up and operation of liquid-solid fluidized systems in industrial applications. Previous experimental studies reported that the onset velocity measured by bed empty time test is pertinent to the physical properties of particle and liquid, and independent of solid inventory, bed geometry, and configuration. The effect of liquid properties is, however, not clear yet since tap water is usually used as the working fluid in laboratory test for convenience. To address this problem, a Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) model is applied to numerically measure the onset velocity of particles in different liquid media. The model is first validated in terms of the experimental data of bed expansion in literature. Then several kinds of particles and three kinds of liquid media of different densities and viscosities are further simulated. We find that the Reynolds number based on onset velocity and Archimedes number follows a power-law relationship. With the decrease of Archimedes number, the discrepancy of Reynolds number based on onset velocity and that on particle terminal velocity becomes highly significant. The ratio of the onset velocity to particle terminal velocity for various particles in one liquid medium is not a constant, and instead dependent on both the particle and liquid properties. (C) 2020 Elsevier B.V. All rights reserved.