Two-phase flows in gas stirred steel bath and their effects on the carbon concentration distribution in the ladle during a decarburisation vacuum treatment are analysed. For this purpose, a numerical tool based on the coupling between the CFD Fluent software and kinetic model developed at IRSID has been used. The CFD Fluent software, equipped with a Lagrangian approach, provides a local description of the ladle hydrodynamics (molten steel flow and turbulence). The kinetic model is based on the physical description of the different reaction sites in the vacuum tank degasser. It is used to determine the evolution of the global carbon concentration in ladle. The coupling gives a complete mapping of the carbon concentration in the ladle according to time treatment. Numerical results applied to a water model configuration reveal good agreement in terms of velocity and turbulence profiles with experimental data provided from laser devices. In industrial configuration, the effect of the bath surface pressure and of the gas flow rate is analysed in terms of generated hydrodynamics and carbon distribution in the ladle.