Abstract

We consider a problem of skin tissue of the human head subjected to thermal diffusion with relaxation time. The bounding 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.

1. 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, 1420]. 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].

2. Formulation of the Problem

We consider a small bounded cross section of the skin tissue of the human head occupies the region subjected to thermal diffusion (such as a gas) in the upper plane of the skin tissue 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, is the external heat source, is the heat source due to metabolic heating, is a measure of thermodiffusion effect, 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 is the diffusion relaxation time, is the diffusion coefficient, and is the chemical potential which satisfieswhere is a measure of diffusive effect.

To simplify the equations in (1)–(3), we will use the following nondimension variables:After applying the dimensionless variables, (1)–(3) take the formswhere and .

We dropped the prime for simplicity.

We will consider the tissue has no external heat source and the heat source due to metabolic heating is constant; then we getBy eliminating between (6) and (7), we obtainWe will apply Laplace transform on (6) and (8) which is defined aswhere .

Hence, we getBy eliminating between the equations in (11), we obtainand by eliminating between the equations in (11), we getwhere .

We can write the general solutions of (12) and (13) in the formswhere are parameters and , are the roots of the characteristic equationwhere .

To get the parameters , we will consider the upper plane of the tissue is subjected to thermal and chemical potential shock which depends on time only as follows:Applying Laplace transform on conditions (16), we getBy using (7) and (17), we getApplying the conditions in (17) and (18) on (14), we obtain the following system of linear equations:By solving the system above, we getThus, we have the solution of the problem in the Laplace transform domain in the form

3. 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 .

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

4. The Numerical Values

The following properties of a skin tissue are used in the numerical analysis [15]:

Figures 2 and 3 represent the temperature increment distribution and the concentration distribution, respectively, with a wide range of the thickness of the layer of the tissue when and at different values of time . Those figures show that the value of the time has a significant effect 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. 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 and at different values of distance . 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.

Conflicts of Interest

The funding mentioned in Acknowledgments did not lead to any conflicts of interest regarding the publication of this manuscript.

Acknowledgments

The author is grateful for the support and the funding for this work provided by Deanship of Scientific Research at Najran University, Saudi Arabia (Grant no. NU/ESCI/15/042).