A thermodynamic model has been developed to predict the distribution behavior of Ni, Cu, Co, Fe, S, As, Sb, and Bi in the Outokumpu flash-smelting process, the Outokumpu direct high-grade matte smelting process, and the INCO flash-smelting process. In this model, as many as 16 elements (Ni, Cu, Co, Fe, As, Sb, Bi, S, O, Al, Ca, Mg, Si, N, C, and H) are considered, and two nickel sulfide species are used to allow for modeling of sulfur-deficient mattes. The compositions of the matte, slag, and gaseous phases in equilibrium are calculated using Gibbs free energies of formation and the activity coefficients of the components derived from the experimental data. The model predictions are compared with the known industrial data from the Kalgoorlie Nickel Smelter (Kalgoorlie, Australia), the Outokumpu Harjavalta Nickel Smelter (Harjavalta, Finland), the INCO Metals Company (Sudbury, Canada), and from a number of experimental data. An excellent agreement is obtained. It was found that the distribution behaviors of Ni, Co, Cu, Fe, S, As, Sb, and Bi in the nickel smelting furnace depend on process parameters such as the smelting temperature, matte grade, and partial pressure of oxygen in the process.