Термофизические основы радиометрических измерений температуры
Оптические методы регистрации светового излучения для измерения температуры [3-7] могут применяться как в диапазоне сравнительно невысоких температур (здесь меньше температуры плавления хрома), так и высоких температур. Вполне актуальным является развитие термопирометрии [1] для измерения действительных температур. Под термической пирометрией здесь понимается «радиометрическая» термометрия спектрального и термического отношения, для которой ниже изложены вычислительные формулы.
Одними из возможных устройств измерения электромагнитного излучения являются спектральные фотометры или радиометры [3].
Цель работы заключается в теоретическом представлении возможности измерения температуры с использованием калиброванных устройств основанных на измерении теплового электромагнитного излучения, регистрируемого в полихроматическом или спектральном фотометре или радиометре.
1. Тепловое излучение
Действительные температуры измеряются по оптическому излучению с использованием свойств функции Планка, описывающей тепловую спектральную мощность абсолютно черного тела [5]. Спектральная мощность теплового излучения черного тела в единице спектра длин волн описывается обобщенной формулой Планка:
Iλ = ӕC1/{λ5[exp(C/(λT)) - 1]};
k = 1/λ;
ӕ(λ, Т) = ε(λ, Т)·α(λ, Т);
ӕ ≈ ӕb, (1)
где Iλ - удельная спектральная мощность по длинам волн (Вт·мкм4·см-2); С1 = 37417,7 Вт·мкм4/см2; С = С2 = 14387,8 мкм·К; T - температура; λ - длина волны; k - волновое число; e - коэффициент черноты излучения; α - коэффициент передачи (пропускания); ӕ - коэффициент преобразования («радиоматический» коэффициент) для регистрируемого излучения, определяемый произведением коэффициента излучения наблюдаемого объекта и коэффициента передачи (пропускания) измеряемого электромагнитного излучения через среду наблюдения. ӕb - ниже изложенное модельное приближенное значение коэффициента ӕ.
Измерение приближенной «радиометрической» температуры допускает использование некоторого модельного значения ӕb для неизвестного «радиоматического» коэффициента ӕ.
Преобразуем уравнение (1) к следующей формуле для измеряемой температуры:
T = C(1/λ - 1/λo)/{ln[C1ӕ(λ, Т)/(λ5Iλ(λ)) + 1] - ln[C1ӕ(λo, Т)/(λo5Iλ(λo)) + 1]}; (2)
T = C/{λ·ln[C1ӕ(λ, Т)/(λ5Iλ(λ)) + 1]},
где λo - некоторая отсчетная длина волны для спектрального диапазона.
Здесь вторая формула есть частный случай для первого выражения, если выполняется приближение большой отсчетной длины волны в сравнении с используемыми λo >> λ; ko << k.
Измерение температур по тепловым излучениям наблюдаемого объекта с использованием фотоприемников на основе из кремния или арсенида индия обычно проводят в диапазонах инфpaкрасного, видимого и ультрафиолетового излучения на длинах волн λ ≈ 0,2-2,5 мкм (k ≈ 0,4-5 мкм-1), которые соответствует области фоточувствительности кремния или арсенида индия. Максимум спектральной мощности Iλ находится на длине волны равной 2899,9 мкм·К/T, где T - температура объекта. Наиболее значительные изменения спектральной мощности теплового излучения происходят в коротковолновой области спектра.
2. Коэффициент излучения
Спектры теплового излучения реальных объектов отличаются от спектра черного тела, эти отличия описываются дополнительным коэффициентом излучения 0 £ e(λ, Т) £ 1, являющимся, в общем случае, функцией зависящей от длины волны и температуры и называемым коэффициентом черноты излучения (излучательной способностью или коэффициентом черноты). Коэффициент излучения e(λ, Т) измеряют экспериментально, и при этом в общем случае для измерения коэффициента излучательной способности требуется дополнительный метод измерения. Экспериментальные данные по коэффициентам излучательной способности разных материалов содержатся в справочнике [2].
Если предположить единичность «радиоматического» (radiomatical) коэффициента для регистрируемого излучения (ӕ = 1), то вместо формулы (2) для измеряемой температуры получим цветовой и яркостной аналог для регистрируемой температуры.
Преобразуем уравнение (2) к следующему определению цветовой температуры:
Tc = C(1/λ - 1/λo)/{ln[C1/(λ5Iλ(λ)) + 1] - - ln[C1 /(λo5Iλ(λo)) + 1]}. (3)
где Tc - цветовая температура.
Решение (3) преобразуем к следующему определению для яркостной температуры:
Tr = C/{λ·ln[C1/(λ5Iλ(λ)) + 1]}, (4)
где Tr - яркостная температура.
В простой оптической пирометрии измерения интенсивности излучения проводят обычно на одной (4) или двух (3) эффективных длинах волн [5-7], соответственных некоторым участкам спектра, используемых в пирометрах частичного излучения. В пирометрах полного излучения измеряется «излучательная» температура.
3. Вычисление температуры
Рассмотрим исходную формулу (1) в преобразованных координатах (B, λ).
B = λP/C; P = ln[1 + (ӕC1)/(λ5Iλ(λ))], (5)
где B = 1/T - значение обратной величины действительной температуры.
Также рассмотрим полученную формулу (5) в преобразованных координатах (B, λ) в коротковолновом приближении (коротковолновой асимптотике).
B = λP/C; (6)
P = ln{(ӕC1)/(λ5Iλ(λ))} - - Σj (-1)j{(λ5Iλ(λ))/(ӕC1)}j/j;
1 ≤ j ≤ ∞ ;
(λ5Iλ(λ))/(ӕC1) << 1,
где Σj - знак суммирования по натуральным числам.
Если провести измерения интенсивности на многих длинах волн и задать решение для многих температур для каждой длины волны λ1, λ2 и так далее, то результирующая средняя величина по спектру для обратных температур есть интегральное решение для величины действительной обратной температуры. Значение функции B(λ) в указанных координатах является постоянной величиной, не зависящей от длины волны.
Эта постоянная равна обратной величине температуры. В реальных измерениях график этой функции регистрируется с дополнительными систематическими и статистическими погрешностями и неопределенностями. Для вычисления среднего значения измеряемой температуры используется метод наименьших квадратов [3] или метод усреднения.
4. Метод измерения температуры
Для калибровки пирометров необходимо откалибровать спектральные радиометры и фотометры, которые используются для измерения спектральной мощности. Эта калибровка относится к аппаратному и программному обеспечению спектрометров для обеспечения требуемой зависимости выходного сигнала от интенсивности регистрируемого излучения.
Для калибровки спектральных пирометров можно использовать термостат с тепловым электромагнитным излучателем с применением эквивалента абсолютно черного тела (АЧТ). Способ предполагает сравнение показаний поверяемого измерительного прибора с образцовым.
Перепишем уравнение (2) в ином представлении от обратной температуры.
B = {ln[C1ӕ(λ,B)/(λ5Iλ(λ)) + 1] - ln[C1ӕ(λo,B)/(λo5Iλ(λo)) + 1]}/{C(1/λ - 1/λo)};
B = 1/T. (7)
С целью перехода к коротковолновому приближению преобразуем выражение (7) к следующему уравнению.
B = Bs + {ln[1 + λ5Iλ(λ)/(C1ӕ(λ, B))] - ln[1 + λo5Iλ(λo)/(C1ӕ(λo,B))]}/{C(1/λ - 1/λo)};
Bs = {ln[λo5Iλ(λo)/λ5Iλ(λ)] + ln[ӕ(λ,B)/ӕ(λo,B)]}/{C(1/λ - 1/λo)}; B = 1/T, (8)
где Bs - значение обратной величины некоторой температуры, которая равна действительной обратной температуре в коротковолновом приближении.
Рассмотрим простейшую коротковолновую асимптотику в соответствии с формулой (6) для вычислений температуры.
B = (P - Po)/{C(1/λ - 1/λo)};
P = ln{(ӕC1)/(λ5Iλ)};
B = Bd + [ln ӕӕ]/[C(1/λ - 1/λo)];
ӕӕ = ӕ(λ, B)/ӕ(λo,B); (9)
Bd = {ln[(λo5Iλ(λo))/(λ5Iλ(λ))]}/{C(1/λ - 1/λo)},
где Bd - значение обратной величины простой «радиометрической» температуры; ӕӕ - коэффициент спектрального отношения для радиоматического коэффициента; |ln[ӕӕ(B)]/[C(1/λ - 1/λo)]| - модуль неопределенности логарифма радиоматического коэффициента нормированного на ширину спектра волновых чисел.
Здесь предполагается усреднение по спектру для простой «радиометрической» обратной температуры и для действительной обратной температуры.
В модельном приближении формула (9) приводит к упрощенному вычислению «радиометрической» обратной температуры с усреднением по спектру длин волн.
B = (P - Po)/{C(1/λ - 1/λo)};
P = ln{(ӕC1)/(λ5Iλ)};
B = Bb + [ln ӕӕb]/[C(1/λ - 1/λo)];
ӕӕb = ӕb(λ,B)/ӕb(λo,B); (10)
Bb = {ln[(λo5Iλ(λo))/(λ5Iλ(λ))]}/{C(1/λ - 1/λo)};
|ln[ӕӕ(B)/ӕӕb(B)]/[C(1/λ - 1/λo)]| << B,
где Bb - значение обратной величины «модельной радиометрической» температуры; |ln[ӕӕ(B)/ӕӕb(B)]/[C(1/λ - 1/λo)]| - модуль неопределенности логарифма коэффициента спектрального отношения нормированного на ширину спектра волновых чисел.
Неопределенность коэффициента спектрального отношения уменьшается при использовании усиления коэффициента излучения [7]. Дополнительно при использовании более широкой области волновых чисел регистрируемого спектра, неопределенность коэффициента спектрального отношения приводит к меньшей относительной неопределенности обратной температуры [4], и в пределе получаем действительное значение обратной температуры [3].
Результирующая формула вычислений (10), соответствующая модельному представлению, позволяет вычислять температуру с некоторой точностью и относительной погрешностью, определяемой выбором параметров модели.
Рассмотрим вопрос о методе измерения температуры с использованием некоторого дополнительного образцового термометра. Допустим, что определение действительной температуры может быть проведено в два этапа: вначале калибровочное измерение контрольной температуры контактным или альтернативным способом (например, бесконтактным аналогом контактного) в некоторой контрольной точке, а затем бесконтактное измерение термическим пирометром разности радиометрической и контрольной температур в области визирования. В соответствии с решением (9) должна использоваться следующая формула.
B = (Bp - Bn) + Bm +
+ln[ӕӕ(B)/ӕӕ(Bm)]/[C(1/λ - 1/λo)], (11)
где B - действительное значение обратной температуры; Bp - измеренное значение обратной величины температуры; Bn - измеренное контрольное значение обратной величины температуры; Bm - калибровочное контрольное значение обратной величины температуры; ӕӕ(B)/ӕӕ(Bm) - коэффициент термического отношения коэффициентов спектрального отношения. Здесь предполагается некоторое усреднение по спектру и по выборке значений контрольных температур с некоторыми весовыми коэффициентами. Действительное значение температуры определяется интерполяцией по выборке контрольных температур.
Такой метод измерения температур называется произведенным в операциях, которые происходят с предварительным или дополнительным регистрированием состояния системы.
Допустим, что коэффициент термического отношения, входящий в определение температуры (11), близок к модельной величине, или измерен достаточно широкий интервал волновых чисел. Тогда различие между реальным и модельным значением последнего слагаемого в (11) близко к нулю в сравнении с обратной температурой и получаем следующую упрощенную формулу.
B = (Bp - Bn) + Bm + ln[ӕӕb(B)/ӕӕb(Bm)]/[C(1/λ - 1/λo)];
|ln[{ӕӕ(B)/ӕӕ(Bm)}/{ӕӕb(B)/ӕӕb(Bm)}]/[C(1/λ - 1/λo)]| << B. (12)
Здесь возможно проведение в модельном приближении (10) измерений радиометрических обратных температур, усредняемых по спектру и выборке температур.
Результирующая относительная погрешность измерений действительной температуры может уменьшиться при использовании дополнительной штатной калибровки (12). Следовательно, в общем случае результирующее штатное измерение с калибровкой (12) соответствует псевдомодельной «коллекционной» радиометрической температуре, где для вычисления температуры посредством калибровки учтена некоторая относительная неопределенность радиоматического коэффициента в формуле (9).
Для примера рассмотрим оценку неравенства в последней формуле в (12) для двух способов (модель {ӕӕb(B)/ӕ?ӕb(Bm) = 1} «единичная» или {ӕӕ(B)/ӕӕ(Bm) ≠ 1} справочная, измеренная [2]) при температурах 1400 К и 2000 К при нагревании циркония и измерении спектра светового излучения по нормали к поверхности образца в диапазоне длин волн 500-1000 нм. Вычисления произведем при отсчетной длине волны равной 750 нм. Данный выбор обусловлен справочными данными для коэффициента излучения циркония, представленными в [2] и изображенными на рисунке.
Коэффициент излучения по нормали с поверхности циркония для двух температур [2]: 1 - 1400 К и 2 - 2000 К
Для упрощения задачи предположим, что поглощение света в оптическом тpaкте измерения равно нулю. Кроме того предположим, что статистическая погрешность измерений близка к нулю. Также предположим, что тpaкт измерения света имеет требуемую аппаратную зависимость выходного сигнала от интенсивности света.
После расчетов неравенства в последней формуле в (12) получаем величину различия между справочным и модельным значениями «радиометрической» обратной температуры ≈ 1/1400 К-1. Это неравенство выполнено и усреднено равно ≈ 1/3700000 К-1 << 1/1400 К-1, что является приемлемым для «коллекционного» пирометрического метода в данном примере.
Вывод
В радиометрической термометрии возможно измерение температуры с использованием коэффициентов термического отношения. Для меняющихся «радиоматических» коэффициентов результат вычислений, соответствует приближенной «радиометрической» температуре.
Список литературы
- Дубас Л.Г. Термопирометрия для измерения действительной температуры. // Бесконтактная термопирометрия для плотного вещества. - URL: http://econf.rae.ru/article/6558 (дата обращения 14.02.2012) - 6 с.
- Латыев Л.Н., и др. Излучательные свойства твердых материалов; под ред. А.Е. Шейндлина. - М.: Энергия, 1974. - 471 с.
- Daniel L.P.Ng. Многоволновой пирометр для серых и не серых поверхностей в присутствии мешающего излучения. Патент США № 5 326 172 МПК : G01J5/00; G01J5/52; G01J5/62.
- Felice R.A. Expert System Spectropyrometer Results for non-black, non-grey, or Changing Emissivity and Selectively Absorbing Environments. Electro-Techno- Exposition. - M., 2003. - Nov. 2-4, - 11 p.
- Quinn T.J. Temperature // London: Academic Press Ltd., 1983. - P. 384. (Куинн Т. Температура. - М.: Мир, 1985. - 448 с.)
- Radiometric temperature measurements. I. Fundamentals. Z.M.Zhang, etc., Eds. / Experimental Methods in the Physical Sciences. - Vol. 42. Amsterdam: Elsevier, 2009.
- Radiometric temperature measurements. II. Applications. Z.M. Zhang, etc., Eds. / Experimental Methods in the Physical Sciences. - Vol. 43. Amsterdam: Elsevier, 2010.
Статья в формате PDF 210 KB...
11 12 2024 1:48:25
Статья в формате PDF 153 KB...
10 12 2024 8:10:26
Инженерная рационализация лесопользования предполагает активное применение достижений древесиноведения. Фундаментальные достижения в этой области вполне могут быть применены в исследованиях свойств живой древесины растущих деревьев. Доказательство биотехнического принципа в данной статье выполнено на основе моделирования экспериментальных данных профессора Б.Н.Уголева по деформативности древесины при действии усилий поперек волокон. ...
09 12 2024 5:51:23
Статья в формате PDF 123 KB...
08 12 2024 19:31:38
Статья в формате PDF 138 KB...
07 12 2024 6:16:32
06 12 2024 16:59:35
Статья в формате PDF 119 KB...
05 12 2024 0:20:24
Статья в формате PDF 221 KB...
03 12 2024 8:46:49
Статья в формате PDF 110 KB...
02 12 2024 13:39:56
30 11 2024 1:35:23
Статья в формате PDF 116 KB...
29 11 2024 10:29:51
Статья в формате PDF 102 KB...
28 11 2024 17:38:20
Статья в формате PDF 111 KB...
27 11 2024 9:57:23
Обследовано 65 больных шизофренией, направленных в психиатрический стационар в недобровольном порядке. Установлено, что только 16% из них госпитализируются по Постановлению суда (что составляет 0,46% от всех поступивших больных в течение года), остальные 84% дают согласие на госпитализацию в первые дни пребывания в стационаре. У большинства больных указанного контингента наблюдается выраженная социально-трудовая дезадаптация и низкая комплаентность. Приводятся рекомендации по стратегии и тактике лечения в условиях стационара и внебольничного звена психиатрической помощи, требования к преемственности и особенности диспансерного наблюдения. ...
25 11 2024 2:42:23
Статья в формате PDF 164 KB...
24 11 2024 11:36:22
Статья в формате PDF 307 KB...
23 11 2024 12:51:23
Статья в формате PDF 124 KB...
22 11 2024 15:26:32
Статья в формате PDF 126 KB...
21 11 2024 13:46:59
Статья в формате PDF 204 KB...
20 11 2024 6:54:13
Статья в формате PDF 106 KB...
19 11 2024 9:39:47
Статья в формате PDF 493 KB...
18 11 2024 19:12:48
Статья в формате PDF 115 KB...
17 11 2024 9:35:23
16 11 2024 7:49:17
Статья в формате PDF 142 KB...
15 11 2024 14:42:41
Статья в формате PDF 250 KB...
14 11 2024 5:28:43
13 11 2024 12:35:59
Статья в формате PDF 109 KB...
12 11 2024 22:40:15
Статья в формате PDF 303 KB...
11 11 2024 4:45:32
Статья в формате PDF 143 KB...
09 11 2024 11:18:28
Статья в формате PDF 110 KB...
08 11 2024 15:25:37
Статья в формате PDF 232 KB...
07 11 2024 18:52:38
Замачивание семян и опрыскивание вегетирующих растений томата растворами сочетаний витаминов: пантотеновая кислота–тиамин и фитогормонов: цитокинин–гибберелловая кислота, и совместным их сочетанием снижает токсическое действие хлоридного засоления, повышая всхожесть семян, рост проростков, стeбля, размеры листьев, интенсивность фотосинтеза и накопление общего белка. Наиболее эффективно во всех случаях комплексное сочетание витаминов с фитогормонами. ...
06 11 2024 1:37:33
Статья в формате PDF 112 KB...
05 11 2024 0:23:59
Статья в формате PDF 112 KB...
02 11 2024 7:18:17
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::