In this study, numerical simulation of natural gas combustion with highly preheated air in a gas-fired furnace is performed. The predictions by using k-ε turbulence model, Reynolds stress model (RSM), fast chemistry model, and DO radiation model are investigated with comparison to recent measurements. As validated by experimental data, results show that Reynolds stress model in association with DO radiation model and fast chemistry model can perform better simulation than those of other models in out-flame flow-field. On the other hand, k-ε turbulence model in association with DO radiation model and fast chemistry model can perform better simulation in in-flame flow region.