Modeling of co-axial bubbles coalescence under moderate Reynolds regimes: A Bi-phase SPH approach