The numerical modeling of a binary solidification with a mushy layer mechanism is considered in this manuscript. The nonlinear coupled system of equations describes the heat and mass diffusions of a one-dimensional spatial variable in the semiinfinite interval. Also formulated is a transformed system in a finite interval. We propose numerical methods for solving the nonlinear system using a threshold strategy based on fixed computation-domain approach. Our calculated results and those from the LeadEx field experiment are well-matched in their tendencies. ? 2012 Global-Science Press.