The fluidization behavior of cohesive particles was investigated using an Euler–Euler approach. To do so, a two-fluid model (TFM) platform was developed to account for the cohesivity of particles. Specifically, the kinetic theory of granular flow (KTGF) was modified based on the solid rheology developed by Gu et al. J. Fluid Mech.2019. The results of our simulations demonstrated that the modified TFM approach can successfully predict the formation of particle agglomerates and clusters in the fluidized bed, induced by the negative (tensile-dominant) pressure. The formation of such granules and clusters highly depended on the particle Bond number and the tensile pressure prefactor. To evaluate fluidization regimes, a set of simulations was conducted for a wide range of particle cohesivity (e.g., Bond number and tensile pressure prefactor) at two different fluidization numbers of 2 and 5. Our simulation results reveal the formation of four different regimes of fluidization for cohesive particles: (i) bubbling, (ii) bubbling–clustering, (iii) bubble-less fluidization, and (iv) stagnant bed. Comprehensive analysis of the shear-to-yield ratio reveals that the observed regime map is attributed to the competition between the shear stress and yield stress acting on the particles. The obtained regime map can be extended to incorporate the effect of dimensionless velocity and dimensionless diameter as a comprehensive fluidization chart for cohesive particles. Such fluidization charts can facilitate the design of fluidized beds by predicting the conditions under which the formation of particle agglomeration and clustering is likely in fluidized beds.