Local heating of human skin by millimeter waves: effect of blood flow

We investigated the influence of blood perfusion on local heating of the forearm and middle finger skin following 42.25 GHz exposure with an open ended waveguide (WG) and with a YAV mm wave therapeutic device. Both sources had bell-shaped distributions of the incident power density (IPD) with peak intensities of 208 and 55 mW/cm(2), respectively. Blood perfusion was changed in two ways: by blood flow occlusion and by externally applied vasodilator (nonivamide/nicoboxil) cream to the skin. For thermal modeling, we used the bioheat transfer equation (BHTE) and the hybrid bioheat equation (HBHE) which combines the BHTE and the scalar effective thermal conductivity equation (ETCE). Under normal conditions with the 208 mW/cm(2) exposure, the cutaneous temperature elevation (DeltaT) in the finger (2.5 +/- 0.3 degrees C) having higher blood flow was notably smaller than the cutaneous DeltaT in the forearm (4.7 +/- 0.4 degrees C). However, heating of the forearm and finger skin with blood flow occluded was the same, indicating that the thermal conductivity of tissue in the absence of blood flow at both locations was also the same. The BHTE accurately predicted local hyperthermia in the forearm only at low blood flow. The HBHE made accurate predictions at both low and high perfusion rates. The relationship between blood flow and the effective thermal conductivity (k(eff)) was found to be linear. The heat dissipating effect of higher perfusion was mostly due to an apparent increase in k(eff). It was shown that mm wave exposure could result in steady state heating of tissue layers located much deeper than the penetration depth (0.56 mm). The surface DeltaT and heat penetration into tissue increased with enlarging the irradiating beam area and with increasing exposure duration. Thus, mm waves at sufficient intensities could thermally affect thermo-sensitive structures located in the skin and underlying tissue.     

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.     

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.     

Finite-element solution of thermal conductivity of muscle during cold water immersion

The in vivo or effective thermal conductivity (keff) of muscle tissue of the human forearm was determined through a finite-element (FE) model solution of the bioheat equation. Data were obtained from steady-state temperatures measured in the forearm after 3 h of immersion in water at temperatures (Tw) of 15 (n = 6), 20 (n = 5), and 30 degrees C (n = 5). Temperatures were measured every 0.5 cm from the longitudinal axis of the forearm to the skin approximately 9 cm distal from the elbow. Heat flux was measured at two sites on the skin adjacent to the temperature probe. The FE model is comprised of concentric annular compartments with boundaries defined by the location of temperature measurements. Through this approach, it was possible to include both the metabolic heat production and the convective heat transfer between blood and tissue at two levels of blood flow, one perfusing the compartment and the other passing through the compartment. Without heat exchange at the passing blood flow level, the arterial blood temperature would be assumed to have a constant value everywhere in the forearm muscles, leading to a solution of the bioheat equation that greatly underpredicts keff. The extent of convective heat exchange at the passing blood flow level is estimated to be approximately 60% of the total heat exchange between blood and tissue. Concurrent with this heat exchange is a decrease in the temperature of the arterial blood as it flows radially from the axis to the skin of the forearm, and this decrease is enhanced with a lowered Tw.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

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.     

Thermal models for microwave hazards and their role in standards development

We consider the thermal response of the body to radiofrequency (RF) energy, with emphasis on partial-body exposure, to assess potential thermal hazards. The thermal analysis is based on Pennes' bioheat equation. In this model, the thermal response is governed by two time constants. One (τ₁) pertains to heat convection by blood flow and is (for physiologically normal perfusion rates) on the order of 3 min. The second (τ₂) characterizes heat conduction, and varies as the square of a distance that characterizes the spatial extent of the heating. We examine three idealized cases. The first is a region of tissue with an insulated surface, subject to irradiation with an exponentially decreasing SAR, which models a large surface area of tissue exposed to microwaves. The second is a region of tissue in contact with a hemispherical electrode that passes current into it, which models exposure from contact with a conductor. The third is a region of tissue with an insulated surface, subject to heating from a dipole located close to it. In all three cases, we estimate the maximum steady-state temperature increase as a function of the relevant electrical and thermal parameters and the thresholds for thermal hazard. We conclude that thermal models are a potentially fruitful but underutilized means of analyzing thermal hazards from RF fields. A quantitative analysis of such hazards enables the development of data-based uncertainty factors, which can replace arbitrary “safety factors” in developing exposure limits. Finally, we comment on the need to marry quantitative modeling of data and risk assessment, and to incorporate contemporary approaches to risk assessment into RF standards development. Bioelectromagnetics 20:52–63, 1999. © 1999 Wiley-Liss, Inc.     

Reflection and penetration depth of millimeter waves in murine skin

Millimeter (mm) wave reflectivity was used to determine murine skin permittivity. Reflection was measured in anesthetized Swiss Webster and SKH1-hairless mice in the 37-74 GHz frequency range. Two skin models were tested. Model 1 was a single homogeneous skin layer. Model 2 included four skin layers: (1) the stratum corneum, (2) the viable epidermis plus dermis, (3) fat layer, and (4) muscle which had infinite thickness. We accepted that the permittivity of skin in the mm wave frequency range results from the permittivity of cutaneous free water which is described by the Debye equation. Using Fresnel equations for reflection we determined the skin parameters best fitting to the reflection data and derived the permittivity of skin layers. The permittivity data were further used to calculate the power density and specific absorption rate profiles, and the penetration depth of mm waves in the skin. In both murine models, mm waves penetrate deep enough into tissue to reach muscle. In human skin, mm waves are mostly absorbed within the skin. Therefore, when extrapolating the effects of mm waves found in animals to humans, it is important to take into account the possible involvement of muscle in animal effects.     

In vivo thermal conductivity of the human forearm tissues

The effective thermal conductivities of the skin + subcutaneous (keff skin + fat) and muscle (keff muscle) tissues of the human forearm at thermal steady state during immersion in water at temperatures (Tw) ranging from 15 to 36 degrees C were determined. Tissue temperature (Tt) was continuously monitored by a calibrated multicouple probe during a 3-h immersion of the resting forearm. Tt was measured every 5 mm from the longitudinal axis of the forearm (determined from computed-tomography scanning) to the skin surface. Skin temperature (Tsk), heat loss (Hsk), and blood flow (Q) of the forearm, as well as rectal temperature (Tre) and arterial blood temperature at the brachial artery (Tbla), were measured during the experiments. When the keff values were calculated from the finite-element (FE) solution of the bioheat equation, keff skin + fat ranged from 0.28 +/- 0.03 to 0.73 +/- 0.14 W.degrees C-1.m-1 and keff muscle varied between 0.56 +/- 0.05 and 1.91 +/- 0.19 W.degrees C-1.m-1 from 15 to 36 degrees C. The values of keff skin + fat and keff muscle, calculated from the FE solution for Tw less than or equal to 30 degrees C, were not different from the average in vitro values obtained from the literature. The keff values of the forearm tissues were linearly related (r = 0.80, P less than 0.001) to Q for Tw greater than or equal to 30 degrees C. It was found that the muscle tissue could account for 92 +/- 1% of the total forearm insulation during immersion in water between 15 and 36 degrees C.     

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.     

Enhanced absorption of millimeter wave energy in murine subcutaneous blood vessels

The aim of the present study was to determine millimeter wave (MMW) absorption by blood vessels traversing the subcutaneous fat layer of murine skin. Most calculations were performed using the finite-difference time-domain (FDTD) technique. We used two types of models: (1) a rectangular block of multilayer tissue with blood vessels traversing the fat layer and (2) cylindrical models with circular and elliptical cross-sections simulating the real geometry of murine limbs. We found that the specific absorption rate (SAR) in blood vessels normally traversing the fat layer achieved its maximal value at the parallel orientation of the E-field to the vessel axis. At 42 GHz exposure, the maximal SAR in small blood vessels could be more than 30 times greater than that in the skin. The SAR increased with decreasing the blood vessel diameter and increasing the fat thickness. The SAR decreased with increasing the exposure frequency. When the cylindrical or elliptical models of murine limbs were exposed to plane MMW, the greatest absorption of MMW energy occurred in blood vessels located on the lateral areas of the limb model. At these areas the maximal SAR values were comparable with or were greater than the maximal SAR on the front surface of the skin. Enhanced absorption of MMW energy by blood vessels traversing the fat layer may play a primary role in initiating MMW effects on blood cells and vasodilatation of cutaneous blood vessels.     

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.     

Numerical study to establish relationship between coagulation volume and target tip temperature during temperature-controlled radiofrequency ablation

The present study aims at proposing a relationship between the coagulation volume and the target tip temperature in different tissues (viz., liver, lung, kidney, and breast) during temperature-controlled radiofrequency ablation (RFA). A 20-min RFA has been modelled using commercially available monopolar multi-tine electrode subjected to different target tip temperatures that varied from 70°C to 100°C with an increment of 10°C. A closed-loop feedback proportional-integral-derivative (PID) controller has been employed within the finite element model to perform temperature-controlled RFA. The coagulation necrosis has been attained by solving the coupled electric field distribution, the Pennes bioheat and the first-order Arrhenius rate equations within the three-dimensional finite element model of different tissues. The computational study considers temperature-dependent electrical and thermal conductivities along with the non-linear piecewise model of blood perfusion. The comparison between coagulation volume obtained from the numerical and in vitro experimental studies has been done to evaluate the aptness of the numerical models. In the present study, a total of 20 numerical simulations have been performed along with 12 experiments on tissue-mimicking phantom gel using RFA device. The study revealed a strong dependence of the coagulation volume on the pre-set target tip temperature and ablation time during RFA application. Further, the effect of target tip temperature on the applied input voltage has been studied in different tissues. Based on the results attained from the numerical study, statistical correlations between the coagulation volume and treatment time have been developed at different target tip temperatures for each tissue.     

A new simplified bioheat equation for the effect of blood flow on local average tissue temperature

A new simplified three-dimensional bioheat equation is derived to describe the effect of blood flow on blood-tissue heat transfer. In two recent theoretical and experimental studies [1, 2] the authors have demonstrated that the so-called isotropic blood perfusion term in the existing bioheat equation is negligible because of the microvascular organization, and that the primary mechanism for blood-tissue energy exchange is incomplete countercurrent exchange in the thermally significant microvessels. The new theory to describe this basic mechanism shows that the vascularization of tissue causes it to behave as an anisotropic heat transfer medium. A remarkably simple expression is derived for the tensor conductivity of the tissue as a function of the local vascular geometry and flow velocity in the thermally significant countercurrent vessels. It is also shown that directed as opposed to isotropic blood perfusion between the countercurrent vessels can have a significant influence on heat transfer in regions where the countercurrent vessels are under 70-micron diameter. The new bioheat equation also describes this mechanism.     

Elderly bioheat modeling: changes in physiology,thermoregulation, and blood flow circulation

A bioheat model for the elderly was developed focusing on blood flow circulatory changes that influence their thermal response in warm and cold environments to predict skin and core temperatures for different segments of the body especially the fingers. The young adult model of Karaki et al. (Int J Therm Sci 67:41–51, 2013) was modified by incorporation of the physiological thermoregulatory and vasomotor changes based on literature observations of physiological changes in the elderly compared to young adults such as lower metabolism and vasoconstriction diminished ability, skin blood flow and its minimum and maximum values, the sweating values, skin fat thickness, as well as the change in threshold parameter related to core or skin temperatures which triggers thermoregulatory action for sweating, maximum dilatation, and maximum constriction. The developed model was validated with published experimental data for elderly exposure to transient and steady hot and cold environments. Predicted finger skin temperature, mean skin temperature, and core temperature were in agreement with published experimental data at a maximum error less than 0.5 °C in the mean skin temperature. The elderly bioheat model showed an increase in finger skin temperature and a decrease in core temperature in cold exposure while it showed a decrease in finger skin temperature and an increase in core temperature in hot exposure.     

Millimeter wave dosimetry of human skin

To identify the mechanisms of biological effects of mm waves it is important to develop accurate methods for evaluating absorption and penetration depth of mm waves in the epidermis and dermis. The main characteristics of mm wave skin dosimetry were calculated using a homogeneous unilayer model and two multilayer models of skin. These characteristics included reflection, power density (PD), penetration depth (delta), and specific absorption rate (SAR). The parameters of the models were found from fitting the models to the experimental data obtained from measurements of mm wave reflection from human skin. The forearm and palm data were used to model the skin with thin and thick stratum corneum (SC), respectively. The thin SC produced little influence on the interaction of mm waves with skin. On the contrary, the thick SC in the palm played the role of a matching layer and significantly reduced reflection. In addition, the palmar skin manifested a broad peak in reflection within the 83-277 GHz range. The viable epidermis plus dermis, containing a large amount of free water, greatly attenuated mm wave energy. Therefore, the deeper fat layer had little effect on the PD and SAR profiles. We observed the appearance of a moderate SAR peak in the therapeutic frequency range (42-62 GHz) within the skin at a depth of 0.3-0.4 mm. Millimeter waves penetrate into the human skin deep enough (delta = 0.65 mm at 42 GHz) to affect most skin structures located in the epidermis and dermis.     

Effects of fibroblasts and microenvironment on epidermal regeneration and tissue function in long-term skin equivalents

In vitro generated skin models find growing interest as promising tools in basic research and clinical application in regenerative medicine. Here, we present further details of an improved long-term skin equivalent (SE) enabling mechanistic studies on skin reconstruction and epidermal function. Growth conditions of fibroblasts in a 3D scaffold were analysed to optimise the dermal microenvironment by providing an authentic dermal matrix for regular tissue reconstruction and function of cocultured keratinocytes. These SEs demonstrate sustained epidermal viability - over 12 weeks - with regular differentiation as substantiated by in vivo-like patterns of all differentiation products, exemplified here by the cornified envelope components loricrin and repetin. The continuous expression of all major tight junction components in the granular layer, shown here for ZO-1 in coherence with the presence of epidermal barrier lipids, and ultrastructural accumulation of lamellar bodies, collectively indicate proper epidermal barrier structures. Remarkably, cocultured keratinocytes exerted an ongoing proliferation-stimulating effect on fibroblasts colonising the scaffold comparable to a cocktail of fibroblast growth factors. Consequently, precultivation of dermal equivalents (DEs) in basal or growth factor-enriched media had only minor effects on the quality of epidermal regeneration in cocultures. As to the role of fibroblast numbers, complete absence of dermal cells resulted in atrophic epithelia but the effect of cell numbers as low as 5 x 10(4)cells/cm(2) on epidermal tissue quality equalled that of the standard density (2 x 10(5)cells/cm(2)). Surprisingly, precultivation of fibroblasts in the DEs for 7 days (standard) showed no better effect on epidermal tissue reformation as compared to 2 days whereas a precultivation period of 14 days resulted in atrophic epidermal and dermal tissue development. These data demonstrate, (i) the strict dependence of epidermal tissue regeneration on the presence of fibroblasts, (ii) the mutual keratinocyte-fibroblast interactions for cell proliferation and organogenesis, and (iii) the importance of the proper microenvironment for epidermal tissue function and supposedly for establishment of a stem cell niche in vitro.     

The effects of large blood vessels on temperature distributions during simulated hyperthermia.

Several three-dimensional vascular models have been developed to study the effects of adding equations for large blood vessels to the traditional bioheat transfer equation of Pennes when simulating tissue temperature distributions. These vascular models include "transiting" vessels, "supplying" arteries, and "draining" veins, for all of which the mean temperature of the blood in the vessels is calculated along their lengths. For the supplying arteries this spatially variable temperature is then used as the arterial temperature in the bioheat transfer equation. The different vascular models produce significantly different locations for both the maximum tumor and the maximum normal tissue temperatures for a given power deposition pattern. However, all of the vascular models predict essentially the same cold regions in the same locations in tumors: one set at the tumors' corners and another around the inlets of the large blood vessels to the tumor. Several different power deposition patterns have been simulated in an attempt to eliminate these cold regions; uniform power in the tumor, annular power in the tumor, preheating of the blood in the vessels while they are traversing the normal tissue, and an "optimal" power pattern which combines the best features of the above approaches. Although the calculations indicate that optimal power deposition patterns (which improve the temperature distributions) exist for all of the vascular models, none of the heating patterns studied eliminated all of the cold regions. Vasodilation in the normal tissue is also simulated to see its effects on the temperature fields.(ABSTRACT TRUNCATED AT 250 WORDS)     

人胎儿真皮乳头和皮纹发生的研究

本文用扫描电镜法研究了入胎儿的皮纹发生过程,包括初级真皮嵴、次级真皮嵴、真皮乳头和表皮隆线的发生。为研究人皮纹的发生和皮肤的异常提供了皮肤正常发育的形态学依据。共观察111例从第6周到第9个月胎儿的皮纹区皮肤,表明第3个月末胎儿开始形成初级真皮嵴,以后逐渐加深,至第16周嵴的顶端中央产生纵沟形成两条平行的次级真皮嵴;自19周后,次级真皮嵴局部隆起,由波浪形逐渐形成乳头。至30周乳头呈犬牙状。表皮隆线于第4—5月形成,随真皮乳头的增高而渐趋明显。至第6个月,全部皮纹图样已可辨认。本文还讨论了真皮乳头发生的过程。