A new Lagrangian Particle Tracking algorithm for LES of bubble plumes is presented.
Delta functions are applied for the first time to the knowledge of the authors to this kind of problem.
A new mapping approach based on the use of the bubble to mesh size ratio is suggested.
We successfully validate the time-averaged velocities and second-order turbulence statistics.
Slip velocity, interfacial forces and confinement effect are also analysed.