Vasant A. Salgaonkar

Passive cavitation imaging with ultrasound arrays

A method is presented for passive imaging of cavitational acoustic emissions using an ultrasound array, with potential application in real-time monitoring of ultrasound ablation. To create such images, microbubble emissions were passively sensed by an imaging array and dynamically focused at multiple depths. In this paper, an analytic expression for a passive image is obtained by solving the Rayleigh-Sommerfield integral, under the Fresnel approximation, and passive images were simulated. A 192-element array was used to create passive images, in real time, from 520-kHz ultrasound scattered by a 1-mm steel wire. Azimuthal positions of this target were accurately estimated from the passive images. Next, stable and inertial cavitation was passively imaged in saline solution sonicated at 520 kHz. Bubble clusters formed in the saline samples were consistently located on both passive images and B-scans. Passive images were also created using broadband emissions from bovine liver sonicated at 2.2 MHz. Agreement was found between the images and source beam shape, indicating an ability to map therapeutic ultrasound beams in situ. The relation between these broadband emissions, sonication amplitude, and exposure conditions are discussed.

Acoustic emissions during 3.1 MHz ultrasound bulk ablation in vitro

Acoustic emissions associated with cavitation and other bubble activity have previously been observed during ultrasound (US) ablation experiments. Because detectable bubble activity may be related to temperature, tissue state and sonication characteristics, these acoustic emissions are potentially useful for monitoring and control of US ablation. To investigate these relationships, US ablation experiments were performed with simultaneous measurements of acoustic emissions, tissue echogenicity and tissue temperature on fresh bovine liver. Ex vivo tissue was exposed to 0.9-3.3-s bursts of unfocused, continuous-wave, 3.10-MHz US from a miniaturized 32-element array, which performed B-scan imaging with the same piezoelectric elements during brief quiescent periods. Exposures used pressure amplitudes of 0.8-1.4 MPa for exposure times of 6-20 min, sufficient to achieve significant thermal coagulation in all cases. Acoustic emissions received by a 1-MHz, unfocused passive cavitation detector, beamformed A-line signals acquired by the array, and tissue temperature detected by a needle thermocouple were sampled 0.3-1.1 times per second. Tissue echogenicity was quantified by the backscattered echo energy from a fixed region-of-interest within the treated zone. Acoustic emission levels were quantified from the spectra of signals measured by the passive cavitation detector, including subharmonic signal components at 1.55 MHz, broadband signal components within the band 0.3-1.1 MHz and low-frequency components within the band 10-30 kHz. Tissue ablation rates, defined as the thermally ablated volumes per unit time, were assessed by quantitative analysis of digitally imaged, macroscopic tissue sections. Correlation analysis was performed among the averaged and time-dependent acoustic emissions in each band considered, B-mode tissue echogenicity, tissue temperature and ablation rate. Ablation rate correlated significantly with broadband and low-frequency emissions, but was uncorrelated with subharmonic emissions. Subharmonic emissions were found to depend strongly on temperature in a nonlinear manner, with significant emissions occurring within different temperature ranges for each sonication amplitude. These results suggest potential roles for passive detection of acoustic emissions in guidance and control of bulk US ablation treatments.

Ultrasound-enhanced penetration of topical riboflavin into the corneal stroma.

PURPOSE To determine whether ultrasound treatment can promote the permeation of topical riboflavin into the corneal stroma. METHODS Fresh cadaveric rabbit eyes with intact epithelium were left for 45 minutes in riboflavin 0.1% solution and divided in the following groups: A--untreated, epithelium-on; B--ultrasound-treated (1 W/cm(2) at 880 kHz for 6 minutes) with epithelium-on; and C--epithelium-off (no ultrasound). Eyes were removed from the riboflavin solution, corneas were excised, and group B was divided into B1 (with epithelium maintained) and B2 (epithelium removed for the fluorescence analysis). Confocal microscopy was performed to quantify the fluorescence intensity in the cornea according to the distance from the surface (with epithelium in groups A and B1; without epithelium in groups B2 and C). RESULTS The average fluorescence intensity of riboflavin at a depth of 100, 150, 200, and 250 μm was 69.97, 58.83, 49.23, and 41.72 arbitrary units (A.U.) in group A, respectively; 255.26, 206.01, 159.81, 124.20 A.U. in group B1; 218.90, 177.90, 141.43, 110.45 A.U. in group B2; and 677.64, 420.10, 250.72 and 145.07 A.U. in group C. The difference in fluorescence was statistically significant between groups A and B1 (P = 0.001) and groups B2 and C (P < 0.0001). CONCLUSIONS Ultrasound treatment increased the entry of topical riboflavin into the corneal stroma despite the presence of a previously intact epithelial barrier. This approach may offer a means of achieving clinically useful concentrations of riboflavin within the cornea with minimum epithelial damage, thereby improving the risk profile of corneal cross-linking procedures.

Interstitial ultrasound ablation of vertebral and paraspinal tumours: Parametric and patient-specific simulations

Abstract Purpose: Theoretical parametric and patient-specific models are applied to assess the feasibility of interstitial ultrasound ablation of tumours in and near the spine and to identify potential treatment delivery strategies. Methods: 3D patient-specific finite element models (n = 11) of interstitial ultrasound ablation of tumours associated with the spine were generated. Gaseous nerve insulation and various applicator configurations, frequencies (3 and 7 MHz), placement trajectories, and tumour locations were simulated. Parametric studies with multilayered models investigated the impacts of tumour attenuation, tumour dimension, and the thickness of bone insulating critical structures. Temperature and thermal dose were calculated to define ablation (>240 equivalent minutes at 43 °C (EM43 °C)) and safety margins (<45 °C and <6 EM43 °C), and to determine performance and required delivery parameters. Results: Osteolytic tumours (≤44 mm) encapsulated by bone could be successfully ablated with 7 MHz interstitial ultrasound (8.1–16.6 W/cm2, 120–5900 J, 0.4–15 min). Ablation of tumours (94.6–100% volumetric) 0–14.5 mm from the spinal canal was achieved within 3–15 min without damaging critical nerves. 3 MHz devices provided faster ablation (390 versus 930 s) of an 18 mm diameter osteoblastic (high bone content) volume than 7 MHz devices. Critical anatomy in proximity to the tumour could be protected by selection of appropriate applicator configurations, active sectors, and applied power schemas, and through gaseous insulation. Preferential ultrasound absorption at bone surfaces facilitated faster, more effective ablations in osteolytic tumours and provided isolation of ablative energies and temperatures. Conclusions: Parametric and patient-specific studies demonstrated the feasibility and potential advantages of interstitial ultrasound ablation treatment of paraspinal and osteolytic vertebral tumours.

Applicators for magnetic resonance-guided ultrasonic ablation of benign prostatic hyperplasia.

ObjectivesThe aims of this study were to evaluate in a canine model applicators designed for ablation of human benign prostatic hyperplasia (BPH) in vivo under magnetic resonance imaging (MRI) guidance, including magnetic resonance thermal imaging (MRTI), determine the ability of MRI techniques to visualize ablative changes in prostate, and evaluate the acute and longer term histologic appearances of prostate tissue ablated during these studies. Materials and MethodsAn MRI-compatible transurethral device incorporating a tubular transducer array with dual 120° sectors was used to ablate canine prostate tissue in vivo, in zones similar to regions of human BPH (enlarged transition zones). Magnetic resonance thermal imaging was used for monitoring of ablation in a 3-T environment, and postablation MRIs were performed to determine the visibility of ablated regions. Three canine prostates were ablated in acute studies, and 2 animals were rescanned before killing at 31 days postablation. Acute and chronic appearances of ablated prostate tissue were evaluated histologically and were correlated with the MRTI and postablation MRI scans. ResultsIt was possible to ablate regions similar in size to enlarged transition zone in human BPH in 6 to 18 minutes. Regions of acute ablation showed a central “heat-fixed” region surrounded by a region of more obvious necrosis with complete disruption of tissue architecture. After 31 days, ablated regions demonstrated complete apparent resorption of ablated tissue with formation of cystic regions containing fluid. The inherent cooling of the urethra using the technique resulted in complete urethral preservation in all cases. ConclusionsProstatic ablation of zones of size and shape corresponding to human BPH is possible using appropriate transurethral applicators using MRTI, and ablated tissue may be depicted clearly in contrast-enhanced magnetic resonance images. The ability accurately to monitor prostate tissue heating, the apparent resorption of ablated regions over 1 month, and the inherent urethral preservation suggest that the magnetic resonance–guided techniques described are highly promising for the in vivo ablation of symptomatic human BPH.

Model-based feasibility assessment and evaluation of prostate hyperthermia with a commercial MR-guided endorectal HIFU ablation array

PURPOSE Feasibility of targeted and volumetric hyperthermia (40-45 °C) delivery to the prostate with a commercial MR-guided endorectal ultrasound phased array system, designed specifically for thermal ablation and approved for ablation trials (ExAblate 2100, Insightec Ltd.), was assessed through computer simulations and tissue-equivalent phantom experiments with the intention of fast clinical translation for targeted hyperthermia in conjunction with radiotherapy and chemotherapy. METHODS The simulations included a 3D finite element method based biothermal model, and acoustic field calculations for the ExAblate ERUS phased array (2.3 MHz, 2.3 × 4.0 cm(2), ∼1000 channels) using the rectangular radiator method. Array beamforming strategies were investigated to deliver protracted, continuous-wave hyperthermia to focal prostate cancer targets identified from representative patient cases. Constraints on power densities, sonication durations and switching speeds imposed by ExAblate hardware and software were incorporated in the models. Preliminary experiments included beamformed sonications in tissue mimicking phantoms under MR temperature monitoring at 3 T (GE Discovery MR750W). RESULTS Acoustic intensities considered during simulation were limited to ensure mild hyperthermia (Tmax < 45 °C) and fail-safe operation of the ExAblate array (spatial and time averaged acoustic intensity ISATA < 3.4 W/cm(2)). Tissue volumes with therapeutic temperature levels (T > 41 °C) were estimated. Numerical simulations indicated that T > 41 °C was calculated in 13-23 cm(3) volumes for sonications with planar or diverging beam patterns at 0.9-1.2 W/cm(2), in 4.5-5.8 cm(3) volumes for simultaneous multipoint focus beam patterns at ∼0.7 W/cm(2), and in ∼6.0 cm(3) for curvilinear (cylindrical) beam patterns at 0.75 W/cm(2). Focused heating patterns may be practical for treating focal disease in a single posterior quadrant of the prostate and diffused heating patterns may be useful for heating quadrants, hemigland volumes or even bilateral targets. Treatable volumes may be limited by pubic bone heating. Therapeutic temperatures were estimated for a range of physiological parameters, sonication duty cycles and rectal cooling. Hyperthermia specific phasing patterns were implemented on the ExAblate prostate array and continuous-wave sonications (∼0.88 W/cm(2), 15 min) were performed in tissue-mimicking material with real-time MR-based temperature imaging (PRFS imaging at 3.0 T). Shapes of heating patterns observed during experiments were consistent with simulations. CONCLUSIONS The ExAblate 2100, designed specifically for thermal ablation, can be controlled for delivering continuous hyperthermia in prostate while working within operational constraints.

Modelling of endoluminal and interstitial ultrasound hyperthermia and thermal ablation: Applications for device design, feedback control and treatment planning

Abstract Endoluminal and catheter-based ultrasound applicators are currently under development and are in clinical use for minimally invasive hyperthermia and thermal ablation of various tissue targets. Computational models play a critical role in device design and optimisation, assessment of therapeutic feasibility and safety, devising treatment monitoring and feedback control strategies, and performing patient-specific treatment planning with this technology. The critical aspects of theoretical modelling, applied specifically to endoluminal and interstitial ultrasound thermotherapy, are reviewed. Principles and practical techniques for modeling acoustic energy deposition, bioheat transfer, thermal tissue damage, and dynamic changes in the physical and physiological state of tissue are reviewed. The integration of these models and applications of simulation techniques in identification of device design parameters, development of real time feedback-control platforms, assessing the quality and safety of treatment delivery strategies, and optimisation of inverse treatment plans are presented.

Approaches for modelling interstitial ultrasound ablation of tumours within or adjacent to bone: Theoretical and experimental evaluations

Abstract Purpose: The objectives of this study were to develop numerical models of interstitial ultrasound ablation of tumours within or adjacent to bone, to evaluate model performance through theoretical analysis, and to validate the models and approximations used through comparison to experiments. Methods: 3D transient biothermal and acoustic finite element models were developed, employing four approximations of 7-MHz ultrasound propagation at bone/soft tissue interfaces. The various approximations considered or excluded reflection, refraction, angle-dependence of transmission coefficients, shear mode conversion, and volumetric heat deposition. Simulations were performed for parametric and comparative studies. Experiments within ex vivo tissues and phantoms were performed to validate the models by comparison to simulations. Temperature measurements were conducted using needle thermocouples or magnetic resonance temperature imaging (MRTI). Finite element models representing heterogeneous tissue geometries were created based on segmented MR images. Results: High ultrasound absorption at bone/soft tissue interfaces increased the volumes of target tissue that could be ablated. Models using simplified approximations produced temperature profiles closely matching both more comprehensive models and experimental results, with good agreement between 3D calculations and MRTI. The correlation coefficients between simulated and measured temperature profiles in phantoms ranged from 0.852 to 0.967 (p-value < 0.01) for the four models. Conclusions: Models using approximations of interstitial ultrasound energy deposition around bone/soft tissue interfaces produced temperature distributions in close agreement with comprehensive simulations and experimental measurements. These models may be applied to accurately predict temperatures produced by interstitial ultrasound ablation of tumours near and within bone, with applications toward treatment planning.

Multiple applicator hepatic ablation with interstitial ultrasound devices: Theoretical and experimental investigation

PURPOSE To evaluate multiple applicator implant configurations of interstitial ultrasound devices for large volume ablation of liver tumors. METHODS A 3D bioacoustic-thermal model using the finite element method was implemented to assess multiple applicator implant configurations for thermal ablation with interstitial ultrasound energy. Interstitial applicators consist of linear arrays of up to four 10 mm-long tubular ultrasound transducers, each under separate and dynamic power control, enclosed within a water-cooled delivery catheter (2.4 mm OD). The authors considered parallel implants with two and three applicators (clustered configuration), spaced 2-3 cm apart, to simulate open surgical placement. In addition, the authors considered two applicator implants with applicators converging and diverging at angles of ∼20°, 30°, and 45° to simulate percutaneous placement. Heating experiments (10-15 min) were performed and compared against simulations employing the same experimental parameters. To estimate the performance of parallel, multiple applicator configurations in an in vivo setting, simulations were performed taking into account a range of blood perfusion levels (0, 5, 12, and 15 kg m(-3) s(-1)) that may occur in tumors of varying vascularity. The impact of tailoring the power supplied to individual transducer elements along the length of applicators is explored for applicators inserted in non-parallel (converging and diverging) configurations. Thermal dose (t(43) > 240 min) and temperature thresholds (T > 52 °C) were used to define the ablation zones, with dynamic changes to tissue acoustic and thermal properties incorporated within the model. RESULTS Experiments in ex vivo bovine liver yielded ablation zones ranging between 4.0-5.6 cm × 3.2-4.9 cm, in cross section. Ablation zone dimensions predicted by simulations with similar parameters to the experiments were in close agreement (within 5 mm). Simulations of in vivo heating showed that 15 min heating and interapplicator spacing less than 3 cm are required to obtain contiguous, complete ablation zones. The ability to create complete ablation zone profiles for nonparallel implants was illustrated by tailoring applied power levels along the length of applicators. CONCLUSIONS Parallel implants consisting of three interstitial ultrasound applicators in a triangular configuration yield complete ablation zones measuring up to 6.2 cm × 5.7 cm after 15 min heating. At larger interapplicator spacing, the level of blood perfusion in the tumor may yield indentations along the periphery of the ablation zone. Tailoring applied power along the length of the applicator can accommodate for nonparallel implants, without compromising safety.

Catheter-based ultrasound technology for image-guided thermal therapy: Current technology and applications

Abstract Catheter-based ultrasound (CBUS) is applied to deliver minimally invasive thermal therapy to solid cancer tumours, benign tissue growth, vascular disease, and tissue remodelling. Compared to other energy modalities used in catheter-based surgical interventions, unique features of ultrasound result in conformable and precise energy delivery with high selectivity, fast treatment times, and larger treatment volumes. We present a concise review of CBUS technology being currently utilized in animal and clinical studies or being developed for future applications. CBUS devices have been categorised into interstitial, endoluminal and endovascular/cardiac applications. Basic applicator designs, site-specific evaluations and possible treatment applications have been discussed in brief. Particular emphasis has been given to ablation studies that incorporate image guidance for applicator placement, therapy monitoring, feedback control, and post-procedure assessment. Examples of devices included here span the entire spectrum of the development cycle from preliminary simulation-based design studies to implementation in clinical investigations. The use of CBUS under image guidance has the potential for significantly improving precision and applicability of thermal therapy delivery.