Simulation of nonlinear Westervelt equation for the investigation of acoustic streaming and nonlinear propagation effects

Process, dynamics and bioeffects of acoustic droplet vaporization induced by dual-frequency focused ultrasound

Acoustic droplet vaporization (ADV) plays a crucial role in ultrasound-related biomedical applications. While previous models have examined the stages of nucleation, growth, and oscillation in isolation, which may limit their ability to fully describe the entire ADV process. To address this, our study developed an integrated model that unifies these three stages of ADV, stimulated by a continuous nonlinear dual-frequency ultrasound wave. Using this integrated model, we investigated the influence of nonlinear dual-frequency ultrasound parameters on ADV dynamics and bioeffects by incorporating tissue viscoelasticity through parametric studies. Our results demonstrated that the proposed model accurately captured the entire ADV process, ensuring continuous vapor bubble formation and evolution throughout the phase transition process. Moreover, the applied unified theory for bubble dynamics can simulate intense bubble collapse with high Mach Number as a result of the nonlinear effects of dual-frequency ultrasound. In addition, cavitation-associated mechanical and thermal damage appeared to be more strongly correlated with rapid bubble collapse than with maximum bubble size. Our research also revealed that the mechanical and thermal effects could be regulated independently to some extent by adjusting dual-frequency ultrasound parameters, as they presented differing sensitivities to frequency and acoustic power. Importantly, dual-frequency combinations such as 1.5 MHz + 3 MHz (fundamental and second harmonic), which exhibit a higher Degree of Nonlinearity (DoN) can extend bubble lifespan, offering a potential pathway to the efficacy of ultrasound treatments.

The effect of temperature constraints on the treatment of tumors using focused ultrasound-induced acoustic streaming

The transport of drugs into tumor cells near the center of the tumor is known to be severely hindered due to the high interstitial pressure and poor vascularization. The aim of this work is to investigate the possibility to induce acoustic streaming in a tumor. Two tumor cases (breast and abdomen) are simulated to find the acoustic streaming and temperature rise, while varying the focused ultrasound transducer radius, frequency, and power for a constant duty cycle (1%). In the absence of perfusion, the simulated rise in temperature, despite the low duty cycle, never reaches a steady state and is fitted to a logarithmic equation, enabling predictions of the temperature for long treatment times. Higher frequencies and larger probe radii are found to result in shorter treatment times relative to the temperature rise, at the cost of a smaller treated area. Results from the simulations indicate that it may be possible to achieve reasonable acoustic streaming values in tumor without the temperature exceeding 50 °C. Treatment times for streaming a distance of 50 μm in the breast case are shown to range from less than one and a half hour to 93 h, depending on the probe settings.

A new strategy for modelling sonochemical reactors: Coupling of the non-linear Louisnard model with mass and heat transport equations with applications to cavitating viscous fluids

In this work, novel numerical models were developed and validated to offer new strategies in modelling sonochemical reactors. More specifically, in our original approach the non-linear Louisnard model was coupled with heat and mass transport equations to predict gradients in temperature and species concentration in a sonicated reactor. Additionally, a new operating window was investigated by modelling mixtures of increasing viscosity on both micro- and macroscale sonochemical effects. On the microscale, the effects of increasing viscosity on bubble dynamics were determined by solving the Keller-Miksis equation. Various cavitation threshold definitions were evaluated. The bubble collapse temperature was determined for all investigated mixtures and the energy dissipation of a single bubble was calculated. On the macroscale, different acoustic attenuation models were compared accounting for either linear or non-linear equations. Specifically, viscous losses were implemented in the non-linear Louisnard model, and model predictions were validated against experimental data. The model was able to predict multiple zones of cavitation in the reactor, as observed experimentally, and to estimate the dissipated energy for the different mixtures. Moreover, it was demonstrated that the cavitation-based attenuation dominates the other dissipation phenomena even for the most viscous solutions. The Louisnard model was coupled with heat transport equations, and using this extended version of the model, the temperature profiles were predicted for mixtures of increasing viscosity during sonication. Using a regression formula available in literature, radical production was related to the acoustic pressure field. By including reactions and mass transport in the acoustic model, for the first time in modelling ultrasonic reactors, the full distribution of light in the reactor during sonochemiluminescence (SCL) experiments for water was quantified.

Experimental Evaluation of an Ultrasound-Guided High-Intensity-Focused Ultrasound Probe for Sonication of Artery

OBJECTIVES

This study aimed to develop an ultrasound-guided high-intensity-focused ultrasound (USgHIFU) probe for arterial sonication and to evaluate vascular contraction.

METHODS

The USgHIFU probe comprised two confocal spherical transducers for sonication and a US color Doppler flow imaging probe for guidance. A vessel-mimicking phantom was sonicated in two directions. In the vascular radial direction, an isolated rabbit aorta embedded in ex vivo pork liver was sonicated at different acoustic powers (245 and 519 W), flow rates (25, 30, and 50 mL/minute), and sonication energies (519, 980, and 1038 J). Changes in the postsonication vessels were evaluated using US imaging, microscopic observation, and histopathological analysis.

RESULTS

Beam focusing along the vascular radial direction caused significant deformation of both tube walls (n = 4), whereas focusing along the axial direction only affected the contraction of the anterior wall (n = 4). The contraction index (Dc) of the vessel sonicated at 245 W and 980 J was 56.2 ± 9.7% (n = 12) with 25 mL/minute. The Dc of the vessel sonicated at 519 W and 1038 J was 56.5 ± 7.8% (n = 17). The Dc of the vessel sonicated at 519 J total energy was 18.3 ± 5.1% (n = 12).

CONCLUSION

The developed USgHIFU probe induced greater vascular contractions by covering a larger area of the vessel wall in the radial direction. Sonication energy affects vascular contraction through temperature elevation of the vessel wall. When the acoustic power was high, an increase in acoustic power, even with comparable sonication energy, did not result in greater vessel contraction.

Parameter effects on arterial vessel sonicated by high-intensity focused ultrasound: an ex vivo vascular phantom study

Objective. This study is aimed to explore the effects of vascular and sonication parameters on ex vivo vessel sonicated by high-intensity focused ultrasound. Approach. The vascular phantom embedding the polyolefin tube or ex vivo vessel was sonicated. The vascular phantom with 1.6 and 3.2 mm tubes was sonicated at three acoustic powers (2.0, 3.5, 5.3 W). The occlusion level of post-sonication tubes was evaluated using ultrasound imaging. The vascular phantom with the ex vivo abdominal aorta of rabbit for three flow rates (0, 5, 10 cm s−1) was sonicated at two acoustic powers (3.5 and 5.3 W). Different distances between focus and posterior wall (2, 4, 6 mm) and cooling times (0 and 10 s) were also evaluated. The diameter of the sonicated vessel was measured by B-mode imaging and microscopic photography. Histological examination was performed for the sonicated vessels. Main results. For the 5 cm s−1 flow rate, the contraction index of vascular diameter (Dc) with 5.3 W and 10 s cooling time at 2 mm distance was 39 ± 9% (n = 9). With the same parameters except for 0 cm s−1 flow rate, the Dc was increased to 45 ± 7% (n = 4). At 3.5 W, the Dc with 5 cm s−1 flow rate was 23 ± 15% (n = 4). The distance and cooling time influenced the lesion along the vessel wall. Significance. This study has demonstrated the flow rate and acoustic power have the great impact on the vessel contraction. Besides, the larger lesion covering the vessel wall would promote the vessel contraction. And the in vivo validation is required in the future study.

GPU-accelerated study of the inertial cavitation threshold in viscoelastic soft tissue using a dual-frequency driving signal

Inertial cavitation thresholds under two forms of ultrasonic excitation (the single- and dual-frequency ultrasound modes) are studied numerically. The Gilmore-Akulichev model coupled with the Zener viscoelastic model is used to model the bubble dynamics. The threshold pressures are determined with two criteria, one based on the bubble radius and the other on the bubble collapse speed. The threshold behavior is investigated for different initial bubble sizes, acoustic signal modes, frequencies, tissue viscosities, tissue elasticities, and all their combinations. Due to the large number of parameters and their many combinations (around 1.5 billion for each threshold criterion), all simulations were executed on graphics processing units to speed up the calculations. We used our own code written in the C++ and CUDA C languages. The results obtained demonstrate that using the dual-frequency signal mode can help to reduce the inertial cavitation threshold (in comparison to the single-frequency mode). The criterion based on the bubble size gives a lower threshold than the criterion using the bubble collapse speed. With an increase of the elasticity, the threshold pressure also increases, whereas changing the viscosity has a very small impact on the optimal threshold, unlike the elasticity. A detailed analysis of the optimal ultrasound frequencies for a dual-frequency driving signal found that for viscosities less than 0.02 Pa·s, the first optimal frequency, in general, is much smaller than the second optimal frequency, which can reach 1 MHz. However, for high viscosities, both optimal frequencies are similar and varied in the range 0.01-0.05 MHz. Overall, this study presents a detailed analysis of inertial cavitation in soft tissue under dual-frequency signal excitation. It may be helpful for the further development of different applications of biomedical ultrasound.

Finite-amplitude acoustic responses of large-amplitude vibration objects with complex geometries in an infinite fluid

High-intensity acoustic waves existing commonly in aeronautical and aerospace vehicles manifest nonlinear propagation behaviors. Large-amplitude vibration and irregular shape of the aerospace vehicles further complicate the acoustic responses. This paper is concerned with numerical analysis of finite-amplitude acoustic responses of complex-shaped vibration objects. The time-dependent effect of the solid boundary position due to the large-amplitude vibration of the objects is considered. A set of first-order differential equations is derived to govern the finite-amplitude acoustic wave. A fourth-order dispersion-relation-preserving finite difference formulation is employed to solve the nonlinear acoustic equations on a fixed Cartesian grid. Acoustic responses of the fluid and the vibration of the complex-shaped object are coupled by considering the compatibility conditions on the fluid-solid interface. A ghost-cell sharp-interface immersed boundary method is utilized to relax the conformity requirement between the computational grid and solid boundary. Numerical filters are employed in the computational procedure to suppress numerical oscillations. The present method is validated through several numerical tests. Numerical analysis of finite-amplitude acoustic responses of a complex-shaped object is performed. The nonlinear effect of a finite-amplitude acoustic wave, the time-dependent effect of solid boundary position, and the coupling effect between them on the propagation behaviors of nonlinear acoustic waves are discussed.

"HIFU Beam:" A Simulator for Predicting Axially Symmetric Nonlinear Acoustic Fields Generated by Focused Transducers in a Layered Medium

"HIFU beam" is a freely available software tool that comprises a MATLAB toolbox combined with a user-friendly interface and binary executable compiled from FORTRAN source code (HIFU beam. (2021). Available: http://limu.msu.ru/node/3555?language=en). It is designed for simulating high-intensity focused ultrasound (HIFU) fields generated by single-element transducers and annular arrays with propagation in flat-layered media that mimic biological tissues. Numerical models incorporated in the simulator include evolution-type equations, either the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation or one-way Westervelt equation, for radially symmetric ultrasound beams in homogeneous and layered media with thermoviscous or power-law acoustic absorption. The software uses shock-capturing methods that allow for simulating strongly nonlinear acoustic fields with high-amplitude shocks. In this article, a general description of the software is given along with three representative simulation cases of ultrasound transducers and focusing conditions typical for therapeutic applications. The examples illustrate major nonlinear wave effects in HIFU fields including shock formation. Two examples simulate propagation in water, involving a single-element source (1-MHz frequency, 100-mm diameter, 90-mm radius of curvature) and a 16-element annular array (3-MHz frequency, 48-mm diameter, and 35-mm radius of curvature). The third example mimics the scenario of a HIFU treatment in a "water-muscle-kidney" layered medium using a source typical for abdominal HIFU applications (1.2-MHz frequency, 120-mm diameter, and radius of curvature). Linear, quasi-linear, and shock-wave exposure protocols are considered. It is intended that "HIFU beam" can be useful in teaching nonlinear acoustics; designing and characterizing high-power transducers; and developing exposure protocols for a wide range of therapeutic applications such as shock-based HIFU, boiling histotripsy, drug delivery, immunotherapy, and others.

Accelerating frequency-domain numerical methods for weakly nonlinear focused ultrasound using nested meshes

The numerical simulation of weakly nonlinear ultrasound is important in treatment planning for focused ultrasound (FUS) therapies. However, the large domain sizes and generation of higher harmonics at the focus make these problems extremely computationally demanding. Numerical methods typically employ a uniform mesh fine enough to resolve the highest harmonic present in the problem, leading to a very large number of degrees of freedom. This paper proposes a more efficient strategy in which each harmonic is approximated on a separate mesh, the size of which is proportional to the wavelength of the harmonic. The increase in resolution required to resolve a smaller wavelength is balanced by a reduction in the domain size. This nested meshing is feasible owing to the increasingly localised nature of higher harmonics near the focus. Numerical experiments are performed for FUS transducers in homogeneous media to determine the size of the meshes required to accurately represent the harmonics. In particular, a fast volume potential approach is proposed and employed to perform convergence experiments as the computation domain size is modified. This approach allows each harmonic to be computed via the evaluation of an integral over the domain. Discretising this integral using the midpoint rule allows the computations to be performed rapidly with the FFT. It is shown that at least an order of magnitude reduction in memory consumption and computation time can be achieved with nested meshing. Finally, it is demonstrated how to generalise this approach to inhomogeneous propagation domains.

Experimental validation of acoustic and thermal modeling in heterogeneous phantoms using the hybrid angular spectrum method

PURPOSE

The aim was to quantitatively validate the hybrid angular spectrum (HAS) algorithm, a rapid wave propagation technique for heterogeneous media, with both pressure and temperature measurements.

METHODS

Heterogeneous tissue-mimicking phantoms were used to evaluate the accuracy of the HAS acoustic modeling algorithm in predicting pressure and thermal patterns. Acoustic properties of the phantom components were measured by a through-transmission technique while thermal properties were measured with a commercial probe. Numerical models of each heterogeneous phantom were segmented from 3D MR images. Cylindrical phantoms 30-mm thick were placed in the pre-focal field of a focused ultrasound beam and 2D pressure measurements obtained with a scanning hydrophone. Peak pressure, full width at half maximum, and normalized root mean squared difference (RMSDn) between the measured and simulated patterns were compared. MR-guided sonications were performed on 150-mm phantoms to obtain MR temperature measurements. Using HAS-predicted power density patterns, temperature simulations were performed. Experimental and simulated temperature patterns were directly compared using peak and mean temperature plots, RMSDn metrics, and accuracy of heating localization.

RESULTS

The average difference between simulated and hydrophone-measured peak pressures was 9.0% with an RMSDn of 11.4%. Comparison of the experimental MRI-derived and simulated temperature patterns showed RMSDn values of 10.2% and 11.1% and distance differences between the centers of thermal mass of 2.0 and 2.2 mm.

CONCLUSIONS

These results show that the computationally rapid hybrid angular spectrum method can predict pressure and temperature patterns in heterogeneous models, including uncertainties in property values and other parameters, to within approximately 10%.

Microstructure-based non-Fourier heat transfer modeling of HIFU treatment for thyroid cancer

BACKGROUND AND OBJECTIVES

High intensity focused ultrasound is an emerging non-invasive technique for the thermal ablation of cancer. Modeling of high intensity focused ultrasound as a method to induce hyperthermia, by considering non-equilibrium convective heat transfer has been under-represented in the previous studies. Therefore, in the present study, we aimed to study the effect of blood vessels during high intensity focused ultrasound ablation of thyroid cancer. In addition, high intensity focused ultrasound modeling was greatly improved by considering non-Fourier heat transfer.

METHODS

The modified dual-phase-lag model was used for the modeling of heat transfer in thyroid cancer during the ultrasound irradiation. The model parameters were linked with the tissue's microstructure parameters. Meanwhile, an interfacial convective heat transfer was considered between the blood vessels and the extravascular matrix. The extent of the vascular region was determined using the field emission scanning electron microscopy images. The non-linear Westervelt equation was solved for the sound wave to determine the heat source for the induced hyperthermia treatment.

RESULTS

Referring to the acoustic results, sharp-wave ripples were observed due to the inclusion of notable amplitudes of excited harmonics. The thermal results showed a maximum temperature rise of 25.08°C and 51.47°C at the powers of 5 W and 10 W using the modified dual-phase-lag model, while the Pennes model predicted a temperature rise of 28.77°C and 55.5°C at the same powers. It was also concluded that a constant blood temperature, overestimates the dissipated energy and the temperature reduction during the cooling period, as a 15% deviation in the tumor temperature was observed from the non-equilibrium state at 10.65 s exposure and 10 W power. Eventually, the calculation of the ablated volumes indicated that the volumes were up to 4.5 times larger by the Pennes model compared to the modified dual-phase-lag model.

CONCLUSIONS

It can be concluded from the results that there should be a serious concern on the high intensity focused ultrasound modeling based on the parameters of blood vessels. Based on the thermal maps, the cancerous tissue should be exposed to a higher energy level of ultrasound waves in order to cause the desired damage against the estimated energy level predicted by the Pennes model.

Nonlinear ultrasound simulation in an axisymmetric coordinate system using a k-space pseudospectral method

A full-wave model for nonlinear ultrasound propagation through a heterogeneous and absorbing medium in an axisymmetric coordinate system is developed. The model equations are solved using a nonstandard or k-space pseudospectral time domain method. Spatial gradients in the axial direction are calculated using the Fourier collocation spectral method, and spatial gradients in the radial direction are calculated using discrete trigonometric transforms. Time integration is performed using a k-space corrected finite difference scheme. This scheme is exact for plane waves propagating linearly in the axial direction in a homogeneous and lossless medium and significantly reduces numerical dispersion in the more general case. The implementation of the model is described, and performance benchmarks are given for a range of grid sizes. The model is validated by comparison with several analytical solutions. This includes one-dimensional absorption and nonlinearity, the pressure field generated by plane-piston and bowl transducers, and the scattering of a plane wave by a sphere. The general utility of the model is then demonstrated by simulating nonlinear transcranial ultrasound using a simplified head model.

Numerical Study on the Possible Scanning Pathways to Optimize Thermal Impacts During Multiple Sonication of HIFU

Findings from numerical investigation of thermal aspects of tissue phantoms subjected to high intensity focused ultrasound (HIFU) during the process of multiple sonications have been presented. Various scanning pathways are chosen to deliver the acoustic energy and to place multiple lesions at the desired locations within the tissue. A new scanning method based on the calculation of the largest distance between any two consecutive lesions is proposed and developed, which can possibly reduce the total treatment time by eliminating the cooling period. This new approach is compared with conventional scanning methods, namely, raster scan, spiral scan from the center outwards and spiral scan from the outside to the center. Various parameters that affect the thermal response of the tissue due to the generation of multiple lesions are defined and studied so that the energy that is given to each spot can be modulated. It is found that modulation of thermal energy provides better control over HIFU exposure and heating time modulation results into a reduced treatment time. The study highlighted the potential advantages of the proposed scanning method vis-à-vis the other conventional approaches. In particular, compared to the other methods, the proposed approach resulted in the realization of almost uniform distribution of thermal energy over the entire ROI leading to almost simultaneous destruction of the affected lesions. The present study can provide the requisite guidance for HIFU-based treatment planning and also highlights the need for optimizing the scanning pathway on the basis of total number of lesions and time parameters.

Simulation of cavitation enhanced temperature elevation in a soft tissue during high-intensity focused ultrasound thermal therapy

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³ 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.

Dynamics of bubble-bubble interactions experiencing viscoelastic drag

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.

Numerical and Experimental Evaluation of High-Intensity Focused Ultrasound-Induced Lesions in Liver Tissue Ex Vivo

OBJECTIVES

Recent advances in the field of acoustics and piezoelectric and ultrasound transducers have led to new approaches to the diagnosis and treatment of certain diseases. One method of treatment with ultrasonic waves is high-intensity focused ultrasound (HIFU) treatment, which is a thermal therapeutic method used to treat malignant tumors. Although a variety of treatment-planning strategies using ultrasonic waves have been investigated, little clinical success has been achieved. Computational modeling is a powerful tool for predicting device performance.

METHODS

The heating induced by a concave transducer with operating powers of 85 and 135 W was studied, and the experimental results presented in this article verify its applicability. Numerical simulations of the nonlinear acoustic field were performed by using the Westervelt and Khokhlov-Zabolotskaya-Kuznetsov equations. Heat transfer was measured for the 2 operational powers, and the results were compared with ex vivo experimental results. In addition, thermal dose contours for both the simulation and experimental results were calculated to investigate the ablated area.

RESULTS

Good agreement was found between the experimental and numerical results. The results show that the average temperature deviations calculated at the focal point were 12.8% and 4.3% for transducer powers of 85 and 135 W, respectively.

CONCLUSIONS

This study provides guidance to HIFU practitioners in determining lesion size and identifying nonlinear effects that should be considered in HIFU procedures.

Multi-resolution simulation of focused ultrasound propagation through ovine skull from a single-element transducer

Transcranial focused ultrasound (tFUS) is emerging as a non-invasive brain stimulation modality. Complicated interactions between acoustic pressure waves and osseous tissue introduce many challenges in the accurate targeting of an acoustic focus through the cranium. Image-guidance accompanied by a numerical simulation is desired to predict the intracranial acoustic propagation through the skull; however, such simulations typically demand heavy computation, which warrants an expedited processing method to provide on-site feedback for the user in guiding the acoustic focus to a particular brain region. In this paper, we present a multi-resolution simulation method based on the finite-difference time-domain formulation to model the transcranial propagation of acoustic waves from a single-element transducer (250 kHz). The multi-resolution approach improved computational efficiency by providing the flexibility in adjusting the spatial resolution. The simulation was also accelerated by utilizing parallelized computation through the graphic processing unit. To evaluate the accuracy of the method, we measured the actual acoustic fields through ex vivo sheep skulls with different sonication incident angles. The measured acoustic fields were compared to the simulation results in terms of focal location, dimensions, and pressure levels. The computational efficiency of the presented method was also assessed by comparing simulation speeds at various combinations of resolution grid settings. The multi-resolution grids consisting of 0.5 and 1.0 mm resolutions gave acceptable accuracy (under 3 mm in terms of focal position and dimension, less than 5% difference in peak pressure ratio) with a speed compatible with semi real-time user feedback (within 30 s). The proposed multi-resolution approach may serve as a novel tool for simulation-based guidance for tFUS applications.

Bubble dynamics in viscoelastic soft tissue in high-intensity focal ultrasound thermal therapy

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.

High-intensity focused ultrasound ablation around the tubing

High-intensity focused ultrasound (HIFU) has been emerging as an effective and noninvasive modality in cancer treatment with very promising clinical results. However, a small vessel in the focal region could be ruptured, which is an important concern for the safety of HIFU ablation. In this study, lesion formation in the polyacrylamide gel phantom embedded with different tubing (inner diameters of 0.76 mm and 3 mm) at varied flow speeds (17–339 cm/s) by HIFU ablation was photographically recorded. Produced lesions have decreased length (~30%) but slightly increased width (~6%) in comparison to that without the embedded tubing. Meanwhile, bubble activities during the exposures were measured by passive cavitation detection (PCD) at the varied pulse repetition frequency (PRF, 10–30 Hz) and duty cycle (DC, 10%-20%) of the HIFU bursts. High DC and low flow speed were found to produce stronger bubble cavitation whereas no significant influence of the PRF. In addition, high-speed photography illustrated that the rupture of tubing was produced consistently after the first HIFU burst within 20 ms and then multiple bubbles would penetrate into the intraluminal space of tubing through the rupture site by the acoustic radiation force. Alignment of HIFU focus to the anterior surface, middle, and posterior surface of tubing led to different characteristics of vessel rupture and bubble introduction. In summary, HIFU-induced vessel rupture is possible as shown in this phantom study; produced lesion sizes and shapes are dependent on the focus alignment to the tubing, flow speed, and tubing properties; and bubble cavitation and the formation liquid jet may be one of the major mechanisms of tubing rupture as shown in the high-speed photography.

Focused ultrasound actuation of shape memory polymers; acoustic-thermoelastic modeling and testing

Controlled drug delivery (CDD) technologies have received extensive attention recently.

Modeling of wave propagation for medical ultrasound: a review

Numerical modeling of medical ultrasound has advanced tremendously in the past two decades. This opens up a great number of opportunities for medical ultrasound and associated technologies. Numerous new governing equations and algorithms have emerged and been applied to studying various medical ultrasound applications, including ultrasound imaging, photo-acoustic imaging, and therapeutic ultrasound. In addition, thanks to the rapid development of computers, modeling acoustic wave propagation in three-dimensional, large-scale domains has become a reality. This article will provide an indepth literature and technical review of recent progress on numerical modeling of medical ultrasound. Future challenges will also be discussed.