Abstract

This research deals with the temperature increment and responsiveness of skin tissue to a continuous flow of surface heat induced by a constant-voltage electrical current. The exact analytical solution for the dual-phase-lag (DPL) of bioheat transfer has been obtained. It is used to confine the variables to a limited domain to solve the governing equations. The transition temperature reactions have been measured and investigated. The figures provide a comparison of the Pennes, Tzou models, and Vernotte–Cattaneo models. The numerical results demonstrate the values of the voltage, resistance, electric shock time, and dual-phase-lag time parameters which have significant influences on the distributions of the dynamic and conductive temperature rise through the skin tissue.

1. Introduction

Many methods of treatment and application rely substantially on temperature transfer through living skin tissue [1]. Lasers, microwaves, and other technological breakthroughs have also contributed to the progress of biothermal transfer. Several experts, including Pennes’s works, weigh in on bioheat transmission research in living biological tissues [2]. He refined the parabolic equation of the model for living tissues and introduced the first bioheat transfer model. Pennes’s biological transmission technique has also been used to describe the consistency and nature of temperature action in live bodies and biological tissues. On the exterior, such irregular results demonstrate a hyperbolic conduction behaviour which is called the non-Fourier model of heat conduction. Vernotte and Cattaneo’s (V-C) multiple modulations of the heat conduction law were given as a linear thermal extension form of the well-known Fourier law to describe this kind of hyperbolic differential equation. To analyze the influence of microwave and flux thermal activity, they created a thermal wave model [3, 4]. Many methods were tried to cure different types of skin tissue problems without hurting the good tissue around or near them. Xu et al. obtained the solution of Pennes’s bioheat transfer equation (PBTE) analytically. In addition, they investigated skin bioheat transfer, skin biomechanics, thermal injury, and skin structure [5]. The dual-phase-lag, hyperbolic, and parabolic biomass transport models were used to investigate the thermomechanical activity based on the non-Fourier skin tissues under different surface thermal loading constraints. Xu and colleagues found substantial discrepancies between Pennes’s thermal wave and dual-phase-lag (DPL) anticipation models [6]. Furthermore, Rossmanna and Haemmerich investigated how the thermal characteristics, perfusion, and dielectric properties of biological tissues vary with temperature at hyperthermic and ablation temperatures [7]. Moradi et al. have focused on the temperature change in skin tissue caused by heating with a time-dependent surface [8]. The model of Tzou was expanded by using the dual-phase-lag (DPL) technique, which considers the delayed activity in a high rate of reaction. Despite the lengthy operation, the small-scale response was discovered in good time [9, 10]. O¨zis¸ik and Tzou have introduced phase lag behaviour for a temperature gradient [11]. Askarizadeh and Ahmadikia used the dual-phase-lag (DPL) model for solving transient heat transmission difficulties in skin tissue [12].

Non-Fourier thermal activity was used in the Liu and Chen heat conduction of the dual-phase-lag (DPL) model to assess the thermal transit in biological tissue and to diagnose hyperthermia [13]. Zhang defined phase lag or relaxation periods in terms of tissue and blood properties, the interphase charge of the parameter of heat transfer, and perfusion rate [14]. He discovered that the lap durations for living tissues are identical. Dutta and Kundu conducted a study on thermal wave propagation for a changing and constant thermal flow through a skin tissue surface to detect hyperthermia [15].

To analyze the increment distribution of deep skin tissues temperature in a human leg, Agrawal and Pardasani suggested utilizing a finite element approach [16]. Kumar and Rai investigated the thermal behaviour using the dual-phase-lag (DPL) model based on time-fractional consideration [17]. Moradi et al. used the dual-phase-lag (DPL) bioheat transfer model with an analytical approach on skin tissue as a finite domain of cosine, continuous, and pulse thermal flow conditions on the skin’s bounding surface [8].

Liu and Xu improved a different analytical approach for Pennes’s bioheat equation solution for fluctuations in the temperature of skin tissue due to the sinusoidal thermal flow [18]. Shih et al. overcame the impacts of studying the temperature reaction on the biological skin tissue under sine heat flux [19]. The exact analytical solution of Pennes’s bioheat and hyperbolized bioheat transfer models for periodic constant and pulsed train heat flux boundary conditions has been provided by Ahmadikia et al. [20]. To characterize the impacts of the pulsatile blood flow in thermal disruption, Horng et al. used a statistical model [21]. He observed that the thermal region of the tumor had a pure little variation on the continuum range of blood from usual to parabolic blood flow velocity profiles. Shih et al. investigated how the thermal relaxation time interacts with the pulsatile blood flow during thermal treatments in living tissues [22]. Youssef and Alghamdi obtained the solution to one-dimensional problems of thermoelastic dual-phase-lag (DPL) skin tissue based on specific temperature stress [23]. Kundu and Dewanjee pioneered the response of a non-Fourier thermal approach in one layer of skin tissue [24]. Ezzat and Alabdulhadi introduced a novel mathematical model of a generalized thermoviscoelasticity theory based on Pennes’s bioheat transfer equation with a dual phase-lag to investigate the biothermomechanics behaviour in the skin tissue [25]. Ezzat constructed the model of the thermoviscoelasticity theory of fractional dual-phase-lag heat conduction law with rheological properties of the volume to investigate one-dimensional bioheat transfer and heat-induced mechanical response in human skin tissue [26]. In Pennes’s bioheat transmission and heat-induced mechanical response in human living tissue with variable thermal conductivity based on rheological parameters of the volume, Ezzat presented a novel notion of memory-dependent derivative [27]. Many applications of analytical analysis of the dual-phase-lag model of bioheat transfer have been constructed and solved [2832].

In the current work, the exact solutions of a novel mathematical model of skin tissue will be obtained in the framework of the dual-phase-lag (DPL) heat conduction model when the bounding plane of the surface of the skin tissue is exposed to a continuous heat flux due to a constant voltage of the thermoelectric effect. To obtain the exact analytical solutions, variables will be separated into a finite domain. The influence of the dual-phase-lag time parameters, voltage value, resistance value, and electric shock time value will be investigated and discussed. Comparing the three different models of bioheat conduction is the main goal of this work and is a novel effort.

2. Basic Equations

The model of bioheat transmission was created to investigate the time-dependent temperature increment caused by a heat source or thermal loading. Pennes created the first biological tissue model in the context of the classical Fourier’s law of heat conduction as follows [23, 3336]:where , , , and are the blood density, blood perfusion rate, specific heat of the blood, and blood temperature, respectively. , , and are the thermal conductivity, density, and specific heat of the skin tissue. denotes the absolute temperature function. is the metabolic heat generated by the chemical reaction inside the skin tissue, and it is assumed to be a constant, while is the external heat source and it could be a function. Finally, is the well-known Laplace operator.

Vernotte and Cattaneo (V-C) have updated the classic Fourier law of thermal transfer by positioning the assumption of the finite and limited speed of the propagation of the thermal wave and utilizing the following hyperbolic heat conduction form [23, 3336]:where is defined as the relaxation time parameter, gives the thermal wave’s speed inside the medium, and gives the thermal diffusivity.

The dual-phase-lag (DPL) heat conduction equation is based on the dual response between the temperature gradient and heat flow, which modifies the heat conduction equation of Tzou’s classical Fourier’s.

Then, in this case, the heat conduction equation takes the following form [23, 3336]:where is the second relaxation time parameter, which gives the phase lag time of the temperature gradient.

We assume the temperature increment’s function to take the following form [35]:

Then, we obtain the heat conduction form as follows:

3. Problem Formulation

We assume that the region of a skin tissue obeys the dual-phase-lag formulation as in equation (5) (see Figure 1).

The medium is initially quiescent and has no external heat source i.e., , while is a constant.

We consider the surface of the skin tissue to be subjected to constant heat flux , while the surface of the other side has a zero-heat flux.

Thus, the heat conduction model takes the following form [35]:

The initial conditions take the following form:

The boundary conditions take the following form:

The boundary value problem (B. V. P.) in equations (6)–(8) contains nonhomogeneous partial differential equations with nonhomogeneous boundary conditions on the surface of the skin tissue. Hence, the differential equations must be formulated in two parts, i.e., a steady part and a transient part as follows [35, 37, 38]:

Then, the transient part takes the following form:where the second part of the steady state is as follows:where .

The steady-state part has the following boundary conditions [35]:

Then, according to the boundary conditions of equation (8), the solution of equation (6) is as follows:

The initial conditions and the boundary conditions of the transient part of the differential equation (7) take the following form:and

To get the solution of equation (10), we write the expansion of the function in the Fourier series expansion as follows [35, 37, 38]:where it must satisfy the boundary conditions in equation (15).

By substituting equation (16) into equation (10), we obtain the following equation:where , and .

The solutions of equation (17) in general forms take the following form [35, 37, 38]:and from this, we havewhere are the solutions of the characteristic equation represented as

The roots of the characteristic equation (20) take the following form:where .

To use the initial condition of equation (14), we must expand in a Fourier series expansion as follows:

The initial condition of equation (14) lead to the following linear system of algebraic equations:and

In the case of , we have

Thus, when , we have

Then, equations (23) and (24) introduce the following systems:and

When , we have

Then, the system of equations (27) and (28) can be reduced to a system of equations as follows:and

The equations (30) and (31) give two systems of algebraic equations as follows:

The solution of the abovementioned two systems gives the following equation:

In this case, the exact solution of equation (6) takes the following form:

For this case , so, we have

By using the initial condition of equation (14), we get

Solving the abovementioned system gives

Thus, in this case, the exact solution of equation (6) is [35]

Now, we will consider that the skin tissue has been considered to be subjected to a thermal shock for a small value of time due to the thermal effect of an electrical shock with a constant voltage V (V) and a constant resistance of the skin tissue [36].

According to Joule’s equation of electrical heating, we have [36]where is the time interval parameter of the electrical shock.

4. Results and Discussions

The temperature distribution of skin tissue is examined by using three bioheat transfer models (V-C, Pennes, and Tzou). Table 1 shows the values of the pertinent thermal parameters that were used in the numerical results and calculations for this project [17, 23, 24, 3032, 35, 39].

In Figures 26, we plot the solutions of the heat conduction temperature given in equation (38) versus the length. All figures from 26 seem to have the same behaviour, which is that the temperature decreases rapidly as the length increases. Specifically, in Figure 2, we compare the temperature increment distribution for the voltage values It is clear from the figure that as the value of the voltage increases, the temperature will rise when the length is less than 0.008. However, there is no big difference shown in temperature after this point.

We also studied the influence of resistance when in the temperature increment distribution, as shown in Figure 3. It clearly shows that the high resistance plays a significant role in reducing the temperature when the length is less than 0.008. This result seems to agree with what is expected intuitively.

Figure 4 shows the temperature increment distribution for the two different values of electric shock time, . We noted that as the electric shock time increases, the temperature increment distribution will also rise which intuitively makes sense.

We can compare this figure with Figure 8 in [26], where both figures represent the same result and the same attitude of the temperature increment based on the different values of the ramp-time heat parameter.

In Figure 5, we compare the temperature increment distribution considering the three cases of relaxation time parameters, . When , the temperature declines slightly and uniformly over the x-axis whereas the temperature will start to increase somewhat if the value of exceeds or is equal to the value of . The values of the temperature increment take the following order:

Figure 6 compares the distribution of the temperature change for three different bioheat models, which are Pennes, Vernotte–Cattaneo, and Tzou. The Tzou model plays a considerable role in reducing the temperature. The values of the temperature increment based on the three studied models take the following order:

Equation (41) shows that the profile of the speed of the thermal wave through the skin tissue is decreasing from Pennes’s model, Vernotte–Cattaneo, to Tzou. In other words, Tzuo’s model ensures that the thermal wave transfer through the skin tissues at finite or limited speeds is better than other models, thus making this model the best and closest to the reality of the physical behaviour of the skin than other models. So, depending on the lag-time parameters of temperature gradient and heat flux makes the model of Tzou more successful and more useful than the other models.

Then, we studied the behaviour of temperature increment distribution over time to understand as to what extent can other factors such as voltage, resistance, electric shock time, and relaxation parameters influence the temperature, and this is shown in Figures 711.

Figures 711 show the deep fluctuation in temperature at a high altitude in the first few seconds. Then, it gradually stabilizes into a periodic motion as time increases.

In Figure 7, the distribution of the temperature increment is illustrated for two values of voltages, . It indicates that the voltage has a significant influence on the temperature increment distribution. Moreover, increasing resistance and shock time seem to have a non-neglected effect on the temperature increment distribution either by increasing or decreasing the temperature, respectively, and this is explained in Figures 8 and 9.

To study the influences of time relaxation parameters on the temperature, we plotted the temperature increment distribution over time considering three cases: It is clear from Figure 10 that the temperature fluctuates from the highest to lowest temperature over time when at a high altitude whereas the temperature zigzags fast at a high altitude at the beginning and then declines significantly over time as . Moreover, the attitude of the temperature increment in Figure 10 is the same as in Figures 3 and 7 in Ezzat [26].

Figure 11 compares the behaviour of temperature increment distribution over time between three different models: Pennes, Vernotte−Cattaneo, and Tzou. Among these models, Tzou seems to illustrate the lowest temperature increment with the smallest altitude.

For validation and to know which model is more close to the real behaviour of the skin tissue, Figure 12 has been performed with conditions similar to the experimental condition in [40] to compare the behaviour of temperature increment distribution over a long interval of time between the three different models: Pennes, Vernotte−Cattaneo, Tzou, and the experimental and cited results in [40]. Thus, this figure shows that Tzou is the best model where the temperature increment based on this model is closer to the experimental results more than the other models.

5. Conclusion

The temperature interaction and the response of the skin tissue caused by a continuous flow of surface heat caused by a constant voltage electrical current was the goal of this study. The bioheat transfer of the dual-phase-lag (DPL) model has a precise analytical solution. It is used to separate the variables into a finite domain for the governing equations. The reactions at the transition temperature have been measured and studied. Figures have been provided that compare the Pennes, Vernotte−Cattaneo, and Tzou models.

The results concluded the following:(i)When the value of the voltage increases, the temperature rises.(ii)When the value of the resistance decreases, the temperature rises.(iii)When the electric shock time increases, the temperature increment distribution also rises.(iv)When the value of the temperature relaxation time increases more than the value of the gradient temperature relaxation time , the temperature increment distribution rises.(v)The temperature increment based on Penne’s model is greater than its values based on the Vernotte–Cattaneo model, Tzou’s model.(vi)The current results show that the thermal waves based on Tzou spread at a finite speed in the skin tissue which removes the Vernotte–Cattaneo and Pennes’s shortcomings.(vii)Depending on the lag-time parameters of temperature gradient and heat flux makes the model of Tzou more successful and more useful than the other models.

Data Availability

No data were used to support the study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.