In the drop-on-demand 3D printing technique, molten drop impact is crucial to the understanding of the formation of the building block of a part. However, existing studies cast less due attention to the rebounding dynamics of molten drops. This paper thus develops a phase field method to that end. The model consists of the Cahn-Hilliard equation, coupled with the Navier-Stokes equation to track the evolving liquid-gas interface. A heat source is added to the energy balance equation to distinguish between solid and fluid. The model is discretized using a finite difference method on a semi-staggered grid, and is paralleled based on the shared memory parallelism, OpenMP. The model was first compared with an experiment before being applied to various impacts. The agreement is reasonable. Typical impacts in 3D printing conditions were considered. The major findings are that the Tin drop rebounds off the substrate with a dimpled bottom when the impact velocity is 1 m/s, but with a rounded one when the impact velocity is reduced by half, and that the dimpled bottom is caused by an entrapped bubble during the receding process. Besides, effect of solidification on drop rebounding was examined as well. [ABSTRACT FROM AUTHOR]