Mesoscale bubbles exist inherently in bubbling fluidized beds and hence should be considered in the constitutive modeling of the drag force. The energy minimization multiscale bubbling (EMMS/bubbling) drag model takes the effects of mesoscale structures (

i.e

.

bubbles) into the modeling of drag coefficient and thus improves the coarse-grid simulation of bubbling and turbulent fluidized beds. However

its dependence on the bubble diameter correlation has not been thoroughly investigated. The hydrodynamic disparity between homogeneous and heterogeneous fluidization is accounted for by the heterogeneity index

H

d

which can be affected by choice of bubble diameter correlation. How this choice of bubble diameter correlation influences the model prediction calls for further fundamental research. This article incorporated seven different bubble diameter correlations into EMMS/bubbling drag model and studied their effects on

H

d

. The performance of these correlations has been compared with the correlation used previously by EMMS/bubbling drag model. We found that some of the correlations predicted lower

H

d

by order of a magnitude than the correlation used by th

e original EMMS/bubbling drag. Based on such analysis

we proposed a modification in the EMMS drag model for bubbling and turbulent fluidized beds. A computational fluid dynamics (CFD) simulation using two-fluid model with the modified EMMS/bubbling drag model was performed for two bubbling and one turbulent fluidized beds. Voidage distribution

time averaged solid concentration and axial solid concentration profiles were studied and compared with the previous version of the EMMS/bubbling drag model and experimental data. We found that the right choice of bubble diameter correlations can significantly improve the results for CFD simulations.