Compare commits
18 Commits
f153939ecc
...
master
Author | SHA1 | Date | |
---|---|---|---|
12f04dbab4 | |||
d703e433d1 | |||
218db7ead6 | |||
96030ac4ff | |||
8737f39b1a | |||
b89355111e | |||
7ed730357c | |||
c1e099321d | |||
d445d36f0d | |||
ba71eeb29c | |||
553a14fdc7 | |||
802a5f97a1 | |||
cf4ed75dfe | |||
2ff9fc5839 | |||
99b7d16746 | |||
113065d23b | |||
970b6679f5 | |||
f538cca79c |
@@ -85,7 +85,7 @@ den = coeffs(D, z);
|
||||
```
|
||||
|
||||
```text
|
||||
num = (Ts^2*wn^2 - 4*Ts*gc*wn*xi + 4) + (2*Ts^2*wn^2 - 8) * z + (Ts^2*wn^2 + 4*Ts*gc*wn*xi + 4) * z^2
|
||||
org_babel_eoe
|
||||
```
|
||||
|
||||
```text
|
||||
@@ -272,7 +272,7 @@ And the transfer function is equal to `gain * num/den`.
|
||||
|
||||
Once the analytical formula of a discrete transfer function is obtained, it is possible to vary some parameters in real time.
|
||||
|
||||
This is easily done in Simulink (see Figure [1](#figure--fig:variable-controller-simulink)) where a `Discrete Varying Transfer Function` block is used.
|
||||
This is easily done in Simulink (see [Figure 1](#figure--fig:variable-controller-simulink)) where a `Discrete Varying Transfer Function` block is used.
|
||||
The coefficients are simply computed with a Matlab function.
|
||||
|
||||
<a id="figure--fig:variable-controller-simulink"></a>
|
||||
|
@@ -12,10 +12,11 @@ Tags
|
||||
|
||||
General purpose / PLC:
|
||||
|
||||
| Manufacturer |
|
||||
|-------------------------------------------------------------------------------------|
|
||||
| [Bechoff](https://www.beckhoff.com/fr-fr/products/i-o/ethercat-terminals/) |
|
||||
| [Wago](https://www.wago.com/global/i-o-systems/fieldbus-coupler-ethercat/p/750-354) |
|
||||
| Manufacturer |
|
||||
|----------------------------------------------------------------------------------------------|
|
||||
| [Bechoff](https://www.beckhoff.com/fr-fr/products/i-o/ethercat-terminals/) |
|
||||
| [Wago](https://www.wago.com/global/i-o-systems/fieldbus-coupler-ethercat/p/750-354) |
|
||||
| [Rexroth](https://apps.boschrexroth.com/microsites/ctrlx-automation/en/portfolio/ctrlx-i-o/) |
|
||||
|
||||
Acquisition systems:
|
||||
|
||||
|
215
content/zettels/heat_transfer.md
Normal file
@@ -0,0 +1,215 @@
|
||||
+++
|
||||
title = "Heat Transfer"
|
||||
author = ["Dehaeze Thomas"]
|
||||
draft = false
|
||||
+++
|
||||
|
||||
Tags
|
||||
:
|
||||
|
||||
|
||||
## Electrical Analogy - Lumped Mass Modeling {#electrical-analogy-lumped-mass-modeling}
|
||||
|
||||
The difference in temperature \\(\Delta T\\) is driving potential energy flow \\(Q\\) (in watts):
|
||||
|
||||
\begin{equation}
|
||||
\Delta T = R\_{th} \cdot Q
|
||||
\end{equation}
|
||||
|
||||
\\(R\_{th}\\) is the analogy of a "thermal resistance", and is expressed in K/W.
|
||||
|
||||
|
||||
## Conduction (diffusion) {#conduction--diffusion}
|
||||
|
||||
The _conduction_ corresponds to the heat transfer \\(Q\\) (in watt) through molecular agitation within a material.
|
||||
|
||||
\begin{equation}
|
||||
R\_{th} = \frac{d}{\lambda A}
|
||||
\end{equation}
|
||||
|
||||
with:
|
||||
|
||||
- \\(\lambda\\) the thermal conductivity in \\([W/m \cdot K]\\)
|
||||
- \\(A\\) the surface area in \\([m^2]\\)
|
||||
- \\(d\\) the length of the barrier in \\([m]\\)
|
||||
|
||||
|
||||
## Convection {#convection}
|
||||
|
||||
The convection corresponds to the heat transfer \\(Q\\) through flow of a fluid.
|
||||
It can be either _natural_ or _forced_.
|
||||
|
||||
\begin{equation}
|
||||
R\_{th} = \frac{1}{h A}
|
||||
\end{equation}
|
||||
|
||||
with:
|
||||
|
||||
- \\(h\\) the convection heat transfer coefficient in \\([W/m^2 \cdot K]\\).
|
||||
\\(h \approx 10.5 - v + 10\sqrt{v}\\) with \\(v\\) the velocity of the object through the fluid in \\([m/s]\\)
|
||||
Typically:
|
||||
- \\(h = 5 - 10\ W/m^2/K\\) for free convection with air
|
||||
- \\(h = 500 - 5000\ W/m^2/K\\) for forced water cooling in a tube of 5mm diameter
|
||||
- \\(A\\) the surface area in \\([m^2]\\)
|
||||
|
||||
Note that clean-room air flow should be considered as forced convection, and \\(h \approx 10 W/m^2/K\\).
|
||||
|
||||
|
||||
## Radiation {#radiation}
|
||||
|
||||
_Radiation_ corresponds to the heat transfer \\(Q\\) (in watt) through the emission of electromagnetic waves from the emitter to its surroundings.
|
||||
|
||||
In the general case, we have:
|
||||
\\[ Q = \epsilon \cdot \sigma \cdot A \cdot (T\_r^4 - T\_s^4) \\]
|
||||
with:
|
||||
|
||||
- \\(\epsilon\\) the emissivity which corresponds to the ability of a surface to emit energy through radiation relative to a black body surface at equal temperature.
|
||||
It is between 0 (no emissivity) and 1 (maximum emissivity)
|
||||
- \\(\sigma\\) the Stefan-Boltzmann constant: \\(\sigma = 5.67 \cdot 10^{-8} \\, \frac{W}{m^2 K^4}\\)
|
||||
- \\(T\_r\\) the temperature of the emitter in \\([K]\\)
|
||||
- \\(T\_s\\) the temperature of the surrounding in \\([K]\\)
|
||||
|
||||
In order to use the lumped mass approximation, the equations can be linearized to obtain:
|
||||
|
||||
\begin{equation}
|
||||
R\_{th} = \frac{1}{h\_{rad} A}
|
||||
\end{equation}
|
||||
|
||||
with:
|
||||
|
||||
- \\(h\_{rad}\\) the effective heat transfer coefficient for radiation in \\(W/m^2 \cdot K\\)
|
||||
- \\(A\\) the surface in \\([m^2]\\)
|
||||
|
||||
|
||||
### Practical Cases {#practical-cases}
|
||||
|
||||
Two parallel plates:
|
||||
|
||||
\begin{equation}
|
||||
h\_{rad} = \frac{\sigma}{1/\epsilon\_1 + 1/\epsilon\_2 - 1} (T\_1^2 + T\_2^2)(T\_1 + T\_2)
|
||||
\end{equation}
|
||||
|
||||
Two concentric cylinders:
|
||||
|
||||
\begin{equation}
|
||||
h\_{rad} = \frac{\sigma}{1/\epsilon\_1 + r\_1/r\_2 (1/\epsilon\_2 - 1)} (T\_1^2 + T\_2^2)(T\_1 + T\_2)
|
||||
\end{equation}
|
||||
|
||||
A small object enclosed in a large volume:
|
||||
|
||||
\begin{equation}
|
||||
h\_{rad} = \epsilon\_1 \sigma (T\_1^2 + T\_2^2)(T\_1 + T\_2)
|
||||
\end{equation}
|
||||
|
||||
|
||||
### Emissivity {#emissivity}
|
||||
|
||||
The emissivity of materials highly depend on the surface finish (the more polished, the lower the emissivity).
|
||||
Some examples are given in <tab:emissivity_examples>.
|
||||
|
||||
Gold coating gives also a very low emissivity and is typically used in cryogenic applications.
|
||||
|
||||
<a id="table--tab:emissivity-examples"></a>
|
||||
<div class="table-caption">
|
||||
<span class="table-number"><a href="#table--tab:emissivity-examples">Table 1</a>:</span>
|
||||
Some examples of emissivity (specified at 25 degrees)
|
||||
</div>
|
||||
|
||||
| Substance | Emissivity |
|
||||
|----------------------------|------------|
|
||||
| Silver (polished) | 0.005 |
|
||||
| Silver (oxidized) | 0.04 |
|
||||
| Stainless Steel (polished) | 0.02 |
|
||||
| Aluminium (polished) | 0.02 |
|
||||
| Aluminium (oxidized) | 0.2 |
|
||||
| Aluminium (anodized) | 0.9 |
|
||||
| Copper (polished) | 0.03 |
|
||||
| Copper (oxidized) | 0.87 |
|
||||
|
||||
<div class="exampl">
|
||||
|
||||
Let's take a polished aluminum plate (20 by 20 cm) at 125K (temperature of zero thermal expansion coefficient of silicon) surrounded by elements are 25 degrees (300 K):
|
||||
\\[ P = \epsilon \cdot \sigma \cdot A \cdot (T\_r^4 - T\_s^4) = 0.36\\, J \\]
|
||||
|
||||
</div>
|
||||
|
||||
|
||||
## Heat {#heat}
|
||||
|
||||
The _heat_ \\(Q\\) (in Joules) corresponds to the energy necessary to change the temperature of the mass with a certain material specific heat capacity:
|
||||
\\[ Q = m \cdot c \cdot \Delta T \\]
|
||||
with:
|
||||
|
||||
- \\(m\\) the mass in \\([kg]\\)
|
||||
- \\(c\\) the specific heat capacity in \\([J/kg \cdot K]\\)
|
||||
- \\(\Delta T\\) the temperature different \\([K]\\)
|
||||
|
||||
<div class="exampl">
|
||||
|
||||
Let's compute the heat (i.e. energy) necessary to increase a 1kg granite by 1 degree.
|
||||
The specific heat capacity of granite is \\(c = 790\\,[J/kg\cdot K]\\).
|
||||
The required heat is then:
|
||||
\\[ Q = m\cdot c \cdot \Delta T = 790 \\,J \\]
|
||||
|
||||
</div>
|
||||
|
||||
<a id="table--tab:specific-heat-capacity"></a>
|
||||
<div class="table-caption">
|
||||
<span class="table-number"><a href="#table--tab:specific-heat-capacity">Table 2</a>:</span>
|
||||
Some examples of specific heat capacity
|
||||
</div>
|
||||
|
||||
| Substance | Specific heat capacity [J/kg.K] |
|
||||
|---------------------|---------------------------------|
|
||||
| Air | 1012 |
|
||||
| Aluminium | 897 |
|
||||
| Copper | 385 |
|
||||
| Granite | 790 |
|
||||
| Steel | 466 |
|
||||
| Water at 25 degrees | 4182 |
|
||||
|
||||
|
||||
## Heat Transport (i.e. Water cooling) {#heat-transport--i-dot-e-dot-water-cooling}
|
||||
|
||||
<a id="figure--fig:heat-transfer-fluid"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/heat_transfer_fluid.png" caption="<span class=\"figure-number\">Figure 1: </span>Heat transfered to the fluid" >}}
|
||||
|
||||
\begin{equation}
|
||||
Q\_{in} = h \cdot A \cdot (T\_{wall} - T\_{mean})
|
||||
\end{equation}
|
||||
|
||||
<a id="figure--fig:heat-transport"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/heat_transport.png" caption="<span class=\"figure-number\">Figure 2: </span>Heat Transport in the fluid" >}}
|
||||
|
||||
\begin{equation}
|
||||
Q\_{out} = \phi \rho c\_p (T\_{mean,in} - T\_{mean,out})
|
||||
\end{equation}
|
||||
|
||||
with:
|
||||
|
||||
- \\(Q\_{out}\\) the transported heat in W
|
||||
- \\(\phi\\) the flow in \\(m^3/s\\)
|
||||
- \\(\rho\\) the fluid density in \\(kg/m^3\\)
|
||||
- \\(c\_p\\) the specific heat capacity of the fluid in \\(J/(kg \cdot K)\\)
|
||||
- \\(T\_{mean}\\) the mean incoming and outgoing fluid temperature
|
||||
|
||||
Because of energy balance, we have in the stationary condition: \\(Q\_{in} = Q\_{out}\\)
|
||||
|
||||
|
||||
## Heat flow {#heat-flow}
|
||||
|
||||
The heat flow \\(P\\) (in watt) is the derivative of the heat:
|
||||
\\[ P = \cdot{Q} = \frac{dQ}{dt} = \frac{dT}{R\_T} = C\_T \cdot dT \\]
|
||||
with:
|
||||
|
||||
- \\(Q\\) the heat in [W]
|
||||
- \\(R\_T\\) the thermal resistance in \\([K/W]\\)
|
||||
- \\(C\_T\\) the thermal conductance in \\([W/K]\\)
|
||||
|
||||
|
||||
## Bibliography {#bibliography}
|
||||
|
||||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||||
</div>
|
@@ -8,10 +8,10 @@ Tags
|
||||
:
|
||||
|
||||
|
||||
## Commercial products {#commercial-products}
|
||||
|
||||
<https://www.tcdirect.fr/product-2-300-1/Cartouche-chauffante>
|
||||
|
||||
<https://www.thorlabs.com/newgrouppage9.cfm?objectgroup_id=305>
|
||||
- <https://www.tcdirect.fr/product-2-300-1/Cartouche-chauffante>
|
||||
- <https://www.thorlabs.com/newgrouppage9.cfm?objectgroup_id=305>
|
||||
|
||||
|
||||
### Cryogenic temperature (~120K, UHV compatible) {#cryogenic-temperature--120k-uhv-compatible}
|
||||
|
@@ -9,6 +9,12 @@ Tags
|
||||
: [Slip Rings]({{< relref "slip_rings.md" >}})
|
||||
|
||||
|
||||
## Spindle Metrology {#spindle-metrology}
|
||||
|
||||
- <https://ingeniqs.be/engipedia/spindle-metrology/>
|
||||
- <https://ingeniqs.be/engipedia/error-separation-techniques/>
|
||||
|
||||
|
||||
## Manufacturers {#manufacturers}
|
||||
|
||||
| Manufacturers | Country |
|
||||
@@ -19,6 +25,59 @@ Tags
|
||||
| [PIC](https://www.airbearings.com/) | USA |
|
||||
|
||||
|
||||
## Spindle Errors {#spindle-errors}
|
||||
|
||||
|
||||
### Air Bearing Spindle (RT150U) {#air-bearing-spindle--rt150u}
|
||||
|
||||
<a id="figure--fig:rotation-stage-rt150u-xyz-errors"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/rotation_stage_rt150u_xyz_errors.png" caption="<span class=\"figure-number\">Figure 1: </span>RT150U - XYZ errors (eccentricity removed)" >}}
|
||||
|
||||
```matlab
|
||||
%% Frequency Analysis
|
||||
Nfft = floor(10.0/Ts);
|
||||
win = hanning(Nfft);
|
||||
Noverlap = floor(Nfft/2);
|
||||
|
||||
[pxx_ab_x, f_ab] = pwelch(x, win, Noverlap, Nfft, 1/Ts);
|
||||
[pxx_ab_y, ~ ] = pwelch(y, win, Noverlap, Nfft, 1/Ts);
|
||||
[pxx_ab_z, ~ ] = pwelch(z, win, Noverlap, Nfft, 1/Ts);
|
||||
```
|
||||
|
||||
<a id="figure--fig:rotation-stage-rt150u-asd-errors"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/rotation_stage_rt150u_asd_errors.png" caption="<span class=\"figure-number\">Figure 2: </span>Amplitude Spectral Density of the Spindle errors" >}}
|
||||
|
||||
<a id="figure--fig:rotation-stage-rt150u-cas-error"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/rotation_stage_rt150u_cas_error.png" caption="<span class=\"figure-number\">Figure 3: </span>Cumulative Amplitude Spectrum - RT150U" >}}
|
||||
|
||||
|
||||
### Mechanical Bearing Spindle (RT100BS) {#mechanical-bearing-spindle--rt100bs}
|
||||
|
||||
Rotation at 60deg/s (10RPM).
|
||||
|
||||
<a id="figure--fig:rotation-stage-rt100bs-asd-errors"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/rotation_stage_rt100bs_asd_errors.png" caption="<span class=\"figure-number\">Figure 4: </span>Amplitude Spectral Density of the Spindle errors" >}}
|
||||
|
||||
<a id="figure--fig:rotation-stage-rt100bs-cas-error"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/rotation_stage_rt100bs_cas_error.png" caption="<span class=\"figure-number\">Figure 5: </span>Cumulative Amplitude Spectrum - RT100BS" >}}
|
||||
|
||||
|
||||
### Ultrasonic Piezo Motor Spindle (XRT-U-30) {#ultrasonic-piezo-motor-spindle--xrt-u-30}
|
||||
|
||||
<https://xeryon.com/products/precision-rotation-stages/xrt-u-30/>
|
||||
|
||||
Rotation at few degrees per second.
|
||||
|
||||
```matlab
|
||||
length(axial_turns)/5 % Number of point per turn
|
||||
```
|
||||
|
||||
|
||||
## Bibliography {#bibliography}
|
||||
|
||||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||||
|
@@ -19,9 +19,11 @@ An active temperature control system generally consists of:
|
||||
- <https://www.thinksrs.com/products/tempcontrol.html>
|
||||
- <https://www.thorlabs.com/thorproduct.cfm?partnumber=TC300B>
|
||||
- <https://www.vescent.com/manuals/doku.php?id=slice:qt>
|
||||
- CN16DPT-305-EIP-DC, OMEGA Inc. (<a href="#citeproc_bib_item_1">Holler et al. 2022</a>)
|
||||
|
||||
|
||||
## Bibliography {#bibliography}
|
||||
|
||||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||||
<div class="csl-entry"><a id="citeproc_bib_item_1"></a>Holler, Mirko, Tomas Aidukas, Lars Heller, Christian Appel, Nicholas W. Phillips, Elisabeth Müller-Gubler, Manuel Guizar-Sicairos, Jörg Raabe, and Johannes Ihli. 2022. “Environmental Control for X-Ray Nanotomography.” <i>Journal of Synchrotron Radiation</i> 29 (5): 1223–31. doi:<a href="https://doi.org/10.1107/s1600577522006968">10.1107/s1600577522006968</a>.</div>
|
||||
</div>
|
||||
|
@@ -10,29 +10,13 @@ Tags
|
||||
|
||||
## Temperature sensors types {#temperature-sensors-types}
|
||||
|
||||
There are three main types of temperature sensors:
|
||||
There are mainly three main types of temperature sensors:
|
||||
|
||||
- [Thermocouples](#org-target--sec-temperature-sensor-thermocouple) with are based on the Seebeck effect
|
||||
- [RTD](#org-target--sec-temperature-sensor-rtd) (Resistance Temperature Detectors): made of pure metals (Pt, Ni or Cu)
|
||||
They are all PTC (Positive Temperature Coefficient): PT100, PT1000, Ni100, Ni1000, ...
|
||||
- [Thermistor](#org-target--sec-temperature-sensor-thermistor): made of metal oxide mixtures (semiconductor materials).
|
||||
It can have a NTC (Negative Temperature Coefficient) or a PTC (Positive Temperature Coefficient).
|
||||
- [Thermocouple](#org-target--sec-temperature-sensor-thermocouple) (Seebeck effect):
|
||||
- Type K (-180 to 1200 degC)
|
||||
- Type T (-250 to 350 degC)
|
||||
|
||||
|
||||
### RTD {#rtd}
|
||||
|
||||
<span class="org-target" id="org-target--sec-temperature-sensor-rtd"></span>
|
||||
|
||||
Sensitivity of PT100 is typically around 0.004 Ohm/Ohm/deg.
|
||||
|
||||
|
||||
### Thermistor {#thermistor}
|
||||
|
||||
<span class="org-target" id="org-target--sec-temperature-sensor-thermistor"></span>
|
||||
|
||||
Sensitivity of NTC is typically around -0.05 Ohm/Ohm/deg.
|
||||
They typically have a NTC (Negative Temperature Coefficient).
|
||||
|
||||
|
||||
### Thermocouple {#thermocouple}
|
||||
@@ -40,17 +24,134 @@ Sensitivity of NTC is typically around -0.05 Ohm/Ohm/deg.
|
||||
<span class="org-target" id="org-target--sec-temperature-sensor-thermocouple"></span>
|
||||
|
||||
|
||||
#### Working Principle {#working-principle}
|
||||
|
||||
Based on Seebeck effect: in a conductor, a temperature difference \\(\Delta T\\) creates an electric field \\(V\\) defined by the Seebeck coefficient \\(\epsilon\\).
|
||||
|
||||
Consider a material with a Seebeck coefficient \\(\epsilon\_A\\), the voltage difference at its send is:
|
||||
|
||||
\begin{equation}
|
||||
V\_A = \int\_{T\_1}^{T\_2} \epsilon\_A dT = \epsilon\_A (T\_2 - T\_1)
|
||||
\end{equation}
|
||||
|
||||
<a id="figure--fig:temperature-seebeck-material-A"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_seebeck_material_A.png" caption="<span class=\"figure-number\">Figure 1: </span>Seebeck effect in a material" >}}
|
||||
|
||||
Now if two materials are "chained", the overall voltage will be:
|
||||
|
||||
\begin{equation}
|
||||
V = (\epsilon\_{A} - \epsilon\_{B})(T\_{2} - T\_{1})
|
||||
\end{equation}
|
||||
|
||||
<a id="figure--fig:temperature-thermocouple-two-materiasl"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_thermocouple_two_materiasl.png" caption="<span class=\"figure-number\">Figure 2: </span>Combination of two materials" >}}
|
||||
|
||||
Several materials can be combined:
|
||||
|
||||
<a id="figure--fig:temperature-termocouples-sensitivities"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_termocouples_sensitivities.png" caption="<span class=\"figure-number\">Figure 3: </span>Some types of thermocouples and associated sensitivities" >}}
|
||||
|
||||
|
||||
#### Measurement of the temperature {#measurement-of-the-temperature}
|
||||
|
||||
The output voltage of a thermocouple is given by \\(\Delta T\\) (between \\(T\_1\\), the temperature at the measured voltage side, and \\(T\_2\\) at the thermocouple side).
|
||||
Therefore, typically \\(T\_1\\) is measured with an RTD (typically what is done inside the Agilent 34970A) and \\(T\_2\\) is therefore estimated from \\(T\_1\\) and \\(\Delta T\\).
|
||||
|
||||
<a id="figure--fig:temperature-thermocouple-meas"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_thermocouple_meas.png" caption="<span class=\"figure-number\">Figure 4: </span>Thermocouple measured voltage" >}}
|
||||
|
||||
|
||||
#### Choice of Thermocouple {#choice-of-thermocouple}
|
||||
|
||||
- J - Iron / Constantan
|
||||
- \\(55\ \mu V / K\\)
|
||||
- Easy to solder
|
||||
- Can act as galvanic element
|
||||
- K - Chromel / Alumel
|
||||
- \\(40\ \mu V / K\\), almost linear
|
||||
- Difficult to solver, welding is better
|
||||
- Low thermal conductivity of wires
|
||||
- T - Copper / Constantan
|
||||
- \\(40\ \mu V / K\\)
|
||||
- Easy to solder
|
||||
- High thermal conductivity is source of errors
|
||||
|
||||
|
||||
#### Summary {#summary}
|
||||
|
||||
Advantages:
|
||||
|
||||
- Simple to use
|
||||
- Standard acquisition systems
|
||||
- Small (thin wires down to 0.08mm)
|
||||
- Fast response
|
||||
- Suitable for high and low temperatures
|
||||
- Can be used in vacuum systems
|
||||
- No self heating
|
||||
|
||||
Disadvantages:
|
||||
|
||||
- Relative expensive
|
||||
- Uncertainty about 0.1K, not feasible to reach mK
|
||||
|
||||
|
||||
### RTD {#rtd}
|
||||
|
||||
<span class="org-target" id="org-target--sec-temperature-sensor-rtd"></span>
|
||||
|
||||
Sensitivity of PT100 is typically around .
|
||||
|
||||
Advantages:
|
||||
|
||||
- Very high stability (better than 1mK per year)
|
||||
- Very high linearity (0.1% over -40 to 125 degrees)
|
||||
- This makes them very useful as calibration reference sensor, which is linked to the international standard
|
||||
|
||||
Disadvantages:
|
||||
|
||||
- Expensive
|
||||
- Low sensitivity (typically 0.004 Ohm/Ohm/deg)
|
||||
- The measurement is sensitivity to lead wire resistance, but four-wire technique may be used
|
||||
- Self heating due to electrical dissipation.
|
||||
Typically, for a Pt100, \\(P = 0.1 mW\\) (the source voltage is typically 0.1 V)
|
||||
This corresponds to approximately 0.1 degree of self heating in "still" air
|
||||
|
||||
|
||||
### Thermistor {#thermistor}
|
||||
|
||||
<span class="org-target" id="org-target--sec-temperature-sensor-thermistor"></span>
|
||||
|
||||
Sensitivity of NTC is .
|
||||
|
||||
Advantages:
|
||||
|
||||
- Highest sensitivity (typically around -0.05 Ohm/Ohm/deg)
|
||||
- Because the resistance is typically high (100k Ohm), no Four wire configuration is necessary, and long wires may be used
|
||||
- Lower heat dissipation than Pt100
|
||||
- Very high stability, especially for glass encapsulated
|
||||
- Very small, and available in all kinds of shapes
|
||||
|
||||
Disadvantages:
|
||||
|
||||
- Non-linear, so compensation is necessary (but not really an issue with software compensation)
|
||||
- Self-heating if mK accuracy is wanted
|
||||
- Cover is necessary for use in vacuum
|
||||
|
||||
|
||||
### Comparison of sensor types {#comparison-of-sensor-types}
|
||||
|
||||
<https://www.ni.com/fr/shop/data-acquisition/sensor-fundamentals/measuring-temperature-with-thermocouples-rtds-and-thermistors.html?srsltid=AfmBOoqrR2VCWVvkCTjzTIYHlyiKVpt6Ket1xfhU1yFPTtXsGKg0RILE>
|
||||
|
||||
| | RTD | Thermistor | Thermocouple |
|
||||
|---------------|------|------------|--------------|
|
||||
| Accuracy | Good | Non-Linear | |
|
||||
| Stability | | | |
|
||||
| Sensitivity | Bad | Good | |
|
||||
| Response time | | | |
|
||||
| Self heating | | | |
|
||||
| | RTD | Thermistor | Thermocouple | |
|
||||
|---------------|-------------------------|------------|--------------|---|
|
||||
| Accuracy | Good | Non-Linear | | A |
|
||||
| Linearity | 0.1% over -40..125 degC | | | |
|
||||
| Stability | better than 1mK/year | | | |
|
||||
| Sensitivity | 0.4%/K | 5%/K | 50uV/K | |
|
||||
| Response time | | | | |
|
||||
| Self heating | | | None | |
|
||||
|
||||
|
||||
## Accuracy of Temperature measurement {#accuracy-of-temperature-measurement}
|
||||
@@ -73,7 +174,7 @@ For a Pt100 at \\(0^oC\\), this corresponds to an accuracy of \\(< \pm 0.04\\,K\
|
||||
#### 2, 3 and 4 wires sensors {#2-3-and-4-wires-sensors}
|
||||
|
||||
The measured resistance is the sum of the resistance of the sensitive element and the resistance of the wires.
|
||||
This corresponds to the 2-wire measurement ([Figure 1](#figure--fig:temperature-sensor-rtd-2-wires)).
|
||||
This corresponds to the 2-wire measurement ([Figure 5](#figure--fig:temperature-sensor-rtd-2-wires)).
|
||||
|
||||
The errors associated with this effect are large when the resistance of the sensitive element is small and then the resistance of all cables and connectors are large.
|
||||
For instance, the effect of contact/wire resistance less important for the PT1000 than for the PT100.
|
||||
@@ -81,18 +182,27 @@ The use of 2 wire PT1000 is possible (whereas for PT100, 4 wire is more accurate
|
||||
|
||||
<a id="figure--fig:temperature-sensor-rtd-2-wires"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_2_wires.png" caption="<span class=\"figure-number\">Figure 1: </span>2-wire measurement" >}}
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_2_wires.png" caption="<span class=\"figure-number\">Figure 5: </span>2-wire measurement" >}}
|
||||
|
||||
The effect of the resistance of the wires (cables, connectors, etc..) can be mitigated by using the 4-wire configuration ([Figure 2](#figure--fig:temperature-sensor-rtd-4-wires)).
|
||||
The effect of the resistance of the wires (cables, connectors, etc..) can be mitigated by using the 4-wire configuration ([Figure 6](#figure--fig:temperature-sensor-rtd-4-wires)).
|
||||
|
||||
<a id="figure--fig:temperature-sensor-rtd-4-wires"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_4_wires.png" caption="<span class=\"figure-number\">Figure 2: </span>4-wire measurement" >}}
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_4_wires.png" caption="<span class=\"figure-number\">Figure 6: </span>4-wire measurement" >}}
|
||||
|
||||
|
||||
### Temperature {#temperature}
|
||||
|
||||
|
||||
#### Effect of conductivity through the wires {#effect-of-conductivity-through-the-wires}
|
||||
|
||||
It is better to use thin wires, of the fix the wires to the part that is to be measured.
|
||||
|
||||
<a id="figure--fig:temperature-effect-wires"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_effect_wires.png" caption="<span class=\"figure-number\">Figure 7: </span>Measured effect of wires. When in \"air\", it conducts the heat from the air to the sensor which can lead to measurement errors." >}}
|
||||
|
||||
|
||||
#### Thermal contact and response time {#thermal-contact-and-response-time}
|
||||
|
||||
The measured temperature is the temperature of the sensitive element.
|
||||
@@ -100,6 +210,10 @@ It may not be equal to the temperature of the element on which the sensor is fix
|
||||
|
||||
It depends on the thermal contact and the response time in play.
|
||||
|
||||
The sensor contact may be improved by using "soft" (i.e. plastically deformable) metals at the contact interface such as indium.
|
||||
|
||||
However, it seems that having too much pressure in the sensor may induce stress in the NTC that can induce measurement errors.
|
||||
|
||||
|
||||
#### Self heating effect {#self-heating-effect}
|
||||
|
||||
@@ -113,12 +227,14 @@ Typically, a constant voltage is applied, such that the generated current is low
|
||||
|
||||
<div class="exampl">
|
||||
|
||||
The applied voltage is typically in the order of 1V.
|
||||
For a Pt100 (\\(R \approx 100\\,\Omega\\)), this would lead a heat generation of \\(P \approx 10\\,mW\\).
|
||||
For a NTC with \\(R\approx 10\\,k\Omega\\), the heat generation will me much lower \\(P\approx 100\\,\mu W\\).
|
||||
The applied voltage is typically in the order of 0.1V to 1V.
|
||||
For a Pt100 (\\(R \approx 100\\,\Omega\\)), this would lead a heat generation of \\(P \approx 1 \text{ to } 10\\,mW\\).
|
||||
For a NTC with \\(R\approx 10\\,k\Omega\\), the heat generation will me much lower \\(P\approx 10 \text{ to } 100\\,\mu W\\).
|
||||
|
||||
</div>
|
||||
|
||||
In order to lower the self heating effect, _intermitted_ currents may be used as is the case with the Agilent 34970A.
|
||||
|
||||
|
||||
### Converting Resistance to Temperature {#converting-resistance-to-temperature}
|
||||
|
||||
@@ -150,26 +266,26 @@ Manufacturers usually provides a lookup table that links the resistance and the
|
||||
|
||||
### PT100 {#pt100}
|
||||
|
||||
A PT100 resistance is quite linear with respect to the temperature as shown in [Figure 3](#figure--fig:temperature-sensor-pt100-resistance).
|
||||
A PT100 resistance is quite linear with respect to the temperature as shown in [Figure 8](#figure--fig:temperature-sensor-pt100-resistance).
|
||||
|
||||
<a id="figure--fig:temperature-sensor-pt100-resistance"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_pt100_resistance.png" caption="<span class=\"figure-number\">Figure 3: </span>Resistance of a PT100 as a function of the temperature" >}}
|
||||
{{< figure src="/ox-hugo/temperature_sensor_pt100_resistance.png" caption="<span class=\"figure-number\">Figure 8: </span>Resistance of a PT100 as a function of the temperature" >}}
|
||||
|
||||
The coefficient of resistance \\(\alpha\\) is defined as the ratio of the rate of change of resistance with temperature to the resistance of the thermistor at a specified temperature:
|
||||
\\[ \alpha(T) = \frac{1}{R(T)}\frac{dR(T)}{dT} \\]
|
||||
|
||||
For a PT100, it is displayed in [Figure 4](#figure--fig:temperature-sensor-pt100-sensitivity).
|
||||
For a PT100, it is displayed in [Figure 9](#figure--fig:temperature-sensor-pt100-sensitivity).
|
||||
At \\(0^oC\\), \\(\alpha(0^oC) \approx 0.004\\,\Omega/\Omega/{}^oC\\).
|
||||
|
||||
<a id="figure--fig:temperature-sensor-pt100-sensitivity"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_pt100_sensitivity.png" caption="<span class=\"figure-number\">Figure 4: </span>Sensitivity of a PT100 as a function of the temperature" >}}
|
||||
{{< figure src="/ox-hugo/temperature_sensor_pt100_sensitivity.png" caption="<span class=\"figure-number\">Figure 9: </span>Sensitivity of a PT100 as a function of the temperature" >}}
|
||||
|
||||
|
||||
### NTC {#ntc}
|
||||
|
||||
A NTC is much more non-linear than a PT100 as shown in [Figure 5](#figure--fig:temperature-sensor-rtd-resistance).
|
||||
A NTC is much more non-linear than a PT100 as shown in [Figure 10](#figure--fig:temperature-sensor-rtd-resistance).
|
||||
|
||||
The NTC used here is "Type F" from Amphenol Thermometrics.
|
||||
|
||||
@@ -180,18 +296,241 @@ R_rtd = 1e4*[68.60 48.16 34.23 24.62 17.91 13.17 9.782 7.339 5.558 4.247 3.274 2
|
||||
|
||||
<a id="figure--fig:temperature-sensor-rtd-resistance"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_resistance.png" caption="<span class=\"figure-number\">Figure 5: </span>Resistance of a RTD as a function of the temperature" >}}
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_resistance.png" caption="<span class=\"figure-number\">Figure 10: </span>Resistance of a RTD as a function of the temperature" >}}
|
||||
|
||||
The huge advantage of RTD compared to PT100 is that the sensitivity is much larger than Pt100 as shown in [Figure 6](#figure--fig:temperature-sensor-rtd-sensitivity).
|
||||
The huge advantage of RTD compared to PT100 is that the sensitivity is much larger than Pt100 as shown in [Figure 11](#figure--fig:temperature-sensor-rtd-sensitivity).
|
||||
|
||||
<a id="figure--fig:temperature-sensor-rtd-sensitivity"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_sensitivity.png" caption="<span class=\"figure-number\">Figure 6: </span>Sensitivity of a RTD as a function of the temperature" >}}
|
||||
{{< figure src="/ox-hugo/temperature_sensor_rtd_sensitivity.png" caption="<span class=\"figure-number\">Figure 11: </span>Sensitivity of a RTD as a function of the temperature" >}}
|
||||
|
||||
|
||||
## Compute temperature from the measured resistance {#compute-temperature-from-the-measured-resistance}
|
||||
|
||||
|
||||
### Pt100 and Pt1000 {#pt100-and-pt1000}
|
||||
|
||||
The resistance as a function of temperature is approximated by the Callendar–Van Dusen equation:
|
||||
\\[ R(T) = \begin{cases}
|
||||
R\_0 (1 + A \cdot T + B \cdot T^2), & \text{for } T>0^oC \\\\
|
||||
R\_0 (1 + A\cdot T + B \cdot T^2 + C \cdot (T - 100) \cdot T^3), & \text{for } T<0^oC
|
||||
\end{cases} \\]
|
||||
with \\(R\_0\\) the resistance value at 0 degrees (\\(100\\,\Omega\\) for a Pt100 and \\(1000\\,\Omega\\) for a Pt1000).
|
||||
|
||||
Values for A, B, C and D are depending on the exact model (summarized in [Table 1](#table--tab:pt100-values)).
|
||||
|
||||
<a id="table--tab:pt100-values"></a>
|
||||
<div class="table-caption">
|
||||
<span class="table-number"><a href="#table--tab:pt100-values">Table 1</a>:</span>
|
||||
Values of the Callendar-Van Dusen equations
|
||||
</div>
|
||||
|
||||
| TCR | A | B | C |
|
||||
|------------|----------------------------|------------------------------|------------------------------|
|
||||
| 3850 ppm/K | \\(3.9083 \cdot 10^{-3}\\) | \\(-5.775 \cdot 10^{-7}\\) | \\(-4.183 \cdot 10^{-12}\\) |
|
||||
| 3911 ppm/K | \\(3.9692 \cdot 10^{-3}\\) | \\(-5.829 \cdot 10^{-7}\\) | \\(-4.3303 \cdot 10^{-12}\\) |
|
||||
| 3750 ppm/K | \\(3.8102 \cdot 10^{-3}\\) | \\(-6.01888 \cdot 10^{-7}\\) | \\(-6 \cdot 10^{-12}\\) |
|
||||
| 3770 ppm/K | \\(3.8285 \cdot 10^{-3}\\) | \\(-5.85 \cdot 10^{-7}\\) | |
|
||||
|
||||
```matlab
|
||||
%% Pt100 (3850 ppm/K)
|
||||
R0 = 100; % [Ohm]
|
||||
|
||||
A = 3.9083e-3; % [degC^-1]
|
||||
B = -5.775e-7; % [degC^-2]
|
||||
C = -4.183e-12; % [degC^-4]
|
||||
|
||||
T1 = -200:0; % [degC]
|
||||
T2 = 0:850; % [degC]
|
||||
T = [T1,T2]; % [degC]
|
||||
|
||||
R = [R0*(1 + A*T1 + B*T1.^2 + C*(T1-100).*T1.^3), R0*(1 + A*T2 + B*T2.^2)]; % [Ohm]
|
||||
```
|
||||
|
||||
<a id="figure--fig:temperature-sensor-pt100-curve"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_pt100_curve.png" caption="<span class=\"figure-number\">Figure 12: </span>Resistance as a function of the temperature for a Pt100" >}}
|
||||
|
||||
For temperatures above 0 degrees, the temperature \\(T\\) can be easily computed from the measured resistance \\(R\\) using:
|
||||
\\[ T = \frac{-A + \sqrt{A^2 - 4 B ( 1 - R/R\_0 )}}{2 B} \\]
|
||||
|
||||
For temperatures below 0 degrees, the equation is harder to solve analytically, and a lookup table is more appropriate.
|
||||
|
||||
Let's compare the temperature given by a Loopup table and the temperature given by the analytical formula in two cases:
|
||||
|
||||
- linear interpolation with one point every degree
|
||||
- cubic interpolation with one point every 10 degrees
|
||||
|
||||
The error is less than 0.1mK over the full range, validating the use of a lookup table to convert the resistance to temperature ([Figure 14](#figure--fig:temperature-sensor-lut-errors)).
|
||||
|
||||
|
||||
### NTC thermistor {#ntc-thermistor}
|
||||
|
||||
The resistance of the NTC thermistor as a function of the temperature can be well approximated with the following equation:
|
||||
\\[ R\_t = R\_{25} \cdot e^{A + B/T + C/T^2 + D/T^3 \\]
|
||||
where \\(T\\) is the temperature in kelvins, \\(R\_{25}\\) the nominal resistance at \\(25^oC\\), \\(A\\), \\(B\\), \\(C\\) and \\(D\\) are coefficients which are specific for a given thermistor.
|
||||
|
||||
Typically, coefficients A, B, C and D are varying with temperature as shown in [Table 2](#table--tab:temperature-sensor-ntc-coefs).
|
||||
|
||||
<a id="table--tab:temperature-sensor-ntc-coefs"></a>
|
||||
<div class="table-caption">
|
||||
<span class="table-number"><a href="#table--tab:temperature-sensor-ntc-coefs">Table 2</a>:</span>
|
||||
Example of A, B, C and D coeficients for an NTC thermistor (DC95F202VN)
|
||||
</div>
|
||||
|
||||
| | A | B | C | D |
|
||||
|------------|----------------|---------------|----------------|----------------|
|
||||
| -50 to 0 | -1.4122478E+01 | 4.4136033E+03 | -2.9034189E+04 | -9.3875035E+06 |
|
||||
| 0 to 50 | -1.4141963E+01 | 4.4307830E+03 | -3.4078983E+04 | -8.8941929E+06 |
|
||||
| 50 to 100 | -1.4202172E+01 | 4.4975256E+03 | -5.8421357E+04 | -5.9658796E+06 |
|
||||
| 100 to 150 | -1.6154078E+01 | 6.8483992E+03 | -1.0004049E+06 | 1.1961431E+08 |
|
||||
|
||||
```matlab
|
||||
%% Compute the resistance as a function of the temperature for a given NTC (DC95F202VN)
|
||||
R0 = 2e3; % Resistance at 25deg
|
||||
|
||||
T1 = 273.15+[-50:0]; % [degK]
|
||||
T2 = 273.15+[1:50]; % [degK]
|
||||
T3 = 273.15+[51:100]; % [degK]
|
||||
T4 = 273.15+[101:150]; % [degK]
|
||||
|
||||
R = R0*exp([[-1.4122478E+01 + 4.4136033E+03./T1 - 2.9034189E+04./T1.^2 - 9.3875035E+06./T1.^3]';
|
||||
[-1.4141963E+01 + 4.4307830E+03./T2 - 3.4078983E+04./T2.^2 - 8.8941929E+06./T2.^3]';
|
||||
[-1.4202172E+01 + 4.4975256E+03./T3 - 5.8421357E+04./T3.^2 - 5.9658796E+06./T3.^3]';
|
||||
[-1.6154078E+01 + 6.8483992E+03./T4 - 1.0004049E+06./T4.^2 + 1.1961431E+08./T4.^3]'])'; % [Ohm]
|
||||
|
||||
T = -273.15+[T1,T2,T3,T4]; % [degC]
|
||||
```
|
||||
|
||||
<a id="figure--fig:temperature-sensor-ntc-curve"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_ntc_curve.png" caption="<span class=\"figure-number\">Figure 13: </span>Resistance as a function of the temperature for a given NTC" >}}
|
||||
|
||||
To calculate the actual thermistor temperature as a function of the measured thermistor resistance, use the following equation:
|
||||
\\[ T = \frac{1}{a + b \ln(R\_t/R\_{25}) + c (Ln Rt/R25)^2 + d (Ln Rt/R25)^3) \\]
|
||||
|
||||
<a id="table--tab:temperature-sensor-ntc-equation"></a>
|
||||
<div class="table-caption">
|
||||
<span class="table-number"><a href="#table--tab:temperature-sensor-ntc-equation">Table 3</a>:</span>
|
||||
Coefficients used to compute the temperature as a function of the resistance
|
||||
</div>
|
||||
|
||||
| Rt/R25 range | a | b | c | d |
|
||||
|--------------------|---------------|---------------|----------------|----------------|
|
||||
| 68.600 to 3.274 | 3.3538646E-03 | 2.5654090E-04 | 1.9243889E-06 | 1.0969244E-07 |
|
||||
| 3.274 to 0.36036 | 3.3540154E-03 | 2.5627725E-04 | 2.0829210E-06 | 7.3003206E-08 |
|
||||
| 0.36036 to 0.06831 | 3.3539264E-03 | 2.5609446E-04 | 1.9621987E-06 | 4.6045930E-08 |
|
||||
| 0.06831 to 0.01872 | 3.3368620E-03 | 2.4057263E-04 | -2.6687093E-06 | -4.0719355E-07 |
|
||||
|
||||
|
||||
### Approximation of formulas using lookup tables {#approximation-of-formulas-using-lookup-tables}
|
||||
|
||||
First, let's compare the analytical formula with a LUT for a Pt100 ([Figure 14](#figure--fig:temperature-sensor-lut-errors)).
|
||||
The error (accuracy) is bellow 0.1mK for relatively small LUT.
|
||||
|
||||
```matlab
|
||||
%% "Perfect" temperature and resistance
|
||||
R0 = 100; % [Ohm]
|
||||
A = 3.9083e-3; % [degC^-1]
|
||||
B = -5.775e-7; % [degC^-2]
|
||||
C = -4.183e-12; % [degC^-4]
|
||||
|
||||
T1 = -200:0.1:0; % [degC]
|
||||
T2 = 0.1:0.1:850; % [degC]
|
||||
T_true = [T1,T2]; % [degC]
|
||||
R_true = [R0*(1 + A*T1 + B*T1.^2 + C*(T1-100).*T1.^3), R0*(1 + A*T2 + B*T2.^2)]; % [Ohm]
|
||||
|
||||
%% Lookup table for Pt100 (3850 ppm/K) - Linear
|
||||
dT = 1;
|
||||
interp_method = 'linear';
|
||||
|
||||
T1 = -200:dT:0; % [degC]
|
||||
T2 = dT:dT:850; % [degC]
|
||||
T_lut_linear = [T1,T2]; % [degC]
|
||||
R_lut_linear = [R0*(1 + A*T1 + B*T1.^2 + C*(T1-100).*T1.^3), R0*(1 + A*T2 + B*T2.^2)]; % [Ohm]
|
||||
|
||||
T_meas_linear = interp1(R_lut_linear,T_lut_linear,R_true,interp_method); % interpolate the resistance using the LUT to find the corresponding temperature
|
||||
|
||||
%% Lookup table for Pt100 (3850 ppm/K) - Makima
|
||||
dT = 10;
|
||||
interp_method = 'makima';
|
||||
|
||||
T1 = -200:dT:0; % [degC]
|
||||
T2 = dT:dT:850; % [degC]
|
||||
T_lut_makima = [T1,T2]; % [degC]
|
||||
R_lut_makima = [R0*(1 + A*T1 + B*T1.^2 + C*(T1-100).*T1.^3), R0*(1 + A*T2 + B*T2.^2)]; % [Ohm]
|
||||
|
||||
T_meas_makima = interp1(R_lut_makima,T_lut_makima,R_true,interp_method); % interpolate the resistance using the LUT to find the corresponding temperature
|
||||
```
|
||||
|
||||
<a id="figure--fig:temperature-sensor-lut-errors"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_lut_errors.png" caption="<span class=\"figure-number\">Figure 14: </span>Interpolation errors in two cases when using a LUT for a Pt100" >}}
|
||||
|
||||
NTC thermistors are more non-linear and therefore require finer LUT to have low accuracy errors.
|
||||
In order to have less than 0.1mK of accuracy, a LUT with linear interpolation requires approximately one point every 0.1 degree ([Figure 15](#figure--fig:temperature-sensor-lut-errors-ntc)).
|
||||
|
||||
```matlab
|
||||
%% "Perfect" temperature and resistance of NTC (DC95F202VN)
|
||||
R0 = 2e3; % Resistance at 25deg
|
||||
dT_true = 0.01;
|
||||
|
||||
T1 = 273.15+[-50:dT_true:0]; % [degK]
|
||||
T2 = 273.15+[0+dT_true:dT_true:50]; % [degK]
|
||||
T3 = 273.15+[50+dT_true:dT_true:100]; % [degK]
|
||||
T4 = 273.15+[100+dT_true:dT_true:150]; % [degK]
|
||||
|
||||
R_true = R0*exp([[-1.4122478E+01 + 4.4136033E+03./T1 - 2.9034189E+04./T1.^2 - 9.3875035E+06./T1.^3]';
|
||||
[-1.4141963E+01 + 4.4307830E+03./T2 - 3.4078983E+04./T2.^2 - 8.8941929E+06./T2.^3]';
|
||||
[-1.4202172E+01 + 4.4975256E+03./T3 - 5.8421357E+04./T3.^2 - 5.9658796E+06./T3.^3]';
|
||||
[-1.6154078E+01 + 6.8483992E+03./T4 - 1.0004049E+06./T4.^2 + 1.1961431E+08./T4.^3]'])'; % [Ohm]
|
||||
|
||||
T_true = -273.15+[T1,T2,T3,T4]; % [degC]
|
||||
|
||||
%% Lookup table for NTC (DC95F202VN) - Linear
|
||||
dT = 0.1;
|
||||
interp_method = 'linear';
|
||||
|
||||
T1 = 273.15+[-50:dT:0]; % [degK]
|
||||
T2 = 273.15+[0+dT:dT:50]; % [degK]
|
||||
T3 = 273.15+[50+dT:dT:100]; % [degK]
|
||||
T4 = 273.15+[100+dT:dT:150]; % [degK]
|
||||
T_lut_linear = -273.15+[T1,T2,T3,T4]; % [degC]
|
||||
R_lut_linear = R0*exp([[-1.4122478E+01 + 4.4136033E+03./T1 - 2.9034189E+04./T1.^2 - 9.3875035E+06./T1.^3]';
|
||||
[-1.4141963E+01 + 4.4307830E+03./T2 - 3.4078983E+04./T2.^2 - 8.8941929E+06./T2.^3]';
|
||||
[-1.4202172E+01 + 4.4975256E+03./T3 - 5.8421357E+04./T3.^2 - 5.9658796E+06./T3.^3]';
|
||||
[-1.6154078E+01 + 6.8483992E+03./T4 - 1.0004049E+06./T4.^2 + 1.1961431E+08./T4.^3]'])'; % [Ohm]
|
||||
|
||||
T_meas_linear = interp1(R_lut_linear,T_lut_linear,R_true,interp_method); % interpolate the resistance using the LUT to find the corresponding temperature
|
||||
|
||||
%% Lookup table for Pt100 (3850 ppm/K) - Makima
|
||||
dT = 1;
|
||||
interp_method = 'makima';
|
||||
|
||||
T1 = 273.15+[-50:dT:0]; % [degK]
|
||||
T2 = 273.15+[0+dT:dT:50]; % [degK]
|
||||
T3 = 273.15+[50+dT:dT:100]; % [degK]
|
||||
T4 = 273.15+[100+dT:dT:150]; % [degK]
|
||||
T_lut_makima = -273.15+[T1,T2,T3,T4]; % [degC]
|
||||
R_lut_makima = R0*exp([[-1.4122478E+01 + 4.4136033E+03./T1 - 2.9034189E+04./T1.^2 - 9.3875035E+06./T1.^3]';
|
||||
[-1.4141963E+01 + 4.4307830E+03./T2 - 3.4078983E+04./T2.^2 - 8.8941929E+06./T2.^3]';
|
||||
[-1.4202172E+01 + 4.4975256E+03./T3 - 5.8421357E+04./T3.^2 - 5.9658796E+06./T3.^3]';
|
||||
[-1.6154078E+01 + 6.8483992E+03./T4 - 1.0004049E+06./T4.^2 + 1.1961431E+08./T4.^3]'])'; % [Ohm]
|
||||
|
||||
T_meas_makima = interp1(R_lut_makima,T_lut_makima,R_true,interp_method); % interpolate the resistance using the LUT to find the corresponding temperature
|
||||
```
|
||||
|
||||
<a id="figure--fig:temperature-sensor-lut-errors-ntc"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_sensor_lut_errors_ntc.png" caption="<span class=\"figure-number\">Figure 15: </span>Interpolation errors in two cases when using a LUT for a NTC" >}}
|
||||
|
||||
|
||||
## Commercial Temperature Sensors {#commercial-temperature-sensors}
|
||||
|
||||
|
||||
### 20 degC {#20-degc}
|
||||
|
||||
|
||||
### 20 degC, Vacuum compatible {#20-degc-vacuum-compatible}
|
||||
|
||||
From (<a href="#citeproc_bib_item_2">Neto et al. 2022</a>), UHV compatible:
|
||||
@@ -200,7 +539,7 @@ From (<a href="#citeproc_bib_item_2">Neto et al. 2022</a>), UHV compatible:
|
||||
> The part-though-hole (PTH) sensors were soldered to thin, 30 AWG, varnish insulated copper wires with small amounts of tin-lead (70/30) alloy.
|
||||
|
||||
|
||||
### Cryogenic temperatures {#cryogenic-temperatures}
|
||||
### Cryogenic temperatures (77K / -200degC) {#cryogenic-temperatures--77k-200degc}
|
||||
|
||||
- <https://shop.allectra.com/products/343-pt100-c2>
|
||||
- <https://www.lakeshore.com/products/categories/overview/temperature-products/cryogenic-temperature-sensors/cernox>
|
||||
@@ -215,10 +554,28 @@ From (<a href="#citeproc_bib_item_2">Neto et al. 2022</a>)
|
||||
> Furthermore, the thin platinum wire of the 10 kΩ RTDs presented bad solderability and its assembly process was too laborious, resulting in unreliable mechanical bonds and a failure rate beyond acceptable for a robust beamline instrumentation.
|
||||
> The alternative was to use **2 kΩ IST RTDs (P2K0.232.3FW.B.007)** with custom-made flat gold-platted terminals, resulting in a full range sensor with better solderability and temperature **resolution below 0.4 mK** over the entire measurable range
|
||||
|
||||
**NTC, Amphenol DC95F**, measured with an Agilent 34970A at 10K:
|
||||
|
||||
- leads to 0.2mK resolution (22 bits)
|
||||
- High interchangeability: offset of 0.01K and sensitivity of 3mK/K
|
||||
|
||||
<a id="figure--fig:temperature-ntc-dc95f-results"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_ntc_dc95f_results.png" caption="<span class=\"figure-number\">Figure 16: </span>6 (un-calibrated) DC95F sensors fixed to the same mass with homogeneous temperature" >}}
|
||||
|
||||
**NTC, Betatherm 10K3A1**, measured at 10K:
|
||||
|
||||
- Resolution of 0.2mK
|
||||
- Low noise and high repeatability
|
||||
|
||||
<a id="figure--fig:temperature-betatherm-10K3A1-results"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/temperature_betatherm_10K3A1_results.png" caption="<span class=\"figure-number\">Figure 17: </span>Measured temperature of two BetaTherm 10K3A1 compared to a reference sensor, at 10K" >}}
|
||||
|
||||
|
||||
## Bibliography {#bibliography}
|
||||
|
||||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||||
<div class="csl-entry"><a id="citeproc_bib_item_1"></a>Ebrahimi-Darkhaneh, Hadi. 2019. “Measurement Error Caused by Self-Heating in Ntc and Ptc Thermistors.” <i>Tex. Instrum. Analog. Des. J. Q</i> 3: 001–7.</div>
|
||||
<div class="csl-entry"><a id="citeproc_bib_item_2"></a>Neto, Joao Brito, Renan Geraldes, Francesco Lena, Marcelo Moraes, Antonio Piccino Neto, Marlon Saveri Silva, and Lucas Volpe. 2022. “Temperature Control for Beamline Precision Systems of Sirius/Lnls.” <i>Proceedings of the 18th International Conference on Accelerator and Large Experimental Physics Control Systems</i> ICALEPCS2021 (nil): China. doi:<a href="https://doi.org/10.18429/JACOW-ICALEPCS2021-WEPV001">10.18429/JACOW-ICALEPCS2021-WEPV001</a>.</div>
|
||||
<div class="csl-entry"><a id="citeproc_bib_item_2"></a>Neto, Joao Brito, Renan Geraldes, Francesco Lena, Marcelo Moraes, Antonio Piccino Neto, Marlon Saveri Silva, and Lucas Volpe. 2022. “Temperature Control for Beamline Precision Systems of Sirius/Lnls.” <i>Proceedings of the 18th International Conference on Accelerator and Large Experimental Physics Control Systems</i> ICALEPCS2021: China. doi:<a href="https://doi.org/10.18429/JACOW-ICALEPCS2021-WEPV001">10.18429/JACOW-ICALEPCS2021-WEPV001</a>.</div>
|
||||
</div>
|
||||
|
51
content/zettels/thermoelectric_cooler.md
Normal file
@@ -0,0 +1,51 @@
|
||||
+++
|
||||
title = "Thermoelectric cooler"
|
||||
author = ["Dehaeze Thomas"]
|
||||
draft = false
|
||||
+++
|
||||
|
||||
Tags
|
||||
: [Temperature Control]({{< relref "temperature_control.md" >}})
|
||||
|
||||
|
||||
## First principles {#first-principles}
|
||||
|
||||
From (<a href="#citeproc_bib_item_1">Evers et al. 2021</a>):
|
||||
|
||||
<a id="figure--fig:thermoelectric-cooler-schematic"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/thermoelectric_cooler_schematic.svg" caption="<span class=\"figure-number\">Figure 1: </span>Schematic of a Peltier module" >}}
|
||||
|
||||
The thermoelectric dynamics is described by 3 phenomena:
|
||||
|
||||
1. the Fourier effect
|
||||
2. Joule heating
|
||||
3. the Peltier effect
|
||||
|
||||
|
||||
### The Fourier effect {#the-fourier-effect}
|
||||
|
||||
The Fourier effect \\(Q\_f\\) describes the energy transfer through **conduction** between the two sides of the Peltier module:
|
||||
\\[ Q\_f^{1 \rightarrow 2} = \frac{K\_m \cdot A}{d} (T\_1 - T\_2) \\]
|
||||
for conduction from temperature \\(T\_1\\) to \\(T\_2\\) with \\(K\_m\\) the conductivity of the Peltier module in \\(W/m \cdot K\\), \\(A\\) the area in \\(m^2\\) and \\(d\\) the thickness in \\(m\\).
|
||||
|
||||
|
||||
### Joule heating {#joule-heating}
|
||||
|
||||
Joule heating \\(Q\_j\\) occurs when an electrical current flows through a resistive element:
|
||||
\\[ Q\_j = R\_m I^2 \\]
|
||||
where \\(R\_m\\) is the electrical resistance in \\(\Omega\\) of the Peltier module and \\(I\\) is the electrical current in \\(A\\).
|
||||
|
||||
|
||||
### The Peltier effect {#the-peltier-effect}
|
||||
|
||||
The Peltier effect describes the occurrence of a heat flow over a semi-conductor in the presence of an electrical potential difference and resulting current:
|
||||
\\[ Q\_p = S\_m T I \\]
|
||||
where \\(S\_m\\) is the Seebeck coefficient of the Peltier module, and \\(T\\) is the temperature at the cold/hot side.
|
||||
|
||||
|
||||
## Bibliography {#bibliography}
|
||||
|
||||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||||
<div class="csl-entry"><a id="citeproc_bib_item_1"></a>Evers, Enzo, Rens Slenders, Rob van Gils, Bram de Jager, and Tom Oomen. 2021. “Thermoelectric Modules in Mechatronic Systems: Temperature-Dependent Modeling and Control.” <i>Mechatronics</i> 79 (nil): 102647. doi:<a href="https://doi.org/10.1016/j.mechatronics.2021.102647">10.1016/j.mechatronics.2021.102647</a>.</div>
|
||||
</div>
|
25
content/zettels/wheatstone_bridge.md
Normal file
@@ -0,0 +1,25 @@
|
||||
+++
|
||||
title = "Wheatstone Bridge"
|
||||
author = ["Dehaeze Thomas"]
|
||||
draft = false
|
||||
+++
|
||||
|
||||
Tags
|
||||
:
|
||||
|
||||
Wheatstone Bridge are used to measure an electrical resistance.
|
||||
|
||||
They are used to read various sensors:
|
||||
|
||||
- [Strain Sensors]({{< relref "strain_sensors.md" >}})
|
||||
- [Temperature Sensors]({{< relref "temperature_sensors.md" >}})
|
||||
|
||||
<a id="figure--fig:wheatstone-bridge"></a>
|
||||
|
||||
{{< figure src="/ox-hugo/wheatstone_bridge.jpg" caption="<span class=\"figure-number\">Figure 1: </span>Electrical schematic of a Wheatstone bridge" >}}
|
||||
|
||||
|
||||
## Bibliography {#bibliography}
|
||||
|
||||
<style>.csl-entry{text-indent: -1.5em; margin-left: 1.5em;}</style><div class="csl-bib-body">
|
||||
</div>
|
BIN
static/ox-hugo/heat_transfer_fluid.png
Normal file
After Width: | Height: | Size: 237 KiB |
BIN
static/ox-hugo/heat_transport.png
Normal file
After Width: | Height: | Size: 71 KiB |
BIN
static/ox-hugo/rotation_stage_rt100bs_asd_errors.png
Normal file
After Width: | Height: | Size: 83 KiB |
BIN
static/ox-hugo/rotation_stage_rt100bs_cas_error.png
Normal file
After Width: | Height: | Size: 58 KiB |
BIN
static/ox-hugo/rotation_stage_rt150u_asd_errors.png
Normal file
After Width: | Height: | Size: 91 KiB |
BIN
static/ox-hugo/rotation_stage_rt150u_cas_error.png
Normal file
After Width: | Height: | Size: 55 KiB |
BIN
static/ox-hugo/rotation_stage_rt150u_xyz_errors.png
Normal file
After Width: | Height: | Size: 53 KiB |
BIN
static/ox-hugo/temperature_betatherm_10K3A1_results.png
Normal file
After Width: | Height: | Size: 569 KiB |
BIN
static/ox-hugo/temperature_effect_wires.png
Normal file
After Width: | Height: | Size: 353 KiB |
BIN
static/ox-hugo/temperature_ntc_dc95f_results.png
Normal file
After Width: | Height: | Size: 237 KiB |
BIN
static/ox-hugo/temperature_seebeck_material_A.png
Normal file
After Width: | Height: | Size: 27 KiB |
BIN
static/ox-hugo/temperature_sensor_lut_errors.png
Normal file
After Width: | Height: | Size: 37 KiB |
BIN
static/ox-hugo/temperature_sensor_lut_errors_ntc.png
Normal file
After Width: | Height: | Size: 28 KiB |
BIN
static/ox-hugo/temperature_sensor_ntc_curve.png
Normal file
After Width: | Height: | Size: 28 KiB |
BIN
static/ox-hugo/temperature_sensor_pt100_curve.png
Normal file
After Width: | Height: | Size: 20 KiB |
BIN
static/ox-hugo/temperature_termocouples_sensitivities.png
Normal file
After Width: | Height: | Size: 189 KiB |
BIN
static/ox-hugo/temperature_thermocouple_meas.png
Normal file
After Width: | Height: | Size: 52 KiB |
BIN
static/ox-hugo/temperature_thermocouple_two_materiasl.png
Normal file
After Width: | Height: | Size: 41 KiB |
2356
static/ox-hugo/thermoelectric_cooler_schematic.svg
Normal file
After Width: | Height: | Size: 686 KiB |
BIN
static/ox-hugo/wheatstone_bridge.jpg
Normal file
After Width: | Height: | Size: 10 KiB |