In this study, a constitutive model that is applicable to both dilute and dense regimes for simulating gas-particle flows with polydisperse particles is proposed. In dilute regime, solid stress is closed by kinetic theory of granular flow (KTGF). In dense regime, solid stress is closed by the η(Is) rheology. The transition from dilute to dense regimes is realized smoothly by the dynamic blending function χmix, which is a function of the inertial number of mixture Imix. The dynamic blending function χmix is also used to link the solid–solid drag force at different regimes. This new model was validated with the experimental data of polydisperse bubbling fluidized bed of Girimonte et al. (2018). When compared with the conventional baseline model, this new model improves the transition from dilute to dense regimes and the predictions in both density and size segregations.