In this work, a new numerical model for the treatment of phase change phenomena is developed. The LBM method is used to install this new model based on determining the image of distribution functions during the phase transition. The classical Stefan problem modeling the solid / liquid phase change is studied with this approach and validated by comparison with analytical results.