Research Papers

Methods for Characterizing Convective Cryoprobe Heat Transfer in Ultrasound Gel Phantoms

[+] Author and Article Information
Michael L. Etheridge

Department of Mechanical Engineering,
Department of Biomedical Engineering,
University of Minnesota,
Minneapolis, MN 55455

Jeunghwan Choi

Department of Mechanical Engineering,
University of Minnesota,
Minneapolis, MN 55455

Satish Ramadhyani

Galil Medical Inc.,
Arden Hills, MN 55112

John C. Bischof

Department of Mechanical Engineering,
Department of Biomedical Engineering,
Department of Urologic Surgery,
University of Minnesota,
Minneapolis, MN 55455

Contributed by the Bioengineering Division of ASME for publication in the JOURNAL OF BIOMECHANICAL ENGINEERING. Manuscript received September 4, 2012; final manuscript received December 6, 2012; accepted manuscript posted December 22, 2012; published online February 7, 2013. Assoc. Editor: Michael Sacks.

J Biomech Eng 135(2), 021002 (Feb 07, 2013) (10 pages) Paper No: BIO-12-1391; doi: 10.1115/1.4023237 History: Received September 04, 2012; Revised December 06, 2012

While cryosurgery has proven capable in treating of a variety of conditions, it has met with some resistance among physicians, in part due to shortcomings in the ability to predict treatment outcomes. Here we attempt to address several key issues related to predictive modeling by demonstrating methods for accurately characterizing heat transfer from cryoprobes, report temperature dependent thermal properties for ultrasound gel (a convenient tissue phantom) down to cryogenic temperatures, and demonstrate the ability of convective exchange heat transfer boundary conditions to accurately describe freezing in the case of single and multiple interacting cryoprobe(s). Temperature dependent changes in the specific heat and thermal conductivity for ultrasound gel are reported down to −150 °C for the first time here and these data were used to accurately describe freezing in ultrasound gel in subsequent modeling. Freezing around a single and two interacting cryoprobe(s) was characterized in the ultrasound gel phantom by mapping the temperature in and around the “iceball” with carefully placed thermocouple arrays. These experimental data were fit with finite-element modeling in COMSOL Multiphysics, which was used to investigate the sensitivity and effectiveness of convective boundary conditions in describing heat transfer from the cryoprobes. Heat transfer at the probe tip was described in terms of a convective coefficient and the cryogen temperature. While model accuracy depended strongly on spatial (i.e., along the exchange surface) variation in the convective coefficient, it was much less sensitive to spatial and transient variations in the cryogen temperature parameter. The optimized fit, convective exchange conditions for the single-probe case also provided close agreement with the experimental data for the case of two interacting cryoprobes, suggesting that this basic characterization and modeling approach can be extended to accurately describe more complicated, multiprobe freezing geometries. Accurately characterizing cryoprobe behavior in phantoms requires detailed knowledge of the freezing medium's properties throughout the range of expected temperatures and an appropriate description of the heat transfer across the probe's exchange surfaces. Here we demonstrate that convective exchange boundary conditions provide an accurate and versatile description of heat transfer from cryoprobes, offering potential advantages over the traditional constant surface heat flux and constant surface temperature descriptions. In addition, although this study was conducted on Joule–Thomson type cryoprobes, the general methodologies should extend to any probe that is based on convective exchange with a cryogenic fluid.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Onik, G. M., Cohen, J. K., Reyes, G. D., Rubinsky, B., Chang, Z., and BaustJ., 1993, “Transrectal Ultrasound-Guided Percutaneous Radical Cryosurgical Ablation of the Prostate,” Cancer, 72(4), pp. 1291–1299. [CrossRef] [PubMed]
Uchida, M., Imaide, Y., Sugimoto, K., Uehara, H., and WatanabeH., 1995, “Percutaneous Cryosurgery for Renal Tumours,” Brit. J. Urology, 75(2), pp. 132–137. [CrossRef]
Ravikumar, T. S., Kane, R., Cady, B., Jenkins, R., Clouse, M., and Steele, G.Jr., 1991, “A 5-Year Study of Cryosurgery in the Treatment of Liver Tumors,” Arch. Surgery126(12), pp. 1520–1524. [CrossRef]
Maiwand, M. O., and AsimakopoulosG., 2004, “Cryosurgery for Lung Cancer: Clinical Results and Technical Aspects,” Tech. Cancer Res. Treat., 3(2), pp. 143–150.
Cox, J. L., Ferguson, T. B.Jr., Lindsay, B. D., and CainM. E., 1990, “Perinodal Cryosurgery for Atrioventricular Node Reentry Tachycardia in 23 Patients,” J. Thorac. Cardiov. Sur., 99(3), pp. 440–449.
Avitall, B., Urboniene, D., Rozmus, G., Lafontaine, D., Helms, R., and UrbonasA., 2003, “New Cryotechnology for Electrical Isolation of the Oulmonary Veins,” J. Cardiovac. Electr., 14(3), pp. 281–286. [CrossRef]
Barnard, D., Lloyd, J., and EvansJ., 1981, “Cryoanalgesia in the Management of Chronic Facial Pain,” J. Maxillofac. Surg., 9, pp. 101–102. [CrossRef] [PubMed]
Allen, B. H., Fallat, L. M., and SchwartzS. M., 2007, “Cryosurgery: An Innovative Technique for the Treatment of Plantar Fasciitis,” J. Foot Ankle Surg., 46(2), pp. 75–79. [CrossRef] [PubMed]
RubinskyB., 2000, “Cryosurgery,” Ann. Rev. Biomed. Eng., 2(1), pp. 157–187. [CrossRef]
Gage, A. A., and BaustJ., 1998, “Mechanisms of Tissue Injury in Cryosurgery,” Cryobiology, 37(3), pp. 171–186. [CrossRef] [PubMed]
Hoffmann, N. E., and BischofJ. C., 2002, “The Cryobiology of Cryosurgical Injury,” Urology, 60(2), pp. 40–49. [CrossRef] [PubMed]
Baust, J. G., Gage, A. A., Klossner, D., Clarke, D., Miller, R., Cohen, J., Katz, A., Polascik, T., Clarke, H., and BaustJ. M., 2007, “Issues Critical to the Successful Application of Cryosurgical Ablation of the Prostate.,” Tech. Cancer Res. Treat., 6(2), pp. 97–109.
Jiang, J., Goel, R., Schmechel, S., Vercellotti, G., Forster, C., and BischofJ., 2010, “Pre-Conditioning Cryosurgery: Cellular and Molecular Mechanisms and Dynamics of TNF-Alpha Enhanced Cryotherapy in an In Vivo Prostate Cancer Model System,” Cryobiology, 61(3), pp. 280–288. [CrossRef] [PubMed]
Cooper, T. E., and TrezekG. J., 1970, “Analytical Prediction of the Temperature Field Emanating From a Cryogenic Surgical Cannula,” Cryobiology, 7(2), pp. 79–83. [CrossRef] [PubMed]
Cooper, T. E., and TrezekG. J., 1971, “Rate of Lesion Growth Around Spherical and Cylindrical Cryoprobes,” Cryobiology, 7(4), pp. 183–190. [CrossRef]
Bischof, J. C., Bastacky, J., and RubinskyB., 1992, “An Analytical Study of Cryosurgery in the Lung,” J. Biomech. Eng., 114, pp. 467–472. [CrossRef] [PubMed]
Schweikert, R. J., and KeaniniR. G., 1999, “A Finite Element and Order of Magnitude Analysis of Cryosurgery in the Lung,” Int. Commun. Heat Mass Transf., 26(1), pp. 1–12. [CrossRef]
Wan, R., Liu, Z., Muldrew, K., and RewcastleJ., 2003, “A Finite Element Model for Ice Ball Evolution in a Multi-probe Cryosurgery,” Comp. Meth. Biomech. Biomed. Eng., 6(3), pp. 197–208. [CrossRef]
RabinY., 2000, “The Effect of Temperature-Dependent Thermal Conductivity in Heat Transfer Simulations of Frozen Biomaterials,” Cryo Lett., 21(3), pp. 163–170.
AusG., 2006, “Current Status of HIFU and Cryotherapy in Prostate Cancer—A Review,” Eur. Urology, 50(5), pp. 927–934. [CrossRef]
Baissalov, R., Sandison, G., Donnelly, B., Saliken, J., McKinnon, J., Muldrew, K., and RewcastleJ., 2000, “A Semi-Empirical Treatment Planning Model for Optimization of Multiprobe Cryosurgery,” Phys. Med. Biol., 45, pp. 1085–1098. [CrossRef] [PubMed]
Rewcastle, J. C., Sandison, G. A., Hahn, L. J., Saliken, J. C., McKinnon, J. G., and DonnellyB. J., 1998, “A Model for the Time-Dependent Thermal Distribution Within an Iceball Surrounding a Cryoprobe,” Phys. Med. Biol., 43, pp. 3519–3534. [CrossRef] [PubMed]
Rewcastle, J. C., Sandison, G. A., Muldrew, K., Saliken, J. C., and DonnellyB. J., 2001, “A Model for the Time Dependent Three-Dimensional Thermal Distribution Within Iceballs Surrounding Multiple Cryoprobes,” Med. Phys., 28, pp. 1125–1137. [CrossRef] [PubMed]
Rossi, M. R., and RabinY., 2007, “Experimental Verification of Numerical Simulations of Cryosurgery With Application to Computerized Planning,” Phys. Med. Biol., 52, pp. 4553–4567. [CrossRef] [PubMed]
Rossi, M. R., Tanaka, D., Shimada, K., and RabinY., 2007, “An Efficient Numerical Technique for Bioheat Simulations and Its Application to Computerized Cryosurgery Planning,” Comp. Meth. Prog. Biomed., 85(1), pp. 41–50. [CrossRef]
Magalov, Z., Shitzer, A., and DeganiD., 2007, “Isothermal Volume Contours Generated in a Freezing Gel by Embedded Cryo-Needles With Applications to Cryo-Surgery,” Cryobiology, 55(2), pp. 127–137. [CrossRef] [PubMed]
Kim, C., O'Rourke, A. P., Mahvi, D. M., and WebsterJ. G., 2007, “Finite-Element Analysis of Ex Vivo and In Vivo Hepatic Cryoablation,” Biomed. Eng. IEEE Trans., 54(7), pp. 1177–1185. [CrossRef]
Blezek, D. J., Carlson, D. G., Cheng, L. T., Christensen, J. A., Callstrom, M. R., and EricksonB. J., 2010, “Cell Accelerated Cryoablation Simulation,” Comp. Meth. Prog. Biodmed., 98(3), pp. 241–252. [CrossRef]
Keanini, R. G., and RubinskyB., 1992, “Optimization of Multiprobe Cryosurgery,” ASME J. Heat Trans., 114(4), pp. 796–801. [CrossRef]
He, X., and Bischof, J. C., 2005, “Analysis of Thermal Stress in Cryosurgery of Kidneys,” J. Biomech. Eng., 127, pp. 656–661. [CrossRef] [PubMed]
Zhang, J., Sandison, G. A., Murthy, J. Y., and XuL. X., 2005, “Numerical Simulation for Heat Transfer in Prostate Cancer Cryosurgery,” J. Biomech. Eng., 127, pp. 279–294. [CrossRef] [PubMed]
ZhangX., 2008, “Multi-Probe Cryosurgical Treatments: An Experimental and Computational Study Confirms Probe Cooling Power is Probe Number and Placement Dependent,” M.S., University of Minnesota, Minneapolis, MN.
Young, J. L., Kolla, S. B., Pick, D. L., Sountoulides, P., Kaufmann, O. G., Ortiz-Vanderdys, C. G., Huynh, V. B., Kaplan, A. G., Andrade, L. A., Osann, K. E., Louie, M. K., McDougall, E. M., and Clayman, R. V., 2010, “In Vitro, Ex Vivo and In Vivo Isotherms for Renal Cryotherapy,” J. Urology, 183(2), pp. 752–758. [CrossRef]
Choi, J., and BischofJ. C., 2010, “Review of Biomaterial Thermal Property Measurements in the Cryogenic Regime and Their Use for Prediction of Equilibrium and Non-Equilibrium Freezing Applications in Cryobiology,” Cryobiology, 60(1), pp. 52–70. [CrossRef] [PubMed]
Choi, J. H., and BischofJ. C., 2008, “A Quantitative Analysis of the Thermal Properties of Porcine Liver With Glycerol at Subzero and Cryogenic Temperatures,” Cryobiology, 57(2), pp. 79–83. [CrossRef] [PubMed]
Patel, P. A., Valvano, J. W., Pearce, J. A., Prahl, S. A., and DenhamC. R., 1987, “A Self-Heated Thermistor Technique to Measure Effective Thermal Properties From the Tissue Surface,” J. Biomech. Eng., 109, pp. 330–335. [CrossRef] [PubMed]
Zhang, M., Che, Z., Chen, J., Zhao, H., Yang, L., Zhong, Z., and LuJ., 2011, “Experimental Determination of Thermal Conductivity of Water- Agar Gel at Different Concentrations and Temperatures,” J. Chem. Eng. Data, 56(4), pp. 859–864. [CrossRef]
BergmanT. L., LavineA. S., IncroperaF. P., and DeWittD. P., 2011, Fundamentals of Heat and Mass Transfer, John Wiley & Sons, Hoboken, NJ.
ReedR. P., and HoriuchiT., 1983, Austenitic Steels at Low Temperatures, Plenum Publishing Corp., New York.
Jensen, J. E., Tuttle, W. A., Stewart, R. B., Brechna, H., and ProdellA. G., 1980, Brookhaven National Laboratory Selected Cryogenic Data Notebook: Volume 1, Sections 1–9, Brookhaven National Lab., Upton, NY.
RabinovichV. A., and SeloverT. B., 1988, Thermophysical Properties of Neon, Argon, Krypton, and Xenon, Hemisphere PubCorp, Berlin, New York, Washington, DC.
HarveyA. H., 2012, “Properties of Ice and Supercooled Water,” CRC Handbook of Chemistry and Physics 2012–2013, CRC Press, Boca Raton, FL, pp. 6–12.
LiC., 1976, “Thermal Conductivity of Liquid Mixtures,” AIChE J., 22(5), pp. 927–930. [CrossRef]
Angell, C. A., and SmithD. L., 1982, “Test of the Entropy Basis of the Vogel-Tammann-Fulcher Equation. Dielectric Relaxation of Polyalcohols Near Tg,” J. Phys. Chem., 86(19), pp. 3845–3852. [CrossRef]
Sun, T., and TejaA. S., 2003, “Density, Viscosity, and Thermal Conductivity of Aqueous Ethylene, Diethylene, and Triethylene Glycol Mixtures Between 290 K and 450 K,” J. Chem. Eng. Data, 48(1), pp. 198–202. [CrossRef]
Silverman, S. G., Tuncali, K., Adams, D. F., vanSonnenberg, E., Zou, K. H., Kacher, D. F., Morrison, P. R., and JoleszF. A., 2000, “MR Imaging-Guided Percutaneous Cryotherapy of Liver Tumors: Initial Experience,” Radiology, 217(3), pp. 657–664. [PubMed]
Wansapura, J. P., Daniel, B. L., Vigen, K. K., and ButtsK., 2005, “In Vivo MR Thermometry of Frozen Tissue Using R2* and Signal Intensity,” Acad. Radiol., 12(9), pp. 1080–1084. [CrossRef] [PubMed]
Rossi, M. R., Tanaka, D., Shimada, K., and RabinY., 2008, “Computerized Planning of Cryosurgery Using Bubble Packing: An Experimental Validation on a Phantom Material,” Int. J. Heat Mass Transf., 51(23), pp. 5671–5678. [CrossRef] [PubMed]
MartinH., 1977, “Heat and Mass Transfer Between Impinging Gas Jets and Solid Surfaces,” Adv. Heat Transf., 13, pp. 1–60. [CrossRef]


Grahic Jump Location
Fig. 3

Model setup for 2D, axisymmetric, single-probe case (a). The temperature (T) distribution around the freezing cryosurgical probe is numerically solved using the described boundary conditions and heat transfer equations. The dual-probe case is modeled using the same heat transfer conditions, but utilizes a 3D geometry with symmetry planes ((b), symmetry highlighted in green on near and left-hand faces).

Grahic Jump Location
Fig. 2

Experimental setup. An array of thermocouples was positioned around the cryoprobe(s) during freezing at different radial positions (r) and on different measurement planes (z). The temperature distribution was measured for both single- and dual-probe cases using the same placement jig.

Grahic Jump Location
Fig. 1

Simplified internal exchange flow for the 1.5 mm diameter IceSeed (Galil Medical Inc., Arden Hills, MN) cryoablation needle studied here. Compressed argon gas flows down through a center channel and impinges on the top of a stainless steel tip, where significant heat exchange occurs. Convective exchange continues up the walls of the probe for 12 mm, at which point the cryoprobe is insulated for the remainder of the shaft.

Grahic Jump Location
Fig. 4

Measured ultrasound gel specific heat (a) and thermal conductivity (b). Reference values for water and ice were included for comparison [42].

Grahic Jump Location
Fig. 5

Measured temperature profiles at the z = 4 mm plane and various radial distances for the single- and dual-probe cases (b–d). The single-probe case is plotted as the black curve for each radial distance. Colored curves correspond to the thermocouple grid locations (a), as first described in Fig. 2.

Grahic Jump Location
Fig. 6

Fitted values for the convective heat transfer boundary conditions for case 4 (see Table 2)

Grahic Jump Location
Fig. 7

Comparison of single- and dual-probe model results. Note the difference in isotherm shapes between the left and right sides of the probe. Isotherms on the left side are strongly influenced by the presence of the second probe and interactions between the probes results in expanded isotherms on the right side as well.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In