The discrete element method (DEM) is used to simulate two types of mixing processes (radial mixing and axial mixing) of particles with and without cohesive interparticle forces in a two-dimensional (2D) fluidized bed. By using the Ashton index, the mixing degree is quantified and the differences between the mixing characteristics of particles with and without cohesive interparticle forces are compared. The cohesive interparticle forces, which impair bubble behaviors, can apparently weaken the mixing process and prolong the mixing time. This effect increases with increasing cohesive interparticle force. Radial mixing, which is mainly dependent on the convective mixing caused by bubbles, is relatively sensitive to cohesive interparticle force and shows a long mixing time, while axial mixing, which is mainly induced by the bubble wake, is less sensitive to cohesive interparticle force and presents a fast mixing process. Finally, to weaken the effect of cohesive interparticle force, an uneven air distribution is introduced, which greatly improves the mixing process by changing the bubble movement. This method may be helpful for the reduction of particle agglomeration in high-temperature fluidized beds. (C) 2017 Elsevier B.V. All rights reserved.