A model for the time-dependent thermal distribution within an iceball surrounding a cryoprobe

The optimal cooling parameters to maximize cell necrosis in different types of tissue have yet to be determined. However, a critical isotherm is commonly adopted by cryosurgeons as a boundary of lethality for tissue. Locating this isotherm within an iceball is problematic due to the limitations of MRI, ultrasound and CT imaging modalities. This paper describes a time-dependent two-dimensional axisymmetric model of iceball formation about a single cryoprobe and extensively compares it with experimental data. Thermal histories for several points around a CRYOprobe(TM) are predicted to high accuracy (5 degrees C maximum discrepancy). A realistic three-dimensional probe geometry is specified and cryoprobe temperature may be arbitrarily set as a function of time in the model. Three-dimensional temperature distributions within the iceball, predicted by the model at different times, are presented. Isotherm locations, as calculated with the infinite cylinder approximation, are compared with those of the model in the most appropriate region of the iceball. Infinite cylinder approximations are shown to be inaccurate when applied to this commercial probe. Adopting the infinite cylinder approximation to locate the critical isotherm is shown to lead the user to an overestimate of the volume of target tissue enclosed by this isotherm which may lead to incomplete tumour ablation.