Numerical simulations of a gas‐particle‐droplet system were performed using an Euler‐Lagrange approach. Models accounting for (1) the interaction between droplets and particles, (2) evaporation from the droplet spray, as well as (3) evaporation of liquid from the surface of non‐porous particles were considered. The implemented models were verified for a packed bed, as well as other standard flow configurations. The developed models were then applied for the simulation of flow, as well as heat and mass transfer in a fluidized bed with droplet injection. The relative importance of droplet evaporation vs. evaporation from the particle surface was quantified. It was proved that spray evaporation competes with droplet deposition and evaporation from the particle surface. Moreover, we show that adopting a suitable surface coverage model is vital when attempting to make accurate predictions of the particle's liquid content. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2569–2587, 2017

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.