An extensive algorithm based on both analytical and numerical solution methodologies is proposed to obtain transient temperature distributions in a three-dimensional living tissue subjected to a moving single-point and multi-point laser beam by considering metabolic heat generation and blood perfusion rate. Here, the dual-phase lag/Pennes equation is analytically solved by using the method of Fourier series and the Laplace transform. The ability to model single-point or multi-point laser beams as an arbitrary function of place and time is a significant advantage of the proposed analytical approach, which can be used to solve similar heat transfer problems in other living tissues. Besides, the related heat conduction problem is numerically solved based on the finite element method. The effects of laser beam transitional speed, laser power, and the number of laser points on the temperature distribution within the skin tissue are investigated. Moreover, the temperature distribution predicted by the dual-phase lag model is compared with that of the Pennes model under different working conditions. For the studied cases, it is observed that the maximum tissue temperature decreased about 63% by an increase of 6mm/s in the speed of the laser beam. An increase in the laser power from 0.8W/cm3 to 1.2W/cm3 results in a 28 °C increase in the maximum temperature of the skin tissue. It is observed that the maximum temperature predicted by the dual-phase lag model is always lower than that of the Pennes model and the temperature variations over time are sharper, while their results are entirely consistent over the simulation time. The obtained numerical results indicated that the dual-phase lag model is preferred in heating processes occurring at short intervals. Among the investigated parameters, the laser beam speed has the most considerable effect on the difference between the results of the Pennes and the dual-phase lag models. © 2022 Elsevier Ltd
