In this study, numerical investigation has been performed on the evolution of weld pool geometry with moving free surface during low-energy density laser welding process. The free surface elevates near the weld pool edge and descends at the center of the weld pool if d$sigma$/dT is dominantly negative. It is shown that the predicted weld pool width and depth with moving free surface are a little greater than those with flat weld pool surface. It is also believed that the weld pool surface oscillation during the melting process augments convective heat transfer rate in the weld pool. The present analysis with moving free surface should be considered when We number is very small compared to 1.0 since the deformation of the weld pool surface is noticeable as We number decreases.