A highly resolved large eddy simulation (LES) of the semi-industrial IFRF coal furnace [1,2] employing the steady flamelet model is presented. The flamelet table is based on mixture fractions of volatile and char off-gases as well as on enthalpy and scalar dissipation rate. Turbulence–chemistry interaction is treated with an assumed pdf approach, with the variance obtained from a transport equation. Radiation is computed by the discrete ordinates method and the grey weighted sum of grey gases model. The simulation is conducted with the massively parallel “PsiPhi” code on up to 1.7 billion cells and with 40 million particles. Results are processed and compared against the comprehensive set of experiments to (i) validate the new flamelet model and the simulation method and to (ii) gain further insight into the combustion process that is not available from the experiment. The simulation results show that the flamelet LES approach can successfully describe the flow field and combustion inside the furnace; major species and velocities are found in good agreement with the experiment. The results are further analyzed with a focus on the processes of particle heating, devolatilization, char combustion and flame stabilization in a highly turbulent environment. Additionally, the relative importance of scalar dissipation rate is highlighted, showing a large separation of mixing scales between volatile and char off-gas combustion due to the long residence time and generally much lower scalar dissipation rates than typical for lab-scale experiments.