CFD models have been developed and numerical simulations have been carried out to predict the formation of foam in oxygen steelmaking. Foam was considered as a separate phase comprising a mixture of gas and liquid. Bubble break up and coalescence models have also been incorporated in a CFD model to predict the number density of individual bubble classes. A population balance equation was used to track the number density of each bubble class. Decarburization with heat generation from chemical reactions was integrated in the process. User subroutines were written in FORTRAN to incorporate the foam formation, the bubble break up and coalescence rate and decarburisation in the main program. The model predicted the foam height, bubble number density, velocity of phases, decarburization, and turbulence. The result from the model has been compared with available data from literature and found to be in reasonable agreement with the experimental and plant data.