Numerical investigation of the effect of vessel size and distance on the cryosurgery of an adjacent tumor

Estimation of the injection criteria for magnetic hyperthermia therapy based on tumor morphology

Intratumoral multi-injection strategy enhances the efficacy of magnetic nanoparticle hyperthermia therapy (MNPH). In this study, criteria for the selection of injections and their location depending on the tumor shape/geometry are developed. The developed strategy is based on the thermal dosimetry results of different invasive 3D tumor models during MNPH simulation. MNPH simulations are conducted on physical tumor tissue models encased within healthy tissue. The tumor shapes are geometrically divided into a central tumor region containing maximum tumor volume and a peripheral tumor portion protruding in any random direction. The concepts of core and invasive radius are used to geometrically divide the tumor volume. Primary & secondary injections are used to inject MNP fluid into these respective tumor regions based on the invasiveness of the tumor. The optimization strategy is devised based on the zone of influence of primary & secondary injection. Results indicate that the zone of influence of secondary injection lies between 0.7 and 0.8 times the radial distance between the center of the tumor core and branch node point (extreme far endpoint on the invasive tumor surface). Additionally, the multi-injection strategy is more effective when the protrusion volume exceeds10%of the total volume. The proposed algorithm is used to devise multi-injection strategies for arbitrarily shaped tumors and will assist in pre-planning magnetic nanoparticle hyperthermia therapy.

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.

Numerical Simulation of bio-heat transfer for cryoablation of regularly shaped tumours in liver tissue using multiprobes

The present study focuses upon pre-operative planning strategy for cryosurgical treatment of multiple regularly shaped tumours inside a three dimensional liver tissue. Numerical simulations provide an optimal framework to predict the number of cryo-probes, their placement, operation time and thermal necrosis to the tumour and surrounding healthy tissues. An efficient cryosurgery process requires keeping the tumour cell under lethal temperature which is between -40 °C to -50 °C. The freezing process of undesired tumour tissues involves phase transition from liquid phase to solid phase, the accurate capturing of transition front and size or location of ice balls generated in the process are the important factors of cryosurgery. In this study, fixed domain heat capacity method has been utilized to take into account the latent heat of phase change in bio-heat transfer equation. The ice balls generated with different number of probes haven been analysed. Numerical simulations have been carried out using standard Finite Element Method with COMSOL 5.5 and results obtained are validated with previous studies.

Effects of injection rates and tissue diffusivity in magnetic nano-particle hyperthermia

Effects of injection rate and tumor physiology on the diffusion of magnetic nano-particles (MNPs) and temperature profile during magnetic hyperthermia are investigated in this work. The study considers three injection rates (2.5 μL/min, 10 μL/min, and 40 μL/min), and two MNP diffusion coefficients (10^-9 m²/s and 10^-11 m²/s). The simulation of this physics has been done on 3D tumor surrounded by healthy tissue. Transient MNP distribution in tissue is evaluated using Darcy's flow model and the MNP transport (convection-diffusion) equation. The temperature profile in the tumor model is computed by solving Penne's bioheat transfer equation (PBHTE). Results show tumors with high collagen content (with low MNP diffusivity) are more restrictive towards MNP transport than tumors having low collagen content. Thus, tumors with low MNP diffusivity need a higher injection rate to increase the homogeneity of MNP concentration as well as temperature profile during thermo-therapy. Results also show that, MNP fluid injected with a higher injection rate produces a more uniform MNP concentration up to greater depth than the lower injection rate.

In Silico Study on Tumor-Size-Dependent Thermal Profiles inside an Anthropomorphic Female Breast Phantom Subjected to Multi-Dipole Antenna Array

Electromagnetic hyperthermia as a potent adjuvant for conventional cancer therapies can be considered valuable in modern oncology, as its task is to thermally destroy cancer cells exposed to high-frequency electromagnetic fields. Hyperthermia treatment planning based on computer in silico simulations has the potential to improve the localized heating of breast tissues through the use of the phased-array dipole applicators. Herein, we intended to improve our understanding of temperature estimation in an anatomically accurate female breast phantom embedded with a tumor, particularly when it is exposed to an eight-element dipole antenna matrix surrounding the breast tissues. The Maxwell equations coupled with the modified Pennes’ bioheat equation was solved in the modelled breast tissues using the finite-difference time-domain (FDTD) engine. The microwave (MW) applicators around the object were modelled with shortened half-wavelength dipole antennas operating at the same 1 GHz frequency, but with different input power and phases for the dipole sources. The total input power of an eight-dipole antenna matrix was set at 8 W so that the temperature in the breast tumor did not exceed 42 °C. Finding the optimal setting for each dipole antenna from the matrix was our primary objective. Such a procedure should form the basis of any successful hyperthermia treatment planning. We applied the algorithm of multi for multi-objective optimization for the power and phases for the dipole sources in terms of maximizing the specific absorption rate (SAR) parameter inside the breast tumor while minimizing this parameter in the healthy tissues. Electro-thermal simulations were performed for tumors of different radii to confirm the reliable operation of the given optimization procedure. In the next step, thermal profiles for tumors of various sizes were calculated for the optimal parameters of dipole sources. The computed results showed that larger tumors heated better than smaller tumors; however, the procedure worked well regardless of the tumor size. This verifies the effectiveness of the applied optimization method, regardless of the various stages of breast tumor development.

GPU-based 3D iceball modeling for fast cryoablation simulation and planning

PURPOSE

The elimination of abdominal tumors by percutaneous cryoablation has been shown to be an effective and less invasive alternative to open surgery. Cryoablation destroys malignant cells by freezing them with one or more cryoprobes inserted into the tumor through the skin. Alternating cycles of freezing and thawing produce an enveloping iceball that causes the tumor necrosis. Planning such a procedure is difficult and time-consuming, as it is necessary to plan the number and cryoprobe locations and predict the iceball shape which is also influenced by the presence of heating sources, e.g., major blood vessels and warm saline solution, injected to protect surrounding structures from the cold.

METHODS

This paper describes a method for fast GPU-based iceball modeling based on the simulation of thermal propagation in the tissue. Our algorithm solves the heat equation within a cube around the cryoprobes tips and accounts for the presence of heating sources around the iceball.

RESULTS

Experimental results of two studies have been obtained: an ex vivo warm gel setup and simulation on five retrospective patient cases of kidney tumors cryoablation with various levels of complexity of the vascular structure and warm saline solution around the tumor tissue. The experiments have been conducted in various conditions of cube size and algorithm implementations. Results show that it is possible to obtain an accurate result within seconds.

CONCLUSION

The promising results indicate that our method yields accurate iceball shape predictions in a short time and is suitable for surgical planning.

Generalized solution and estimation method for cooling performance of downscaled cryoprobe

In cryosurgery, downscaling of cryoprobes is important to minimize surgical invasion. In this study, a set of analytical solutions to the freezing phenomenon around a cryoprobe in a dimensionless form is derived and the general trend is discussed to clarify the relationship between the freezing ability of a biological tissue and the cooling power of a cryoprobe. A one-dimensional axisymmetric model in the steady-state condition is considered. The relationship between the size of the frozen region, fluid temperature in the cryoprobe, and heat transfer coefficient on the wall of the cryoprobe in the dimensional form is derived under the condition mentioned above. The fluid temperature and heat transfer coefficient are eliminated from the solutions by introducing the steady-state cryoprobe surface temperature. This transformation indicates that the steady-state surface temperature directly affects the size of the frozen region and combination of fluid temperature and heat transfer coefficient occurs, which has the same cooling effect. The derived solutions are transformed into a dimensionless form using the characteristic length of bioheat transfer and normalizing the temperature distribution in an unfrozen tissue. The applicability of these analytical solutions is evaluated by comparing them with numerical simulation results from existing studies. The dimensionless solutions describe the general trend of the relationship between the frozen region and the cooling power of a cryoprobe, which is independent of the type of organ, fluid temperature, and heat transfer coefficient. Finally, the concept of freezing limit is established using the derived solutions. The freezing limit describes the minimum requirements to freeze a tissue, and it can be used as guideline to design future downscaled cryoprobes with a suitable cooling mechanism.