Physics:Steinhart–Hart equation

From HandWiki
Short description: Semiconductor resistance model

The Steinhart–Hart equation is a model of the resistance of a semiconductor at different temperatures. The equation is

[math]\displaystyle{ \frac{1}{T} = A + B \ln R + C (\ln R)^3, }[/math]

where

[math]\displaystyle{ T }[/math] is the temperature (in kelvins),
[math]\displaystyle{ R }[/math] is the resistance at [math]\displaystyle{ T }[/math] (in ohms),
[math]\displaystyle{ A }[/math], [math]\displaystyle{ B }[/math], and [math]\displaystyle{ C }[/math] are the Steinhart–Hart coefficients, which vary depending on the type and model of thermistor and the temperature range of interest.

Uses of the equation

The equation is often used to derive a precise temperature of a thermistor, since it provides a closer approximation to actual temperature than simpler equations, and is useful over the entire working temperature range of the sensor. Steinhart–Hart coefficients are usually published by thermistor manufacturers.

Where Steinhart–Hart coefficients are not available, they can be derived. Three accurate measures of resistance are made at precise temperatures, then the coefficients are derived by solving three simultaneous equations.

Inverse of the equation

To find the resistance of a semiconductor at a given temperature, the inverse of the Steinhart–Hart equation must be used. See the Application Note, "A, B, C Coefficients for Steinhart–Hart Equation".

[math]\displaystyle{ R = \exp\left(\sqrt[3]{y - x/2} - \sqrt[3]{y + x/2}\right), }[/math]

where

[math]\displaystyle{ \begin{align} x &= \frac{1}{C}\left(A - \frac{1}{T}\right), \\ y &= \sqrt{\left(\frac{B}{3C}\right)^3 + \frac{x^2}{4}}. \end{align} }[/math]

Steinhart–Hart coefficients

To find the coefficients of Steinhart–Hart, we need to know at-least three operating points. For this, we use three values of resistance data for three known temperatures.

[math]\displaystyle{ \begin{bmatrix} 1 & \ln R_1 & \ln^3 R_1 \\ 1 & \ln R_2 & \ln^3 R_2 \\ 1 & \ln R_3 & \ln^3 R_3 \end{bmatrix}\begin{bmatrix} A \\ B \\ C \end{bmatrix} = \begin{bmatrix} \frac{1}{T_1} \\ \frac{1}{T_2} \\ \frac{1}{T_3} \end{bmatrix} }[/math]

With [math]\displaystyle{ R_1 }[/math], [math]\displaystyle{ R_2 }[/math] and [math]\displaystyle{ R_3 }[/math] values of resistance at the temperatures [math]\displaystyle{ T_1 }[/math], [math]\displaystyle{ T_2 }[/math] and [math]\displaystyle{ T_3 }[/math], one can express [math]\displaystyle{ A }[/math], [math]\displaystyle{ B }[/math] and [math]\displaystyle{ C }[/math] (all calculations):

[math]\displaystyle{ \begin{align} L_1 &= \ln R_1, \quad L_2 = \ln R_2, \quad L_3 = \ln R_3 \\ Y_1 &= \frac{1}{T_1}, \quad Y_2 = \frac{1}{T_2}, \quad Y_3 = \frac{1}{T_3} \\ \gamma_2 &= \frac{Y_2 - Y_1}{L_2 - L_1}, \quad \gamma_3 = \frac{Y_3 - Y_1}{L_3 - L_1} \\ \Rightarrow C &= \left( \frac{ \gamma_3 - \gamma_2 }{ L_3 - L_2} \right) \left(L_1 + L_2 + L_3\right)^{-1} \\ \Rightarrow B &= \gamma_2 - C \left(L_1^2 + L_1 L_2 + L_2^2\right) \\ \Rightarrow A &= Y_1 - \left(B + L_1^2 C\right) L_1 \end{align} }[/math]

Developers of the equation

The equation is named after John S. Steinhart and Stanley R. Hart who first published the relationship in 1968.[1] Professor Steinhart (1929–2003), a fellow of the American Geophysical Union and of the American Association for the Advancement of Science, was a member of the faculty of University of Wisconsin–Madison from 1969 to 1991.[2] Dr. Hart, a Senior Scientist at Woods Hole Oceanographic Institution since 1989 and fellow of the Geological Society of America, the American Geophysical Union, the Geochemical Society and the European Association of Geochemistry,[3] was associated with Professor Steinhart at the Carnegie Institution of Washington when the equation was developed.

Derivation and alternatives

The most general form of the equation can be derived from extending the B parameter equation to an infinite series:

[math]\displaystyle{ R = R_0 e^{B\left(\frac{1}{T} - \frac{1}{T_0}\right)} }[/math]
[math]\displaystyle{ \frac{1}{T} = \frac{1}{T_0} + \frac{1}{B} \left(\ln \frac{R}{R_0}\right) = a_0 + a_1 \ln \frac{R}{R_0} }[/math]
[math]\displaystyle{ \frac{1}{T} = \sum_{n=0}^\infty a_n \left(\ln \frac{R}{R_0}\right)^n }[/math]

[math]\displaystyle{ R_0 }[/math] is a reference (standard) resistance value. The Steinhart–Hart equation assumes [math]\displaystyle{ R_0 }[/math] is 1 ohm. The curve fit is much less accurate when it is assumed [math]\displaystyle{ a_2=0 }[/math] and a different value of [math]\displaystyle{ R_0 }[/math] such as 1 kΩ is used. However, using the full set of coefficients avoids this problem as it simply results in shifted parameters.[4]

In the original paper, Steinhart and Hart remark that allowing [math]\displaystyle{ a_2 \neq 0 }[/math] degraded the fit.[1] This is surprising as allowing more freedom would usually improve the fit. It may be because the authors fitted [math]\displaystyle{ 1/T }[/math] instead of [math]\displaystyle{ T }[/math], and thus the error in [math]\displaystyle{ T }[/math] increased from the extra freedom.[5] Subsequent papers have found great benefit in allowing [math]\displaystyle{ a_2 \neq 0 }[/math].[6]

The equation was developed through trial-and-error testing of numerous equations, and selected due to its simple form and good fit.[1] However, in its original form, the Steinhart–Hart equation is not sufficiently accurate for modern scientific measurements. For interpolation using a small number of measurements, the series expansion with [math]\displaystyle{ n=4 }[/math] has been found to be accurate within 1 mK over the calibrated range. Some authors recommend using [math]\displaystyle{ n=5 }[/math].[6] If there are many data points, standard polynomial regression can also generate accurate curve fits. Some manufacturers have begun providing regression coefficients as an alternative to Steinhart–Hart coefficients.[7]

References

  1. 1.0 1.1 1.2 John S. Steinhart, Stanley R. Hart, Calibration curves for thermistors, Deep-Sea Research and Oceanographic Abstracts, Volume 15, Issue 4, August 1968, Pages 497–503, ISSN 0011-7471, doi:10.1016/0011-7471(68)90057-0.
  2. "Memorial resolution of the faculty of the University of Wisconsin-Madison on the death of professor emeritus John S. Steinhart". University of Wisconsin. 5 April 2004. http://www.secfac.wisc.edu/senate/2004/0405/1775(mem_res).pdf. 
  3. "Dr. Stan Hart". Woods Hole Oceanographic Institution. http://www.whoi.edu/science/GG/people/shart/index.htm. Retrieved 2 July 2015. 
  4. Matus, Michael (October 2011). "Temperature Measurement in Dimensional Metrology – Why the Steinhart–Hart Equation works so well". MacroScale 2011. Wabern, Switzerland. https://oar.ptb.de/files/download/810.20130620D.pdf. 
  5. Hoge, Harold J. (1 June 1988). "Useful procedure in least squares, and tests of some equations for thermistors". Review of Scientific Instruments 59 (6): 975–979. doi:10.1063/1.1139762. ISSN 0034-6748. https://aip.scitation.org/doi/10.1063/1.1139762. 
  6. 6.0 6.1 Rudtsch, Steffen; von Rohden, Christoph (1 December 2015). "Calibration and self-validation of thermistors for high-precision temperature measurements" (in en). Measurement 76: 1–6. doi:10.1016/j.measurement.2015.07.028. ISSN 0263-2241. https://www.sciencedirect.com/science/article/pii/S0263224115003632. Retrieved 8 July 2020. 
  7. "Comments on the Steinhart–Hart Equation". Building Automation Products Inc.. 11 November 2015. http://www.bapihvac.com/wp-content/uploads/SHH_Equation_Comments.pdf. Retrieved 8 July 2020. 

External links