Finite element analysis of hepatic radiofrequency ablation probes using temperature-dependent electrical conductivity

Safety and Effectiveness of Triple-Antenna Hepatic Microwave Ablation

Microwave ablation is becoming a standard procedure for treating tumors based on heat generation, causing an elevation in the tissue temperature level from 50 to 60 °C, causing tissue death. Microwave ablation is associated with uniform cell killing within ablation zones, multiple-antenna capability, low complication rates, and long-term survival. Several reports have demonstrated that multiple-antenna microwave ablation is a promising strategy for safely, rapidly, and effectively treating large tumors. The key advantage of multi-antenna tumor microwave ablation is the creation of a large, well-defined ablation zone without excessively long treatment times or high power that can damage healthy tissue. The strategic positioning of multiple probes provides a fully ablated volume, even in regions where individual probe damage is incomplete. Accurate modeling of the complex thermal and electromagnetic behaviors of tissue is critical for optimizing microwave ablation because material parameters and tissue responses can change significantly during the procedure. In the case of multi-antenna microwave ablation, the calculation complexity increases significantly, requiring significant computational resources and time. This study aimed to evaluate the efficacy and safety of liver percutaneous microwave ablation using the simultaneous activation of three antennas for the treatment of lesions larger than 3 cm. Based on the known results from a single-probe setup, researchers can estimate and evaluate various spatial configurations of the three-probe array to identify the optimal arrangement. Due to the synergistic effects of the combined radiation from the three antennas, the resulting ablation zone can be significantly larger, leading to better outcomes in terms of treatment time and effectiveness. The obtained results revealed that volumetric damage and the amount of damaged healthy tissue are smaller for a three-antenna configuration than for microwave ablation using a single-antenna and two-antenna configurations.

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.

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.

Bioheat Model of Spinal Column Heating During High-Density Spinal Cord Stimulation

INTRODUCTION

High-density (HD) spinal cord stimulation (SCS) delivers higher charge per time by increasing frequency and/or pulse duration, thus increasing stimulation energy. Previously, through phantom studies and computational modeling, we demonstrated that stimulation energy drives spinal tissue heating during kHz SCS. In this study, we predicted temperature increases in the spinal cord by HD SCS, the first step in considering the potential impact of heating on clinical outcomes.

MATERIALS AND METHODS

We adapted a high-resolution computer-aided design-derived spinal cord model, both with and without a lead encapsulation layer, and applied bioheat transfer finite element method multiphysics to predict temperature increases during SCS. We simulated HD SCS using a commercial SCS lead (eight contacts) with clinically relevant intensities (voltage-controlled: 0.5-7 Vrms) and electrode configuration (proximal bipolar, distal bipolar, guarded tripolar [+-+], and guarded quadripolar [+--+]). Results were compared with the conventional and 10-kHz SCS (current-controlled).

RESULTS

HD SCS waveform energy (reflecting charge per second) governs joule heating in the spinal tissues, increasing temperature supralinearly with stimulation root mean square. Electrode configuration and tissue properties (an encapsulation layer) influence peak tissue temperature increase-but in a manner distinct for voltage-controlled (HD SCS) compared with current-controlled (conventional/10-kHz SCS) stimulation. Therefore, depending on conditions, HD SCS could produce heating greater than that of 10-kHz SCS. For example, with an encapsulation layer, using guarded tripolar configuration (500-Hz, 250-μs pulse width, 5-Vpeak HD SCS), the peak temperature increases were 0.36 °C at the spinal cord and 1.78 °C in the epidural space.

CONCLUSIONS

As a direct consequence of the higher charge, HD SCS increases tissue heating; voltage-controlled stimulation introduces special dependencies on electrode configuration and lead encapsulation (reflected in impedance). If validated with an in vivo measurement as a possible mechanism of action of SCS, bioheat models of HD SCS serve as tools for programming optimization.

Tissue Temperature Increases by a 10 kHz Spinal Cord Stimulation System: Phantom and Bioheat Model

OBJECTIVE

A recently introduced Spinal Cord Stimulation (SCS) system operates at 10 kHz, faster than conventional SCS systems, resulting in significantly more power delivered to tissues. Using a SCS heat phantom and bioheat multi-physics model, we characterized tissue temperature increases by this 10 kHz system. We also evaluated its Implanted Pulse Generator (IPG) output compliance and the role of impedance in temperature increases.

MATERIALS AND METHODS

The 10 kHz SCS system output was characterized under resistive loads (1-10 KΩ). Separately, fiber optic temperature probes quantified temperature increases (ΔTs) around the SCS lead in specially developed heat phantoms. The role of stimulation Level (1-7; ideal pulse peak-to-peak of 1-7mA) was considered, specifically in the context of stimulation current Root Mean Square (RMS). Data from the heat phantom were verified with the SCS heat-transfer models. A custom high-bandwidth stimulator provided 10 kHz pulses and sinusoidal stimulation for control experiments.

RESULTS

The 10 kHz SCS system delivers 10 kHz biphasic pulses (30-20-30 μs). Voltage compliance was 15.6V. Even below voltage compliance, IPG bandwidth attenuated pulse waveform, limiting applied RMS. Temperature increased supralinearly with stimulation Level in a manner predicted by applied RMS. ΔT increases with Level and impedance until stimulator compliance was reached. Therefore, IPG bandwidth and compliance dampen peak heating. Nonetheless, temperature increases predicted by bioheat multi-physic models (ΔT = 0.64°C and 1.42°C respectively at Level 4 and 7 at the cervical segment; ΔT = 0.68°C and 1.72°C respectively at Level 4 and 7 at the thoracic spinal cord)-within ranges previously reported to effect neurophysiology.

CONCLUSIONS

Heating of spinal tissues by this 10 kHz SCS system theoretically increases quickly with stimulation level and load impedance, while dampened by IPG pulse bandwidth and voltage compliance limitations. If validated in vivo as a mechanism of kHz SCS, bioheat models informed by IPG limitations allow prediction and optimization of temperature changes.

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.

Kilohertz-frequency stimulation of the nervous system: A review of underlying mechanisms

BACKGROUND

Electrical stimulation in the kilohertz-frequency range has gained interest in the field of neuroscience. The mechanisms underlying stimulation in this frequency range, however, are poorly characterized to date.

OBJECTIVE/HYPOTHESIS

To summarize the manifold biological effects elicited by kilohertz-frequency stimulation in the context of the currently existing literature and provide a mechanistic framework for the neural responses observed in this frequency range.

METHODS

A comprehensive search of the peer-reviewed literature was conducted across electronic databases. Relevant computational, clinical, and mechanistic studies were selected for review.

RESULTS

The effects of kilohertz-frequency stimulation on neural tissue are diverse and yield effects that are distinct from conventional stimulation. Broadly, these can be divided into 1) subthreshold, 2) suprathreshold, 3) synaptic and 4) thermal effects. While facilitation is the dominating mechanism at the subthreshold level, desynchronization, spike-rate adaptation, conduction block, and non-monotonic activation can be observed during suprathreshold kilohertz-frequency stimulation. At the synaptic level, kilohertz-frequency stimulation has been associated with the transient depletion of the available neurotransmitter pool - also known as synaptic fatigue. Finally, thermal effects associated with extrinsic (environmental) and intrinsic (associated with kilohertz-frequency stimulation) temperature changes have been suggested to alter the neural response to stimulation paradigms.

CONCLUSION

The diverse spectrum of neural responses to stimulation in the kilohertz-frequency range is distinct from that associated with conventional stimulation. This offers the potential for new therapeutic avenues across stimulation modalities. However, stimulation in the kilohertz-frequency range is associated with distinct challenges and caveats that need to be considered in experimental paradigms.

Radiofrequency ablation with four electrodes as a building block for matrix radiofrequency ablation: Ex vivo liver experiments and finite element method modelling. Influence of electric and activation mode on coagulation size and geometry

PURPOSE

Radiofrequency ablation (RFA) is increasingly being used to treat unresectable liver tumors. Complete ablation of the tumor and a safety margin is necessary to prevent local recurrence. With current electrodes, size and shape of the ablation zone are highly variable leading to unsatisfactory local recurrence rates, especially for tumors >3 cm. In order to improve predictability, we recently developed a system with four simple electrodes with complete ablation in between the electrodes. This rather small but reliable ablation zone is considered as a building block for matrix radiofrequency ablation (MRFA). In the current study we explored the influence of the electric mode (monopolar or bipolar) and the activation mode (consecutive, simultaneous or switching) on the size and geometry of the ablation zone.

MATERIALS AND METHODS

The four electrode system was applied in ex vivo bovine liver. The electric and the activation mode were changed one by one, using constant power of 50 W in all experiments. Size and geometry of the ablation zone were measured. Finite element method (FEM) modelling of the experiment was performed.

RESULTS

In ex vivo liver, a complete and predictable coagulation zone of a 3 × 2 × 2 cm block was obtained most efficiently in the bipolar simultaneous mode due to the combination of the higher heating efficacy of the bipolar mode and the lower impedance by the simultaneous activation of four electrodes, as supported by the FEM simulation.

CONCLUSIONS

In ex vivo liver, the four electrode system used in a bipolar simultaneous mode offers the best perspectives as building block for MRFA. These results should be confirmed by in vivo experiments.

Determination of Tissue Thermal Conductivity as a Function of Thermal Dose and Its Application in Finite Element Modeling of Electrosurgical Vessel Sealing

OBJECTIVE

Electrosurgical vessel sealing is a process commonly used to control bleeding during surgical procedures. Finite element (FE) modeling is often performed to obtain a better understanding of thermal spread during this process. The accuracy of the FE model depends on the implemented material properties. Thermal conductivity is one of the most important properties that affect temperature distribution. The goal of this study is to determine the tissue thermal conductivity as a function of thermal dose. Methods: We developed an iterative approach to correlating tissue thermal conductivity to more accurately calculated thermal dose, which cannot be experimentally measured. The resulting regression model was then implemented into an electrosurgical vessel sealing FE model to examine the accuracy of this FE model. Results: The results show that with the regression model, more reasonable temperature and thermal dose prediction can be achieved at the center of the sealed vessel tissue. The resulting electrical current and impedance from the FE model match with the experimental results. Conclusion: The developed approach can be used to determine the correlation between thermal dose and thermal conductivity. Describing the thermal conductivity as a function of thermal dose allows modeling of irreversible changes in tissue properties. Significance: By having a more accurate temperature estimation at the center of the sealed vessel, more insight is provided into how the tissue reacts during the vessel sealing process.

Computer modeling and in vitro experimental study of water-cooled microwave ablation array

INTRODUCTION

Microwaves (MWs) quickly deliver relatively high temperatures into tumors and cover a large ablation zone. We present a research protocol for using water-cooled double-needle MW ablation arrays for tumor ablation here.

MATERIAL AND METHODS

Our research program includes computer modeling, tissue-mimicking phantom experiments, and in vitro swine liver experiments. The computer modeling is based on the finite element method (FEM) to evaluate ablation temperature distributions. In tissue-mimicking phantom and in vitro swine liver ablation experiments, the performances of the new device and the single-needle MW device currently used in clinical practice are compared.

RESULTS

FEM shows that the maximum transverse ablation diameter (MTAD) is 4.2 cm at 100 W output and 300 s (assessed at the 50 °C isotherm). In the tissue-mimicking phantom, the MTDA is 2.6 cm at 50 W and 300 s in single-needle MW ablation, and 4 cm in double needle MW ablation array. In in vitro swine liver experiments, the MTAD is 2.820 ± 0.127 cm at 100 W and 300 s in single-needle MW ablation, and 3.847 ± 0.103 cm in MW ablation array.

CONCLUSION

A new type of water-cooled MW ablation array is designed and tested, and has potential advantages over currently used devices.

Real-time monitoring radiofrequency ablation using tree-based ensemble learning models

BACKGROUND

Radiofrequency ablation is a minimally-invasive treatment method that aims to destroy undesired tissue by exposing it to alternating current in the 100 kHz-800 kHz frequency range and heating it until it is destroyed via coagulative necrosis. Ablation treatment is gaining momentum especially in cancer research, where the undesired tissue is a malignant tumor. While ablating the tumor with an electrode or catheter is an easy task, real-time monitoring the ablation process is a must in order to maintain the reliability of the treatment. Common methods for this monitoring task have proven to be accurate, however, they are all time-consuming or require expensive equipment, which makes the clinical ablation process more cumbersome and expensive due to the time-dependent nature of the clinical procedure.

METHODS

A machine learning (ML) approach is presented that aims to reduce the monitoring time while keeping the accuracy of the conventional methods. Two different hardware setups are used to perform the ablation and collect impedance data at the same time and different ML algorithms are tested to predict the ablation depth in 3 dimensions, based on the collected data.

RESULTS

Both the random forest and adaptive boosting (adaboost) models had over 98% R² on the data collected with the embedded system-based hardware instrumentation setup, outperforming Neural Network-based models.

CONCLUSIONS

It is shown that an optimal pair of hardware setup and ML algorithm (Adaboost) is able to control the ablation by estimating the lesion depth within a test average of 0.3mm while keeping the estimation time within 10ms on a ×86-64 workstation.

Waveform Dependent Electrosurgical Effects on Soft Hydrated Tissues

Electrosurgical procedures are ubiquitously used in surgery. The commonly used power modes, including the coagulation and blend modes, utilize non-sinusoidal or modulated current waveforms. For the same power setting, the coagulation, blend and pure cutting modes have different heating and thermal damage outcomes due to the frequency dependence of electrical conductivity of soft hydrated tissues. In this paper, we propose a multi-physics model of soft tissues to account for the effects of multi-frequency electrosurgical power modes within the framework of a continuum thermomechanical model based on mixture theory. Electrical and frequency spectrum results from different power modes at low and high power settings are presented. Model predictions are compared with in vivo electrosurgical heating experiments on porcine liver tissue. The accuracy of the model in predicting experimentally observed temperature profiles is found to be overall greater when frequency-dependence is included. An Arrhenius type model indicates that more tissue damage is correlated with larger duty cycles in multi-frequency modes.

A continuum thermomechanical model of in vivo electrosurgical heating of hydrated soft biological tissues

Radio-frequency (RF) heating of soft biological tissues during electrosurgical procedures is a fast process that involves phase change through evaporation and transport of intra- and extra-cellular water, and where variations in physical properties with temperature and water content play significant role. Accurately predicting and capturing these effects would improve the modeling of temperature change in the tissue allowing the development of improved instrument design and better understanding of tissue damage and necrosis. Previous models based on the Pennes' bioheat model neglect both evaporation and transport or consider evaporation through numerical correlations, however, do not account for changes in physical properties due to mass transport or phase change, nor capture the pressure increase due to evaporation within the tissue. While a porous media approach can capture the effects of evaporation, transport, pressure and changes in physical properties, the model assumes free diffusion of liquid and gas without a careful examination of assumptions on transport parameters in intact tissue resulting in significant under prediction of temperature. These different approaches have therefore been associated with errors in temperature prediction exceeding 20% when compared to experiments due to inaccuracies in capturing the effects of evaporation losses and transport. Here, we present a model of RF heating of hydrated soft tissue based on mixture theory where the multiphase nature of tissue is captured within a continuum thermomechanics framework, simultaneously considering the transport, deformation and phase change losses due to evaporation that occur during electrosurgical heating. The model predictions are validated against data obtained for in vivo ablation of porcine liver tissue at various power settings of the electrosurgical unit. The model is able to match the mean experimental temperature data with sharp gradients in the vicinity of the electrode during rapid low and high power ablation procedures with errors less than 7.9%. Additionally, the model is able to capture fast vaporization losses and the corresponding increase in pressure due to vapor buildup which have a significant effect on temperature prediction beyond 100 °C.

Temperature increases by kilohertz frequency spinal cord stimulation

INTRODUCTION

Kilohertz frequency spinal cord stimulation (kHz-SCS) deposits significantly more power in tissue compared to SCS at conventional frequencies, reflecting increased duty cycle (pulse compression). We hypothesize kHz-SCS increases local tissue temperature by joule heat, which may influence the clinical outcomes.

METHODS

To establish the role of tissue heating in KHZ-SCS, a decisive first step is to characterize the range of temperature changes expected during conventional and KHZ-SCS protocols. Fiber optic probes quantified temperature increases around an experimental SCS lead in a bath phantom. These data were used to verify a SCS lead heat-transfer model based on joule heat. Temperature increases were then predicted in a seven-compartment (soft tissue, vertebral bone, fat, intervertebral disc, meninges, spinal cord with nerve roots) geometric human spinal cord model under varied parameterization.

RESULTS

The experimentally constrained bio-heat model shows SCS waveform power (waveform RMS) determines tissue heating at the spinal cord and surrounding tissues. For example, we predict temperature increased at dorsal spinal cord of 0.18-1.72 °C during 3.5 mA peak 10 KHz stimulation with a 40-10-40 μs biphasic pulse pattern, 0.09-0.22 °C during 3.5 mA 1 KHz 100-100-100 μs stimulation, and less than 0.05 °C during 3.5 mA 50 Hz 200-100-200 μs stimulation. Notably, peak heating of the spinal cord and other tissues increases superlinearly with stimulation power and so are especially sensitive to incremental changes in SCS pulse amplitude or frequency (with associated pulse compression). Further supporting distinct SCS intervention strategies based on heating; the spatial profile of temperature changes is more uniform compared to electric fields, which suggests less sensitivity to lead position.

CONCLUSIONS

Tissue heating may impact short and long-term outcomes of KHZ-SCS, and even as an adjunct mechanism, suggests distinct strategies for lead position and programming optimization.

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.

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.

Energy Dissipation in Ex Vivo Porcine Liver During Electrosurgery

This paper explores energy dissipation in ex vivo liver tissue during radiofrequency current excitation with application in electrosurgery. Tissue surface temperature for monopolar electrode configuration is measured using infrared thermometry. The experimental results are fitted to a finite-element model for transient heat transfer taking into account energy storage and conduction in order to extract information about "apparent" specific heat, which encompasses storage and phase change. The average apparent specific heat determined for low temperatures is in agreement with published data. However, at temperatures approaching the boiling point of water, apparent specific heat increases by a factor of five, indicating that vaporization plays an important role in the energy dissipation through latent heat loss.

Thermal analysis of induced damage to the healthy cell during RFA of breast tumor

Effective pre-clinical computational modeling strategies have been demonstrated in this article to enable risk free clinical application of radiofrequency ablation (RFA) of breast tumor. The present study (a) determines various optimal regulating parameters required for RFA of tumor and (b) introduces an essential clinical monitoring scheme to minimize the extent of damage to the healthy cell during RFA of tumor. The therapeutic capabilities offered by RFA of breast tumor, viz., the rise in local temperature and induced thermal damage have been predicted by integrating the bioheat transfer model, the electric field distribution model and the thermal damage model. The mathematical model has been validated with the experimental results available in the literature. The results revealed that, the effective damage of tumor volume sparing healthy tissue essentially depends on the voltage, the exposure time, the local heat distribution, the tumor stage and the electrode geometric configuration. It has been confirmed that, the assessment of damage front can accurately determine the extent of damage as compared to the thermal front. The study further evaluates the damaged healthy and tumor volumes due to RFA of different stages of breast cancer. The assessment of cell survival and damage fractions discloses the propensity of reappearance/healing of tumor cells after treatment.

Cell death, perfusion and electrical parameters are critical in models of hepatic radiofrequency ablation

Purpose: A sensitivity analysis has been performed on a mathematical model of radiofrequency ablation (RFA) in the liver. The purpose of this is to identify the most important parameters in the model, defined as those that produce the largest changes in the prediction. This is important in understanding the role of uncertainty and when comparing the model predictions to experimental data. Materials and methods: The Morris method was chosen to perform the sensitivity analysis because it is ideal for models with many parameters or that take a significant length of time to obtain solutions. A comprehensive literature review was performed to obtain ranges over which the model parameters are expected to vary, crucial input information. Results: The most important parameters in predicting the ablation zone size in our model of RFA are those representing the blood perfusion, electrical conductivity and the cell death model. The size of the 50 °C isotherm is sensitive to the electrical properties of tissue while the heat source is active, and to the thermal parameters during cooling. Conclusions: The parameter ranges chosen for the sensitivity analysis are believed to represent all that is currently known about their values in combination. The Morris method is able to compute global parameter sensitivities taking into account the interaction of all parameters, something that has not been done before. Research is needed to better understand the uncertainties in the cell death, electrical conductivity and perfusion models, but the other parameters are only of second order, providing a significant simplification.

Fourier, hyperbolic and relativistic heat transfer equations: a comparative analytical study

Parabolic heat equation based on Fourier's theory (FHE), and hyperbolic heat equation (HHE), has been used to mathematically model the temperature distributions of biological tissue during thermal ablation. However, both equations have certain theoretical limitations. The FHE assumes an infinite thermal energy propagation speed, whereas the HHE might possibly be in breach of the second law of thermodynamics. The relativistic heat equation (RHE) is a hyperbolic-like equation, whose theoretical model is based on the theory of relativity and which was designed to overcome these theoretical impediments. In this study, the three heat equations for modelling of thermal ablation of biological tissues (FHE, HHE and RHE) were solved analytically and the temperature distributions compared. We found that RHE temperature values were always lower than those of the FHE, while the HHE values were higher than the FHE, except for the early stages of heating and at points away from the electrode. Although both HHE and RHE are mathematically hyperbolic, peaks were only found in the HHE temperature profiles. The three solutions converged for infinite time or infinite distance from the electrode. The percentage differences between the FHE and the other equations were larger for higher values of thermal relaxation time in HHE.

The effect of the shape and size of gold seeds irradiated with ultrasound on the bio-heat transfer in tissue

The aim of this report is to propose a new methodology to treat prostate cancer with macro-rod-shaped gold seeds irradiated with ultrasound and develop a new computational method for temperature and thermal dose control of hyperthermia therapy induced by the proposed procedure. A computer code representation, based on the bio-heat diffusion equation, was developed to calculate the heat deposition and temperature elevation patterns in a gold rod and in the tissue surrounding it as a result of different therapy durations and ultrasound power simulations. The numerical results computed provide quantitative information on the interaction between high-energy ultrasound, gold seeds and biological tissues and can replicate the pattern observed in experimental studies. The effect of differences in shapes and sizes of gold rod targets irradiated with ultrasound is calculated and the heat enhancement and the bio-heat transfer in tissue are analyzed.

Flu-like symptoms following radiofrequency liver transection: a new variety of the post-radiofrequency syndrome

BACKGROUND/AIMS

The aim of our study was to determine whether post-radiofrequency syndrome may also develop following hepatectomy using saline-cooled radiofrequency coagulation.

METHODS

We retrospectively reviewed 95 consecutive patients who underwent 110 liver resections between May 2000 and September 2012. We stated that 80.9% of the resections were carried out employing the saline-cooled radiofrequency device. All medical records were searched for the occurrence of flu-like symptoms, without evidence of sepsis or infection, in the first two postoperative weeks.

RESULTS

Eleven patients (11.5%) developed flu-like symptoms after hepatectomy without evidence of sepsis or infection. All their hepatectomies were performed employing the saline-cooled radiofrequency probe (p = .089), and all cases but one appeared following colorectal liver metastases surgery (p = .042). Eight of them were readmitted to the hospital because of their symptoms. In all 11 cases, a fluid collection was present, 8 of them with gas presence. Nine patients underwent a percutaneous drainage whose cultures were negative. Ten patients recovered without treatment or with the intake of nonsteroidal anti-inflammatory drugs within 1 week, but one patient developed a secondary infection with gram-positive bacteria after percutaneous drainages that prolonged his hospital stay.

CONCLUSION

Liver splitting using saline-cooled radiofrequency coagulation may cause postoperative symptoms that may mimic surgical site infection. Surgeons employing this device should keep this in mind to avoid potentially unwarranted treatments that may be unnecessary, expensive, and even harmful.

Histological and finite element analysis of cell death due to irreversible electroporation

Irreversible electroporation (IRE) has been shown to be an effective method of killing cells locally. In contrast to radiofrequency ablation, the mechanism by which cells are thought to die via IRE is the creation of pores in cell membranes, without substantial increase in tissue temperature. To determine the degree to which cell death is non-thermal, we evaluated IRE in porcine hepatocytes in vivo. Using pulse widths of 10μs, bursts of 3 kV square-wave pulses were applied through a custom probe to the liver of an anesthetized pig. Affected tissue was evaluated histologically via stainings of hematoxylin & eosin (H&E), nitroblue tetrazolium (NBT) to monitor cell respiration and TUNEL to gauge apoptosis. Temperature was measured during the application of electroporation, and heat transfer was modeled via finite element analysis. Cell death was calculated via Arrhenius kinetics. Four distinct zones were observed within the ring return electrode; heat-fixed tissue, coagulation, necrotic, and viable. The Arrhenius damage integral estimated complete cell death only in the first zone, where the temperature exceeded 70°C, and partial or no cell death in the other zones, where maximum temperature was approximately 45°C. Except for a limited area near the electrode tip, cell death in IRE is predominantly due to a non-thermal mechanism.

Validation of accuracy of liver model with temperature-dependent thermal conductivity by comparing the simulation and in vitro RF ablation experiment

Radiofrequency (RF) ablation is increasingly used to treat cancer because it is minimally invasive. However, it is difficult for operators to control precisely the formation of coagulation zones because of the inadequacies of imaging modalities. To overcome this limitation, we previously proposed a model-based robotic ablation system that can create the required size and shape of coagulation zone based on the dimensions of the tumor. At the heart of such a robotic system is a precise temperature distribution simulator for RF ablation. In this article, we evaluated the simulation accuracy of two numerical simulation liver models, one using a constant thermal conductivity value and the other using temperature-dependent thermal conductivity values, compared with temperatures obtained using in vitro experiments. The liver model that reflected the temperature dependence of thermal conductivity did not result in a large increase of simulation accuracy compared with the temperature-independent model in the temperature range achieved during clinical RF ablation.

Temperature control at DBS electrodes using a heat sink: experimentally validated FEM model of DBS lead architecture

There is a growing interest in the use of deep brain stimulation (DBS) for the treatment of medically refractory movement disorders and other neurological and psychiatric conditions. The extent of temperature increases around DBS electrodes during normal operation (joule heating and increased metabolic activity) or coupling with an external source (e.g. magnetic resonance imaging) remains poorly understood and methods to mitigate temperature increases are being actively investigated. We developed a heat transfer finite element method (FEM) simulation of DBS incorporating the realistic architecture of Medtronic 3389 leads. The temperature changes were analyzed considering different electrode configurations, stimulation protocols and tissue properties. The heat-transfer model results were then validated using micro-thermocouple measurements during DBS lead stimulation in a saline bath. FEM results indicate that lead design (materials and geometry) may have a central role in controlling temperature rise by conducting heat. We show how modifying lead design can effectively control temperature increases. The robustness of this heat-sink approach over complimentary heat-mitigation technologies follows from several features: (1) it is insensitive to the mechanisms of heating (e.g. nature of magnetic coupling); (2) it does not interfere with device efficacy; and (3) can be practically implemented in a broad range of implanted devices without modifying the normal device operations or the implant procedure.

Tissue temperature feedback control of power: the key to successful ablation

Multiple ablation technologies are used to treat atrial fibrillation during cardiac operations. All such ablation technologies use locally induced temperature extremes (>50°C or <-20°C) to kill tissue and create a lesion pattern in the atria which blocks activation pathways that initiate and sustain atrial fibrillation. The technologies used to heat tissue have included radiofrequency (RF), microwave, high-intensity focused ultrasound, and infrared laser. RF accounts for more than 95% of the heating-based ablation technology used by cardiac surgeons. Energy delivery with RF is easier to control than with some other technologies, the heating produced by the energy source is well understood, and manufacturing costs are not excessive. Whichever heating technology is used, control of energy delivery is required to ensure both safe and effective heating of the targeted tissue. All targeted tissue needs to be heated above 50°C to achieve cell death. However, the targeted tissue should not be heated above 100°C, as this can cause perforation due to a steam pop. In addition, adjacent noncardiac tissues must not be damaged during the ablation procedure. The best method to achieve this control uses direct measurement of tissue temperature, because the tissue temperature defines both the safe and effective limits for the ablative process.

Treatment of advanced hepatocellular carcinoma with very low levels of amplitude-modulated electromagnetic fields

Background:

Therapeutic options for patients with advanced hepatocellular carcinoma (HCC) are limited. There is emerging evidence that the growth of cancer cells may be altered by very low levels of electromagnetic fields modulated at specific frequencies.

Methods:

A single-group, open-label, phase I/II study was performed to assess the safety and effectiveness of the intrabuccal administration of very low levels of electromagnetic fields amplitude modulated at HCC-specific frequencies in 41 patients with advanced HCC and limited therapeutic options. Three-daily 60-min outpatient treatments were administered until disease progression or death. Imaging studies were performed every 8 weeks. The primary efficacy end point was progression-free survival ⩾6 months. Secondary efficacy end points were progression-free survival and overall survival.

Results:

Treatment was well tolerated and there were no NCI grade 2, 3 or 4 toxicities. In all, 14 patients (34.1%) had stable disease for more than 6 months. Median progression-free survival was 4.4 months (95% CI 2.1–5.3) and median overall survival was 6.7 months (95% CI 3.0–10.2). There were three partial and one near complete responses.

Conclusion:

Treatment with intrabuccally administered amplitude-modulated electromagnetic fields is safe, well tolerated, and shows evidence of antitumour effects in patients with advanced HCC.

Temperature dependence of thermal conductivity of liver based on various experiments and a numerical simulation for RF ablation

Radiofrequency ablation (RFA) for liver cancer has increasingly been used over the past few years because RFA is minimally invasive treatment for patients. However, precise control of the formation of coagulation zones is difficult for operators due to inadequate imaging modalities. With this in mind, we have proposed a model-based robotic ablation system using numerical simulation to analyze temperature distributions in the organ to overcome this deficiency. The objective of our work is to develop a temperature-dependent thermophysical organ model to construct a precise numerical simulator for RFA. However, no standard methods exist for obtaining the thermophysical properties of biological tissues, as detailed evaluations of the accuracy of properties obtained from various experiments have not been completed. The purpose of this study was thus to measure and model the temperature dependence of thermal conductivity in hog liver from three representative methods, and to compare these results using our developed numerical simulator to reveal differences in temperature distributions stemming from differences in thermal conductivities.

Targeted tissue ablation with nanosecond pulses

In-vivo porcine studies on the effect of nanosecond high voltage pulses on liver tissue have shown that cell death can be induced in well-defined tissue volumes without damaging collagen-predominant structures. Comparison of the experimental results with the results of a three-dimensional finite element model allowed us to determine the threshold electric field for cell death. For 30, 100 nanosecond long pulses this was found to be in the range from 12 to 15 kV/cm. Modelling of the temperature distribution in the tissue using Pennes' bioheat equation showed that the lethal effect of nanosecond pulses on cells is non-thermal. Muscle contractions, generally caused by high voltage pulses, were significantly reduced for the 100 nanosecond pulses compared to microsecond long pulses. The results of these studies indicate that high voltage nanosecond pulses reliably kill normal liver cells in vivo and therefore may be useful for liver tumor treatments.

A theoretical model of the application of RF energy to the airway wall and its experimental validation

Background

Bronchial thermoplasty is a novel technique designed to reduce an airway's ability to contract by reducing the amount of airway smooth muscle through controlled heating of the airway wall. This method has been examined in animal models and as a treatment for asthma in human subjects. At the present time, there has been little research published about how radiofrequency (RF) energy and heat is transferred to the airways of the lung during bronchial thermoplasty procedures. In this manuscript we describe a computational, theoretical model of the delivery of RF energy to the airway wall.

Methods

An electro-thermal finite-element-analysis model was designed to simulate the delivery of temperature controlled RF energy to airway walls of the in vivo lung. The model includes predictions of heat generation due to RF joule heating and transfer of heat within an airway wall due to thermal conduction. To implement the model, we use known physical characteristics and dimensions of the airway and lung tissues. The model predictions were tested with measurements of temperature, impedance, energy, and power in an experimental canine model.

Results

Model predictions of electrode temperature, voltage, and current, along with tissue impedance and delivered energy were compared to experiment measurements and were within ± 5% of experimental averages taken over 157 sample activations.

The experimental results show remarkable agreement with the model predictions, and thus validate the use of this model to predict the heat generation and transfer within the airway wall following bronchial thermoplasty.

Conclusions

The model also demonstrated the importance of evaporation as a loss term that affected both electrical measurements and heat distribution. The model predictions showed excellent agreement with the empirical results, and thus support using the model to develop the next generation of devices for bronchial thermoplasty. Our results suggest that comparing model results to RF generator electrical measurements may be a useful tool in the early evaluation of a model.

A randomized comparison of a new flexible bipolar hemostasis forceps designed principally for NOTES versus a conventional surgical laparoscopic bipolar forceps for intra-abdominal vessel sealing in a porcine model

BACKGROUND

Current devices for hemostasis in flexible endoscopy are inferior to methods used during open or laparoscopic surgery and might be ineffective for natural orifice transluminal endoscopic surgery.

OBJECTIVE

To compare new flexible bipolar forceps (FBF), designed principally for natural orifice transluminal endoscopic surgery, with laparoscopic bipolar forceps (LBF) for hemostasis of intra-abdominal porcine arteries.

SETTING

Surgical laboratories in Europe and the United States.

DESIGN AND INTERVENTIONS

New FBF for hemostasis (3.7-mm diameter), featuring electrode isolation, were compared with rigid 5-mm LBF (ERBE BiClamp LAP forceps) at recommended settings. A porcine model of acute hemostasis was prepared by suturing the uterine horns and cecum to the abdominal wall, exposing uterine arteries, ovarian pedicles, cecal mesenteric bundles, and the inferior mesenteric artery. This allowed access to 10 vessels in each pig by transabdominal laparoscopic devices or a transgastric double-channel gastroscope. Vessels were measured, coagulated at 4 and more points, and transected. Blood pressure was increased to more than 200 mm Hg for 10 minutes by administering phenylephrine. Delayed bleeding was identified.

MAIN OUTCOME MEASUREMENTS

In 7 pigs, a total of 65 vessels (1.5-6.0 mm) were randomly allocated to FBF (n = 32) or LBF (n = 33). Successful hemostasis both before and after blood pressure increase was equivalent between the 2 groups (before: 88% FBF vs 88% LBF, not significant [NS]; after: 97% FBF vs 94% LBF, NS). With FBF, the number of seals per vessel was 4.8 vs 4.4 with LBF (NS). The energy used to create FBF seals was 19.8 J vs 38.2 J for LBF (P < .05).

LIMITATIONS

Results from porcine studies may not reflect patient outcomes.

CONCLUSIONS

In a porcine model, transgastric FBF endoscopic hemostasis was as effective as conventional laparoscopic hemostasis using LBF across a wide range of vessels.

Conductive interstitial thermal therapy (CITT) inhibits recurrence and metastasis in rabbit VX2 carcinoma model

PURPOSE

To investigate the potential of conductive interstitial thermal therapy (CITT) to inhibit recurrence and metastasis in a partially resected tumour model.

METHOD

Fifteen New Zealand white rabbits were implanted with VX2 tumour intramuscularly in the rear thigh. Once the tumour size reached 20-25 mm in diameter, three animals were randomly selected to serve as controls, while the remaining animals were designated as the study group and treated with CITT. In the CITT group, the partially resected tumour and margins were thermally ablated. In the control group the tumour was partially resected to simulate positive margins. The animals were monitored for up to 12 weeks. At the endpoint, the animals were sacrificed, and whole-body diagnostic necropsy was conducted immediately.

RESULTS

Recurrences and metastatic lesions were observed in iliac and popliteal lymph nodes and abdomens of all control animals. In contrast, the observed rate of recurrence and metastatic lesion was 0% among CITT-treated animals, significantly less than the >or=50% null-hypothesis rate expected upon treatment failure (exact binomial P = 0.0002). Complete histopathological healing was obtained in 2 of 12 rabbits, and residual inflammation remained at the ablation site up to 12 weeks post-ablation in 10 of 12 rabbits. This pattern of necrosis and inflammatory response was not observed in any of the control rabbits.

CONCLUSIONS

The CITT device effectively ablated partially resected VX2 carcinoma in a rabbit model, and inhibited recurrence and metastasis in this model. CITT evoked an inflammatory response that may be linked to the mechanism involved in reduced metastatic spread.

Conductive interstitial thermal therapy device for surgical margin ablation:In vivoverification of a theoretical model

PURPOSE

To demonstrate the efficacy and predictability of a new conductive interstitial thermal therapy (CITT) device to ablate surgical margins.

METHOD

The temperature distributions during thermal ablation of CITT were calculated with finite element modelling in a geometrical representation of perfused tissue. The depth of ablation was derived using the Arrhenius and the Sapareto and Dewey (S&D) models for the temperature range of 90 to 150 degrees C. The female pig animal model was used to test the validity of the mathematical model. Breast tissues were ablated to temperatures in the range of 79-170 degrees C, in vivo. Triphenyltetrazolium chloride viability stain was used to delineate viable tissue from ablated regions and the ablation depths were measured using digital imaging.

RESULTS

The calculations suggest that the CITT can be used to ablate perfused tissues to a 10-15 mm width within 20 minutes. The measured and calculated depths of ablation were statistically equivalent (99% confidence intervals) within +/- 1mm at 170 degrees C. At lower temperatures the equivalence between the model and the observations was within +/- 2 mm.

CONCLUSION

The CITT device can reliably and uniformly ablate a 10-15 mm wide region of soft tissue. Thus, it can be used to secure negative margins following the resection of a primary tumor, which could impede local recurrences in the treatment of local diseases such as early staged, non-metastatic, breast cancer.

Modeling the temperature dependence of thermophysical properties: Study on the effect of temperature dependence for RFA

Radio frequency ablation (RFA) has increasingly been used over the past few years and RFA treatment is minimally invasive for patients. However, it is difficult for operators to control the precise formation of coagulation zones due to inadequate imaging modalities. With this in mind, an ablation system using numerical simulation to analyze the temperature distribution of the organ is needed to overcome this deficiency. The objective of our work is to develop a temperature dependent thermophysical liver model. First, an overview is given of the development of the thermophysical liver model. Second, a simulation to evaluate the effect of temperature dependence of the thermophysical properties of the liver is explained. Finally, the result of the simulation, which indicated that the temperature dependence of thermophysical properties accounts for temperature differences influencing the accuracy of RFA treatment is described.