Thermal ablation of biological tissues in disease treatment: A review of computational models and future directions

Percutaneous thermal ablation has proved to be an effective modality for treating both benign and malignant tumors in various tissues. Among these modalities, radiofrequency ablation (RFA) is the most promising and widely adopted approach that has been extensively studied in the past decades. Microwave ablation (MWA) is a newly emerging modality that is gaining rapid momentum due to its capability of inducing rapid heating and attaining larger ablation volumes, and its lesser susceptibility to the heat sink effects as compared to RFA. Although the goal of both these therapies is to attain cell death in the target tissue by virtue of heating above 50 oC, their underlying mechanism of action and principles greatly differs. Computational modelling is a powerful tool for studying the effect of electromagnetic interactions within the biological tissues and predicting the treatment outcomes during thermal ablative therapies. Such a priori estimation can assist the clinical practitioners during treatment planning with the goal of attaining successful tumor destruction and preservation of the surrounding healthy tissue and critical structures. This review provides current state-of- the-art developments and associated challenges in the computational modelling of thermal ablative techniques, viz., RFA and MWA, as well as touch upon several promising avenues in the modelling of laser ablation, nanoparticles assisted magnetic hyperthermia and non- invasive RFA. The application of RFA in pain relief has been extensively reviewed from modelling point of view. Additionally, future directions have also been provided to improve these models for their successful translation and integration into the hospital work flow

Effective treatment of solid tumors via Cryosurgery

Ph.DDOCTOR OF PHILOSOPH

Image-guidance and computational modeling to develop and characterize microwave thermal therapy platforms

Doctor of PhilosophyDepartment of Electrical and Computer EngineeringPunit PrakashThis dissertation focuses on the development of magnetic resonance imaging (MRI)-guided microwave thermal therapy systems for driving experimental studies in small animals, and to experimentally validate computational models of microwave ablation, which are widely employed for device design and characterization. MRI affords noninvasive monitoring of spatial temperature profiles, thereby providing a means to to quantitatively monitor and verify delivery of prescribed thermal doses in experimental studies and clinical use, as well as a means to validate thermal profiles predicted by computational models of thermal therapy. A contribution of this dissertation is the development and demonstration of a system for delivering mild hyperthermia to small animal targets, thereby providing a platform for driving basic research studies investigating the use of heating as part of cancer treatment strategies. An experimentally validated 3D computational model was employed to design and characterize a non-invasive directional water-cooled microwave hyperthermia applicator for MRI guided delivery of hypethermia in small animals. Following a parametric model-based design approach, a reflector aperture angle of 120°, S-shaped monopole antenna with 0.6 mm displacement, and a coolant flow rate of 150 ml/min were selected as applicator parameters that enable conformal delivery of mild hyperthermia to tumors in experimental animals. The system was integrated with real-time high-field 14.1 T MRI thermometry and feedback control to monitor and maintain target temperature elevations in the range of 4 – 5 °C (hypethermic range). 2 - 4 mm diameter targets positioned 1 – 3 mm from the applicator surface were heated to hyperthermic temperatures, with target coverage ratio ranging between 76 - 93 % and 11 – 26 % of non-targeted tissue heated. Another contribution of this dissertation is using computational models to determine how the fibroids altered ablation profile of a microwave applicator for global endometrial ablation. Uterine fibroids are benign pelvic tumors located within the myometrium or endometrium,and may alter the profile of microwave ablation applicators deployed within the uterus for delivering endometrial ablation. A 3D computational model was employed to investigate the effect of 1 – 3 cm diameter uterine fibroids in different locations around the uterine cavity on endometrial ablation profiles of microwave exposure with a 915 MHz microwave triangular loop antenna. The maximum change in simulated ablation depths due to the presence of fibroids was 1.1 mm. In summary, this simulation study suggests that 1 – 3 cm diameter uterine fibroids can be expected to have minimal impact on the extent of microwave endometrial ablation patterns achieved with the applicator studied in this dissertation. Another contribution of this dissertation is the development of a method for experimental validation of 3D transient temperature profiles predicted by computational models of MWA. An experimental platform was developed integrating custom designed MR-conditional MWA applicators for use within the MR environment. This developed platform was employed to conduct 30 - 50 W, 5 - 10 min MWA experiments in ex vivo tissue. Microwave ablation computational models, mimicking the experimental setting in MRI, were implemented using the finite element method, and incorporated temperature-dependent changes in tissue physical properties. MRI-derived Arrhenius thermal damage maps were compared to Model-predicted ablation zone extents using the Dice similarity coefficient (DSC). Mean absolute error between MR temperature measurements and fiber-optic temperature probes, used to validate the accuracy of MR temperature measurements, during heating was in the range of 0.5 – 2.8 °C. The mean DSC between model-predicted ablation zones and MRI-derived Arrhenius thermal damage maps for 13 experimental set-ups was 0.95. When comparing simulated and experimentally (i.e. using MRI) measured temperatures, the mean absolute error (MAE %) relative to maximum temperature change was in the range 5 % - 8.5 %

Investigation of Heat Therapies using Multi-Scale Models and Statistical Methods

Ph.DDOCTOR OF PHILOSOPH

Magnetic resonance imaging and navigation of ferromagnetic thermoseeds to deliver thermal ablation therapy

Minimally invasive therapies aim to deliver effective treatment whilst reducing off-target burden, limiting side effects, and shortening patient recovery times. Remote navigation of untethered devices is one method that can be used to deliver targeted treatment to deep and otherwise inaccessible locations within the body. Minimally invasive image-guided ablation (MINIMA) is a novel thermal ablation therapy for the treatment of solid tumours, whereby an untethered ferromagnetic thermoseed is navigated through tissue to a target site within the body, using the magnetic field gradients generated by a magnetic resonance imaging (MRI) system. Once at the tumour, the thermoseed is heated remotely using an alternating magnetic field, to induce cell death in the surrounding cancer tissue. The thermoseed is then navigated through the tumour, heating at pre-defined locations until the entire volume has been ablated. The aim of this PhD project is to develop MINIMA through a series of proof-of-concept studies and to assess the efficacy of the three key project components: imaging, navigation, and heating. First, an MR imaging sequence was implemented to track the thermoseeds during navigation and subsequently assessed for precision and accuracy. Secondly, movement of the thermoseeds through a viscous fluid was characterised, by measuring the effect of different navigation parameters. This was followed by navigation experiments performed in ex vivo tissue. To assess thermoseed heating, a series of in vitro experiments were conducted in air, water, and ex vivo liver tissue, before moving onto in vivo experiments in the rat brain and a murine subcutaneous tumour model. These final experiments allowed the extent of cell death induced by thermoseed heating to be determined, in both healthy and diseased tissue respectively

Computer modeling and experimentation in radiofrequiency-based minimally invasive therapies

Tesis por compendio[ES] La ablación por radiofrecuencia (RF) se ha convertido en una técnica ablativa importante, ampliamente utilizada en el área de las terapias mínimamente invasivas de la medicina moderna. El avance en el campo de las tecnologías basadas en RF a lo largo de los años ha llevado a un número creciente de aplicaciones en diferentes áreas terapéuticas tales como arritmias cardíacas, epilepsia, oncología, resección asistida, apnea, dolor o cirugía estética. Sin embargo, existe una constante necesidad de desarrollar estudios computacionales y experimentales para mejorar el rendimiento de estas técnicas. El enfoque principal de esta tesis doctoral está centrado en examinar los efectos térmicos y eléctricos de ablación por radiofrecuencia de tejidos para mejorar la eficacia y la seguridad de las terapias y dispositivos basados en energía de radiofrecuencia. Las dos áreas principales de interés han sido el tratamiento del dolor y la cirugía hepática oncológica, que se han organizado en tres estudios independientes. La metodología de los estudios se ha basado en modelos computacionales y estudios experimentales sobre phantom de agar, modelos ex vivo e in vivo y ensayos clínicos. El estudio focalizado en el tratamiento del dolor ha incluido el análisis de los efectos eléctricos y térmicos del tratamiento con radiofrecuencia pulsada (PRF) y el riesgo relacionado con el daño térmico al tejido. Se han estudiado diferentes protocolos pulsados empleados en la práctica clínica utilizando modelos computacionales. La exactitud del modelo se ha validado mediante un modelo en phantom de agar. Se han propuesto también modelos computacionales adicionales para los protocolos pulsados alternativos en los cuales se reduciría el efecto térmico sin afectar al efecto eléctrico. En el estudio se ha discutido también el concepto de electroporación leve como el resultado de PRF. En el área de la cirugía hepática oncológica se han analizado dos técnicas diferentes. El primer estudio se ha centrado en examinar la hidratación del tejido durante la ablación por RF con un nuevo electrodo ICW. El nuevo diseño ha incluido dos agujas de perfusión expandibles integradas en el catéter. El objetivo principal ha sido mejorar la precisión del modelo computacional de ablación por RF de tumor utilizando una geometría realista de la distribución de solución salina en el tejido y evaluar el rendimiento del catéter de RF. Se han modelado diferentes casos de tumor infundido con solución salina y los resultados simulados se han comparado con los datos clínicos de un ensayo en 17 pacientes con cáncer hepático. Con el fin de obtener una distribución espacial realista de la solución salina infundida, se ha empleado un estudio in vivo sobre el modelo de hígado de cerdo. El segundo estudio se ha centrado en el desarrollo de una nueva técnica de sellado endoluminal basada en catéter, como una alternativa más efectiva para el manejo del remanente pancreático. El método ha consistido en una ablación por radiofrecuencia guiada por impedancia con la técnica de pullback. El ajuste del tipo de catéter de RF y del protocolo de ablación se ha realizado mediante modelos porcinos ex vivo. Posteriormente, la efectividad del sellado se ha evaluado sobre un modelo de cerdo in vivo.[CA] L'ablació per radiofreqüència (RF) s'ha convertit en una tècnica ablativa important, àmpliament utilitzada en l'àrea de les teràpies mínimament invasives de la medicina moderna. L'avanç en el camp de les tecnologies basades en RF al llarg dels anys ha portat a un número creixent d'aplicacions en diferents àrees terapèutiques com ara arítmies cardíaques, epilèpsia, oncologia, resecció assistida, apnea, dolor o cirurgia estètica. No obstant això, hi ha una constant necessitat de desenvolupar estudis computacionals i experimentals per a millorar el rendiment d'aquestes tècniques. Aquesta tesi doctoral ha estat centrada en examinar els efectes tèrmics i elèctrics de l'ablació per radiofreqüència de teixits per tal de millorar l'eficàcia i la seguretat de les teràpies i dispositius basats en energia de radiofreqüència. Dos àrees principals són el tractament del dolor i la cirurgia hepàtica. Aquestos han sigut organitzats en tres estudis independents. La metodologia dels estudis ha estat basada en models computacionals i experimentals sobre phantom d'agar, models ex vivo i in vivo i assajos clínics. L'estudi enfocat en el tractament del dolor ha inclòs l'anàlisi dels efectes elèctrics i tèrmics del tractament amb radiofreqüència polsada (PRF) i el risc relacionat amb el dany tèrmic al teixit. S'han estudiat diferents protocols polsats emprats en la pràctica clínica utilitzant models computacionals. L'exactitud del model ha estat validada per mitjà d'un model de phantom d'agar. S'han proposat també models computacionals addicionals per a protocols polsats alternatius en els quals es reduiria l'efecte tèrmic sense afectar l'efecte elèctric. En aquest estudi s'ha discutit també el concepte d'electroporació lleu com el resultat de PRF. A l'àrea de la cirugía hepàtica han sigut analitzades dos tècniques diferents. El primer estudi s'ha centrat en la hidratació del teixit durant l'ablació per RF amb un nou elèctrode ICW. El nou disseny ha inclòs dos agulles de perfusió expandibles integrades en el catèter. L' objetiu principal ha sigut millorar la precisió del model computacional d' ablació de tumors per RF utilitzant una geometria realista per a la distribució de sèrum salií en el teixit i evaluar el rendiment del catèter de RF. S'han modelat diferents casos de tumor infundit amb sèrum salí i els resultats simulats han sigut comparats amb les dades clíniques d'un assaig dut a terme sobre 17 pacients amb càncer hepàtic. Amb l'objetiu d'obtenir una distribució espacial realista del sèrum salí injectat, s'ha du a terme un estudi in vivo basat en un model de fetge de porc. El segon estudi s'ha centrat en el desenvolupament d'una nova tècnica de tancament endoluminal bassat en catèter, com una alternativa més efectiva per a gestionar el romanent pancreàtic. El mètode ha consistit en una ablació per radiofreqüència guiada per impedància amb la tècnica de pullback. L'ajust del tipus de catèter de RF i del protocol d'ablació ha sigut realitzat per mitjà de models porcins ex vivo. Posteriorment, l'efectivitat del tancament ha sigut avaluada sobre un model de porc in vivo.[EN] Radiofrequency (RF) ablation has become an important ablative technique widely used in the area of minimally invasive therapies of the modern medicine. The advancement in the field of RF-based technologies over the years has led to a growing number of applications in different therapeutic areas such as cardiac arrhythmias, epilepsy, oncology, assisted resection, apnea, pain or aesthetic surgery. There is, however, a constant need for the development of computer and experimental studies, which would enhance the performance and safety of these techniques. The main focus of this PhD Thesis was on examining the thermal and electrical phenomena behind tissue radiofrequency ablation in order to improve the efficacy and safety of the RF-based therapies and applicators. Two main areas of interest were pain management and oncology, which were organized into three independent studies. The research methodology was based on computer modeling and experimental studies on phantoms, ex vivo and in vivo models, and clinical trials. The research on pain management involved the analysis of electrical and thermal effects of the pulsed radiofrequency (PRF) treatment and the related risk of tissue thermal damage. Different pulse protocols used in clinical practice were studied using computer modeling and the study accuracy was validated by means of agar phantom model. Additional computer models for alternative pulse protocols were also proposed, in which thermal effect would be reduced but the electrical effect would remain unchanged. The study also discussed the concept of a mild electroporation from PRF. In the area of oncology, two different techniques were analyzed. First study focused on examining tissue hydration technique during RF ablation with a novel internally cooled wet (ICW) electrode. The new design involved two expandable perfusion needles built into the catheter. The main aim was to improve the accuracy of computer model of tumor RF ablation using a realistic geometry of saline distribution in tissue, and to assess the performance of the RF catheter. Different cases of saline-infused tumor were modeled and the simulated results were compared with the clinical data from a trial on 17 hepatic cancer patients. An in vivo study on pig liver model was used to obtain a realistic spatial distribution of the infused saline. The second study focused on the development of a new catheter-based endoluminal sealing technique as more effective alternative for management of the pancreatic stump. The method consisted of the impedance-guided radiofrequency ablation with pullback. Fine-tuning involving RF catheter type and ablation protocol was performed using ex vivo porcine models, and posteriorly, sealing effectiveness was assessed on an in vivo pig model.The completion of this work would have not been possible without the financial support of the Spanish Ministerio de Economía, Industría y Competitividad that provided funding for the development of this research project, my Predoctoral scholarship, and also Travel Grant for the research stay in The Wellman Center for PhotomedicineEwertowska, E. (2019). Computer modeling and experimentation in radiofrequiency-based minimally invasive therapies [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/134057TESISCompendi

Thermal-mechanical response modelling and thermal damage prediction of soft tissues during thermal ablation

During thermal ablation, target soft tissue responses both thermally and mechanically simultaneously. However, current thermal ablation treatment mainly relies on the quantitative temperature indication to evaluate tissue behaviours and control the delivered thermal energy, which is ineffective and inaccurate. Based on these, our research study focuses on: bioheat transfer theory, linear and nonlinear elasticity of soft tissues at varied temperatures, as well as thermal damage prediction theory, and the whole program was developed in Netbeans IDE 8.1. The main contributions of our research work lie in the following aspects: Firstly, considering a situation where soft tissue’s mechanical deformation during thermal ablation is only caused by thermal loading, it is reasonable to assume that the generated strain value is within the linear range of stress-strain relationship characterisation which is also thermal stable (nearly temperature independent). Therefore, we propose our first model by integrating the heating process with thermally-induced mechanical deformations of soft tissues for simulation and analysis of the thermal ablation process. This method combines classical Fourier based bioheat transfer and constitutive elastic mechanics derived from the method of multiplicative decomposition of thermal mechanical deformation gradient, as well as non-rigid motion dynamics. The 3D governing equations are discretised spatially using finite difference scheme and temporally using implicit time integration scheme and the obtained linear system of equations are subsequently solved using a Gauss-Seidel iterative solver. Simulation implement based on proposed method can serve as a visible assistance for relevant surgeons on analysing soft tissue’s behaviours from both thermal and mechanical deformation fields rather than from just determined temperature distribution. Secondly, we present a method to characterize soft tissue thermal damage by taking into account of thermal mechanical interactions during thermal ablation, concerning stored energy by both thermal and mechanical effects can affect the energy barrier for macromolecular transitions, leading to further or the reverse damage to treated biological tissues. To do this, traditional tissue damage model of Arrhenius integration is improved by including the thermally and mechanically induced strain energy term. Simulations and comparison analysis based on different types of soft tissues are also performed to study its influences. Our findings may provide more reliable guidelines for relevant surgeons to control the tissue damage zone during thermal ablation practice. Thirdly, thermal relaxation time used to describe heating process in homogeneous substance is usually referred to as the characteristic time in non-homogeneous biological materials, which is needed to accumulate enough energy to transfer to the nearest point. Such non-Fourier thermal behaviour has also been experimentally observed in biological tissues. Our second model is presented by integrating non-Fourier bioheat transfer and constitutive elastic mechanics derived from the method of multiplicative decomposition of thermal mechanical deformation gradient, as well as non-rigid motion of dynamics to predict and analyse thermal distribution, thermal-induced mechanical deformation and tissue damage under purely thermal loads. The simulation performances are compared between two numerical methods: Finite Difference Method and Finite Element Method, from perspectives of accuracy and computing efficiency, and also against available existed experimental data and other commercialized analysis tools. Finally, our research moves on to nonlinear range characterization of tissue deformation under combined thermal and mechanical loads. Basically, the contribution of our proposed nonlinear thermal mechanical model is by extending the finite strain framework of Neo-Hookean energy function to the heating process of soft tissues during thermal ablation. Meanwhile, our nonlinear thermal mechanical model also considers the effect of collagen fibre bundles as embedded in many biological tissues. Separating free energy density modelling into isotropic and anisotropic parts, it is assumed that the anisotropy is due to the collagen fibre bundles behaviour, while the ground substance, behaves in an isotropic manner can be modelled using selected nonlinear biomaterial model. The necessary ingredients for the finite element method implementation including: weak form and time integration are also included in this chapter. Keywords: Thermal ablation, soft tissue, non-Fourier bioheat transfer, thermal mechanical deformation, anisotropic nonlinear, tissue damage