3-D modeling of the thermal coagulation necrosis induced by an interstitial ultrasonic transducer

The Agreement of the Nomogram Tool and Ultrasound Biomicroscopy Images in Calculating Ultrasound Cycloplasty Probe Model in Chinese Patients

PURPOSE

The aim of the study was to compare and explore the agreement between the nomogram tool and ultrasound biomicroscopy (UBM) images method to calculate the ultrasound cycloplasty (UCP) probe model in Chinese glaucoma patients.

METHODS

Retrospective analysis of Chinese glaucoma patients who visited Zhongshan Ophthalmic Center in Guangzhou from January to December 2019 and were eligible for UCP surgery. Visual acuity, intraocular pressure (IOP), ocular axial length (AL), and horizontal corneal diameter (white to white [WTW]) were measured. UBM images with clear ciliary body imaging and AL and WTW data were sent to trained personnel for probe model measurements. The data calculated by both methods were analyzed using unweighted and weighted κ statistics. The level of agreement refers to Landis and Koch's guideline for the strength of agreement indicated with weighted κ values.

RESULTS

1,061 eyes of 642 patients were involved, with a mean age of 61.66 ± 11.66 years. Their best-corrected visual acuity converted to logarithm of minimal-angle-of-resolution (logMAR) scores of -0.18-3.00 with a mean value of 0.69 ± 0.77. IOP was 22.0-60.0 mm Hg with a mean of 27.97 ± 5.66 mm Hg. The mean AL and WTW were 22.88 ± 1.33 (19.15-32.14) mm and 11.52 ± 0.49 (10.00-12.90) mm, respectively. The agreement between the two methods was fair (weighted κ = 0.299), matching in 62.86% of eyes (weighted κ = 0.299, κ = 0.264). The agreement in primary open angle glaucoma, acute primary angle-closure glaucoma, chronic primary angle-closure glaucoma, and secondary glaucoma patients was 60.85% (weighted κ = 0.336, κ = 0.301), 65.06% (weighted κ = 0.146, κ = 0.127), 62.26% (weighted κ = 0.204, κ = 0.184), and 57.97% (weighted κ = 0.332, κ = 0.280) of eyes, respectively.

CONCLUSION

The agreement between UBM images and the nomogram tool to calculate the UCP probe model of Chinese patients is at a fair level. The nomogram tool prefers to use larger probes. Improvements to the nomogram tool, such as including data from more ethnic groups and being able to calculate separately for different types of glaucoma, are needed to improve accuracy. The inclusion of parameters or images from more directions of the eye may help measure probe models more accurately for both the nomogram tool and the UBM image measurement.

Development of a Numerical Model of High-Intensity Focused Ultrasound Treatment in Mobile and Elastic Organs: Application to a Beating Heart

High-intensity focused ultrasound (HIFU) is a promising method used to treat cardiac arrhythmias, as it can induce lesions at a distance throughout myocardium thickness. Numerical modeling is commonly used for ultrasound probe development and optimization of HIFU treatment strategies. This study was aimed at describing a numerical method to simulate HIFU thermal ablation in elastic and mobile heart models. The ultrasound pressure field is computed on a 3-D orthonormal grid using the Rayleigh integral method, and the attenuation is calculated step by step between cells. The temperature distribution is obtained by resolution of the bioheat transfer equation on a 3-D non-orthogonally structured curvilinear grid using the finite-volume method. The simulation method is applied on two regions of the heart (atrioventricular node and ventricular apex) to compare the thermal effects of HIFU ablation depending on deformation, motion type and amplitude. The atrioventricular node requires longer sonication than the ventricular apex to reach the same lesion volume. Motion considerably influences treatment duration, lesion shape and distribution in cardiac HIFU treatment. These results emphasize the importance of considering local motion and deformation in numerical studies to define efficient and accurate treatment strategies.

A three-dimensional visualization preoperative treatment planning system for microwave ablation in liver cancer: a simulated experimental study

PURPOSE

To evaluate the application value of three-dimensional (3D) visualization preoperative treatment planning system (VPTPS) for microwave ablation (MWA) in liver cancer.

MATERIALS AND METHODS

The study was a simulated experimental study using the CT imaging data of patients in DICOM format in a model. Three students (who learn to interventional ultrasound for less than 1 year) and three experts (who have more than 5 years of experience in ablation techniques) in MWA performed the preoperative planning for 39 lesions (mean diameter 3.75 ± 1.73 cm) of 32 patients using two-dimensional (2D) image planning method and 3D VPTPS, respectively. The number of planning insertions, planning ablation rate, and damage rate to surrounding structures were compared between2D image planning group and 3D VPTPS group.

RESULTS

There were fewer planning insertions, lower ablation rate and higher damage rate to surrounding structures in 2D image planning group than 3D VPTPS group for both students and experts. When using the 2D ultrasound planning method, students could carry out fewer planning insertions and had a lower ablation rate than the experts (p < 0.001). However, there was no significant difference in planning insertions, the ablation rate, and the incidence of damage to the surrounding structures between students and experts using 3D VPTPS.

CONCLUSIONS

3DVPTPS enables inexperienced physicians to have similar preoperative planning results to experts, and enhances students' preoperative planning capacity, which may improve the therapeutic efficacy and reduce the complication of MWA.

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

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.

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.

Phase estimation for a phased array therapeutic interstitial ultrasound probe

This paper deals about high intensity ultrasound interstitial therapy simulation. The simulated phased array ultrasound probe allows a dynamic electronic focusing of the therapeutic beam. In order to maximize the power deposit at the focal point we propose a method which allows to optimally defining the phase shift of the electrical control signal for each individual element.

Considerations for theoretical modelling of thermal ablation with catheter-based ultrasonic sources: implications for treatment planning, monitoring and control

PURPOSE

To determine the impact of including dynamic changes in tissue physical properties during heating on feedback controlled thermal ablation with catheter-based ultrasound. Additionally, we compared the impact of several indicators of thermal damage on predicted extents of ablation zones for planning and monitoring ablations with this modality.

METHODS

A 3D model of ultrasound ablation with interstitial and transurethral applicators incorporating temperature-based feedback control was used to simulate thermal ablations in prostate and liver tissue. We investigated five coupled models of heat dependent changes in tissue acoustic attenuation/absorption and blood perfusion of varying degrees of complexity. Dimensions of the ablation zone were computed using temperature, thermal dose, and Arrhenius thermal damage indicators of coagulative necrosis. A comparison of the predictions by each of these models was illustrated on a patient-specific anatomy in the treatment planning setting.

RESULTS

Models including dynamic changes in blood perfusion and acoustic attenuation as a function of thermal dose/damage predicted near-identical ablation zone volumes (maximum variation < 2.5%). Accounting for dynamic acoustic attenuation appeared to play a critical role in estimating ablation zone size, as models using constant values for acoustic attenuation predicted ablation zone volumes up to 50% larger or 47% smaller in liver and prostate tissue, respectively. Thermal dose (t(43) ≥ 240 min) and thermal damage (Ω ≥ 4.6) thresholds for coagulative necrosis are in good agreement for all heating durations, temperature thresholds in the range of 54°C for short (<5 min) duration ablations and 50°C for long (15 min) ablations may serve as surrogates for determination of the outer treatment boundary.

CONCLUSIONS

Accounting for dynamic changes in acoustic attenuation/absorption appeared to play a critical role in predicted extents of ablation zones. For typical 5-15 min ablations with this modality, thermal dose and Arrhenius damage measures of ablation zone dimensions are in good agreement, while appropriately selected temperature thresholds provide a computationally cheaper surrogate.

Dual-focus therapeutic ultrasound transducer for production of broad tissue lesions

In noninvasive high-intensity focused ultrasound (HIFU) treatment, formation of a large tissue lesion per sonication is desirable for reducing the overall treatment time. The goal of this study is to show the feasibility of enlarging tissue lesion size with a dual-focus therapeutic ultrasound transducer (DFTUT) by increasing the depth-of-focus (DOF). The proposed transducer consists of a disc- and an annular-type element of different radii of curvatures to produce two focal zones. To increase focal depth and to maintain uniform beamwidth of the elongated DOF, each element transmits ultrasound of a different center frequency: the inner element at a higher frequency for near field focusing and the outer element at a lower frequency for far field focusing. By activating two elements at the same time with a single transmitter capable of generating a dual-frequency mixed signal, the overall DOF of the proposed transducer may be extended considerably. A prototype transducer composed of a 4.1 MHz inner element and a 2.7 MHz outer element was fabricated to obtain preliminary experimental results. The feasibility the proposed technique was demonstrated through sound field, temperature and thermal dose simulations. The performance of the prototype transducer was verified by hydrophone measurements and tissue ablation experiments on a beef liver specimen. When several factors affecting the length and the uniformity of elongated DOF of the DFTUT are optimized, the proposed therapeutic ultrasound transducer design may increase the size of ablated tissues in the axial direction and, thus, decreasing the treatment time for a large volume of malignant tissues especially deep-seated targets.

Fast FFT-based bioheat transfer equation computation

This paper describes a modeling method of the tissue temperature evolution over time in hyper or hypothermia. The tissue temperature evolution over time is classically described by Pennes' bioheat transfer equation which is generally solved by a finite difference method. In this paper we will present a method where the bioheat transfer equation can be algebraically solved after a Fourier transformation over the space coordinates. As an example, we implemented this method for the simulation of a percutaneous high intensity ultrasound hepatocellular carcinoma curative treatment and compared it with the finite difference method and experimental data.

Interstitial thermal ablation with a fast rotating dual-mode transducer

Interstitial ultrasound applicators can be a minimally invasive alternative for treating targets that are unresectable or are inaccessible by extracorporeal methods. Dual-mode transducers for ultrasound imaging and therapy were developed to address the constraints of a miniaturized applicator and real-time treatment monitoring. We propose an original treatment strategy that combines ultrasound imaging and therapy using a dual-mode transducer rotating at 8 revolutions per second. Real-time B-mode imaging was interrupted to emit high-intensity ultrasound over a selected therapy aperture. A full 360 degrees image was taken every 8th rotation to image the therapy aperture. Numerical simulations were performed to study the effect of rotation on tissue heating, and to study the effect of the treatment sequence on transducer temperature. With the time-averaged transducer surface intensity held at 12 W/cm(2) to maintain transducer temperature below 66 degrees C, higher field intensities and deeper lesions were produced by narrower therapy apertures. A prototype system was built and tested using in vitro samples of porcine liver. Lesions up to 8 mm were produced using a time-averaged transducer surface intensity of 12 W/cm(2) applied for a period of 240 s over a therapy aperture of 40 degrees. Apparent strain imaging of the therapy aperture improved the contrast between treated and spared tissues, which could not be differentiated on B-mode images. With appropriate limits on the transducer output, real-time imaging and deep thermal ablation are feasible and sustainable using a rotating dual-mode transducer.

In vivo evaluation of a mechanically oscillating dual-mode applicator for ultrasound imaging and thermal ablation

Unresectable liver tumors are often treated with interstitial probes that modify tissue temperature, and efficacious treatment relies on image guidance for tissue targeting and assessment. Here, we report the in vivo evaluation of an interstitial applicator with a mechanically oscillating five-element dual-mode transducer. After thoroughly characterizing the transducer, tissue response to high-intensity ultrasound was numerically calculated to select parameters for experimentation in vivo. Using perfused porcine liver, B-mode sector images were formed before and after a 120-s therapy period, and M-mode imaging monitored the therapy axis during therapy. The time-averaged transducer surface intensity was 21 or 27 W/cm (2). Electroacoustic conversion efficiency was maximally 72 +/- 3% and impulse response length was 295 +/- 1.0 ns at -6 dB. The depth of thermal damage measured by gross histology ranged from 10 to 25 mm for 13 insertion sites. For six sites, M-mode data exhibited a reduction in gray-scale intensity that was interpreted as the temporal variation of coagulation necrosis. Contrast ratio analysis indicated that the gray-scale intensity dropped by 7.8 +/- 3.3 dB, and estimated the final lesion depth to an accuracy of 2.3 +/- 2.4 mm. This paper verified that the applicator could induce coagulation necrosis in perfused liver and demonstrated the feasibility of real-time monitoring.

FFT-based computation of the bioheat transfer equation for the HCC ultrasound surgery therapy modeling

This paper describes a modeling method of the tissue temperature evolution over time in hyperthermia. More precisely, this approach is used to simulate the hepatocellular carcinoma curative treatment by a percutaneous high intensity ultrasound surgery. The tissue temperature evolution over time is classically described by Pennes' bioheat transfer equation which is generally solved by a finite difference method. In this paper we will present a method where the bioheat transfer equation can be algebraically solved after a Fourier transformation over the space coordinates. The implementation and boundary conditions of this method will be shown and compared with the finite difference method.

Preoperative surgery planning for percutaneous hepatic microwave ablation

A novel preoperative surgery planning method is proposed for percutaneous hepatic microwave ablation. An iterative framework for necrosis field simulation and 3D necrosis zone reconstruction is introduced here, and the necrosis model is further superimposed to patient anatomy structures using advanced GPU-accelerated visualization techniques. The full surgery planning is performed by the surgeon in an interactively way, until the optimal surgery plan is achieved. Experiments have been performed on realistic patient with hepatic cancer and the actual necrosis zone are measured in postoperative CT images for patient. Results show that this method is relative accurate for preoperative trajectory plan and could be used as an assistant to the clinical practice.