Research Papers

50 Years of Computer Simulation of the Human Thermoregulatory System

[+] Author and Article Information
George M. Netscher

Department of Biomedical Engineering,
The University of Texas at Austin,
Austin, TX 78712

Eugene H. Wissler

Department of Chemical Engineering,
The University of Texas at Austin,
Austin, TX 78712

Kenneth R. Diller

Department of Biomedical Engineering,
The University of Texas at Austin,
Austin, TX 78712
e-mail: kdiller@mail.utexas.edu

1Corresponding author.

Contributed by the Bioengineering Division of ASME for publication in the JOURNAL OF BIOMECHANICAL ENGINEERING. Manuscript received November 10, 2012; final manuscript received December 25, 2012; accepted manuscript posted January 18, 2013; published online February 7, 2013. Editor: Victor H. Barocas.

J Biomech Eng 135(2), 021006 (Feb 07, 2013) (9 pages) Paper No: BIO-12-1548; doi: 10.1115/1.4023383 History: Received November 10, 2012; Revised December 25, 2012; Accepted January 18, 2013

This paper presents an updated and augmented version of the Wissler human thermoregulation model that has been developed continuously over the past 50 years. The existing Fortran code is translated into C with extensive embedded commentary. A graphical user interface (GUI) has been developed in Python to facilitate convenient user designation of input and output variables and formatting of data presentation. Use of the code with the GUI is described and demonstrated. New physiological elements were added to the model to represent the hands and feet, including the unique vascular structures adapted for heat transfer associated with glabrous skin. The heat transfer function and efficacy of glabrous skin is unique within the entire body based on the capacity for a very high rate of blood perfusion and the novel capability for dynamic regulation of blood flow. The model was applied to quantify the absolute and relative contributions of glabrous skin flow to thermoregulation for varying levels of blood perfusion. The model also was used to demonstrate how the unique features of glabrous skin blood flow may be recruited to implement thermal therapeutic procedures. We have developed proprietary methods to manipulate the control of glabrous skin blood flow in conjunction with therapeutic devices and simulated the effect of these methods with the model.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Wissler, E. H., 1961, “Steady-State Temperature Distribution in Man,” J. Appl. Physiol., 16(4), pp. 734–740. [PubMed]
Wissler, E. H., 1964, “A Mathematical Model of the Human Thermal System,” Bull. Math. Biophys., 26, pp. 147–158. [CrossRef] [PubMed]
Wissler, E. H., 1970, “Comparison of Results Obtained from Two Mathematical Models – A Simple 14-Node Model and a Complex 250-Node Model,” J. Physiol. (Paris), 63, pp. 455–458.
Wissler, E. H., 1985, “Mathematical Simulation of Human Thermal Behavior Using Whole Body Models,” in Heat Transfer in Medicine and Biology, A.Shitzer and R. C.Eberhart, eds., Plenum Press, New York, Vol. 1, pp. 325–373.
Wissler, E., 1986, “Simulation of Fluid-Cooled or Heated Garments that Allow Man to Function in Hostile Environments,” Chem. Eng. Sci., 41, pp. 1689–1698. [CrossRef]
Nyberg, K. L., Diller, K. R., and Wissler, E. H., 2000, “Automatic Control of Thermal Neutrality for Space Suit Applications Using a Liquid Cooling Garment,” Aviat. Space Environ. Med., 71, pp. 904–915. [PubMed]
Nyberg, K. L., Diller, K. R., and Wissler, E. H., 2001, “Model of Human/Liquid Cooling Garment Interaction for Space Suit Automatic Thermal Control,” J. Biomech. Eng., 123, pp. 114–120. [CrossRef] [PubMed]
Wissler, E. H., 2003, “Probability of Survival during Accidental Immersion in Cold Water,” Aviat. Space Environ. Med., 74, pp. 47–55. [PubMed]
Pavlinic, D. Z., Wissler, E., and Mekjavic, I. B., 2009, “Modeling Thermophysiological Responses of Military Personnel Conducting a Variety of Activities during Simulated Field Operations in a Cold Environment,” NATO OTAN, 28, pp. 1–19.
Wissler, E. H., 2009, “Whole-Body Human Thermal Models,” Adv. Numer. Heat Transf., 3, pp. 306–344.
Hensley, D. W., Mark, A. E., Wissler, E. H., and Diller, K. R., 2012, “Quantitative Analysis of Glabrous Skin Blood Flow and its Role in Human Thermoregulation,” ASME SBCSBC2012-80854, pp. 1–2.
Wissler, E. H., 2012, “Whole-Body Human Thermal Modeling, an Alternative to Immersion in Cold Water and Other Unpleasant Endeavors,” J. Heat Transfer, 134, 031019, pp. 1–11. [CrossRef]
Stolwijk, J. A. J., and Hardy, J. D., 1966, “Temperature Regulation in Man: A Theoretical Study,” Pfluegers Arch., 291, pp. 126–162.
Gagge, P. A., Stolwijk, J. A. J., and Nishi, Y., 1970, “An Effective Temperature Scale Based on a Simple Model of Human Physiological Regulatory Response,” ASHRAE Trans., 70, pp. 247–262.
Kraning, K. K., and Gonzalez, R. R., 1977, “A Mechanistic Computer Simulation of Human Work in Heat That Accounts for Physical and Physiological Effects of Clothing, Aerobic Fitness, and Progressive Dehydration,” J. Therm. Biol., 22, pp. 331–342. [CrossRef]
Fiala, D., Havenith, G., Brode, P., Kampmann, B., and Jendritzky, G., 2012, “UTCI-Fiala Multi-Node Model of Human Heat Transfer and Temperature Regulation,” Int. J. Biometeorol., 56, pp. 429–441. [CrossRef] [PubMed]
Xu, Q., and Luo, X., 2006, “Dynamic Thermal Comfort Numerical Simulation Model on 3D Garment CAD,” Appl. Math. Comput., 182, pp. 106–118. [CrossRef]
Gao, N., Niu, J., and Zhang, H., 2006, “Coupling CFD and Human Body Thermo-Regulation Model for the Assessment of Personalized Ventilation,” HVAC&R Res., 12, pp. 497–518. [CrossRef]
Bazett, H. C., and McGlone, B., 1927, “Temperature Gradient in the Tissues in Man,” Am. J. Physiol., 82, pp. 415–451.
Burton, A. C., 1934, “The Application of the Theory of Heat Flow to the Study of Energy Metabolism,” J. Nutrition, 7, pp. 497–533.
Pennes, H. H., 1948, “Analysis of Tissue and Arterial Blood Temperature in the Resting Human Forearm,” J. Appl. Physiol., 1, pp. 93–122.
Stolwijk, J. A. J., 1970, “A Mathematical Model of Physiological Temperature Regulation in Man,” NASA-9-9531.
Wissler, E. H., 1998, “Pennes' 1948 Paper Revisited,” J. Appl. Physiol., 85, pp. 35–41.
Berlinerblau, F., 1875, “Ueber den direkten Uebergang von Arterien in Venen,” Arch. Anat. Physiol. (cited and discussed in [26]).
Hoyer, H., 1877, “Ueber unmittelbare Einmeundung kleinstrer Arterien in Gefassaste venosen Charackters,” Arch. Mikr. Anat., 13, p. 603 (cited and discussed in [21]). [CrossRef]
Clark, E. R., 1938, “Arterio-Venous Anastomoses,” Physiol. Rev., 18, pp. 229–245.
Sherman, J. L., Jr., 1963, “Normal Arteriovenous Anastomoses,” Medicine, 42, pp. 247–268. [CrossRef]
Zanick, D. C., and Delaney, J. P., 1973 “Temperature Influences on Arteriovenous Anastomoses,” Proc. Soc. Exp. Biol. Med., 144, pp. 616–620.
Charkoudian, N., 2003 “Skin Blood flow in Adult Human Thermoregulation: How It Works, When It Does Not, and Why,” Mayo Clin. Proc., 78, pp. 603–612. [CrossRef]
Grahn, D. A., Cao, V. H., and Heller, H. C., 2005, “Heat Extraction Through the Palm of One Hand Improves Aerobic Exercise Endurance in a Hot Environment,” J. Appl. Physiol., 99, pp. 972–978. [CrossRef]
Taylor, N. A. S., Machado-Moreira, C., van den Heuvel, A., Caldwell, J., Taylor, W. A., and Tipton, M. J., 2009, “The Roles of Hands and Feet in Temperature Regulation in Hot and Cold Environments,” 13th International Conference on Environmental Ergonomics, Boston, MA, pp. 405–409.
Bergerson, T. K., 1992, “A Search for Arteriovenous Anastomoses in Human Skin Using Ultrasound Doppler,” Acta Physiol. Scand., 147, pp. 195–201. [CrossRef]
Tikuisis, P., Meunier, P., and Jubenville, C. E., 2001, “Human Body Surface Area: Measurement and Prediction Using Three Dimensional Body Scans,” Eur. J. Appl. Physiol., 85, pp. 264–271. [CrossRef]
Lee, J.-Y., Choi, J.-W., and Kim, H., 2007, “Determination of Hand Surface Area by Sex and Body Shape using Alginate,” J. Physiol. Anthrop., 26, pp. 475–483. [CrossRef]
Tipton, M. J., Allsopp, A., Balmi, P. J., and House, J. R., 1993, “Hand Immersion as a Method of Cooling and Rewarming: A Short Review,” J. R.l Naval Med. Serv., 79(3), pp. 125–131.
Taylor, N. A. S., Caldwell, J. N., van den Heuvel, A. M. J., and Patterson, M. J., 2008, “To Cool, But Not Too Cool: That Is the Question – Immersion Cooling for Hyperthermia,” Med. Sci. Sports Exer., 40(11), pp. 1962–1969. [CrossRef]
Casa, D. J., and Kenny, G. P., 2009, “In Defense of Cold Water,” Med. Sci. Sports Exer., 41(5), p. 64. [CrossRef]
Taylor, N. A. S., Caldwell, J. N., van den Heuvel, A. M. J., and Patterson, M. J., 2009, “Response,” Med. Sci. Sports Exer., 41(5), p. 65. [CrossRef]
Taylor, N. A. S., Casa, D. J., and Kenny, G. P., 2010, “Immersion Treatment for Exertional Hyperthermia; Cold or Temperature Water? Contrasting Perspectives,” Med. Sci. Sports Exer., 42(5), pp. 1246–1252. [CrossRef]
Wilson, T. E., Zhang, R., Levine, B. D., and Crandall, C. G., 2005, “Dynamic Autoregulation of Cutaneous Circulation: Differential Control in Glabrous Versus Nonglabrous Skin,” Am. J. Physiol. Heart Circ. Physiol., 289, pp. H385–H391. [CrossRef]
Kellogg, D. L., 2006, “In Vivo Mechanisms of Cutaneous Vasodilation and Vasoconstriction in Humans During Thermoregulatory Challenges,” J. Appl. Physiol., 100, pp. 1709–1718. [CrossRef]
Johnson, J. M., and Proppe, D. W., 1996, “Cardiovascular Adjustments to Heat Stress,” in Handbook of Physiology, Sec 4: Environmental Physiology, Vol. 1, Oxford University Press, pp. 215–243.
Hales, J. R. S., and Molyneux, G. S., 1988, “Control of Cutaneous Arteriovenous Anastomoses,” in Vasodilation: Vascular Smooth Muscle, Peptides, Autonomic Nerves, and Endothelium, P. M.Vanhoutte, ed., Raven Press, New York, pp. 321–332.
Wissler, E. H., 2008, “A Quantitative Assessment of Skin Blood Flow in Humans,” Eur. J. Appl. Physiol., 104, pp. 145–157. [CrossRef]
Roddie, I. C., 1983, “Circulation to Skin and Adipose Tissue,” Handbook of Physiology, Sec 2 Cardiovas System, Vol III. Peripheral Circulation Organ Blood Flow, Part I, pp. 285–317.
Colemont, L., and Decoutere, P., 1981, “The In Vivo Determination of the Blood Volume in the Foot Using a Mercury Starin Gauge Plethysmograph and I-131 Labeled Serum Albumin,” Cardiovasc. Res., 15, pp. 262–267. [CrossRef]
Nagasaka, R., Cabanac, M., Hirata, K., and Nunomura, T., 1987, “Control of Local Heat Grain by Vasomotor Response of the Hand,” J. Appl. Physiol., 63(4), pp. 1335–1338.
Wenger, C. B., Roberts, M. F., Nadel, E. R., and Stolwijk, J. A. J., 1975, “Thermoregulatory Control of Finger Blood Flow,” J. Appl. Physiol., 38(6), pp. 1078–1082.
Grahn, D. A., Davidson, H. L., and Heller, H. C., 2011, “Incorporating Variable Vascular Heat Exchangers Into Models of Human Thermoregulation,” American Association on Artificial Intelligence Spring Symposium Series, pp. 1–3.
Fiala, D., Lomas, K. J., and Stohrer, M., 1999, “A Computer Model of Human Thermoregulation for a Wide Range of Environmental Conditions: The Passive System,” J. Appl. Physiol., 87, pp. 1957–1972.
Lee, Y., and Hwang, K., 2002, “Skin Thickness of Korean Adults,” Surg. Radiol. Anat., 24, pp. 183–189. [CrossRef]
Manelli, A., Sangiorgi, S., Ronga, M., Reguzzoni, M., Mini, A., and Raspanti, M., 2005, “Plexiform Vascular Structures in the Human Digital Dermal Layer: A SEM-Corrosion Casting Morphological Study,” Eur. J. Morph., 42, pp. 173–177. [CrossRef]
Bergersen, T. K., Eriksen, M., and Walloe, L., 1997, “Local Constriction of Arteriovenous Anastomoses in the Cooled Finger,” Am. J. Physiol., 273, pp. R880–R886.
Savastano, D. M., Gorback, A. M., Eden, H. S., Brady, S. M., Reynolds, J. C., and Yanovski, J. A., 2009, “Adiposity and Human Regional Body Temperature,” Am. J. Clin. Nutr., 90, pp. 1124–1131. [CrossRef]


Grahic Jump Location
Fig. 1

Graphical representation of the classic Wissler model that terminates at the wrists and ankles. Each of the 25 elements consists of 21 radial layers and 12 angular segments. The new model has four additional elements for two hands and two feet that are highlighted in red.

Grahic Jump Location
Fig. 2

Graphical representation of the control algorithm for normalized glabrous vasodilation (GVD) of Eqs. (1) and (2) that regulates glabrous skin blood flow (GBF). There is a region of intermediate level perfusion associated with thermoneutral control inputs. Increasing core temperature above this region causes a sigmoidal elevation in glabrous perfusion, and decreasing core and mean skin temperatures likewise depresses the perfusion. A minimum level of 5% nutritive blood flow is assumed to be always maintained.

Grahic Jump Location
Fig. 3

GUI screens for specifying input parameters for a thermoregulation process simulation. (a) Physical and physiological data for the subject. (b) AVA blood perfusion levels and control temperatures. (c) Physical activity sequence, types, and durations.

Grahic Jump Location
Fig. 4

Example GUI screen to input clothing type and thermal properties and boundary temperature for each of the 25 model elements

Grahic Jump Location
Fig. 5

Simulations by the Wissler and BME models for a male walking for 30 min following equilibration at rest in a thermoneutral environment. (a) Mean skin temperature; (b) skin temperature on thigh; (c) skin temperature on hand (BME model) and end of arm (Wissler model); (d) esophageal temperature.

Grahic Jump Location
Fig. 6

Effect on modulation of core temperature as a function of the magnitude of glabrous skin blood flow as a percent of cardiac output (CO) to both hands and both feet with surface cooling at 22  °C. The metabolic protocol consisted of during continuous exposure to air at 23  °C.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In