Определение локальных значений параметров плазмы: концентрации электронов, атомов в основном и возбужденных состояниях, атомной и электронной температуры в некоторых случаях возможно, если определены интенсивности и формы контуров спектральных линий, излучаемых различными элементарными объемами плазмы [1], что и является задачей локальной спектроскопии.
В решении этой задачи можно выделить следующие основные этапы:
- сбор массивов экспериментальных данных, связанных с пространственным распределением спектральной энергетической яркости поверхности источника.
- определение на основе этих данных пространственного распределения спектральной энергетической яркости поверхности источника с учетом кривой спектральной чувствительности использованных фотоприемников и возможных искажений полученных распределений из-за конечного разрешения спектральной аппаратуры.
- переход от массива яркости поверхности источника к массиву локальных значений коэффициентов излучения ε (λ, r) (λ- длина волны, r-пространственная координата) на основе принятой модели источника. В данной работе предполагается осевая симметрия источника, тогда такой переход (радиальное преобразование) состоит в решении интегрального уравнения Абеля .
Принципы и разнообразные технические и программные средства для реализации этих этапов известны давно: сбор данных осуществляется автоматизированной системой (как правило, уникальной) под управлением ЭВМ, для обработки данных с целью радиальных преобразований и исключения аппаратных искажений применяются различные методы статистической регуляризации [2] (с использованием не всегда обоснованной априорной информации об искомых распределениях) .
Отличительная черта предлагаемой технологии состоит в применении программно аппаратного комплекса, построенного по модульному принципу на основе стандартных, промышленно выпускаемых технических и программных средств и обработки данных с помощью алгоритмов, обеспечивающих устойчивость результата к шумам эксперимента на основе внутренней корреляции обрабатываемых массивов.
Измерительный комплекс использует дифракционный спектрометр ДФС-12 с фотоэлектрической регистрацией спектра и систему пространственного стробирования [3] В качестве приборных интерфейсов используются модули КАМАК или платы LabPC-1200AI, PCI-1802L, PCI-1202L, а также сетевой интерфейс (для удаленного управления экспериментальным оборудованием). Все программное обеспечение для управления процессом сбора данных и их обработки разработано в среде графического программирования LabView [4]. Комплекс работает под управлением IBM-совместимого компьютера с ОС Windows-98 , однако легко может быть адаптирован к любой операционной системе , поддерживающей LabView.
Программное обеспечение включает следующие основные функциональные модули:
- а) градуировки спектрального прибора по длинам волн путем регистрации известного спектра,
- б) градуировки чувствительности фотоприемной системы путем регистрации спектра аттестованной температурной лампы.
- в) автоматизированного пространственно- спектрального стробирования в заданных спектральном и пространственном интервалах и в заданных фазах периода тока ( для источников , работающих на переменном токе).
- г) загрузки ранее зарегистрированных данных для последующего анализа,
Каждый модуль оформлен в виде "виртуального прибора", на панели которого имеются окна для ввода необходимой информации, отражается результат работы модуля и имеются кнопки управления (запуск модуля, сохранение результата).
Модуль "г" позволяет автоматически перевести отсчеты регистрирующей системы в абсолютные значения яркости поверхности источника, извлечь различную информацию из спектра (например, измерить положение максимумов и ширины линий) непосредственно используя средства Labview для работы с графиками и выполнить обработку данных с использованием описанного ниже алгоритма.
Результаты работы любого модуля автоматически протоколируются в информационом файле, в который автоматически вводится дата эксперимента, вся информация об условиях измерений, которая была на панелях виртуальных приборов, а также вспомогательная информация, которую пользователь может ввести (например, тип источника, ток разряда).
При работе модуля "в" система осуществляет сканирование по пространственной координате х - расстоянию от оси источника и для определенных положений хk (k=1,2...m) регистрирует распределение отсчетов фотоприемного устройства Fik (относящиеся к определенной фазе тока), соответствующих различным длинам волн (λi) (i=1,2... n)
(1)
Здесь b(λ,хk) - спектральная энергетическая яркость поверхности источника на расстоянии от оси хk, g - аппаратная функция спектрального прибора, q - коэффициент, зависящий от чувствительности фотоприемника, практически постоянен в пределах контура одной спектральной линии, поэтому в дальнейшем описании алгоритма опускается. Для получения искомых значений коэффициентов излучения ε (λ, r) в определенном сечении источника (теперь r- расстояние от его оси) надо исключить аппаратные искажения из каждого профиля Fik при фиксированном k (m раз решить уравнение 1) , с учетом значения q определить b(λ,хk) и найти e (λ, r), решив интегральное уравнение Абеля для каждой спектральной компоненты (n раз)
=R{bλ} (2)
Вместо этого проводится "коллективная" обработка исходного массива Fik, основанная на использовании метода главных компонент [5], при которой существенно сокращается число операций, и главное, каждая операция совершается над распределением, полученным путем усреднения значительного числа отсчетов, что повышает устойчивость результата к шумам.
Сначала находится средний вектор отсчетов, в зависимости от пространственной координаты (его компоненты усреднены по всем длинам волн):
(3)
Затем вычисляется ковариационная матрица массива отсчетов:
(4)
Массив Fik раскладывается по собственным векторам матрицы А, причем учитываются только вектора Up, отвечающие собственным значениям матрицы А, превосходящим величину оценки дисперсии воспроизводимости одного отсчета F (оценивается при повторных измерениях одного и того же массива). Вследствие сильной корреляции отсчетов внутри массива число таких векторов v всегда существенно меньше m.
Вычисляются проекции спектрального распределения на вектор с номером р- Мi,p:
(5)
только они зависят от длин волн и из них, в случае необходимости, надо исключать аппаратные искажения ( v операций), результат такого исключения обозначим Li,p.
F¯k и Up зависят только от координаты х и над ними осуществляются радиальные преобразования (2) ( всего v+1 раз)
В результате искомые контура в различных точках плазмы получаются согласно соотношению:
e (λ, r)= R{F¯ }+ ∑V P=1 Lp(λ) R{ Up} (6)
Комплекс применен для исследования конкретных плазменных объектов. В металлогалоидной лампе (закрытая дуга в парах ртути высокого давления с добавкой иодида таллия, питаемая переменным током) зарегистрированы поперечные профили яркости 12 спектральных линий таллия и ртути в двух фазах тока (соответствующих максимальной и минимальной яркости линий) в m=10 пространственных точках. Число спектральных точек n в каждой линии в зависимости от ее ширины составляло от 85 до 280. Описанный алгоритм обработки применен к линиям, которые по оценкам не испытывали заметного поглощения внутри источника. Ни в одном случае число значимых собственных векторов не превысило v=2. Таким образом, обработка свелась к двум решениям уравнения (1) вместо 10 и трем радиальным преобразованиям (вместо более чем 80) Эффективность коллективной обработки иллюстрируется на рисунке 1.
Рисунок 1. Нормированные на максимум контура линии 434.7нм (наиболее слабой и, соответственно, зашумленной из линий ртути) для центра разряда, полученные при "коллективной" (1) и "индивидуальной" (2) обработке исходного массива в двух фазах тока:
а) максимальная, b) минимальная интенсивность излучения.
По построенным в результате обработки исходных массивов отсчетов определены контурам e (λ, r) найдены радиальные распределения концентрации атомов ртути в основном и возбужденном состояниях, возбужденных атомов таллия и концентрации электронов. При этом обнаружен ряд ранее неизвестных эффектов (например, немонотонный ход определенных по штарковскому уширению линий [1] электронных концентраций), исследование которых с использованием описанной технологии продолжается.
Работа выполнена при поддержке гранта PZ-013-02 CRDF, Мин.обр. РФ и Правительства Карелии.
Литература
- Энциклопедия низкотемпературной плазмы / Под. ред В.Е.Фортова. Вводный том II, М: Наука, 2000, 655с.
- Тихонов А. И., Арсенин В. Я. Методы решения некорректных задач. М.: Наука.- 1979. - 285 с.
- А.с. № 545954; СССР, Оптико-механическое линейносканирующее устройство:/ Кюльмясу И.И., Хахаев А.Д. Щербина А.И. (СССР) ; Заяв.20.11.76; Опубл. 05.02.1977, Бюл.№5
- Жарков Ф.П., Каратаев В.В., Никифоров В.Ф., Панов В.С. Использование виртуальных инструментов LabView.-М: Радио и связь, 1999, 268 с. ( см. также http://www.labview.ru/ )
- J. Edward Jackson. A User´s Guide To Principal Component. 1991. John Willey &Sons , New York , 567 pp.