Issue 
Int. J. Metrol. Qual. Eng.
Volume 8, 2017



Article Number  18  
Number of page(s)  13  
DOI  https://doi.org/10.1051/ijmqe/2017007  
Published online  05 June 2017 
Research Article
Characterization of a spherical heat source for measuring thermal conductivity and water content of ethanol and water mixtures
^{1}
National Institute of Metrology, Quality and Technology (Inmetro),
Av. Nossa Senhora das Graças, 50, Xerém,
Duque de Caxias,
RJ,
CEP: 25250020, Brazil
^{2}
Pontifical Catholic University of Rio de Janeiro (PUCRio),
Rua Marquês de São Vicente, 225, Gávea,
Rio de Janeiro,
RJ,
CEP: 22451900, Brazil
^{3}
Georges Bonnier Consultant,
Rue Charles Perrault, 14,
Montsoult
95560, France
^{⁎} Corresponding author: jdbrionizio@inmetro.gov.br
Received:
4
July
2016
Accepted:
2
May
2017
The study and the development of measuring methods of thermal conductivity are essential in several engineering applications, since as a consequence of the current justified demands on saving and rational use of thermal energy, the heat transfer with the maximum efficient as possible is of great relevance. The measurement of the water content is also a relevant parameter in several research areas and industrial sectors, since the quantity of water in the substances influences several biological, chemical and physical processes. The aim of this paper is to present an experimental and theoretical study, following the good metrological practices, of a method based on a spherical heat source in order to measure the thermal conductivity of liquids, focusing on water, ethanol and their mixtures, with later determination of the water content of the binary samples.
Key words: thermal conductivity / water content / spherical heat source / metrology
© J.D. Brionizio et al., published by EDP Sciences, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Due to the environmental needs, the rational use and saving of thermal energy have become increasingly relevant in the last decades, especially in the industrial sectors. Thus, the heat transfer with the best efficiency is as important as minimizing the heat loss through the use of appropriate materials. Hence, the knowledge of the thermal conductivity of materials is of great relevance in several engineering applications where heat transfer plays a fundamental role. The water content is also another important quantity for many industrial sectors and research areas, since the amount of water in the substances influences several biological, chemical and physical processes. Nevertheless, the number of equipment available on the market for the measurement of these quantities is not vast when compared to other quantities, mainly for liquid substances. For this reason, researches on the development of alternative measuring techniques have considerable significance.
The thermal conductivity is a transport property that furnishes an indication of the rate at which energy is transferred by the conduction process (net transfer of energy by random molecular motion). The measuring methods of the thermal conductivity are often classified as being steadystate and transient. In the case of fluids, its measurement is more difficult than that of solids due to the possibility of the onset of natural convection, which consists in energy transfer due to random molecular motion and the global motion of the fluid. This fluid motion is due to the buoyancy forces induced by density differences caused by temperature variations in the fluid.
The thermal conductivity depends on several parameters, and the water content is one of them. Consequently, the water content can be indirectly determined by means of the thermal conductivity value of the material, which is higher as the water content increases, since the thermal conductivity of water is higher than that of many nonmetallic solids, gases and liquids.
For a long while researchers have been working on new approaches for theoretical estimation of thermal conductivity of liquids, by linking it to other physical property, usually using temperature as a variable, or by deriving a model from the existing ones. Also, several experiments have been developed, using different techniques based on steadystate and transient methods, in order to measure it properly.
The line heat source probe (hotwire method) is usually recognized as the most accurate technique for measuring thermal conductivity of liquids. The concept and the initial experiments with heated wires started about 1780, and the first transient hotwire instrument was proposed by Stâlhane and Pyk in 1931 to measure the thermal conductivity of solids, powders and some liquids. Since 1780, the method has been studied by several researchers [1]. The standard test method for the determination of thermal conductivity of nonmetallic liquids of the American Society for Testing and Materials (ASTM) is based on this method [2]. An extensive uncertainty assessment elaborated for the hotwire method obtained the value of 5.8% for the thermal conductivity [3], although other literatures claim for uncertainty values better than 5%.
Models based on spherical heat sources have also been studied to measure thermal conductivity of different materials (solid, powder, paste, etc.). In this technique, the spherical device simultaneously generates heat and monitors the temperature response with time. Heat sources of spherical symmetry are free of lateral thermal effects and they yield to the steadystate regime at long times.
In 1933, Gibbs [4] described a device primarily designed to determine qualitatively the changes in blood flow through a tissue or organ by means of an electrically heated thermocouple. From 1951, Grayson [5–7] discussed the use of the apparatus for measuring the thermal conductivity, taking into consideration a steadystate regime in which the relation between heat production and heat loss is represented by equation (2), described by Carslaw [8] in 1921 for an electrically heated spherical source in an infinite mass of material. In 1968, Chato [9] was to first to use a thermistor as temperature sensor and heating element for measuring thermal conductivity and thermal diffusivity of biological materials. Several other researchers have also applied the thermistor method for measuring thermal properties of biomaterials [10–15]. Woodbury [16] used the thermistor method to measure thermal conductivity of building insulation with varying degree of wetness. Fujii et al. [17] employed the technique to measure thermal conductivities of liquid mixtures of water–ethanol, water–methanol and R113oil. Dougherty [18] used the thermistor method to perform thermal conductivity measurements in materials ranging from low viscosity fluids to insulation ones. Kravets [19] used the technique to measure thermal conductivity of milk and cream over the range of 25 °C to 125 °C. Holeschovsky et al. [20] used the thermistor method to measure the thermal conductivity of liquids and silica gel. Radhakrishnan [21] measured the thermal properties of ten different seafood in the range from 5 °C to 30 °C. Gelder [22] employed the thermistor technique to measure thermal properties of moist food materials at high temperatures. In 2003, Zhang et al. [23], based on a transient thermal model of a point heat pulse, developed a dualthermistor probe, in which two thermistor beads serve as point heater and temperature sensor, in order to measure thermal diffusivity and thermal conductivity. From 2005, Kubicar et al. [24–26] presented a sensor to measure thermal conductivity, in which a heat source and a thermometer were unified in a single component in the form of a small ball.
In principle, a model based on a spherical heat source is an absolute measuring method of thermal conductivity, which means that the sensor can furnish a result without being calibrated against a standard or a reference material. Nevertheless, parameters of the model need to be determined or obtained by means of calibration.
The research was focused on ethanol because it is a strategic biofuel in Brazil, and also because it is widely used as feedstock (pure or mixed with water) in several industrial sectors, such as pharmaceutical, food processing, hygiene and others. Since the early days, ethanol has helped fuel motorized transport in Brazil. In 1903, Brazil's First National Congress on Industrial Applications of Alcohol recommended the development of infrastructure to produce automotive ethanol. Several ethanol programs were implemented by the Brazilian government since then [27]. As a consequence, many ethanol production plants are spread all over the country, which makes Brazil the world's second largest producer of ethanol fuel. So, studies on ethanol are of great relevance for the industries, especially those related to measurement processes.
2 Theoretical model
In the measuring principle of the thermal properties by spherical geometry, the sensor is inserted in the medium of interest. Electrical power is supplied to the sensor and a constant heat flow is generated inside. Through the radius of the sphere, the heat propagates through the surrounding medium by a radial distance, resulting in temperature rise of the medium until its stabilization. The heat dissipated by the ball to maintain this temperature rise depends on the thermal properties of the medium.
The heat transfer model of a spherical heat source is generally based on the following hypothesis: the sensor and the medium are in thermal equilibrium before heating; the surrounding medium is homogeneous, isotropic and infinite; the sensor is homogeneous and isotropic; the sensor has spherical shape; the heat is generated uniformly within the sensor; there is no thermal contact resistance between the sensor and the environment; and the only heat transfer mode acting in the process is conduction.
The temperature rise of the spherical heat source and the measured medium in the steadystate regime can be described by the following initial and boundary conditions: where the subscripts s and m represent the sphere and the medium, respectively; T is the temperature; t is the time; r is the radial distance; r_{s} is the radius of the sphere; and k is the thermal conductivity.
The temperature distribution of a region bounded internally by a sphere, with initial temperature as zero and constant heat flux in the interface, is given by [28]: (1) where q″ is the heat transfer per unit area (heat flux) and erfc(x) is the complementary error function.
For a temperature measured at the surface of the sphere (r = r_{s}) at long times (t → ∞), an equation for determining the thermal conductivity of the medium that surrounds the sensor in the steadystate regime can be obtained, as previously shown by Carslaw [8]: (2) where is the heat transfer rate.
The temperature distribution within a sphere, with initial temperature as zero, considering its heat capacity and the thermal contact resistance between the heated sphere and the surrounding medium is given by [28]: (3) where ; ϕ = H_{c}/k_{m} = 1/k_{m}R_{c}; ρ is the density; c is the specific heat; m is the mass; and H_{c} is the thermal contact conductance, which is the inverse of the thermal contact resistance (R_{c}).
For the steadystate regime, the temperature distribution within a sphere, with initial temperature as zero, considering the thermal contact resistance between the heated sphere and the surrounding medium is given by [28]: (4) where is the heat transfer rate per unit volume.
3 Materials and experimental methods
Two commercial measuring systems were used in the experiments. Each instrument is composed by a remote spherical sensor and an electronic unit. The sensor consists of a NTC (Negative Temperature Coefficient) thermistor of 47 kΩ and an electrical resistance of 100 Ω fixed together by means of epoxy resin, resulting in a sphere with nominal diameter of 2 mm. The operating range of the thermistor is −40 °C to 150 °C, the temperature response is 0.1 °C to 5 °C, and the power generation range of the device is 1 mW to 30 mW. Figure 1 shows one measuring system (left) and the spherical sensors (right).
The rate of heat transfer from the heated spherical sensor to the medium results from the generating process of internal power, which comes from the conversion of electrical energy into thermal energy. This phenomenon occurs due to resistive heating, i.e., an electric current (I) passes through a resistance (R) inside the sensor (control volume), electric power is then dissipated at a rate equal to I^{2}R which corresponds to the rate at which thermal energy is released.
The thermistors were calibrated in the range from 15 °C to 30 °C, in steps of 2.5 °C, by comparison against a platinum resistance thermometer of 100 Ω (at 0 °C) calibrated by the Thermometry Laboratory of Inmetro. The electrical power values adjusted in the devices were compared to those calculated by means of voltage and current measurements performed by precise digital multimeters (of 6½ and 7½ digits of resolution) calibrated by the Voltage and Electric Current Laboratory of Inmetro. The diameters of the sensors were measured by the Dimensional Metrology Laboratory of Inmetro by means of an optical microscope traceable to a standard ruler, which in turn was calibrated by a HeNe laser interferometer. The measurement results of the diameter of the sensors were 2.32 mm (sensor 01) and 2.08 mm (sensor 02), both having expanded measurement uncertainty of 0.03 mm, with coverage probability of 95.45%.
For the measurements, the following samples were used: distilled water, anhydrous ethanol (purity around 99.3%), ethyl alcohol (purity around 95%), solution of 75% of ethanol and 25% of water, solution of 50% of ethanol and 50% of water, and solution of 25% of ethanol and 75% of water. The water content of the binary samples was measured by the Organic Analysis Laboratory of Inmetro by means of a Karl Fischer coulometric titration. All the samples were stored in dark borosilicate glass bottles with the same characteristics and size. The samples had approximately volume of 1 L, although lower volumes could be used. In previous works [15,29,30], the volume of the sample was estimated from 5 to 12 times the radius of the sensor, resulting in values in the order of 1 cm^{3}.
When measuring the thermal conductivity of fluids, it is desirable that the heat transfer from the measuring system to the medium should occur only by conduction. So, care must be taken to prevent the onset of natural convection. In order to detect its onset, the experimental method developed by Schmidt and Milverton [31] was used. The researchers noted that in the heat transfer by conduction, the temperature difference imposed in the fluid increased linearly with the power supplied to the heating system of the apparatus, and this linearity was interrupted when convection started to act. Therefore, the spherical sensors were immerged in distilled water and ethyl alcohol at 20 °C, and thirty nominal values of electrical power (with increments of approximately 0.45 mW) were set in the devices in the range from 1 mW to 14 mW. For both sensors, it was clearly observed that the relation between the electrical power and the temperature increment was linear for the whole range. Nevertheless, for ethyl alcohol, this linearity was interrupted from 8 mW and 6 mW for the sensors numbered as 01 and 02, respectively, which is an indication of the convection onset.
The spherical sensor of the measuring system was placed approximately in the middle of the bottle, which in turn was conditioned at 20 °C in a thermostatic bath with stability lower than 0.02 °C. The sample within the bottle was completely covered by the fluid of the bath. The sensor was guided to the middle of the bottle by means of a very small hole prepared in the removable lid of the container. In order to avoid infiltration of air to the inner of the bottle, the tip and part of the sensor cable were wrapped with polyvinyl chloride (PVC) film. After removal of the spherical sensor from a sample, it was exhaustively rinsed with distilled water and dried with paper towel and hot air before submerging it on another sample. This was done to prevent crosscontamination.
Before starting the measurements, the samples were left for at least 12 h in the bath so as to guarantee their completely stabilization at 20 °C. The nominal electrical powers set in the electronic units were, respectively, 8 mW and 6 mW for the spherical sensors 01 and 02. Several measurement cycles were performed for each sensor, where each one was composed by temperature monitoring before, during and after the heating of the spherical sensor. Five measurements were acquired before the heating of the sphere, ninetyfive measurements were obtained during it, and fifty measurements were collected when the heating process was stopped. The interval between measurements was 1 s, so that each cycle took 2.5 min from the beginning to the end. An interval from 10 min to 20 min between cycles was given so as to prevent that the temperature rise of a cycle could affect the following one. This was confirmed by comparing the first measurement of a cycle with the last one of the previous cycle.
The thermal conductivity of the samples was determined in the steadystate regime. For this reason, it was very important to establish when this regime was reached. During the calibration of the thermometers, where several measurements were taken during 30 min after the stabilization of the liquid of the bath, an amplitude value of temperature oscillation was defined for each sensor. For both temperature sensors, the amplitude values were lower than 0.05 °C. The temperature measurements of the samples varied within the amplitude value 22 s after the beginning of the heating of the spherical sensor. Figure 2 shows the temperature monitoring during a measurement cycle.
Fig. 1 Measuring system (left) and spherical sensors (right). 
Fig. 2 Temperature monitoring during a measurement cycle. 
4 Method validation
The water content of the binary samples (ethanol and water mixtures) was compared to the reference measurements performed at the Organic Analysis Laboratory of Inmetro by means of a Karl Fischer coulometric titration. Table 1 presents for each sample the nominal ethanol content EC_{N}, the nominal water content WC_{N}, the reference water content WC_{R} and their expanded uncertainties, with a coverage factor of 95.45%.
The thermal conductivity of the samples was compared to literature values, since the institute does not have reference standards for measuring this quantity in liquids. A mean value was used as reference, once different values for the thermal conductivity of the samples at 20 °C were found in the literatures.
Using a hotwire equipment, Assael et al. [32] measured the thermal conductivity of the liquid mixtures of water–ethanol, water–methanol and water–propanol in the range from 26.85 °C to 71.85 °C. The thermal conductivity measurements were adjusted as a function of temperature, and the coefficients ϕ have been determined for pure components and for the mixtures, as shown in Table 2. The researchers compared their thermal conductivity values to five other works at 16.85 °C and 56.85 °C for the liquid mixtures of water–methanol and water–ethanol. The maximum difference between their measurements and those of the five works were 15%, although this difference could be reduced to 5% when threes works were considered (5)
Filippov [33] studied the thermal conductivity of liquid solutions as a function of their concentration. He analyzed twelve systems of unassociated components, ten systems with a single associated component, aqueous solutions of three alcohols and solutions of some salts and acids. From the thermal conductivities of the pure substances, an equation was proposed for estimating the thermal conductivity of binary mixtures: (6) where k_{1} and k_{2} are the thermal conductivities of the pure substances 1 and 2; w_{1} and w_{2} are the mass fractions of these substances; and λ is a coefficient that depends on the substances of the solution.
For solutions of unassociated substances, λ is equal to 0.72; for solutions of associated substances, λ ranges from 0.5 to 1; and for aqueous solutions it ranges from 0.3 to 0.7. According to Khan [34], based on 120 data, the value of λ for Filippov's equation for liquid mixtures of water–ethanol in the range from −70 °C to 60 °C is 0.571066.
Reid et al. have also proposed an equation for estimating the thermal conductivity of binary mixtures as a function of the concentrations and thermal conductivities of the pure substances [34]: (7)
According to Khan [34], the value of λ for the equation of Reid et al. for liquid mixtures of water–ethanol in the range from −70 °C to 60 °C is 0.051892.
From several experimental measurements performed by the hotwire method, Ramires et al. [35] proposed an equation for estimating the thermal conductivity of water in the range from 0.85 °C to 96.85 °C at atmospheric pressure: (8)
Miller and Yaws proposed an equation for estimating the thermal conductivity of ethanol in the range from −114 °C to 190 °C [36]: (9)
Touloukian et al. proposed an equation for estimating the thermal conductivity of pure liquids, where for ethanol in the range from −123 °C to 127 °C the coefficients λ_{1} and λ_{2} are 609.512 and −0.70924, respectively [34]: (10)
The Korea Thermophysical Properties Data Bank (KDB), which provides information and estimating methods for the thermophysical properties of several substances, presents the following equations for estimating the thermal conductivities of water (in the range from −0.15 °C to 349.85 °C) and ethanol (in the range from −113.15 °C to 189.85 °C), respectively [37]: (11) (12)
The International Association for the Properties of Water and Steam (IAPWS), a nonprofit association of national organizations concerned with the properties of water and steam, provides in its website an online calculator, prepared by Moscow Power Engineering Institute and Russian National Committee of IAPWS, for estimating the thermal conductivity of water as a function of temperature and density [38]. The calculator was based on the “Release on the IAPWS Formulation 2011 for the Thermal Conductivity of Ordinary Water Substance” [39]. The density of water was estimated by the equation presented by Tanaka et al. [40]. The output value of thermal conductivity is presented in Table 3.
For estimating the thermal conductivity values of the binary samples by means of the equations (6) and (7), the thermal conductivities of the pure substances were necessary. In these cases, the thermal conductivities of ethanol and water were obtained, respectively, from KDB [37] and IAPWS [38]. Besides the thermal conductivity values estimated by the aforementioned equations for water, ethanol and their mixtures, Table 3 also presents values which were interpolated or directly obtained in tables available in literatures [41–43].
Concentrations and expanded uncertainties of the samples.
Thermal conductivity values of the samples from literatures (in W/m °C) at 20 °C.
5 Preliminary results
For each sample, three measurement cycles were performed, and for each cycle a temperature step was determined as: (13) where, ΔT is the temperature step; is the mean temperature of the five readings before the heating of the sphere; and is the last five readings during the heating of sphere.
As explained before, during the calibration of the thermometers it was observed that the measurements varied within the range of 0.05 °C in a stabilized medium. This parameter was adopted as criterion to determine the steadystate condition, which was reached 22 s after the heating of the sphere. The temperature step could then be determined in approximately 30 s. However, a slightly longer measuring time was used.
Knowing the radius of the sphere, the heat transfer rate (electrical power) and the temperature step (measured as detailed in Sect. 3), the thermal conductivities of the samples could then be determined by means of equation (2). Nevertheless, the estimated results were highly unsatisfactory, since they were very diverging from the literature values shown in Table 3. The percentage differences from the estimated values to the literature ones increased from approximately 10–130% (sensor 01) and 30–160% (sensor 02) as the ethanol content increased.
Another option could be to use an effective radius instead of the geometrical radius of the sphere. An effective radius, which can be estimated by the calibration of the spherical sensor in certain substances (selected by the user according to his measurement needs), consists of a virtual radius that compensates the lack of knowledge of some parameters involved in the measurement process, so as to accurately reproduce the thermal conductivities of the substances taken as reference. Using the thermal conductivities of the samples as reference (Tab. 3), the heat transfer rate and the temperature steps, the effective radii for each sensor could be estimated. Nevertheless, for both sensors, the effective radii varied considerably from one sample to another. The effective radius of the sensors increased linearly as the ethanol content increased. The estimated effective radii for sensors 01 and 02 were respectively 1.28 mm and 1.34 mm for distilled water, and 2.67 mm and 2.73 mm for anhydrous ethanol. Due to these considerable variations, the results were not useful for the measurement of the ethanol and water mixtures, because it would be necessary to use an effective radius value according to the concentration of the sample to be measured, which is usually unknown information in practice. Except for the half and half sample, the use of a mean effective radius is not convenient either. The percentage differences from the estimated thermal conductivities to the literature values increased as the difference of the concentrations of the solutions increased, reaching approximately 35% in the pure substances (distilled water and anhydrous ethanol).
The determination of the thermal conductivity values of the samples by the previous methods did not succeed, because equation (2) should only be used in ideal models, i.e., free from error sources. Nevertheless, the real process presents errors. So, a new model needed to be developed.
6 Development of an experimental model
In an ideal model, the temperatures of the sphere and the medium are the same in the interface. However, as shown in [25], the temperature measured in the sphere is shifted due to the thermal gradient within the sphere and the thermal contact resistance between the heated sphere and the surrounding medium. It can be stated as: (14) where δ_{TG} is the temperature shift due to the thermal gradient within the sphere; δ_{TR} is the temperature shift due to the thermal contact resistance between the sphere and the medium; and δ_{T} is the temperature shift due to δ_{TG} and δ_{TR}.
Rearranging equation (14) as δ_{T} = T_{s} − T_{m}, and using equations (2) and (4), δ_{T} can be deduced as: (15)
Rearranging equation (14) as T_{m} = T_{s} − δ_{T}, equation (2) can be rewritten as: (16)
Two parameters were experimentally measured for the determination of the thermal conductivities of the samples: the heat transfer rate and the temperature step ΔT, which means that k should be expressed as a function of them. Equation (16) is valid when the initial temperature is zero; for the cases where it is not zero, ΔT is used instead of T_{s}. So, equation (16) can be rearranged as: (17)
Equation (17) shows that the thermal conductivity of the medium can be determined by means of the product of the experimental ratio and a term, named as ξ, that takes into consideration the possible temperature errors. The term ξ was calculated for each sample through the calibration of the spherical sensors by means of the following equation: (18) where k_{L} is the mean thermal conductivity of the sample obtained from literature.
Table 4 presents for each sample the experimental ratio , which is a mean value of the three measurement cycles, k_{L} and the terms ξ.
An equation relating the experimental ratios and the terms ξ could be adjusted for each sensor. This equation is the calibration curve of the instrument. It means that, knowing and measuring ΔT, the equation compensates the errors of the sensor and gives the thermal conductivity value for water, ethanol and any unknown mixture between these substances. Nevertheless, as the thermal conductivities of the binary samples do not vary linearly with the water content, a polynomial function is the best option. Thus, a third degree equation (with coefficient of determination equals 1) was then chosen for each sensor: (19) (20)
Table 5 presents the thermal conductivities of the samples, determined by means of the equations (19) and (20), for the three measurement cycles, and the percentage differences Δ% from the estimated values to the literature ones.
Unlike the previous cases, the determination of the thermal conductivity values of the samples by the developed model was successful. The thermal conductivities estimated for the measured samples clearly present small divergences from the literature values, which show that the equations set for the sensors by means of calibration are fairly good.
The developed experimental model worked because the errors of the real measurement process were compensated by the calibration. As shown by equation (15), some errors of the measurement process come from the lack of some information, such as the thermal conductivity of the sphere, the right positioning of the temperature sensor within the sphere and the thermal resistance between the sphere and the medium. In addition, other phenomena may also cause errors in the process, such as the influence of the sensor cable and the possibility of some convective effects. That is why it is very important to calibrate each heated sphere sensor individually, so that a compensation for the unknown errors may be included in the ideal model.
The determination of the water content WC of the binary samples was done by means of the thermal conductivity measurements. Using the literature values (Tab. 3), a regression equation (a polynomial of third degree) was fitted correlating the water content to the thermal conductivity: (21)
Due to the measurement uncertainty of the water content, the measurement of alcohol containing low concentrations of water, such as anhydrous ethanol (about 0.7% water) is not suitable. Thus, a sample of ethyl alcohol (about 5% water) was used. Table 6 presents for the three measurement cycles the thermal conductivity and the water content estimated for the binary mixtures of water and ethanol.
Experimental ratios and ξ for each sample obtained at 20 °C.
Thermal conductivities of the samples by the developed model (in W/m °C) at 20 °C.
Water content of the binary mixtures of water and ethanol at 20 °C.
7 Uncertainty analysis
Based on the Guide to the Expression of Uncertainty in Measurement (GUM) [44], the combined standard uncertainties of the values of thermal conductivity and water content of the samples were respectively estimated as: (22) where

u(k) is the combined standard uncertainty of k;

u(r) is the standard uncertainty due to the thermal conductivity resolution adopted for the equipment (based on a rectangular distribution);

u(s_{k}) is the standard uncertainty of the mean thermal conductivity of the sample obtained from literatures (based on a normal distribution);

u(f_{k}) is the standard uncertainty related to the polynomial equations (19) and (20) (based on a normal distribution);

is the standard uncertainty related to the experimental ratio (based on a normal distribution), which combined the uncertainties of the individual components and ΔT;

c_{i}(k) is the coefficient sensitivity of the thermal conductivity to a particular input, estimated as the first order derivate of k with respect to the experimental ratio .

u(WC) is the combined standard uncertainty of WC;

u(r_{WC}) is the standard uncertainty due to the water content resolution adopted for the equipment (based on a rectangular distribution;

u(f_{WC}) is the standard uncertainty related to the polynomial equation (21) (based on a normal distribution);

u(k) is the combined standard uncertainty of k;

c_{i}(WC) is the coefficient sensitivity of the water content to a particular input, estimated as the first order derivate of WC with respect to k.
The expanded uncertainties of the values of thermal conductivity and water content of the samples were respectively estimated as: (24) (25) where U(k) is the expanded uncertainty of the thermal conductivity; κ is the coverage factor (κ equals two was used, which gives a coverage probability of 95.45%); and U(WC) is the expanded uncertainty of the water content.
It is worth pointing out that the reproducibility of ΔT is due to the reproducibility of the temperature measurements, which was estimated during the calibration of the thermometer (by means of the repetition of the calibration point 20 °C) and included in its uncertainty. The repeatability of ΔT was estimated during the temperature measurements of the samples as the standard deviation of the mean values (before and during the heating of the sphere). In the case of , the reproducibility consists on how the equipment reproduces the heat generation. It was estimated during the characterization of the equipment (by means of the repetition of three calibration points) and included in its uncertainty. The repeatability of was also estimated during the characterization of the equipment by means of the standard deviation of the mean values of voltage and current and included in its uncertainty.
Table 7 presents the standard uncertainties, the combined standard uncertainty, the expanded uncertainty and the percentage expanded uncertainty of the thermal conductivity determined for each sample.
Table 8 presents the standard uncertainties, the combined standard uncertainty and the expanded uncertainty of the water content determined for each sample.
Uncertainties for the thermal conductivities of the samples.
Uncertainties for the water content of the binary mixtures.
8 Comparing the results
The thermal conductivity values reported in some literatures were not provided with their uncertainty statements or error estimations. So, the agreement of the measured values with the literature ones was checked by means of the percentage expanded uncertainties of the thermal conductivities determined for the samples. The percentage differences between all the values of thermal conductivity given in the literatures (by means of tables or by the equations from (6) to (12)) and the mean thermal conductivity measured for each sample (calculated by means of the values of the three measurement cycles) were smaller than its percentage expanded uncertainty, which confirm the agreement of the measured values with the literature ones. Figures 3 (sensor 01) and 4 (sensor 02) show the difference between each literature value and the mean thermal conductivity, and the uncertainty of the thermal conductivity measurement of each sample.
It can be clearly seen by means of Figures 3 and 4 that there is a good agreement between the measured thermal conductivities and the literature ones for all the samples, which validates the developed method for measuring the thermal conductivity of water, ethanol and their mixtures.
In order to validate the water content determined for the binary samples, these were compared to the reference measurements performed at Inmetro and presented in Table 1. The compatibility of the measurements was checked by means of the normalized error (E_{n}), which is calculated according to the following equation [45]: (26) Table 9 shows the reference measurements of the water content of the binary samples, the values determined by the developed method in this study and the normalized errors.
A comparison between two measurements is satisfactory when E_{n} ≤ 1. Consequently, the water content measurements of this work and those performed by the reference laboratory are clearly compatible, since all the E_{n} values were lower than one, validating the developed method for measuring the water content of binary samples of water and ethanol.
Fig. 3 Differences between literature values and the mean thermal conductivity for sensor 01. 
Fig. 4 Differences between literature values and the mean thermal conductivity for sensor 02. 
Comparison of the measurements of WC of the binary samples (in % water).
9 Conclusions
The method of the spherical heat source, in principle, is an absolute measuring method of thermal conductivity, which means that the sensor can provide an output without being calibrated against a standard or a reference material. However, some parameters of the model need to be carefully taken into consideration. Thus, to compensate for the lack of some theoretical evaluations and the difficulty for obtaining accurately some practical information, the heated sphere sensors needed to be calibrated by means of mediums with known thermal properties. As the devices have different constructive characteristics from each other, the calibration must be done individually.
The use of spherical heat sources in the industrial sectors for measuring the thermal conductivity presents considerable advantages, such as wide measuring range, relatively fast measurements, measurement uncertainty compatible with other techniques and the possibility of using the sensor in situ. The method has also a great potential to be employed in research institutes and laboratories that provide calibration and testing services.
The applicability of the method of the spherical heat source for measuring the thermal conductivity of water, ethanol and their mixtures proved to be quite satisfactory, since the measurements of this study showed excellent agreement with the values proposed by several researchers. This agreement occurred with values obtained by other measurement methods, such as the traditional hotwire technique, as with those obtained by means of equations for estimating the thermal conductivity. The applicability of the method for determining the water content of the binary samples was also quite satisfactory, since the results of the proposed method showed good agreement with those performed by the reference equipment from the national institute (Inmetro).
10 Implications and influences
The paper presents the metrological aspects and the calibration procedure of a spherical heat source to measure thermal conductivity and water content of some liquid samples. This will strongly contribute and stimulate future works on the investigation of the applicability of the spherical heat source method and its metrological aspects for measuring thermal conductivity and water content of other liquids and other mediums.
References
 M.J. Assael, K.D. Antoniadis, W.A. Wakeham, Int. J. Thermophys. 31, 1051 (2010) [Google Scholar]
 ASTM D271715, Standard Test for Thermal Conductivity of Liquids (American Society for Testing and Materials, Philadelphia, PA, 2009) [Google Scholar]
 U. Hammerschmidt, W. Sabuga, Int. J. Thermophys. 21, 1255 (2000) [Google Scholar]
 F.A. Gibbs, Proc. Soc. Exp. Biol. Med. 31, 141 (1933) [Google Scholar]
 J. Grayson, J. Physiol. 114 (Suppl.), 29 (1951) [Google Scholar]
 W. Chester, J. Grayson, Nature 4273, 521 (1951) [Google Scholar]
 J. Grayson, J. Physiol. 118, 54 (1952) [CrossRef] [PubMed] [Google Scholar]
 H.S. Carslaw, in Introduction to the Mathematical Theory of the Conduction of Heat in Solids (Macmillan, London, 1921), 2nd ed. [Google Scholar]
 J.C. Chato, Therm. Prob. Biotechnol. Trans. ASME 16 (1968) [Google Scholar]
 T.A. Balasubramaniam, H.F. Bowman, J. Heat Trans. Trans. ASME 296 (1974) [CrossRef] [Google Scholar]
 R.K. Jain, J. Biomech. Eng. Trans. ASME 101, 82 (1979) [CrossRef] [Google Scholar]
 M.M. Chen, K.R. Holmes, V. Rupinskas, J. Biomech. Eng. Trans. ASME 103, 253 (1981) [Google Scholar]
 J.W. Valvano, The use of thermal diffusivity to quantify tissue perfusion, Ph.D. thesis, Harvard University, 1981 [Google Scholar]
 G. Hamilton, Investigation of the thermal properties of human and animal tissues, Ph.D. thesis, University of Glasgow, 1998 [Google Scholar]
 H. Zhang, S. Cheng, L. He, A. Zhang, Y. Zheng, D. Gao, Cell Preserv. Technol. 1, 141 (2002) [Google Scholar]
 K.A. Woodbury, Experimental and analytical investigation of liquid moisture distribution in roof insulating systems, Ph.D. thesis, Virginia Polytechnic Institute and State University, 1984 [Google Scholar]
 M. Fujii, H. Takamatsu, T. Fujii, in Proceedings of the 1st Asian Thermophysical Properties Conference, Beijing, 1986, edited by W. Buxuan et al. (China Academic Publishers, Beijing, 1986), pp. 462–467 [Google Scholar]
 B.P. Dougherty, An automated probe for thermal conductivity measurements, M.Sc. thesis, Virginia Polytechnic Institute and State University, 1987, cited in [22] [Google Scholar]
 R.R. Kravets, Determination of thermal conductivity of food materials using a bead thermistor, Ph.D. thesis, Virginia Polytechnic Institute and State University, 1988, cited in [22] [Google Scholar]
 U.B. Holeschovsky, G.T. Martin, J.W. Tester, Int. J. Heat Mass Transf. 39, 1135 (1996) [Google Scholar]
 S. Radhakrishnan, Measurement of thermal properties of seafood, M.Sc. thesis, Virginia Polytechnic Institute and State University, 1997 [Google Scholar]
 M.F. Gelder, A thermistor based method for the measurement of thermal conductivity and thermal diffusivity of moist food materials at high temperatures, Ph.D. thesis, Virginia Polytechnic Institute and State University, 1998 [Google Scholar]
 H. Zhang, L. He, S. Cheng, Z. Zhai, D. Gao, Meas. Sci. Technol. 14, 1396 (2003) [Google Scholar]
 L. Kubicar, L. Bagel, V. Vretenar, V. Stofanik, in Proceedings of the Meeting of the Thermophysical Society – Thermophysics 2005, Kocovce, 2005, edited by P. Matiasovsky, O. Koronthalyova (ICA SAS, Bratislava, 2005), pp. 38–42 [Google Scholar]
 L. Kubicar, V. Vretenar, V. Stofanik, V. Bohac, Int. J. Thermophys. (2008), doi:10.1007/s1076500805444 [Google Scholar]
 L. Kubicar, U. Hammerschmidt, D. Fridrikova, P. Dieska, V. Vretenar, in Proceedings of Thermophysics 2010, Valtice, 2010 (Brno University of Technology, Brno, 2010), pp. 166–171 [Google Scholar]
 J. Goldemberg, L.A.H. Nogueira, http://bioenergyconnection.org/article/sweeteningbiofuelsectorhistorysugarcaneethanolbrazil, accessed on May 2016 [Google Scholar]
 H.S. Carslaw, J.C. Jaeger, in Conduction of Heat in Solids (Oxford University Press, London, 1959), 2nd ed. [Google Scholar]
 H.F. Bowman, T.A. Balasubramaniam, Cryobiology 13, 572 (1976) [CrossRef] [PubMed] [Google Scholar]
 J.W. Valvano, J.R. Cochran, K.R. Diller, Int. J. Thermophys. 6, 301 (1985) [Google Scholar]
 R.J. Schmidt, S.W. Milverton, Proc. Roy. Soc. Lond. A: Mater. 152, 586 (1935) [CrossRef] [Google Scholar]
 M.J. Assael, E. Charitidou, W.A. Wakeham, Int. J. Thermophys. 10, 793 (1989) [Google Scholar]
 L.P. Filippov, Int. J. Heat Mass Transf. 11, 331 (1968) [Google Scholar]
 M.H. Khan, Modeling, simulation and optimization of ground source heat pump systems, M.Sc. thesis, Oklahoma State University, 2004 [Google Scholar]
 M.L.V. Ramires, C.A.N. de Castro, Y. Nagasaka, A. Nagashima, M.J. Assael, W.A. Wakeham, J. Phys. Chem. Ref. Data 24, 1377 (1995) [Google Scholar]
 S. Henke, P. Kadlec, Z. Bubník, J. Food Eng. 99, 497 (2010) [Google Scholar]
 Korea Thermophysical Properties Data Bank, KDB: Pure Components Properties (2012), available on http://www.cheric.org, accessed on 2012 [Google Scholar]
 IAPWS, Online Calculation for General and Scientific Use (2012), available on http://www.iapws.org/relguide/ThCond.html, accessed on 2012 [Google Scholar]
 IAPWS, Release on the IAPWS Formulation 2011 for the Thermal Conductivity of Ordinary Water Substance (2012), available on http://www.iapws.org/release.html, accessed on 2012 [Google Scholar]
 M. Tanaka, G. Girard, R. Davis, A. Peuto, N. Bignell, Metrologia 38, 301 (2001) [Google Scholar]
 L.Q. Fang, L.R. Sen, N.D. Yan, H.Y. Chun, J. Chem. Eng. Data 42, 971 (1997) [Google Scholar]
 A. Melinder, Thermophysical properties of aqueous solutions used as secondary working fluids, Ph.D. thesis, Royal Institute of Technology, 2007 [Google Scholar]
 B.E. Polin, J.M. Prausnitz, J.P. O'Connell, in The Properties of Gases and Liquids (McGrawHill, New York, 2001), 5th ed. [Google Scholar]
 INMETRO, Avaliação de Dados de Medição: Guia para a Expressão da Incerteza de Medição − GUM 2008 (INMETRO, Rio de Janeiro, 2012), 1st ed. [Google Scholar]
 Conformity Assessment − General Requirements for Proficiency Testing, ISO/IEC 17043 (International Organization for Standardization and International Electrotechnical Commission, 2010) [Google Scholar]
Cite this article as: Júlio Dutra Brionizio, Alcir de Faro Orlando, Georges Bonnier, Characterization of a spherical heat source for measuring thermal conductivity and water content of ethanol and water mixtures, Int. J. Metrol. Qual. Eng. 8, 18 (2017)
All Tables
Thermal conductivity values of the samples from literatures (in W/m °C) at 20 °C.
Thermal conductivities of the samples by the developed model (in W/m °C) at 20 °C.
All Figures
Fig. 1 Measuring system (left) and spherical sensors (right). 

In the text 
Fig. 2 Temperature monitoring during a measurement cycle. 

In the text 
Fig. 3 Differences between literature values and the mean thermal conductivity for sensor 01. 

In the text 
Fig. 4 Differences between literature values and the mean thermal conductivity for sensor 02. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.