The discrete element method-large eddy simulation (DEM-LES) is used to model coal combustion at the particle level in a bubbling fluidized bed. The gas phase is modelled as a continuum and the solid phase is modeled by DEM. Chemical reactions consist in the heterogeneous reactions of char with O-2, CO, CO2, N2O, and N2O, and in the homogeneous reactions involving CO, O-2, NO, and N2O. The colliding particle-particle heat transfer is based on the analysis of the elastic deformation of the spheres during their contact. The model predicts the effects of the particle heterogeneous flow structure on the thermal characteristics of coal particles when heating and burning, and the gaseous emissions from a fluidized sand-coal binary mixture. The heating rates are 1627 and 432 K/s for, respectively, 0.8 and 2.0 mm diameter coal particles fed into the fluidized bed. The instantaneous contribution of the collision heat transfer is weak, less than 5.0% of the total power exchanges (coal combustion, radiation, convection and collision) during the heating and 1.5% during the combustion. The temperature of the coal particles exceeds the bed temperature, which is in qualitative agreement with experimental data from literature. The effects of the diameter of coal particles, of the bed temperature, and of the inlet gas velocity on the thermal characteristics are also studied. (C) 2004 Elsevier Ltd. All rights reserved.