Premature melt solidification is the solidification of a melt during mold filling. In this study, a numerical model is used to analyze the influence of the pouring process on the premature solidification. The numerical model considers three phases, namely, air, melt, and equiaxed crystals. The crystals are assumed to have originated from the heterogeneous nucleation in the undercooled melt resulting from the first contact of the melt with the cold mold during pouring. The transport of the crystals by the melt flow, in accordance with the so-called “big bang” theory, is considered. The crystals are assumed globular in morphology and capable of growing according to the local constitutional undercooling. These crystals can also be remelted by mixing with the superheated melt. As the modeling results, the evolutionary trends of the number density of the crystals and the volume fraction of the solid crystals in the melt during pouring are presented. The calculated number density of the crystals and the volume fraction of the solid crystals in the melt at the end of pouring are used as the initial conditions for the subsequent solidification simulation of the evolution of the as-cast structure. A five-phase volume-average model for mixed columnar-equiaxed solidification is used for the solidification simulation. An improved agreement between the simulation and experimental results is achieved by considering the effect of premature melt solidification during mold filling. Finally, the influences of pouring parameters, namely, pouring temperature, initial mold temperature, and pouring rate, on the premature melt solidification are discussed.