A study of heat transfer during cryosurgery of lung cancer

In this study, a mathematical model describing two-dimensional bio-heat transfer during cryosurgery of lung cancer is developed. The lung tissue is cooled by a cryoprobe by imposing its surface at a constant temperature or a constant heat flux or a constant heat transfer coefficient. The freezing starts and the domain is distributed into three stages namely: unfrozen, mushy and frozen regions. In stage I where the only unfrozen region is formed, our problem is an initial-boundary value problem of the hyperbolic partial differential equation. In stage II where mushy and unfrozen regions are formed, our problem is a moving boundary value problem of parabolic partial differential equations and in stage III where frozen, mushy, and unfrozen regions are formed, our problem is a moving boundary value problem of parabolic partial differential equations. The solution consists of the three-step procedure: (i) transformation of problem in non-dimensional form, (ii) by using finite differences, the problem converted into ordinary matrix differential equation and moving boundary problem of ordinary matrix differential equations, (iii) applying Legendre wavelet Galerkin method the problem is transferred into the generalized system of Sylvester equations which are solved by applying Bartels-Stewart algorithm of generalized inverse. The complete analysis is presented in the non-dimensional form. The consequence of the imposition of boundary conditions on moving layer thickness and temperature distribution are studied in detail. The consequence of Stefan number, Kirchoff number and Biot number on moving layer thickness are also studied in specific.