A Combined Conduction-Radiation Model for analyzing the role of radiation on freezing of a Biological Tissue

[+] Author and Article Information
Siladitya Sukumar

School of Mechanical Engineering KIIT Deemed to be University Bhubaneswar, Odisha 751024 India siladitya.sukumar@rocketmail.com

Satya Prakash Kar

School of Mechanical Engineering, KIIT Bhubaneswar Bhubaneswar, Odisha 751024 India satyapkar@gmail.com

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the Journal of Thermal Science and Engineering Applications. Manuscript received April 4, 2019; final manuscript received July 24, 2019; published online xx xx, xxxx. Assoc. Editor: Sandip Mazumder.

ASME doi:10.1115/1.4044428 History: Received April 04, 2019; Accepted July 24, 2019


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.

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





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