A numerical algorithm for the time-independent Stefan problem is presented. The heat conducting plate is assumed to be infinitely broad and the heat source of constant intensity is assumed to move with a constant velocity in the plate.
The shape of molten pool is determined by correcting successively the provisional shape of the molten pool given beforehand until it satisfies the Stefan condition. The isotherms showing the temperature field are determined successively by starting from the obtained fusion boundary. The algorithm is given for the two-dimensional case and for the axi-symmetric three-dimensional case. Some examples of the results are given for the case of a linear heat source and a point heat source of various intensities.