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

Experimental study for in vitro prostate cancer treatment with microwave ablation and pulsed electromagnetic field

This paper presents the findings of a comprehensive study exploring the synergistic effects arising from the combination of microwave ablation and pulsed electromagnetic field (PEMF) therapy on prostate cancer cells. The research encompassed five distinct experimental groups, with continuous electric field measurements conducted during the entire treatment process. Group 1 and Group 2, subjected to microwave power below 350 W, exhibited specific electric field values of 72,800 V/m and 56,600 V/m, respectively. In contrast, Group 3 and Group 4, exposed to 80 W microwave power, displayed electric field levels of approximately 1450 V/m, while remaining free from any observable electrical discharges. The migratory and invasive capacities of PC3 cells were assessed through a scratch test in all groups. Notably, cells in Group 3 and Group 4, subjected to the combined treatment of microwave ablation and PEMF, demonstrated significantly accelerated migration in comparison to those in Groups 1 and 2. Additionally, Group 5 cells, receiving PEMF treatment in isolation, exhibited decreased migratory ability. These results strongly suggest that the combined approach of microwave ablation and PEMF holds promise as a potential therapeutic intervention for prostate cancer, as it effectively reduced cell viability, induced apoptosis, and impeded migration ability in PC3 cells. Moreover, the isolated use of PEMF demonstrated potential in limiting migratory capacity, which could hold critical implications in the fight against cancer metastasis.

Particle swarm optimization solution for roll-off control in radiofrequency ablation of liver tumors: Optimal search for PID controller tuning

The study investigates the efficacy of a bioinspired Particle Swarm Optimization (PSO) approach for PID controller tuning in Radiofrequency Ablation (RFA) for liver tumors. Ex-vivo experiments were conducted, yielding a 9th order continuous-time transfer function. PSO was applied to optimize PID parameters, achieving outstanding simulation results: 0.605% overshoot, 0.314 seconds rise time, and 2.87 seconds settling time for a unit step input. Statistical analysis of 19 simulations revealed PID gains: Kp (mean: 5.86, variance: 4.22, standard deviation: 2.05), Ki (mean: 9.89, variance: 0.048, standard deviation: 0.22), Kd (mean: 0.57, variance: 0.021, standard deviation: 0.14) and ANOVA analysis for the 19 experiments yielded a p-value ≪ 0.05. The bioinspired PSO-based PID controller demonstrated remarkable potential in mitigating roll-off effects during RFA, reducing the risk of incomplete tumor ablation. These findings have significant implications for improving clinical outcomes in hepatocellular carcinoma management, including reduced recurrence rates and minimized collateral damage. The PSO-based PID tuning strategy offers a practical solution to enhance RFA effectiveness, contributing to the advancement of radiofrequency ablation techniques.

Using Patient-Specific 3D Modeling and Simulations to Optimize Microwave Ablation Therapy for Liver Cancer

Microwave ablation (MWA) of liver tumors presents challenges like under- and over-ablation, potentially leading to inadequate tumor destruction and damage to healthy tissue. This study aims to develop personalized three-dimensional (3D) models to simulate MWA for liver tumors, incorporating patient-specific characteristics. The primary objective is to validate the predicted ablation zones compared to clinical outcomes, offering insights into MWA before therapy to facilitate accurate treatment planning. Contrast-enhanced CT images from three patients were used to create 3D models. The simulations used coupled electromagnetic wave propagation and bioheat transfer to estimate the temperature distribution, predicting tumor destruction and ablation margins. The findings indicate that prolonged ablation does not significantly improve tumor destruction once an adequate margin is achieved, although it increases tissue damage. There was a substantial overlap between the clinical ablation zones and the predicted ablation zones. For patient 1, the Dice score was 0.73, indicating high accuracy, with a sensitivity of 0.72 and a specificity of 0.76. For patient 2, the Dice score was 0.86, with a sensitivity of 0.79 and a specificity of 0.96. For patient 3, the Dice score was 0.8, with a sensitivity of 0.85 and a specificity of 0.74. Patient-specific 3D models demonstrate potential in accurately predicting ablation zones and optimizing MWA treatment strategies.

Effect of high blood flow on heat distribution and ablation zone during microwave ablation-numerical approach

Microwave ablation has become a viable alternative for cancer treatment for patients who cannot undergo surgery. During this procedure, a single-slot coaxial antenna is employed to effectively deliver microwave energy to the targeted tissue. The success of the treatment was measured by the amount of ablation zone created during the ablation procedure. The significantly large blood vessel placed near the antenna causes heat dissipation by convection around the blood vessel. The heat sink effect could result in insufficient ablation, raising the risk of local tumor recurrence. In this study, we investigated the heat loss due to large blood vessels and the relationship between blood velocity and temperature distribution. The hepatic artery, with a diameter of 4 mm and a height of 50 mm and two branches, is considered in the computational domain. The temperature profile, localized tissue contraction, and ablation zones were simulated for initial blood velocities 0.05, 0.1, and 0.16 m/s using the 3D Pennes bio-heat equation, temperature-time dependent model, and cell death model, respectively. Temperature-dependent blood velocity is modeled using the Navier-Stokes equation, and the fluid-solid interaction boundary is treated as a convective boundary. For discretization, we utilizedH curl Ω $$ H\left(\operatorname{curl},\Omega \right) $$ elements for the wave propagation model,H 1 Ω $$ {H}^1\left(\Omega \right) $$ elements for the Pennes bio-heat model, andH 1 Ω 3 × L 0 2 Ω $$ {\left({H}^1\left(\Omega \right)\right)}^3\times {L}_0^2\left(\Omega \right) $$ elements for the Navier-Stokes equation, whereΩ $$ \Omega $$ represents the computational domain. The simulated results show that blood vessels and blood velocity have a significant impact on temperature distribution, tissue contraction, and the volume of the ablation zone.

The influence of microwave ablation parameters on the positioning of trocar in different cancerous tissues: a numerical study

The present study analyzed the microwave ablation of cancerous tumors located in six major cancer-prone organs and estimated the significance of input power and treatment time parameters in the apt positioning of the trocar into the tissue during microwave ablation. The present study has considered a three-dimensional two-compartment tumour-embedded tissue model. FEA based COMSOL Multiphysics software with inbuilt bioheat transfer, electromagnetic waves, heat transfer in solids and fluids, and laminar flow physics has been used to obtain the numerical results. Based on the mortality rates caused by cancer, the present study has considered six major organs affected by cancer, viz. lung, breast, stomach/gastric, liver, liver (with colon metastasis), and kidney for MWA analysis. The input power (100 W) and ablation times (4 minutes) with apt and inapt positioning of the trocar have been considered to compare the ablation volume of various cancerous tissues. The present study addresses one of the major problems clinicians face, i.e. the proper placement of the trocar due to poor imaging techniques and human error, resulting in incomplete tumor ablation and increased surgical procedures. The highest values of the ablation region have been observed for the liver, colon metastatic liver and breast cancerous tissues compared with other organs at the same operating conditions.

Advances in screening hyperthermic nanomedicines in 3D tumor models

Hyperthermic nanomedicines are particularly relevant for tackling human cancer, providing a valuable alternative to conventional therapeutics. The early-stage preclinical performance evaluation of such anti-cancer treatments is conventionally performed in flat 2D cell cultures that do not mimic the volumetric heat transfer occurring in human tumors. Recently, improvements in bioengineered 3D in vitro models have unlocked the opportunity to recapitulate major tumor microenvironment hallmarks and generate highly informative readouts that can contribute to accelerating the discovery and validation of efficient hyperthermic treatments. Leveraging on this, herein we aim to showcase the potential of engineered physiomimetic 3D tumor models for evaluating the preclinical efficacy of hyperthermic nanomedicines, featuring the main advantages and design considerations under diverse testing scenarios. The most recent applications of 3D tumor models for screening photo- and/or magnetic nanomedicines will be discussed, either as standalone systems or in combinatorial approaches with other anti-cancer therapeutics. We envision that breakthroughs toward developing multi-functional 3D platforms for hyperthermia onset and follow-up will contribute to a more expedited discovery of top-performing hyperthermic therapies in a preclinical setting before their in vivo screening.

MATLAB-based innovative 3D finite element method simulator for optimized real-time hyperthermia analysis

BACKGROUND AND OBJECTIVE

Owing to the significant role of hyperthermia in enhancing the efficacy of chemotherapy or radiotherapy for treating malignant tissues, this study introduces a real-time hyperthermia simulator (RTHS) based on the three-dimensional finite element method (FEM) developed using the MATLAB App Designer.

METHODS

The simulator consisted of operator-defined homogeneous and heterogeneous phantom models surrounded by an annular phased array (APA) of eight dipole antennas designed at 915 MHz. Electromagnetic and thermal analyses were conducted using the RTHS. To locally raise the target temperature according to the tumor's location, a convex optimization algorithm (COA) was employed to excite the antennas using optimal values of the phases to maximize the electric field at the tumor and amplitudes to achieve the required temperature at the target position. The performance of the proposed RTHS was validated by comparing it with similar hyperthermia setups in the FEM-based COMSOL software and finite-difference time-domain (FDTD)-based Sim4Life software.

RESULTS

The simulation results obtained using the RTHS were consistent, both for the homogeneous and heterogeneous models, with those obtained using commercially available tools, demonstrating the reliability of the proposed hyperthermia simulator. The effectiveness of the simulator was illustrated for target positions in five different regions for both homogeneous and heterogeneous phantom models. In addition, the RTHS was cost-effective and consumed less computational time than the available software. The proposed method achieved 94% and 96% accuracy for element sizes of λ/26 and λ/36, respectively, for the homogeneous model. For the heterogeneous model, the method demonstrated 93% and 95% accuracy for element sizes of λ/26 and λ/36, respectively. The accuracy can be further improved by using a more refined mesh at the cost of a higher computational time.

CONCLUSIONS

The proposed hyperthermia simulator demonstrated reliability, cost-effectiveness, and reduced computational time compared to commercial software, making it a potential tool for optimizing hyperthermia treatment.

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.

Effects of Pulsed Radiofrequency Source on Cardiac Ablation

Heart arrhythmia is caused by abnormal electrical conduction through the myocardium, which in some cases, can be treated with heat. One of the challenges is to reduce temperature peaks-by still guaranteeing an efficient treatment where desired-to avoid any healthy tissue damage or any electrical issues within the device employed. A solution might be employing pulsed heat, in which thermal dose is given to the tissue with a variation in time. In this work, pulsed heat is used to modulate induced temperature fields during radiofrequency cardiac ablation. A three-dimensional model of the myocardium, catheter and blood flow is developed. Porous media, heat conduction and Navier-Stokes equations are, respectively, employed for each of the investigated domains. For the electric field, solved via Laplace equation, it is assumed that the electrode is at a fixed voltage. Pulsed heating effects are considered with a cosine time-variable pulsed function for the fixed voltage by constraining the product between this variable and time. Different dimensionless frequencies are considered and applied for different blood flow velocity and sustained voltages. Results are presented for different pulsed conditions to establish if a reasonable ablation zone, known from the obtained temperature profiles, can be obtained without any undesired temperature peaks.

Computer simulations of consecutive radiofrequency pulses applied at the same point during cardiac catheter ablation: Implications for lesion size and risk of overheating

BACKGROUND AND OBJECTIVES

To study temperature distribution and lesion size during two repeated radiofrequency (RF) pulses applied at the same point in the context of RF cardiac ablation (RFCA).

METHODS

An in-silico RFCA model accounting for reversible and irreversible changes in myocardium electrical properties due to RF-induced heating. Arrhenius damage model to estimate lesion size during the application of two 20 W pulses at intervals (INT) of from 5 to 70 s. We considered two pulse durations: 20 s and 30 s.

RESULTS

INT has a significant effect on lesion size and maximum tissue temperature (TMAX). The shorter the INT the greater the increase in lesion size after the second pulse but also the greater the TMAX. If the second pulse is applied almost immediately (INT=5 s), depth increases 1.4 mm and 1.5 mm for pulses of 20 s and 30 s, respectively. If INT is longer than 30 s it increases 1.1 mm and 1.3 mm for pulses of 20 s and 30 s, respectively. While a single 20 s pulse causes T_MAX=79 ºC, a second pulse produces values of from 92 to 96 ºC (the higher the temperature the shorter the INT). For 30 s pulses, T_MAX=93 ºC for a single pulse, and varied from 98 to 104 ºC for a second pulse.

CONCLUSIONS

Applying a second RF pulse at the same ablation site increases lesion depth by 1 - 1.5 mm more than a single pulse and could lead to higher temperatures (up to 17 ºC). Both lesion depth and maximum tissue temperature increased at shorter inter-pulse intervals, which could cause clinical complications from overheating such as steam pops.

Mathematical modeling of microwave liver ablation with a variable-porosity medium approach

BACKGROUND AND OBJECTIVES

Thermal ablation of tumors plays a key role to fight cancer, since it is a minimally invasive treatment which involves some advantages compared to surgery and chemotherapy, such as shorter hospital stays and consequently lower costs, along with minor side effects. In this context, computational modeling of heat transfer during thermal ablation is relevant to accurately predict the obtained ablation zone in order to avoid tumor recurrence risk caused by incomplete ablation, and the same time to save the surrounding healthy tissue. The aim of this work is to develop a more realistic porous media-based mathematical model to simulate a microwave thermal ablation (MWA) of an in vivo liver tumor surrounded by healthy tissue.

METHODS

The domain is made up of a spherical tumor bounded by a cylindrical healthy liver tissue. The simulated microwave antenna is a 14 G HS Amica-Gen Probe, and the supplied power of 60 W is applied for 300 s and 600 s. The model consists in coupling modified Local Thermal Non Equilibrium (LTNE) equations with the electromagnetic equations. The LTNE equations include a variable porosity function which fits the porosity changing from the tumor core to the rim based on experimental measures in in vivo cases. Moreover, four different blood vessels' uniform distributions are investigated to compare the effects of different vascularizations of the considered target tissue.

RESULTS

The results are shown in terms of temperature fields, ablation diameters and volumes based on the Arrhenius thermal damage model with 99% of cell death probability. The outcomes show a very good agreement with a clinical study on human patients with hepatocellular carcinoma using the same antenna and energy setting, when terminal arteries distribution is included.

CONCLUSIONS

In this work, an in vivo microwave ablation of liver tumor surrounded by healthy tissue is modeled with a variable-porosity medium approach based on experimental measures. The outcomes shown for distinct vascularizations underline the key relevance of modeling more and more accurately tumor MWA, by considering increasingly realistic features, avoiding tumor recurrence, and improving both medical protocols and devices.

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.

A Computational Study on Magnetic Nanoparticles Hyperthermia of Ellipsoidal Tumors

The modelling of magnetic hyperthermia using nanoparticles of ellipsoid tumor shapes has not been studied adequately. To fill this gap, a computational study has been carried out to determine two key treatment parameters: the therapeutic temperature distribution and the extent of thermal damage. Prolate and oblate spheroidal tumors, of various aspect ratios, surrounded by a large healthy tissue region are assumed. Tissue temperatures are determined from the solution of Pennes’ bio-heat transfer equation. The mortality of the tissues is determined by the Arrhenius kinetic model. The computational model is successfully verified against a closed-form solution for a perfectly spherical tumor. The therapeutic temperature and the thermal damage in the tumor center decrease as the aspect ratio increases and it is insensitive to whether tumors of the same aspect ratio are oblate or prolate spheroids. The necrotic tumor area is affected by the tumor prolateness and oblateness. Good comparison is obtained of the present model with three sets of experimental measurements taken from the literature, for animal tumors exhibiting ellipsoid-like geometry. The computational model enables the determination of the therapeutic temperature and tissue thermal damage for magnetic hyperthermia of ellipsoidal tumors. It can be easily reproduced for various treatment scenarios and may be useful for an effective treatment planning of ellipsoidal tumor geometries.

Fluid–Structure Interaction and Non-Fourier Effects in Coupled Electro-Thermo-Mechanical Models for Cardiac Ablation

In this study, a fully coupled electro-thermo-mechanical model of radiofrequency (RF)-assisted cardiac ablation has been developed, incorporating fluid–structure interaction, thermal relaxation time effects and porous media approach. A non-Fourier based bio-heat transfer model has been used for predicting the temperature distribution and ablation zone during the cardiac ablation. The blood has been modeled as a Newtonian fluid and the velocity fields are obtained utilizing the Navier–Stokes equations. The thermal stresses induced due to the heating of the cardiac tissue have also been accounted. Parametric studies have been conducted to investigate the effect of cardiac tissue porosity, thermal relaxation time effects, electrode insertion depths and orientations on the treatment outcomes of the cardiac ablation. The results are presented in terms of predicted temperature distributions and ablation volumes for different cases of interest utilizing a finite element based COMSOL Multiphysics software. It has been found that electrode insertion depth and orientation has a significant effect on the treatment outcomes of cardiac ablation. Further, porosity of cardiac tissue also plays an important role in the prediction of temperature distribution and ablation volume during RF-assisted cardiac ablation. Moreover, thermal relaxation times only affect the treatment outcomes for shorter treatment times of less than 30 s.