The Skin Tissue of the Human Head Subjected to Thermal Diffusion

We consider a problemof skin tissue of the humanhead subjected to thermal diffusionwith relaxation time.Thebounding surface of the tissue of the human head is subjected to a thermal shock.The chemical potential is also assumed to be a known function of time on the bounding plane. Laplace transform techniques are used and the results are obtained numerically using a numerical method for the inversion of the Laplace transform based on Tzou method.The temperature distribution and the concentration distribution are obtained and represented in figures. According to the results, the values of the time and the distance have significant effects on the temperature increment and the concentration. The values of the temperature increment and the concentration increase when the time increases and decrease when the distance increases.


Introduction
Biomechanics is the area of science where the laws, principles, and methods of mechanics are applied to the structure and function of the human body.Mechanics can be divided into two categories: statics, which is the study of stationary objects, and dynamics, which is the study of moving objects [1].
Some human trial studies have expressed that, in the short term, being subjected to cellular phones frequencies can cause such problems like exciting the brain cells when subjected to a wave [1], the possible occurrence of cataracts with continuous exposure weakening the neural behavioral and cognitive functioning, and possibly a disruption in the functionality of the nervous system [1].Lindholm et al. [2] conducted a human trial to see the effects of short term radio frequency when the electromagnetic fields are subjected to preadolescents.
Ferreri et al. studied the mobile phone emissions and human brain excitability [3] and Keetley et al. studied the neuropsychological squeal of digital mobile phone exposure in humans [4], while Lin discussed the cataracts and personal communication radiation [5].The effects of electromagnetic radiation of mobile phones on the central nervous system have been studied by Hossmann and Hermann [6] and the thermal effects of radiation from cellular telephones have been studied by Wainwright [7].
The most important work regarding specific absorption rate and temperature distributions in human head subjected to mobile phone radiation at different frequencies was introduced by Wessapan et al. [8].Elwassif et al. constructed the bioheat transfer model of deep brain stimulation-induced temperature changes [9].Ley and Bayazitoglu analyzed the brain temperature distribution during deep hypothermic circulatory arrest in humans [10] and got the effect of physiology on the temperature distribution of a layered head with external convection [1].
The thermal wave model of bioheat transfer (TWMBT) equation is an energy balance equation that includes the heat caused by conduction, convection, metabolic heat generation, and the heat produced by radiation due to the cellular frequency [11].The model incorporates a relaxation time to account for the lag time it takes for the radio frequency wave to respond to the heat penetrating through the tissue, that is, heat flux associated with the tissues [12,13].Biological experimental studies have shown that this value can range from 15 to 30 s [12,[14][15][16][17][18][19][20].An additional relaxation time accommodating the phase lag of the temperature gradient is analytically studied.Morimoto et al. [21] studied the time constants of the temperature elevations in human head and body models exposed to simulated radiation from dipole antennas, electromagnetic beams, and plane waves within frequency range considered to be from 1 to 30 GHz.Tang et al. discussed the thermal and vibrational characterization of human skin to find suitable and reproducible biomarkers for further clinical studies on intrinsic and extrinsic aging of dermis [22].

Formulation of the Problem
We consider a small bounded cross section of the skin tissue of the human head occupies the region 0 ≤  ≤ ℓ subjected to thermal diffusion (such as a gas) in the upper plane of the skin tissue  = 0 as in Figure 1, and this tissue is initially quiescent.
The heat transfer in living tissue is accompanied by blood perfusion and metabolic heat generation; the bioheat equation takes the following form [4]: where   =   (, ) is the temperature of the tissue,   is the blood temperature (constant),  is the thermal conductivity,  is the density of the skin tissue,   is the blood mass density,   is the specific heat of the blood,   is the specific heat of the skin tissue,   is the blood perfusion rate,  out is the external heat source,  met is the heat source due to metabolic heating,  is a measure of thermodiffusion effect,  1 is the relaxation time of the skin tissue, and  = (, ) is the concentration of the diffusive material in the skin tissue.The diffusion equations take the form where  2 is the diffusion relaxation time,  is the diffusion coefficient, and  = (, ) is the chemical potential which satisfies where  is a measure of diffusive effect.

Inverse Laplace Transforms
To get solutions in the time domain, the Riemann-sum approximation method is used to obtain the numerical results.In this method, any function in Laplace domain can be inverted to the time domain as follows: where Re is the real part and  is the imaginary number unit.For faster convergence, many numerical experiments have shown that the value of  satisfies the relation  ≈ 4.7.
We considered that the functions on the boundary in ( 16) will be defined as a shock in the following forms: where  is constant (the intensity of the shock) and () is the Heaviside unit step function.By using Laplace transform, (23) takes the form

The Numerical Values
The following properties of a skin tissue are used in the numerical analysis [1][2][3][4][5]: the concentration.The values of the temperature increment and the concentration increase when the time increases and decrease when the distance increases.The thermal and the concentration waves have finite speeds of wave propagation which agrees with the real physical behavior.
Figures 4 and 5 represent the temperature increment distribution and the concentration distribution, respectively, with a wide range of the time when 0 ≤  ≤ 1.0 and at different values of distance  = (0.2, 0.4, 0.6).Those figures show that the value of the distance has a significant effect on the temperature increment and the concentration.The values of the temperature increment and the concentration increase when the distance decreases and go to the steady state position when the time increases.Moreover, the steady state positions start first in the nearest distance.The steady states of the thermal wave and the concentration wave occur almost at the same instance.

7 Figure 2 :Figures 2 Figure 3 :Figure 4 :
Figure 2: The temperature increment distribution at different values of time.

Figure 5 :
Figure 5: The concentration distribution at different values of distance.