Computer modelling of an impedance-controlled pulsing protocol for RF tumour ablation with a cooled electrode

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).

Microwave ablation modeling with AMICA antenna: Validation by means a numerical analysis

BACKGROUND AND OBJECTIVES

Percutaneous microwave thermal ablation is based on electromagnetic waves that generate dielectric heating, and it is widely recognized as one of the mostly used techniques for tumor treatment. The aim of this work is to validate a predictive model capable of providing physicians with guidelines to be used during thermal ablation procedures avoiding collateral damage.

METHODS

A finite element commercial software, COMSOL Multiphysics, is employed to implement a tuning-parameter approach. Governing equations are written with reference to variable-porosity and Local Thermal Non-Equilibrium (LTNE) equations are employed. The simulations results are compared with available ex-vivo and in-vivo data with the help of regression analysis. For in-vivo data simulations, velocity vector modulus and direction are varied between 0.0007 and 0.0009 m/s and 90-270°, respectively, in order to use this parameter as a tuning one to simulate - and lately optimize with respect to the differences from experimental outcomes - all the possible directions of the blood flow with respect to the antenna, whose insertion angle is not registered in the dataset.

RESULTS

The model is validated using reference data provided by the manufacturer (AMICA), which is obtained from ex-vivo bovine liver. The model accurately predicts the size and shape of the ablated area, resulting in an overestimation lesser than 10 %. Additionally, predictive data are compared to an in-vivo dataset. The ablated volume is accurately predicted with a mean underestimation of 6 %. The sphericity index is calculated as 0.75 and 0.62 for the predictions and in-vivo data, respectively.

CONCLUSION

This study developed a predictive model for microwave ablation of liver tumors that showed good performance in predicting ablation dimensions and sphericity index for ex-vivo bovine liver and for in-vivo human liver data with the tuning technique. The study emphasizes the necessity for additional development and validation to enhance the accuracy and reliability of in-vivo application.

Gold nanorods assisted photothermal therapy of bladder cancer in mice: A computational study on the effects of gold nanorods distribution at the centre, periphery, and surface of bladder cancer

BACKGROUND AND OBJECTIVES

Gold nanorod-assisted photothermal therapy (GNR-PTT) is a cancer treatment whereby GNRs incorporated into the tumour act as photo-absorbers to elevate the thermal destruction effect. In the case of bladder, there are few possible routes to target the tumour with GNRs, namely peri/intra-tumoural injection and intravesical instillation of GNRs. These two approaches lead to different GNR distribution inside the tumour and can affect the treatment outcome.

METHODOLOGY

The present study investigates the effects of heterogeneous GNR distribution in a typical setup of GNR-PTT. Three cases were considered. Case 1 considered the GNRs at the tumour centre, while Case 2 represents a hypothetical scenario where GNRs are distributed at the tumour periphery; these two cases represent intratumoural accumulation with different degree of GNR spread inside the tumour. Case 3 is achieved when GNRs target the exposed tumoural surface that is invading the bladder wall, when they are delivered by intravesical instillation.

RESULTS

Results indicate that for a laser power of 0.6 W and GNR volume fraction of 0.01%, Case 2 and 3 were successful in achieving complete tumour eradication after 330 and 470 s of laser irradiation, respectively. Case 1 failed to form complete tumour damage when the GNRs are concentrated at the tumour centre but managed to produce complete tumour damage if the spread of GNRs is wider. Results from Case 2 also demonstrated a different heating profile from Case 1, suggesting that thermal ablation during GNR-PTT is dependant on the GNRs distribution inside the tumour. Case 3 shows similar results to Case 2 whereby gradual but uniform heating is observed. Cases 2 and 3 show that uniformly heating the tumour can reduce damage to the surrounding tissues.

CONCLUSIONS

Different GNR distribution associated with the different methods of introducing GNRs to the bladder during GNR-PTT affect the treatment outcome of bladder cancer in mice. Insufficient spreading during intratumoural injection of GNRs can render the treatment ineffective, while administered via intravesical instillation. GNR distribution achieved through intravesical instillation present some advantages over intratumoural injection and is worthy of further exploration.

An in silico assessment on the potential of using saline infusion to overcome non-confluent coagulation zone during two-probe, no-touch bipolar radiofrequency ablation of liver cancer

No-touch bipolar radiofrequency ablation (bRFA) is known to produce incomplete tumour ablation with a 'butterfly-shaped' coagulation zone when the interelectrode distance exceeds a certain threshold. Although non-confluent coagulation zone can be avoided by not implementing the no-touch mode, doing so exposes the patient to the risk of tumour track seeding. The present study investigates if prior infusion of saline into the tissue can overcome the issues of non-confluent or butterfly-shaped coagulation. A computational modelling approach based on the finite element method was carried out. A two-compartment model comprising the tumour that is surrounded by healthy liver tissue was developed. Three cases were considered; i) saline infusion into the tumour centre; ii) one-sided saline infusion outside the tumour; and iii) two-sided saline infusion outside the tumour. For each case, three different saline volumes were considered, i.e. 6, 14 and 22 ml. Saline concentration was set to 15% w/v. Numerical results showed that saline infusion into the tumour centre can overcome the butterfly-shaped coagulation only if the infusion volume is sufficient. On the other hand, one-sided infusion outside the tumour did not overcome this. Two-sided infusion outside the tumour produced confluent coagulation zone with the largest volume. Results obtained from the present study suggest that saline infusion, when carried out correctly, can be used to effectively eradicate liver cancer. This presents a practical solution to address non-confluent coagulation zone typical of that during two-probe bRFA treatment.

An in silico derived dosage and administration guide for effective thermochemical ablation of biological tissues with simultaneous injection of acid and base

BACKGROUND AND OBJECTIVES

Thermochemical ablation (TCA) is a thermal ablation technique involving the injection of acid and base, either sequentially or simultaneously, into the target tissue. TCA remains at the conceptual stage with existing studies unable to provide recommendations on the optimum injection rate, and reagent concentration and volume. Limitations in current experimental methodology have prevented proper elucidation of the thermochemical processes inside the tissue during TCA. Nevertheless, the computational TCA framework developed recently by Mak et al. [Mak et al., Computers in Biology and Medicine, 2022, 145:105494] has opened new avenues in the development of TCA. Specifically, a recommended safe dosage is imperative in driving TCA research beyond the conceptual stage.

METHODS

The aforesaid computational TCA framework for sequential injection was applied and adapted to simulate TCA with simultaneous injection of acid and base at equimolar and equivolume. The developed framework, which describes the flow of acid and base, their neutralisation, the rise in tissue temperature and the formation of thermal damage, was solved numerically using the finite element method. The framework will be used to investigate the effects of injection rate, reagent concentration, volume and type (weak/strong acid-base combination) on temperature rise and thermal coagulation formation.

RESULTS

A higher injection rate resulted in higher temperature rise and larger thermal coagulation. Reagent concentration of 7500 mol/m³ was found to be optimum in producing considerable thermal coagulation without the risk of tissue overheating. Thermal coagulation volume was found to be consistently larger than the total volume of acid and base injected into the tissue, which is beneficial as it reduces the risk of chemical burn injury. Three multivariate second-order polynomials that express the targeted coagulation volume as functions of injection rate and reagent volume, for the weak-weak, weak-strong and strong-strong acid-base combinations were also derived based on the simulated data.

CONCLUSIONS

A guideline for a safe and effective implementation of TCA with simultaneous injection of acid and base was recommended based on the numerical results of the computational model developed. The guideline correlates the coagulation volume with the reagent volume and injection rate, and may be used by clinicians in determining the safe dosage of reagents and optimum injection rate to achieve a desired thermal coagulation volume during TCA.

Radiofrequency ablation for liver tumors abutting complex blood vessel structures: treatment protocol optimization using response surface method and computer modeling

OBJECTIVE

To achieve a result of a large tumor ablation volume with minimal thermal damage to the surrounding blood vessels by designing a few clinically-adjustable operating parameters in radiofrequency ablation (RFA) for liver tumors abutting complex vascular structures.

METHODS

Response surface method (RSM) was employed to correlate the ablated tumor volume (Ra) and thermal damage to blood vessels (Dt) based on RFA operating parameters: ablation time, electrode position, and insertion angle. A coupled electric-thermal-fluid RFA computer model was created as the testbed for RSM to simulate RFA process. Then, an optimal RFA protocol for the two conflicting goals, namely (1) large tumor ablation and (2) small thermal damage to the surrounding blood vessels, has been achieved under a specific ablation environment.

RESULTS

Linear regression analysis confirmed that the RFA protocol significantly affected Ra and Dt (the adjusted coefficient of determination Radj2 = 93.61% and 95.03%, respectively). For a proposed liver tumor scenario (liver tumor with a dimension of 4×3×2.9 cm³ abutting a complex vascular structure), an optimized RFA protocol was found based on the regression results in RSM. Compared with a reference RFA protocol, in which the electrode was centered in the tumor with a 12-min ablation time, the optimized RFA protocol has increased Ra  from 98.1% to 99.6% and decreased Dt from 4.1% to 0.4%, achieving nearly the complete ablation of proposed liver tumor and ignorable thermal damages to vessels.

CONCLUSION

This work showed that it is possible to design a few clinically-adjustable operating parameters of RFA for achieving a large tumor ablation volume while minimizing thermal damage to the surrounding blood vessels.

A New Thermal Damage-Controlled Protocol for Thermal Ablation Modeled with Modified Porous Media-Based Bioheat Equation with Variable Porosity

Thermal ablation of tumors is a minimally invasive technique more and more employed in cancer treatments. The main shortcomings of this technique are, on the one hand, the risk of an incomplete ablation, and on the other hand, the destruction of the surrounding healthy tissue. In this work, thermal ablation of a spherical hepatocellular carcinoma tumor (HCC) surrounded by healthy tissue is modeled. A modified porous media-based bioheat model is employed, including porosity variability from tumor core to healthy tissue, following experimental in vivo measures. Moreover, three different protocols are investigated: a constant heating protocol, a pulsating protocol, and a new developed damage-controlled protocol. The proposed damage-controlled protocol changes the heating source from constant to pulsating according to the thermal damage probability on the tumor rim. The equations are numerically solved by means of the commercial software COMSOL Multiphysics, and the outcomes show that the new proposed protocol is able to achieve the complete ablation in less time than the completely pulsating protocol, and to reach tissue temperature on the tumor rim 10 °C smaller than the constant protocol. These results are relevant to develop and improve different patient-based and automated protocols which can be embedded in medical devices’ software or in mobile applications, supporting medical staff with innovative technical solutions.

How does saline backflow affect the treatment of saline-infused radiofrequency ablation?

BACKGROUND AND OBJECTIVE

Saline infusion is applied together with radiofrequency ablation (RFA) to enlarge the ablation zone. However, one of the issues with saline-infused RFA is backflow, which spreads saline along the insertion track. This raises the concern of not only thermally ablating the tissue within the backflow region, but also the loss of saline from the targeted tissue, which may affect the treatment efficacy.

METHODS

In the present study, 2D axisymmetric models were developed to investigate how saline backflow influence saline-infused RFA and whether the aforementioned concerns are warranted. Saline-infused RFA was described using the dual porosity-Joule heating model. The hydrodynamics of backflow was described using Poiseuille law by assuming the flow to be similar to that in a thin annulus. Backflow lengths of 3, 4.5, 6 and 9 cm were considered.

RESULTS

Results showed that there is no concern of thermally ablating the tissue in the backflow region. This is due to the Joule heating being inversely proportional to distance from the electrode to the fourth power. Results also indicated that larger backflow lengths led to larger growth of thermal damage along the backflow region and greater decrease in coagulation volume. Hence, backflow needs to be controlled to ensure an effective treatment of saline-infused RFA.

CONCLUSIONS

There is no risk of ablating tissues around the needle insertion track due to backflow. Instead, the risk of underablation as a result of the loss of saline due to backflow was found to be of greater concern.

A numerical study to investigate the effects of tumour position on the treatment of bladder cancer in mice using gold nanorods assisted photothermal ablation

Gold nanorods assisted photothermal therapy (GNR-PTT) is a new cancer treatment technique that has shown promising potential for bladder cancer treatment. The position of the bladder cancer at different locations along the bladder wall lining can potentially affect the treatment efficacy since laser is irradiated externally from the skin surface. The present study investigates the efficacy of GNR-PTT in the treatment of bladder cancer in mice for tumours growing at three different locations on the bladder, i.e., Case 1: closest to skin surface, Case 2: at the bottom half of the bladder, and Case 3: at the side of the bladder. Investigations were carried out numerically using an experimentally validated framework for optical-thermal simulations. An in-silico approach was adopted due to the flexibility in placing the tumour at a desired location along the bladder lining. Results indicate that for the treatment parameters considered (laser power 0.3 W, GNR volume fraction 0.01% v/v), only Case 1 can be used for an effective GNR-PTT. No damage to the tumour was observed in Cases 2 and 3. Analysis of the thermo-physiological responses showed that the effectiveness of GNR-PTT in treating bladder cancer depends not only on the depth of the tumour from the skin surface, but also on the type of tissue that the laser must pass through before reaching the tumour. In addition, the results are reliant on GNRs with a diameter of 10 nm and an aspect ratio of 3.8 - tuned to exhibit peak absorption for the chosen laser wavelength. Results from the present study can be used to highlight the potential for using GNR-PTT for treatment of human bladder cancer. It appears that Cases 2 and 3 suggest that GNR-PTT, where the laser passes through the skin to reach the bladder, may be unfeasible in humans. While this study shows the feasibility of using GNRs for photothermal ablation of bladder cancer, it also identifies the current limitations needed to be overcome for an effective clinical application in the bladder cancer patients.

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.

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.

Quasi-distributed fiber optic sensor-based control system for interstitial laser ablation of tissue: theoretical and experimental investigations

This work proposes the quasi-distributed real-time monitoring and control of laser ablation (LA) of liver tissue. To confine the thermal damage, a pre-planning stage of the control strategy based on numerical simulations of the bioheat-transfer was developed to design the control parameters, then experimentally assessed. Fiber Bragg grating (FBG) sensors were employed to design the automatic thermometry system used for temperature feedback control for interstitial LA. The tissue temperature was maintained at a pre-set value, and the influence of different sensor locations (on the direction of the beam propagation and backward) on the thermal outcome was evaluated in comparison with the uncontrolled case. Results show that the implemented computational model was able to properly describe the temperature evolution of the irradiated tissue. Furthermore, the realized control strategy allowed for the accurate confinement of the laser-induced temperature increase, especially when the temperature control was actuated by sensors located in the direction of the beam propagation, as confirmed by the calculated fractions of necrotic tissues (e.g., 23 mm³ and 53 mm³ for the controlled and uncontrolled LA, respectively).

Numerical Investigation of a Thermal Ablation Porous Media-Based Model for Tumoral Tissue with Variable Porosity

Thermal ablation is a minimally or noninvasive cancer therapy technique that involves fewer complications, shorter hospital stays, and fewer costs. In this paper, a thermal-ablation bioheat model for cancer treatment is numerically investigated, using a porous media-based model. The main objective is to evaluate the effects of a variable blood volume fraction in the tumoral tissue (i.e., the porosity), in order to develop a more realistic model. A modified local thermal nonequilibrium model (LTNE) is implemented including the water content vaporization in the two phases separately and introducing the variable porosity in the domain, described by a quadratic function changing from the core to the rim of the tumoral sphere. The equations are numerically solved employing the finite-element commercial code COMSOL Multiphysics. Results are compared with the results obtained employing two uniform porosity values (ε = 0.07 and ε = 0.23) in terms of coagulation zones at the end of the heating period, maximum temperatures reached in the domain, and temperature fields and they are presented for different blood vessels. The outcomes highlight how important is to predict coagulation zones achieved in thermal ablation accurately. In this way, indeed, incomplete ablation, tumor recurrence, or healthy tissue necrosis can be avoided, and medical protocols and devices can be improved.

Pennes' bioheat equation vs. porous media approach in computer modeling of radiofrequency tumor ablation

The objective of this study was to compare three different heat transfer models for radiofrequency ablation of in vivo liver tissue using a cooled electrode and three different voltage levels. The comparison was between the simplest but less realistic Pennes' equation and two porous media-based models, i.e. the Local Thermal Non-Equilibrium (LTNE) equations and Local Thermal Equilibrium (LTE) equation, both modified to take into account two-phase water vaporization (tissue and blood). Different blood volume fractions in liver were considered and the blood velocity was modeled to simulate a vascular network. Governing equations with the appropriate boundary conditions were solved with Comsol Multiphysics finite-element code. The results in terms of coagulation transverse diameters and temperature distributions at the end of the application showed significant differences, especially between Pennes and the modified LTNE and LTE models. The new modified porous media-based models covered the ranges found in the few in vivo experimental studies in the literature and they were closer to the published results with similar in vivo protocol. The outcomes highlight the importance of considering the three models in the future in order to improve thermal ablation protocols and devices and adapt the model to different organs and patient profiles.

Numerical analysis of the pulsating heat source effects in a tumor tissue

BACKGROUND AND OBJECTIVES

Hyperthermia treatment is nowadays recognized as the fourth additional cancer therapy technique following surgery, chemotherapy, and radiation; it is a minimally or non-invasive technique which involves fewer complications, a shorter hospital stay, and fewer costs. In this paper, pulsating heat effects on heat transfer in a tumor tissue under hyperthermia are analyzed. The objective of the paper is to find and quantify the advantages of pulsatile heat protocols under different periodical heating schemes and for different tissue morphologies.

METHODS

The tumor tissue is modeled as a porous sphere made up of a solid phase (tissue, interstitial space, etc.) and a fluid phase (blood). A Local Thermal Non-Equilibrium (LTNE) model is employed to consider the local temperature difference between the two phases. Governing equations with the appropriate boundary conditions are solved with the finite-element code COMSOL Multiphysics®. The pulsating effect is modeled with references to a cosine function with different frequencies, and such different heating protocols are compared at equal delivered energy, i. e. different heating times at equal maximum power.

RESULTS

Different tissue properties in terms of blood vessels sizes and blood volume fraction in tissue (porosity) are investigated. The results are shown in terms of tissue temperature and percentage of necrotic tissue obtained. The most powerful result achieved using a pulsating heat source instead of a constant one is the decreasing of maximum temperature in any considered case, even reaching about 30% lower maximum temperatures. Furthermore, the evaluation of tissue damage at the end of treatment shows that pulsating heat allows to necrotize the same tumoral tissue area of the non-pulsating heat source.

CONCLUSIONS

Modeling pulsating heat protocols in thermal ablation under different periodical heating schemes and considering different tissues morphologies in a tumor tissue highlights how the application of pulsating heat sources allows to avoid high temperature peaks, and simultaneously to ablate the same tumoral area obtained with a non-pulsating heat source. This is a powerful result to improve medical protocols and devices in thermal ablation of tumors.

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.

Radiofrequency ablation combined with conductive fluid-based dopants (saline normal and colloidal gold): computer modeling and ex vivo experiments

BACKGROUND

The volume of the coagulation zones created during radiofrequency ablation (RFA) is limited by the appearance of roll-off. Doping the tissue with conductive fluids, e.g., gold nanoparticles (AuNPs) could enlarge these zones by delaying roll-off. Our goal was to characterize the electrical conductivity of a substrate doped with AuNPs in a computer modeling study and ex vivo experiments to investigate their effect on coagulation zone volumes.

METHODS

The electrical conductivity of substrates doped with normal saline or AuNPs was assessed experimentally on agar phantoms. The computer models, built and solved on COMSOL Multiphysics, consisted of a cylindrical domain mimicking liver tissue and a spherical domain mimicking a doped zone with 2, 3 and 4 cm diameters. Ex vivo experiments were conducted on bovine liver fragments under three different conditions: non-doped tissue (ND Group), 2 mL of 0.9% NaCl (NaCl Group), and 2 mL of AuNPs 0.1 wt% (AuNPs Group).

RESULTS

The theoretical analysis showed that adding normal saline or colloidal gold in concentrations lower than 10% only modifies the electrical conductivity of the doped substrate with practically no change in the thermal characteristics. The computer results showed a relationship between doped zone size and electrode length regarding the created coagulation zone. There was good agreement between the ex vivo and computational results in terms of transverse diameter of the coagulation zone.

CONCLUSIONS

Both the computer and ex vivo experiments showed that doping with AuNPs can enlarge the coagulation zone, especially the transverse diameter and hence enhance sphericity.

Role of saline concentration during saline-infused radiofrequency ablation: Observation of secondary Joule heating along the saline-tissue interface

Infusion of saline prior to radiofrequency ablation (RFA) is known to enlarge the thermal coagulation zone. The abundance of ions in saline elevate the electrical conductivity of the saline-saturated region. This promotes greater electric current flow inside the tissue, which increases the amount of RF energy deposition and subsequently enlarges the coagulation zone. In theory, infusion of higher concentration of saline should lead to larger coagulation zone due to the greater number of ions. Nevertheless, existing studies on the effects of concentration on saline-infused RFA have been conflicting, with the exact role of saline concentration yet to be fully elucidated. In this paper, computational models of saline-infused RFA were developed to investigate the role of saline concentration on the outcome of saline-infused RFA. The elevation in tissue electrical conductivity was modelled using the microscopic mixture model, while RFA was modelled using the coupled dual porosity-Joule heating model. Results obtained indicated that the presence of a concentration threshold to which no further elevation in tissue electrical conductivity and enlargement in thermal coagulation can occur. This threshold was determined to be at 15% NaCl. Analysis of the Joule heating distribution revealed the presence of a secondary Joule heating site located along the interface between wet and dry tissue. This secondary Joule heating was responsible for the enlargement in coagulation volume and its rapid growth phase during ablation.

How large is the periablational zone after radiofrequency and microwave ablation? Computer-based comparative study of two currently used clinical devices

Purpose:

To compare the size of the coagulation (CZ) and periablational (PZ) zones created with two commercially available devices in clinical use for radiofrequency (RFA) and microwave ablation (MWA), respectively.

Methods:

Computer models were used to simulate RFA with a 3-cm Cool-tip applicator and MWA with an Amica-Gen applicator. The Arrhenius model was used to compute the damage index (Ω). CZ was considered when Ω> 4.6 (>99% of damaged cells). Regions with 0.6<Ω< 2.1 were considered as the PZ (tissue that has undergone moderate sub-ablative hyperthermia). The ratio of PZ volume to CZ volume (PZ/CZ) was regarded as a measure of performance, since a low value implies achieving a large CZ while keeping the PZ small.

Results:

Ten-min RFA (51 W) created smaller periablational zones than 10-min MWA (11.3 cm³ vs. 17.2 22.9 cm³, for 60 100 W MWA, respectively). Prolonging duration from 5 to 10 min increased the PZ in MWA more than in RFA (2.7 cm³ for RFA vs. 8.3–11.9 cm³ for 60–100 W MWA, respectively). PZ/CZ for RFA were relatively high (65–69%), regardless of ablation time, while those for MWA were highly dependent on the duration (increase of up to 25% between 5 and 10 min) and on the applied power (smaller values as power was raised, 102% for 60 W vs. 81% for 100 W, both for 10 min). The lowest PZ/CZ across all settings was 56%, obtained with 100 W-5 min MWA.

Conclusions:

Although RFA creates smaller periablational zones than MWA, 100 W-5 min MWA provides the lowest PZ/CZ.

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.

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

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

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

How coagulation zone size is underestimated in computer modeling of RF ablation by ignoring the cooling phase just after RF power is switched off

All the numerical models developed for radiofrequency ablation so far have ignored the possible effect of the cooling phase (just after radiofrequency power is switched off) on the dimensions of the coagulation zone. Our objective was thus to quantify the differences in the minor radius of the coagulation zone computed by including and ignoring the cooling phase. We built models of RF tumor ablation with 2 needle-like electrodes: a dry electrode (5 mm long and 17G in diameter) with a constant temperature protocol (70°C) and a cooled electrode (30 mm long and 17G in diameter) with a protocol of impedance control. We observed that the computed coagulation zone dimensions were always underestimated when the cooling phase was ignored. The mean values of the differences computed along the electrode axis were always lower than 0.15 mm for the dry electrode and 1.5 mm for the cooled electrode, which implied a value lower than 5% of the minor radius of the coagulation zone (which was 3 mm for the dry electrode and 30 mm for the cooled electrode). The underestimation was found to be dependent on the tissue characteristics: being more marked for higher values of specific heat and blood perfusion and less marked for higher values of thermal conductivity.