Mathematical Simulation and Numerical Computation of the Temperature Profiles in the Peripherals of Human Brain during the Tepid Sponge Treatment to Fever
Background. Fever is one of the frequently occurring diseases in human beings, and the body is said to have befallen in fever if the arterial or internal body temperature rises to 38°C. The patient who suffers from fever is either given paracetamol or tepid sponging or both. Objective. This paper is aimed at studying the effects of the tepid sponge in normalizing the high temperature of the human body during fever. Among the various available methods for tepid sponging, the impact of holding a cool wet cloth on the forehead for reducing the fever is analyzed and pictured graphically. Method. For analyzing the effects of tepid sponge on the temperature distribution of the domain consisting of scalp, skull, and cerebrospinal fluid (CSF), a cool wet cloth is brought in contact with the skin allowing the heat to transfer from the brain to the wet cloth through these layers. The heat transfer in living biological tissues is different from ordinary heat transfer in other nonliving materials. Therefore, a model based on the bioheat equation has been constructed. The model has been solved by numerical methods for both steady- and unsteady-state cases. The domain, which consists of the scalp, skull, and CSF layers of the human head, has been discretized into four equal parts along the axes of the three-dimensional coordinate system. The forward difference and forward time centered space approximations were employed for numerical temperature distribution results at the nodal points. Results. The effects of tepid sponge in reducing the body temperature with fever at 38°C, 39.5°C, and 41°C have been numerically calculated, and the results were pictured graphically. For transient cases, the corresponding calculations have been carried out at times , 4 minutes, and 6 minutes. Conclusion. Among all the available remedies to fever, tepid sponging has shown a significant effect in controlling fever.
Fever is one of the common and universally known diseases from which human beings frequently suffer. Fever is generally the body’s response to infection . Fever is said to occur when the oral temperature exceeds 38°C. The normal range of oral temperature of the human body is 37.2–37.7°C as asserted by Mackowiak et al. . The findings of Mackowiak et al.  conflicted with Wunderlich’s definition of the normal temperature of the human body as the latter had ranged it between 36.8 and 37°C. The authors exhibited comparable thermal variability and observed experimentally that the mean oral temperature of our body varies with the time of day (morning, noon, or evening) and sex (men or women). Morrison et al.  demonstrated that the central neural circuits orchestrate a homeostatic repertoire to maintain the body temperature during abnormal environmental conditions and to alter body temperature during the inflammatory response.
Man has always been striving for remedies to various diseases and has been experimenting with drug placebo and physical methods for optimum outcomes. Treatment of fever by tepid sponging at home, especially to children by caregivers prior to a doctor’s consultation, is a common practice, and it changes the temperature of the body significantly, as analyzed by Hendrawati and Elvira . Tepid or cold sponging is a remedy where the application of moist cold liquid reduces the excess temperature of the body. Although, this treatment at home is generally successful but sometimes may lead to complications due to mishandling by a person with negligible medical competency. The treatment techniques employed include exposing the patient to air, tepid sponging, and using drugs such as paracetamol and other traditional remedies . Among all these treatments, tepid sponging is most commonly practiced. The sponging action is a process where the water film by cool wet cloth is constantly pressed and moved over a part of the body, especially the forehead, resulting in the replacement of the warm water with a fresh layer of cooler water on the skin surface, thus helping in maximizing the outward heat transfer through conduction, convection, and perfusion .
Aluka et al.  compared the treatment by cold water sponging with paracetamol on children and analyzed the effects of both for the reduction of fever temperature. They substantiated that as compared to oral paracetamol, cold water sponging yields a more speedy reduction in body temperature of febrile children in a tropical environment, and the results are more effective when the sponging is ongoing. Meremikwu et al.  also observed the effect of physical methods versus drug placebo or no treatment for managing fever in children. Further, the comparisons of the various treatments like alcohol versus tepid sponging, iced water versus tepid sponging, and physical methods versus drug injection were carried and analyzed. Agbolosu et al.  studied the efficacy of tepid sponging over paracetamol on a sample of 80 children, among which 40 were enrolled for paracetamol treatment and the remaining 40 for tepid sponging. They concluded that in a tropical climate, paracetamol is more effective in reducing fever than tepid sponging. However, tepid sponging together with paracetamol is highly effective in serving the purpose of reducing the patient’s excess temperature due to fever. Eman  considered the problem of heat transfer in the skin tissue of the human head subjected to thermal diffusion with relaxation time and concluded that the values of the time and the distance have significant effects on the temperature increment and temperature concentration. The values of the temperature increment and the concentration increase when the time increases, and their values decrease when the distance increases. Shen et al.  examined the heat transfer and elastic deformation in soft tissues by using Penne’s bioheat transfer equation and the modified Duhamel-Neumann equations. A 19-point finite-difference scheme was used to solve the three-dimensional governing equation. Ferreira and Yanagihara  also solved the model based on 3D heat conduction and used the elliptical cylinders to adequately approximate the body geometry.
The present paper explores the three-dimensional heat transfer and corresponding temperature distribution during tepid sponging treatment to fever. The water, being a compound with high specific heat capacity, is an excellent and readily available material for extracting the heat from the patient’s body. Thus, when a piece of wet cloth is kept on the forehead, it immediately cools down the temperature of the scalp and hence of the skull and CSF. This paper derives a mathematical estimation of temperature distribution in the scalp, skull, and CSF (the surrounding layers of the human brain) during tepid sponging treatment in fever. The temperature of the skin decreases due to its contact with a piece of wet cloth which is repeatedly damped in cold water, wrung for any extra moisture, and applied to the forehead. The heat (due to fever) from the body is taken away by the water resulting in cooling down the temperature of the forehead and, in turn, reducing the fever in subsequent regions. Aijaz et al.  elaborated how the water that diffuses into the skin plays a significant role in the fluctuation of temperature profiles near the skin surface and subsequent regions. Khanday et al.  estimated the fluid distribution patterns in human skin at various values of metabolic heat generation. Further, Tsunetsugu and Sugiyama  compared the physiological changes in the human body which is in direct contact with various materials especially wood. In a similar way, during tepid sponging treatment, the skin is in direct contact with wet cloth; therefore, heat flows out through biological tissues involving conduction, convection, and perfusion processes. Further, various physiological parameters have a remarkable impact on the heat transfer through tissue as assessed by Ley and Bayazitoglu . Therefore, for the formulation of the mathematical model, Pennes’ bioheat equation  is the most appropriate equation as discussed by Aijaz et al. . Therefore, the rate of heat flow through the human head is also dependent on the layer under study. Hence, for better results, it becomes imperative to discretize the domain into small elements. In this direction, the model has been solved for both steady- and unsteady-state cases by the forward difference approximation and the forward time centered space (FTCS), respectively. The uniform grid of dimensions is introduced, and then, the differential equations for internal node have been solved as demonstrated by Majchrzak and Turchan .
2. Materials and Methods
We intend to formulate a model for estimating the heat transfer in the domain consisting of scalp, skull, and CSF during tepid sponging. Therefore, to make the domain compatible with the three-dimensional coordinate system, we take both a piece of wet cloth and the part of the forehead which is in contact with it as cuboidal shaped. Let the domain be of length (along the +ve -axis), height (along the -axis), and width (along the -axis). , , , and are five points on the surface of wet cloth equidistant from the origin as shown in Figure 1(a), and is the extreme point (the interface between CSF and brain) up to which heat transfer is to be estimated. The diffusion equation applicable to the present in vivo heat transfer is the three-dimensional bioheat equation given by where , , , and are, respectively, the temperature, the thermal conductivity, the density, and the specific heat capacity of the tissue. is the blood perfusion rate, is the specific heat of the blood, and is the arterial blood temperature, but in the present paper, it is fever temperature which is 38°C, 39.5°C, and 41°C for three different cases. and are, respectively, the metabolic heat generation and effect by external cooling source. In this paper, the external cooling source is the wet cloth used for tepid sponging at a uniform temperature of 7°C.
It is to be noted that , , and have specific constant values in each layer of the domain and, hence, take distinct values in scalp, skull, and CSF as demonstrated by Aijaz et al. . Thus, in the numerical computation, their numerical values will be taken according to the positions of nodal points pertaining to in the numerical discretization of the domain (see Figure 1(b)). The appropriate boundary conditions for Equation (1) are where is the heat transfer coefficient and is the temperature of the wet cloth. are the dimensions of the body part which is beneath tepid sponging and is the surface area of the skin under tepid sponging.
In the biological system of human beings, the heat not only transfers through conduction but also by perfusion through blood circulation; see, for example, Aijaz et al. . Further, the heat transfer in humans also depends upon the metabolic heat generation and the impact of the external heat source as well as ambient temperature. The heat loss due to tepid sponging of the fevered body involves the heat transfer through human tissue; therefore, the bioheat equation is the most appropriate equation for the formulation of the model. The cubical domain of the model is discretized into small elements having meshes and nodes. This paper studies the effects of a tepid sponge on a small cuboidal portion of the forehead of the human body. Therefore, this assumption does not depart the model from the actual situation. The values of the temperature at the nodes on the boundary of the domain are known. To get the needed results and make the numerical calculations simple, finite difference methods (forward/backward/central difference approximations) are the most appropriate methods. However, in this paper, we shall opt for the forward difference method as it is frequently used in numerical methods and hence simple. For more generalized results, Equation (1) can be solved for both steady-state case and transient case.
3. Steady-State Case
The steady-state case of Equation (1) reduces to
For numerical approximations, the unit of discretization along all the axes is the same and is denoted by (numerically equal to 2 mm). So, the domain discretizes into 64 small cubes each of size , as shown in Figure 2. Thus, the forward difference approximation at the point is
The outer skin is in touch with the wet cloth. Therefore, the heat transfer coefficient is incorporated at the outer layer of the skin. We consider the cuboidal domain; hence, the temperature gradient at the five boundaries (excluding plane, i.e., where ) can be taken as zero. Thus, the boundary conditions are
Further, assuming that the wet cloth is dipped in cold water making the temperature of cloth 7°C uniformly, i.e., . Fever is generally said to occur when the internal body temperature exceeds 38°C. Therefore, we take the situations when the internal body temperature due to fever is 38°C, 39.5°C, and 41°C for the three different cases as shown in Tables 1–3, respectively, where
Discretizing the domain into four subintervals along the -axis, resulting in three unknown points on the -axis, so that, for three-dimensional discretization, there are unknown points. Giving in (4), we obtain 27 linear equations for 27 points of with matrix form of these equations as where
The constants and , , are defined as
MATLAB software has been used for the numerical solutions of . Then, for more accurate values, the domain discretization has been refined and numerical values have been enlisted in Tables 1–3 and displayed in Figures 3–5.
4. Transient Case
The corresponding numerical approximation of Equation (1) is based on the forward difference approximation for the time derivative and centered difference approximation for the space derivatives applied at the point . Thus, by FTCS (forward time centered space) method, we have where (numerically equal to 2 mm) is the uniform length of elements between two adjacent nodes and (numerically equal to 2 minutes) is the time interval between two consecutive readings.
The boundary conditions are the same as discussed in the steady-state case (5), and the initial conditions are
Furthermore, the various values of fever temperature for the transient case are the same as given by (6). The temperature of the surface of wet cloth covering the skin of the forehead is reset to 7°C after a time interval of 2 minutes, and the results are recomputed for numerical approximations of the temperature profiles as shown in Tables 4–6.
5. Domain Discretization and Numerical Scheme
The domain considered for tepid sponging is of the cuboidal shape having its outer skin surface as plane, and the -axis is in the normal direction to the skin surface. Excluding boundaries, three equidistant points , , and are taken along the -axis in the domain. Similarly, take the equally spaced points , , and along the -axis and , , and along the -axis inside the cuboidal domain as shown in Figure 2. For austerity, the three-dimensional frame is chosen in such a way that the point represents origin, i.e., . The domain so formed has known values of temperatures at peripheral boundaries. The forward difference approximation (in steady-state case) and forward time centered space (in transient case) are employed to find the numerical values of the temperatures at the different mesh points of the gridded domain. The discretization of the domain is based on the variation of numerical values of the physiological parameters of the forehead as elucidated by Ley and Bayazitoglu .
6. Results and Discussions
A mathematical model based on the bioheat equation has been formulated to study the temperature distribution and heat transfer in the scalp, skull, and CSF of the forehead during tepid sponging treatment to fever. The bioheat equation is preferred because it incorporates those parameters which make a significant contribution to heat transfer. The model has been solved for the numerical results by the forward difference method for the steady-state case and by the forward time centered space for the transient case separately. To display the results graphically, it becomes indispensable to get the numerical approximations by using the numerical values of the parameters. The numerical values of all the associated parameters (see Table 7) have been taken from the various research papers which include Khanday et al. , Janssen et al. , Aijaz et al. , Shirkavand and Nazif , Melo et al. , de Dear et al. , and Giering et al. . It is to be noted that there is a small variation in the values of the thickness of scalp, skull, and CSF, depending on the age and geographical factors. In this model, we shall take the average thickness of these layers as , , and .
Temperature distribution in the domain which lies under the wet cloth during tepid sponging is computed numerically. The temperature distribution for the steady-state case is shown in Table 1 (at fever 38°C), Table 2 (at fever 39.5°C), and Table 3 (at fever 41°C). For the transient case, the temperature distribution in the domain has been computed at , 4 min, and 6 min as shown in Table 4 (at fever 38°C), Table 5 (at fever 39.5°C), and Table 6 (at fever 41°C). These results have also been displayed graphically in Figures 3–5 (steady-state case) and Figures 6–8 (transient case).
In Figure 9, a comparison of the temperature distribution obtained before and after tepid sponging treatment in the patient’s forehead has been pictured. Figure 10 shows the comparison of the temperature distribution, recorded before and after tepid sponging at , , and . Based on the results obtained for the heat transfer and the temperature distribution during tepid sponging, it is evident from the tables and graphs that tepid sponging has shown a positive effect on the reduction of fever. The fevered body temperature regulates to normal temperature in a few minutes giving the patient relief and comfort.
Tepid sponging is one of the old practices used for the treatment of fever. Just at the beginning of the tepid sponging, the fever at the outer surface seems to have been normalized, but the internal parts of the body and the exterior layers of the brain are still in fever. Therefore, for diagnosing the effects of tepid sponging treatment for reducing temperature due to fever, continuous monitoring of temperature distribution needs to be tracked by some means. The various processes like temperature regulation and mechanism of fever control are also playing their part in reducing fever. In the present paper, an effort has been made to show the temperature variation numerically in the various layers of the forehead during tepid sponging. Therefore, the present paper can be a good tool in analyzing the significance of tepid sponging treatment for fever. Hence, or otherwise, a combination therapy (tepid sponging and paracetamol) could be effective in the treatment of fever. From Figures 3–8, it can be concluded that the reduction of temperature due to tepid sponging at outer surface drops quickly and the process of heat transfer starts for the rest of the layers of the domain. Moreover, Figures 9 and 10 show a blueprint picture of the comparison (before and after tepid sponging) of the temperature distribution in the domain.
Among the various treatments for fever, tepid sponging is useful, result-oriented, readily available, and without side effects. Generally, tepid sponging works nicely for some time, and then, it needs to be combined with paracetamol. The present study helps in determining whether the treatment by the tepid sponge is sufficient or not. If not, after how much time, tepid sponging is to be discontinued or supplemented with paracetamol. The diffusion of fluids and the heat transfer through the skin (and hence the other regions under it) of the forehead is different from the skin at the other parts of the human body as estimated by Eman . Therefore, the results cannot be directly applied to other parts of the body. However, there is a scope for improving this paper for the other parts of the human body. The novelty of this paper over the experimental-based papers is that the needed results can be directly obtained from this model by giving the numerical values to the parameter used in the formulation of the model.
No data has been used in this paper. The tables are numerical calculations and can be calculated by the reader himself with the help of equations and numerical values as given by Table 7.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
The authors thank and extend the appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through the research groups program under grant number R.G.P.2/82/42.
M. Meremikwu and A. OyoIta, Physical Methods for Treating Fever in Children, no. article CD004264, The Cochrane Library, Oxford, 2003.
P. A. Mackowiak, S. Wasserman, and M. Levine, “A critical appraisal of F, the upper limit of the normal body temperature, and other legacies of Carl Reinhold August Wunderlich,” JAMA, vol. 268, no. 12, pp. 1578–1580, 1992.View at: Publisher Site | Google Scholar
S. F. Morrison and K. Nakamura, “Central neural pathways for thermoregulation,” Frontiers in Bioscience, vol. 16, no. 1, pp. 74–104, 2011.View at: Publisher Site | Google Scholar
Hendrawati and M. Elvira, “Effect of tepid sponge on changes in body temperature in children under five who have fever in Dr. Achmad Mochtar Bukittinggi Hospital,” Enfermería Clínica, vol. 29, no. S1, pp. 91–93, 2019.View at: Publisher Site | Google Scholar
N. C. Emeka, “Treatment received by under fives having fever before presenting at the childrens outpatient clinic of a tertiary health facility in Owerri, Nigeria,” Annals of African Medicine, vol. 4, pp. 68–71, 2005.View at: Google Scholar
L. Y. Luchakov and A. D. Nozdrachev, “Mechanism of heat transfer in different regions of human body,” Physiology, vol. 36, no. 1, pp. 53–57, 2009.View at: Publisher Site | Google Scholar
T. M. Aluka, U. E. Asibong, A. N. Gyuse, M. M. Meremikwu, A. Oyo-Ita, and N. E. Udonwa, “Comparison of cold water sponging and acetaminophen in control of fever among children attending a tertiary hospital in South Nigeria,” Journal of Family Medicine and Primary Care, vol. 2, no. 2, pp. 153–158, 2013.View at: Publisher Site | Google Scholar
M. M. Meremikwu, A. Oyo-Ita, and Cochrane Infectious Diseases Group, “Physical methods versus drug placebo or no treatment for managing fever in children,” Cochrane Database of Systematic Reviews, no. 2, article CD004264, 2003.View at: Publisher Site | Google Scholar
N. B. Agbolosu, L. E. Cuevas, P. Milligan, R. L. Broadhead, D. Brewster, and S. M. Graham, “Efficacy of tepid sponging versus paracetamol in reducing temperature in febrile children,” Annals of Tropical Paediatrics, vol. 17, no. 3, pp. 283–288, 1997.View at: Publisher Site | Google Scholar
A. L. Eman, “The skin tissue of the human head subjected to thermal diffusion,” Mathematical Problems in Engineering, vol. 2018, Article ID 8781950, 6 pages, 2018.View at: Publisher Site | Google Scholar
W. Shen, J. Zhang, and F. Yang, “Modeling and numerical simulation of bioheat transfer and biomechanics in soft tissue,” Mathematical and Computer Modelling, vol. 41, no. 11-12, pp. 1251–1265, 2005.View at: Publisher Site | Google Scholar
M. S. Ferreira and J. I. Yanagihara, “A transient three-dimensional heat transfer model of the human body,” International Communications in Heat and Mass Transfer, vol. 36, no. 7, pp. 718–724, 2009.View at: Publisher Site | Google Scholar
M. Aijaz, I. M. Almanjahie, and J. G. Dar, “Mathematical estimation of fluid concentration in human skin during water immersion,” Journal of Advanced Research, vol. 28, pp. 1–6, 2021.View at: Publisher Site | Google Scholar
M. A. Khanday, M. Aijaz, and A. Rafiq, “Numerical estimation of the fluid distribution pattern in human dermal regions with heterogeneous metabolic fluid generation,” Journal of Mechanics in Medicine and Biology, vol. 15, no. 1, article 1550001, 2015.View at: Publisher Site | Google Scholar
Y. Tsunetsugu and M. Sugiyama, “Heat transfer, physiological responses, and subjective perceptions during short contact time with wood or other materials,” Journal of Wood Science, vol. 67, no. 1, p. 27, 2021.View at: Publisher Site | Google Scholar
O. Ley and Y. Bayazitoglu, “Effect of physiology on the temperature distribution of a layered head with external convection,” International Journal of Heat and Mass Transfer, vol. 46, no. 17, pp. 3233–3241, 2003.View at: Publisher Site | Google Scholar
H. H. Pennes, “Analysis of tissue and arterial blood temperatures in the resting human forearm,” Journal of Applied Physiology, vol. 1, no. 2, pp. 93–122, 1948.View at: Publisher Site | Google Scholar
M. Aijaz, M. A. Khanday, and A. Rafiq, “Variational finite element approach to study the thermal stress in multi-layered human head,” International Journal of Biomathematics, vol. 7, no. 6, article 1450073, 2014.View at: Publisher Site | Google Scholar
E. Majchrzak and L. Turchan, “Numerical analysis of tissue heating using the bioheat transfer porous model,” Computer Assisted Methods in Engineering and Science, vol. 20, pp. 123–131, 2013.View at: Google Scholar
M. Aijaz, I. M. Almanjahie, and J. G. Dar, “Spherical simulation of temperature profiles in tumour tissue and its peripherals during targeted hyperthermia,” Chiang Mai Journal of Science, vol. 48, no. 4, pp. 1–12, 2021.View at: Google Scholar
E. M. Janssen, M. J. Leeuwen, and A. A. Steenhoven, “Numerical simulation of scalp cooling to prevent hemotherapy–induced alopecia,” Journal of Mechanical Engineering, vol. 51, pp. 386–390, 2005.View at: Google Scholar
A. Shirkavand and H. R. Nazif, “Numerical study on the effects of blood perfusion and body metabolism on the temperature profile of human forearm in hyperthermia conditions,” Journal of Thermal Biology, vol. 84, pp. 339–350, 2019.View at: Publisher Site | Google Scholar
A. R. Melo, M. M. S. Loureiro, and F. Loureiro, “Blood perfusion parameter estimation in tumors by means of a genetic algorithm,” Procedia Computer Science, vol. 108, pp. 1384–1393, 2017.View at: Publisher Site | Google Scholar
R. J. de Dear, E. Arens, Z. Hui, and M. Oguro, “Convective and radiative heat transfer coefficients for individual human body segments,” International Journal of Biometeorology, vol. 40, no. 3, pp. 141–156, 1997.View at: Publisher Site | Google Scholar
K. Giering, I. Lamprecht, and O. Minet, “Specific heat capacities of human and animal tissues,” in Laser-Tissue Interaction and Tissue Optics, International Society for Optics and Photonics, 1996.View at: Google Scholar