Next Article in Journal
Design of the MEMS Piezoresistive Electronic Heart Sound Sensor
Next Article in Special Issue
Smart Toys Designed for Detecting Developmental Delays
Previous Article in Journal
Simultaneous Life Detection and Localization Using a Wideband Chaotic Signal with an Embedded Tone
Previous Article in Special Issue
When Ultrasonic Sensors and Computer Vision Join Forces for Efficient Obstacle Detection and Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calorimetry Minisensor for the Localised Measurement of Surface Heat Dissipated from the Human Body

by
Fabiola Socorro
,
Pedro Jesús Rodríguez de Rivera
and
Manuel Rodríguez de Rivera
*
Departamento de Física, Universidad de Las Palmas de Gran Canaria, Las Palmas de Gran Canaria E-35017, Spain
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(11), 1864; https://doi.org/10.3390/s16111864
Submission received: 21 September 2016 / Revised: 28 October 2016 / Accepted: 3 November 2016 / Published: 6 November 2016
(This article belongs to the Special Issue Sensing Technology for Healthcare System)

Abstract

:
We have developed a calorimetry sensor that can perform a local measurement of the surface heat dissipated from the human body. The operating principle is based on the law of conductive heat transfer: heat dissipated by the human body passes across a thermopile located between the individual and a thermostat. Body heat power is calculated from the signals measured by the thermopile and the amount of power dissipated across the thermostat in order to maintain a constant temperature. The first prototype we built had a detection area measuring 6 × 6 cm2, while the second prototype, which is described herein, had a 2 × 2 cm2 detection area. This new design offers three advantages over the initial one: (1) greater resolution and three times greater thermal sensitivity; (2) a twice as fast response; and (3) it can take measurements from smaller areas of the body. The sensor has a 5 mW resolution, but the uncertainty is greater, up to 15 mW, due to the measurement and calculation procedure. The order of magnitude of measurements made in healthy subjects ranged from 60 to 300 mW at a thermostat temperature of 28 °C and an ambient room temperature of 21 °C. The values measured by the sensor depend on the ambient temperature and the thermostat’s temperature, while the power dissipated depends on the individual’s metabolism and any physical and/or emotional activity.

1. Introduction

The main developments in calorimetry instrumentation have occurred in the area of thermal analysis, a field that has observed the construction of a myriad of devices that are only differentiated by the specific characteristics of the thermal process being studied [1,2]. In calorimetry, in order to measure the heat energy produced during a process, it must take place in an enclosed, controlled area where the external factors influencing the thermodynamic process can be minimised. This work presents a calorimetric sensor designed for localised measurement of surface heat dissipated by the human body. We therefore attempt to measure the power transmitted by conduction using a thermopile placed between the human body and a thermostat maintained at a constant temperature. In this particular calorimetry application, the dissipation being measured does not occur inside the measuring device but rather outside it; this instrument does not, therefore, comply with normal calorimetry standards. The working principle for this small instrument is based on the law of conductive heat transfer and, as its operation requires a thermostat to reach a constant temperature, it can be classed as an isothermal heat conduction calorimeter [3].
The instrument has a wide variety of applications as it can be used to study different energetic processes occurring in the human body in terms of power and energy. One potential application is to control the power dissipated while using different methods to burn malignant cells. It also offers the opportunity to study dissipation from different areas of the body, to identify abnormal dissipation, or to examine body heat dissipation according to metabolism, and physical and/or emotional activity. We have already presented our first prototype with a 6 × 6 cm2 detection area in previous studies [4,5,6,7]. In this work, we present a second, smaller prototype with a detection area measuring just 2 × 2 cm2. The advantages of this second prototype are: (1) greater resolution; (2) quicker response; and (3) it can take measurements from smaller areas of the body.
Body heat measurements are currently taken for the entire body and primarily through indirect calorimetry techniques based on CO2 production (VCO2) or O2 consumption (VO2) measurements. These measurements and some empirical formulae are used to calculate the total heat dissipated by the subject in function of the physical activity being performed [8,9]. However, the instrument presented in this work provides the means for a “direct calorimetry” technique. While the modern instruments for CO2 and O2 measurement used in indirect calorimetry are actually very well developed, direct calorimetry applied to living things has seen very little development because the organism being studied must be located within an enclosed and controlled chamber [10,11]. The sensor presented here can take direct measurements and in any environment.
We should mention that a thermopile acting as a heat flow meter has been used to measure solar energy incident on a free surface; however, pyrometers are more efficient for these applications [12]. Thermopile-based systems have also been created for use in other applications designed to harness electrical energy from body heat emissions [13,14].
In previous publications [4,5,6,7], we described different aspects of the first prototype’s performance. With the present study we aim to provide a comprehensive account of the second prototype’s performance, emphasising the most important characteristics and publishing some of the experimental data obtained from measurements on the human body. We begin with a detailed description of the experimental setup, then we model the device and finally present the experimental verification of its validity. Modelling consists in treating the instrument as a linear time invariant system with two inputs, two outputs and four transfer functions. This model can be used to simulate the performance, define the sensor’s operating range and establish a method to measure the power dissipated by the human body.

2. Experimental Section

2.1. Sensor and Measurement and Control Instrumentation

The calorimetry sensor was primarily based around a 13.2 × 13.2 × 2.2 mm3 HOT20-65-F2A Laird thermopile (Part c in Figure 1) placed between a thermostat and a 20 × 20 × 1 mm3 aluminium plate (Part a in Figure 1), which was attached to the area where the heat flow is measured. The thermopile produced a calorimetric signal thanks to the Seebeck effect.
The thermostat comprised a small (10 × 12 × 3 mm3) aluminium block (Part d in Figure 1) containing a heating resistor (R = 16.74 Ω) and an RTD sensor (1PT100KN1515, Omega, Manchester, UK). The thermostat also included a cooling system (Parts e, f and g in Figure 1) based around a thermopile (HOT20-65-F2A, Laird, London, UK, 13.2 × 13.2 × 2.2 mm3), which absorbed heat from the thermostat through the Peltier effect, and an aluminium heat sink (with its corresponding fan) attached to the cooling thermopile’s hot surface.
The sides of the entire device were thermally insulated by covering them with expanded polystyrene and reflective aluminium foil to minimise conductive and radiative heat transfer (Part b in Figure 1). The device was also equipped with a second RTD, positioned externally to measure the room temperature associated with each experimental measurement. Figure 1 shows a diagram of the sensor and its constituents: a (aluminium plate), b (thermal insulation), c (measurement thermopile), d (thermostat), e (cooling thermopile), f (aluminium heat sink) and g (fan). The device included a resistor (R = 27 Ω) contained within a 10 × 10 mm2 copper plate which was placed on the calibration base (Figure 2).
A programmable power supply (Agilent E3631A, 80 W Triple Output Power Supply, 6 V, 5 A and ±25 V, 1 A, Keysight, Santa Clara, CA, USA) was used to power the Joule calibration resistor, the heating resistor located inside the thermostat and the Peltier cooling module. The fan was powered by an auxiliary power supply. The calorimetric signal, thermostat temperature and ambient temperature were all recorded using a data acquisition system (Agilent 34970A Data Acquisition/Switch Unit and an Agilent 34901A 20-Channel Multiplexer). We used an Agilent 82357B USB/GPIB Interface to provide a direct connection between laptop USB ports and GPIB instruments. The experimental setup was controlled with a program in MatLab®; the program used a PID controller to ensure the thermostat maintained a constant temperature and stored all the variables (calorimetric signal, temperatures and the power dissipated by each resistor) for subsequent data analysis. Baseline and data collection sampling periods were both set to one second. Figure 3 shows the instruments forming the measurement and control system.
The experimental setup also included a third power supply to run another resistor placed on the laboratory bench. This additional resistor was located underneath a copper plate with a larger surface area (60 × 35 mm2) than the sensor (20 × 20 mm2) and was used to simulate a measurement made on the surface of the human body (see Figure 3).

2.2. Sensor Operating Diagram

To explain how the sensor works, Figure 4 is a block diagram where each block represents a different part of the device and features its corresponding input and output variables. The diagram includes the temperature control loop for the thermostat, where Tref is the desired thermostat temperature and Tpid is the thermostat’s actual temperature. The PID controller was used to calculate the power Wpid dissipated in the thermostat according to the following equation:
W p i d = K p e ( t ) + K i 0 t e ( t ) d t + K d d e ( t ) d t
where the error e(t) = (Tref − Tpid). Having determined the power to be dissipated, a power limiter was introduced so that 0 < Wpid < (Wpid)max. In the present case, (Wpid)max = 734 mW. Parameters for the PID controller were determined experimentally using the Ziegler-Nichols tuning rules [15,16] and were as follows: Kp = 0.8 W·K−1; Ki = 0.14 W·K−1·s−1; Kd = 2.4 W·K−1·s. Hence we managed to reach a steady-state temperature (with oscillations of ±5 mK) in 200 s for a programmed 5 K temperature step.
The block representing the calorimetry sensor comprises three sub-blocks: The first sub-block is the thermal insulation block which served to attenuate the ambient room temperature Troom. The second sub-block is the Peltier element whose function was to decrease the thermostat temperature, if necessary. Although the thermopile can be calibrated separately [17], we decided to perform calibration on the actual configuration. In this case, the functional relationship between the applied voltage (VPeltier) and the temperature difference (ΔTPeltier), with the fan turned on, is as follows:
Δ T P e l t i e r = T p i d T r o o m * = 1.7 9.3 V P e l t i e r + 1.5 V P e l t i e r 2
The third sub-block is the calorimetric model. Its inputs were the powers passing across the sensor (W) and dissipated in the thermostat (Wpid), while the outputs were the calorimetric signals taken from the measurement thermopile (ycal) and any fluctuations recorded in the thermostat’s actual temperature (ΔTpid). Therefore, it is a multiple-input, multiple-output (MIMO) system.
Finally, and in accordance with the schematic in Figure 4, the thermostat’s actual temperature is given by:
T p i d = T r o o m * + Δ T p e l t i e r + Δ T p i d

3. Sensor Modelling and Calibration

3.1. Modelling

As shown schematically in Figure 4, the model chosen for sensor identification purposes consists in treating the instrument as a linear time invariant system with two inputs, two outputs and four transfer functions. The relationship between the inputs and outputs in Laplace domain is as follows:
( Δ Y c a l ( s ) Δ T p i d ( s ) ) = ( T F 1 T F 2 T F 3 T F 4 ) ( Δ W ( s ) Δ W p i d ( s ) )
where W(s), Wpid(s), ΔY(s) and ΔTpid(s) are the Laplace transforms of the inputs and outputs, and each transfer function (TFi) can be calculated from:
T F i ( s ) = K i ( 1 + s τ i ) ( 1 + s τ 1 ) ( 1 + s τ 2 )
where Ki is the sensitivity or steady-state response to a unit step, and τi = −1/si and τi* = −1/si* (si represents the poles and si* the zeros for each TFi). This is equivalent to a total of four parameters for each TFi; if we force all the TFi to take the same pole, this will give a total of 10 parameters. To further reduce the number of parameters, we selected a model based on what is known as the localised-constant model, a model widely used in calorimetry [17,18,19,20]. Taking into account that the calorimetric signal depends on the temperature recorded on each surface of the measurement thermopile, we can consider two separate bodies with these temperatures. Hence, the first body represents the domain encompassing the dissipation site and one of the thermopile’s surfaces, while the second body represents the domain formed by the thermostat and the other side of the thermopile. Each body was assumed to have an infinite thermal conductivity and heat capacity Ci. In addition, each body was connected to the other and to the exterior using thermal couplings with a thermal conductivity of Pik (see Figure 5).
The power Wi in each domain is equal to the power required to change its temperature plus the power transmitted by conduction to the neighbouring domains, where each domain’s power can be calculated from:
W 1 = C 1 d T 1 d t + P 12 ( T 1 T 2 ) + P 1 ( T 1 T 0 ) W 2 = C 2 d T 2 d t + P 12 ( T 2 T 1 ) + P 2 ( T 2 T 0 )
If we assume that T0 and T′0 are constants and correct the baselines for all the temporal variables, then we get:
Δ W 1 = C 1 d Δ T 1 d t + P 12 ( Δ T 1 Δ T 2 ) + P 1 Δ T 1 Δ W 2 = C 2 d Δ T 2 d t + P 12 ( Δ T 2 Δ T 1 ) + P 2 Δ T 2
Performing a Laplace transform, with null initial conditions having corrected the baselines, gives:
( C 1 s + P 1 + P 12 P 12 P 12 C 2 s + P 2 + P 12 ) ( Δ W 1 Δ W 2 ) = ( Δ T 1 Δ T 2 )
Considering that the calorimetric output is proportional to the difference in temperatures:
Δ Y c a l = k ( T 1 T 2 ) = k ( Δ T 1 Δ T 2 )
If we proceed with W = W1, Wpid = W2 and ΔTpid = ΔT2, then we obtain the relationship between the inputs and outputs given in Equation (4), where each TFi is given by:
T F 1 ( s ) = k ( C 2 s + P 2 ) Δ T F 2 ( s ) = k ( C 1 s + P 1 ) Δ T F 3 ( s ) = P 12 Δ T F 4 ( s ) = C 1 s + P 1 + P 12 Δ w h e r e Δ = ( C 1 s + P 1 + P 12 ) ( C 2 s + P 2 + P 12 ) P 12 2
The four TFi have the same poles, but different sensitivities and zeros. With this approach we managed to reduce the number of model parameters to just six: C1, C2, P1, P12, P2 and k. The sensitivities, poles and zeros for each TFi were calculated from these parameters using Equation (10).

3.2. Calibration

The model parameters (Table 1) were determined by minimising a certain error criteria between the experimental ( Δ y c a l , Δ T p i d ) and theoretical curves ( Δ y ^ c a l , Δ T ^ p i d ) calculated with equations from the model (Equation (7)). We achieved this end using the Nelder–Mead simplex search algorithm [21] and MatLab software [22]. The error criterion selected (σ) was the mean squared error given by the equations:
σ y = 1 N i = 1 N ( Δ y c a l [ i ] Δ y ^ c a l [ i ] ) 2 σ T = 1 N i = 1 N ( Δ T p i d [ i ] Δ T ^ p i d [ i ] ) 2 σ = σ y + γ σ T
To carry out model identification, the sensor was placed on the calibration base and the thermostat set to 24 °C. Once the system reached a steady state, the thermostat was reprogrammed to 28 °C and two 117 mW pulses were simultaneously dissipated across the calibration resistor, thereby producing some significant variations in the input (W and Wpid) and output (ycal and Tpid) curves. These curves and the room temperature (Troom) during measurement are shown in Figure 6. Figure 7 shows the excellent fit between the experimental output curves and those calculated with the model. Table 1 presents the result of the calibration for the minisensor described in this work (which we called S4, given the 4 cm2 detection area), as well as the parameters for the first prototype, S36 (with a 36 cm2 detection area) [5]. The static response sensitivity of the new sensor (S4) was checked and found to be higher (K1 is about three times higher). Figure 8 compares the dynamic responses of the different TFi for the S4 and S36 sensors; the magnitudes are presented in dB corresponding to each frequency, and assuming unit sensitivity (Ki = 1) in all cases. TF1, TF2 and TF4 decreased by 20 dB per decade; they could then be approximated with TF and just one time constant, with results of 17.5, 7.3 and 48.0 s, respectively. The S36 sensor, however, had greater time constants: 33, 15 and 85 s. As expected, a smaller sensor involves lower heat capacity elements (Ci) that decrease the time constants and, in our case, the response was twice as fast because the time constants of the prototype decreased by half.

4. Results

4.1. Method for Measuring and Calculating the Mean Power Dissipated from the Surface of the Human Body

The instrument was operated as follows. Before taking a reading from a human body, the sensor was situated on the base until it reached the temperature set on the thermostat. The sensor was then placed on the area of the body to be measured until the thermostat again recorded a steady-state temperature (Figure 9), at which point the sensor was returned to its base. A calibration measurement by the Joule effect can be taken before and/or after taking a measurement on the human body, i.e., when the sensor is on its base.
According to the equations from the model (Equation (10)), the sensitivities (Ki = TFi(0)) are independent of C1 and the calibration base. This is confirmed experimentally by changing the size of the calibration resistor and observing that, while the time constants vary when this resistor is changed, the Ki sensitivities remain constant. Thus, Equation (4) for a steady-state condition is transformed into an equation for the sensor that is independent of the size of surface measured:
( Δ Y c a l Δ T p i d ) = ( K 1 K 2 K 3 K 4 ) ( Δ W Δ W p i d )
The calculation of the mean power Wmean flowing across the sensor during the sensor application period td consists of determining the areas under the calorimetric curves and the curves for the power developed in the thermostat after correcting the baselines:
W m e a n = 1 t d W ( t ) d t = 1 t d K 1 ( y c a l ( t ) d t K 2 W p i d ( t ) d t )
For the dissipation in the calibration base, as shown in Figure 10a,b, the calculated Wmean = 303.4 mW, while the experimental Wmean = 305.1 mW; this represents a difference of −0.6% (measurement made with Tpid = 28.0 °C and Troom = 21.8 °C).
To check the method’s validity we took a measurement using a large resistor placed on the laboratory workbench (Figure 3); the resistor was placed in the centre, under a thin, 20 cm2 sheet of copper. This resistor was used to simulate a measurement on the surface of the human body, where the dissipation occurs across the whole surface yet the measurement is made on an area of just 4 cm2. Figure 10c,d shows three consecutive measurements taken using the following method: the sensor was placed on the calibration base and once it reached a steady-state temperature it was placed on the large resistor; the sensor was then returned to the calibration base and the process repeated three times. The results were: 191.6, 184.0 and 187.2 mW, in other words 187.8 ± 3.8 mW. In this case the large resistor dissipated 404 mW, but it must be taken into account that the surface of this resistor was 20 cm2, while the detection surface was 4 cm2 (measurement made with Tpid = 28.0 °C and Troom = 21.3 °C). The measured power was not proportional to the area because the dissipation was not uniformly distributed over copper surface (it was higher in the central zone).
We performed several measurements on the human body, mainly on the left wrist. The results of this type of measurement depend on the individual’s condition; however, in this work, we did not study the measurements while correlating them against the test subject’s metabolic or emotional status, we instead focused on the reproducibility of the results recorded with the sensor. As an example, Figure 11 shows seven consecutive, short measurements (50–90 s) made on the left wrist of a healthy, 58-year-old male subject; the powers recorded were 147, 162, 173, 177, 167, 162 and 150 mW, giving a mean value of 162 ± 15 mW, measured at Troom = 21.0 °C and Tpid = 28.0 °C. Recordings made on the same subject but a different day gave results of 219 mW on the left wrist, 245 mW on the forehead, and 282 mW on the left pectoral region, where Troom = 21.4 °C and Tpid = 28.0 °C. Figure 12 displays three consecutive, long measurements (180 s) made on the left wrist of a healthy, 22-year-old male subject; the results were 199, 214 and 219 mW (209 ± 10 mW), measured at Troom = 21.4 °C and Tpid = 28.0 °C. Note that this last measurement took a total time of 30 minutes including the initial thermostat temperature stabilisation period. We conclude that the device and measurement method are valid, and the results obtained lie within an uncertainty of ±15 mW.
We also included in this work measurements of the left wrist of a healthy 58-year-old male at different thermostat temperatures ranging from 24 °C to 36 °C and at an ambient temperature of 26 °C. Figure 13 shows three sets of measurements made on three different days. Each series consists of seven consecutive measurements, one for each thermostat temperature. The values obtained can be linearly adjusted, the place the slope cuts the horizontal axis depends on the physical and/or emotional state of the subject.

4.2. Sensor Resolution and Operating Domain

Operating domain refers to the thermostat’s viable temperature range (Tpid) and to the minimum and maximum measurable powers (W). The temperature range depends on the ambient room temperature (Troom). It can be calculated using Equations (3) and (12), arriving at the following equation:
T p i d = T r o o m + Δ T P e l t i e r + K 3 W + K 4 W p i d
The temperature decrease caused by the Peltier element is given by Equation (2). The sensitivities are K3 = 10.1 K/W and K4 = 12.0 K/W (Table 1). Initially, we shall make the following assumptions: there is no dissipation (W = 0), (ΔTPeltier)min = −12.7 K, (ΔTPeltier)max = 1.7 K, and (Wpid)max = 734 mW. Then, we shall calculate the minimum and maximum possible values for the thermostat temperature: (Tpid)min = Troom − 12.7; (Tpid)max = Troom + 10.5. However, the permissible temperature range depends on the maximum measurable power as described below.
With respect to the sensor’s operating range, we must also consider the minimum and maximum measurable powers. The minimum power is related to the noise from the signals. At steady state, the oscillations in the calorimetric signal and in the power dissipated by the thermostat are ±0.2 mV and ±10 mW, respectively; hence, the oscillations in the calculated power will be ±4.3 mW. Therefore, we can conclude that the sensor has a resolution of 5 mW. However, the measurements have a greater uncertainty because the sensor is applied to the skin manually; this translates into a varying degree of contact in function of the pressure exerted by the person taking the reading. The measurement shown in Figure 11 reveals the experimental uncertainty was ±15 mW.
The maximum measurable power is related to the possible saturation of temperature control. The sensor becomes saturated when it can no longer maintain the preset temperature, either because there is insufficient power in the thermostat (upper saturation) or the measured power is too high (lower saturation). Temperature control ensures the power dissipated in the thermostat Wpid decreases, in order to maintain a constant thermostat temperature, when the input power W increases. Equation (14) can be used to determine a feasible value for ΔTPeltier (between −12.7 and 1.7 K) and Tpid (between 26 and 30 °C) over a room temperature range of 20–30 °C; under these conditions powers of up to 855 mW can be measured without saturating the sensor. Table 2 shows the voltages that must be applied to the cooling thermopile to produce the ΔTPeltier necessary to reach the preset thermostat temperature at a given ambient temperature; in the examples presented in Table 2, the maximum thermostat power (Wpid)max was taken as 720 mW. In practice, the measurement program calculates and applies the voltage VPeltier based on the room temperature and the temperatures set for the thermostat.

5. Conclusions

We have developed a calorimetric sensor that can perform a localised measurement of surface heat dissipated from the human body. The prototype was designed to measure a square surface area covering 4 cm2. The model derived from the conductive heat transfer equations is valid and reliably reproduces and simulates the sensor’s operation. It therefore establishes the operating domain and a method for calculating the mean power passing over the sensor. Heat power is hard to measure, as the thermal process being studied must be insulated to minimise the effect of all external influences; in calorimetry, any measurement system operating with an uncertainty of less than 2% can be considered a success. In the present study, the working conditions were very unfavourable as the sensor was used to take measurements in an open environment, hence thermal results with a 10% uncertainty are acceptable. The sensor provides true values for surface heat dissipated from localised areas of the human body; it can measure powers of 5–855 mW for thermostat temperatures of 26–30 °C and over a wide range of room temperatures (20–30 °C). It is easy to operate and has a very quick response. It takes a total of 30 min to perform the measurement in triplicate. In conclusion, the sensor is ready for use in clinical trials.

Author Contributions

All authors contributed equally to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wendlandt, W.W. Thermal Methods of Analysis; John Wiley and Sons: New York, NY, USA, 1974. [Google Scholar]
  2. Brown, M.E. (Ed.) Handbook of Thermal Analysis and Calorimetry: Principles and Practice; Elsevier Science: Amsterdam, The Netherlands, 1998; Volume 1.
  3. Hansen, L.D. Toward a standard nomenclature for calorimetry. Thermochim. Acta 2001, 371, 19–22. [Google Scholar] [CrossRef]
  4. Socorro, F.; de Rivera, M.R. Development of a calorimetric sensor for medical application. Part I. Operating model. J. Therm. Anal. Calorim. 2010, 99, 799–802. [Google Scholar] [CrossRef]
  5. Jesús, C.; Socorro, F.; de Rivera, M.R. Development of a calorimetric sensor for medical application. Part II. Identification and Simulation. J. Therm. Anal. Calorim. 2013, 113, 1003–1007. [Google Scholar] [CrossRef]
  6. Jesús, C.; Socorro, F.; de Rivera, M.R. Development of a calorimetric sensor for medical application. Part III. Operating methods and applications. J. Therm. Anal. Calorim. 2013, 113, 1009–1013. [Google Scholar] [CrossRef]
  7. Jesús, C.; Socorro, F.; de Rivera, H.J.R.; de Rivera, M.R. Development of a calorimetric sensor for medical application. Part IV. Deconvolution of the calorimetric signal. J. Therm. Anal. Calorim. 2014, 116, 151–155. [Google Scholar] [CrossRef]
  8. Brooks, A.G.; Withers, R.T.; Gore, C.J.; Vogler, A.J.; Plummer, J.; Cormack, J. Measurement and prediction of METs during household activities in 35- to 45-year-old females. Eur. J. Appl. Physiol. 2004, 91, 638–648. [Google Scholar] [CrossRef]
  9. Kozey, S.; Lyden, K.; Staudenmayer, J.; Freedson, P. Errors in MET estimates of physical activities using 3.5 mL × kg(−1) × min(−1) as the baseline oxygen consumption. J. Phys. Act. Health 2010, 7, 508–516. [Google Scholar] [CrossRef]
  10. McLean, J.L.; Tobin, G. Animal and Human Calorimetry; Cambridge University Press: Cambridge, UK, 1987. [Google Scholar]
  11. Lamprecht, I. Calorimetry and thermodynamics of living systems. Thermochim. Acta 2003, 405, 1–13. [Google Scholar] [CrossRef]
  12. Hukseflux. Thermal Sensors. Available online: http://www.hukseflux.com/ (accessed on 4 November 2016).
  13. Leonov, V. Thermoelectric energy harvester on heated human machine. J. Micromech. Microeng. 2011, 21, 125013. [Google Scholar] [CrossRef]
  14. Carmo, J.P.; Antunes, J.; Silva, M.F.; Ribeiro, J.F.; Goncalves, L.M.; Correia, J.H. Characterization of thermoelectric generators by measuring the load-dependence behavior. Measurement 2011, 44, 2194–2199. [Google Scholar] [CrossRef]
  15. Ogata, K. Modern Control Engineering; Prentice Hall: Upper Saddle River, NJ, USA, 2010. [Google Scholar]
  16. Goodwin, G.C.; Graebe, S.F.; Salgado, M.E. Control System Design; Prentice Hall: Upper Saddle River, NJ, USA, 2001. [Google Scholar]
  17. Martin, J. Protocols for the high temperature measurement of the Seebeck coefficient in thermoelectric materials. Meas. Sci. Technol. 2013, 34, 085601. [Google Scholar] [CrossRef]
  18. Isalgue, A.; Ortin, J.; Torra, V.; Viñals, J. Heat flux calorimeters: Dynamical model localized time constants. An. Fis. 1980, 76, 192–196. [Google Scholar]
  19. Socorro, F.; de Rivera, M.R.; Jesús, C. A thermal model of a flow calorimeter. J. Therm. Anal. Calorim. 2001, 64, 357–366. [Google Scholar] [CrossRef]
  20. Kirchner, R.; de Rivera, M.R.; Seidel, J.M.; Torra, V. Identification of micro-scale calorimetric devices. Part VI. An approach by RC-representative model to improvements in TAM microcalorimeters. J. Therm. Anal. Calorim. 2005, 82, 179–184. [Google Scholar] [CrossRef]
  21. Nelder, J.A.; Mead, C. A simplex method for function minimization. Comput. J. 1965, 7, 308–313. [Google Scholar] [CrossRef]
  22. The MathWorks, Inc. Optimization ToolboxTM User’s Guide; 5th printing; Revised for Version 3.0 (Release 14); The MathWorks, Inc.: Natick, MA, USA, 2004. [Google Scholar]
Figure 1. Cross-sectional diagram and exploded view of the calorimetry sensor: a (aluminium plate), b (thermal insulation), c (measurement thermopile), d (thermostat), e (cooling thermopile), f (aluminium heat sink) and g (fan).
Figure 1. Cross-sectional diagram and exploded view of the calorimetry sensor: a (aluminium plate), b (thermal insulation), c (measurement thermopile), d (thermostat), e (cooling thermopile), f (aluminium heat sink) and g (fan).
Sensors 16 01864 g001
Figure 2. Base with the calibration resistor (a) and the sensor placed on the calibration base (b).
Figure 2. Base with the calibration resistor (a) and the sensor placed on the calibration base (b).
Sensors 16 01864 g002
Figure 3. Measurement and control instrumentation.
Figure 3. Measurement and control instrumentation.
Sensors 16 01864 g003
Figure 4. Calorimetry sensor operating diagram.
Figure 4. Calorimetry sensor operating diagram.
Sensors 16 01864 g004
Figure 5. Calorimetry sensor model: C1 represents the power dissipation element and one surface of the measurement thermopile, C2 represents the thermostat and the thermopile’s other surface. T0 and T′0 represent temperatures outside the sensor. W1 and W2 are the powers dissipated in each domain.
Figure 5. Calorimetry sensor model: C1 represents the power dissipation element and one surface of the measurement thermopile, C2 represents the thermostat and the thermopile’s other surface. T0 and T′0 represent temperatures outside the sensor. W1 and W2 are the powers dissipated in each domain.
Sensors 16 01864 g005
Figure 6. Input and output curves used in the calibration: (a) Calorimetric signal; (b) Power dissipated in the thermostat resistor; (c) Power dissipated into the calibration base by the Joule effect; (d) Thermostat temperature; and (e) Room temperature.
Figure 6. Input and output curves used in the calibration: (a) Calorimetric signal; (b) Power dissipated in the thermostat resistor; (c) Power dissipated into the calibration base by the Joule effect; (d) Thermostat temperature; and (e) Room temperature.
Sensors 16 01864 g006
Figure 7. Calibration. Fit between the experimental curves (lines) and those calculated using the model (points): (a) Calorimetric signal; and (b) thermostat temperature. Curves feature a corrected baseline.
Figure 7. Calibration. Fit between the experimental curves (lines) and those calculated using the model (points): (a) Calorimetric signal; and (b) thermostat temperature. Curves feature a corrected baseline.
Sensors 16 01864 g007
Figure 8. Log-magnitude curves for the transfer function given by Equations (4) and (5). S4, 4 cm2 minisensor described in this study. S36, an earlier 36 cm2 sensor (Reference [5]).
Figure 8. Log-magnitude curves for the transfer function given by Equations (4) and (5). S4, 4 cm2 minisensor described in this study. S36, an earlier 36 cm2 sensor (Reference [5]).
Sensors 16 01864 g008
Figure 9. Placing the S36 (a) and S4 (b) sensors on the human body.
Figure 9. Placing the S36 (a) and S4 (b) sensors on the human body.
Sensors 16 01864 g009
Figure 10. Measurement performed on the calibration base: (a) calorimetric signal; and (b) power dissipated in the thermostat. Measurement performed on the large resistor placed on the laboratory workbench: (c) calorimetric signal; and (d) power dissipated in the thermostat. The baseline was corrected for the four curves displayed.
Figure 10. Measurement performed on the calibration base: (a) calorimetric signal; and (b) power dissipated in the thermostat. Measurement performed on the large resistor placed on the laboratory workbench: (c) calorimetric signal; and (d) power dissipated in the thermostat. The baseline was corrected for the four curves displayed.
Sensors 16 01864 g010
Figure 11. Consecutive measurements made on the left wrist of a healthy, 58-year-old male subject. Calculated mean powers: 147, 162, 173, 177, 167, 162 and 150 mW (162 ± 15 mW), at Troom = 21 °C, Tpid = 28 °C. The baselines of the curves for the calorimetric signal and the thermostat’s power were corrected. Thermostat temperature is also shown (Tpid). VPeltier = 0.2 V.
Figure 11. Consecutive measurements made on the left wrist of a healthy, 58-year-old male subject. Calculated mean powers: 147, 162, 173, 177, 167, 162 and 150 mW (162 ± 15 mW), at Troom = 21 °C, Tpid = 28 °C. The baselines of the curves for the calorimetric signal and the thermostat’s power were corrected. Thermostat temperature is also shown (Tpid). VPeltier = 0.2 V.
Sensors 16 01864 g011
Figure 12. Three consecutive measurements made on the left wrist of a healthy, 22-year-old male subject. Calculated mean powers: 199, 214 and 219 mW (209 ± 10 mW) at Troom = 21.4 °C, Tpid = 28.0 °C. The baselines of the curves for the calorimetric signal and the thermostat’s power were corrected. Thermostat temperature is also shown (Tpid). VPeltier = 0.2 V.
Figure 12. Three consecutive measurements made on the left wrist of a healthy, 22-year-old male subject. Calculated mean powers: 199, 214 and 219 mW (209 ± 10 mW) at Troom = 21.4 °C, Tpid = 28.0 °C. The baselines of the curves for the calorimetric signal and the thermostat’s power were corrected. Thermostat temperature is also shown (Tpid). VPeltier = 0.2 V.
Sensors 16 01864 g012
Figure 13. Three measurements made on the left wrist of a healthy, 58-year-old male subject on three different days (squares, circles and diamonds) for different thermostat temperatures Tpid (Troom = 26 °C).
Figure 13. Three measurements made on the left wrist of a healthy, 58-year-old male subject on three different days (squares, circles and diamonds) for different thermostat temperatures Tpid (Troom = 26 °C).
Sensors 16 01864 g013
Table 1. Calibration results. RC model (Equation (7)) and TFi (Equation (5)) parameters. S4, the minisensor described in this study; S36, an earlier 36 cm2 model (Reference [5]).
Table 1. Calibration results. RC model (Equation (7)) and TFi (Equation (5)) parameters. S4, the minisensor described in this study; S36, an earlier 36 cm2 model (Reference [5]).
RC ModelTF Model
S4 Sensor S4 SensorS36 Sensor
C12.9146 J·K−1K1148.7 mV·W−146.0 mV·W−1
C23.9735 J·K−1K2−45.4 mV·W−1−9.3 mV·W−1
P10.0203 W·K−1K310.1 K·W−11.76 K·W−1
P20.0663 W·K−1K412.0 K·W−11.70 K·W−1
P120.1116 W·K−1τ182 s147 s
K24.6678 mV·K−1τ213 s27 s
Errors (Equation (11))τ1*60 s91 s
σy29.5 µVτ2*144 s243 s
σT1.5 mKτ3*0 s0 s
N1100 pointsτ4*22 s49 s
Table 2. Sensor operating domain for a thermostat temperature (Tpid) range of 26–30 °C and at an ambient temperature (Troom) of 20–30 °C. The table shows the voltage VPeltier across the cooling thermopile required to produce the temperature decrease ΔTPeltier, and the maximum measurable power (Wmax) for each example.
Table 2. Sensor operating domain for a thermostat temperature (Tpid) range of 26–30 °C and at an ambient temperature (Troom) of 20–30 °C. The table shows the voltage VPeltier across the cooling thermopile required to produce the temperature decrease ΔTPeltier, and the maximum measurable power (Wmax) for each example.
TpidTroomVPeltier (V)ΔTPeltier (°C)Wmax (mW)
26 °C20 °C0.51−2.64855
25 °C1.26−7.64855
30 °C2.88−12.64855
28 °C20 °C0.26−0.64855
25 °C0.93−5.64855
30 °C1.92−10.64855
30 °C20 °C0.04+1.36855
25 °C0.64−3.64855
30 °C1.45−8.64855

Share and Cite

MDPI and ACS Style

Socorro, F.; Rodríguez de Rivera, P.J.; Rodríguez de Rivera, M. Calorimetry Minisensor for the Localised Measurement of Surface Heat Dissipated from the Human Body. Sensors 2016, 16, 1864. https://doi.org/10.3390/s16111864

AMA Style

Socorro F, Rodríguez de Rivera PJ, Rodríguez de Rivera M. Calorimetry Minisensor for the Localised Measurement of Surface Heat Dissipated from the Human Body. Sensors. 2016; 16(11):1864. https://doi.org/10.3390/s16111864

Chicago/Turabian Style

Socorro, Fabiola, Pedro Jesús Rodríguez de Rivera, and Manuel Rodríguez de Rivera. 2016. "Calorimetry Minisensor for the Localised Measurement of Surface Heat Dissipated from the Human Body" Sensors 16, no. 11: 1864. https://doi.org/10.3390/s16111864

APA Style

Socorro, F., Rodríguez de Rivera, P. J., & Rodríguez de Rivera, M. (2016). Calorimetry Minisensor for the Localised Measurement of Surface Heat Dissipated from the Human Body. Sensors, 16(11), 1864. https://doi.org/10.3390/s16111864

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop