As a surface finishing technique for rapid remelting and re-solidification, laser polishing can effectively eliminate the asperities so as to approach the feature size. Nevertheless, the polished surface quality is significantly sensitive to the processing parameters, especially with respect to melt hydrodynamics. In this paper, a transient two-dimensional model was developed to demonstrate the molten flow behavior for different surface morphologies of the Ti6Al4V alloy. It is illustrated that the complex evolution of the melt hydrodynamics involving heat conduction, thermal convection, thermal radiation, melting and solidification during laser polishing. Results show that the uniformity of the distribution of surface peaks and valleys can improve the molten flow stability and obtain better smoothing effect. The high cooling rate of the molten pool resulting in a shortening of the molten lifetime, which prevents the peaks from being removed by capillary and thermocapillary forces. It is revealed that the mechanism of secondary roughness formation on polished surface. Moreover, the double spiral nest Marangoni convection extrudes the molten to the outsides. It results in the formation of expansion and depression, corresponding to nearby the starting position and at the edges of the polished surface. It is further found that the difference between the simulation and experimental depression depths is only about 2 μm. Correspondingly, the errors are approximately 8.3%, 14.3% and 13.3%, corresponding to Models 1, 2 and 3, respectively. The aforementioned results illustrated that the predicted surface profiles agree reasonably well with the experimentally measured surface height data.