Temperature-controlled radiofrequency ablation of different tissues using two-compartment models

Non-Fourier Bioheat Transfer Analysis in Brain Tissue During Interstitial Laser Ablation: Analysis of Multiple Influential Factors

This work presents the dual-phase lag-based non-Fourier bioheat transfer model of brain tissue subjected to interstitial laser ablation. The finite element method has been utilized to predict the brain tissue's temperature distributions and ablation volumes. A sensitivity analysis has been conducted to quantify the effect of variations in the input laser power, treatment time, laser fiber diameter, laser wavelength, and non-Fourier phase lags. Notably, in this work, the temperature-dependent thermal properties of brain tissue have been considered. The developed model has been validated by comparing the temperature obtained from the numerical and ex vivo brain tissue during interstitial laser ablation. The ex vivo brain model has been further extended to in vivo settings by incorporating the blood perfusion effects. The results of the systematic analysis highlight the importance of considering temperature-dependent thermal properties of the brain tissue, non-Fourier behavior, and microvascular perfusion effects in the computational models for accurate predictions of the treatment outcomes during interstitial laser ablation, thereby minimizing the damage to surrounding healthy tissue. The developed model and parametric analysis reported in this study would assist in a more accurate and precise prediction of the temperature distribution, thus allowing to optimize the thermal dosage during laser therapy in the brain.

Sublethal thermal stress promotes migration and invasion of thyroid cancer cells

OBJECTIVE

Radiofrequency ablation is a viable option in the treatment of benign thyroid nodules. Some reports suggest that thermal ablation may also be safe for the management of low-risk thyroid cancer. In this study, we applied transient heat treatment to thyroid cancer cells to mimic clinical scenarios in which insufficient ablation leads to incomplete eradication of thyroid cancer.

METHODS

Differentiated thyroid cancer cell lines B-CPAP, TPC-1, and FTC-133 were subjected to heat treatment at different temperatures for 10 min. Effects on cell growth, clonogenicity, wound healing assay, and Transwell invasion were determined.

RESULTS

Heat treatment at 45°C or higher reduced cell growth, whereas viability of thyroid cancer cells was not changed after heat treatment at 37, 40, or 42°C. Heat treatment at 40°C increased the number of colony formations by 16% to 39%. Additionally, transient heat treatment at 40°C resulted in a 1.75-fold to 2.56-fold higher migratory activity than treatment at 37°C. Invasive capacity was increased after heat treatment, ranging from 115% to 126%. Expression of several epithelial-mesenchymal transition markers, including ZEB1, N-cadherin, and MMP2, was upregulated following heat treatment at 40°C.

CONCLUSION

We for the first time demonstrate that sublethal thermal stress may increase clonogenicity, migration, and invasion of thyroid cancer cells.

Predicting spatio-temporal radiofrequency ablation temperature using deep neural networks

Radiofrequency ablation (RFA) of the medial branch nerve is a widely used therapeutic intervention for facet joint pain. However, denervation of the multifidus muscle is an inevitable consequence of RFA. New ablation techniques with the potential to prevent muscle denervation can be designed using computational simulations. However, depending on the complexity of the model, they could be computationally expensive. As an alternative approach, deep neural networks (DNNs) can be used to predict tissue temperature during RFA procedure. The objective of this paper is to predict the tissue spatial and temporal temperature distributions during RFA using DNNs. First, finite element (FE) models with a range of distances between the probes were run to obtain the temperature readings. The measured temperatures were then used to train the DNNs that predict the spatio-temporal temperature distribution within the tissue. Finally, a separate data obtained from FE simulations were used to test the efficacy of the network. The results presented in this paper demonstrate that the network can achieve an error rate as low as 0.05%, accompanied by a 92% reduction in time compared to FE simulations. The approach proposed in this study will play a major role in the design of new RFA treatments for facet joint pain.

Enlarging the thermal coagulation volume during thermochemical ablation with alternating acid-base injection by shortening the injection interval: A computational study

BACKGROUND AND OBJECTIVES

Thermochemical ablation (TCA) is a cancer treatment that utilises the heat released from the neutralisation of acid and base to raise tissue temperature to levels sufficient to induce thermal coagulation. Computational studies have demonstrated that the coagulation volume produced by sequential injection is smaller than that with simultaneous injection. By injecting the reagents in an ensuing manner, the region of contact between acid and base is limited to a thin contact layer sandwiched between the distribution of acid and base. It is hypothesised that increasing the frequency of acid-base injections into the tissue by shortening the injection interval for each reagent can increase the effective area of contact between acid and base, thereby intensifying neutralisation and the exothermic heat released into the tissue.

METHODS

To verify this hypothesis, a computational model was developed to simulate the thermochemical processes involved during TCA with sequential injection. Four major processes that take place during TCA were considered, i.e., the flow of acid and base, their neutralisation, the release of exothermic heat and the formation of thermal damage inside the tissue. Equimolar acid and base at 7.5 M was injected into the tissue intermittently. Six injection intervals, namely 3, 6, 15, 20, 30 and 60 s were investigated.

RESULTS

Shortening of the injection interval led to the enlargement of coagulation volume. If one considers only the coagulation volume as the determining factor, then a 15 s injection interval was found to be optimum. Conversely, if one places priority on safety, then a 3 s injection interval would result in the lowest amount of reagent residue inside the tissue after treatment. With a 3 s injection interval, the coagulation volume was found to be larger than that of simultaneous injection with the same treatment parameters. Not only that, the volume also surpassed that of radiofrequency ablation (RFA); a conventional thermal ablation technique commonly used for liver cancer treatment.

CONCLUSION

The numerical results verified the hypothesis that shortening the injection interval will lead to the formation of larger thermal coagulation zone during TCA with sequential injection. More importantly, a 3 s injection interval was found to be optimum for both efficacy (large coagulation volume) and safety (least amount of reagent residue).

Model-Optimizing Radiofrequency Parameters of 3D Finite Element Analysis for Ablation of Benign Thyroid Nodules

Radiofrequency (RF) ablation represents an efficient strategy to reduce the volume of thyroid nodules. In this study, a finite element model was developed with the aim of optimizing RF parameters, e.g., input power and treatment duration, in order to achieve the target volume reduction rate (VRR) for a thyroid nodule. RF ablation is modelled as a coupled electro-thermal problem wherein the electric field is applied to induce tissue heating. The electric problem is solved with the Laplace equation, the temperature distribution is estimated with the Pennes bioheat equation, and the thermal damage is evaluated using the Arrhenius equation. The optimization model is applied to RF electrode with different active tip lengths in the interval from 5 mm to 40 mm at the 5 mm step. For each case, we also explored the influence of tumour blood perfusion rate on RF ablation outcomes. The model highlights that longer active tips are more efficient as they require lesser power and shorter treatment time to reach the target VRR. Moreover, this condition is characterized by a reduced transversal ablation zone. In addition, a higher blood perfusion increases the heat dispersion, requiring a different combination of RF power and time treatment to achieve the target VRR. The model may contribute to an improvement in patient-specific RF ablation treatment.

Spacer-Supported Thermal Ablation to Prevent Carbonisation and Improve Ablation Size: A Proof of Concept Study

Thermal ablation offers a minimally invasive alternative in the treatment of hepatic tumours. Several types of ablation are utilised with different methods and indications. However, to this day, ablation size remains limited due to the formation of a central non-conductive boundary layer. In thermal ablation, this boundary layer is formed by carbonisation. Our goal was to prevent or delay carbonisation, and subsequently increase ablation size. We used bovine liver to compare ablation diameter and volume, created by a stand-alone laser applicator, with those created when utilising a spacer between laser applicator and hepatic tissue. Two spacer variants were developed: one with a closed circulation of cooling fluid and one with an open circulation into hepatic tissue. We found that the presence of a spacer significantly increased ablation volume up to 75.3 cm³, an increase of a factor of 3.19 (closed spacer) and 3.02 (open spacer) when compared to the stand-alone applicator. Statistical significance between spacer variants was also present, with the closed spacer producing a significantly larger ablation volume (p < 0.001, M_Diff = 3.053, 95% CI[1.612, 4.493]) and diameter (p < 0.001, M_Diff = 4.467, 95% CI[2.648, 6.285]) than the open spacer. We conclude that the presence of a spacer has the potential to increase ablation size.

Three-Phase-Lag Bio-Heat Transfer Model of Cardiac Ablation

Significant research efforts have been devoted in the past decades to accurately modelling the complex heat transfer phenomena within biological tissues. These modeling efforts and analysis have assisted in a better understanding of the intricacies of associated biological phenomena and factors that affect the treatment outcomes of hyperthermic therapeutic procedures. In this contribution, we report a three-dimensional non-Fourier bio-heat transfer model of cardiac ablation that accounts for the three-phase-lags (TPL) in the heat propagation, viz., lags due to heat flux, temperature gradient, and thermal displacement gradient. Finite element-based COMSOL Multiphysics software has been utilized to predict the temperature distributions and ablation volumes. A comparative analysis has been conducted to report the variation in the treatment outcomes of cardiac ablation considering different bio-heat transfer models. The effect of variations in the magnitude of different phase lags has been systematically investigated. The fidelity and integrity of the developed model have been evaluated by comparing the results of the developed model with the analytical results of the recent studies available in the literature. This study demonstrates the importance of considering non-Fourier lags within biological tissue for predicting more accurately the characteristics important for the efficient application of thermal therapies.

Depth of thermal dispersion of monopolar radiofrequency heating in the vaginal wall

The use of energy-based devices to treat genitourinary syndrome of menopause, termed vaginal thermotherapy (VTT), has gained significant interest in recent years. Among the primary safety concerns of this relatively new procedure is the possibility of unintentionally heating tissues adjacent to the vaginal wall, i.e., heating too deeply. Herein we use numerical simulations to evaluate monopolar radiofrequency-based (RF) VTT specifically focusing on the resultant depth of heating through a range of input parameters. Varying RF power, exposure time, and the simulated rate of blood perfusion, we map the parameter space identifying which combinations of input parameters are likely to heat past the depth of the vaginal wall and affect adjacent tissue. We found that the device parameters commonly used in the literature are likely to heat past the vaginal wall and merit further investigation. In addition, we found that the parameter typically used to describe VTT devices, total energy delivered, does not reliably indicate the resultant depth of heat dispersion.

Unidirectional ablation minimizes unwanted thermal damage and promotes better thermal ablation efficacy in time-based switching bipolar radiofrequency ablation

Switching bipolar radiofrequency ablation (bRFA) is a thermal treatment modality used for liver cancer treatment that is capable of producing larger, more confluent and more regular thermal coagulation. When implemented in the no-touch mode, switching bRFA can prevent tumour track seeding; a medical phenomenon defined by the deposition of cancer cells along the insertion track. Nevertheless, the no-touch mode was found to yield significant unwanted thermal damage as a result of the electrodes' position outside the tumour. It is postulated that the unwanted thermal damage can be minimized if ablation can be directed such that it focuses only within the tumour domain. As it turns out, this can be achieved by partially insulating the active tip of the RF electrodes such that electric current flows in and out of the tissue only through the non-insulated section of the electrode. This concept is known as unidirectional ablation and has been shown to produce the desired effect in monopolar RFA. In this paper, computational models based on a well-established mathematical framework for modelling RFA was developed to investigate if unidirectional ablation can minimize unwanted thermal damage during time-based switching bRFA. From the numerical results, unidirectional ablation was shown to produce treatment efficacy of nearly 100%, while at the same time, minimizing the amount of unwanted thermal damage. Nevertheless, this effect was observed only when the switch interval of the time-based protocol was set to 50 s. An extended switch interval negated the benefits of unidirectional ablation.

Finite element modeling of Non-Fourier heat transfer in a cancerous tissue with an injected fat layer during hyperthermia treatment

Hyperthermia technique has received much attention over the last decade being less invasive among the others. Laser therapy is among the most commonly investigated types of ablative hyperthermia for treatment of cancer. In this method an external heat source provided by a laser fiber leads the cancerous tissue to the necrosis stage. For its simulation a cylindrical geometry of a breast tissue containing a tumor is acted upon by a Gaussian form of laser radiation. Then the feasibility of a fat layer injection around the tumor during the therapy is investigated numerically. In order to consider the finite speed of heat transfer, dual phase lag (DPL) model is implemented for prediction of the thermal results. The therapy is addressed with and without the presence of a fat layer around the breast tumor. Results show that the temperature in the tumor increases up to 15 % by the injection of a fat layer. Also, the presence of a fat layer around the tumor shows that the irreversible ablation happens at a faster rate.

Comparisons between impedance-based and time-based switching bipolar radiofrequency ablation for the treatment of liver cancer

Switching bipolar radiofrequency ablation (bRFA) is a cancer treatment technique that activates multiple pairs of electrodes alternately based on a predefined criterion. Various criteria can be used to trigger the switch, such as time (ablation duration) and tissue impedance. In a recent study on time-based switching bRFA, it was determined that a shorter switch interval could produce better treatment outcome than when a longer switch interval was used, which reduces tissue charring and roll-off induced cooling. In this study, it was hypothesized that a more efficacious bRFA treatment can be attained by employing impedance-based switching. This is because ablation per pair can be maximized since there will be no interruption to RF energy delivery until roll-off occurs. This was investigated using a two-compartment 3D computational model. Results showed that impedance-based switching bRFA outperformed time-based switching when the switch interval of the latter is 100 s or higher. When compared to the time-based switching with switch interval of 50 s, the impedance-based model is inferior. It remains to be investigated whether the impedance-based protocol is better than the time-based protocol for a switch interval of 50 s due to the inverse relationship between ablation and treatment efficacies. It was suggested that the choice of impedance-based or time-based switching could ultimately be patient-dependent.

Bipolar radiofrequency ablation treatment of liver cancer employing monopolar needles: A comprehensive investigation on the efficacy of time-based switching

Radiofrequency ablation (RFA) is a thermal ablative treatment method that is commonly used to treat liver cancer. However, the thermal coagulation zone generated using the conventional RFA system can only successfully treat tumours up to 3 cm in diameter. Switching bipolar RFA has been proposed as a way to increase the thermal coagulation zone. Presently, the understanding of the underlying thermal processes that takes place during switching bipolar RFA remains limited. Hence, the objective of this study is to provide a comprehensive understanding on the thermal ablative effects of time-based switching bipolar RFA on liver tissue. Five switch intervals, namely 50, 100, 150, 200 and 300 s were investigated using a two-compartment 3D finite element model. The study was performed using two pairs of RF electrodes in a four-probe configuration, where the electrodes were alternated based on their respective switch interval. The physics employed in the present study were verified against experimental data from the literature. Results obtained show that using a shorter switch interval can improve the homogeneity of temperature distribution within the tissue and increase the rate of temperature rise by delaying the occurrence of roll-off. The coagulation volume obtained was the largest using switch interval of 50 s, followed by 100, 150, 200 and 300 s. The present study demonstrated that the transient thermal response of switching bipolar RFA can be improved by using shorter switch intervals.

Pre-operative Assessment of Ablation Margins for Variable Blood Perfusion Metrics in a Magnetic Resonance Imaging Based Complex Breast Tumour Anatomy: Simulation Paradigms in Thermal Therapies

BACKGROUND AND OBJECTIVES

Image-guided medical interventions facilitates precise visualization at treatment site. The conformal prediction for sparing healthy tissue fringes precisely in the vicinity of irregular tumour anatomy remains clinically challenging. Pre-clinical image-based computational modelling is imperative as it helps in enhancement of treatment quality, augmenting clinical-decision making, while planning, targeting, controlling, monitoring and assessing treatment response with an effective risk assessment before the onset of treatment in clinical settings. In this study, the influence of heat deposition rate (SAR), exposure duration, and variable blood perfusion metrics for a patient-specific breast tumour is quantified considering the tumour margins thereby suggesting need of geometrically accurate models.

METHODS

A three-dimensional realistic model mimicking dimensions of a female breast, comprising ~1.7 cm irregular tumour, was generated from patient specific two-dimensional DICOM format MRI images through image segmentation tools MIMICS 19.0® and 3-Matic 11.0® which is finally exported to COMSOL Multiphysics 5.2® as a volumetric mesh for finite element analysis. The Pennes bioheat transfer model and Arrhenius thermal damage model of cell-death are integrated to simulate a coupled biophysics problem. A comparative blood perfusion analysis is done to evaluate the response of tumour during heating considering thermal damage extent, including the tumour margins while sparing critical adjoining healthy tissues.

RESULTS

The evaluated thermal damage zones for 1 mm, 2 mm and 3 mm fringe heating region (beyond tumour boundary) reveals 0.09%, 0.21% and 0.34% thermal damage to the healthy tissue (which is <1%) and thus successful necrosis of the tumour. The iterative computational experiments suggests treatment margins < 5 mm are sufficient enough as heating beyond 3 mm fringe layer leads to higher damage surrounding the tumour approximately 1.5 times the tumour volume. Further, the heat-dosage requirements are 22% more for highly perfused tumour as compared to moderately perfused tumour with an approximate double time to ablate the whole tumour volume.

CONCLUSIONS

Depending on the blood perfusion characteristics of a tumour, it is a trade-off between heat-dosage (SAR) and exposure/treatment duration to get desired thermal damage including the irregular tumour boundaries while taking into account, the margin of healthy tissue. The suggested patient-specific integrated multiphysics-model based on MRI-Images may be implemented for pre-treatment planning based on the tumour blood perfusion to evaluate the thermal ablation zone dimensions clinically and thereby avoiding the damage of off-target tissues. Thus, risks involving underestimation or overestimation of thermal coagulation zones may be minimised while preserving the surrounding normal breast parenchyma.

Computational Modeling of Cardiac Ablation Incorporating Electrothermomechanical Interactions

The application of radio frequency ablation (RFA) has been widely explored in treating various types of cardiac arrhythmias. Computational modeling provides a safe and viable alternative to ex vivo and in vivo experimental studies for quantifying the effects of different variables efficiently and reliably, apart from providing a priori estimates of the ablation volume attained during cardiac ablation procedures. In this contribution, we report a fully coupled electrothermomechanical model for a more accurate prediction of the treatment outcomes during the radio frequency cardiac ablation. A numerical model comprising of cardiac tissue and the cardiac chamber has been developed in which an electrode has been inserted perpendicular to the cardiac tissue to simulate actual clinical procedures. Temperature-dependent heat capacity, electrical and thermal conductivities, and blood perfusion rate have been considered to model more realistic scenarios. The effects of blood flow and contact force of the electrode tip on the treatment outcomes of a fully coupled model of RFA have been systematically investigated. The numerical study demonstrates that the predicted ablation volume of RFA is significantly dependent on the blood flow rate in the cardiac chamber and also on the tissue deformation induced due to electrode insertion depth of 1.5 mm or higher.

Thermal and thermal damage responses during switching bipolar radiofrequency ablation employing bipolar needles: A computational study on the effects of different electrode configuration, input voltage and ablation duration

Recent studies have demonstrated the effectiveness of switching bipolar radiofrequency ablation (bRFA) in treating liver cancer. Nevertheless, the clinical use of the treatment remains less common than conventional monopolar RFA - likely due to the lack of understanding of how the tissues respond thermally to the switching effect. The problem is exacerbated by the numerous possible switching combinations when bRFA is performed using bipolar needles, thus making theoretical deduction and experimental studies difficult. This article addresses this issue via computational modelling by examining if significant variation in the treatment outcome exists amongst six different electrode configurations defined by the X-, C-, U-, N-, Z- and O-models. Results indicated that the tissue thermal and thermal damage responses varied depending on the electrode configuration and the operating conditions (input voltage and ablation duration). For a spherical tumour, 30 mm in diameter, complete ablation could not be attained in all configurations with 70 V input voltage and 5 minutes ablation duration. Increasing the input voltage to 90 V enlarged the coagulation zone in the X-model only. With the other configurations, extending the ablation duration to 10 minutes was found to be the better at enlarging the coagulation zone.

Optimal Localization of a Novel Shifted 1T-Ring Based Microwave Ablation Probe in Hepatocellular Carcinoma

OBJECTIVE

Hepatocellular carcinoma (HCC) is one of the most dangerous, and fatal cancers. Thermal ablation proved its power as the best treatment method for HCC. In microwave thermal ablation, microwave probes are used to generate electromagnetic waves (EMW) at microwave (MW) frequency 2.45 GHz. In this paper, the design/model of a novel microwave ablation probe, namely a single slot with a shifted 1T-ring probe is presented for HCC therapy.

METHODS

A Finite Element Method (FEM) is employed to model the probe and the hepatic tumor liver tissues. The relation between the tip of probe position and the center of the hepatic tumor was studied to determine the best probe location at which a minimum MW power is required to ablate the entire tumor tissues with the smallest damage in the nearby healthy tissues to the tumor.

RESULTS

The results indicated that the ablated part of the tissues varies depending on the MW probe type, the amount of used power, the location of the probe, and the exposure time. Hepatic tumors' diameters from 2-5cm were studied.

CONCLUSION

It was shown that the proposed SSS 1T-ring (single slot with shifted 1T-ring) probe provided the best ablation performance when the probe's tip placed below the tumor's center by 11 mm, which achieved 100% damage in the tumor tissues using 6 W power for 10 minutes.

SIGNIFICANCE

When the probe's tip is located at the center of the tumor, the ablation rate was 73.45% in the tumor tissues under the same conditions.

Domain Heterogeneity in Radiofrequency Therapies for Pain Relief: A Computational Study with Coupled Models

The objective of the current research work is to study the differences between the predicted ablation volume in homogeneous and heterogeneous models of typical radiofrequency (RF) procedures for pain relief. A three-dimensional computational domain comprising of the realistic anatomy of the target tissue was considered in the present study. A comparative analysis was conducted for three different scenarios: (a) a completely homogeneous domain comprising of only muscle tissue, (b) a heterogeneous domain comprising of nerve and muscle tissues, and (c) a heterogeneous domain comprising of bone, nerve and muscle tissues. Finite-element-based simulations were performed to compute the temperature and electrical field distribution during conventional RF procedures for treating pain, and exemplified here for the continuous case. The predicted results reveal that the consideration of heterogeneity within the computational domain results in distorted electric field distribution and leads to a significant reduction in the attained ablation volume during the continuous RF application for pain relief. The findings of this study could provide first-hand quantitative information to clinical practitioners about the impact of such heterogeneities on the efficacy of RF procedures, thereby assisting them in developing standardized optimal protocols for different cases of interest.

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.

Low-Level Radiofrequency Exposure Does Not Induce Changes in MSC Biology: An in vitro Study for the Prevention of NIR-Related Damage

Background

The ubiquitous diffusion of radiofrequency (RF) radiation across human living environments has attracted the attention of scientists. Though the adverse health effects of RF exposure remain debatable, it has been reported that the interaction of such radiation with biological macromolecular structures can be deleterious for stem cells, inducing impairment of their main functions involving self-renewal and differentiation.

Purpose

The purpose of this study was to determine whether exposure to RF of 169 megahertz (MHz) that is part of very high radiofrequency (VHF) range 30–300 MHz, could cause damage to stem cells by inducing senescence and loss of regenerative and DNA repair capacity.

Methods

The study was conducted on mesenchymal stromal cells (MSCs) containing a subpopulation of stem cells. The MSCs were exposed to RFs of 169 MHz administered via an open meter 2G “Smart Meter” for different durations of time.

Result

We did not observe modifications in MSC biology as a result of the RF exposure conducted in our experiments.

Conclusion

We concluded that MSCs are insensitive to RF radiation exposure at 169 MHz for various time intervals, including longer durations.

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.

A Novel Method to Increase Tumor Ablation Zones With RFA by Injecting the Cationic Polymer Solution to Tissues: In Vivo and Computational Studies

OBJECTIVE

This study aims to examine, for the first time, the introduction of cationic polymer solutions to improve radiofrequency ablation (RFA) in terms of a potentially enlarged ablation zone.

METHODS

By using in vivo and computational RFA studies, two cationic polymers, Chitooligosaccharides (COS) and carboxymethyl chitosan (CMC), diluted in deionized water, were injected into tissues separately surrounding the RF bipolar electrode prior to power application. A total of 9 rabbits were used to 1) measure the increase in electrical conductivity of tissues injected with the cationic polymer solutions, and 2) explore the enhancement of the ablation performance in RFA trials. A computer model of RFA comprising a model of the solution diffusion with an RF thermal ablation model was also built, validated by the in vivo experiment, to quantitatively study the effect of cationic polymer solutions on ablation performances.

RESULTS

Compared to the control group, the electrical conductivity of rabbit liver tissues was increased by 42.20% (0.282 ± 0.006 vs. 0.401 ± 0.048 S/m, P = 0.001) and 43.97% (0.282 ± 0.006 vs. 0.406 ± 0.042 S/m, P = 0.001) by injecting the COS and CMC solution at the concentration of 100 mg/mL into the tissues, denoted COS_DW100 and CMC_DW100, respectively. Consequently, the in vivo experiments show that the ablation zone was enlarged by 95% (47.6 ± 6.3 vs. 92.6 ± 11.5 mm², P < 0.001) and 87% (47.6± 6.3 vs. 88.8 ± 9.6 mm², P < 0.001) by COS_DW100 and CMC_DW100, respectively. The computer simulation shows that the ablation zone was enlarged by 71% (51.9 vs. 88.7 mm²) and 63% (51.9 vs. 84.7 mm²) by COS_DW100 and CMC_DW100, respectively.

CONCLUSION

The injection of the cationic solution can greatly improve the performance of RFA treatment in terms of enlarging the ablation zone, which is due to the increase in the electrical conductivity of liver tissues surrounding the RF electrode.

SIGNIFICANCE

This study contributes to the improvement of RFA in the treatment of large tumors.

A computational model to investigate the influence of electrode lengths on the single probe bipolar radiofrequency ablation of the liver

BACKGROUND AND OBJECTIVES

Recently, there have been calls for RFA to be implemented in the bipolar mode for cancer treatment due to the benefits it offers over the monopolar mode. These include the ability to prevent skin burns at the grounding pad and to avoid tumour track seeding. The usage of bipolar RFA in clinical practice remains uncommon however, as not many research studies have been carried out on bipolar RFA. As such, there is still uncertainty in understanding the effects of the different RF probe configurations on the treatment outcome of RFA. This paper demonstrates that the electrode lengths have a strong influence on the mechanics of bipolar RFA. The information obtained here may lead to further optimization of the system for subsequent uses in the hospitals.

METHODS

A 2D model in the axisymmetric coordinates was developed to simulate the electro-thermophysiological responses of the tissue during a single probe bipolar RFA. Two different probe configurations were considered, namely the configuration where the active electrode is longer than the ground and the configuration where the ground electrode is longer than the active. The mathematical model was first verified with an existing experimental study found in the literature.

RESULTS

Results from the simulations showed that heating is confined only to the region around the shorter electrode, regardless of whether the shorter electrode is the active or the ground. Consequently, thermal coagulation also occurs in the region surrounding the shorter electrode. This opened up the possibility for a better customized treatment through the development of RF probes with adjustable electrode lengths.

CONCLUSIONS

The electrode length was found to play a significant role on the outcome of single probe bipolar RFA. In particular, the length of the shorter electrode becomes the limiting factor that influences the mechanics of single probe bipolar RFA. Results from this study can be used to further develop and optimize bipolar RFA as an effective and reliable cancer treatment technique.

Design of a Novel Electrode of Radiofrequency Ablation for Large Tumors: In Vitro Validation and Evaluation

In the present study, a monopolar expandable electrode (MEE) in radiofrequency ablation (RFA) proposed in our previous study was validated and evaluated using the in vitro experiment and computer simulation. Two commercial RF electrodes (conventional electrode, CE and umbrella electrode, UE) was used to compare the ablation results with MEE using the in vitro egg white model (experiment and computer simulation) and in vivo liver tumor model (computer simulation) to verify the efficacy of MEE in the large tumor ablation. The sharp increase in impedance during RFA procedures was taken as the termination of RFA protocols. The volume and sphericity of ablation zone generated by MEE, CE, and UE in the in vitro egg white experiment were 75.3 1.6 cm3, 2.7 0.4 cm3, 12.4 1.8 cm3 (P <0.001), and 88.1 0.9%, 12.9 1.3%, 62.0 3.0% (P <0.001), respectively. Correspondingly, a similar result was obtained in the egg white simulation. In the liver tumor simulation, the volume and sphpericity of ablation zone generated by MEE, CE, and UE were 35.4 cm3 and 86.8%, 3.7 cm3 and 17.7%, and 12.7 cm3 and 59.6%, respectively. In summary, MEE has the potential to achieve complete ablation in the treatment of large tumors (>3 cm in diameter) compared with CE and UE due to the larger electrode-tissue interface and more round shape of hooks.

Quantification of Thermal Injury to the Healthy Tissue Due to Imperfect Electrode Placements During Radiofrequency Ablation of Breast Tumor

Radiofrequency ablation (RFA) has emerged as an alternative treatment modality for treating various tumors with minimum intervention. The application of RFA in treating breast tumor is still in its infancy stage. Nevertheless, promising results have been obtained while treating early stage localized breast cancer with RFA procedure. The outcome of RFA is tremendously dependent on the precise insertion of the electrode into the geometric center of the tumor. However, there remains plausible chances of inaccuracies in the electrode placement that can result in slight displacement of the electrode tip from the actual desired location during temperature-controlled RFA application. The present numerical study aims at capturing the influence of inaccuracies in electrode placement on the input energy, treatment time and damage to the surrounding healthy tissue during RFA of breast tumor. A thermo-electric analysis has been performed on three-dimensional heterogeneous model of multilayer breast with an embedded early stage spherical tumor of 1.5 cm. The temperature distribution during the RFA has been obtained by solving the coupled electric field equation and Pennes bioheat transfer equation, while the ablation volume has been computed using the Arrhenius cell death model. It has been found that significant variation in the energy consumption, time required for complete tumor necrosis, and the shape of ablation volume among different positions of the electrode considered in this study are prevalent.

Design of a Novel Electrode of Radiofrequency Ablation for Large Tumors: A Finite Element Study

The aim of the study was to design a novel radiofrequency (RF) electrode for larger and rounder ablation volumes and its ability to achieve the complete ablation of liver tumors larger than 3 cm in diameter using finite element method. A new RF expandable electrode comprising three parts (i.e., insulated shaft, changing shaft, and hooks) was designed. Two modes of this new electrode, such as monopolar expandable electrode (MEE) and hybrid expandable electrode (HEE), and a commercial expandable electrode (CEE) were investigated using liver tissue with (scenario I) and without (scenario II) a liver tumor. A temperature-controlled radiofrequency ablation (RFA) protocol with a target temperature of 95 °C and an ablation time of 15 min was used in the study. Both the volume and shape of the ablation zone were examined for all RF electrodes in scenario I. Then, the RF electrode with the best performance in scenario I and CEE were used to ablate a large liver tumor with the diameter of 3.5 cm (scenario II) to evaluate the effectiveness of complete tumor ablation of the designed RF electrode. In scenario I, the ablation volumes of CEE, HEE, and MEE were 12.11 cm3, 33.29 cm3, and 48.75 cm3, respectively. The values of sphericity index (SI) of CEE, HEE, and MEE were 0.457, 0.957, and 0.976, respectively. The best performance was achieved by using MEE. In scenario II, the ablation volumes of MEE and CEE were 71.59 cm3 and 19.53 cm3, respectively. Also, a rounder ablation volume was achieved by using MEE compared to CEE (SI: 0.978 versus 0.596). The study concluded that: (1) compared with CEE, both MEE and HEE get larger and rounder ablation volumes due to the larger electrode–tissue interface and rounder shape of hook deployment; (2) MEE has the best performance in getting a larger and rounder ablation volume; and (3) computer simulation result shows that MEE is also able to ablate a large liver tumor (i.e., 3.5 cm in diameter) completely, which has at least 0.785 cm safety margin.

Comparison of intraluminal radiofrequency ablation and stents vs. stents alone in the management of malignant biliary obstruction

PURPOSE

To retrospectively evaluate the added benefit of adding intraluminal radiofrequency ablation (RFA) to biliary metal stent placement for patients with malignant biliary obstruction (MBO).

METHODS

From November 2013 to December 2015, 89 patients with MBO who had undergone percutaneous intraluminal RFA and stent placement (RFA-stent group, n = 50) or stent placement only (stent group, n = 39) were included. Outcomes were compared according to the type of tumour: cholangiocarcinoma or non-cholangiocarcinoma.

RESULTS

Primary and secondary stent patency (PSP, SSP) were significantly higher for the RFA-stent group than the stent group (PSP: 7.0 months vs. 5.0 months, p = 0.006; SSP: 10.0 months vs. 5.6 months, p < 0.001), with overall survival being comparable (5.0 months vs. 4.7 months, p = 0.068). In subgroup analysis, RFA-stent showed significant PSP benefits compared to stent alone in patients with cholangiocarcinoma (7.4 months vs. 4.3 months; p = 0.009), but with comparable outcomes in patients with non-cholangiocarcinoma (6.3 months vs. 5.2 months; p = 0.266). The SSP was improved in both subgroups (cholangiocarcinoma, 12.6 months vs. 5.0 months, p < 0.001; non-cholangiocarcinoma, 10.3 months vs. 5.5 months, p = 0.013). Technical success and clinical success were not significantly different between the two groups. The rate of complication was higher for the RFA-stent group, but tolerable when compared to the stent group.

CONCLUSIONS

Although survival was comparable between the groups, RFA-stent confers therapeutic benefits to patients with MBO in terms of stent patency compared to stent placement alone, especially in those with cholangiocarcinoma.