A numerical Euler–Lagrange method for bubble tower CO2 dissolution modeling