Issue 
Int. J. Metrol. Qual. Eng.
Volume 10, 2019



Article Number  17  
Number of page(s)  8  
DOI  https://doi.org/10.1051/ijmqe/2019017  
Published online  12 December 2019 
Research Article
Modelling of thermal effects in a slim tin fixedpoint cell
^{1}
Temperature and Humidity Laboratory, Instituto Nacional de Metrología de Colombia − INM, Bogotá D.C 111321, Colombia
^{2}
Departamento de Ingeniería Química y Ambiental, Universidad Nacional de Colombia, Bogotá 11001, Colombia
^{*} Corresponding author: sacarvajal@inm.gov.co
Received:
21
October
2019
Accepted:
18
November
2019
Slim fixedpoint cells are an economic and suitable alternative for calibration of industrial platinum resistance thermometers, which are widely used in research and industrial applications. The most relevant disadvantage of slim cells is related with thermal effects produced by interactions between the cell and the furnace used for the realization of the fixedpoints. This work presents the implementation of a twodimensional numerical model that considers the transient heat transfer due to conduction, radiation and latent heat. The main goal was to evaluate different geometrical parameters in the design of slim fixedpoint cells. Results show that thermal effects are mitigated when the annular distance between interior and exterior walls of the crucible is maximum, while the length of the cell and the thermowell must be setup according to the temperature gradients in the furnace to avoid undesirable heat fluxes.
Key words: Slim fixedpoint / modeling / thermal effects / ITS90
© S.A. Carvajal et al., published by EDP Sciences, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://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
The objective of the National Metrology Institute of Colombia (INM) is to disseminate traceable measurements to the SI. In thermometry, this is done through calibration services using the fixedpoints of the ITS90. These fixedpoints are realized through the phase change of high purity materials such as mercury, water, gallium, indium, tin, zinc, etc. In the interval between the triple point of mercury (−38.8344 °C) and the freezing point of zinc (419.527 °C), the realization of the ITS90 requires the use of furnaces with uniformities of the order of tens of millikelvins, as well as considerable amounts of high purity metals (∼1 kg) in addition to highlydeveloped facilities and expert personnel. Most of the services required by industry and academy are in this interval, therefore, it is important to develop working standards to cover the traceability demands.
In order to increase the operating capacity and to obtain a better understanding of the phenomena involved in the realization of the ITS90, the Temperature and Humidity Laboratory of the INM has decided to construct a tin slim cell (231.928 °C). The main difference with traditional fixedpoint cells is that slim cells contain about a quarter of the typical amount of high purity metal, which implies that the realization in slim cells can be done in less demanding facilities. The consequence of the use of slim cells instead of traditional ones, is a poor thermal contact between the thermometer and the material undergoing the phase transition, which is evidenced as a high contributions of the thermal effects in the measurement uncertainty [1]. Due to this limitation, Standard Platinum Resistance Thermometers (SPRTs), which have the higher accuracy available are not calibrated using slim cells. However, in applications with lower accuracy requirements such as calibration of Industrial Platinum Resistance Thermometers (IPRTs) [2], calibration of thermocouples [3] or study of impurity effects [4], slim cells fit the purpose.
The thermal effects are a consequence of dynamic and static behaviors. Dynamic thermal effects are related with the heat capacity of the cellthermometer system and its geometry, while static thermal effects are related with the heat fluxes generated between the sensing element of the thermometer and the colder and warmer zones of the cell [5]. Regarding the static thermal effects, its two main contributions are the heat transfer between the sensing element and the furnace and the heat transfer between the sensing element and the external surroundings [6]. While the effect of the external surrounding can be reduced using a proper immersion of the thermometer in the cell, the interaction between the furnace and the cell can be more complex to manage, this is especially important in slim cells designed for tin due to the influence of the supercooling behavior in freezing curves [7] and impurity distributions in melting curves [8].
In order to find the optimal geometrical parameters that allow the minimization of the aforementioned effects, a twodimensional numerical simulation was implemented using the finiteelement method by means of the commercial software FLUENT^{®}. The coupled effects of phase transition, heat conduction and radiation heat transfer on transient mode were evaluated in the cell assembly as well as in the thermometer, mainly considering the influence on the sensing element.
2 Methodology
2.1 Geometry of the cell
The main components and the structure of a slim fixedpoint cell are presented in Figure 1a. The cell consists of a highpurity metal contained in a crucible, which is designed to enhance the formation of a liquid–solid interface around the thermowell, where a temperature sensor can be introduced. The crucible and the thermowell are assembled into a cell holder, which is sealed at the top, to ensure a constant pressure of 101 325 Pa inside the cell.
Figure 1b presents the geometrical parameters used in the design of the crucible. The parameter a corresponds to the annular distance between interior and exterior walls of the crucible, b is the length of high purity material contained in the crucible and d is the distance between the bottom of the thermowell and the bottom of the crucible.
The fixedpoint cell was modeled as a twodimensional system with axial symmetry, as shown in Figure 2b and c. The cell was divided into two sections and eight subdomains distributed as follows:
Section 1 : Cell (Fig. 2b)

Filling gas

Cell holder

Thermowell

High purity tin

Crucible
Section 2: Thermometer (Fig. 2c)

Thermometer cover

Thermometer filling gas

Sensing element
Fig. 1 Fixedpoint cell scheme, (a) main components, (b) geometrical parameters of the crucible. 
Fig. 2 Fixedpoint cell assembly model, (a) main components, (b) cell model, (c) thermometer model. 
2.2 Materials properties
The crucible material was simulated as highpurity graphite, the filling gas as highpurity argon and the thermowell and the cell holder as stainless steel 303. The thermometer was modeled as a platinum wire inside a hollow Inconel sheathing containing argon as filling gas. The thermophysical properties of the materials used are presented in Table 1, potential impurities were not considered.
Thermophysical properties of the materials used.
2.3 Numerical model
The transient energy balance over the measurement system can be written as:(1)where ρ is the density, H is the enthalpy, k is the thermal conductivity, T is the temperature and q_{R} is the radiative heat flux. The enthalpy in equation (1) can be expressed as: (2)
where the term h represents the contribution due to the sensible heat and ΔH represents the contribution due to the latent heat. The last one is associated with the highpurity metal and was computed through the enthalpyporosity technique proposed by Brent et al. [19]. This technique is especially useful for solving diffusion phasechange problems through fixedgrid methods. The cells of the fixedgrid in the region undergoing the phase transition are numerically modeled as a pseudo porous media, where the porosity is equal to zero in the solid phase and equal to unity in the liquid phase. The porosity β relates ΔH to the latent heat of the material L through equation (3), and can be expressed as a linear function of temperature, according to equation (4):(3) (4)where T_{s} is the temperature when all the material is solid and T_{l} is the temperature when all the material is molten.
The radiative heat transfer was estimated through the diffusegray approximation and assuming nonparticipating media. In this case, the radiative component can be calculated with equation (5).(5)where ε is the emissivity, σ is the SteffanBoltzman constant and J is the total heat flux leaving the surface, often called radiosity. The radiosity is calculated from the view factor F at each surface i, according to equation (6) [20].(6)
2.4 Thermal effects
Ideally, in the realization of a fixedpoint, the sensing element of the thermometer is completely surrounded by the liquid–solid interface, as shown in Figure 3a. Due to the heat transfer between cell, furnace and environment, the interface can be moved from the ideal position or even disappear in some regions, promoting the development of thermal bridges between the sensing element and the furnace, as shown in Figure 3b. These phenomena are known as “thermal effects” and along with impurities, they represent the major contributions to the uncertainty when realizing the fixedpoints [21].
In order to evaluate the influence of thermal effects on different geometrical configurations of the crucible, the penetration ratio µ was calculated using the methodology suggested by Fahr and Rudtsch [6] through equation (7):(7)where ΔT_{meas} is the change in temperature measured in the cell due to the change in the furnace temperature ΔT_{i} . If the parameter µ is equal to 0, then there is a proper fixedpoint realization, i.e. a total absence of thermal effects. If µ is equal to 1, it means the worst possible realization, where the effect of the wall temperature dominates the measurement. The penetration ratio as a function of time is presented in Figure 4, where the red curve is the profile assigned to the furnace temperature and the blue curve is the temperature measured in the cell.
Fig. 3 Representation of thermal effects in a fixedpoint cell. (a) Proper liquid–solid interface. (b) Discontinuous liquid–solid interface. 
Fig. 4 Penetration ratio µ in a fixedpoint realization. Cell temperature (blue) and furnace temperature of the cell (red). 
2.5 Boundary conditions
The boundary conditions were imposed on the cell holder and on the lid, as depicted in Figure 5. In order to calculate the penetration ratio µ, the wall temperature on the cell holder was established as a timedependent rectangular function, as shown in Figure 4 (red curve) varying the temperatures between 504.778 K and 504.978 K. On the lid of the cell a convective boundary condition was imposed.
For different configuration of the parameters a, b and d of Figure 1, the coupled equations (1)–(6) were solved using the finiteelement method implemented in the software FLUENT^{®} with the boundary conditions defined in Figure 5.
Fig. 5 Fixed point cell boundary conditions. 
3 Results and discussion
Figures 6–8 present the penetration ratio for different geometrical configurations for solid fraction F between 0.6 and 1. The different combinations of parameters a, b and d have the restriction of producing the same amount of highpurity metal inside the crucible in order for the results to be comparable. As pointed out by Krizmanic et al. [22] the coupled calculations of conductive and radiative heat transfer implemented in FLUENT^{®} can produce errors in the order of 0.1 mK in the relative position of the plateau, however, the results of the simulation are useful in analyzing tendencies. It is also important to note that the model represents only the macroscopic effects and it did not consider microscopic effects related with crystal growth and nucleation, which are important in the analysis of initiation techniques [23].
Figure 6 shows the variation of the penetration ratio with the solid fraction at different b/d ratios, keeping a constant. All the penetration ratios present an increasing trend at b/d = 8 and b/d = 10, forming a flat region from F = 0.82 and F = 0.95, respectively. At b/d = 8, the penetration ratio is larger with respect to the other b/d ratios from F = 0.7. At b/d = 10, the penetration ratio is larger than b/d ratios between 18 and 35 from F = 0.85. The flat regions indicate that thermal effects dominate the process, therefore it can be concluded that thermal effects begin to be significant at lower F values when the b/d ratio decreases between 8 and 10, however, between 18 and 35 b/d ratios, the thermal effects are slightly larger for the highest b/d ratio. This indicates that there is an optimal b/d ratio that minimizes the penetration ratio. The marked effect of the b/d ratio on the realization of the fixedpoint is related to the heat flux through the sensing element of the thermometer. In short cells with excessive volume of highpurity metal below the bottom of the thermowell and a small volume in the surrounding area of the thermometer, i.e. low b/d ratios, the liquid–solid interface along the thermometer is not enough to guarantee a uniform temperature distribution in the sensing element.
Figure 9 illustrates the temperature profiles across the sensing element of the thermometer in the fixedpoint assembly for the b/d ratios of Figure 6 with solid fractions of 0.6, 0.8 and 0.95. At F = 0.6, there are no temperature gradients around the sensing element, indicating that thermal effects from the furnace are low. When F = 0.8, there is a displacement of the isothermal zone around the sensing element. At lower b/d ratios, the liquid–solid interface goes to the bottom of the crucible, uncovering the tip of the sensing element, which increases the thermal effects. At the end of the plateau, using an F = 0.95, the previous pattern is more evident. At b/d = 8 and b/d = 10, the sensing element is far from the isothermal zone, while at b/d = 35, the liquid–solid interface moves to the top of the sensing element uncovering the tip. At b/d = 18 and F = 0.95, the isothermal zone is smaller than b/d = 35, but the liquid–solid interface is better distributed around the sensing element, which explains the behavior of the penetration ratio when F is close to 1. It is important to note that, at b/d = 35 and F = 0.95, a higher temperature on the bottom of the furnace will produce an isothermal zone more uniform around the sensing element. As pointed out by Pearce et al. [24] there is an optimum temperature distribution in the axial direction of the furnace that minimizes thermal effects, which means that the appropriate b/d ratio will depend on the vertical profile of the furnace.
Figure 7 shows the variation of the penetration ratio with the solid fraction at different b/a ratios, keeping d constant. As indicated by the presence of flat regions, thermal effects become dominant at lower F values when the b/a ratio increases, for instance, at b/a = 17 the flat region begins at F = 0.94 while at b/a = 44 the flat region begins at F = 0.86. It also shows that even at high b values, the effect of the annular distance a is relevant for a proper realization of the fixedpoint, since at low annular distances there is not enough space for the formation of a uniform phase boundary, which in turn promotes the development of thermal bridges between the furnace and the sensing element of the thermometer, as mentioned by Fahr et al. [5].
Figure 8 shows the variation of the penetration ratio as a function of the solid fraction at different d/a ratios, keeping b constant. Once again, low annular distances impair the formation of a uniform phase boundary as explained previously. At d/a = 8.7, the thermal effects dominate the realization even for a solid fraction as low as 0.6. Consequently, at lower values of the d/a ratio, the thermal effects become less important, as shown at d/a = 1.7.
As observed in Figures 6–8, especially with an F > 0.8, lower penetration ratios, i.e. lower thermal effects, are obtained at higher b/d ratios and at lower b/a and d/a ratios; this means that the design of slim cells should have the maximum a distance, keeping a b/d ratio high enough to assure a uniform liquid–solid interface near to the sensing element of the thermometer, considering the vertical temperature profile of the furnace.
Fig. 6 Penetration ratio µ as a function of the solid fraction F for different b/d ratios at a constant. 
Fig. 7 Penetration ratio µ as a function of the solid fraction F for different b/a ratios at d constant. 
Fig. 8 Penetration ratio µ against solid fraction F for different d/a ratios at b constant. 
Fig. 9 Temperature profiles of different b/d ratios for solid fraction F between 0.6 and 1 (white sections represent temperatures lower than 504.9 K). 
4 Conclusions
The modeling of thermal effects due to furnacecell interaction in a slim cell was used to evaluate different geometrical designs for the realization of the freezing point of tin. The fixedpoint cell assembly was modeled as a twodimensional axisymmetric system that considers transient conductiveradiative heat transfer and the latent heat released during the freezing process through the diffusegray approximation and the enthalpyporosity method, respectively. The evaluation of thermal effects was done through the penetration ratio, which relates the variation in the temperature of the furnace to the temperature measured in the cell.
The results indicate that the geometrical parameters of the crucible can be optimized in order to minimize the penetration ratio and consequently reduce thermal effects. It was found that the annular distance between interior and exterior walls of the crucible should be as long as it can be in order to guarantee the proper formation of a phase boundary that avoids generation of thermal bridges between the furnace and the sensing element of the thermometer. The length of the cell and the thermowell should be established in accordance with the vertical temperature profile of the furnace, since the form and location of the liquid–solid interface depends on the thermal profile developed inside the cell.
References
 G.F. Strouse, Small fixedpoint cells for use in dry well block calibrators, in Proceedings of TEMPMEKO , 2001, pp. 783–788 [Google Scholar]
 X. Li, D. Farley, D. Chen, M. Zhao, M. Hirst, Mini metalcased fixedpoint cells, in Proceedings of 8th International Symposium on Temperature and Thermal Measurements in Industry and Science , 2001, 2, pp. 777–782 [Google Scholar]
 F. Edler, S. Haupt, S.A. Mokdad, G. Failleau, M. Sadli, Investigation of selfvalidating thermocouples with integrated fixedpoint units, Int. J. Metrol. Qual. Eng. 6 , 103 (2015) [CrossRef] [EDP Sciences] [Google Scholar]
 M. Fahr, S. Rudtsch, Oxides in metal fixed points of the ITS90, Metrologia 46, 423–438 (2009) [Google Scholar]
 M. Fahr, S. Rudtsch, A.K. Gerlitzke, Fixedpoint comparison uncertainties for two cell geometries, Int. J. Thermophys. 32, 2269–2280 (2011) [Google Scholar]
 M. Fahr, S. Rudtsch, A new method for the quantification and correction of thermal effects on the realization of fixed points, Int. J. Thermophys. 29, 126–138 (2008) [Google Scholar]
 M. RagayEnot, Y.H. Lee, Y.G. Kim, Fabrication of a mini multifixedpoint cell for the calibration of industrial platinum resistance thermometers, Meas. Sci. Technol. 28 , 075007 (2017) [Google Scholar]
 E.V. Vasil'ev, A.A. Ignatov, M.I. Kalinin, L.D. Markin, Investigation of the reproducibility of melting and freezing curves for tin in a miniature cell, Meas. Tech. 61, 494–500 (2018) [CrossRef] [Google Scholar]
 G.G. Graf, Tin, Tin alloys, and Tin Compounds, Ullmann's encyclopedia of industrial chemistry (WileyVCH Verlag GmbH & Co. KGaA, Weinheim, 2012) [Google Scholar]
 Y.S. Touloukian, R.W. Powell, C.Y. Ho, P.G. Klemens, Thermal conductivity: metallic elements and alloys (IFI/Plenum, New York, NY, 1970) [CrossRef] [Google Scholar]
 H.O. Pierson, Handbook of carbon, graphite, diamonds and fullerenes: processing, properties and applications (Noyes Publications, New Jersey, 1993) [Google Scholar]
 R. Taylor, The thermal conductivity of pyrolytic graphite, Philos. Mag. 13, 157–166 (1966) [CrossRef] [Google Scholar]
 A.T.D. Butland, R.J. Maddison, The specific heat of graphite: an evaluation of measurements, J. Nucl. Mater. 49, 45–56 (1973) [Google Scholar]
 P.D. Harvey, Engineering properties of steel (American Society for Metals, Ohio, 1982) [Google Scholar]
 Y.S. Touloukian, D.P. DeWitt, Thermal radiative properties: metallic elements and alloys (Springer, New York, 1970) [CrossRef] [Google Scholar]
 H. Renner et al., Platinum group metals and compounds, Ullmann's encyclopedia of industrial chemistry (WileyVCH GmbH, Weinheim, 2000) [Google Scholar]
 H. Yokokawa, Y. Takahashi, Laserflash calorimetry II. Heat capacity of platinum from 80 to1000 K and its revised thermodynamic functions, J. Chem. Thermodyn. 11, 411–420 (1979) [Google Scholar]
 W.F. Gale, T.C. Totemeier, Smithells metals reference book (Elsevier Science, Amsterdam, 2003) [Google Scholar]
 A.D. Brent, V.R. Voller, K.T.J. Reid, Enthalpyporosity technique for modeling convectiondiffusion phase change: Application to the melting of a pure metal, Numer. Heat Transf. Part A Appl. 13, 297–318 (1988) [Google Scholar]
 M.F. Modest, Radiative heat transfer (Elsevier Science, Amsterdam, 2013) [Google Scholar]
 J.V. Pearce et al., Guide to the Realization of the ITS90 Metal Fixed Points for Contact Thermometry, BIPM, 2018 [Google Scholar]
 S. Krizmanic, T. Veliki, D. Zvizdic, Modeling of transient heat transfer in zinc fixedpoint cell, Int. J. Thermophys. 32, 326–336 (2011) [Google Scholar]
 D.R. White, R.S. Mason, Improved initiation technique for the metal fixed points, Int. J. Thermophys. 32, 348–359 (2011) [Google Scholar]
 J.V. Pearce, R.I. Veltcheva, D.H. Lowe, Z. Malik, J.D. Hunt, Optimization of SPRT measurements of freezing in a zinc fixedpoint cell, Metrologia 49 , 359 (2012) [Google Scholar]
Cite this article as: Sergio A. Carvajal, J. Herney Ramírez, Andrés J. Bohórquez, Ciro A. Sánchez, Modelling of thermal effects in a slim tin fixedpoint cell, Int. J. Metrol. Qual. Eng. 10, 17 (2019)
All Tables
All Figures
Fig. 1 Fixedpoint cell scheme, (a) main components, (b) geometrical parameters of the crucible. 

In the text 
Fig. 2 Fixedpoint cell assembly model, (a) main components, (b) cell model, (c) thermometer model. 

In the text 
Fig. 3 Representation of thermal effects in a fixedpoint cell. (a) Proper liquid–solid interface. (b) Discontinuous liquid–solid interface. 

In the text 
Fig. 4 Penetration ratio µ in a fixedpoint realization. Cell temperature (blue) and furnace temperature of the cell (red). 

In the text 
Fig. 5 Fixed point cell boundary conditions. 

In the text 
Fig. 6 Penetration ratio µ as a function of the solid fraction F for different b/d ratios at a constant. 

In the text 
Fig. 7 Penetration ratio µ as a function of the solid fraction F for different b/a ratios at d constant. 

In the text 
Fig. 8 Penetration ratio µ against solid fraction F for different d/a ratios at b constant. 

In the text 
Fig. 9 Temperature profiles of different b/d ratios for solid fraction F between 0.6 and 1 (white sections represent temperatures lower than 504.9 K). 

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.