Thermocouple emf reference functions
Python module containing calibration data and lookup functions for standard thermocouples of types B, C, D, E, G, J, K, M, N, P, R, S, T, and some less standard types too.
Usage and examples
Below, the first computation shows that the type K thermocouple emf at 42 °C, with reference junction at 0 °C, is 1.694 mV (compare to NIST table); the second calculation shows how passing in an array applies the function for each element, in the style of numpy:
>>> from thermocouples_reference import thermocouples >>> typeK = thermocouples['K'] >>> typeK <Type K thermocouple reference (-270.0 to 1372.0 degC)> >>> typeK.emf_mVC(42, Tref=0) 1.6938477049901346 >>> typeK.emf_mVC([-3.14159, 42, 54], Tref=0) array([-0.12369326, 1.6938477 , 2.18822176])
An inverse lookup function is provided that you can use to get a temperature out of a measured voltage, including cold junction compensation effects. If we put our type K thermocouple into a piece of spam and we read 1.1 mV, using our voltmeter at room temperature (23 °C), then the spam is at 50 °C. 
>>> typeK.inverse_CmV(1.1, Tref=23.0) 49.907928030075773 >>> typeK.emf_mVC(49.907928030075773, Tref=23.0) # check result 1.1000000000000001
The functions are called emf_mVC and inverse_CmV just to remind you about the units of voltage and temperature. Other temperature units are supported as well:
|Temperature unit||EMF lookup||Inverse lookup|
You can also compute derivatives of the emf function. These are functional derivatives, not finite differences. The Seebeck coefficients of chromel and alumel differ by 42.00 μV/°C, at 687 °C:
>>> typeK.emf_mVC(687,derivative=1) 0.041998175982382979
Readers may be familiar with thermocouple lookup tables (example table). Such tables are computed from standard reference functions, generally piecewise polynomials.  This module contains the source polynomials directly, and so in principle it is more accurate than any lookup table. Lookup tables also often also include approximate polynomials for temperature lookup based on a given compensated emf value. Such inverse polynomials are not included in this module; rather, the inverse lookup is based on numerically searching for a solution on the exact emf function.
For any thermocouple object, information about calibration and source is available in the repr() of the .func attribute:
>>> typeK.func <piecewise polynomial+gaussian, domain -270.0 to 1372.0 in degC, output in mV; ITS-90 calibrated, from NIST SRD 60, type K>
The data sources are:
- Types B, E, J, K, N, R, S, T use coefficients from NIST’s website, and are calibrations to the ITS-90 scale. 
- Types G, M, P, and non-lettered types Au/Pt, Au/Pd, AuFe 0.07, IrRh 40/0, PtMo 5/0.1, PtRh 40/20 use coefficients from ASTM E 1751-00 and are calibrations to ITS-90.
- Types C, D  use coefficients found from a publication of OMEGA Engineering Inc., and are calibrations to IPTS-68 scale. 
- scipy (optional, only needed for inverse lookup)
- python2 or python3 languages
Recommended installation is via pip. First, install pip. Then:
pip install thermocouples_reference --user
(Remove the --user option if you are superuser and want to install system-wide.)
This module is provided for educational purposes. For any real-world process, I strongly recommend that you check the output of this module against a known good standard.
I make no warranties as to the accuracy of this module, and shall not be liable for any damage that may result from errors or omissions.
|||This is the optimal temperature for spam. Always make sure your spam reads around 1.1 millivolt and you’ll have a tasty treat.|
|||A notable exception is NIST’s type K curve which uses a polynomial plus gaussian. The gaussian conveniently captures a wiggle in the Seebeck coefficient of alumel, that happens around 130 °C.|
|||The ITS-90 value T90 is believed to track the true thermodynamic temperature T very closely. The error T − T90 is quite small, of order 0.01 K for everyday conditions (up to about 200 °C), rising to around 0.05 K up at 1000 °C, and increasing even further after that. See Supplementary Information for the ITS-90. Generally your thermocouple accuracy will be more limited by manufacturing variations and by degradation of the metals in the thermal gradient region.|
|||An extra type G IPTS68 curve from the same source is available in thermocouples_reference.source_OMEGA.thermocouples. The type G in the main thermocouples_reference.thermocouples contains the ASTM curve which is ITS-90 calibrated.|
|||IPTS-68 reads higher than ITS-90 by about 1 °C at high temperatures around 2000 °C. See Supplementary Information for the ITS-90 (specifically Fig. 5 in the Introduction) for more information about the difference.|