Boiling is a complex phenomenon where different non-linear physical interactions take place and for which the quantitative modeling of the mechanism involved is not fully developed yet. In the last years, many works have been published focusing on the numerical analysis of this problem. However, a lack of numerical works assessing quantitatively the sensitivity of these numerical simulations to grid parameters can be identified, especially for the Lattice Boltzmann method (LBM). The main goal of this work is to propose a mesh refinement methodology for simulating phase-change heat transfer problems by means of the pseudopotential LBM. This methodology was based on relating the physical parameters to their lattice counterparts for an arbitrary mesh under the viscous regime (where Δ ∝ Δ 2 ). A suitable modification of the EOS parameters and the adjusting of thermodynamic consistency and surface tension for a certain Δ were the main steps of the proposed methodology. A first ensemble of simple simulations including the droplet vaporization and the Stefan problems was performed to validate the proposed method and to assess the influence of some physical mechanisms. Global norms in space and time were used to evaluate the variations of both the density and temperature fields for pool boiling simulations when the lattice discretization is refined. It was observed that the proposed methodology provides convergent results for all the problems considered, and the convergence orders depend on the complexity of the simulated phenomena.