Tissue freezing has significant applications in cryopreservation and cryosurgery process. The freezing rate is an important factor during cryopreservation process. To improve and make accurate estimation of the freezing rate, radiation heat transfer plays a major role. To analyze in details, a two-dimensional coupled conduction-radiation model is developed. The tissue is frozen from the left side while other sides are at the initial temperature. Finite volume method (FVM) is used to discretize both the Radiative transfer equation (RTE) and Energy equation. The algebraic equation after discretization is solved by the Tri-Diagonal Matrix Algorithm (TDMA). The radiative heat flux is calculated by solving the RTE. The Energy equation provides the temperature field. The Enthalpy-Porosity Method is used to update the liquid volume fraction and thus the freezing front is captured. A linearization technique is proposed to linearize the radiative source term in the Energy equation to avoid chances of divergence of the solution. The present model is first validated with the results of the existing literature and a good agreement is found. The effects of different parameters such as conduction-radiation parameter, scattering albedo, extinction coefficient and Stefan number on the prediction of temperature field and the position of the freezing front are studied in details. It is found that decrease in the values of conduction-radiation parameter, scattering albedo and increase in the values of extinction coefficient and Stefan number cause more radiative heat loss. Hence freezing is improved and more frozen region is observed.

This content is only available via PDF.

Article PDF first page preview

Article PDF first page preview
You do not currently have access to this content.