The present study investigates the hydrodynamics and heat transfer behavior of municipal solid waste (MSW) gasification in a pilot scale bubbling fluidized bed reactor. A multiphase 2-D numerical model following an Eulerian-Eulerian approach within the FLUENT framework was implemented. User defined functions (UDFs) were coupled to improve hydrodynamics and heat transfer phenomena, and to minimize deviations between the experimental and numerical results. A grid independence study was accomplished through comparison of the bed volume fraction profiles and by reasoning the grid accuracy and computational cost. The standard deviation concept was used to determine the mixing quality indexes. Simulated results showed that UDFs improvements increased the accuracy of the mathematical model. Smaller size ratio of the MSW-dolomite mixture revealed a more uniform mixing, and larger ratios enhanced segregation. Also, increased superficial gas velocity promoted the solid particles mixing. Heat transfer within the fluidized bed showed strong dependence on the MSW solid particles sizes, with smaller particles revealing a more effective process.