Analysis of thermal stress in cryosurgery of kidneys

In this study, the thermal stress distribution in cryosurgery of kidney was investigated using a multi-physics finite element model developed in ANSYS (V7.0). The thermal portion of the model was verified using experimental data and the mechanics portion of the model (elastic) was verified using classic analytical solutions. Temperature dependent thermal and mechanical properties were used in the model. Moreover, the model accounts for thermal expansion due to both temperature change and volumetric expansion associated with phase change of tissue water to ice. For a clinical cylindrical cryoprobe inserted into the renal cortex from the top-middle renal capsule, it was found that the thermal stress distributions along the radial position are very different at different depths from the top renal capsule. The thermal stress is much higher at both ends than in the middle of the cryoprobe surface. It was found that there might be more tissue next to the top renal capsule which undergoes microcrack formation or plastic deformation. Furthermore, it was found that macrocrack formation is more likely to occur in tissue adjacent to the cryoprobe surface (especially on the sharp point tip) and during the thawing phase of cryosurgery. Because the thermal stress adjacent to the cryoprobe is much higher than the yield stress of frozen renal tissue, a plastic stress model is required for better modeling of the thermal stress distribution in cryosurgery of kidney in future although the computational effort will be drastically increased due to the strong nonlinear nature of the plastic model and more studies are required to understand the mechanical properties of frozen tissue.