Modeling on Metabolic Rate and Thermoregulation in Three Layered Human Skin during Carpentering, Swimming and Marathon

Metabolisms play a vital role in thermoregulation in the human body. The metabolic rate varies with the activity levels and has different behaviors in nature depending on the physical activities of the person. During the activity, metabolic rate increases rapidly at the beginning and then increases slowly to become almost constant after a certain time. So, its behavior is as logistics in nature. The high metabolic rate during activity causes the increase of body core temperature up to 39˚C [1] [2]. The logistic model of metabolic rate is used to re-model Pennes’ bioheat equation for the study of temperature distribution in three layered human dermal parts during carpentering, swimming and marathon. The finite element method is used to obtain the solution of the model equation. The results demonstrate that there is a significant change in tissue temperature due to sweating and ambient temperature variations.


Thermoregulation
Thermoregulation is the process of transporting thermal energy through the biochemical process in the human body. It maintains its required internal temperature. It keeps the body temperature in the equilibrium position. This process is applicable for diagnostic and therapeutic applications and involves either mass or heat transfer. In the human body, the heat transfer is affected by blood vessel geometry, local blood flow rates, the thermal capacity of blood and produces metabolic energy. The Body regulates its temperature through internal metabolic processes and internally maintains a narrow range of internal temperature. Heat is usually generated by the metabolic process but under the condition of excessive cold, the body generates heat by shivering. Heat is lost and gained through the process of convection, radiation, blood perfusion and conduction, while evaporation contributes only heat loss from the body to the environment. Total heat loss from body surface depends on the temperature difference between skin and environment.
Hypothalamus is the main part of the brain that controls thermo-regulation. When it senses the internal temperature becoming too low or too high, it sends signals to different organs such as muscles, glands, and nervous. The body mechanism responds and helps to maintain the temperature to normal position. Different people have different thermal behaviors even in the same environment.
The biochemical processes have been divided into three broad categories: hypothermia, hyperthermia and cryobiology. Hypothermia is the phenomena in which the body core temperature falls to 35˚C or below it. This phenomenon can potentially lead to cardiac arrest, brain damage, or even death. In Hyperthermia, body core temperature rises to 42˚C or above it. Due to this phenomenon, it can suffer brain damage or even death, where cryobiology is the blood subfreezing temperature period. The average healthy person has a normal body core temperature between 36.1˚C to 37.8˚C [3]. But temperature rising up to 39.5˚C during the high intensity exercises doing on long period [1] [2]. However, if the body temperature is extreme, it can affect the body's ability to work. Thermo-regulation phenomena are affected by either environment or biological factors such as metabolic rate, dehydration, gender, etc. The internal temperature of the body rises due to fever, physical exercise and digestion of the food. Besides the internal temperature decreases due to use of alcohol, drug and some functioning of thyroid glands, etc. [4].

Circulatory System
The circulatory system also plays a vital role to maintain body temperature. The vasodilation and vasoconstriction are the two main processes which maintain heat and balanced the temperature of the body. If the body has a high temperature, the body controls the temperature and keeps in normal by occurring the mechanism of vasodilation and sweating. The blood flow rate to the skin increases by expanding the wider blood vessels under skin and heat energy dissipated by the radiation process. On the other hand, sweat glands lying under the skin surface release sweat in the form of water and vapour from the skin surface and maintain the body temperature is normal. If the body has a low temperature, it maintains the body temperature by the mechanism of vasoconstriction and thermo-genesis. The blood flow rate to the skin reduces by contracting the peripheral blood vessels lying under the skin surface and conserving heat energy near the warm inner body and keeps the body temperature normal. Besides body muscles, organs, glands and hormones also produce heat by their mechanism. Figure 1 represents the thermoregulation system in the human body. It shows that if the body has a low temperature, it maintains by the process of vasoconstriction and thermogenesis and if body temperature is high, it controls by the process of vasodilation and sweating to keeps the body temperature is normal.

Metabolism
The metabolic rate of a person increases due to increases in the kinetic energy and helps on increasing the temperature of the body. Some of the hormones and thyroid gland releases hormones to increase metabolism. This increased metabolism creates energy in the form of heat and maintains body temperature as normal. Basal metabolism is the minimum amount of energy release in the body to sustain life in the rest position. Body uses this energy for the circulation of blood, respiration and functioning of cells and different organs. The active person has a more metabolic rate than a less active person [5]. In general, the metabolic rate is proportional to the body weight and depends upon the type of exercises, body surface area, health, mental state, sex, thermal conductivity, age and atmospheric conditions, etc. Due to the hormones, the metabolic rate is highly increased in fever [6]. The metabolic heat produced by the exercising muscles is transported by the circulating blood to the surface of the body and released to the environment either by radiation and convection or by evaporation of sweet appears in a hot environment. According to the conservation law of energy in a healthy body, the amount of heat gain is equal to the amount of heat loss, and is given by [7] M W R C E H = ± ± + +  when the body is at rest, there is normal blood flow and metabolic activity however, during exercise there is abnormal blood flow and metabolic activity. When a person engaged in different types of physical activities, the body requires additional fuel which increases the metabolic rate and the rate of heat production [8].
The body must use its mechanism to remove the additional heat produced to keep the internal temperature at a healthy level. The body always produces heat, so the metabolic heat (M) is always positive, varying with the activity level. When environment temperature is lower than the body temperature, a thermal gradient is created which favors heat loss from the body core to the environment. That heat energy is transformed into the skin by convective blood flow and emitted by radiation, convection and sweat evaporation. The driving force of radiation and convective heat loss depends on the maintains of a large temperature gradient between the environment and body. When the environment temperature exceeds than body temperature, the gradient for heat exchange is reversed and the body gains the heat. Different people have different thermal behavior even in the same environment. However, as environment temperature rises, sweat rate also increases. In a resting lightly clad man, the sweet begins at 29˚C environment while in the resting nude man sweating begins at the environment temperature at 31˚C [9]. Sweat evaporation becomes the primary means of heat dissipation when convective and radiative heat exchange is minimal. In a hot moist environment, a healthy man approximately 1.94 ml sweat losses from the body [10]. Being lost of heat through evaporation even an entirely high metabolic rate can also be limited by the body mechanism. If the combined heat loss from the body is less than the metabolic heat rate production of the body, the excess heat must be store in the body, but storage heat energy (H) is always small because the body has a limited thermal storage capacity. This storage heat energy helps on bring the body temperature rise and when the metabolic heat is less than combine heat loss, the body temperature falls.
During the mild exercises, the muscles are the main source of metabolic heat, and during intense exercise, it may occur up to 90% [11]. In heavy exercises, the respiration process plays the role of increasing metabolism. In low-intensity exercises as sitting, typing, cooking, etc there is very little difference in the vascular system so the metabolic rate is similar to the basal metabolic rate. While in high-intensity exercises, as marathon, wrestling, etc most of the organs come in movement with high kinetic energy so, there is a loss of energy. In moderate-intensity exercise as carpentering, driving, etc. helps on decreasing body weight by increasing cardiovascular exercise.
From the experimental data, the average basal metabolic rate (BMR) of (30 -60) years old healthy man, of body weight 74 kg has 1737 kcal/day [12]. This shows that the BMR of the healthy human body has 1114 w/m 3 , but it may vary from person to person. On the contrary to the basal metabolic energy, the body produces additional metabolic energy during the activity period. From the next experimental data, the energy consumption during different activities carpen-tering, swimming and marathon (athletics) by a man of 65 kg weight consumed 16.7 kJ/min, (21 -31.5) kJ/min, and 31.5 kJ/min energy respectively [13]. This data provides that the estimation value of the metabolic rate of carpentry person is 4198.08 w/m 3 and the metabolic rate of the swimming person is 6598.79 w/m 3 . The scale of metabolic rate during various activities is shown in Table 1.

Blood Perfusion
Blood perfusion is the physiological term that refers to the process of delivery of arteries blood to a capillaries bed in the biological tissue. The arterioles and venous blood temperature may be different from the local tissue temperature. The rate of heat transfer between blood and local tissue is proportional to the product of a volumetric perfusion rate and the difference between the arterial blood temperature and the local tissue temperature [14]. This temperature may vary as a function of many transient and physiological and physical parameters. Pennes' mathematical model describes the effect of blood perfusion and metabolism on the energy balance within the tissue. These effects were incorporated into the standard thermal diffusion equation and modeled the bio-heat equation Several computer simulated methods are developed for the estimation of temperature distribution in the human body. Saxena et al. [15] used a variational finite element approach with linear shape function to find one dimensional unsteady temperature distribution in epidermis, dermis and subcutaneous tissues (SST) region assuming a rate of blood flow and rate of metabolic heat generation as a variable in the dermis part. Chao et al. [16] assumed two simple models and studied the temperature distribution curve in skin and subcutaneous tissue for certain paramaters values are constant. Agrawal [13] and Kenefick et al. [2] experimentally studied the amount of energy expenditure and blood flow rate during swimming and marathon. They estimated the maximum core temperature occurring 39.5˚C in the human body during the exercise period. Acharya et al. [17] divided the dermal part into six skin layers and studied the metabolic effect in thermoregulation on human males and females. The authors suggested that skin temperature in males has more than females. Khanday and Saxena [18] assumed five layered skin in the human body and used one dimensional steady state model to the estimation of cold effect on the human dermal part. Saxena and Gupta [15] and Saxena and Arya [19] contributed the papers on the effect of blood flow and heat flow in human skin and subcutaneous tissue by using variational finite element method. Kumari and Adlakha [20] developed a numerical model to study the temperature distribution in human peripheral regions incorporated the blood mass flow rate, thermal conductivity and metabolic heat generation rate were constants during and after the exercise.
Gurung et al. [21] investigated the temperature distribution in three layers of a human dermal part in one dimensional unsteady state with quadratic shape function assuming that the outer skin is exposed to the atmosphere. Agrawal et al. [22] developed a model for the temperature distribution in a human limb by assuming an irregular tapered shaped limb with a variable radius and eccentricity. Khanday [23] explained the appearance of thermal stress on the human brain tissue in hypothermic conditions. Khanday and Sexana [18] studied the thermoregulation and fluid regulation in the human head and dermal region at cold environmental conditions by using variational finite element method. K.C. et al. [24] investigated the thermal effect of eyelid based on the properties of ambient temperature, evaporation rate, blood temperature and lens thermal conductivity in the human eye. Previously, developed models have not studied on the temperature distribution in the human dermal part during the exercise. So this mathematical model has presented to estimate the metabolic energy produces differently in various exercises. The main objective of this study is to investigate the temperature profiles of epidermis, dermis and subcutaneous tissues during exercise. Since the body is an irregular geometry, the finite element method is appropriate to handle such an inhomogeneous discretized problem to get realistic values of temperature of different layers.

Model for Metabolic Energy
During physical exercise, metabolic rate increases due to an increase in the rate of blood flow. The continuous increases in blood flow are controlled by the mechanism of the body by rapidly producing metabolic energy at the beginning of the exercise and become constant after a certain time so it is plausible to consider the metabolic rate, increasing logistically in the form similar to the logistic curve. We consider the metabolic rate S(t) equation based on exercise as ( ) ( )

Discretization of Domain
Skin is the main organ that keeps helping the temperature balance in the human body. If any illness occurs in the body, the first symptom is changing body temperature [8]. In mathematical treatments of temperature distribution in the human dermal part, the skin layers can be regarded as a physical and physiological barrier with complex structures. The three natural layers of skin are epidermis, dermis and subcutaneous tissues (SST). The schematic diagram of the temperature distribution model in dermal parts of the human body is shown as in Figure  3. T be the temperature function of epidermis, dermis and subcutaneous tissues respectively.

Solution of the Model
The governing equation that characterized the heat regulation in in-vivo tissue of human body during exercises is given by the partial differential Equation (1)  ( ) where, The outer surface of the skin layer is exposed to the environment during the exercise period. Therefore heat flux is dissipated for 0 x = and heat loss takes place from the outer surface due to conduction, convection, radiation and evaporation methods. So the net heat flux is calculated by mixed boundary, 2) Boundary condition at 3 x l = (body core) During Exercise, due to the rapid movement of muscle mass and an increase in heart rate, the body produces a large amount of metabolic heat energy and cannot dissipate all energy instantaneously at that time. The rate of heat loss does not offset the rate of heat gain so some heat energy stores in the body. That excess (stored) heat energy helps to increase the body core temperature up to 39˚C. So, the Dirichlet's inner boundary condition during exercise period is taken as Using Euler-Lagrange formula, the variational integral form of Equation (3) and Equation (4) is given by In the model, the physical and physiological parameters depending on the layers of dermal part and are considered as given in Table 2 A

B T C T D T E T F T T T T T T t
A

B T C T D T E T F T T T T T T t
constants whose values depend upon physical and physiological parameters of dermal part as given in Table 4. As a next step to finite element method, we differentiate 1 I , 2 I and 3 I with respect to 0 T , 1 T and 2 T and set d 0 On simplification, we obtain system of equations in matrix form where,

Numerical Results
The threshold values of metabolic rate during physical activities: carpentering, swimming and marathon running are 4198.08 w/m 3 , 6598.79 w/m 3 and 7918 w/m 3 respectively. The values of physical and physiological parameters used for numerical simulation are taken as shown in Table 3 and Table 4.
In normal condition when the atmospheric temperature is below body core At normal atmospheric temperature, initial skin temperature is considered 21˚C. The use of µ in Equation (8) is constant, whose numerical value is determined by taking the known value We use the iterative method and the Crank-Nicolson method to solve the Equation (10). The Crank -Nicolson method gives Here, t ∆ is the time interval and ( ) 0 T is the initial nodal temperature in 3 1 × matrix form.
For the steady case of the model, we obtain the following matrix form of the system of algebraic equations

Steady State Results
The results of the analysis for dermal layers temperature distribution during marathon, swimming and carpentering in steady state are presented through Figures 4 to 6 and Table 5.    From the graph in Figure 4 it is observed that there are significant variations in temperature distribution in the layers of skin due to change in atmospheric temperature in carpentering, swimming and marathon. But no significant changes in temperature distribution are observed during these activities at a given atmospheric temperature. From the graph in Figure 5(a), it is observed that the epidermis layer has more temperature by 1.90˚C, the dermis layer has more temperature by 1.95˚C and subcutaneous tissue has more temperature by 2.02˚C at From the graph in Figure 5(b), the result shows that the epidermis layer has more temperature by 1.96˚C, the dermis layer has more temperature by 2.02˚C and subcutaneous tissue has more temperature by 2.06˚C at From the graph in Figure 5(c), it is observed that the epidermis layer has more temperature by 2.00˚C, the dermis layer has more temperature by 2.05˚C and subcutaneous tissue has more temperature by 2.09˚C at ⋅ at fixed environment temperature 25˚C in carpentering, swimming and marathon. This shows that temperature of dermal layer control by producing the essential sweat from the body when the environment temperature is high during the activities.
In the present Table 5, the table shows the steady state temperature of three skin layered during different activities at various ambient temperatures. The temperature of each dermal layer during the marathon is more than swimming and carpentering at each ambient temperature 15˚C, 25˚C, and 35˚C. The temperature of each dermal layer during swimming is more than carpentering due  Table 6 and Table 7. Figure 7 represents the tissue temperatures during the different activities at    On observing Figure 8, at 25 C T ∞ =  and 2 0.00004 kg m s E = ⋅ , the subcutaneous tissue temperature is greater by 3.20˚C than epidermis layer and is greater by 1.97˚C than dermis layer temperatures during carpentering. The subcutaneous tissue temperature is greater by 3.18˚C than the epidermis and is greater by 1.96˚C than dermis layer temperatures during swimming. During the marathon, the subcutaneous tissue temperature is greater by 3.17˚C than the epidermis and is greater by 1.94˚C than dermis layer temperatures. These results exhibit that subcutaneous tissue temperature is higher than dermis and epidermis layers because subcutaneous tissue is nearer from the body core with constant body core temperature. Figure 9 reveals the temperatures variation between the dermis and subcutaneous tissue is higher than the temperatures variations between the epidermis and dermis layers at 35 C T ∞ =  and 2 0.00004 kg m s E = ⋅ during carpentering, swimming and marathon.
In present  Table 6 and Table 7, the results exhibit that the temperature of the skin layers fall due to loss of more heat energy in the form of sweat at high sweat rate.

Discussion and Conclusion
The metabolic rate increases due to increases of the oxygen consumption by the body during physical activities. It helps to increase the temperature of the body.
When the environment temperature increases, the metabolic rate in the body decreases. In the human body, the normal core temperature is 37˚C so, on or above 37˚C, the body produces negligible metabolic heat [17]. In this model, the ambient temperatures 15˚C, 25˚C and 35˚C are taken.
More active physical activity produces more sweat due to the fast movement of muscle mass. The analysis also shows, carpentering produces less metabolic energy than swimming and swimming produces less metabolic energy than the marathon. In this model sweat release rates on the exercises are taken as 0.00004 kg/m 2 s and 0.00008 kg/m 2 s. temperatures. This is due to increase of the metabolic rate on increasing the lean muscle mass and lose adipose fat tissue in subcutaneous tissue during activities.
These results also demonstrate that the steady-state temperature of each dermal layer is achieved more during the marathon than swimming and carpentering due to more metabolic rate during the marathon. The unsteady and steady both results execute that the environment temperature gives significant changes in temperature in human dermal part.
The above result suggests that the temperature of each node increases by decreasing the sweat evaporation rate and vice versa at the same ambient temperature. This shows that the sweat evaporation rate is the catalyst in the thermoregulation process. All of the previous researchers developed models have not determined the temperature distribution in human dermal part by the computational method during the activity period. So this model is developed for the temperature distribution in dermal parts of human body during activities. The thermal effect of physical activities is quite significant. It also uses to develop models regarding physical work for labourers, military, other sports persons, etc.
based on their physical and physiological characters.