CC BY
352
УДК 536.531
ИЗМЕРЕНИЕ ТЕМПЕРАТУРЫ ТЕРМИСТОРАМИ
А. В. Серяков
ОАО «Специальное конструкторско-технологическое бюро по релейной технике»
Россия, 173000, Великий Новгород, ул. Нехинская, 55. E-mail: [email protected]
Рассмотрен вопрос измерения температуры полупроводниковыми термисторами и повышения точности этих измерений. Повышение точности осуществлено за счет привязки температуры к точке перегиба (inflection point) Tinf функциональной зависимости логарифма сопротивления термистора и учета временного дрейфа коэффициентов разложения.
Измерение температуры, термисторы, точка перегиба функциональной зависимости логарифма сопротивления термистора, дрейф коэффициентов разложения.
THERMISTOR THERMOMETRY
A. V. Seryakov
"Special Design and Technological Bureau of Relay Technique"
55 Nekhinskaya street, Veliky Novgorod, 173000, Russia. E-mail: [email protected]
In the article we consider the problem of thermistor thermometry and increase of accuracy of such measurements. Such accuracy increase is achieved by means of temperature being attached to the inflection point T|nf of functional dependence of the logarithm of thermistor resistance, and the account of the timing drift of expansion coefficients.
Temperature measurement, thermistors, inflection point of the functional dependence of the logarithm of thermistor resistance, drift of expansion coefficients.
Проблема измерения температуры термисторами, применяемыми при исследовании тепловых процессов в устройствах охлаждения, например в тепловых трубах среднетемпературного диапазона, исключительно актуальна. При разработке новых тепловых труб сравнение и выбор тех или иных конструктивных вариантов проводится с помощью определения теплопередающих характеристик труб, и точность измерений температуры (разностей температур) при этом является решающим фактором.
С этой целью была выполнена калибровка термистора СТ3-19 в интервале температур 0^200 °С, путем сличения показаний в равновесных термодинамических условиях с образцовым платиновым термометром сопротивления [1; 2].
Экспериментальная часть. Термисторы СТ3-19 представляют собой термочувствительные элементы [3-5], выполненные из керамических оксидных материалов на основе никеля, марганца и кобальта, по составу определяемые формулой (№0,2Мп0,7Со0,і)3О4, с отрицательным ТКС и сопротивлением при комнатной температуре около 10 кОм. Бусинка термистора СТ3-19 тонкими платиновыми проволочками (0,01 мм) приварена к выводным траверсам (0-3 мм) и покрыта тонким слоем молибденового стекла. Термисторы имеют радиационную стойкость не хуже 2 Ус и длительную выдержку при больших g.
В качестве термометра сопротивления использован чувствительный элемент образцового платинового термометра сопротивления ПТС-10. Кварцевый
геликоид с платиновой спиралью помещен в выполненный из молибденового стекла тонкостенный цилиндрический чехол диаметром 4 мм. Перед запаиванием чехол с геликоидом был заполнен гелием при давлении, несколько меньшим атмосферного. После изготовления термометр сопротивления заново проградуирован в Сибирском институте метрологии (СНИИМ), г. Новосибирск, в интервале температур 0^200 °С. Абсолютная погрешность измерения температуры термометром оценена в 0,02 К.
Термометр с термистором размещены в массивном медном блоке цилиндрической формы, длиной 100 мм и диаметром 70 мм, установленном в вакуумной камере на специальных жестких подвесах, закрепленных на внутренней поверхности верхнего фланца камеры. В верхний фланец вакуумной камеры вварены две трубы, через одну из которых осуществляется вакуумная откачка камеры, а через другую выведены все измерительные провода. Вакуум в камере поддерживался не хуже 1,3 • 10-3 Па (10-5тор).
Вакуумная камера с медным блоком помещена в жидкостный термостат вместимостью 40 л, температуру в котором устанавливали и длительно поддерживали с помощью регулятора ВРТ-3 и платинового термометра ТСП-Ш в интервале 0 + 200 °С с минимальной дискретностью 5 -10-3 К. Охлаждающим устройством являлся кондиционер БК-1500, штатный испаритель-теплообменник которого был заменен разветвленным медным змеевиком. Рабочей жидко-
стью в термостате служило силиконовое масло ПЭС-В2, позволяющее работать при температурах 0-200 °С. Колебания температуры внутри термостата на протяжении нескольких часов наблюдений не превышали 110-3 К, дрейф температуры был менее 110-4 К/час.
Проведение измерений. На протяжении нескольких лет в в ИТ СО РАН, Новосибирск, была осуществлена калибровка термисторов СТ3-19. Измерения проводились в стационарном режиме с изотермической оболочкой [6; 7], с постепенным увеличением температуры от 0 °С до 200 °С с шагом 10 °С. Длительность одного непрерывного цикла подъема температуры и проведения измерений достигала 48 часов. Всего был проведен 21 цикл измерений.
Собственно калибровка термистора состояла в точном измерении стандартным потенциометрическим методом сопротивления термистора Кс, Ом, в стационарном состоянии при заданной температуре Т, К, определяемой по платиновому термометру сопротивления ПТС-10. В качестве измерителя использован компаратор напряжений Р3003 класса точности 0,0005, вместе с образцовой катушкой сопротивления Р 321 класса точности 0,01, размещенной в пассивном термостате. Источником тока служил набор батарей «Бакен» в заземленном металлическом кожухе.
Всего получено более 500 экспериментальных точек, из которых был сформирован массив исходных данных температур Т и логарифмов сопротивлений термистора 1пКс. Максимальная случайная погрешность измерения температуры термометром ПТС-10 не превышала (2-3) • 10-3 К, сопротивлений термистора 5 10-4 Ом.
Все измерения температуры, выполненные на термометре сопротивления ПТС-10, включая калибровку в СНИИМ, сделаны при величине измерительного тока 1 та (1 • 10-3 А), и выделяемая тепловая мощность равна Жптс = (10-17)10-6 —.
На рис. 1 приведены результаты одного цикла измерений сопротивления термистора ИС в зависимости от температуры, при этом для удобства дальнейшего анализа график представлен в виде функции обратной
температуры 10 /Т, 1/К, от 1пЯС. Эта зависимость близка к линейной, но в области повышенных температур Т ~ 473К, где 1пКс ~ 4-5, и в области пониженных температур Т ~ 273К, где 1пКс ~ 9-10, наблюдаются отклонения. При средних температурах, в интервале значений 1пКс ~ 6-8, возможно существование точки перегиба.
Рассматривая датчик температуры как систему со сосредоточенными параметрами, уравнение термистора записывают следующим образом [1; 2]:
Тс +тс
Кс
Кс
(1)
где Тс - температура термистора, К; Тс - производная по времени температуры термистора при динамическом режиме измерений, тс - характерное время запаздывания термистора, 5; Жс - выделяемая на термисторе электрическая мощность, Ж; Кс - коэффициент теплопередачи между термистором и медным блоком, Ж/К; Т0 - измеряемая термометром ПТС-10 температура медного блока или, например, тепловой трубы т -время, 5; сс - теплоемкость термистора, I /К.
Решением уравнения (1) является выражение (2) [8], в котором в момент времени т = т
(
Тс (т) — Тс • ехР
т-т
Ес
V Кс
- + Тп
1 - ехр
Ґ
т-т
(2)
температура термистора считается известной и рави /7"т ^
ной Тс .
Электрическая мощность Жс, выделяемая на термисторе измерительным током, при всех калибровочных испытаниях постоянна и равна 20 • 10-6 —, и перегрев термистора ДТс вычисляют по формуле
ДТс —
Ес
К
— тс
Ес
с
(3)
Рис. 1. Экспериментальная зависимость логарифма сопротивления термистора СТ3-19 от обратной температуры 104/Т, 1/К. Прямая линия - расчет по уравнению (6); график кривой - расчет по уравнению (9)
Рис. 2. Экспериментально определенные времена Tc запаздывания термистора СТ3-19
Для определения перегрева термистора были проведены специальные опыты, в которых в стационарных состояниях при температурах от 2 до 195 °С, измерялся темп нагрева и остывания термистора (и термометра) при ступенчатом изменении выделяемой в них электрической мощности Wc. С помощью измерительной системы жидкостного калориметра, описание которой приведено в [6; 7], были проведены подробные измерения релаксационных характеристик датчиков. Расчет времени запаздывания термистора тС по уравнению (2) с помощью стандартной процедуры метода наименьших квадратов (МНК) [9-12] дал результаты, приведенные на рис. 2.
Аппроксимирующая кривая на рис. 2 задана полиномиальным выражением
тС (t) = 3,3963027 -10-5 • t2 -
-1,39779 -10-2 • t + 2,261388 , мин, (4)
где t - температура в градусах Цельсия, среднеквадратичное отклонение экспериментальных точек с ~ 0,18 мин.
Теплоемкость термистора Cc ~ (0,3-0,35) Дж/К, коэффициент теплопередачи Kc ~ (2-5) 10-3 W/K, величина перегрева измерительным током не превышает ATC ~ (10-4) • 10-3 К и учитывается при всех измерениях.
Разность температур между медным блоком и окружающей ее изотермической оболочкой - термостатом - не превышала 0,05-0,1 К, величина температурного хода при проведении калибровки была меньше 10-7 К/s, поэтому термодинамическое состояние датчиков температуры внутри медного блока в этот период являлось квазистационарным.
Обработка результатов. Функциональная зависимость электрического сопротивления оксидного полупроводникового термистора Rc от температуры Т достаточно сложна, и в первом приближении ее представляют в виде сопротивления идеального полупроводника со строго одинаковым количеством дырок и носителей заряда, в экспоненциальном виде
RC = А • exp ГB1, (5)
где Яс - электрическое сопротивление термистора, Ом, при температуре Т, К; А - величина сопротивления термистора, Ом, при бесконечно большой температуре; В - параметр чувствительности термистора, в общем случае также зависящий от температуры, 1/К. При температуре / = 1 °С сопротивление термистора СТ3-19 ИС ~ 30 кОм, при температуре / = 200 °С Яс ~ 50 Ом. Величина В ~ 4000 1/К, А ~ 0,013 Ом. После логарифмирования (5) получают
1
1
1
— =-----------------ln А +----------------ln R
T
B
B
(6)
Для уточнения вопроса о перегибе, была вычислена производная d(1/T)/d(lnRC), при этом анализировался весь массив экспериментальных точек. Значения производной d(1/T)/d(lnRC), рассчитанные по результатам одного цикла измерений сопротивления термистора в зависимости от lnRC, приведены на рис. 3.
По результатам численного дифференцирования [13] всего массива экспериментальных точек, минимум кривой определен нами при значении ln Rc min =
7,63 ± 0,01, что соответствует температуре точки перегиба (inflection point) Tinf = 336,34 К или 63,19 °С. Разлагая производную d(1/T)/d(lnRC) в ряд в окрестности точки минимума, было получено выражение
d(1/T^ = a1 + a2 •(lnRC -7,63)2 +
d ln RC
+ a3 • (lnRc -7,63) +...
(7)
где а - коэффициенты разложения.
После интегрирования разложение (7) приводят к следующему виду
Т = А0 )+ А1 (^ )-(1п Яс - 7,63) + А2 (^ )х
;(lnRc - 7,63)3 + A3 (Xd )ln Rc - 7,63)4
(8)
где A,(xd) - коэффициенты разложения; Td - время дрейфа (drift) коэффициентов.
Рис. 3. Расчетная зависимость производной —С(1/Т)/(С(1пЛс)10 ,1/К, в зависимости от логарифма сопротивления 1пЯс термистора СТ3-19
Таким образом, калибровка термистора СТ3-19 по существу, состоит в определении численного значения коэффициентов полинома (8) А,(тй) и учета их дрейфа во времени.
Основная особенность расчета по уравнению (8) заключается в привязке температуры к точке перегиба Тм функциональной зависимости логарифма сопротивления термистора.
Расчет коэффициентов А,(тй) был проведен с помощью стандартной процедуры метода наименьших квадратов (МНК) [9-12]. Абсолютные погрешности расчета коэффициентов А,(тй) следующие: 5А0 —1 • 10-3, §А —1 • 103, §А2 —1 • 10-4, ЪАъ —1 • 105.
Для контроля качества калибровок и их долговременной стабильности после каждого из 21 определения набора коэффициентов А,(тй) из выражения (8) были вычислены производные й(1/Т)/й(1пЯс). С разбросом не более 5 10-7 все значения производных лежат на кривой рис. 3. В области пониженных температур при значениях логарифмов 1пЯс — 9-10, разброс калибровок несколько больше и достигает X — (5-7) •Ш-7. Связано это с тем, что там, где больше сопротивление термистора Яс и круче его температурная зависимость йЯс/йТ, там меньшая плотность экспериментальных точек, поскольку калибровки термистора преимущественно проводились равномерно по температуре с шагом 10 °С.
Динамика изменения коэффициентов А,(тй), временной дрейф калибровки термистора представлена на рис. 4. Первые восемь циклов измерений проводились с нагревом термистора до температуры 200 °С. Следствием такого нагрева был значительный дрейф калибровки и изменение коэффициентов А,(тй), например йА0/йт — 6 • 10-4 1/мес.
После того как верхний предел температуры был снижен до 190 °С, временной дрейф коэффициентов заметно уменьшился, к примеру скорость изменения коэффициента А0(тй) стала всего йА^йт — 1,5 10-4 1/мес. Значения временного дрейфа коэффициентов А(тй) приведены ниже, тй - месяцы.
Рис. 4. Динамика изменений коэффициентов А1 во времени при нагревах до 200 °С - левые части графиков, до 190 °С - правые части графиков
Коэффициент А0(тй):
А0 = - 6,21997 10-4 та + 29,824488 при нагреве до 200 °С;
А0 = - 2,3075444 10-4 та + 29,819432 при нагреве до 190 °С.
Коэффициент А^):
А! = 4,051207 10-5 тй + 2,489278 при нагреве до 200 °С;
Aj = 1,5876991 10 5-Td + 2,48958 при нагреве до 190 °С.
Коэффициент A2(xd):
A2 = - 2,2266277-10-5-Td + 0,0022733 при нагреве до 200 °С;
A2 = - 1,0559017-10-5-Td + 0,0021054 при нагреве до 190 °С.
Коэффициент A3(xd):
A3 = - 3,98635-10-8 • (xd)2 + 1,771915• 10-6 •id + 6.3241 • 10-5 при нагреве до 190 °С и до 200 °С.
Подстановка коэффициентов Ai(Td), вычисленных при нагреве термистора до 190 °С в уравнение (8), снижает разброс производных d(1IT)Id(\nRC) на кривой рис. 3 до величины X ~ (2-3)-10-7, и позволяет более точно определить точку минимума \nRCmin = 7,632 ±0,01, тем самым уменьшить погрешность расчета температуры при проведении измерений термистором СТ3-19.
Таким образом, рекомендуемое уравнение для расчета температуры термистором СТ3-19, учитывающее как точку перегиба (inflection point)
Tinf = 336, 34 К функциональной зависимости логарифма \nR сопротивления термистора, так и временной дрейф полиномиальных коэффициентов разложения A,(xd) при периодических нагревах термистора до 190 °С, выглядит следующим образом:
T=Л (Td)+ л (Td )-(\nRc - 7,632) +
A2 (Td ) • (\n Rc - 7,632 )3 + A3 (Td ) • (\nRc - 7,632 )4 . (9)
Обсуждение результатов. Известно полиномиальное уравнение третьей степени J. S. Steinhart, S. R. Hart [14], предназначенное для расчета температуры термистором и содержащее линейный и кубический по логарифму сопротивления \nRC члены
1 3
----= A + B • \n RC + C • \n RC. (10)
Диапазон применения по температуре данного уравнения с одним набором коэффициентов не превышает 50 -75 К [14; 15]. Используя опубликованные в современных интернет-изданиях [15] наборы численных значений коэффициентов уравнения Steinhart'a и Hart'а, было проведено сопоставление результатов расчета температуры по уравнению Steinhart'a и Hart'a и по нашему уравнению четвертой степени логарифмов сопротивления термистора (9).
Относительные разности температур
5 = (TSH - T) IT-100 %,
в зависимости от \nRC представлены на рис. 5. Проведенное сопоставление результатов показывает, что разности температур 5 в зависимости от \nRC имеют знакопеременный характер.
Минимальные значения разностей температур 5 порядка 0,25 %, наблюдаются при величине \nRC ~ 6 термистора СТ3-19. Максимальные значения разностей 5, достигающие 0,8 -1 %, имеют место на краях интервала применимости уравнения (9).
Погрешность определения температуры (разностей температур) термистором СТ3-19 с помощью уравнения четвертой степени (9) с учетом точки перегиба при \nRC = 7,632, не превышает (3-5)-10-4 К.
Применение полиномиального уравнения четвертой степени (9) с привязкой температуры к точке перегиба функциональной зависимости логарифма сопротивления термистора, позволяет расширить диапазон и повысить точность определения температуры.
Учет временного дрейфа коэффициентов разложения полинома (9) A,(xd) позволяет повысить точность определения точки минимума \nRCmin = 7,632 ±0,01
характеристик термистора и тем самым точность определения температуры.
Применение полиномиального уравнения третьей степени Steinhart'a и Hart'a увеличивает погрешности расчета температуры термистором на краях интервала 273 К - 473 К до 1,5-2 К.
Рис. 5. Относительные разности температур 5, вычисленные по уравнению 81ет!1аг1'а и НагТа, и по уравнению четвертой степени относительно точки перегиба характеристик термистора (9)
SH
Библиографические ссылки
1. Линевег Ф. Измерение температур в технике. Справочник. М. : Металлургия, 1980.
2. Кондратьев Г. М. Тепловые измерения. М. : Машгиз, 1957.
3. Мэклин Э. Д. Терморезисторы. М. : Радио и связь, 1983.
4. Гендин Г. С. Все о резисторах : справочное издание. М. : Горячая Линия-Телеком, 2000.
5. Шефтель И. Т. Терморезисторы. М. : Наука, 1973.
6. Серяков А. В. Автоматизированный калориметр постоянного объема для исследования растворов электролитов // Теплофизические свойства растворов, расплавов и композитов : сб. научн. тр. ; Институт теплофизики СО РАН. Новосибирск, 1991. С. 139-153.
7. Груздев В. А., Серяков А. В. Жидкостный калориметр // Теплофизические свойства веществ : тр. 8-й Всесоюзн. конф. Новосибирск, 1989. Ч. 1. С. 256-261.
8. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М. : Наука, 1976.
9. Демидович Б. П., Марон И. А. Основы вычислительной математики. М. : Наука, 1966.
10. Амосов А. А., Дубинский Ю. А., Копченова Н. В. Вычислительные методы для инженеров. М. : Высшая школа, 1994.
11. Хемминг Р. В. Численные методы. М. : Наука, 1968.
12. Дэвенпорт Дж., Сирэ И., Турнье Э. Компьютерная алгебра: Системы и алгоритмы алгебраических вычислений. М. : Мир, 1991.
13. Мудров А. Е. Численные методы для ПЭВМ на языках Бейсик, Фортран и Паскаль. Томск : Раско, 1991.
14. Steinhart J. S., Hart S. R. Ca\ibration curves for termistors II Deep Sea Research and Oceanographic Abstracts. 1968. Vo\. 15. № 4. Р. 497-503.
15. Cornerstone Sensors. A, B, С Coefficients for Steinhart-Hart Equation. Temperature Sensors for Hea\th, Science and Industry. 2010. 2 p. Ca\ifornia 92083. USA.
© А. В. Серяков, 2013
УДК 628.16.087
ОЦЕНКА ВЛИЯНИЯ КОАГУЛЯЦИИ НА СТЕПЕНЬ ОЧИСТКИ ВОДЫ ОТ ИОНОВ МЕТАЛЛОВ
И. Я. Шестаков1, О. В. Раева2
'Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева 660014, г. Красноярск, просп. им. газ. «Красноярский рабочий», 31. E-mai\: yakov\evish@mai\.ru 2Сибирский федеральный университет, институт цветных металлов и материаловедения 660025, Красноярск, просп. им. газ. «Красноярский рабочий», 95. E-mai\: O.V.Raeva@mai\.ru
Представлены результаты экспериментальных исследований очистки воды, заключающейся в пропускании через нее переменного асимметричного тока с использованием нерастворимых разнородных электродов, введении коагулянта и отстаивании.
Ключевые слова: вода, электрохимический способ, переменный ток, коагулянт, отстаивание.
ASSESSMENT OF COAGULATION INFLUENCE ON INDEX OF PURIFICATION OF WATER FROM METAL IONS
I. Ya. Shestakov1, O. V. Raeva2
1Siberian State Aerospace University named after academician M. F. Reshetnev 31 “Krasnoyarskiy Rabochiy” prospect, Krasnoyarsk,660014, Russia. E-mai\: yakov\evish@mai\.ru
2Siberian Federa\ University, Institute of Non-Ferrous Meta\s and Materia\ Science 95 “Krasnoyarskiy Rabochiy” prospect, Krasnoyarsk, 660025, Russia. E-mai\: O.V.Raeva@mai\.ru
The article presents the results of experimental studies of water purification, consisting ofpassing through it asymmetric alternating current with application of insoluble heterogeneous electrodes, introducing coagulating agent and settling.
Keywords: water, electrochemical method, alternating current, coagulating agent, settling.
Загрязнение воды является одной из острейших источники. В большой мере причиной тому являются
экологических проблем в мире. Более 90 % сточных сточные воды гальванических производств, содержа-
вод сбрасываются в открытые водоемы без предвари- щие целый ряд ионов металлов, пагубно воздейст-
тельной очистки. Загрязнены реки, озера, подземные вующих на здоровье людей и окружающую среду.
