A mathematical model has been developed to describe the various processes occurring in a flash furnace shaft. The model incorporates turbulent fluid dynamics, chemical reaction kinetics, and heat and mass transfer. The key features include the use of the k-epsilon turbulence model, incorporating the effect of particles on the turbulence, and the four-flux model for radiative heat transfer. The model predictions were compared with measurements obtained in a laboratory flash furnace and a pilot plant flash furnace. Good agreement was obtained between the predicted and measured data in terms of the SO2 and O2 concentrations, the amount of sulfur remaining in the particles, and the gas temperature. Model predictions show that the reactions of sulfide particles are mostly completed within about 1 m of the burner, and the double-entry burner system with radial feeding of the concentrate particles gives better performance than the single-entry burner system. The model thus verified was used to further predict various aspects of industrial flash furnace operation. The results indicate that from the viewpoint of sulfide oxidation, smelting rate can be substantially increased in most existing industrial flash furnaces.