Finite-element analysis of hepatic cryoablation around a large blood vessel

AI-Based Isotherm Prediction for Focal Cryoablation of Prostate Cancer

RATIONALE AND OBJECTIVES

Focal therapies have emerged as minimally invasive alternatives for patients with localized low-risk prostate cancer (PCa) and those with postradiation recurrence. Among the available focal treatment methods for PCa, cryoablation offers several technical advantages, including the visibility of the boundaries of frozen tissue on the intraprocedural images, access to anterior lesions, and the proven ability to treat postradiation recurrence. However, predicting the final volume of the frozen tissue is challenging as it depends on several patient-specific factors, such as proximity to heat sources and thermal properties of the prostatic tissue.

MATERIALS AND METHODS

This paper presents a convolutional neural network model based on 3D-Unet to predict the frozen isotherm boundaries (iceball) resultant from a given a cryo-needle placement. Intraprocedural magnetic resonance images acquired during 38 cases of focal cryoablation of PCa were retrospectively used to train and validate the model. The model accuracy was assessed and compared against a vendor-provided geometrical model, which is used as a guideline in routine procedures.

RESULTS

The mean Dice Similarity Coefficient using the proposed model was 0.79±0.08 (mean+SD) vs 0.72±0.06 using the geometrical model (P<.001).

CONCLUSION

The model provided an accurate iceball boundary prediction in less than 0.4second and has proven its feasibility to be implemented in an intraprocedural planning algorithm.

Calorimetric analysis of ice onset temperature during cryoablation: a model approach to identify early predictors of effective applications

Aim of the present study is to analyze thermal events occurring during cryoablation. Different bovine liver samples underwent freezing cycles at different cooling rate (from 0.0075 to 25 K/min). Ice onset temperature and specific latent heat capacity of the ice formation process were measured according to differential scanning calorimetry signals. A computational model of the thermal events occurring during cryoablation was compiled using Neumann’s analytical solution. Latent heat (#1 = 139.8 ± 7.4 J/g, #2 = 147.8 ± 7.9 J/g, #3 = 159.0 ± 4.1 J/g) of all liver samples was independent of the ice onset temperature, but linearly dependent on the water content. Ice onset temperature was proportional to the logarithm of the cooling rate in the range 5 ÷ 25 K/min (#3a = − 12.2 °C, #3b = − 16.2 °C, #3c = − 6.6 °C at 5K/min; #3a = − 16.5 °C, #3b = − 19.3 °C, #3c = − 11.6 °C at 25 K/min). Ice onset temperature was associated with both the way in which the heat involved into the phase transition was delivered and with the thermal gradient inside the tissue. Ice onset temperature should be evaluated in the early phase of the ablation to tailor cryoenergy delivery. In order to obtain low ice trigger temperatures and consequent low ablation temperatures a high cooling rate is necessary.

Thermodynamic properties of atrial fibrillation cryoablation: a model-based approach to improve knowledge on energy delivery

The objective of this study is to describe a suitable model of atrial fibrillation cryoablation thermodynamic properties. Three different thermal loads were applied to a cylindrical copper element simulating the cryoprobe, thermally coupled with a Peltier stack producing the freezing effect, and in contact with a bovine liver sample. Thermal events occurring inside the samples were measured using mirror image technique. Heat subtracted flux during ice formation and minimum temperature measured at probe-tissue interface were, respectively, 1.33 W cm^-2 and -27.8°C for Sample#0, 1.88 W cm^-2 and -35.6°C for Sample#1 and 1.82 W cm^-2 and 1.44 W cm^-2 before and after the ice trigger, respectively, and -29.3°C for Sample#2. Ice trigger temperature was around -8.5°C for Sample#0 and Sample#2, and -10.4°C for Sample#1. In all the investigated samples, ice front penetration was proportional to the square root of time and its velocity depended on the heat flux subtracted. The fraction of the useful energy spent for ice formation was less than 60% for Sample#0, and about 80% for Sample#1 and for Sample#2, before the reduction of the removed heat flux. Freezing time exceeding a cut-off, according to the heat subtracted flux, does not improve the procedure effectiveness and is detrimental to the surrounding tissues.

Computational study of the effects of arterial bifurcation on the temperature distribution during cryosurgery

BACKGROUND

Thermally significant blood flows into locally cooled diseased tissues and warm them during cryosurgery so that the iceball is often hard to cover the whole diseased volume. This paper is aimed at investigating the effects of large arterial bifurcation on the temperature distribution during cryosurgery through simulation method.

METHODS

A parametric geometry model is introduced to construct a close-to-real arterial bifurcation. The three-dimensional transient conjugate heat transfer between bifurcated artery and solid tissues with phase change during cryosurgery is performed by finite volume method.

RESULTS

The discussion was then made on the effects of the relative position between cryoprobe and artery bifurcation, the inlet velocity of root artery and the layout of multiple cryoprobes on the temperature distribution and iceball evolution. The results show that the thermal interaction between blood flow and iceball growth near bifurcation is considerable complex. The thermal effects of bifurcation could modulate the iceball morphology, severely weaken its freezing volume and prevent the blood vessel from being frozen.

CONCLUSION

The present work is expected to be valuable in optimizing cryosurgery scheme of the situation that the bifurcated artery is embedded into the disease tissue.

Mirror image technique for the thermal analysis in cryoablation: Experimental setup and validation

The paper presents a set of experiments that were performed to characterize the freezing front propagation in water first, and in an agar-gel solution afterwards. The experimental setup made of Peltier devices, to emulate the cryogenic effect, and a copper cold finger, to mimic the cold probe interface, are described. We claim that by monitoring some temperatures at the generating cryodevice, several pieces of information can be derived through the cold interface to assess the outside thermodynamic changes. The employed technique, known as mirror image, allows determining the occurrence of the initial ice formation outside the cryo-probe and in the surrounding material, also with different magnitudes of the thermal contact resistance at the cold interface. For both water and agar the ice penetration was found to be non linear versus time, and proportional to the square root of time in the performed experiments. The ice drift velocity decreases according to its penetration inside the tested materials. At the beginning of ice formation, the measured drift velocities are approximately 0.11 mm/s and 0.06 mm/s for water and agar, respectively, and after the ice penetrates 2 mm, they become approximately 0.03 mm/s for both materials.

Incorporating an immersed boundary method to study thermal effects of vascular systems during tissue cryo-freezing

In this paper, the three-dimensional thermal effects of a clinically-extracted vascular tissue undergoing cryo-freezing are numerically investigated. Based on the measured experimental temperature field, the numerical results of the Pennes bioheat model combined with the boundary condition-enforced immersed boundary method (IBM) agreed well with experimental data with a maximum temperature discrepancy of 2.9°C. For simulating the temperature profile of a tumor sited in a dominantly vascularized tissue, our model is able to capture with ease the thermal effects at specified junctions of the blood vessels. The vascular complexity and the ice-ball shape irregularity which cannot be easily quantified via clinical experiments are also analyzed and compared for both two-dimensional and three-dimensional settings with different vessel configurations and developments. For the three-dimensional numerical simulations, a n-furcated liver vessels model from a three-dimensional segmented volume using hole-making and subdivision methods is applied. A specific study revealed that the structure and complexity of the vascular network can markedly affect the tissue's freezing configuration with increasing ice-ball irregularity for greater blood vessel complexity.

Augmenting Surgery via Multi-scale Modeling and Translational Systems Biology in the Era of Precision Medicine: A Multidisciplinary Perspective

In this era of tremendous technological capabilities and increased focus on improving clinical outcomes, decreasing costs, and increasing precision, there is a need for a more quantitative approach to the field of surgery. Multiscale computational modeling has the potential to bridge the gap to the emerging paradigms of Precision Medicine and Translational Systems Biology, in which quantitative metrics and data guide patient care through improved stratification, diagnosis, and therapy. Achievements by multiple groups have demonstrated the potential for (1) multiscale computational modeling, at a biological level, of diseases treated with surgery and the surgical procedure process at the level of the individual and the population; along with (2) patient-specific, computationally-enabled surgical planning, delivery, and guidance and robotically-augmented manipulation. In this perspective article, we discuss these concepts, and cite emerging examples from the fields of trauma, wound healing, and cardiac surgery.

Computational Modeling for Enhancing Soft Tissue Image Guided Surgery: An Application in Neurosurgery

With the recent advances in computing, the opportunities to translate computational models to more integrated roles in patient treatment are expanding at an exciting rate. One area of considerable development has been directed towards correcting soft tissue deformation within image guided neurosurgery applications. This review captures the efforts that have been undertaken towards enhancing neuronavigation by the integration of soft tissue biomechanical models, imaging and sensing technologies, and algorithmic developments. In addition, the review speaks to the evolving role of modeling frameworks within surgery and concludes with some future directions beyond neurosurgical applications.

An efficient thermal evolution model for cryoablation with arbitrary multi-cryoprobe configuration

Cryoablation has been demonstrated powerful in treating of a variety of diseases, especially for the tumor ablation, which destroys the target tissue through the controlled freezing of cryoprobe. The prediction of temperature evolution during cryoablation is of great importance for developing and improving clinical procedure. This paper presented an efficient thermal model to characterize the freezing effect of cryoprobe with arbitrary layout including its size, orientation and number. The key step of the presented model is to establish a boundary heat source method to implicitly characterize the heat transfer from cryoprobe with fixed temperature or convective heat transfer boundary condition, which is furthermore incorporated to a fast parallel alternating direction explicit (PADE) finite difference method for computation acceleration. A novel dynamical and conformal computational region is designed through the shortest distance definition to balance the thermal effect of tissue and computational efficiency. The detailed test cases including a real head tissue demonstrated that the current model can accurately predict the temperature field evolution induced by arbitrary multi-cryoprobe configuration, and achieve significant computational ability due to allowable large time step (100-fold compared with the explicit finite difference method), compact computational region (at least reducing 40% number of voxels) and high parallel efficiency (speedup ratio about 8 for 12 threads) for complex tissue structure.

Image-guided tumor ablation: standardization of terminology and reporting criteria--a 10-year update

Image-guided tumor ablation has become a well-established hallmark of local cancer therapy. The breadth of options available in this growing field increases the need for standardization of terminology and reporting criteria to facilitate effective communication of ideas and appropriate comparison among treatments that use different technologies, such as chemical (eg, ethanol or acetic acid) ablation, thermal therapies (eg, radiofrequency, laser, microwave, focused ultrasound, and cryoablation) and newer ablative modalities such as irreversible electroporation. This updated consensus document provides a framework that will facilitate the clearest communication among investigators regarding ablative technologies. An appropriate vehicle is proposed for reporting the various aspects of image-guided ablation therapy including classification of therapies, procedure terms, descriptors of imaging guidance, and terminology for imaging and pathologic findings. Methods are addressed for standardizing reporting of technique, follow-up, complications, and clinical results. As noted in the original document from 2003, adherence to the recommendations will improve the precision of communications in this field, leading to more accurate comparison of technologies and results, and ultimately to improved patient outcomes.

Characterization of irreversible electroporation ablation in in vivo porcine liver

OBJECTIVE

The purpose of this study was to prospectively characterize and optimize irreversible electroporation ablation to determine the best parameters to achieve the largest target zones of coagulation for two electrodes.

MATERIALS AND METHODS

Ultrasound-guided irreversible electroporation ablation (n=110) was performed in vivo in 25 pig livers using two 18-gauge electroporation electrodes and an irreversible electroporation generator. Five variables for energy deposition and electrode configuration were sequentially studied: number of electrical pulses (n=20-90), length of pulses (20-100 microseconds), generator voltage (2250-3000 V), interelectrode spacing (1.5-2.5 cm), and length of active electrode exposure (1.0-3.0 cm). Zones of ablation were determined at gross pathology and histopathology 2-3 hours after irreversible electroporation. Dimensions were compared and subjected to statistical analysis.

RESULTS

For 1.5-cm spacing and 2-cm electrode exposure at 2250 V, there was no statistical difference in the size of coagulation when varying the number or length of pulses from 50 to 90 repetitions or 50-100 microseconds, respectively, with each parameter combination yielding 3.0±0.4×1.7±0.4×3.0±0.6 cm (width, depth, and height, respectively). Yet, increasing the pulse width or number over 70 caused increased hyperechogenic or gas and coagulation around the electrode. Increasing the voltage from 2250-3000 V for 70 pulses of 70 microseconds increased coagulation to 3.1±0.4×2.0±0.2 cm (p<0.01 for depth). Greater coagulation width of 3.9±0.5 cm (p<0.01) was achieved at 2-cm interelectrode spacing (with similar depth of 1.9±0.4 cm). However, consistent results required 90 repetitions and a 100-microsecond pulse width; 2.5-cm spacing resulted in two separate zones of ablation. Although electrode exposure did not influence width or depth, a linear correlation (r2=0.77) was noted for height, which ranged from 2.0±0.2-5.0±0.8 cm (for 1- and 3-cm exposures, respectively).

CONCLUSION

Predictable zones of tissue destruction can be achieved for irreversible electroporation. Ablation dimensions are sensitive to multiple parameters, suggesting that precise technique and attention to detail will be particularly important when using this modality.