A three-dimensional mathematical model of the combustion of pulverized coal and coke is developed. The model is applied to the region of lance-blowpipe-tuyere-raceway-coke bed to simulate in-furnace phenomena of pulverized coal injection in an ironmaking blast furnace. The model integrates not only pulverized coal combustion model in the blowpipe-tuyere-raceway-coke bed but also coke combustion model in the coke bed. The model is validated against the measurements under different conditions. The comprehensive in-furnace phenomena are investigated in the raceway and coke bed, in terms of flow, temperature, gas composition, and coal burning characteristics. The underlying mechanisms for the in-furnace phenomena are also analysed. The simulation results indicate that it is important to include recirculation region in the raceway and the coke bed reactions for better understanding in-furnace phenomena. The model provides a cost-effective tool for understanding and optimizing the in-furnace flow-thermo-chemical characteristics of the PCI operation in full-scale blast furnaces. (C) 2010 Elsevier Ltd. All rights reserved.