Modeling the evolution of the melt front under gravity in the presence of a horizontal thermal gradient is a challenging issue, hitherto tackled exclusively with the concepts and tools of computational continuum thermomechanics, too phenomenologically driven to have satisfactory predictive capabilities. Here, we show that this complex phenomenon is amenable to treatment by the methods and tools of Non-Equilibrium Molecular Dynamics (NEMD). To do so, we addressed all the difficulties caused by the necessity of applying suitable boundary conditions and minimizing surface effects so that the bulk behavior of the system in non-equilibrium conditions can be detected. Sufficient adiabatic separation of the time scales permits us to use macroscopically relatively short-but microscopically long enough-time averages to get the macroscopic bulk behavior of the system accurately. To get an adequate signal-to-noise ratio, we had to use an unphysically large value of the gravity. However, we know from NEMD simulations in transport studies that the phenomena produced are stable over many orders of magnitude. In conclusion, our work proves that molecular simulation can be a good tool to study this family of non-equilibrium phenomena, although further work is needed to achieve quantitative predictive capabilities.