Steel refining is a complex phenomenon which depends on numerous variables, so, a kinetic approach is necessary for precise understanding of the refining process. In this study, based on a previously proposed model for hot metal dephosphorization, a new simulation model for the steel refining process in BOF is presented. In most cases, steelmaking slag is saturated with dicalcium-silicate (C 2 S) and it is well known that C 2 S forms solid solution with tricalcium-phosphate (C 3 P) in a wide composition range and the partition ratio of phosphorus between C 2 S and liquid slag is large. On the other hand, C 2 S formed around the lime surface is known as a barrier to lime dissolution into liquid slag. In this simulation model not only the effect of solid phase in slag is considered but also the effects of temperature dependence of variables as well as top and bottom blowing and scrap melting are taken into account.The calculation results are compared with industrial data and the good agreement between experimental and simulation results evidence the validity of this kinetic approach to steel refining process in BOF. Moreover, by using this model the influence of various parameters on the reaction efficiency is discussed. steel research int. 81 (2010) No. 8 Process Metallurgy Figure 8. a) Effect of Ar flow rate and b) effect of oxygen flow rate on the relation between carbon content and free oxygen.