The three-dimensional code ESTET developed at the LNH has been used to predict the reactive flow in a 600 MW coal fired boiler. Assuming a no-slip condition between the gas and the coal, the equations for a gas-particle mixture can be written. The pulverized coal particle size distribution is represented by a discrete number of particle size groups determined by the measured fineness distribution. The combustion models taking into account the pyrolysis of the particle and the heterogeneous combustion of char have been validated using intensive measurements performed on the 600 MW utility boiler. Heat fluxes were measured along the walls of the furnace and satisfactory agreement between computation and measurements has been achieved in terms of maximum flux location and heat flux intensity. Local measurements of velocities using LDV probe, gas temperature and gas species concentrations were performed in the vicinity of one burner and compared with the computed variables. Again we have observed a good agreement between the computations and the measurements in terms of jet penetration, temperature distribution, oxygen concentration and ash content. The three-dimensional code ESTET developed at the LNH has been used to predict the reactive flow in a 600 MW coal fired boiler. Assuming a no-slip condition between the gas and the coal, the equations for a gas-particle mixture can be written The pulverized coal particle size distribution is represented by a discrete number of particle size groups determined by the measured fineness distribution. The combustion models taking into account the pyrolysis of the particle and the heterogeneous combustion of char have been validated using intensive measurements performed on the 600 MW utility boiler. Heat fluxes were measured along the walls of the furnace and satisfactory agreement between computation and measurements has been achieved in terms of maximum flux location and heat flux intensity. Local measurements of velocities using LDV probe, gas temperature and gas species concentrations were performed in the vicinity of one burner and compared with the computed variables. Again we have o bserved a good agreement between the computations and the measurements in terms of jet penetration, temperature distribution, oxygen concentration and ash content.