The present study aims to investigate temperature distribution caused by bubble oscillations in a soft tissue during focused ultrasound therapy by introducing a coupled temperature-cavitation model. The proposed model is capable of describing bubble dynamics, viscoelastic properties of surrounding tissue-like medium, temperature distribution inside and outside the bubble, vapor diffusion within the bubble and vapor flux through the bubble wall to the exterior. The continuous temperature distribution inside and outside the oscillating bubble in soft tissue subject to ultrasound wave with high acoustic pressure is presented. The temperature close to the bubble wall can reach the value of about 10 3 K. The elasticity of soft tissue reduces temperature values. The relaxation time effect strongly depends on the period of the ultrasound wave. If the vapor mass flux effect is taken into account in the simulations, the rectified growth effect can be observed, which can lead to the decrease of the temperature values. Due to the growth of the bubble, the effects of elasticity and relaxation time on the temperature become less prominent during several bubble oscillation cycles. The impact of cavitation heat source terms on the exterior temperature was examined and led us to draw conclusion that, even though these heat sources can increase the outside temperature values, they can not be treated as main mechanisms for the temperature elevation during a few microseconds. The performed comparison with uncoupled conventional model for the outside temperature calculation revealed that coupling with inside temperature model delivers incomparably higher values to the bubble's exterior and, therefore, it is essential for the accurate description of the treatment process.
