The bioheat transfer with phase change in biological tissues during the freezing process is simulated by the dual phase lag conduction heat transfer model. A numerical algorithm based on the enthalpy method is established to solve the solidification of biological tissues. The linearly temperature-dependent enthalpy (non-isothermal phase change) is considered here. The results of the parabolic heat conduction model for a slice of cucumber are compared with the experimental data. A comparison between dual phase lag and hyperbolic solutions with small values of relaxation times is applied in order to verify the corresponding parabolic solutions accuracy of the dual phase lag and hyperbolic solutions. The heating source effect owing to blood perfusion and metabolic heat on the heat transfer in a biological tissue subject to freezing process is studied. The relaxation time has an important influence on the transient temperature and temperature gradient. A major discrepancy among bioheat transfer models is found for zones closer to the cooling boundary. The heat source term, owing to blood flow and metabolism in a phase change problem in the biological tissue, has a significant influence on thermal effects of the subject tissue.
AkhtarTPeggDEForemanJ. The effect of cooling and warming rates on survival of cryopreserved L-cell. Cryobiology1979; 16: 224–229.
2.
FahyGM. Analysis of solution effect injury: cooling rate dependence of the functional and morphological sequelae of freezing in rabbit renal cortex protected with dimethy solsuxide. Cryobiology1981; 18: 550–570.
3.
OrpwoodRD. Biophysical and engineering aspect of cryosurgery. Phys Med Biol1981; 24(6): 555–575.
4.
OnikGRubinskyBZemelR. Ultrasound guided hepatic cryosurgery in the treatment of metastatic colon carcinoma. Cancer1991; 67(4): 901–907.
5.
WhittinghamDGLeiboSPMazurP. Survival of mouse embryos frozen to −196°C and −269°C. Science1972; 178: 411–414.
6.
LudwigMAl-HasaniSFelderbaumR. New aspects of cryopreservation of oocytes and embryos in assisted reproduction and future perspectives. Hum Reprod1999; 14: 162–185.
7.
AgcaY. Cryopreservation of oocyte and ovarian tissue. ILAR J2000; 4: 207–220.
8.
FowlerAJTonerM. Cryopreservation of cells using ultrarapid freezing. J Heat Transfer1997; 37: 179–183.
9.
KandraD. Tissue interactions with lasers and liquid nitrogen – a novel cryopreservation method. MSc Thesis, Department of Mechanical Engineering, Bangalore University, 2001.
10.
DengZSLiuJ. Modeling of multidimensional freezing problem during cryosurgery by the dual reciprocity boundary element method. Engng Analysis with Boundary Elements2004; 28(2): 97–108.
11.
RubinskyBShitzerA. Analysis of Stefan-like problem in a biological tissue around a cryosurgical probe. ASME J Heat Transfer1976; 98: 514–519.
12.
RabinYShitzerA. Exact solution for the inverse Stefan problem in non-ideal biological tissues. ASME J Heat Transfer1995; 117: 425–431.
13.
SaddMHDidlakeJE. Non-Fourier melting of a semi-infinite solid. ASME J Heat Transfer1977; 99: 25–28.
YooJRubinskyB. A finite element method for study of solidification processes in presence of natural convection. Int J Numerical Methods Engng1986; 23: 1785–1805.
16.
LazaridisA. A numerical solution for multidimensional solidification (or melting) problem. Int J Heat Mass Transfer1970; 13: 1459–1477.
17.
ShamsundarNSporrowEM. Analysis of multidimensional conduction phase change via the enthalpy method. ASME J Heat Transfer1975; 97(3): 333–340.
18.
SanzPDde ElviraC. Freezing rate simulation as an aid to reducing crystallization damage in foods. Meat Science1998; 52: 275–278.
19.
DengZSLiuJ. Non-Fourier heat conduction effect on prediction of temperature transients and thermal stress in skin cryopreservation. J Thermal Stresses2003; 26: 779–798.
20.
De SocioLMGualtieriG. A hyperbolic Stefan problem. Quart Appl Math1983; 41: 253–259.
21.
WangZWuHZhaoG. One-dimensional finite-difference modeling on temperature history and freezing time of individual food. J Food Engng2006; 79: 502–510.
22.
KumarSKatiyarVK. Mathematical modeling of freezing and thawing process in tissues: A porous media approach. Int J Appl Mech2010; 2(3): 617–633.
23.
DengZSLiuJ. Numerical simulation of 3-D freezing and heating problems for combined cryosurgery and hyperthermia therapy. Numerical Heat Transfer, Part A2004; 46: 587–611.
24.
ThakurHSinghKMSahooPK. Phase change problems using the MLPG method. Numerical Heat Transfer, Part A2011; 59: 438–458.
25.
DengZSLiuJ. Numerical study of the effects of large blood vessels on three-dimensional tissue temperature profile during cryosurgery. Numerical Heat Transfer, Part A2006; 49: 47–67.
26.
ZhaoGZhangHFGuoXJ. Effect of blood flow and metabolism on multidimensional heat transfer during cryosurgery. Medical Engng Physics2007; 29: 205–215.
27.
VollerVRSwaminathanCR. General source-based method for solidification phase change. Numerical Heat Transfer, Part B1991; 19: 175–189.
28.
DengZSLiuJ. Analytical study on bioheat transfer problems with spatial or transient heating on skin surface or inside biological bodies. ASME J Biomech Eng2002; 124: 638–649.
29.
ThomasBGSamarasekeraIVBrimacombeJK. Comparison of numerical modeling techniques for complex, two dimensional, transient heat-conduction problems. Metall Trans B1984; 15: 307–318.
30.
ClyneTW. Modeling of heat flow in solidification. Mater Sci Engng1984; 65: 111–124.
31.
BasuBSekharJA. Modeling of multidimensional solidification of an alloy. Metall Trans A1989; 20: 1833–1845.
32.
DroouxJJ. Three-dimensional numerical simulation of solidification by an improved explicit scheme. Comp Methods Appl Mech Engng1991; 85: 57–74.
33.
TannehillJCAndersonDAPletcherRH. Computational fluid mechanics and heat transfer.2nd ed.Washington DC: Taylor & Francis, 1997.
34.
KunadianIMc DonoughJMTagaviKA. Numerical simulation of heat transfer mechanism during femtosecond laser heating of nano-film using 3D dual phase lag model. In: Proceedings of the HTFED04 2004 ASME heat transfer/fluids engineering summer conference, Charlotte, North Carolina, USA, 2004; 4: 661–667.
35.
WesslingFCBlackshearPL. The thermal properties of human blood during freezing process. ASME J Heat Transfer1973; 95: 246–249.
36.
DengZSLiuJ. Numerical simulation of selective freezing of target biological tissues following injection of solutions with specific thermal properties. Cryobiology2005; 50: 183–192.
37.
TuJXLiuBL. Numerical simulation and experimental study on the freezing process of cucumber. J University of Shanghai Sci Technol2000, 22(4): 304–307.
