We report about a numerical algorithm based on the lattice Boltzmann method and its applications for simulations of turbulent convection in multi-phase flows. We discuss the issue of ’latent heat’ definition using a thermodynamically consistent pseudo-potential on the lattice. We present results of numerical simulations in 3D with and without boiling, showing the distribution of pressure, density and temperature fluctuations inside a convective cell.