This paper deals with an extended analytical solution for the mixture solidification problem, in which temperature is inherently coupled with the solute transport due to the presence of volume contraction induced flow. A new exact solution to the energy equation accounting for the convection effect in the melt is successfully derived, which allows the present analysis to cover a high initial superheating. Difference in properties between the solid and liquid phases is rigorously incorporated into the model equations in the solid fraction weighted form. Taking advantage of linearized correction factors, a systematic and easy-to-implement algorithm for determining the solidus and liquidus positions is introduced, which proves not only to converge stably but also to be very efficient. For a specific case, the present results show excellent agreements with the existing solution. The effect of convection in the melt becomes appreciable with increasing the initial superheating. It is revealed that variable properties in the mushy region significantly affect the solidification behaviors. The present study is also capable of resolving the interaction between microsegregation and macrosegregation.