A semianalytical method based on the separation of variables is presented to solve the problem of heat transfer in multilayer skin tissues. This method is used to predict the temperature field and thermal damage of skin tissues, which are subjected to the time-varying laser heating and fluid cooling. The effect of convective heat transfer coefficient of the skin surface on temperature field and thermal damage is studied. The temperature field and thermal damage of skin tissue with different laser heating cycles are calculated. The relationship between the laser heating period and the temperature fluctuation is obtained. The effects of the peak time of laser heating on temperature and the thermal damage of skin tissues are investigated. The relationship between the peak time of laser and the temperature fluctuation of skin tissue is drawn. The results calculated in this paper can provide some guidance for hyperthermia of skin tissue. The presented semianalytical method can be applied to solve the general heat transfer problems of multilayer structures.
Department of Education of Guangdong Province2017KTSCX2181. Introduction
Thermotherapy has been widely used in the treatment of skin tissue diseases. Laser, microwave, and radio frequency are often used to irradiate the surface of skin tissue during thermotherapy [1–3]. Accurate prediction of skin tissue temperature and thermal damage during hyperthermia plays an important role in thermotherapy [4]. More and more researchers have studied the heat transfer and thermal damage of the skin tissue. Important research can be summarized as follows.
Deng and Liu gave several closed-form analytical solutions to the bioheat transfer problems which are based on Green’s function method [5]. Xu et al. studied the heat transfer of skin tissue by using pence model; they established a mathematical model of heat transfer of three-layer skin tissue and solved the model by using finitely different method [6]. Liu et al. established DPL model for heat transfer of multilayer skin tissue; the Laplace transform method was used to solve this model [7]. Tian et al. established a skin-fabric five-layer heat transfer model; a semianalytical method was employed to solve this model [8]. Ahmadikia et al. used analytical methods based on Laplace transform to study the distribution of temperature field in skin tissue under constant or transient heat flux conditions; parabolic and hyperbolic skin tissue heat transfer equations were studied [9]. In the study of Kumar et al., a dual-phase-lag model of bioheat transfer had been studied by using Gaussian distribution source term with the most generalized boundary condition during hyperthermia treatment [10]. Lin solved the analytical solution of heat transfer of single skin layer which was based on the separation of variable methods; Fourier and non-Fourier cases were studied [11]. Hooshmand et al. developed an analytical solution method for solving the problem of heat transfer in biological tissue, generalized dual-phase-lagging (DPL) model was studied, and the calculation results were compared with the classical DPL model [12]. Kumar et al. dealt with the study of heat transfer and thermal damage in triple layer skin tissue using fractional bioheat model. An implicit finite difference scheme is obtained by a new method [13]. Kumar and Rai studied the thermal characteristics in skin tissue using time fractional dual-phase-lag bioheat transfer (DPLBHT) model subjected to Dirichlet boundary condition. They solved this problem using the finite element Legendre wavelet Galerkin method. They treated skin tissue as a single layer [14]. Li et al. predicted the temperature of different biological tissues using Fourier and non-Fourier heat transfer models. And the temperature response of different tissues was measured. The experimental results were compared with those predicted by the model [15]. Hobiny and Abbas used the Laplace transform method to solve the problem of biological tissues subjected to moving heat source [16]. Ghanmi and Abbas derived an analytical solution for the fractional transient heating within the skin tissue during hyperthermia [17]. Sur et al. used Laplace transform to study the thermal behavior of skin tissue when it was subjected to moving heat sources [18]. Sharma and Kumar established a nonlinear dual-phase-lag heat transfer model of biological tissue. The finite element Runge–Kutta method was used to solve this model. The temperature field and thermal damage of skin tissue were studied using the numerical solution of the model [19].
Skin tissue is the largest organ of the human body, including epidermis, dermis, and subcutaneous tissue [20]. Thermophysical parameters of different layers of skin tissue are different. So it is difficult to solve the heat transfer model of skin tissue by analytical method [9, 11, 21, 22]. Numerical methods were often used to solve this problem [23]. Several papers used analytical methods, but all of them treated skin tissue as a single layer [12]. In this paper, a semianalytical method is developed to solve heat transfer in multilayer skin tissues. A special method of separation of variables is used to derive the analytical solution. In the derivation process, the time domain is firstly divided into several small time periods. Then the temperature function at the interface of different layers of skin tissue is assumed to be linear in a small period of time. After that, the analytical solution of separated variables of skin tissue heat transfer is deduced in a small period of time. According to the condition that the temperature at the interface is equal and the heat flux at the interface is continuous, the temperature at the interface is calculated. Temperature field of skin tissue during laser heating and fluid cooling is obtained by the semianalytical solution. The effects of several parameters such as convective heat transfer coefficient, period, and peak time of laser heating on temperature field are studied. These parameters can be used to control the temperature of the skin tissue and the thermal damage in the appropriate range.
2. Mathematical Model
Human skin tissue consists of epidermis, dermis, and subcutaneous. The geometric model of skin tissue is shown in Figure 1. The corresponding mathematical model can be established by combining the boundary conditions. The governing equation can be expressed as (1) [2, 24–27]:
The boundary conditions can be expressed as follows:
At x = 0,(2)k1∂T0,t∂x=h0T−T∞+qt.
At x = L3,(3)TL3,t=Ta.
The initial condition is(4)Tx,0=Ta.
The corresponding dimensionless governing equation is(5)rki∂2θ∂ξ2−ω¯biθ−rρirci∂θ∂τ+q¯met,i+q¯ext,i=0,ξi−1<x<ξi,τ>0.
The dimensionless boundary conditions are as follows:
At ξ = 0,(6)∂θ0,τ∂ξ−Nuθ0,τ=−Nu0θ∞+f¯1,1τ=f1,1τ.
At ξ = 1,(7)θ1,τ=f3,2τ.
The dimensionless initial condition is(8)θξ,0=0.
3. Solution Method
If the value of Δτ is small enough, in time domain [0 Δτ], f1,1τ can be expressed as [28](9)f1,1τ=f1,1Δτ−f1,10Δττ=v1,1τ,(10)θξ1,τ=v1,2τ=f1,2τ=f2,1τ,(11)θξ2,τ=v2,3τ=f2,2τ=f3,1τ.
Equations (7) and (9)–(11) show that the boundary conditions on both sides of each layer of skin tissue are known in a short time interval [0 Δτ]. A extended shifting variable method is developed by Lee and Lin, which can be used to solve the temperature field of each layer of skin tissue [29].
Assume(12)θiξ,τ=viξ,τ+∑m=12gi,mξfi,mτ.
Here, gi,m (ξ), i = 1, 2, 3; m = 1, 2 are the shifting functions to be specified and vi (ξ, τ) is the transformed function. Substituting (12) into (5), the two following subsystems are obtained. The first subsystem is expressed in terms of the transforming variables viξ,τ as follows; the transform governing equation is(13)rki∂2vi∂ξ2−ω¯bvi−rρirci∂vi∂τ+q¯ts,i=0,(14)q¯ts,i=−ω¯b∑m=12gi,mf1,m−rρirci∑m=12gi,mdfi,mdτ+rki∑m=12d2gi,mdξ2fi,mτ+q¯met,i+q¯ext,i.
If i = 1, substitute (12) into (6) and (10):(15)∂v10,τ∂ξ−Nuvi0,τ+d∑m=12gi,m0fi,mτdξ−Nu∑m=12gi,m0fi,mτ=f1,1τ,(16)v1ξ1,τ+∑m=12g1,mξ1f1,mτ=f1,2τ.
In order to homogenize the boundary conditions when separating variables, (17) and (18) need to be satisfied:(17)∂v10,τ∂ξ−Nuvi0,τ=0,(18)v1ξ1,τ=0.
The characteristic equation of (13) is(29)rki∂2vi∂ξ2−rρirci∂vi∂τ=0.
Assume(30)vi=ψiξϕiτ.
Substitute (30) into (29):(31)d2ψiξdξ2−λψiξ=0,(32)dϕiτdτ−λrkirρirciϕiτ=0.
If i = 1, substitute equation (30) into (17) and (18):(33)rk1dψ10dξ−Nuψ110=0,(34)ψ1ξ1=0.
Equation (31) is a second-order homogeneous ordinary differential equation. Its characteristic function is as follows:(35)ψ1,jξ=cosλjξ1ξ+Nu0rk1λjsinλjξ1ξ,where λj is the positive solution of the following equation:(36)rk1λNu+tanλξ1=0.
If i = 2, 3, substitute (30) into (24) and (25):(37)ψi,jξ=cosλjξsinλjξi−1−cosλjξi−1sinλjξ,where λj is the positive solution of the following equation:(38)cosλsinλξi−1−cosλξi−1sinλ=0.
Obviously, ψi,jξ satisfies the following equation:(39)∫ξi−1ξiψi,jξψi,kξdξ=0,for j≠k;δi,j,for j=k.
Assume(40)viξ,τ=∑j=1∞ψi,jξγi,jτ.
Substitute (40) into (13) and integrate it in [ξi−1,ξi]:(41)dγi,jdτ+rkiλj+ηirρirciγi,j=qts,i,(42)qts,i=1rρirci1δi,j∫ξi−1ξiψi,jq¯ts,idξ,(43)ηi=Nui=1;ω¯bii=2,3.
The solution of (41) is(44)γi,jτ=e−αiτγi,j0+e−αiτ∫0τeα1τqtsidτ,(45)αi=rkiλj+ηirρirci,(46)θiξ,τ=∑j=1∞λi,jτψi,iξ+∑m=12gi,mξfi,mτ.
4. Thermal Damage
In addition to predicting the temperature field of skin tissue during thermotherapy, it is often necessary to predict its thermal damage. Formula (47) is developed by Henriques and Moritz which can be used to calculate thermal damage of skin tissue [1]:(47)Ωt=∫0tAexp−EaRTdt,where Ω represents the thermal damage factor of skin tissue. Ea represents the activation energy. T represents the temperature of the skin tissue. R represents the universal gas constant. According to the thermal damage factors, the thermal damage is divided into three levels in the literatures. Ω=0.53, the first level of thermal damage; Ω=1, the second level of thermal damage; Ω=104, the third level of thermal damage.
5. Numerical Results
When the surface of skin tissue is irradiated by pulsed laser, the variation of laser energy with time can be shown in Figure 2. Assume the initial temperature of the skin tissue T (x, 0) = 37°C. At x = L3, T (L3, 0) = 37°C. The thermophysical parameters of skin tissue can be shown in Table 1 [4].
Schematic diagram of pulsed laser heating and fluid cooling at skin tissue surface.
Thermophysical parameters of skin tissue.
Parameters
i = 1
i = 2
i = 3
References
Li (mm)
0. 1
1.6
6
[30]
ρi (kg/m3)
1190
1116
971
[31]
ki (W/m-K)
0.235
0.445
0.185
[26]
ci (J/kg-K)
3600
3300
2700
[26]
ϖb,i (ml/ml/s)
0
0.1
0
[32]
qmet,i (W/m3)
368.1
368.1
368.3
[9]
qext,i (W/m3)
0
0
0
[9]
Figure 3 shows the temperature history at the skin surface, ED interface, and DS interface in a laser heating cycle. It is found that the temperature rises first and then decreases with time at x = 0, x = 0.1 mm. If h = 7 W/m2°C, the temperature peaks appear at t = 2 s and t = 2.4 s, respectively. If h = 700 W/m2°C, the temperature peaks appear at t = 1.4 s and t = 1.7 s, respectively. So we can speculate that the larger the value of x, the shorter the time required for temperature peak. If h = 7 W/m2°C, at x = 1.6 mm, the temperature first remains constant and then continues to rise. When t = 10 s, the peak value of temperature still does not appear. If h = 700 W/m2°C, at x = 1.6 mm, because this location is far away from the heating source, temperature varies a little with time. It can be inferred from the comparison of solid and dotted lines that the larger the value of h, the lower the temperature. Investigating its reason we can know that the larger the value of h is, the more heat is taken away by water cooling on the surface of skin tissue, so the lower the temperature of skin tissue is.
Temperature history at skin surface, ED interface, and DS interface (TC = 10 s, tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C).
Figure 4 shows the temperature distribution in the x direction at different times. As shown in the picture, if h = 7 W/m2°C, t = 1 s, the temperature of the skin tissue surface is the highest, which is 65°C. The temperature distribution of skin tissue decreases in the x-axis. Comparing the temperature distribution curves of t = 0.5 s and t = 1 s, it is found that the temperature of T = 1 s is higher than that of t = 1 s. This is because the heat flux of laser heating keeps rising in 0 to 1 s. Because the heat flux of laser heating begins to decrease after t = 1 s, the temperature distribution at t = 5 s and 10 s is much lower than that at t = 1 s. Comparing the temperature distribution of h = 7 W/m2°C and h = 700 W/m2°C, it is found that, at the same time, the temperature distribution of skin tissue at h = 700 W/m2°C is significantly lower than that at h = 7 W/m2°C. Because the cooling effect is greater than the heating effect, the temperature of skin tissue surface is lower than 37°C at t = 5 s and t = 10 s.
Temperature distribution along x direction at different time (TC = 10 s, tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C).
Figure 5 demonstrates the relationship between heat flux and heating time at different locations of skin tissue. As can be seen from the figure, when h = 7 W/m2°C, the heat flux density on the surface of skin tissue (x = 0) increases first and then decreases, which is in accordance with the trend of surface heating heat flux. At the skin-dermis interface (x = 0.1 mm), the heat flux is slightly lower than that at the skin surface. At the epidermal-subcutaneous interface (x = 1.6 mm), the heat flux increases first and then decreases, but the range of variation is very small. On the left side of the skin tissue (x = 6 mm), the heat flux is 0 and remains unchanged. When h = 700 W/m2°C, the heat flux of heating is less than that of cooling in the initial period of time. The heat flux of skin surface (x = 0) and dermal-epidermal interface (x = 0.1 mm) changed from negative to positive, continued to increase, and then gradually decreased. Because the boundary condition we set is that the right side of the skin tissue is adiabatic, on the left side of the skin tissue (x = 6 mm), the heat flux is 0 and remains unchanged too.
Relationship between heat flux and heating time at different locations of skin tissue (TC = 10 s, tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C).
Figure 6 demonstrates the relationship between thermal damage and heating time at different locations of skin tissue. As can be seen from the figure, thermal damage of skin tissue increases with time accumulation. The surface of skin tissue has the greatest thermal damage. If h = 7 W/m2°C, the thermal damage value of skin tissue surface is close to 100, reaching the second-degree thermal damage, and the thermal damage value at the junction of epidermis and dermis exceeds 0.53, reaching the first-degree thermal damage. Because the greater the value of h, the better the heat dissipation and the smaller the thermal damage of skin tissue, the thermal damage at h = 700 W/m2°C is six orders lower than that at h = 7 W/m2°C. It can be seen from the figure that although the heat flux of laser heating increases first and then decreases, the heat damage increases continuously. As the heat flux of laser heating decreases, the trend of increasing skin tissue thermal damage slows down. In order to analyze the effect of R and H on thermal damage, a variance analysis was carried out. Table 2 shows the calculated results. By calculation, R2 (x) = 0.163812, R2 (h) = 0.318886, R2(interaction) = 0.318877. Therefore, it can be inferred that the influence level of x on thermal damage is lower than the average and the influence level of h on thermal damage is above average. These two factors and their interaction have a great influence on thermal damage.
Relationship between thermal damage and heating time at different locations of skin tissue (TC = 10 s, tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C).
Analysis of variance about x and h (α = 0.01).
Source of difference
SS
df
MS
F
P value
F crit
x
109139.8
1
109139.8
490.3853
1.13E − 79
6.677735
h
212457.6
2
106228.8
477.3055
2.5E − 124
4.641058
Interaction
212452
2
106226
477.2929
2.5E − 124
4.641058
Internal
132200.2
594
222.5593
—
—
—
Total
666249.6
599
—
—
—
—
Figure 7 demonstrates the effect of heating cycle Tc on the temperature history at skin surface. As can be seen from the figure, because heating and cooling are the same in the first five seconds, so the temperature curve coincides completely. If TC = 5 s, the first pulse of laser heating is completed after 5 seconds, the heat flux of heating no longer decreases, and the second laser heating pulse starts. Therefore, the temperature on the surface of skin tissue starts to rise when t = 5 s, and the rising process is fluctuant, and the fluctuating period is TC = 5 s. After five seconds, the shorter the heating cycle is, the faster the temperature rises. Temperature fluctuates with heating period. After 30 seconds, the amplitude of temperature fluctuation tends to be balanced. The shorter the heating period TC is, the greater the temperature fluctuation amplitude Ts,max is.
Effect of heating cycle TC on the temperature history at skin surface (tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 0).
Figure 8 shows the effect of heating cycle TC on the temperature history at the interface between epidermis and dermis. Comparing with Figure 7, it can be seen that the temperature fluctuation at the interface of epidermis and dermis is smaller than that at the skin surface. It is the same as the result shown in Figure 6, The shorter the heating period TC is, the greater the temperature fluctuation amplitude Ts,max is. When TC = 5 s, the amplitude of fluctuation reaches the maximum value; at the junction of epidermis and dermis, Ted1,max = 82.5°C, it is 9.5°C lower than the surface temperature of skin tissue. When TC = 15 s, the amplitude of fluctuation reaches minimum value; at the junction of epidermis and dermis, Ted3,max = 65.2°C, it is 7.9°C lower than the surface temperature of skin tissue.
Effect of heating cycle TC on the temperature history at the interface between epidermis and dermis (tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 0.1 mm).
Figure 9 shows the effect of heating cycle TC on the temperature history at the interface between dermis and subcutaneous. Because the location is far from the heat source, the range of temperature fluctuation is much smaller than the skin surface and the interface between epidermis and dermis. It takes a long time for temperature fluctuation to reach a steady state. As can be seen from the figure, when TC values are different, the time required for temperature fluctuation to reach stability is 75 s, 100 s, and 150 s, respectively. It can be inferred that the greater the value of temperature TC is, the less time it takes for temperature fluctuation to reach stability. It can also be seen from the figure that the larger TC value is, the larger temperature fluctuation range is.
Effect of heating cycle TC on the temperature history at the interface dermis and subcutaneous (tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 1.6 mm).
Figure 10 shows the effect of heating cycles on the thermal damage at different locations of skin tissue. Because the surface of skin tissue is close to the heat source, the thermal damage on the surface of skin tissue (x = 0) is much greater than the interface between dermis and epidermis (x = 0.1). With the accumulation of heating time, the thermal damage of skin surface in all three cases exceeds 104, reaching the third-degree burn. The thermal damage with TC = 5 s is four orders higher than the thermal damage with TC = 10 s. It can be inferred from the figure that the smaller the heating period is, the more serious the thermal damage of skin tissue is. The thermal damage at the interface between epidermis (x = 0.1 mm) and dermis exceeds 104, reaching the third-degree burn when TC = 5 s or TC = 10 s. When TC = 15 s, the thermal damage value at the interface between epidermis and dermis exceeds 100 to reach second-degree burn. Because the dermis-subcutaneous interface (x = 1.6 mm) is far away from the heat source, no burn occurs with TC = 10 s and TC = 15 s, and the first-degree burn occurs with TC = 5 s. In order to study the influence of TC and x on thermal damage, variance analysis was carried out. Table 3 shows the calculation results. By calculation, R2 (x) = 0.165575, R2 (TC) = 0.167528, and R2(Interaction) = 0.3308943. Therefore, it can be inferred that the influence level of x and h on thermal damage is lower than the average, respectively. These two factors and their interaction have a significant effect on thermal damage.
Effect of heating cycle TC on the thermal damage history at different locations of skin tissue (tp = 1 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 1.6 mm).
Analysis of variance about TC and x (α = 0.01).
Source of difference
SS
df
MS
F
P value
F crit
x
2.04E + 25
2
1.02E + 25
13309.32
0
4.605563007
h
2.06E + 25
2
1.03E + 25
13466.32
0
4.605563007
Interaction
4.07E + 25
4
1.02E + 25
13305.7
0
3.319522681
Internal
4.13E + 25
53991
7.65E + 20
—
—
—
Total
1.23E + 26
53999
—
—
—
—
Figure 11 indicates the relationship between laser heating cycle and temperature fluctuation amplitude. Ts,max and Ts,min represent the maximum and minimum values of temperature fluctuation amplitude of skin tissue surface, respectively. Ted,max and Ted,min represent the maximum and minimum values of temperature fluctuation amplitude of the interface between epidermis and dermis, respectively. Tds,max and Ted,min represent the maximum and minimum values of temperature fluctuation amplitude at the interface between dermis and subcutaneous, respectively. It can be seen from the figure that the amplitude of temperature fluctuation increases first and then decreases with the increase of laser heating period. The range of temperature fluctuation increases with the increase of laser heating period. Skin tissue surface temperature fluctuation range is the largest, the epidermis and dermis interface is the second, and the dermal and subcutaneous interface is the smallest. With the increase of laser period, the minimum amplitude of temperature fluctuation becomes closer and closer; when TC > 16 s, the lowest temperature of skin tissue is equal at all locations.
Relationship between laser heating period and temperature fluctuation amplitude (tp = 2 s, qpeak = 3 W/cm2, T∞ = 20°C).
Figure 12 shows the effect of time tp on skin surface temperature history. It can be seen from the figure that the skin tissue surface temperature fluctuates and rises at the beginning time and then tends to fluctuate steadily. If tp = 3 s, the fluctuation range of skin surface temperature is 80–107°C. If tp = 1 s, the fluctuation range of skin surface temperature is 45–75°C. If tp = 0.5 s, the fluctuation range of skin surface temperature is 42–63°C. We can infer that the larger the tp value is, the larger the maximum value of skin tissue surface temperature fluctuation is and the smaller the minimum value of skin tissue surface temperature fluctuation is.
Effect of peak time tp on the temperature history at skin tissue surface (TC = 10 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 0.1 mm).
Figure 13 shows the effect of peak time tp on the temperature history at the interface between epidermis and dermis. Comparing with Figure 12, we can see that the temperature of this interface is lower than that of the skin surface. The trend of fluctuation is consistent with Figure 12.
Effect of peak time tp on the temperature history at the interface between epidermis and dermis (TC = 10 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 0.1 mm).
Figure 14 shows the effect of peak time tp on the temperature history at the interface between dermis and subcutaneous. Because this location is far away from the heat source, the temperature of this location is much lower than that of the skin tissue surface and the interface between the epidermis and dermis. The temperature at this position fluctuates with the heat flux heated by laser.
Effect of peak time tp on the temperature history at the interface between dermis and subcutaneous (TC = 10 s, qpeak = 3 W/cm2, T∞ = 20°C, x = 0.1 mm).
Figure 15 indicates the relation between peak time tp and temperature fluctuation Amplitude. The meaning of Ts,max, Ts,min, Ted,max, and Ted,min is the same as shown above. It can be seen from the figure that, with the increase of peak time tp, the amplitude of temperature fluctuation becomes larger and larger at all positions of skin tissue. The temperature fluctuation range of skin surface and epidermis-dermis interface increases with the increase of peak time tp. The temperature fluctuation range of dermal-subcutaneous interface increases first and then remains constant with the increase of tp. The larger the tp value is, the smaller the range of temperature fluctuation is.
Relation between peak time tp and temperature fluctuation amplitude (TC = 2 s, qpeak = 3 W/cm2, T∞ = 20°C).
Figure 16 indicates the effect of peak time tp on the thermal damage history at different locations of skin tissue. It can be seen from the figure that the greater the tp value, the more severe the thermal damage of the skin tissue. With tp = 3 s thermal damage value of skin surface and epidermal-dermal interface exceeds 104, reaching the third-degree burn, and thermal damage value of dermal-subcutaneous interface exceeds 1 to reach a second-degree burn. If tp = 0.5, the thermal damage on the skin surface (x = 0) exceeds 100 to reach a second-degree burn. The thermal damage on the skin surface (x = 0.1 mm) exceeds 1 to reach a first-degree burn. Because the dermal-subcutaneous interface is far away from the laser heating heat source, there is no burn. If tp = 1, the thermal damage on the surface of skin tissue exceeds 106, reaching the third-degree burn. The thermal damage on the epidermis-dermis interface of skin tissue exceeds 106, reaching the second-degree burn. The thermal damage on the dermal-subcutaneous interface of skin tissue is less than 0.53, lower than the first-level burn standard. In order to study the influence of tp and x on thermal damage, variance analysis was carried out. Table 4 shows the calculation results. By calculation, R2 (x) = 0.165721, R2 (tp) = 0.167244, and R2(Interaction) = 0.3325. Therefore, it can be inferred that the influence level of x and tp on thermal damage is lower than the average, respectively. These two factors and their interaction have a significant effect on thermal damage.
Effect of peak time tp on the thermal damage history at different locations of skin tissue (TC = 10 s, qpeak = 3 W/cm2, T∞ = 20°C).
Analysis of variance about tp and x (α = 0.01).
Source of difference
SS
df
MS
F
P value
F crit
x
2E + 32
2
9.98E + 31
13329.93
0
4.605563
tp
2.01E + 32
2
1.01E + 32
13450.85
0
4.605563
Interaction
3.99E + 32
4
9.98E + 31
13329.93
0
3.319523
Internal
4.04E + 32
53991
7.48E + 27
Total
1.2E + 33
53999
—
—
—
—
6. Discussion
It is of great significance to accurately predict the temperature field and control the thermal damage during skin tissue hyperthermia. In this paper, a special separation variable method is used to derive the semianalytical solution of temperature field in skin tissue hyperthermia. Compared with the traditional finite element method, finite difference method, Laplace transform method, and so on, this method has clearer physical meaning and higher calculation efficiency. The temperature field and thermal damage of skin tissue during laser hyperthermia are predicted by the analytical solution derived in this paper. Several parameters that can affect the temperature field and thermal damage of skin tissue during heating have been studied. Through analysis we can find that the convective heat transfer coefficient between the skin tissue surface and the cooling fluid has a greater effect on the temperature field and thermal damage. We can change the injection speed of the cooling fluid to change the convection heat transfer coefficient and then control the temperature and thermal damage of the skin tissue in an appropriate range. In addition, the heating period and peak time of the laser also have a greater impact on the skin tissue temperature field and thermal damage during hyperthermia. Adjusting the parameters such as laser heating period or peak time can control the skin tissue temperature to fluctuate within a stable range. This operation can achieve the desired effect of hyperthermia without thermal damage to skin tissue.
7. Conclusion
In this paper, a semianalytical method is presented to solve the heat transfer problem of multilayer skin tissues. The temperature field and thermal damage of skin tissue under time-varying laser heating and fluid cooling are obtained by the presented method. The effects of several parameters on skin tissue temperature and thermal damage are investigated. The main phenomena are revealed as follows:
The smaller the convective heat transfer coefficient h is, the faster the temperature of skin tissue rises
The larger the convective heat transfer coefficient h is, the smaller the heat flux of skin tissue is
The smaller the heating period is, the more serious the thermal damage of skin tissue is
The range of temperature fluctuation decreases with the increase of laser heating period TC
The greater the tp value, the more severe the thermal damage of the skin tissue
In this paper, the Fourier model of heat transfer in biological tissue is used to study the heat transfer in skin tissue. In recent years, some literatures have used non-Fourier model to simulate the heat transfer in biological tissue. Our future research direction is to use the method of this paper to solve the non-Fourier heat transfer model of biological tissue and then to study the heat transfer of skin tissue.
Nomenclaturei:
Layer number of skin tissue: i = 1, epidermis; i = 2, dermis; i = 3, subcutaneous
ki:
Thermal conductivity of i-th layer
T:
Temperature
ρi:
Density of i-th layer
ci:
Specific heat capacity of i-th layer
t:
Time
δi:
Thickness of i-th layer
Li:
x-axis coordinates, Li=∑i=1iδi
ϖb:
Blood perfusion rate
ρb:
Density of epidermis blood
cb:
Specific heat capacity of blood
qmet,i:
Metabolic heat generation in the skin of i-th layer
qext,i:
Heat source due to external heating i-th layer
h0:
Heat convective coefficient
T∞:
Temperature of cooling fluid
Ta:
Temperatures of arterial blood
qt:
Heat flux function
rki:
Dimensionless kirki=ki/k1
θ:
Dimensionless Tθ=T−Ta/Ta
θ∞:
Dimensionless θ∞=T∞−Ta/Ta
ξ:
Dimensionless ξ=x/L3
ω¯bi:
Dimensionless ϖbω¯b=L2ϖbρbcb/k1
rρi:
Dimensionless rρirρi=ρi/ρ1
rci:
Dimensionless rci=ci/c1
τ:
Dimensionless tτ=k1t/ρcL2
q¯met,i:
Dimensionless q¯met,iq¯met,i=qmet,iL32/k1Ta
q¯ext,i:
Dimensionless q¯ext,iq¯ext,i=qext,iL32/k1Ta
Nu0:
Dimensionless Nu0=h0L3/k1
f¯1,1τ:
Dimensionless f¯1,1=L3qt/k1Ta.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are grateful to the financial body, Department of Education of Guangdong Province (associated grant number: 2017KTSCX218).
HenriquesF. H.MoritzA. R.Studies of thermal injury: I. The conduction of heat to and through skin and the temperatures attained therein. A theoretical and an experimental investigation1947234530549PennesH. H.Analysis of tissue and arterial blood temperatures in the resting human forearm1948129312210.1152/jappl.1948.1.2.93MaJ.SunY.YangJ.Analytical solution of non-Fourier heat conduction in a square plate subjected to a moving laser pulse201711560661010.1016/j.ijheatmasstransfer.2017.07.0492-s2.0-85026262165LiuK.-C.LinC.-N.WangJ.-S.Numerical solutions for the hyperbolic heat conduction problems in a layered solid cylinder with radiation surface2005164380582010.1016/j.amc.2004.06.0532-s2.0-16344396127DengZ.-S.LiuJ.Analytical study on bioheat transfer problems with spatial or transient heating on skin surface or inside biological bodies2002124663864910.1115/1.15168102-s2.0-0036960340XuF.LuT. J.SeffenK. A.Biothermomechanics of skin tissues20085651852188410.1016/j.jmps.2007.11.0112-s2.0-41749114034LiuK.-C.WangY.-N.ChenY.-S.Investigation on the bio-heat transfer with the dual-phase-lag effect201258293510.1016/j.ijthermalsci.2012.02.0262-s2.0-84860716751TianM.ZhuS.PanN.Skin thermal stimulation on touching cool fabric from the transient stage to steady-state stage201253808810.1016/j.ijthermalsci.2011.10.0082-s2.0-84855764886AhmadikiaH.FazlaliR.MoradiA.Analytical solution of the parabolic and hyperbolic heat transfer equations with constant and transient heat flux conditions on skin tissue201239112113010.1016/j.icheatmasstransfer.2011.09.0162-s2.0-84855347055KumarP.KumarD.RaiK. N.A numerical study on dual-phase-lag model of bio-heat transfer during hyperthermia treatment201549-509810510.1016/j.jtherbio.2015.02.0082-s2.0-84923289665LinS.-M.Analytical solutions of bio-heat conduction on skin in fourier and non-fourier models2013134135006310.1142/s02195194135006322-s2.0-84880173384HooshmandP.MoradiA.KhezryB.Bioheat transfer analysis of biological tissues induced by laser irradiation20159021422310.1016/j.ijthermalsci.2014.12.0042-s2.0-84922771837KumarS.DamorR. S.ShuklaA. K.Numerical study on thermal therapy of triple layer skin tissue using fractional bioheat model2018114185005210.1142/s17935245185005232-s2.0-85045474799KumarD.RaiK. N.Numerical simulation of time fractional dual-phase-lag model of heat transfer within skin tissue during thermal therapy201767495810.1016/j.jtherbio.2017.05.0012-s2.0-85019075647LiC.MiaoJ.YangK.Fourier and non-Fourier bio-heat transfer models to predict ex vivo temperature response to focused ultrasound heating20181231710.1063/1.50226222-s2.0-85046642669HobinyA. D.AbbasI. A.Theoretical analysis of thermal damages in skin tissue induced by intense moving heat source20181241011101410.1016/j.ijheatmasstransfer.2018.04.0182-s2.0-85045186290GhanmiA.AbbasI. A.An analytical study on the fractional transient heating within the skin tissue during the thermal therapy20198222923310.1016/j.jtherbio.2019.04.0032-s2.0-85064758401SurA.MondalS.KanoriaM.Influence of moving heat source on skin tissue in the context of two-temperature memory-dependent heat transport law2020431557110.1080/01495739.2019.16602882-s2.0-85073985071SharmaS. K.KumarD.A study on non-linear DPL model for describing heat transfer in skin tissue during hyperthermia treatment202022448110.3390/e22040481DowlingD. F.The thickness of cattle skin19556577610.1071/ar95507762-s2.0-3543060758LiuJ.ChenX.XuL. X.New thermal wave aspects on burn evaluation of skin subjected to instantaneous heating199946442042810.1109/10.7529392-s2.0-0032932502ZhangH.WangW.ZhangS.ZhaoZ.Semi-analytic solution of three-dimensional temperature distribution in multilayered materials based on explicit frequency response functions201811820822210.1016/j.ijheatmasstransfer.2017.10.1182-s2.0-85032820536AbbasI. A.MarinM.Analytical solution of thermoelastic interaction in a half-space by pulsed laser heating20178725426010.1016/j.physe.2016.10.0482-s2.0-85008411860CattaneoC.A form of heat conduction equation which eliminates the paradox of instantaneous propagation1958247431433VernotteP.Les paradoxes de la théorie continue de léquation de la chaleur195840DuckF. A.1990Cambridge, MA, USAAcademic PressVyasR.RustgiM. L.Green’s function solution to the tissue bioheat equation19921951319132410.1118/1.5967672-s2.0-0026785712PalaniG.AbbasI. A.Free convection MHD flow with thermal radiation from an impulsively-started vertical plate2009141738410.15388/na.2009.14.1.14531LeeS. Y.LinS. M.Dynamic analysis of nonuniform beams with time-dependent elastic boundary conditions199663247447810.1115/1.27888922-s2.0-0005081372BraznikovA.KarpychevV.LuikovaA.One engineering method of calculating heat conduction process197528KaminskiW.Hyperbolic heat conduction equation for materials with a nonhomogeneous inner structure1990112355556010.1115/1.29104222-s2.0-0025471339XuF.SeffenK. A.LuT. J.Non-Fourier analysis of skin biothermomechanics2008519-102237225910.1016/j.ijheatmasstransfer.2007.10.0242-s2.0-41649101655