  In this paper, we propose a lattice Boltzmann (LB) model to solve the coupled
Cahn-Hilliard-Navier-Stokes equations. Differently from previous efforts, the
LB equation for the fluid velocity is decomposed in a space of non-orthogonal
central moments where the surface tension force is inserted directly in the
equilibrium state. The present scheme is validated against well-consolidated
benchmark tests, showing very good accuracy. Moreover, it outperforms the BGK
approach in terms of stability, as our algorithm allows us to simulate a very
large range of viscosity contrast.
