The simulation of multiphase flows is an outstanding challenge, due to the inherent complexity of the underlying physical phenomena and to the fact that multiphase flows are very diverse in nature, and so are the laws governing their dynamics. In the last two decades, a new class of mesoscopic methods, based on minimal lattice formulation of Boltzmann kinetic equation, has gained significant interest as an efficient alternative to continuum methods based on the discretisation of the NS equations for non ideal fluids.