Abstract

The physiological mechanisms conduction, convection, and radiation exchange the heat energy in bi-directional routes between the body and the temperature field. Metabolism and evaporation are the uni-directional routes for the exchange of heat energy. In the metabolic process, the body creates internal heat energy, whereas the body loses excess heat energy through the evaporation process and maintains the body temperature. This study has shown steady and unsteady state temperature distribution in three skin layers: epidermis, dermis, and subcutaneous tissue, during walking and marathon. The results have analyzed that each skin layer temperature is higher during a marathon compared with walking due to more metabolic effects. The computation has been carried out for the two-dimensional Pennes’ bio-heat equation using a finite element approach. The generated results have been exhibited graphically.

1. Thermal Regulation Mechanism

Human organisms maintain or control heat energy, primarily through the radiation, convection, metabolism, diffusion, and evaporation processes. Radiation, diffusion, and convection contribute to receiving heat energy when the ambient temperature is higher than skin temperature and losing heat energy when the ambient temperature is lower than skin temperature. Hypothalamus is the central unit and plays a significant role in thermoregulatory responses. It checks the current body’s core temperature and compares it with the natural temperature of about 37°C. If the body’s core temperature is too low, the body generates and maintains heat through the processes of shivering or metabolism to sustain a healthy life. If the current body core temperature is too high, the surplus heat is eliminated by the evaporation process, which cools the skin and maintains the body temperature. Metabolic heat generation, thermal environmental heat load, and heat loss capacity are the key components of the heat balance equation.

During walking, the body creates metabolic heat energy and dissipates it on the body’s surface. There is no storage heat energy. As a result, the body is thermally balanced. The heat gain and heat loss by the body is in equilibrium. Therefore, the heat exchange between the body and the temperature field is minimum, and the body core temperature maintains constant at 37°C. Since the body core temperature is proportional to the metabolic rate and largely independent of a wide range of environmental conditions. In a temperate or cool environment, the body loses heat energy through the convection process and reduces the skin’s blood flow rate. The physiological mechanisms are only incapable to maintain the body temperature. Therefore, additional clothing or external heating, require to better the body temperature in a cold climate. Marathon racers suffer more heat stress than walkers [1]. Furthermore, the fast movement of muscle mass during the marathon releases a large amount of metabolic heat energy, which cannot be dissipated instantly. Therefore, the rate of heat loss is not equivalent to the heat gain rate. The storage of metabolic heat energy causes skin temperature to rise. It also increases the body core temperature up to 39.5°C [2]. If the body mechanism is unable to control the body core temperature higher than 39.5°C, the temperature causes cell death [3]. Sweat evaporation is the primary factor to maintain the body temperature in a hot environment. The sweat vaporization from the skin depends on the surface area exposed to the surroundings, convective airflow around the body, ambient temperature, and relative humidity of ambient air. In a hot environment, during walking, the average sweating rate is 170.83 ml/hour [4]. During the marathon in a hot environment, the skin and the body core temperature increase due to rapid skin blood flow. At the same time, due to the body mechanism, the eccrine sweat glands become active and lose the excess heat energy in the form of sweat vapor from the body [5]. The body loses sweat upto 10 liters per day by occurring the sensible perspiration process and controls the body temperature, preventing the body from hyperthermia disorder [4]. This shows that sweat loss affects significantly by the activity level.

The evaporate rate is independent of the temperature gradient between the skin and the environment. On the other hand, sweat rate is proportional to the water vapor pressure gradient between the skin and the environment.

Convection is the heat exchange process between the fluid and the body. The heat exchange rate depends on the direction and the speed of the fluid flow. The convective heat exchange rate increases by increasing the temperature gradient between the fluid flow and the skin surface. A healthy body loses about 15% of heat energy through this process [6].

Radiation is the electromagnetic heat transfer mechanism from the sun to the body without contact with the surface. This radiation penetrates through to the earth’s surface, increases the body’s temperature, and loses heat energy through infrared radiation. If the body temperature is higher than the ambient temperature, more heat emits than is received. The body loses about 60% of heat energy by the radiation process at rest in a temperate room [6]. The sun rays of a wavelength longer than  m absorb by the atmospheric environment and do not reach the earth’s surface due to having a less penetrating force. On the other hand, the radiation of a wavelength shorter than  m does not affect the skin surface [7]. Thus, the skin’s capability to reflect the light of a shorter wavelength assists the possibilities of their radiating almost perfectly. The driving force of radiation heat loss depends on a large temperature gradient between the body and environment temperatures. The radiated heat energy by the skin surface given by Stefan’s Boltzmann as; where Er is the radiative heat energy, σ is the Stefan Boltzmann constant, is the emissivity of the skin surface depending on the emitter of wavelength and Plank constant, is the absolute temperature of the skin surface, and is the temperature of the temperature field.

Conduction is another process that occurs in the human body to exchange heat from or to the body directly in contact with other objects. If the body is warmer, then it loses heat. If the body is colder, it gains heat from the contact objects. The heat flow rate is the same in each direction if the body is in the thermal equilibrium position [8]. where is the storage heat energy, is the metabolic heat energy, Ww is the work energy, Ee is the evaporative heat energy, Er is the radiative heat energy, is the conductive heat energy, and is the convective heat energy.

Skin blood flow plays a crucial role in maintaining body temperature, whereas the blood flow rate significantly changes in the presence of a magnetic field [9]. The blood flow rate reflects by the temperature gradient between the body core temperature and the skin surface temperature [10]. The hypothalamus regulates skin blood flow and changes the skin temperature. In a cold environment, the skin blood flow decreases by the vasoconstriction process due to narrow blood vessels. It minimizes the loss of heat energy from the body. In a hot environment, the skin blood flow rate rises through the vasodilation process due to expanding the blood vessel. The continuous increase in blood flow controls by maintaining the pressure with the help of active muscle.

The parameter thermal conductivity also affects the tissue temperature. It has various values at different temperatures. The temperature of the tissue, which is above the healthy range, decreases as the thermal conductivity of the blood decreases [11].

Basal metabolic rate (BMR) is the minimum energy required for a person to sustain a healthy life at rest. The BMR depends upon the body weight [12]. In the metabolic process, the fat adipose tissue rapidly breaks down to produce a large amount of heat energy and maintains the body’s core temperature. During walking, the muscles are the primary source of metabolic heat. In the marathon, the inhaled oxygen rate of the body and exhaled carbon dioxide rate from the body also increase. These help in increasing the metabolism to provide the necessary energy.

During walking and marathon, most of the organs come in movement. Therefore, the body requires additional energy that substitutes by the metabolic process. The metabolic rate increases as much as 15 times above the resting level in the marathon [13]. In this process, the body loses the thermal energy generated by the organs. On increasing the metabolic rate, half of the increased metabolic energy uses to maintain the skin temperature and evaporates heat energy from the body in the form of a tear. The remaining metabolic energy is used to increase the blood flow rate and body core temperature. The continuous blood flow controls by the body mechanism by rapidly increasing the metabolic rate at the beginning of exercises and becomes constant after a certain period. Therefore, the metabolic rate S(t) has logistic behavior given by equation (3) [14]. where S0 is the BMR (w/m3), t is activity time (seconds), E is the activities threshold metabolism (w/m3), β is an activity controlled parameter (/seconds), and t0 (seconds) is the Sigmoid’s midpoint of the curve over time t for extensive exercise. The basal metabolism is the minimum required metabolic energy to sustain life at rest. The threshold metabolism is the upper limit of the metabolic value, and the Sigmoid midpoint is the middle of the period, from where the phenomena of the curve change.

The average BMR in a healthy adult human body is 1114 w/m3, but it may vary from person to person. The unsteady to steady behavior of metabolic rate for walking is 3889.43 w/m3 and during a marathon is 7918 w/m3 [14, 15]. The behavior of metabolic rate during the marathon and walking is graphically shown in Figure 1.

Blood perfusion is the nutritious delivery process of arterial blood to a capillary bed in the biological tissue. Changes in blood perfusion may be more complicated to detect in skin layers. A decrease in blood perfusion below normal levels is associated with ill health, and raising skin blood perfusion above normal levels has a beneficial or detrimental effect on different organs of the body [16]. Change in blood perfusion rate has an enhancing effect on thermoregulation in the human body [17]. The difference between the blood perfusion rates in the artery to the vein is negligible. In capillary, the rate of blood flow is slow, so takes place the equilibrium position.

First, the perfusion has been formulated by Pennes’ and proposed a mathematical model equation incorporated into the standard thermal diffusion and metabolism given by [18]. equation (4) is known as Pennes’ bio-heat equation.

Pennes’ analyzed the effect of thermal diffusivity, perfusion, and metabolism components in living tissue and calculated the temperature of the human forearm. He checked the validity of the temperature results obtained by his predicated equation with experimentally measured temperature in the human forearm.

By using equation (3), Pennes’ equation is modified to; where cb is the blood-specific heat capacity (J/kg °C), ρb is the blood density (kg/m3), wb is blood perfusion rate (kg/m3 s), ρ is the tissue density (kg/m3), c is the tissue-specific heat capacity (J/kg °C), K is the tissue thermal conductivity (J/ms °C), TA is the arterial blood temperature (°C), and T is the tissue temperature (°C).

Exercises, such as walking and marathon, are beneficial for health. They play a remarkable role to keeps the body healthy. However, some novelty problem arises due to the long marathon. Training and racing for marathons can also cause muscle damage due to repetitive muscle contractions. During the marathon, kidney cells become damaged by the lack of blood flow to the organs and the loss of fluid volume but typically recover within two days. For those who are participating in a long marathon, and intake insufficient oxygen, it seems the risk of sudden cardiac death is approximately 2.5 times greater than during walking [19].

Gokul et al. [20] studied the effect of blood perfusion and metabolism on temperature distribution in the human eye. They suggested that eye temperature is negligibly affected by blood perfusion and metabolism, but is affected by the parameters of blood temperature, ambient temperature, and evaporation rate. Acharya et al. [21] compared the metabolic effect in thermoregulation on the human male and female skin layers in a two-dimensional finite element method. The study deals that males lead to higher skin temperature than females due to thinner skin layers in males. Agrawal et al. [22] prepared a model for the temperature distribution in a human limb by assuming an irregular tapered shape limb with a variable radius and eccentricity. Gurung [23] prepared a model in two-dimensional temperature distribution in the human dermal region exposed at low ambient temperature with airflow. Kenefick et al. [2] experimentally studied the temperature of skin layers during the marathon and noticed that the maximum body core temperature reaches 39.5°C. This temperature is proportional to the metabolic rate and largely independent of the environmental condition during the exercise period. Khanday and Sexana [24] observed the thermoregulation and fluid regulation in the human head and dermal parts at cold ambient conditions by using the finite element method.

Previous researchers prepared models that have presented the temperature distribution in the human dermal parts with various sweating rates and ambient temperatures with constant metabolic rates. However, this model is prepared for the temperature distribution in the dermal part of the human body using the different metabolic rates produced during walking and marathon. The excess heat energy loses and keeps the body in thermoregulation by the body mechanism, which provides the realistic temperature of the dermal parts. The finite element approach has been used for numerical results and graphs of the temperature profiles.

2. Methods

In two-dimensional discretization of the dermal layers, the skin thickness ‘L’ has been taken along Y-axis, and the width of the skin ‘W’ has been taken along the X-axis on the skin surface. The outer surface of the body has exposed to the environment, so heat loss occurs by convection, radiation, and sweat evaporation. The mixed boundary condition of heat flux from the outer skin surface is given by [21]: where hcv is the convective heat transfer coefficient between the skin and the temperature field (w/m2 °C), T is the atmospheric temperature (°C), σ is the Stefan Boltzmann constant (5.67 × 10−8 w/m2 °C), ε is the emissivity of the skin surface, L is the latent heat of evaporation (J/kg), and E is the evaporative heat loss between the skin surface and the temperature field (kg/m2 s).

The radiation term presented in equation (6) is nonlinear. Therefore, it is difficult to find the solution of equation (6). For the solution of the problem, we introduce a suitable iterative method. The above nonlinear boundary condition stated in equation (6) can be written as;

If the value of the term is known, equation (6) can be viewed as a generalized convection condition between convection and radiation. For this, we introduce the iterative algorithm: with where hcr is the combined heat transfer coefficient of convection and radiation with radiation coefficient .

The term is the temperature sequences for . . . , and is the initial guess temperature of the skin surface. The iteration process has completed when the convergent condition is satisfied. where is the specified interaction tolerance.

The transport of heat within tissue occurs along normal to the skin surface from the body core, and hence we assume negligible heat flux in the x-direction. Therefore, the boundary conditions are assumed to be:

During marathon and walking, the human body maintains its body core temperature uniformly by 39.5°C and 37°C, respectively. Therefore, the inner boundary condition has taken as: and are the Dirichlet’s boundary conditions during walking and marathon, respectively, and is the body core temperature.

3. Skin Geometry and Assumption Parameters

3.1. Skin Geometry

Skin plays a protective role and performs various functions in the thermoregulatory process. The skin thickness in the human body has divided into three primary layers: epidermis, dermis, and subcutaneous tissue. Each skin layer has its own thermo-physical and optical properties [25]. The outer layer epidermis provides mechanical strength and rigidity to the skin’s structure. The epidermis contains melanocytes that produce melanin pigments. Melanin helps to change the skin’s color and protect the body from ultraviolet radiation when the sunlight incident on the skin’s surface. The middle layer dermis contains blood vessels, elastic fibers, and collagen. The primary role of the dermis is to support the epidermis and enable the skin to thrive. It is responsible for the contraction or dilation of the blood vessels to maintain homeostasis body. During walking and marathon in a hot environment, the blood volume increases due to dilating the blood vessel and releasing the heat in the form of sweat.

The subcutaneous tissue is the inner layer of tissue, which has composed of adipose fat cells. It works as insulation to maintain our body’s core temperature and prevent hyperthermia sickness when the body exposes to a hot environment. In the model, the skin thickness has considered as a two-dimensional rectangular skin region. The skin diameter (width) along the X-axis is  m, and the total thickness along the vertical Y-axis is  m. Initially, the skin is divided into 380 elements with a triangular shape and has 220 nodes. The epidermis layer, dermis layer, and subcutaneous tissue have divided into 40, 140, and 200 triangular meshes and have 33, 88, and 121 nodes, respectively, as shown in Figure 2. Figure 3 depicts the discretization of skin thickness into the epidermis, dermis, and subcutaneous tissue with their respective nodal temperatures, , , and .

A triangular element ‘e’ has three global Cartesian coordinates: , , and has taken as shown in Figure 4. Each node temperature has expressed in the form of shape (interpolation) functions [26]. We assume that the temperature field (trial function) over the element ‘e’ is given by

On reducing the nodal temperature at the nodal point 1, 2, and 3, we get:

On solving equation (15) we get,

Let us assume, area of triangular element = and

Then equation (13) becomes as,

This provides, , , and .

Using , , and in equation (14), finally we get, where is a shape function written as, .

3.2. Assumption Parameters

Since the epidermis layer is composed of the dead cell, therefore, the arterial blood temperature TA, blood perfusion rate , and metabolic heat generation rate have been taken to zero in the epidermis layer [21]. The arterial blood temperature TA in the dermis and subcutaneous tissue are assumed equal to Tb. During walking and marathon races, the adipose tissue becomes more active and increases metabolic energy. Therefore, the metabolic rate in the subcutaneous tissue is taken double that of the dermis layer. Since the blood perfusion rate has been considered a function of depth, however in this model, the perfusion rate is considered as same in the dermis and subcutaneous region.

4. Solution of the Model

The governing equation (5) together with boundary conditions is transformed into the variational form: where is the domain of the region of the dermal parts and is the outer boundary of the skin. On expressing the function I as the sum of E elemental quantities as where Here, is the domain of the element ‘e’. Equation (22) is broken into six parts as below. where

For the minimization of the function I, we have

Using equations (22) and (25), we get,

Differentiate equation (26) with respect to each nodal temperature T0, T1, and T2 set the derivative equal to zero for the minimization. Equation (26) leads to a linear system of differential in the form

The system of equation (27) can be solved by using the Crank–Nicolson method concerning time with the following relation where Δt is the time interval and is the interface temperature of the epidermis, dermis, and subcutaneous tissue in the nth time step. On solving equation (28) repeatedly, we get the nodal temperature of each layer.

For the steady case of the model, we obtain the system of algebraic equations in matrix form as;

5. Numerical Results

The threshold value of metabolic rate during walking is 3889.43 w/m3, and during the marathon is 7918.00 w/m3. The value of skin layers’ thickness and physiological parameters used for numerical simulation has been taken as shown in Tables 1 and 2, respectively.

Assuming the ambient temperature is less than 37°C, the tissue temperature increases from the skin surface towards the body core. The increased tissue temperature shows in linear order given by the equation:

The initial skin temperature is considered 24.91°C at normal ambient temperature. The use of in equation (30) is constant, whose numerical value is determined by taking the known value °C for walking and °C for the marathon at .

5.1. Steady State Temperature Results

The steady-state temperature results of the natural skin layers for the epidermis, dermis, and subcutaneous tissue at the various ambient temperatures °C, 25°C, and 35°C, and sweat evaporation rates , 0.00002, and 0.00004 kg/m2 s during walking and marathon are calculated. The obtained results are presented through the graph in Figures 5, 6, 7, and 8 and Tables 3, 4, and 5.

On comparing the temperatures between Figure 5 and Figure 6 during walking and marathon at  = 15°C, 25°C, 35°C, and  kg/m2 s, the skin layers temperature , , and have more by 2.38°C, 2.47°C, and 2.55°C, respectively, during the marathon compared with walking at each ambient temperature. These are due to more metabolic effects during the marathon. These results also have shown that the temperature reached 39.60°C at during the marathon without sweat evaporation. It may cause hyperthermia disorder. Therefore, the body controls the temperature by sensible perspiration process and keeps the body in thermoregulation.

Figure 7 shows the comparison of the skin layers’ temperature during walking and marathon at °C and E = 0.00002 kg/m2 s. The graph of the results delineates that is more by 2.38°C, is more by 2.47°C, and is more by 2.55°C during the marathon compared with walking. These are due to more heat energy produced by the body during the marathon, and excess heat energy release helps to control the body temperature.

Figure 8 has presented the skin layers’ temperature during walking and marathon at °C and  kg/m2 s. These results reveal the temperatures , , and control to 34.93°C, 35.60°C, and 36.61°C during walking and 37.31°C, 38.08°C, and 39.15°C, by releasing sweat during the marathon, respectively. Due to high ambient temperature, the body gets energy from the hot environment and increases the interface temperature. The study indicates that the temperature distribution has an enhancing effect on heat flux at the skin surface due to the exercises. On the other hand, the body plays a vital role, increasing the sweating rate has to reduce the temperature of the living tissues and keeps the body in thermoregulation during exercise.

5.2. Unsteady State Temperature Results

The unsteady state temperatures solution for the epidermis node, for the dermis node, and for the subcutaneous node are carried out by solving the system of equation (27). Figures 9, 10, 11, and 12 and Table 6 illustrate the unsteady state to steady state temperatures during walking and marathon at °C, 25°C, and 35°C at various evaporation rates.

Figures 9 and 10 indicate the unsteady state temperatures, at  kg/m2 s in different ambient temperatures during walking and marathon. The results reveal that the temperature is more by 2.38°C, is more by 2.47°C, and is more by 2.55°C during the marathon compared with walking at °C. The temperatures, , , and is more by 2.37°C, 2.46°C, and 2.54°C, respectively, during a marathon compared with walking at T = 35°C. These are due to the body producing more metabolic heat energy during the marathon compared with walking. Figures 9 and 10 indicate that each temperature increases on increasing the ambient temperatures from 15°C to 35°C, and shows the highest temperature at °C. Figures 9 and 10 also delineate that is less deviated and has reached faster in steady state from the unsteady state than and during both walking and marathon.

In comparing the unsteady state temperatures during walking and marathon as presented in Figures 11(a) and 11(b). Figures 11(a) and 11(b) delineate that each nodal temperature has more during the marathon compared with walking. These are due to the more metabolic effect and conduct more heat energy during the marathon.

Figure 12 indicates the temperatures and during walking and marathon at  kg/m2 s and °C. These results delineate that the temperatures, and , rise higher and reach faster to the steady state during the marathon compared with walking. These are due to the faster movement of muscle mass during the marathon. The temperature is affected more due to the closer the subcutaneous tissue is to the body core.

5.3. Validity of the Temperature Results
5.3.1. Temperature during walking

Procter et al. studied the temperature occurring in the human body during moderate-intensity exercise. They suggested that the body core temperature occurs at °C during moderate-intensity exercise [28].

de Andrade Fernandes et al. suggested that the maximum body core temperature occurs at 37.70°C during moderate-intensity exercise. They indicated this experimental result in the published article [29].

Since the subcutaneous tissue is closer to the body’s core, in our results, we observed that the maximum temperature of subcutaneous tissue during walking (moderate-intensity exercise) is 37.05°C. This result is closely equal to the body’s core temperature during moderate-intensity exercise suggested by the researchers [28, 29]. Therefore, the temperature results obtained in this study during walking are valid.

5.3.2. Temperature During Marathon

Del Coso et al. suggested in their experimental study that during a 5-km marathon race, the body core temperatures rise rapidly at the beginning of the race and control the body core temperature up to 39.50°C [30].

Kenefick et al. suggested that the body core temperature in a human body occurs ranges from 38.50°C to 39.50°C during a marathon [2].

In this study, we obtained the maximum body core temperature during the marathon is 39.60°C. The result is nearly equal to the body core temperature suggested by the researchers during the marathon [2, 30]. Therefore, the temperature result conveyed in this study during the marathon is valid.

6. Stability and Convergence Analysis of the Solution

6.1. Stability Analysis

To show the stability of the solution of an equation in the two-dimensional discretization domain, it is sufficient to show the discretization is stable. For this, we use the following theorem.

6.1.1. Theorem

Let be a continuous finite element space defined on the triangulation . f and b are the linear and bilinear forms of the problem . Then finite element approximation exists and discretization is stable in the norm.

6.1.2. Proof

To prove the theorem, first, we have to show that F and b are continuous.

Where,

It shows that F is continuous.

Using Schwarz inequality of the inner product , the bilinear form b is coercive. Here, the continuity and coercivity constants both are 1, independent of space step size h, therefore, the discretization is stable.

6.2. Convergence of the Solution

Let be the approximation solution with the exact solution T, of the problem. Then, for the convergence of the solution, it is sufficient to show that where k is the degree of Lagrange finite element approximation . If m is the solution to the problem, then, Aubin–Nitsche duality argument provides where m has the same Neumann condition as T.

Since the elliptic regularity provides that

Then, using the testing function v and integrating by parts we get, and so denotes the global interpolating mapping with dividing both sides by we get gives the result. Thus, the solution is convergent with space step size h by using the norm instead of the norm.

7. Conclusion

The results in the model analyze that the steady temperature of skin layers is higher during the marathon than walking due to more metabolic effects. These results indicate that the temperature of the human body increases by increasing the activity level. However, due to the threshold value of body core temperature, the body loses more heat energy during the marathon compared with walking in the form of sweat and frequently controls the body temperature. It indicates that the sweating mechanism plays a vital role in protecting the body from hyperthermia sickness during walking and the marathon.

Earlier researchers developed many mathematical models in the temperature distribution in the human dermal parts, but they have not determined the metabolic rate during the exercise period. Therefore, this model is prepared for the realistic temperature distribution in dermal parts of the human body with the metabolic rate during walking and marathon. The various sweat evaporation rates and ambient temperatures have been used in this model. Therefore, this research study will assist to maintains the physical structure of aging and pregnant women. It also uses to develop the model regarding the different exercises as a sports player, mountain climber, laborer, plumber, typist, and driver based on their physical and physiological parameters.

Data Availability

Data supporting this research article are available from the corresponding author or First author upon reasonable request.

Conflicts of Interest

The author(s) declare(s) that they have no conflicts of interest.