Natural convection heat transfer coefficient for newborn baby Thermal manikin assessed convective heat loses

The energy balance and heat exchange for newborn baby in radiant warmer environment are considered. The present study was performed to assess the body dry heat loss from an infant in radiant warmer, using copper cast anthropomorphic thermal manikin and controlled climate chamber laboratory setup. The total body dry heat losses were measured for varying manikin surface temperatures (nine levels between 32.5 ◦ C and 40.1 ◦ C) and ambient air temperatures (five levels between 23.5 ◦ C and 29.7 ◦ C). Radiant heat losses were estimated based on measured climate chamber wall temperatures. After subtracting radiant part, resulting convective heat loses were compared with computed ones (based on Nu correlations for common geometries). Simplified geometry of newborn baby was represented as: (a) single cylinder and (b) weighted sum of 5 cylinders and sphere. The predicted values are significantly overestimated relative to measured ones by: 28.8% (SD 23.5%) for (a) and 40.9% (SD 25.2%) for (b) . This showed that use of adopted general purpose correlations for approximation of convective heat losses of newborn baby can lead to substantial errors. Hence, new Nu number correlating equation is proposed. The mean error introduced by proposed correlation was reduced to 1.4% (SD 11.97%), i.e. no significant overestimation. The thermal manikin appears to provide a precise method for the noninvasive assessment of thermal conditions in neonatal care.


Introduction
The heat exchange with the environment is of greater importance in the infant than in the adult, explaining the increased risk of body hypo-or hyperthermia [2]. As result, maintaining fluid and heat balance is of vital importance to the newborn infant [3]. The infants are exposed to a cold and dry environment, and are then at risk of dehydration and hypothermia. These conditions may have serious consequences and significantly influence mortality and morbidity. Moreover in recent years there is a growing interest in therapeutic use of hypothermia for treatment of neonatal encephalopathy [4][5][6][7]. To assist both therapy design and control, attempts are made to use mathematical models of human thermal physiology (i.e. heat transfer and active thermoregulation) [8][9][10][11][12].
To take advantage of all above mentioned issues, including the therapeutic hypothermia modeling, it is important to understand and quantify the heat transfer processes between the baby and its environment whilst it is in the radiant warmer.
There are numerous studies on natural convection heat transfer coefficient in adults based on direct measurements and thermal manikins [13,14], as well as based on numerical simulations [15]. Review and comparison of convective heat transfer coefficients of the adult human body in several studies is presented in [16].
However, for a case of newborn babies, the available literature sources are very limited. The only available data are obtained for newborn thermal manikins in a single-walled, air-heated closed incubators. Early reports, like [17], were based on simplified geometry manikins, while recent studies of dry heat losses from the manikin [18,19] account for anthropomorphic premature newborn sized manikins. Yet, expect of initial measurements reported in [20], no detailed reports have been found for open type radiant warmers that are widely used in pediatric centers.
Due to ethical reasons, experiments cannot be performed on living infants, while they are at risk of body cooling and cannot tolerate the experimental procedures. In such situation, the idea of thermal manikins' usage for heat exchange identification arises in natural way. As result, within current research, the heated thermal manikin and controlled climate chamber setup is being proposed, designed and used to study dry heat loses in stagnant air for the infant laying in infant care bed conditions.

Thermal manikin
The anthropomorphic thermal manikin is used in the present study to represent newborn baby. Its size correspond to the dimensions of prematurely born baby (35 week of gestational age) in 50th percentile [21]. The manikin, shown in Fig. 1, was cast in copper and painted matt (graphite) black (see Fig. 6). The model was designed to mimic newborn baby with body surface area S = 0, 13278m 2 (measured by means of 3D manikin scan). The dimensions of manikin are presented in the Table 1.
To maintain and control the constant surface temperature of the manikin, the water heating system was proposed. Stabilized (±0.01 • C) temperature water (Advanced Digital Controller Refrigerated/Heated Circulating Baths, model AD07R-40, VWR International, USA) was supplied into the thermal manikin by means of rubber hoses. To obtain uniform temperature distribution over body surface, the water inlet was divided into 5 separate outflow ports located in limbs (four) and head (one). The outlet port was located on the back of the manikin (see Fig. 2). The mass flow rate of heating water was kept on elevated level, to ensure uniform manikin surface temperature distribution. Uniform temperature distribution was achieved, an exemplary IR thermogram (ThermaCam SC2000, FLIR, USA) of body surface temperature is shown in Fig. 3, mean value of surface temperature with one standard deviation (SD) is 36.42 • C (0.13 • C). Slightly elevated temperature was measured only in self irradiated areas of the manikin surface (and these locations were excluded form above mentioned mean temperature computation).
Personal computer (LabView Signal Express software and 24-Bit Universal Analog Input card NI9219 National Instruments, USA) was used to record the manikin surface temperatures and in/out water supply (thermocouples type K, diameter 0.5mm, CZAH-Pomiar, Poland). Manikin surface temperature was recorded by two contact thermocouples (the same type as above) attached to the manikin surface by adhesive tape and protected by aluminium foil coating to prevent its direct irradiation (see Fig. 4). First was located in anterior (front) trunk, the latter one on posterior (back) tight. The heating power delivered to the manikin  (i.e. dry heat loses) was measured by means of water mass flow rate measurement (rotameter type KM, Z.A. Rotametr sp. z o.o., Poland) and the two thermocouples installed in inlet and outlet water ports of manikin. Therefore the heating power (i.e. dry heat loses) can be calculated as: where: P -heating power (dry heat losses), W; mwater mass flow rate, kg·s −1 ; c p -water heat capacity, J·kg −1 ·K −1 ; T w,in , T w,out -thermal manikin inlet and outlet water temperatures, respectively, • C.

Climate chamber
The measurements were performed in the double wall (chamber-in-chamber type) climate chamber located in the Laboratory Hall of Department of Heating, Ventilation and Dust Removal Technology (Silesian University of Technology, Gliwice, Poland). The climate chamber was made of wooden oriented strand boards (OSB), having inside dimensions: length 334cm, width 232cm and height 220cm. The internal wooden chamber is located in slightly larger room, with thermal insulated walls. Central air conditioning system was used to deliver stabilized temperature air both to in-and outside zone of the chamber. The inside air inflow is provided by means of holes in ceiling panels made of stainless steel, while four air outflow ports are located in corners (floor level). Such setup allows to control precisely the climate chamber wall temperature (by means of circulating air on both sides of chamber wall). There were 45 temperature sensors (Programmable Resolution 1-Wire Digital Thermometer, type 18B20, Maxim Integrated, USA) installed inside the climate chamber to continuous monitor and record both wall and air temperatures. Sensors were located in 11 locations, in

Radiant warmer
The manikin was lying in a spread-supine (with hands heading up) position on a mattress in an infant care bed (Babytherm 8000, Dräger Medical, Germany). Built in radiant warmer was set off. The mattress was not heated, and in current research can be treated as thermal insulation. Under this assumption (i.e. no conductive heat losses on the manikin -mattress interface), the body surface area used in further computations should be corrected. The manikin surface being in contact with mattress has been marked and measured (3D scan). The effective body surface area was measured S e = 0.107m 2 .

Measurements
During each experimental session, steady state was strictly required. First, for given air temperature and lowest surface (manikin) temperature, the steady state was reached (after approximately 4 hours), then the air flow inside the climate chamber was turned off. When the heat transfer between thermal manikin and thermal environment reached equilibrium (after approx. 15 minutes) the chamber air, wall temperatures, manikin surface temperature & heat loses were recorded in steady state conditions. Then the surface Fig. 6 Thermal manikin covered with graphite (high emissivity) spray coating (manikin) temperature was increased up to next test case value, and after reaching equilibrium data were recorded again. Due to low dissipated power (range 4-16 W) and relatively large size of the chamber, the air temperature inside the chamber was constant during experiment, i.e. the overall air temperature increase was negligible in course of taking measurements.
The temperature measurements were recorded with 2Hz frequency for every test case in steady conditions (for approx. 10 minutes each, resulting in approx. 1200 samples for every test case). Mean values of recorded: manikin surface temperatures, air & wall temperatures and dry heat losses, calculated according to Eq. 1, are presented in Table 3. Low levels of standard deviations (SD) in recorded temperatures for all test cases confirm that the data were recorded in the steady state conditions regime, as planned.
Prior to the test sessions all temperature measurement devices (ie. thermocouples and digital sensors) has been calibrated in ice water (0 • C) as well as in dry-well calibration device (Transcat 5853T/model 9100S, Hart Scientific, USA). Device specific offsets has been set accordingly.
During data recording sessions there were almost no vertical temperature gradients inside the chamber (within ±0.1 • C). Air velocity was kept far below 0.1m/s.
To account for radiative heat loses, emissivity of all surfaces exchanging heat by radiation have to be determined. The thermal manikin surface has been covered with graphite high emissivity spray coating (see Fig. 6), with known emissivity of 0.99. Remaining surfaces data has obtained from manufacturers' datasheets, and are all listed in Table 2.

Results
The measurements were carried out under steady state conditions. The manikin was exposed to five different ambient air temperatures (between 23.5 • C and 29.7 • C) and nine  Table 3.

Dry heat loses
Similar to the observations reported in [18] for single walled, air-heated closed incubator, the dissipated power in current experiment increases linearly with the increasing excess temperature of the manikin surface over ambient air (cf. Fig. 7):

Convective heat loses
In order to estimate pure convective heat loses, the radiation part of the manikin heat balance should be estimated first. The radiative part of total dissipated heat flux (P r ) is computed as: where: S e -effective body surface area, m 2 , σ = 5.67×10 8 , W·m −2 ·K −4 -blackbody radiation constant; F ε -transfer factor (to account geometrical view factors for radiative heat exchange between manikin and surrounding), T s and T w are manikin surface and chamber wall temperatures, respectively, K. Convective only heat loses P c are then computed by subtracting radiative heat losses P r (2) from total dry heat losses P (1):

Comparison with available correlating equations
Measured convective heat loses P c (3) are compared with computed ones and are shown in Fig. 8. Computed values are based on general purpose (GP) Nusselt (Nu) correlating formulas for common geometries in free convection regime: horizontal cylinder [23] and sphere [24] (see Appendix).  Detailed measurement data are presented in Table 3 Two cases are considered for comparison, for which convective heat transfer coefficient is estimated by means of simplified geometry of newborn baby represented as: (a) single horizontal cylinder (with diameter as for trunk), (b) cylinders and sphere, where the weighted sum of 5 horizontal cylinders (trunk and limbs) and sphere (head) is used to compute the overall Nu number. Summation weights are taken according to a Rule of nines 1 , see Table 4 (the skin surface areas weights are modified in infants and children because of the proportionately larger head size).
It should be noted that the above approach assumes free (not disturbed) flow of air around modeled surfaces. In case of newborn baby in infant care bed, this assumption can not be fulfilled: newborn baby body is lying on the mattress (i.e. bottom part of the skin surface is insulated) and the presence of the infant care bed itself disturbs the air flow around the body. As result, it can be noticed that usage of classical correlating equations, valid for objects in open environment, leads to overestimation of dissipated energy. Statistical analysis t-test (MatLAB, The MathWorks, Inc., USA) shows that computed values of convective heat loses P c are significantly (p < 0.001) overestimated relative to measured values for both approaches, by: where values are presented as mean (with one standard deviation).
To overcome overestimation, the dissipated power due to convection can be approximated using two-term power series function of excess temperature of the manikin surface over ambient air (cf.

Proposed correlating equation
To generalize natural convection heat transfer coefficient estimation of the newborn baby in infant care bed the nondimensional analysis of experimental data (listed in Table 3) has been performed. Air properties has been determined at film temperature (i.e arithmetic mean of surface and free stream fluid temperatures): using CoolProp libraries [22]. The diameter of the trunk has been chosen as characteristic dimension. Based on experimental data (in range of 2.5 × 10 5 ≤ Ra D ≤ 10 6 ) a Nusselt number (Nu D ) correlating equation for natural convection from newborn baby in air has been worked out. The two-term power series data fit (r 2 = 0.4159) of Nu D = f (Ra D ) reads: Nu D = 9.179 + 1.043 × 10 14 Ra −2.452 D (5) and is shown in Fig. 9.

Comparison of correlating equations
The plot of proposed correlating equation is shown together with experimental data in Fig. 9. To maintain continuity, the computed values using GP equations (a) and (b) are shown as well.
For all test cases, the absolute error for Nu number computed using proposed correlating equation comes from a normal distribution and a zero mean (t-test, p = 0.994) with SD 1.2493.
In Fig. 10 the Nu numbers determined in course of experiment are presented together with data predicted using: the proposed correlation Eq. 5 and GP correlating equations, i.e.  The improvement in accuracy of convective heat losses prediction has been achieved. The computed heat loses P c using proposed correlating Eq. 5 is shown in Fig. 11. The mean relative error of is 1.4% (SD: 11.97%) and what should be noticed -in contrary to errors introduced by approaches (a) and (b) -it shows no significant overestimation (t-test, p = 0.718).

Summary
The anthropomorphic newborn baby thermal manikin was used to carry out measurements of the whole body dry heat loses under free convection regime in the radiant warmer-like environment. The experimental setup and proposed methodology proved to be effective way of estimating dry and convective heat loses for non-standard geometrical cases.
Presented results show that usage of general purpose Nu correlations leads to significant overestimation. This is due to the fact that general purpose formulas are valid for cylinders/spheres surrounded by air (with no obstructions to air flow), while in analyzed geometry, i.e. newborn baby laying on the mattress, part of the body surface is insulated and air flow is obstructed by presence of the bed itself.
Proposed correlating Eq. 5 can be used to model convective heat loses of newborn babies in real hospital conditions (radiant warmer environment) in range of ambient air temperatures 23 • C to 30 • C and skin temperatures 32 • C to 40 • C. The formula has been fitted to measurements of thermal manikin which size corresponds to the prematurely (35th week) born baby, in 50th percentile [21] (see Table  1). This correspond to the Raileigh (Ra D ) number range: 2.5 × 10 5 to 10 6 (trunk diameter taken as characteristic dimension). The extrapolation of the results on other newborn baby sizes has not been tested in current study.