Numerical simulation of mold filling in low pressure die casting is considered in this study. The physical model includes modeling of free surface flow, heat transfer with phase change, surface tension, natural convection together with effect of trapped air in the mold. The governing equations are discretized by control volume finite difference method. The pressure field is computed by two-step projection method and the free surface is tracked by PLIC-VOF method. Water analog model is used for the validation purpose. Good agreement between numerical and experimental results is observed which supports the feasibility of the presented method.