In this study, a three-dimensional gravity casting problem has been examined to investigate a coupled phenomenon of the filling and solidification process. This work simultaneously considers the two key phenomena of metal casting : the fluid flow during mold filling, and solidification process. The VOF method is used to analyze the free surface flow during filling and the equivalent specific heat method is employed to model the latent heat release during solidification. The time-implicit filling algorithm is applied to save the computational time for analyzing the mold filling process. The three-dimensional benchmark problem used in the MCWASP VII has been solved using both the implicit and explicit algorithm, and the present results are compared with the benchmark experimental results and the other numerical results.