Hydrodynamics of multiple rising bubbles as a fundamental two-phase phenomenon is studied numerically by lattice Boltzmann method and using Lee two-phase model. Lee model based on Cahn-Hilliard diffuse interface approach uses potential form of intermolecular forces and isotropic finite difference discretization. This approach is able to avoid parasitic currents and leads to a stable procedure t...