This study investigates the influence of blood flow on temperature distribution during high-intensity focused ultrasound (HIFU) ablation of liver tumors. A three-dimensional acoustic-thermal-hydrodynamic coupling model is developed to compute the temperature field in the hepatic cancerous region. The model is based on the nonlinear Westervelt equation, bioheat equations for the perfused tissue and blood flow domains. The nonlinear Navier-Stokes equations are employed to describe the flow in large blood vessels. The effect of acoustic streaming is also taken into account in the present HIFU simulation study. A simulation of the Westervelt equation requires a prohibitively large amount of computer resources. Therefore a sixth-order accurate acoustic scheme in three-point stencil was developed for effectively solving the nonlinear wave equation. Results show that focused ultrasound beam with the peak intensity 2470 W/cm2 can induce acoustic streaming velocities up to 75 cm/s in the vessel with a diameter of 3 mm. The predicted temperature difference for the cases considered with and without acoustic streaming effect is 13.5 °C or 81% on the blood vessel wall for the vein. Tumor necrosis was studied in a region close to major vessels. The theoretical feasibility to safely necrotize the tumors close to major hepatic arteries and veins was shown.

The present study is aimed to investigate bubble dynamics in a soft tissue, to which HIFU's continuous harmonic pulse is applied by introducing a viscoelastic cavitation model. After a comparison of some existing cavitation models, we decided to employ Gilmore-Akulichev model. This chosen cavitation model should be coupled with the Zener viscoelastic model in order to be able to simulate soft tissue features such as elasticity and relaxation time. The proposed Gilmore-Akulichev-Zener model was investigated for exploring cavitation dynamics. The parametric study led us to the conclusion that the elasticity and viscosity both damp bubble oscillations, whereas the relaxation effect depends mainly on the period of the ultrasound wave. The similar influence of elasticity, viscosity and relaxation time on the temperature inside the bubble can be observed. Cavitation heat source terms (corresponding to viscous damping and pressure wave radiated by bubble collapse) were obtained based on the proposed model to examine the cavitation significance during the treatment process. Their maximum values both overdominate the acoustic ultrasound term in HIFU applications. Elasticity was revealed to damp a certain amount of deposited heat for both cavitation terms.

The subject of the present theoretical study is the dynamics of bubble-bubble interactions in a viscoelastic medium. First, new equations for calculating the viscoelastic drag exerted on bubbles during their translational motion in a viscoelastic medium are derived. The drag equations are incorporated in the bubble-bubble interaction model in which, thereby, both the translational and radial motions of the bubbles are affected by the viscoelastic features of the medium. Second, the derived equations are applied to investigate how the viscoelastic properties of the medium can affect the dynamics of multiple bubbles, as well as how the bubbles can affect each other. It was discovered that the bubble-bubble interaction can significantly influence the dynamics of a single bubble. As the distance between the bubbles increases, their effect on each other decreases, and at a distance of several millimeters, this effect can be neglected. Moreover, it was concluded that with increasing elasticity and viscosity of the medium, as well with decreasing relaxation time, the effects of other bubbles on the current bubble's radial motion can become negligible. The translational motion of the bubbles was investigated for different viscoelastic models. The elasticity resists the motion of bubbles in space, resulting in a dynamical steady state of the distance between the bubbles at high elasticity values. The relaxation time of the medium was also found to be important in terms of the bubbles' translational movement.