The eXtended Discrete Element Method (XDEM) platform which is a Coupled Eulerian-Lagrangian framework with heat and mass transfer, is extended for melting of granular packed beds. In this method, the fluid is simulated by computational fluid dynamics (CFD) and the soft-sphere discrete element approach (DEM) is used for the particle system. A four-way coupling accounts for solid-liquid interaction via drag and buoyancy forces and the collisions between the particles and the walls. The contact forces between the particles and wall-particle contacts have been calculated by the hertz-mindlin model. The particles heat up, melt and shrink due to heat and mass exchange, and the temperature distributions inside the particles are described. In order to validate the method, melting of a single ice particle and of a packed bed of ice in flowing water have been carried out. Very good agreement between the simulation and experiment has been achieved. The effects of the temperature and velocity of flowing water on melting rate are also discussed. (C) 2017 Elsevier Ltd. All rights reserved.