SPH-based numerical technique for modelling of impact of molten drops on a surface with heat transfer and phase transitions effects is proposed. Computational algorithm uses SPH with procedure of restoring of particle consistence and variational approach to calculation of acceleration field. Also, boundary algorithm for free and contact surfaces in 3D setting are developed