The coking behaviour and NOx emission characteristics of the coke oven are of great significance for coke production. Based on the realisation k-epsilon turbulence model, the P-1 radiation model, the species transport model with eddy dissipation, and the thermal NOx model. A comprehensive model for traditional regenerative coke oven was established with reversing mode and external flue gas recirculation model considering temperature-dependent thermophysical properties. The convective term of model equations was the discretised using the second-order upwind scheme, while the others were applied by the central difference scheme. The coking behaviour of the coke oven with external flue gas recirculation was numerically solved through the SIMPLE algorithm, and the effects of reversing period and external flue gas recirculation ratio on the coking behaviour were explored. When increasing the reversing period, the mean temperature of the vertical flue increases, and the coking time decreases, but the temperature uniformity of the carbonisation chamber becomes worse. As the external flue gas recirculation ratio increases, the average temperature of the vertical flue decreases, and the NOx concentration at the outlet decreases, yet the coking time increases. These results would provide useful information for guiding the improvement of coke oven operation and reducing NOx emissions.