In this paper, an algorithm for simulating fluid flow and heat transfer for mold filling of shaped castings is presented. The main features of the algorithm include: 1) a simple but practical technique based on VOF method to determine free surface, 2) an explicit scheme of enthalpy to solve the energy equation more efficiently, and 3) an effective treatment to modify the flux deviation due to pressure iteration. In order to verify these methods, well controlled experiments have been repeatedly done with both water analog and gray iron pouring experiments to record the flow patterns and temperature variations. The calculated results are in accordance with the experimental ones. For the applications, the simulated initial temperature distribution right after mold filling was used to analyse subsequent solidification and to predict shrinkage defects.Actual castings were poured and tested in a foundry plant. The reuslts show that the defects predication with considering fluid flow effects is more precise than that without considering the effects.