This work aims at figuring out the influence of gas bubble size distribution on the ladle stirring process. The work is conducted through three-dimensional (3D) numerical simulation based on the finite volume method. Mesh sensitivity test and the cross-validation are performed to ensure the results are mesh independent and the numerical set-up is correct. Two distributions, uniform and Log-normal function, are investigated under different gas flow rates and number of porous plugs. The results indicate that the results, e.g., the axial velocity and the area of the slag eye, have little difference for low flow rate. The difference becomes dominant whilst the flow rate is increasing, such as 600 NL/min. The Log-normal function bubble size distribution gives a larger axial velocity and a smaller slag eye area compared to the uniform bubble size distribution. This work indicated that, at a higher flow rate, the Log-normal function is a better choice to predict the melt behavior and the slag open eye in the ladle refining process if the bubble interaction is not considered.