A 3D mathematical model on the multiphase flow and splashing phenomenon in a 210 t Basic oxygen furnace converter is established. The turbulent multiphase flow of the molten steel, slag, and oxygen is simulated by the combined realizable k−ε model and volume of fluid model. The trajectory of argon bubbles and the interaction between the molten steel and argon is simulated by the two‐way coupled Lagrangian discrete phase model. Effects of the lance height and bottom‐blowing gas flow rate on the multiphase flow and splashing phenomenon are investigated. The results indicated that the cavity depth is decreased from 0.45 to 0.36 m when the lance height increases from 1700 to 2300 mm. The cavity depth is not significantly affected by the bottom‐blowing argon flow rate with the 44 000 m3 h−1 top‐blowing oxygen flow rate and 1700 mm lance height. The increased lance height resulted in a decrease in splashing. When the oxygen lance height increases from 1700 to 2300 mm, the amount of the splashed molten steel decreases from 3774.3 to 1953.0 kg.