Current research in cryosurgery is concerned with finding a thermal history that will definitively destroy tissue. In this study, we measured and predicted the thermal history obtained during freezing and thawing in a cryosurgical model. This thermal history was then compared to the injury observed in the tissue of the same cryosurgical model (reported in companion paper (Hoffmann and Bischof, 2001)). The dorsal skin flap chamber, implanted in the Copenhagen rat, was chosen as the cryosurgical model. Cryosurgery was performed in the chamber on either normal skin or tumor tissue propagated from an AT-1 Dunning rat prostate tumor. The freezing was performed by placing a ∼1 mm diameter liquid-nitrogen-cooled cryoprobe in the center of the chamber and activating it for approximately 1 minute, followed by a passive thaw. This created a 4.2 mm radius iceball. Thermocouples were placed in the tissue around the probe at three locations (r=2, 3, and 3.8 mm from the center of the window) in order to monitor the thermal history produced in the tissue. The conduction error introduced by the presence of the thermocouples was investigated using an in vitro simulation of the in vivo case and found to be <10°C for all cases. The corrected temperature measurements were used to investigate the validity of two models of freezing behavior within the iceball. The first model used to approximate the freezing and thawing behavior within the DSFC was a two-dimensional transient axisymmetric numerical solution using an enthalpy method and incorporating heating due to blood flow. The second model was a one-dimensional radial steady state analytical solution without blood flow. The models used constant thermal properties for the unfrozen region, and temperature-dependent thermal properties for the frozen region. The two-dimensional transient model presented here is one of the first attempts to model both the freezing and thawing of cryosurgery. The ability of the model to calculate freezing appeared to be superior to the ability to calculate thawing. After demonstrating that the two-dimensional model sufficiently captured the freezing and thawing parameters recorded by the thermocouples, it was used to estimate the thermal history throughout the iceball. This model was used as a basis to compare thermal history to injury assessment (reported in companion paper (Hoffmann and Bischof, 2001)).

1.
Hoffmann
,
N. E.
, and
Bischof
,
J. C.
,
2001
, “
Cryosurgery of Normal and Tumor Tissue in the Dorsal Skin Flap Chamber: Part II — Injury Response
,”
ASME J. Biomech. Eng.
,
123
, this issue, pp.
310
316
.
2.
Cooper
,
T. E.
, and
Petrovic
,
W. K.
,
1974
, “
An Experimental Investigation of the Temperature Field Produced by a Cryosurgical Cannula
,”
ASME J. Heat Transfer
,
96
, pp.
415
420
.
3.
Rewcastle
,
J. C.
,
Sandison
,
G. A.
,
Hahn
,
L. J.
,
Saliken
,
J. C.
,
McKinnon
,
J. G.
, and
Donnelly
,
B. J.
,
1998
, “
A Model for the Time-Dependent Thermal Distribution Within an Iceball Surrounding a Cryoprobe
,”
Phys. Med. Biol.
,
43
, No.
12
, pp.
3519
3534
.
4.
Budman
,
H.
,
Shitzer
,
A.
, and
Del Guidice
,
S.
,
1986
, “
Investigation of Temperature Fields Around Embedded Cryoprobes
,”
ASME J. Biomech. Eng.
,
108
, pp.
42
48
.
5.
Hong
,
J.-S.
,
Wong
,
S.
,
Pease
,
G.
, and
Rubinsky
,
B.
,
1994
, “
MR Imaging Assisted Temperature Calculations During Cryosurgery
,”
Magn. Reson. Imaging
,
12
, No.
7
, pp.
1021
1031
.
6.
Bischof
,
J. C.
,
Bastacky
,
J.
, and
Rubinsky
,
B.
,
1992
, “
An Analytical Study of Cryosurgery in the Lung
,”
ASME J. Biomech. Eng.
,
114
, pp.
467
472
.
7.
Gori, F., 1981, “On Heat Transfer in Non-Living Materials Around a Cryosurgical Probe,” in: Numerical Methods in Thermal Problems, Pineridge Press, Venice, Italy.
8.
Rubinsky
,
B.
, and
Shitzer
,
A.
,
1976
, “
Analysis of a Stefan-Like Problem in a Biological Tissue Around a Cryosurgical Probe
,”
ASME J. Heat Transfer
,
98
, pp.
514
519
.
9.
Smith, D. J., Devireddy, R. V., and Bischof, J. C., 1999, “Prediction of Thermal History and Interface Propagation During Freezing in Biological Systems—Latent Heat and Temperature-Dependent Property Effects,” in: Proc. 54th ASME/JSME Joint Thermal Engineering Conference, San Diego, CA, AJTE 99-6520.
10.
Rabin
,
Y.
, and
Shitzer
,
A.
,
1997
, “
Combined Solution of the Inverse Stefan Problem for Successive Freezing/Thawing in Nonideal Biological Tissues
,”
ASME J. Biomech. Eng.
,
119
, pp.
146
152
.
11.
Rabin
,
Y.
, and
Shitzer
,
A.
,
1998
, “
Numerical Solution of the Multidimensional Freezing Problem During Cryosurgery
,”
ASME J. Biomech. Eng.
,
120
, pp.
32
37
.
12.
Gage
,
A. A.
,
Caruana
,
J. A.
, and
Garamy
,
G.
,
1983
, “
A Comparison of Instrument Methods of Monitoring Freezing in Cryosurgery
,”
J. Dermatol. Surg. Oncol.
,
9
, No.
3
, pp.
209
214
.
13.
Rabin
,
Y.
,
1998
, “
Uncertainty in Temperature Measurements Using Thermocouples
,”
Cryo-Letters
,
19
, No.
4
, pp.
213
224
.
14.
Papenfuss
,
H. D.
,
Gross
,
J. F.
,
Intaglietta
,
M.
, and
Treese
,
F. A.
,
1979
, “
A Transparent Access Chamber for the Rat Dorsal Skin Fold
,”
Microvasc. Res.
,
18
, No.
3
, pp.
311
318
.
15.
Gaber
,
M. H.
,
Wu
,
N. Z.
,
Hong
,
K.
,
Huang
,
S. K.
,
Dewhirst
,
M. W.
, and
Papahadjopoulos
,
D.
,
1996
, “
Thermosensitive Liposomes: Extravasation and Release of Contents in Tumor Microvascular Networks
,”
Int. J. Radiat. Oncol., Biol., Phys.
,
36
, No.
5
, pp.
1177
1187
.
16.
Aggarwal
,
S. J.
,
Shah
,
S. J.
,
Diller
,
K. R.
, and
Baxter
,
C. R.
,
1989
, “
Fluorescence Digital Microscopy of Interstitial Macromelecular Diffusion in Burn Injury
,”
Comput. Biol. Med.
,
19
, No.
4
, pp.
245
261
.
17.
Yuan
,
F.
,
Leunig
,
M.
,
Berk
,
D. A.
, and
Jain
,
R. K.
,
1993
, “
Microvascular Permeability of Albumin, Vascular Surface Area, and Vascular Volume Measured in Human Adenocarcinoma LS174T Using Dorsal Chamber in SCID Mice
,”
Microvasc. Res.
,
45
, No.
3
, pp.
269
89
.
18.
Endrich
,
B.
,
Laprell-Moschner
,
C.
,
Brendel
,
W.
, and
Messmer
,
K.
,
1982
, “
Effects of Prolonged Cold Injury on the Subcutaneous Microcirculation of the Hamster. I. Technique, Morphology and Tissue Oxygenation
,”
Res. Exp. Med.
,
181
, No.
1
, pp.
49
61
.
19.
Diller
,
K. R.
, and
Hayes
,
L. J.
,
1991
, “
Analysis of Tissue Injury by Burning: Comparison of in Situ and Skin Flap Models
,”
Int. J. Heat Mass Transf.
,
34
, No.
6
, pp.
1393
1406
.
20.
Gross
,
J. F.
,
Roemer
,
R.
,
Dewhirst
,
M.
, and
Meyer
,
T.
,
1982
, “
A Uniform Thermal Field in a Hyperthermia Chamber for Microvascular Studies
,”
Int. J. Heat Mass Transf.
,
25
, pp.
1313
1320
.
21.
Bischof
,
J. C.
,
Smith
,
D. J.
,
Pazhayannur
,
P. V.
,
Manivel
,
C.
,
Hulbert
,
J.
, and
Roberts
,
K. P.
,
1997
, “
Cryosurgery of Dunning AT-1 Rat Prostate Tumor: Thermal, Biophysical and Viability Response at the Cellular and Tissue Level
,”
Cryobiology
,
34
, pp.
42
69
.
22.
Smith
,
D. J.
,
Fahssi
,
W. M.
,
Swanlund
,
D. J.
, and
Bischof
,
J. C.
,
1999
, “
A Parametric Study of Freezing Injury in AT-1 Rat Prostate Tumor Cells
,”
Cryobiology
,
39
, pp.
13
28
.
23.
Fennema, O. R., Powrie, W. D., and Marth, E. H., 1973, Low Temperature Preservation of Food and Living Matter, Marcel Dekker, New York.
24.
Devireddy
,
R. V.
,
Smith
,
D. J.
, and
Bischof
,
J. C.
,
1999
, “
Mass Transfer During Freezing in Rat Prostate Tumor Tissue
,”
AIChE J.
,
45
, No.
3
, pp.
639
654
.
25.
Reinoso
,
R. F.
,
Telfer
,
B. A.
, and
Rowland
,
M.
,
1997
, “
Tissue Water Content in Rats Measured by Desiccation
,”
J. Pharmacol. Toxicol. Methods
,
38
, No.
2
, pp.
87
92
.
26.
Alexiades, A., and Solomon, A. D., 1993, “Formulation of the Stefan Problem,” in: Mathematical Modeling of Melting and Freezing Processes, Hemisphere Publishing, Washington, pp. 6–26.
27.
Guide to Plastics, 1987, Modern Plastics Publication, McGraw-Hill New York.
28.
Lunardini, V., 1981, “Finite Difference Methods for Freezing and Thawing,” in: Heat Transfer in Cold Climates, Van Nostrand Reinhold, Co., New York.
29.
Voller
,
V. R.
, and
Swaminathan
,
C. R.
,
1993
, “
Treatment of Discontinuous Thermal Conductivity in Control-Volume Solutions of Phase Change Problems
,”
Numer. Heat Transfer
,
24
, pp.
161
180
.
30.
Burton
,
M. A.
,
Kelleher
,
D. K.
,
Gray
,
B. N.
, and
Morgan
,
C. K.
,
1990
, “
Effect of Temperature on Liver Tumour Blood Flow
,”
Eur. J. Cancer
,
26
, p.
999
999
.
31.
Cooper
,
T. E.
, and
Trezek
,
G. J.
,
1971
, “
Rate of Lesion Growth Around Spherical and Cylindrical Cryoprobes
,”
Cryobiology
,
7
, pp.
183
190
.
32.
Hayes
,
L. J.
,
Diller
,
K. R.
,
Chang
,
H. J.
, and
Lee
,
H. S.
,
1988
, “
Prediction of Local Cooling Rates and Cell Survival During the Freezing of a Cylindrical Specimen
,”
Cryobiology
,
25
, No.
1
, pp.
67
82
.
33.
Rabin
,
Y.
,
2000
, “
The Effect of Temperature-Dependent Thermal Conductivity in Heat Transfer Simulations of Frozen Biomaterials
,”
Cryo-Letters
,
21
, No.
3
, pp.
163
170
.
34.
Bevington, P. R., and Robinson, D. K., 1992, “Error Analysis,” in: Data Reduction and Error Analysis for the Physical Sciences, 2nd ed., McGraw-Hill, New York, pp. 38–48.
35.
Hayes
,
L. J.
, and
Diller
,
K. R.
,
1983
, “
Implementation of Phase Change in Numerical Models of Heat Transfer
,”
J. Energy Resour. Technol.
,
105
, pp.
431
435
.
You do not currently have access to this content.