A numerical approach is presented for predicting the species concentrations and temperature profiles of chemically reacting flow in the non-catalytic partial oxidation of hot coke oven gas (HCOG) in a pilotscale reformer installed on an operating coke oven. A detailed chemical kinetic model consisting of 2216 reactions with 257 species ranging in size from the hydrogen radical to coronene was used to predict the chemistries of HCOG reforming and was coupled with a plug model and one-dimensional (1D) flow with axial diffusion model. The HCOG was a multi-component gas mixture derived from coal dry distillation, and was approximated with more than 40 compounds: H-2, CO, CO2, CH4, C-2 hydrocarbons, H2O, aromatic hydrocarbons such as benzene and toluene, and polycyclic aromatic hydrocarbons up to coronene. The measured gas temperature profiles were reproduced successfully by solving the energy balance equation accounting for the heat change induced by chemical reactions and heat losses to the surroundings. The approach was evaluated critically by comparing the computed results with experimental data for exit products such as H-2, CO, CO2, and CH4, in addition to the total exit gas flow rate. The axial diffusion model slightly improves the predictions of H-2, CO, and CO2, but significantly improves those of CH4 and total exit flow rate. The improvements in the model predictions were due primarily to the improved temperature predictions by accounting for axial diffusion in the flow model. (C) 2014 Elsevier B.V. All rights reserved.