Kinetic modeling is very important to understand steelmaking reactions. The coupled reaction model proposed by Prof. D.G.C Robertson is the best suited program, and therefore, the author has applied it to the various steelmaking reactions, a few examples of which are discussed. Hot metal dephosphorization occurs under non-equilibrium conditions as it takes place between the slag, which has a high oxygen potential, and the hot metal, which has a low oxygen potential. A process simulation model based on the coupled reaction model was constructed. This model considers the solid and liquid slag, liquid metal phases, and the reaction between the solid and liquid slag, besides the reaction between the liquid slag and liquid metal. Using this model, the ruling factor to increase the reaction efficiency of dephosphorization has been verified. In the secondary refining process, the values of inclusion composition calculated by the thermodynamic model, under the assumption of equilibrium with the metal composition, were found to be different from those observed in practice. A kinetic model, incorporating the reactions between slag, metal, inclusion, and refractory has been postulated, and the composition change of inclusion during treatment has been analyzed.