Nonclassical Heat Transfer Models for Laser-Induced Thermal Damage in Biological Tissues

Author(s):  
Jianhua Zhou ◽  
J. K. Chen ◽  
Yuwen Zhang

To ensure personal safety and improve treatment efficiency in laser medical applications, one of the most important issues is to understand and accurately assess laser-induced thermal damage to biological tissues. Biological tissues generally consist of nonhomogeneous inner structures, in which heat flux equilibrates to the imposed temperature gradient via a thermal relaxation mechanism which cannot be explained by the traditional parabolic heat conduction model based on Fourier’s law. In this article, two non-Fourier heat conduction models, hyperbolic thermal wave model and dual-phase-lag (DPL) model, are formulated to describe the heat transfer in living biological tissues with blood perfusion and metabolic heat generation. It is shown that the non-Fourier bioheat conduction models could predict significantly different temperature and thermal damage in tissues from the traditional parabolic model. It is also found that the DPL bioheat conduction equations can be reduced to the Fourier heat conduction equations only if both phase lag times of the temperature gradient (τT) and the heat flux (τq) are zero. Effects of laser parameters and blood perfusion on the thermal damage simulated in tissues are also studied. The result shows that the overall effects of the blood flow on the thermal response and damage are similar to those of the time delay τT. The two-dimensional numerical results indicate that for a local heating with the heated spot being smaller than the tissue bulk, the variations of the non-uniform distributions of temperature suggest that the multi-dimensional effects of thermal wave and diffusion not be negligible.

2012 ◽  
Vol 134 (5) ◽  
Author(s):  
Hai-Dong Wang ◽  
Bing-Yang Cao ◽  
Zeng-Yuan Guo

Fourier’s law is a phenomenological law to describe the heat transfer process. Although it has been widely used in a variety of engineering application areas, it is still questionable to reveal the physical essence of heat transfer. In order to describe the heat transfer phenomena universally, Guo has developed a general heat conduction law based on the concept of thermomass, which is defined as the equivalent mass of phonon gas in dielectrics according to Einstein’s mass–energy relation. The general law degenerates into Fourier’s law when the thermal inertia is neglected as the heat flux is not very high. The heat flux in carbon nanotubes (CNTs) may be as high as 1012 W/m2. In this case, Fourier’s law no longer holds. However, what is estimated through the ratio of the heat flux to the temperature gradient by molecular dynamics (MD) simulations or experiments is only the apparent thermal conductivity (ATC); which is smaller than the intrinsic thermal conductivity (ITC). The existing experimental data of single-walled CNTs under the high-bias current flows are applied to study the non-Fourier heat conduction under the ultrahigh heat flux conditions. The results show that ITC and ATC are almost equal under the low heat flux conditions when the thermal inertia is negligible, while the difference between ITC and ATC becomes more notable as the heat flux increases or the temperature drops.


2016 ◽  
Vol 20 (5) ◽  
pp. 1457-1472 ◽  
Author(s):  
Poor Ziaei ◽  
Hassan Moosavi ◽  
Amir Moradi

This article focuses on temperature response of skin tissue due to time-dependent surface heat fluxes. Analytical solution is constructed for DPL bio-heat transfer equation with constant, periodic and pulse train heat flux conditions on skin surface. Separation of variables and Duhamel?s theorem for a skin tissue as a finite domain are employed. The transient temperature responses for constant and time-dependent boundary conditions are obtained and discussed. The results show that there is major discrepancy between the predicted temperature of parabolic (Pennes bio-heat transfer), hyperbolic (thermal wave) and DPL bio-heat transfer models when high heat flux accidents on the skin surface with a short duration or propagation speed of thermal wave is finite. The results illustrate that the DPL model reduces to the hyperbolic model when ?T approaches zero and the classic Fourier model when both thermal relaxations approach zero. However for ?q = ?T the DPL model anticipates different temperature distribution with that predicted by the Pennes model. Such discrepancy is due to the blood perfusion term in energy equation. It is in contrast to results from the literature for pure conduction material, where the DPL model approaches the Fourier heat conduction model when ?q = ?T . The burn injury is also investigated.


Author(s):  
Hai-Dong Wang ◽  
Bing-Yang Cao ◽  
Zeng-Yuan Guo

Fourier’s law is a phenomenological law to describe the heat transfer process. Although it has been widely used in a variety of engineering application areas, it is still questionable to reveal the physical essence of heat transfer. In order to describe the heat transfer phenomena universally, Guo has developed a general heat conduction law based on the concept of thermomass, which is defined as the equivalent mass of phonon gas in dielectrics according to Einstein’s mass-energy relation. The general law degenerates into Fourier’s law when the thermal inertia is neglected as the heat flux is not very high. The heat flux in carbon nanotubes (CNTs) may be as high as 1012 W/m2. In this case Fourier’s law no longer holds. However, what is estimated through the ratio of the heat flux to the temperature gradient by MD simulations or experiments is only the apparent thermal conductivity (ATC); which is smaller than the intrinsic thermal conductivity (ITC). The existing experimental data of single-walled CNTs under the high-bias current flows are applied to study the non-Fourier heat conduction under the ultra-high heat flux conditions. The results show that ITC and ATC are almost equal under the low heat flux conditions when the thermal inertia is negligible, while the difference between ITC and ATC becomes more notable as the heat flux increases or the temperature drops.


Author(s):  
Wen-Qiang Lu ◽  
Junfeng Lu

The model of non-Fourier heat conduction in a kind of two-phase mediums with great different heat conductivity is deduced by the idea and mathematics of dual phase lag. It is pointed out that the relaxation times to establish heat flux and temperature gradient include both kinds in this model: the relaxation time appeared under the conditions of applied high heat flux and rapid heating, the relaxation time introduced by the non-equilibrium heat exchange between the two-phase mediums. It is very important to distinguish the both kinds of relaxation times for analyzing and explaining the experimental phenomena of non-Fourier heat conduction in this kind of two-phase mediums.


Author(s):  
А.И. Жмакин

The Fourier law adequately describes heat conduction in most practical problems. However for heat transfer in fast processes and at micro/nano scales, in materials with inner structure (porous media, biological tissues) other models are needed that account for the nonlinear effects and both time (memory) and space nonlocality. Such models including phase lag models, phonon and thermodynamics models as well as fractional derivatives models are reviewed. .


2015 ◽  
Vol 20 (4) ◽  
pp. 443-456 ◽  
Author(s):  
Sonalika Singh ◽  
Sushil Kumar

This paper considers hyperbolic heat conduction model for biological tissue freezing during cryosurgery with non ideal property of tissue, metabolism and blood perfusion. Mathematical model is solved numerically using finite difference method to obtain temperature distribution and phase change interfaces in tissue during freezing. The effects of phase lag of heat flux in hyperbolic bio-heat model on freezing process are studied. Comparative study of parabolic and hyperbolic bio-heat models is also made here.


2012 ◽  
Vol 134 (11) ◽  
Author(s):  
Stephen Akwaboa ◽  
Patrick Mensah ◽  
Ebubekir Beyazouglu ◽  
Ravinder Diwan

This paper presents a numerical solution of the hyperbolic heat conduction equation in a thermal barrier coating (TBC) structure under an imposed heat flux on the exterior of the TBC. The non-Fourier heat conduction equation is used to model the heat conduction in the TBC system that predicts the heat flux and the temperature distribution. This study presents a more realistic approach to evaluate in-service performance of thin layers of TBCs typically found in hot sections of land based and aircraft gas turbine engines. In such ultrafast heat conduction systems, the orders of magnitude of the time and space dimensions are extremely short which renders the traditional Fourier conduction law, with its implicit assumption of infinite speed of thermal propagation, inaccurate. There is, therefore, the need for an advanced modeling approach for the thermal transport phenomenon taking place in microscale systems. A hyperbolic heat conduction model can be used to predict accurately the transient temperature distribution of thermal barrier structures of turbine blades. The hyperbolic heat conduction equations are solved numerically using a new numerical scheme codenamed the mean value finite volume method (MVFVM). The numerical method yields minimal numerical dissipation and dispersion errors and captures the discontinuities such as the thermal wave front in the solution with reliable accuracy. Compared with some traditional numerical methods, the MVFVM method provides the ability to model the behavior of the single phase lag thermal wave following its reflection from domain boundary surfaces. In addition, parametric studies of properties of the substrate on the temperature and the heat flux distributions in the TBC revealed that relaxation time of the substrate material, unlike the thermal diffusivity and thermal conductivity has very little effect on the transient thermal response in the TBC. The study further showed that for thin film structures subject to short time durations of heat flux, the hyperbolic model yields more realistic results than the parabolic model.


Author(s):  
Kyunghan Kim ◽  
Zhixiong Guo

Ultrafast laser radiation heat transfer in biological tissues is governed by time-dependent equation of radiative transfer and modeled using the transient discrete ordinates method. The divergence of radiative heat flux is then obtained and used for predicting the local temperature response due to radiation energy absorption within the ultrashort time period. To this end, the lumped method is employed and heat diffusion is negligible. Both single pulse and pulse train irradiations are considered. For the single pulse irradiation, the transient radiation field is obtained and the local temperature keeps rising until a time of about 20 times of the short pulse width; and then a stable local temperature profile is reached and maintained until the start of heat conduction. For the pulse train case (104 ultrashort pulses until 1 ms), the local temperature response is an accumulation of continuous single pulses because the thermal relaxation time of biological tissues was reported in the range of 1-100 sec and is much longer than the pulse train duration (1 ms). After a stable local temperature field is achieved, the hyperbolic heat conduction model is adopted to describe the heat conduction. MacCormark’s scheme is utilized for solving the thermal wave equations. Thermal wave behavior is observed during the heat transfer process. It is found that the hyperbolic wave model predicts a higher temperature rise than the classical heat diffusion model. After several thermal relaxation times the thermal wave behavior is substantially weakened and the predictions between the hyperbolic and diffusion models match.


Sign in / Sign up

Export Citation Format

Share Document