A mathematical model is proposed to study the impingement of multiple supersonic jets onto the free surface of the liquid bath containing molten slag and metal in a steelmaking converter by means of volume of fluid (VOF) approach. The applicability of the proposed model is verified by the good agreement with measured and calculated results. The model is then used to study the jet hydrodynamic behaviours such as jet profile, impact force, penetration depth and impact area, as well as time-dependent evolution of cavity, with respect to temperature, operating pressure and lance height. The results show that the temperature effect on total impact force is negligible but significant for penetration depth and impact area. A semi-empirical correlation is used to describe penetration depth under the steelmaking condition of high temperature. It is also shown that the proposed model can reproduce the splashing phenomenon on the bath surface.