In this work, a numerical model is used to study the flow and coal combustion along the coal plume in a large-scale setting simulating the lance-blowpipe-tuyere-raceway region of a blast furnace. The model formulation is validated against the measurements in terms of burnout for both low and high volatile coals. The typical phenomena related to coal combustion along the coal plume are simulated and analyzed. The effects of some operational parameters on combustion behavior are also investigated. The results indicate that oxygen as a cooling gas gives a higher coal burnout than methane and air. The underlying mechanism of coal combustion is explored. It is shown that under the conditions examined, coal burnout strongly depends on the availability of oxygen and residence time. Moreover, the influences of two related issues, i.e. the treatment of volatile matter (VM) and geometric setting in modeling, are investigated. The results show that the predictions of final burnouts using three different VM treatments are just slightly different, but all comparable to the measurements. However, the influence of the geometric setting is not negligible when numerically examining the combustion of pulverized coal under blast furnace conditions.