Numerical simulation of time fractional dual-phase-lag model of heat transfer within skin tissue during thermal therapy

In this paper, we investigated the thermal behavior in living biological tissues using time fractional dual-phase-lag bioheat transfer (DPLBHT) model subjected to Dirichelt boundary condition in presence of metabolic and electromagnetic heat sources during thermal therapy. We solved this bioheat transfer model using finite element Legendre wavelet Galerkin method (FELWGM) with help of block pulse function in sense of Caputo fractional order derivative. We compared the obtained results from FELWGM and exact method in a specific case, and found a high accuracy. Results are interpreted in the form of standard and anomalous cases for taking different order of time fractional DPLBHT model. The time to achieve hyperthermia position is discussed in both cases as standard and time fractional order derivative. The success of thermal therapy in the treatment of metastatic cancerous cell depends on time fractional order derivative to precise prediction and control of temperature. The effect of variability of parameters such as time fractional derivative, lagging times, blood perfusion coefficient, metabolic heat source and transmitted power on dimensionless temperature distribution in skin tissue is discussed in detail. The physiological parameters has been estimated, corresponding to the value of fractional order derivative for hyperthermia treatment therapy.     

Theoretical analysis of thermal damage in biological tissues caused by laser irradiation

A bioheat transfer approach is proposed to study thermal damage in biological tissues caused by laser radiation. The laser light propagation in the tissue is first solved by using a robust seven-flux model in cylindrical coordinate system. The resulting spatial distribution of the absorbed laser energy is incorporated into the bioheat transfer equation for solving temperature response. Thermal damage to the tissue is assessed by the extent of denatured protein using a rate process equation. It is found that for the tissue studied, a significant protein denaturation process would take place when temperature exceeds about 53 degrees C. The effects of laser power, exposure time and beam size as well as the tissue absorption and scattering coefficients on the thermal damage process are examined and discussed. The laser conditions that cause irreversible damage to the tissue are also identified.     

Analytical study on bioheat transfer problems with spatial or transient heating on skin surface or inside biological bodies

Several closed form analytical solutions to the bioheat transfer problems with space or transient heating on skin surface or inside biological bodies were obtained using Green's function method. The solutions were applied to study several selected typical bioheat transfer processes, which are often encountered in cancer hyperthermia, laser surgery, thermal comfort analysis, and tissue thermal parameter estimation. Thus a straightforward way to quantitatively interpret the temperature behavior of living tissues subject to constant, sinusoidal, step, point or stochastic heatings etc. both in volume and on boundary were established. Further solution to the three-dimensional bioheat transfer problems was also given to illustrate the versatility of the present method. Implementations of this study to the practical problems were addressed.     

A thermal model for human thresholds of microwave-evoked warmth sensations

Human thresholds for skin sensations of warmth were measured at frequencies from 2.45 to 94 GHz. By solving the one-dimensional bioheat equation, we calculated the temperature increase at the skin surface or at a depth of 175 μm at incident power levels corresponding to the observed thresholds. The thermal analysis suggests that the thresholds correspond to a localized temperature increase of about 0.07 °C at and near the surface of the skin. We also found that, even at the highest frequency of irradiation, the depth at which the temperature receptors are located is not a relevant parameter, as long as it is within 0.3 mm of the surface. Over the time range of the simulation, the results of the thermal model are insensitive to blood flow, but sensitive to thermal conduction; and this sensitivity increases strongly with frequency. We conclude with an analysis of the effect of thermal conduction on surface temperature rise, which becomes a dominant factor at microwave frequencies over 10 GHz. Bioelectromagnetics 18:578–583, 1997. © 1997 Wiley-Liss, Inc.     

Experimental validation of an inverse heat transfer algorithm for optimizing hyperthermia treatments

Hyperthermia is a cancer treatment modality in which body tissue is exposed to elevated temperatures to destroy cancerous cells. Hyperthermia treatment planning refers to the use of computational models to optimize the heating protocol with the goal of isolating thermal damage to predetermined treatment areas. This paper presents an algorithm to optimize a hyperthermia treatment protocol using the conjugate gradient method with the adjoint problem. The output of the minimization algorithm is a heating protocol that will cause a desired amount of thermal damage. The transient temperature distribution in a cylindrical region is simulated using the bioheat transfer equation. Temperature and time are integrated to calculate the extent of thermal damage in the region via a first-order rate process based on the Arrhenius equation. Several validation experiments are carried out by applying the results of the minimization algorithm to an albumen tissue phantom. Comparisons of metrics describing the damage region (the height and radius of the volume of thermally ablated phantom) show good agreement between the desired extent of damage and the measured extent of damage. The sensitivity of the bioheat transfer model and the Arrhenius damage model to their constituent parameters is calculated to create a tolerable range of error between the desired and measured extent of damage. The measured height and radius of the ablated region fit well within the tolerable range of error found in the sensitivity analysis.     

Experimental study on characteristic of bioheat transfer and numerical simulation for the temperature fields in transverse and longitudinal sections of pig tongue

The influence factors of the tongue tissue structure on the temperature fields are importantly laid on the diameters and positions of arteries. Then temperature fields of certain sections are simulated by the finite element method (FEM). The results show that temperature gradients of the tongue in length are wider than in depth as a result of the specific vascularity. Thereby the relationship between temperature fields of tongue surface and those in transverse and longitudinal sections can be developed. This work will be further referential to study on the characteristics of heat transfer in human tongue, and will promote the development of bioheat transfer science.     

Heat transfer analysis of skin during thermal therapy using thermal wave equation

Specifying exact geometry of vessel network and its effect on temperature distribution in living tissues is one of the most complicated problems of the bioheat field. In this paper, the effects of blood vessels on temperature distribution in a skin tissue subjected to various thermal therapy conditions are investigated. Present model consists of counter-current multilevel vessel network embedded in a three-dimensional triple-layered skin structure. Branching angles of vessels are calculated using the physiological principle of minimum work. Length and diameter ratios are specified using length doubling rule and Cube law, respectively. By solving continuity, momentum and energy equations for blood flow and Pennes and modified Pennes bioheat equations for the tissue, temperature distributions in the tissue are measured. Effects of considering modified Pennes bioheat equation are investigated, comprehensively. It is also observed that blood has an impressive role in temperature distribution of the tissue, especially at high temperatures. The effects of different parameters such as boundary conditions, relaxation time, thermal properties of skin, metabolism and pulse heat flux on temperature distribution are investigated. Tremendous effect of boundary condition type at the lower boundary is noted. It seems that neither insulation nor constant temperature at this boundary can completely describe the real physical phenomena. It is expected that real temperature at the lower levels is somewhat between two predicted values. The effect of temperature on the thermal properties of skin tissue is considered. It is shown that considering temperature dependent values for thermal conductivity is important in the temperature distribution estimation of skin tissue; however, the effect of temperature dependent values for specific heat capacity is negligible. It is seen that considering modified Pennes equation in processes with high heat flux during low times is significant.     

Uncertainty analysis for temperature prediction of biological bodies subject to randomly spatial heating

An analytical solution to the Pennes bioheat transfer equation in three-dimensional geometry with practical hyperthermia boundary conditions and random heating was obtained in this paper. Uncertainties for the predicted temperatures of tissues due to approximate parameters were studied based on analyzing one-dimensional heat transfer in the biological bodies subject to a spatially decay heating. Contributions from each of the thermal parameters such as heat conductivity, blood perfusion rate, and metabolic rate of the tissues, the scattering coefficient and the surface power flux of the heating apparatus were compared and the uncertainty limit for temperature distribution in this case was estimated. The results are useful in a variety of clinical hyperthermia and biological thermal parameter measurement.     

Conventional and newly developed bioheat transport models in vascularized tissues: A review

Heat transfer in a biological system is a complex process and its analysis is difficult. Heterogeneous vascular architecture, blood flow in the complex network of arteries and veins, varying metabolic heat generation rates and dependence of tissue properties on its physiological condition contribute to this complexity. The understanding of heat transfer in human body is important for better insight of thermoregulatory mechanism and physiological conditions. Its understanding is also important for accurate prediction of thermal transport and temperature distribution during biomedical applications. During the last three decades, many attempts have been made by researchers to model the complex thermal behavior of the human body. These models, viz., blood perfusion, countercurrent, thermal phase-lag, porous-media, perturbation, radiation, etc. have their corresponding strengths and limitations. Along with their biomedical applications, this article reviews various contextual issues associated with these models. After brief discussion of early bioheat models, the newly developed bioheat models are discussed in detail. Dependence of these models on biological properties, viz., thermophysical and optical properties are also discussed.     

A parametric study of thermal therapy of skin tissue

A thermal therapy for cancer in skin tissue is numerically investigated using three bioheat conduction models, namely Pennes, thermal wave and dual-phase lag models. A laser is applied at the surface of the skin for cancer ablation, and the temperature and thermal damage distributions are predicted using the three bioheat models and two different modeling approaches of the laser effect. The first one is a prescribed surface heat flux, in which the tissue is assumed to be highly absorbent, while the second approach is a volumetric heat source, which is reasonable if the scattering and absorption skin effects are of similar magnitude. The finite volume method is applied to solve the governing bioheat equation. A parametric study is carried out to ascertain the effects of the thermophysical properties of the cancer on the thermal damage. The temperature distributions predicted by the three models exhibit significant differences, even though the temperature distributions are similar when the laser is turned off. The type of bioheat model has more influence on the predicted thermal damage than the type of modeling approach used for the laser. The phase lags of heat flux and temperature gradient have an important influence on the results, as well as the thermal conductivity of the cancer. In contrast, the uncertainty in the specific heat and blood perfusion rate has a minor influence on the thermal damage.     

Pulsatile blood flow effects on temperature distribution and heat transfer in rigid vessels

The effect of blood velocity pulsations on bioheat transfer is studied. A simple model of a straight rigid blood vessel with unsteady periodic flow is considered. A numerical solution that considers the fully coupled Navier-Stokes and energy equations is used for the simulations. The influence of the pulsation rate on the temperature distribution and energy transport is studied for four typical vessel sizes: aorta, large arteries, terminal arterial branches, and arterioles. The results show that: the pulsating axial velocity produces a pulsating temperature distribution; reversal of flow occurs in the aorta and in large vessels, which produces significant time variation in the temperature profile. Change of the pulsation rate yields a change of the energy transport between the vessel wall and fluid for the large vessels. For the thermally important terminal arteries (0.04-1 mm), velocity pulsations have a small influence on temperature distribution and on the energy transport out of the vessels (8 percent for the Womersley number corresponding to a normal heart rate). Given that there is a small difference between the time-averaged unsteady heat flux due to a pulsating blood velocity and an assumed nonpulsating blood velocity, it is reasonable to assume a nonpulsating blood velocity for the purposes of estimating bioheat transfer.     

A comparative study of thermal effects of 3 types of laser in eye: 3D simulation with bioheat equation

The application of laser in ophthalmology and eye surgery is so widespread that hardly can anyone deny its importance. On the other hand, since the human eye is an organ susceptible to external factors such as heat waves, laser radiation rapidly increases the temperature of the eye and therefore the study of temperature distribution inside the eye under laser irradiation is crucial; but the use of experimental and invasive methods for measuring the temperature inside the eye is typically high-risk and hazardous. In this paper, using the three-dimensional finite element method, the distribution of heat transfer inside the eye under transient condition was studied through three different lasers named Nd:Yag, Nd:Yap and ArF. Considering the metabolic heat and blood perfusion rate in various regions of the eye, numerical solution of space–time dependant Pennes bioheat transfer equation has been applied in this study. Lambert–Beer's law has been used to model the absorption of laser energy inside the eye tissues. It should also be mentioned that the effect of the ambient temperature, tear evaporation rate, laser power and the pupil diameter on the temperature distribution have been studied. Also, temperature distribution inside the eye after applying each laser and temperature variations of six optional regions as functions of time have been investigated. The results show that these radiations cause temperature rise in various regions, which will in turn causes serious damages to the eye tissues. Investigating the temperature distribution inside the eye under the laser irradiation can be a useful tool to study and predict the thermal effects of laser radiation on the human eye and evaluate the risk involved in performing laser surgery.     

Heat transport mechanisms in vascular tissues: a model comparison

We have conducted a parametric comparison of three different vascular models for describing heat transport in tissue. Analytical and numerical methods were used to predict the gross temperature distribution throughout the tissue and the small-scale temperature gradients associated with thermally significant blood vessels. The models are: an array of unidirectional vessels, an array of countercurrent vessels, and a set of large vessels feeding small vessels which then drain into large vessels. We show that three continuum formulations of bioheat transfer (directed perfusion, effective conductivity, and a temperature-dependent heat sink) are limiting cases of the vascular models with respect to the thermal equilibration length of the vessels. When this length is comparable to the width of the heated region of tissue, the local temperature changes near the vessels can be comparable to the gross temperature elevation. These results are important to the use of thermal techniques used to measure the blood perfusion rate and in the treatment of cancer with local hyperthermia.     

Studies on the three-dimensional temperature transients in the canine prostate during transurethral microwave thermal therapy

Thermal therapy of benign prostatic hyperplasia requires accurate prediction of the temperature distribution induced by the heating within the prostatic tissue. In this study, the Pennes bioheat transfer equation was used to model the transient heat transfer inside the canine prostate during transurethral microwave thermal therapy. Incorporating the specific absorption rate of microwave energy in tissue, a closed-form analytical solution was obtained. Good agreement was found between the theoretical predictions and in-vivo experimental results. Effects of blood perfusion and the cooling at the urethral wall on the temperature rise were investigated within the prostate during heating. The peak intraprostatic temperatures attained by application of 5, 10, or 15 W microwave power were predicted to be 38 degrees C, 41 degrees C, and 44 degrees C. Results from this study will help optimize the thermal dose that can be applied to target tissue during the therapy.     

Application of Different Variants of the BEM in Numerical Modeling of Bioheat Transfer Problems

Heat transfer processes proceeding in the living organisms are described by the different mathematical models. In particular, the typical continuous model of bioheat transfer bases on the most popular Pennes equation, but the Cattaneo-Vernotte equation and the dual phase lag equation are also used. It should be pointed out that in parallel are also examined the vascular models, and then for the large blood vessels and tissue domain the energy equations are formulated separately. In the paper the different variants of the boundary element method as a tool of numerical solution of bioheat transfer problems are discussed. For the steady state problems and the vascular models the classical BEM algorithm and also the multiple reciprocity BEM are presented. For the transient problems connected with the heating of tissue, the various tissue models are considered for which the 1st scheme of the BEM, the BEM using discretization in time and the general BEM are applied. Examples of computations illustrate the possibilities of practical applications of boundary element method in the scope of bioheat transfer problems.     

Heating of tissues by microwaves: A model analysis

We consider the thermal response times for heating of tissue subject to nonionizing (microwave or infrared) radiation. The analysis is based on a dimensionless form of the bioheat equation. The thermal response is governed by two time constants: one(τ₁) pertains to heat convection by blood flow, and is of the order of 20–30 min for physiologically normal perfusion rates; the second (τ₂) characterizes heat conduction and varies as the square of a distance that characterizes the spatial extent of the heating. Two idealized cases are examined. The first is a tissue block with an insulated surface, subject to irradiation with an exponentially decreasing specific absorption rate, which models a large surface area of tissue exposed to microwaves. The second is a hemispherical region of tissue exposed at a spatially uniform specific absorption rate, which models localized exposure. In both cases, the steady-state temperature increase can be written as the product of the incident power density and an effective time constant τ_eff, which is defined for each geometry as an appropriate function of τ₁ and τ₂. In appropriate limits of the ratio of these time constants, the local temperature rise is dominated by conductive or convective heat transport. Predictions of the block model agree well with recent data for the thresholds for perception of warmth or pain from exposure to microwave energy. Using these concepts, we developed a thermal averaging time that might be used in standards for human exposure to microwave radiation, to limit the temperature rise in tissue from radiation by pulsed sources. We compare the ANSI exposure standards for microwaves and infrared laser radiation with respect to the maximal increase in tissue temperature that would be allowed at the maximal permissible exposures. A historical appendix presents the origin of the 6-min averaging time used in the microwave standard. Bioelectromagnetics 19:420–428, 1998. © 1998 Wiley-Liss, Inc.     

Theoretical temperature distributions produced by an annular phased array-type system in CT-based patient models

Theoretical calculations for the specific absorption rate (SAR) and the resulting temperature distributions produced by an annular phased array (APA)-type system are made. The finite element numerical method is used in the formulation of both the electromagnetic (EM) and thermal boundary value problems. A number of detailed two-dimensional patient models based on CT-scan data from the pelvic, visceral, and thoracic regions are generated to simulate a variety of tumor locations and surrounding normal tissues. The SAR values from the EM solution are put into the bioheat transfer equation, and steady-state temperature distributions are calculated for a wide range of blood flow rates. Based on our theoretical modeling, the APA shows no preferential heating of superficial over deep-seated tumors. However, in most cases for all three regions of the human trunk only fair thermal profiles (therapeutic area near 60%) are obtained in tumors with little or no blood flow and poor temperature patterns (therapeutic area less than 50%) are found in tumors with moderate to high perfusion rates. These theoretical calculations should aid the clinician in the evaluation of the effectiveness of APA-type devices in heating tumors located in the trunk region.     

Laser-induced hyperthermia of nanoshell mediated vascularized tissue – A numerical study

Laser-induced hyperthermia treatment of tumor in a 2-D axisymmetric tissue embedded with moderate size (100–150 µm) blood vessels is studied. Laser absorption is enhanced by embedding gold–silica nanoshells in the tumor. Heat transfer in the tissue is modeled using Weinbaum–Jiji bioheat transfer equation. With laser irradiation, the volumetric radiation is accounted in the governing bioheat equation. Radiative information needed in the bioheat equation is calculated using the discrete ordinate method, and the coupled bioheat-radiation equation is solved using the finite volume method. Effects of power density, laser exposure time, beam radius, diameter of blood vessel and volume fractions of nanoshells on temperature spread in the tissue are analyzed.