Finite element modeling of cooled-tip probe radiofrequency ablation processes in liver tissue

Investigate the relationship between pulsed field ablation parameters and ablation outcomes

BACKGROUND

Pulsed field ablation (PFA) is an emerging non-thermal ablation method. The primary challenge is the control of multiple parameters in PFA, as the interplay of these parameters remains unclear in terms of ensuring effective and safe tissue ablation.

PURPOSE

This study employs the response surface method (RSM) to explore the interactions between various PFA parameters and ablation outcomes, and seeks to enhance the efficacy and safety of PFA.

METHODS

In vivo experiments were conducted using rabbit liver for varying PFA parameters: pulse amplitude (PA), pulse interval (PI), number of pulse trains (NT), and number of pulses in a pulse train (NP). Ablation outcomes assessed included three ablation sizes, surface temperature, and muscle contraction strength. Additionally, histological analysis was performed on the ablated tissue. We analyzed the relationship between PFA parameters and ablation outcomes, and results were then compared with those from a simulation using an electric-thermal coupling PFA finite element model.

RESULTS

A linear relationship between ablation outcomes and PFA parameters was established. PA and NT exhibited extremely significant (P < 0.0001) and significant effects (P < 0.05) on all ablation outcomes, respectively. NP showed an extremely significant impact (P < 0.0001) on surface temperature and muscle contraction strength, while PI significantly influenced (P < 0.05) muscle contraction strength alone. Histological analysis revealed that PFA produces controlled, well-defined areas of liver tissue necrosis. Surface temperature results from simulations and experiments were highly consistent (R2 > 0.97).

CONCLUSIONS

This study clarifies the relationship between various PFA parameters and ablation outcomes, and aims to improve the efficacy and safety of PFA.

Analysis of hyperbolic Pennes bioheat equation in perfused homogeneous biological tissue subject to the instantaneous moving heat source

The one-dimensional hyperbolic Pennes bioheat equation under instantaneous moving heat source is solved analytically based on the Eigenvalue method. Comparison with results of in vivo experiments performed earlier by other authors shows the excellent prediction of the presented closed-form solution. We present three examples for calculating the Arrhenius equation to predict the tissue thermal damage analysis with our solution, i.e., characteristics of skin, liver, and kidney are modeled by using their thermophysical properties. Furthermore, the effects of moving velocity and perfusion rate on temperature profiles and thermal tissue damage are investigated. Results illustrate that the perfusion rate plays the cooling role in the heating source moving path. Also, increasing the moving velocity leads to a decrease in absorbed heat and temperature profiles. The closed-form analytical solution could be applied to verify the numerical heating model and optimize surgery planning parameters.

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

Percutaneous thermal ablation has proven to be an effective modality for treating both benign and malignant tumours 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°C, 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 tumour 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.

Shape-shifting thermal coagulation zone during saline-infused radiofrequency ablation: A computational study on the effects of different infusion location

BACKGROUND AND OBJECTIVE

The majority of the studies on radiofrequency ablation (RFA) have focused on enlarging the size of the coagulation zone. An aspect that is crucial but often overlooked is the shape of the coagulation zone. The shape is crucial because the majority of tumours are irregularly-shaped. In this paper, the ability to manipulate the shape of the coagulation zone following saline-infused RFA by altering the location of saline infusion is explored.

METHODS

A 3D model of the liver tissue was developed. Saline infusion was described using the dual porosity model, while RFA was described using the electrostatic and bioheat transfer equations. Three infusion locations were investigated, namely at the proximal end, the middle and the distal end of the electrode. Investigations were carried out numerically using the finite element method.

RESULTS

Results indicated that greater thermal coagulation was found in the region of tissue occupied by the saline bolus. Infusion at the middle of the electrode led to the largest coagulation volume followed by infusion at the proximal and distal ends. It was also found that the ability to delay roll-off, as commonly associated with saline-infused RFA, was true only for the case when infusion is carried out at the middle. When infused at the proximal and distal ends, the occurrence of roll-off was advanced. This may be due to the rapid and more intense heating experienced by the tissue when infusion is carried out at the electrode ends where Joule heating is dominant.

CONCLUSION

Altering the location of saline infusion can influence the shape of the coagulation zone following saline-infused RFA. The ability to 'shift' the coagulation zone to a desired location opens up great opportunities for the development of more precise saline-infused RFA treatment that targets specific regions within the tissue.

Coupled thermo-electro-mechanical models for thermal ablation of biological tissues and heat relaxation time effects

Thermal ablation is a widely applied electrosurgical process in medical treatment of soft biological tissues. Numerical modeling and simulations play an important role in prediction of temperature distribution and damage volume during the treatment planning stage of associated therapies. In this contribution we report a coupled thermo-electro-mechanical model, accounting for heat relaxation time, for more accurate and precise prediction of the temperature distribution, tissue deformation and damage volume during the thermal ablation of biological tissues. Finite element solutions are obtained for most widely used percutaneous thermal ablative techniques, viz., radiofrequency ablation (RFA) and microwave ablation (MWA). Importantly, both tissue expansion and shrinkage have been considered for modeling the tissue deformation in the coupled model of high temperature thermal ablation. The coupled model takes into account the non-Fourier effects, considering both single-phase-lag (SPL) and dual-phase-lag (DPL) models of bio-heat transfer. The temperature-dependent electrical and thermal parameters, damage-dependent blood perfusion rate and phase change effect accounting for tissue vaporization have been accounted for obtaining more clinically relevant model. The proposed model predictions are found to be in good agreement against the temperature distribution and damage volume reported by previous experimental studies. The numerical simulation results revealed that the non-Fourier effects cause a decrease in the predicted temperature distribution, tissue deformation and damage volume during the high temperature thermal ablative procedures. Furthermore, the effects of different magnitudes of phase lags of the heat flux and temperature gradient on the predicted treatment outcomes of the considered thermal ablative modalities are also quantified and discussed in detail.

Optimization of an Endoscopic Radiofrequency Ablation Electrode

Radiofrequency ablation (RFA) is an increasingly used, minimally invasive, cancer treatment modality for patients who are unwilling or unable to undergo a major resective surgery. There is a need for RFA electrodes that generate thermal ablation zones that closely match the geometry of typical tumors, especially for endoscopic ultrasound-guided (EUS) RFA. In this paper, the procedure for optimization of an RFA electrode is presented. First, a novel compliant electrode design is proposed. Next, a thermal ablation model is developed to predict the ablation zone produced by an RFA electrode in biological tissue. Then, a multi-objective genetic algorithm is used to optimize two cases of the electrode geometry to match the region of destructed tissue to a spherical tumor of a specified diameter. This optimization procedure is then applied to EUS-RFA ablation of pancreatic tissue. For a target 2.5 cm spherical tumor, the optimal design parameters of the compliant electrode design are found for two cases. Cases 1 and 2 optimal solutions filled 70.9% and 87.0% of the target volume as compared to only 25.1% for a standard straight electrode. The results of the optimization demonstrate how computational models combined with optimization can be used for systematic design of ablation electrodes. The optimization procedure may be applied to RFA of various tissue types for systematic design of electrodes for a specific target shape.

Model-based optimal planning of hepatic radiofrequency ablation

This article presents a model-based pre-treatment optimal planning framework for hepatic tumour radiofrequency (RF) ablation. Conventional hepatic radiofrequency (RF) ablation methods rely on pre-specified input voltage and treatment length based on the tumour size. Using these experimentally obtained pre-specified treatment parameters in RF ablation is not optimal to achieve the expected level of cell death and usually results in more healthy tissue damage than desired. In this study we present a pre-treatment planning framework that provides tools to control the levels of both the healthy tissue preservation and tumour cell death. Over the geometry of tumour and surrounding tissue, we formulate the RF ablation planning as a constrained optimization problem. With specific constraints over the temperature profile (TP) in pre-determined areas of the target geometry, we consider two different cost functions based on the history of the TP and Arrhenius index (AI) of the target location, respectively. We optimally compute the input voltage variation to minimize the damage to the healthy tissue while ensuring a complete cell death in the tumour and immediate area covering the tumour. As an example, we use a simulation of a 1D symmetric target geometry mimicking the application of single electrode RF probe. Results demonstrate that compared to the conventional methods both cost functions improve the healthy tissue preservation.

Performance of radiofrequency ablation used for metastatic spinal tumor: Numerical approach

Surgical treatment of spine metastases follows only local anatomical and biomechanical objectives. Few cases of actual solitary metastases are rather exceptional, while removal of these metastases and the primary tumor may help to eradicate the process. The aim of our subsequent numerical simulations was to find out the temperature distribution and the volume lesion in a spinal tumor. For this purpose, the parametric three-dimensional numerical model was developed. It was shown that by finite element modeling approach not only the temperature distribution but even the resulted cavity may be estimated. The numerical approach was shown as a strong tool in surgery planning.

Nano-assisted radiofrequency ablation of clinically extracted irregularly-shaped liver tumors

Radiofrequency ablation (RFA) for liver tumors is a minimally invasive procedure that uses electrical energy and heat to destroy cancer cells. One of the critical factors that impedes its successful outcome is the use of inappropriate radiofrequency levels that will not completely destroy the target tumor tissues, resulting in therapy failure. Additionally, the surrounding healthy tissues may suffer from serious damage due to excessive ablation. To address these challenges, this work proposes the employment of injected nanoparticles to thermally promote the ablation efficacy of conventional RFA. A three-dimensional finite difference analysis is employed to simulate the RFA treatment. Based on the data acquired from measured experiments, the simulation results have demonstrated close agreement with experimental data with a maximum discrepancy of within ±8.7%. Several types of nanoparticles were selected to evaluate their influences on liver tissue's thermal and electrical properties. We analysed the effects of nanoparticles on liver RFA via a tumor rending process incorporating several clinically-extracted tumor profiles and vascular systems. Simulations were conducted to explore the temperature difference responses between conventional RFA treatment and one with the inclusion of assisted nanoparticles on several irregularly-shaped tumors. Results have indicated that applying selected nanoparticles with high thermal conductivity and electrical conductivity on the targeted tissue zone promotes heating rate while sustaining a similar ablation zone that experiences lower maximum temperature when compared with the conventional RFA treatment. In sum, incorporating thermally-enhancing nanoparticles promotes heat transfer during the RFA treatment, resulting in improved ablation efficiency.

MRI-based finite element simulation on radiofrequency ablation of thyroid cancer

In order to provide a quantitative disclosure on the RFA (radiofrequency ablation)-induced thermal ablation effects within thyroid tissues, this paper has developed a three-dimensional finite element simulation strategy based on a MRI (magnetic resonance imaging)-reconstructed model. The thermal lesion's growth was predicted and interpreted under two treatment conditions, i.e. single-cooled-electrode modality and two-cooled-electrode system. The results show that the thermal lesion's growth is significantly affected by two factors including the position of RF electrode and thermal-physiological behavior of the breathing airflow. Additional parametric studies revealed several valuable phenomena, e.g. with the electrode's movement, thermal injury with varying severity would happen to the trachea wall. Besides, the changes in airflow mass produced evident effects on the total heat flux of thyroid surface, while the changes in breathing frequency only generated minor effects that can be ignored. The present study provided a better understanding on the thermal lesions of RFA within thyroid domain, which will help guide future treatment of the thyroid cancer.

Image-based 3D modeling and validation of radiofrequency interstitial tumor ablation using a tissue-mimicking breast phantom

PURPOSE

Minimally invasive treatment of solid cancers, especially in the breast and liver, remains clinically challenging, despite a variety of treatment modalities, including radiofrequency ablation (RFA), microwave ablation or high-intensity focused ultrasound. Each treatment modality has advantages and disadvantages, but all are limited by placement of a probe or US beam in the target tissue for tumor ablation and monitoring. The placement is difficult when the tumor is surrounded by large blood vessels or organs. Patient-specific image-based 3D modeling for thermal ablation simulation was developed to optimize treatment protocols that improve treatment efficacy.

METHODS

A tissue-mimicking breast gel phantom was used to develop an image-based 3D computer-aided design (CAD) model for the evaluation of a planned RF ablation. First, the tissue-mimicking gel was cast in a breast mold to create a 3D breast phantom, which contained a simulated solid tumor. Second, the phantom was imaged in a medical MRI scanner using a standard breast imaging MR sequence. Third, the MR images were converted into a 3D CAD model using commercial software (ScanIP, Simpleware), which was input into another commercial package (COMSOL Multiphysics) for RFA simulation and treatment planning using a finite element method (FEM). For validation of the model, the breast phantom was experimentally ablated using a commercial (RITA) RFA electrode and a bipolar needle with an electrosurgical generator (DRE ASG-300). The RFA results obtained by pre-treatment simulation were compared with actual experimental ablation.

RESULTS

A 3D CAD model, created from MR images of the complex breast phantom, was successfully integrated with an RFA electrode to perform FEM ablation simulation. The ablation volumes achieved both in the FEM simulation and the experimental test were equivalent, indicating that patient-specific models can be implemented for pre-treatment planning of solid tumor ablation.

CONCLUSION

A tissue-mimicking breast gel phantom and its MR images were used to perform FEM 3D modeling and validation by experimental thermal ablation of the tumor. Similar patient-specific models can be created from preoperative images and used to perform finite element analysis to plan radiofrequency ablation. Clinically, the method can be implemented for pre-treatment planning to predict the effect of an individual's tissue environment on the ablation process, and this may improve the therapeutic efficacy.

High-fidelity computer models for prospective treatment planning of radiofrequency ablation with in vitro experimental correlation

PURPOSE

To evaluate the accuracy of computer simulation in predicting the thermal damage region produced by a radiofrequency (RF) ablation procedure in an in vitro perfused bovine liver model. The thermal dose end point in the liver model is used to assess quantitatively computer prediction for use in prospective treatment planning of RF ablation procedures.

MATERIALS AND METHODS

Geometric details of the tri-cooled tip electrode were modeled. The resistive heating of a pulsed voltage delivery was simulated in four dimensions using finite element models (FEM) implemented on high-performance parallel computing architectures. A range of physically realistic blood perfusion parameters, 3.6-53.6 kg/sec/m(3), was considered in the computer model. An Arrhenius damage model was used to predict the thermal dose. Dice similarity coefficients (DSC) were the metric of comparison between computational predictions and T1-weighted contrast-enhanced images of the damage obtained from a RF procedure performed on an in vitro perfused bovine liver model.

RESULTS

For a perfusion parameter greater than 16.3 kg/sec/m(3), simulations predict the temporal evolution of the damaged volume is perfusion limited and will reach a maximum value. Over a range of physically meaningful perfusion values, 16.3-33.1 kg/sec/m(3), the predicted thermal dose reaches the maximum damage volume within 2 minutes of the delivery and is in good agreement (DSC > 0.7) with experimental measurements obtained from the perfused liver model.

CONCLUSIONS

As measured by the computed volumetric DSC, computer prediction accuracy of the thermal dose shows good correlation with ablation lesions measured in vitro in perfused bovine liver models over a range of physically realistic perfusion values.

Mathematical spatio-temporal model of drug delivery from low temperature sensitive liposomes during radiofrequency tumour ablation

PURPOSE

Studies have demonstrated a synergistic effect between hyperthermia and chemotherapy, and clinical trials in image-guided drug delivery combine high-temperature thermal therapy (ablation) with chemotherapy agents released in the heating zone via low temperature sensitive liposomes (LTSL). The complex interplay between heat-based cancer treatments such as thermal ablation and chemotherapy may require computational models to identify the relationship between heat exposure and pharmacokinetics in order to optimise drug delivery.

MATERIALS AND METHODS

Spatio-temporal data on tissue temperature and perfusion from heat-transfer models of radiofrequency ablation were used as input data. A spatio-temporal multi-compartmental pharmacokinetic model was built to describe the release of doxorubicin (DOX) from LTSL into the tumour plasma space, and subsequent transport into the extracellular space, and the cells. Systemic plasma and tissue compartments were also included. We compared standard chemotherapy (free-DOX) to LTSL-DOX administered as bolus at a dose of 0.7 mg/kg body weight.

RESULTS

Modelling LTSL-DOX treatment resulted in tumour tissue drug concentration of approximately 9.3 microg/g with highest values within 1 cm outside the ablation zone boundary. Free-DOX treatment produced comparably uniform tissue drug concentrations of approximately 3.0 microg/g. Administration of free-DOX resulted in a considerably higher peak level of drug concentration in the systemic plasma compartment (16.1 microg/g) compared to LTSL-DOX (4.4 microg/g). These results correlate well with a prior in vivo study.

CONCLUSIONS

Combination of LTSL-DOX with thermal ablation allows localised drug delivery with higher tumour tissue concentrations than conventional chemotherapy. Our model may facilitate drug delivery optimisation via investigation of the interplays among liposome properties, tumour perfusion, and heating regimen.

Analytical validation of COMSOL Multiphysics for theoretical models of Radiofrequency ablation including the Hyperbolic Bioheat transfer equation

In this paper we outline our main findings about the differences between the use of the Bioheat Equation and the Hyperbolic Bioheat Equation in theoretical models for Radiofrequency (RF) ablation. At the moment, we have been working on the analytical approach to solve both equations, but more recently, we have considered numerical models based on the Finite Element Method (FEM). As a first step to use FEM, we conducted a comparative study between the temperature profiles obtained from the analytical solutions and those obtained from FEM. Regarding the differences between both methods, we obtain agreement in less than 5% of relative differences. Then FEM is a good alternative to model heating of biological tissues using BE and HBE in, for example, more complex and realistic geometries.