This work is devoted to the development and validation of a subgrid model describing heat transfer between the bulk flow and moving particles undergoing phase change under the influence of forced and free convection. Such kinds of submodels play the role of "scale bridges" between microscale (e.g. interfacial heat transfer) and macroscale (e.g. bulk flow) phenomena. Applied to multiscale modelling of particulate flows with phase change phenomena, our model serves as a coupling between equations describing particle movement in Lagrangian space and the mass and heat conservation equation defining melt flow in Eulerian space. The input parameters in our model are the particulate Reynolds number (Re), the Grashof number (Gr), the Stefan number (Ste), and the Prandlt number (Pr). The model has been validated against experimental data published recently in the literature applied to the melting of ice spheres under different flow conditions. Good agreement between our model predictions and published experimental data (Shukla et al.([16]) and Hao and Tao([14])) is observed.