Based on the two-phase model (Eulerian-Eulerian model), the three dimensional fluid flow in water and that liquid steel systems stirred by one or two multiple gas jets are simulated. In the Eulerian-Eulerian two-phase model, the gas and the liquid phase are considered to be two different continuous fluids interacting with each other through the finite inter-phase areas. The exchange between the phases is represented by source terms in conversation equations. Turbulence is assumed to be a property of the liquid phase. A new turbulence modification k-ε model is introduced to consider the bubbles movement contribution to k and ε. The dispersion of phases due to turbulence is represented by introducing a diffusion term in mass conservation equation. The mathematical simulation agrees well with the experiment results. The study results indicate that the distance of two nozzles has big effect on fluid flow behavior in the vessel. Using two gas injection nozzles at the half radii of one diameter of the bottom generates a much better mixing than with one nozzle under the condition of the same total gas flow rate.