In the present article, we extend the dynamic model of Lee to investigate film boiling numerically based on the Lattice Boltzmann method. The phase-change process is modeled by incorporating a proper source term at the phase interface, and then the extension of the classical convective Cahn–Hilliard’s equation in the presence of phase change is employed in the multiphase LB framework of Lee. In the first step, the two-dimensional Rayleigh–Taylor instability is modeled. Comparing the results with the other computational results reveals the consistency of our method. Then the model is extended to simulate single- and multiple-mode film boiling. The results of simulation are compared with the previous analytical correlations. In addition, the effects of Jacob number on the interface dynamic, Nusselt number, film thickness and maximum velocity are illustrated. In this simulation, the bubble temperature is assumed not to be constant, which makes the simulation more usable and realistic.