The article presents the results of numerical simulation of an experiment on irradiating ex vivo bovine liver sample by the therapeutic array of the MR-HIFU clinical system (Sonalleve V1 3.0T, Profound Medical Corp., Canada). Continuous quasi-linear and pulsed shock-wave exposures with the same time-averaged power are compared. Volumetric thermal lesions were generated by moving the focus of the array in its focal plane along discrete trajectories consisting of two or four concentric circles with a maximum radius of 4 mm. The effect of using the criteria for controlling the thermal dose during treatment and ending the sonication on the shape, volume, and exposure time of generating thermal lesion were analyzed. The acoustic field in tissue was calculated using the Westervelt equation; the temperature field was simulated with the inhomogeneous heat conduction equation; and the lesion boundary was determined according to the thermal dose threshold. In the quasi-linear mode corresponding to the clinical one, thermal diffusion leads to elongation of the lesion by a factor of 2–3 along the beam axis compared to the transverse dimension of the trajectory. The use of pulsed shock-wave exposures with switching off the inner circles of the trajectory as they reach the threshold value of the thermal dose makes it possible to significantly suppress the thermal diffusion effects in the axial direction of the beam and obtain localized thermal lesion of a given shape with a thermal ablation rate comparable to the clinical case.