The high temperature air combustion performance of natural gas in an industrial furnace with a swirling burner was numerically modeled. A Beta function PDF (Probability Density Function) combustion model was selected to simulate the gas combustion combined with the Reynolds Stress Model (RSM) to simulate the turbulent flow. The radiation was simulated by a Discrete Ordinates method. The NO chemistry was simulated by thermal NO model. The simulation was performed at inlet air oxygen fraction 8% and the total air excess ratio 1.1 for natural gas. The effect of preheated air temperature on NO emission, temperature, O2 and CO distribution in the furnace was investigated. Results showed that thermal NO emission increased when the preheated air temperature increased from 1073 K to 1473K. When the preheated air temperature increased, both of the maximum and averaged temperature in the furnace increased. The oxygen was consumed by the formation of thermal NO at higher inlet air temperature and the fuel was not fully burnt out.