“…In comparison to analytical methods, numerical methods could only provide results approximately within an acceptable error tolerance, but they are more flexible to deal with the complicated yet practical situations such as the transient problem of multiple heat sources moving in a complex geometry of the material with time-dependent speeds [3]. However, most of numerical studies, regardless using meshless methods [13,14] or meshbased methods such as the finite element method [6,10], were concerned about problems involving only a heat source moving along a straight line with a constant speed, or multiple heat sources moving along parallel straight lines with the same constant speed. Apart from these, the technique of moving coordinate system, such that the heat source is stationary in the new coordinate system, is often introduced in both analytical and numerical analyses of the quasistationary problem [1,13].…”