A new model is developed to describe the coalescence and breakup between bubbles and inclusions. In this model, the density of bubble attached inclusions is calculated through mass conservation equation. The momentum exchange after breakup or coalescence is derived through momentum conservation equation, which are tracked by discrete particle method (DPM). Three continuous phases (air-slag-steel) are considered in a ladle, and the unsteady turbulent flow is computed through k-e method. What's more, bubble expansion due to decreasing of hydraulic pressure is also taken into account. Results show that due to the bubble expansion, bubble density is mostly decreased to 0.3 to 0.8 kg/m 3 . However, after attaching some inclusions, the bubble density significantly rises to 20 to 60 kg/m 3 . The optimal bubble diameters attaching inclusions are ranged between 1.5 and 10 mm. For traditional slot-slot matched tuyeres (S-S mode), the inclusion removement ratio is 29.48 pct; By comparison, after employing the slot-porous matched tuyeres (S-P mode), the inclusion removements ratio rises to 36.34 pct. The mixing time is also shortened after adopting the S-P mode. The reason for this phenomenon is because the slot tuyere produces a strong asymmetry stream that drives more liquid to flow at the bottom of the ladle. And at the same time, the porous tuyere produces more fine bubbles to entrap more inclusions to the top. Taking advantages of porous and slot tuyeres, the mixing behavior and inclusion removements improves a lot. The result is beneficial for improving ladle refining.