Dominant factors affecting temperature rise in simulations of human thermoregulation during RF exposure

Numerical models of the human thermoregulatory system can be used together with realistic voxel models of the human anatomy to simulate the body temperature increases caused by the power absorption from radio-frequency electromagnetic fields. In this paper, the Pennes bioheat equation with a thermoregulatory model is used for calculating local peak temperatures as well as the body-core-temperature elevation in a realistic human body model for grounded plane-wave exposures at frequencies 39, 800 and 2400 MHz. The electromagnetic power loss is solved by the finite-difference time-domain (FDTD) method, and the discretized bioheat equation is solved by the geometric multigrid method. Human thermoregulatory models contain numerous thermophysiological and computational parameters--some of which may be subject to considerable uncertainty--that affect the simulated core and local temperature elevations. The goal of this paper is to find how greatly the computed temperature is influenced by changes in various modelling parameters, such as the skin blood flow rate, models for vasodilation and sweating, and clothing and air movement. The results show that the peak temperature rises are most strongly affected by the modelling of tissue blood flow and its temperature dependence, and mostly unaffected by the central control mechanism for vasodilation and sweating. Almost the opposite is true for the body-core-temperature rise, which is however typically greatly lower than the peak temperature rise. It also seems that ignoring the thermoregulation and the blood temperature increase is a good approximation when the local 10 g averaged specific absorption rate is smaller than 10 W kg(-1).