Acoustic droplet vaporization (ADV) has been proven to enhance high intensity focused ultrasound (HIFU) thermal ablation of tumor. It has also been demonstrated that triggering droplets before HIFU exposure could be a potential way to control both the size and the shape of the thermal lesion. In this paper, a numerical model is proposed to predict the thermal lesion created in ADV enhanced HIFU treatment. Bubble oscillation was coupled into a viscoelastic medium in the model to more closely represent real applications in tissues. Several physical processes caused by continuous wave ultrasound and elevated temperature during the HIFU exposure were considered, including rectified diffusion, gas solubility variation with temperature in the medium, and boiling. Four droplet concentrations spanning two orders of magnitude were calculated. The bubble cloud formed from triggering of the droplets by the pulse wave ultrasound, along with the evolution of the shape and location of the bubble cloud and thermal lesion during the following continuous wave exposure was obtained. The increase of bubble void fraction caused by continuous wave exposure was found to be consistent with the experimental observation. With the increase of droplet concentration, the predicted bubble cloud shapes vary from tadpole to triangular and double triangular, while the thermal lesions move toward the transducer. The results show that the assumptions used in this model increased the accuracy of the results. This model may be used for parametrical study of ADV enhanced HIFU treatment and be further used for treatment planning and optimization in the future.