In the current work we construct a nonlocal mathematical model describing the phase transition occurs during the resistance spot welding process in the industry of metallurgy. We then consider a time discretization scheme for solving the resulting nonlocal moving boundary problem. The scheme consists of solving at each time step a linear elliptic partial differential equation and then making a correction to account for the nonlinearity. The stability and error estimates of the developed scheme are investigated. Finally some numerical results are presented confirming the efficiency of the developed numerical algorithm.