Convective Heat Transfer Coefficients of the Human Body under Forced Convection from Ceiling

The average convective heat transfer coefficient for a seated human body exposed to downward flow from above was determined. Thermal manikin with complex body shape and size of an average Scandinavian female was used. The surface temperature distribution of the manikin’s body was as the skin temperature distribution of an average person. The measurements were performed in a room with controlled thermal environment. Air temperature was set at 26ºC for cooling and at 20ºC for heating. The radiant temperature asymmetry in horizontal and vertical direction was close to zero, i.e. mean radiant temperature was equal to the air temperature. The air velocity of the isothermal downward flow from the ceiling at height of 1.5 m above the floor (above the top of the head) was set in a range between still air and 0.73 m/s. Based on the analyses of the results relationships for determination of the convective heat transfer coefficient of the whole body (h c [W/(m 2 •K)]) was proposed: h c =4.088+6.592V 1.715 for a seated naked body at 20ºC and h c =2.874+7.427V 1.345 for a seated naked body at 26ºC. Differences in the convective heat transfer coefficient of the whole body in low air velocity range, V<0.3 m/s, due to the natural convection were found. The results may be useful during design of air distribution in rooms, e.g. low impulse ventilation, diffuse ventilation, etc.


Introduction
Air conditioning in office spaces is conventionally based on supply of air from the ceiling.Often occupants' comfort (thermal comfort and inhaled air quality) depends on the interaction of the supplied downward ventilation flow with the thermal plume generated by human body.This interaction is essential for the heat transfer from the human body and thus for its cooling.Convective cooling of the body at warm environment improves peoples' thermal comfort and may cause thermal discomfort due to draught at comfortable or cooler environment.
Draft is defined as local cooling of the body due to air movement.Studies that consider the issues of draft include those by Houghten et al. [1], Fanger et al. [2,3], McIntyre [4], Mayer [5], Berglund and Fobelets [6], Melikov [7] and Lee et al. [8].The studies by Houghten et al. [1], Fanger et al. [2,3], McIntyre [4], Berglund and Fobelets [6], Melikov [7] and Lee et al. [8] target forced convection by confronting the airflow from the frontal coronal plane of the body, and do not target downdraft from the upper horizontal plane of the body.On the other hand, Mayer [5] considers the head as the subject but does not consider the influence of other parts of the body.
When taking the air conditioners from the ceiling into consideration for offices, it is essential to have a grasp on the heat balance of the human body when planning, designing or evaluating the thermal environment.The usual angle factor of the human body and radiative heat transfer area is applicable for thermal radiation.In office spaces, the conductive heat transfer area is small.Moreover, the difference between the floor surface temperature and the skin temperature of the body is little.Therefore, It is possible to neglect the heat conduction [9,10].However, the direction of airflow from the ceiling is different from previous studies.The previous ones are the direction of the confronting airflow from the body's frontal coronal plane.Therefore, it is indispensable to determine whether the conventional research findings are applicable.
Previous researches on convective heat transfer from human body have focused on experimental rather than theoretical approach.In some studies simple heat emitting devices or assemblies has been used to model the convective heat transfer coefficient of the human body, other studies have used real subjects or thermal mannequins to investigate the heat balance of the human body.The mass and heat transport properties of naphthalene have been used as well.Most of these experimental studies on forced convections use wind tunnels or artificial climate chambers, and focus on the confronted airflow from the frontal coronal plane of the body.Under some conditions natural convection of an updraft from the lower horizontal plane of the body exists, but its characteristics differ from forced convection.
Research focused on the impact of direction of airflow on the heat transfer from human body has been reported [11][12][13][14][15]. Colin and Houdas [11] examined the convective heat transfer coefficient of a human body to airflow from the foot direction of the body to the horizontal plane, and the lateral direction of the body to the sagittal plane, for a body in a dorsal position.de Dear et al. [12] considered the convective heat transfer coefficient of human body in standing and seated posture to airflow from 8 horizontal directions, focusing on the frontal coronal plane of the body.Kurazumi et al. [15] determined the convective heat transfer coefficient of the human body to airflow from the front and back of the coronal plane of the body in a sitting with legs out position, and seated position.Mayer and Schwab [13,14] focused on the head area of the body, and examined the convective heat transfer coefficient of the head to the airflow from above, bellow, front and back.However, the convective heat transfer coefficient of the whole body has not been clarified.de Dear et al. [12] and Kurazumi et al. [15] show that the convective heat transfer coefficient of the human body is significantly different depending on the direction of the airflow.
In the case of air conditioning facilities in the ceiling surface of an office space, the direction of airflow is in downward, vertical to the main axis of the human body in a standing or seated position.However, previous research does not cover the convective heat transfer coefficient of the whole body with respect to the orientation of airflow at all.The purpose of this study was to examine the convective heat transfer by forced convection from seated human body exposed to airflow from above.

Calculation of Convective Heat Transfer Coefficient
In previous studies, the total amount of heat transfer was obtained by using physiological measurements and a heat flow meter.The convective heat transfer was determined by subtracting the radiant heat transfer from the total heat transfer.The angle factor value of the human body is essential to calculate the radiation heat transfer.However, the angle factor value between the measurement unit corresponding to the body section and enclosure has not been determined.Moreover, the appropriateness of the estimated value is inaccurate.
The radiative heat transfer coefficient of the whole body can be calculated in the process of determining the convective heat transfer coefficient or heat balance of the human body, but this specific value is rarely mentioned.Stolwijk [16], de Dear et al. [12], Ichihara et al. [17], Ishii [18] and Kurazumi et al. [15] have calculated the specific value.Values provided by Stolwijk [16], de Dear et al. [12], Ichihara et al. [17] and Kurazumi et al. [15] are for a human body model, and values by Ishii and Kurazumi et al. are for a human body.Studies by Stolwijk [16], de Dear et al. [12], Ichihara et al. [17] and Ishii [18] do not directly measure the amount of radiant heat transfer, but determine the angle factor of the space between the surface and each part, and conduct calculations through the Stefan-Boltzmann law.On the other hand, the studies carried out by Kurazumi et al. are on directly measured values [15,[19][20][21].These studies produce results that strongly illustrate the influence of radiative heat transfer area and body surface distribution.
In this paper, the total sensible heat transfer coefficient is obtained through experiments in a thermal environment in which the air temperature and each surrounding wall temperature are equal.The convective heat transfer coefficient, hc can be calculated by using the total amount of sensible heat transfer Q from the body.

Experimental Design
The experiment was carried out in an experimental test room (4.71 m×5.33 m×2.60 m) located in a large room (11.84 m×9.73 m×10.13 m).The test room was isolated and placed 0.75 m above the floor.The air temperature in the climatic room and experimental test room was kept the same, i.e. mean radiant temperature was the same as the air temperature and radiant temperature asymmetry in horizontal and vertical direction was the same.Air was supplied isothermally downward from the ceiling.The supplied airflow was with low initial turbulence intensity.The size of the air supply device, leveled with the ceiling, was 1.08 m×1.00 m.
Human body was resembled by a thermal mannequin with a body surface area of 1.4794 m 2 .The body of the thermal mannequin was divided into 23 segments: top of head and face (left and right), back of head and neck, chest (left and right), back, abdomen, buttocks, upper arm (left and right), forearm (left and right), hands (left and right), anterior thigh (left and right), posterior thigh (left and right), lower leg (left and right), and feet (left and right), making it possible to control the calorific value and surface temperature of each region.Table 1 shows the surface area ratio of each region of the thermal mannequin.The emissivity of the surface of the thermal mannequin is 0.95.

Regin Area ratio [-]
Parietal The mannequin measurements were taken in a seated position.The thermal mannequin was seated in a support base made with support rods that only touch the mannequin in 3 locations, and the experiment was carried out in conditions where the heat transfer amount between the thermal mannequin and the environment could be ignored.Therefore, all body segments were exposed to the airflow.The thermal mannequin was set up so that the soles of the feet were 0.10 m above the floor.The thermal mannequin and thermal environment were kept in a steady state.
Measurements were performed at two constant levels of air temperature: 26ºC (assuming cooling) and 20ºC (assuming heating period).The air velocity, determined from preliminary measurements, was set at six levels: still air, 0.2, 0.3, 0.4, 0.5, 0.7 m/s, where is near the top of head at 1.5 m above the floor.Relative humidity was the same for all conditions.
The air temperature and humidity, vertical air temperature distribution, air velocity, and surface temperature of each side of the room configuration were measured for the thermal environmental conditions.A thermistor thermometer was used to measure the air temperature and humidity at 0.93m vertically above the floor at the chestlevel of the mannequin, and the vertical air temperature distribution and surface temperature of each side of the room configuration was measured with a thermistor thermometer at 10 seconds intervals.The vertical air temperature distribution was measured at the soles (0.12 m), lower leg (0.53 m), chest (0.93 m), top of head (1.55 m), and at 2.00, 2.15 and 2.29 m above the floor.The downdraft air velocity sensor itself is a disturbance factor for the distribution of air velocity in the intervening areas of the body.Therefore, prior to the start of the experiments, the air velocity 1.50 m above the floor (corresponding to the top of the head level) was measured using an omni-directional anemometer (measurement accuracy ± 0.03 m/s).The initial velocity field at the air supply device was identified by measurements at 64 points positioned at 0.15 m intervals.Table 2 shows the air velocity and the measured results.During the experiment, the air velocity was controlled by monitoring the amount of air supplied to the room.
The surface temperature and heat supply to the 23 segments of the thermal mannequin was measured in 10 seconds intervals.
Once thermal equilibrium was established between the thermal mannequin and the surrounding environment, the measurements began.Each region of the thermal mannequin's electric energy was controlled to maintain each region's surface temperature at a constant level of 34ºC.
After confirming that the set thermal environment conditions and the heat loss from the thermal mannequin were in a steady state, the amount of heat supplied to the thermal mannequin and the surface temperature of the thermal mannequin were recorded.The thermal environment conditions 5 minutes before the end of the experiment and the surface temperature and amount of heat supplied to the thermal mannequin were used for the analysis.

Results and Discussion
The measured results of the thermal environment are shown in Table 3.The air temperature in setting conditions was kept in the range of ± 0.2ºC from the set value.The difference between the mean radiant temperature and the air temperature was about 0.1ºC.The calculation of the mean radiant temperature was performed using Fanger et al.'s angle factor value of the human body [22].The vertical air temperature distribution measured at the vicinity of the thermal mannequin identified almost constant temperature with differences in the range of ± 0.1ºC.The supplied airflow rate was stable throughout the experiment.Relative humidity was at around 15%. Therefore the thermal environment during each experiment nearly satisfied all setting conditions.Since the vertical air temperature distribution in the range of the thermal mannequin differed very little (in the range of ± 0.1ºC), and the difference between the mean radiant temperature and the air temperature was around 0.1ºC, it can be concluded that the thermal environment was nearly homogeneous.
The measured results of the thermal mannequin's surface temperature are shown in Table 3.It was determined to be around 34ºC for all conditions.Therefore the thermal mannequin's environment satisfied the setting conditions.
The convective heat transfer coefficient of the whole body of the thermal mannequin as a function of the air velocity shown in Figure 1   a larger convective heat transfer coefficient can be seen at a temperature of 20ºC when compared to 26ºC.Remarkable difference in the heat transfer coefficient at the two temperatures can be seen for velocities up to 0.3 m/s.The buoyancy flow generated by natural convection around mannequin's body, (due to the temperature difference of the thermal mannequin's surface temperature and the air temperature as a driving force) can be thought of as being more dominant than the downdraft from the ceiling.
Colin et al. [11] establish the air velocity of confronted airflow from the frontal coronal plane of a body at 0.15 m/s or less as a natural convection region.Rapp [23] defines the confronted airflow from the frontal coronal plane at 0.05-0.2m/s as a combined convection region, 0.05 m/s or less as a natural convection region, and 0.2 m/s or more as a forced convection region.Mitchell [24] defines the air velocity of confronted airflow from the frontal coronal plane of a body at 0.2 m/s or less as the natural convection region.Mayer and Schwab [13] have clarified the convective heat transfer coefficient for the head in the case of downward flow.At a head surface temperature of 34ºC and air temperature of 23ºC, downward airflow with velocity of up to 0.3 m/s shows a strong impact on the natural convection and thus on the convective heat transfer coefficient.Lee et al. [8] detected that in forced convection by confronted airflow from the frontal coronal plane of the body at less than 0.2 m/s, the updraft from natural convection makes an impact on the convective heat transfer coefficient.Mochida et al. [25] stated that due to a large variation in measured data from experiments using human subjects, separating the effect of natural convection in the low air velocity range of confronted airflow from the frontal coronal plane of the body is difficult.Judging from the experimental results, the air velocity is inferred 0.3 m/s.
In the present experiment downward flow with velocity of 0.3 m/s and higher became dominant over the upward thermal flow generated over the head of the seated heated mannequin and affected the relationships between the convective heat transfer coefficient and air velocity.Therefore, the interaction between the downward flows with the upward thermal flow is rather complex and different than the interaction of a horizontal flow with an upward thermal flow.In addition, in order to isolate and capture the effect of natural convection from the convective heat transfer coefficient by downdraft from the ceiling is extremely difficult.
The results of this study show that for downward flows with velocity 0.3 m/s and higher the same regression equation can be used regardless of the surrounding air temperature, i.e. 20ºC or 26ºC.The regression equation for the convective heat transfer coefficient (h c [W/(m 2 •K)]) for the whole body in forced convection is shown below.
20ºC nude seated position downward flow: h c = 4.088+6.592V 1.715  (4) 26ºC nude seated position downward flow: h c = 2.874+7.427V 1.345  (5) Here, It is difficult to isolate the effect of natural convection from the convective heat transfer coefficient by downdraft from the ceiling, but the regression line slope of the convective heat transfer coefficient at a velocity of 0.2 m/s or less is not strongly affected by air velocity.There are differences in the average values of convective heat transfer coefficient for different temperature conditions.Kurazumi et al. [21] proposed a formula for convective heat transfer coefficient due to natural convection while seated.If the thermal mannequin's surface temperature and the air temperature differ by 8ºC or 14ºC, the calculated convective heat transfer coefficient changes by about 20%.The convective heat transfer coefficient at an air velocity of below 0.2m/s at air temperatures of 20ºC and 26ºC can be expressed respectively at different air temperature conditions at mean values of 4.29W/(m 2 •K) and 3.04W/(m 2 •K).
Figure 2 compares the convective heat transfer coefficient of the seated nude body due to downward airflow from the ceiling with the convective heat transfer coefficient of the seated nude body.Since no comparative research has been conducted on the same direction of airflow a strict comparison cannot be done, but the dependent tendency of the convective heat transfer coefficient at an air velocity of 0.2 m/s or more to the air velocity can be regarded similarly in all experiments.Air temperature conditions at 20ºC where the temperature difference between the skin temperature and air temperature are nearly equal, show very similar calculated values for forced convection regions over 0.2 m/s as studies conducted by Mochida et al. [25].In addition, air temperature conditions at 26ºC show nearly the same calculated value as those of Mitchell [24].
At air velocities of less than 0.2 m/s, where the effect of natural convection on the convective heat transfer coefficient becomes strong, the convective heat transfer coefficient is strongly affected by the temperature difference between the skin temperature and air temperature.Kurazumi et al. [21] focused on the temperature difference between the skin temperature and air temperature during natural convection, and have proposed a formula for the convective heat transfer coefficient.Comparing the results calculated for different

Conclusion
The convective heat transfer coefficient of the whole seated body during forced convection by downward flow from the ceiling has been revealed through measurements using a thermal mannequin.At flow velocity over 0.3 m/s, the relationship between the convective heat transfer coefficient and the air velocity does not depend on the surrounding temperature, while at velocity up to 0.2 m/s, the free convection flow plays the major role for the heat exchange.The following empirical formula is suggested for calculation of the convective heat transfer coefficient of the whole body (h c [W/(m 2 •K)]).

Volume 4 •
Issue 1 • 1000126 J Ergonomics ISSN: 2165-7556 JER, an open access journal sensible heat transfer [W/m 2 ] C: Convective heat transfer [W/m 2 ] h c : Convective heat transfer coefficient [W/(m 2 •K)] T s : Surface temperature of the thermal mannequin [K] T a : Surrounding air temperature [K] T r : Mean radiant temperature of thermal mannequin [K] ε h : Emissivity of thermal mannequin surface [ND] ε w : Emissivity of the interior surface configuration [ND] σ: Stefan-Boltzmann constant [W/(m 2 •K 4 )] F w-h : Angle factor between thermal mannequin and room configuration [ND] f conv : Convective heat transfer area ratio of thermal mannequin [ND] f rad : Radiative heat transfer area ratio of thermal mannequin [ND].
. A trend for Ta is air temperature [ºC].Air velocity is mean ± standard deviation [m/s]

Figure 2 :
Figure 2: Comparison of convective heat transfer coefficients of the human body in previous studies.hc: convective heat transfer coefficient [W/(m 2 •ºC)].V: air velocity [m/s].∆T: difference between mean skin temperature and air temperature [ºC].

Table 1 :
Surface area ratios of segments of thermal mannequin.

Table 2 :
Measured values of air velocity.
Ta is air temperature[ºC].Tr is mean radiant temperature[ºC].∆Th is deviation of vertical air temperature in the vicinity of thermal mannequin[ºC].RH is relative humidty [%].Ts is mean surface temperature of thermal manneqin[ºC]

Table 3 :
Results of thermal conditions.