Документ взят из кэша поисковой машины. Адрес оригинального документа : http://medphys.phys.msu.ru/pract/dopler.doc Дата изменения: Thu Oct 20 18:57:28 2011 Дата индексирования: Mon Oct 1 19:54:20 2012 Кодировка: koi8-r

Задача 6. Лазерный доплеровский анемометр.

Цель задачи
Целью работы является ознакомление студентов с основными принципами
лазерной доплеровской анемометрии на примере измерения скорости вращения
светорассеивающего диска. Использование ЭВМ для регистрации и обработки
экспериментальной информации и для представления результатов дает студентам
возможность почувствовать на практике те новые возможности, которые
представляет автоматизация научных исследований.

§1. Введение
Лазерная доплеровская анемометрия (от греческого "anemos" - ветер,
движение) - метод бесконтактного измерения скорости движения газообразных,
жидких и твёрдых сред, содержащих светорассеивающие неоднородности. Этот
метод в настоящее время широко используется в научных исследованиях и
технических приложениях. Лазерные доплеровские анемометры (ЛДА) позволяют
решать большой круг задач: от исследования медленных направленных движенний
в капиллярах и живых клеткак до дистанционных измерений турбулентной
скорости потоков газа в сверхзнуковых трубах и скорости ветра в атмосфере.
Величины измеряемых скоростей лежат в диапазоне от нескольких мкм/с до
нескольких км/с.
В данной работе разбираются основные принципы лазерной доплеровской
анемометрии на примере так называемой дифференциаьной схемы ЛДА,
получившей в настоящее время наибольшее распространение. Возможности этой
схемы демонстрируются на задаче измерения различных проекций линейной
скорости равномерно врашающегося шероховатого стеклянного диска, как
функции координаты точки измерения. В качестве метода обработки
доплеровского сигнала используется цифровой спектральный анализ. При этом
анемометр сопряжён с персональным компьютером.
Подробно теория доплеровского метода измерения скорости, а также
анализ других оптических схем и систем обработки сигнала ЛДА изложены в [1-
2].

§2. Общие сведения о лазерной доплеровской анемометрии
§2.1. Принцип работы ЛДА

Принцип работы ЛДА основан на эффекте Доплера, суть которого
заключается в зависимости частоты излучения света, отражённого или
рассеянного движущимся объектом, от скорости его движения.
Если на объект, движущийся со скоростью U, направить монохроматическое
излучение с частотой [pic] (рис. 1а), то отражённое (рассеянное) излучение,
вследствие эффекта Доплера, будет иметь частоту [pic], которая отличается
от частоты падающего излучения.
Доплеровский сдвиг частоты (д.с.ч.) определяется выражением [1, 2]:

[pic],
где [pic] - разностный волновой вектор, [pic] волновой вектор
падающего излучения, [pic] - волновой вектор рассеянного излучения.
С учётом введённых на рис.1 обозначений
[pic], (1)
где [pic]- показатель преломления среды, [pic] - длина волны падающего
излучения в вакууме.
Характерной особенностью выражения (1) является линейная зависимость
д.с.ч. от величины скорости, причём коэффициент пропорциональности зависит
от геометрии схемы (углов [pic] и [pic]) и оптических свойств среды.
Рассмотрим теперь схемy с двумя зондирующими пучками, имеющими
волновые векторы [pic] и [pic] и близкие частоты [pic] и [pic] (рис. 16).
Разность д.с.ч. волн рассеянных в одном направлении двумя пучками
определяется соотношением [1,2]:

[pic], (2)
где [pic]- разностная частота, которая по причине, поясняемой ниже,
часто называется частотой модуляции, [pic] - вектор чувствительности,
модуль которого равен

[pic]. (3)
Здесь [pic], [pic] - показатель преломления среды, [pic] - угол между
пучками.

[pic]
Отметим, что разность частот (2) не зависит от направления наблюдения
(покажите это !).
Схемы с двумя зондирующими пучками называются дифференциальными. При
[pic]и с учетом (2), (3), а также обозначений, введённых на рис. 1б, имеем:

[pic]. (4)
Формулы (1)-(3) позволяют рассчитывать д.с.ч. при известной скорости
объекта и заданном направлении падающих и рассеянных пучков. Из них же
можно найти и величину проекции вектора скорости [pic]на выбранное
направление, определяемое векторами [pic]и [pic] при известной величине
[pic]. Принцип работы дифференциальной схемы ЛДА можно объяснить и не
используя явным образом эффект Доплера. Суть этого объяснения заключается в
следующем: движущаяся частица рассеивает излучение в интерференционном
поле, которое образуется в области пересечения двух когерентных пучков. Для
того чтобы показать эквивалентность обоих подходов рассмотрим картину
интерференции. Для простоты будем считать, что пересекаются две плоские
монохроматические волны (рис. 2), которые запишем в следующем виде:
[pic]
[pic]
где [pic], и [pic]- амплитуды падающих волн.
Усреднённая за промежуток времени ? [pic]интенсивность света будет
равна
[pic].
Таким образом, в области пересечения двух когерентных монохроматических
пучков с разными частотами образуется периодически меняющееся во времени
распределение интенсивности, эквивалентное бегущей световой решётке.
Направление движения решётки и ее период определяются так называемым
«вектором чувствительности» [pic], а скорость движения светлых и тёмных
полос определяется разностью частот интерферирующих волн [pic].

[pic]

Если [pic], то есть частоты интерферирующих волн совпадают, то имеет
место стационарная интерференционная картина. Видность этой
интерференционной картины [pic] равна:
[pic]. (5)
При равенстве интенсивностей интерферирующих волн видность максимальна и
равна 1.
Период изменения интерференционных максимумов и минимумов
интенсивности в плоскости, перпендикулярной биссектрисе угла [pic] равен:
[pic] (6)
Чем меньше угол между падающими пучками, тем больше период
интерференционной картины. Так, при [pic] и [pic] мкм период [pic]= 0,63
мкм, а при [pic] [pic]= 109 мкм.
Пусть в области пересечения пучков движется малая и в простейшем
случае сферическая (радиуса r< [pic]) частица со скоростью,
характеризующейся вектором [pic]. Мощность рассеянного этой частицей
излучения зависит от её местоположения, размера и оптических характеристик.
Если частица находится в центре светлой интерференционной полосы, она
рассеивает максимальную мощность, а если в центре тёмной полосы -
минимальную. Особенности работы ЛДА с большими частицами (r [pic])
рассмотрены в работе [2].
Мощность, рассеянную частицами движущимися со скоростью U через
измерительный объем, можно записать в следующем виде:
[pic], (7)
где [pic], [pic] - медленно меняющаяся (по сравнению с [pic]) случайная
компонента. Коэффициент [pic] обычно принимает значения от 1 до 10 -2 в
зависимости от условий согласования амплитуд, фаз и состояния поляризации
рассеянных волн в пределах приёмной апертуры, а также от величины отношения
[pic].
Из уравнения (7) видно, что усредненная по времени [pic] рассеянная
мощность меняется с частотой [pic] равной:

[pic], (8)
и совпадающей с д.с.ч. в выражении (2). В случае [pic] (частоты
зондирующих пучков равны) регистрация д.с.ч, не даёт информации о знаке
скорости (8).
Рассмотрим случай бегущей интерференционной картины [pic]. Если частица
движется в ту же сторону, что и интерференционные полосы, то [pic]; если
частица движется навстречу полосам, то [pic]. Таким образом, ЛДА
чувствителен не только к величине проекции скорости частицы на вектор
чувствительности, но и к знаку проекции скорости.

§2.1. Блок-схема ЛДА. Характеристика основных элементов
Рассмотрим в общем виде основные функциональные блоки ЛДА, которые
позволяют представить в целом работу измерителя скорости, его возможности и
требования, предъявляемые к различным его элементам.
[pic]
Блок-схема ЛДА, построенного по дифференциальной схеме, изображена на
рис. 3. Источннком когерентного излучения является лазер, как правило
непрерывного действия, хотя в некоторых случаях возможно применение и
импульсных лазеров. Излучение лазера в блоке (2) делится на два пучка, один
из которых в знакочувствительных ЛДА пропускается через однополосный
модулятор (3), сдвигающий частоту [pic] излучения на величину [pic]. Для
сдвига частоты часто используется акустооптический эффект. Один из пучков в
некоторых измерителях проходит через линию задержки (4) для выравнивання
длины оптического пути. Блок передаюшей оптики (5) направляет зондирующие
пучки на исследуемый объект (6), содержащий движущиеся светорассеивающие
частицы. Для транспортировки излучения к измерительному объекту часто
используют оптические волноводы. Для увеличения светосилы метода рассеянный
свет собирается приёмной оптикой (7) и направляется на блок выделения
д.с.ч. (8). Выделенный сигнал, содержащий информацию о скорости движения
частиц в области пересечения зондирующих пучков, обрабатывается в блоке 9.
Дадим краткую характеристику основных элементов ЛДА и рассмотрим
требования, предъявляемые к ним.
1) Лазер. Основными параметрами, характеризующими лазер, как источник
излучения в ЛДА, являются мощность, длина волны излучения [pic], спектp и
модовый состав излучения. Выбор конкретного типа лазера обуславливается
задачей, для которой предназначен ЛДА. Например, СО2 лазеры обладают
большой мощностью (до нескольких десятков Вт в режиме непрерывноro
излучения), длина волны излучения которых [pic]=10,6 мкм лежит в
инфракрасной области спектра («окне прозрачности» атмосферы). Такое
излучение сильно поглощается жидкостями, что исключает возможность их
применения для исследования гидродинамических потоков. С другой стороны,
большая мощность CO2 лазеров позволяет создавать ЛДА, которые предназначены
для исследования воздушных потоков, содержащих аэрозоли, на расстояниях до
нескольких километров.
Длина волны [pic] определяет минимальный размер области пересечения
зондирующих пучков, то есть величину измерительного объёма. Даже если пучки
фокусируются и пересекаются в перетяжках, то поперечный размер перетяжки
ограничен дифракционным пределом, который тем меньше, чем меньше [pic].
Пространственная локальность измерений особенно важна при решении таких
задач, как исследования микроструктуры турбулентных потоков, измерение
скорости течения в микрокапиллярах и т.п. Важным параметром лазерного
излучения является его временная когерентность. Для получения хороших
характеристик ЛДА, как по пространственному разрешению, так и по отношению
сигнал/шум, излучение лазера должно иметь самый низкий поперечный тип
колебаний ТЕМ00. При измерении малых скоростей лазеры могут работать на
нескольких типах колебаний, когда д.с.ч. меньше частотного игнтервала между
этими типами колебаний. При измерении больших скоростей, как правило, нужны
одночастотные лазеры.
Мощность лазеров определяет чувствительность установки. В случае
большой концентрации рассеиваюших частиц, что характерно, например, для
двухфазных потоков и твердых тел с шероховатыми nоверхностями, возможно
применение лазеров небольшой мощности (около единиц мВт). Как правило, в
этом случае используются получившие широкое распространение He-Ne или
полупроводниковые лазеры. Если ЛДА применяется для исследования структуры
потоков газа с малой концентрацией светорассеивающих частиц, то здесь
необходимо использовать лазер большой мощности, например, аргоновый или
СО2.
2) Блок формирования зондирующих пучков может быть реализован
различными способами. Исходный пучок можно разделить на два «по амплитуде»
или «по фронту». В первом случае, как правило, используется делительный
кубик, пластинка или полупрозрачное зеркало. Во втором случае используются
призмы или маски той или иной формы.
3) В качестве однополосного модулятора, осуществляющего сдвиг частоты
излучения в одном из двух зондирующих пучков (в знакочувствительных ЛДА),
чаще всего используется акусто- или электрооптические преобразователи и
вращающиеся дифракционные решётки.
4) Линия задержки предназначена для выравнивания длин оптических путей
зондирующих пучков. Допустимая разность хода определяется длиной
когерентности 1k излучения лазера.
5) Блок передающей оптики осуществляет сведение зондирующих пучков в
область, где производятся измерения. В тех случаях, когда важна локальность
измерения, пучки фокусируются. В простейшем случае этот блок представляет
собой фокусирующую линзу или объектив. В некоторых случаях линза или
объектив ставится в каждый пучок, а сведение пучкoв в область измерения
осуществляется поворотными зеркалами.
7) Блок приёмной оптики, как правило, представляет собой линзу или
объектив, иногда с переменной апертурой. В случаях, когда д.с.ч. выделяется
из света, расееянного назад (например, при измерении скорости движения
оптически непрозрачного твёрдого тела), функции приёмной оптики может
выполнять передающая линза. Задача блока приемной оптики состоит в том,
чтобы собрать на фотоприемник рассеянное излучение из измерительного объема
в возможно большем телесном угле.
8) Функции блока выделения д.с.ч. выполняют фотоприёмники (ФП):
фотодиоды (ФД) или фотоумножители (ФЭУ), обладающие большой
чувствительностью и высоким быстродействием. Фэу обычно применяются при
измерениях скорости движения слабо рассеивающих сред, когда регистрация
рассеянного излучения ведётся в режиме счета фотонов. В дифференциальнык
схемах ЛДА ФП осуществляет выделение д.с.ч. путём прямого
фотодетектирования колебаний интенсивности рассеянного света,преобразуя их
в колебания силы фототока. Частотно модулированный выходной ток ФП и
является выходным сигналом ЛДА.
9) устройство обработки сигнала преобразует информацию, содержашуюся в
фототоке, в удобную для пользователя форму. При этом используются либо
аналоговые, либо цифровые методы, либо их комбинация. Более подробно
охарактеризуем эти методы после того, как рассмотрим свойства сигнала ЛДА.

§2.2. Свойства сигнала ЛДА
Если часть света, рассеянногo некоторой частицей при пересечении ею
интерференционного поля, попадает в блок приёмной оптики и далее на
фотоприёмник, то на выходе фотоприёмника появляется импульс фототока,
промодулированный частотой [pic] (см. рис. 4).
Длительность импульса [pic] определяется временем нахождения частицы в
измерительном объеме - области пересечения зондирующих пучков: [pic], где
[pic] - размер этой области в направлении движения частицы (см. рис.2),
[pic] - скорость частицы. Глубина модуляции определяется видностью
интерференционной картины (5). Очевидно, что спектральная плотность
мощности [pic] этого импульса (рис. 5) содержит низкочастатную
составляющую, а также собственно доплеровскую составляющую, сдвинутую
относительно нулевой частоты на [pic] или [pic] (в зависимости от
направления скорости U частицы). Низкочастотная составляющая спектра
обусловлена как постоянной составляющей сигнала с фотоприемника, так и
фликкер-шумом источника излучения и фотоприемника. Ширина доплеровской
составляющей спектра [pic], а, следовательно, и неопределённость в
нахождении частоты и проекции скорости тем больше, чем меньше измерительный
объём, т.е. чем выше локальность измерения. Если интерференционное поле
пересекает поток случайно расположенных частиц различных размеров, то
сигнал на выходе фотоприёмника представляет собой наложение радиоимпульсов,
приходящих в случайные моменты времени и имеющих разную начальную фазу,
длительность и амплитуду, поскольку рассеивающие свойства частиц также
могут силъно отличаться. Реальный сигнал и его выборочный спектр
изображены на рисунках 6 и 7 соответственно.


[pic]

[pic]Рис. 7. Выборочный спектр реального сигнала.

Для того чтобы из такого квазистохастического сигнала извлечь
информацию о скорости движения исследуемогo объекта, разработаны различные
методы его обработки. Самым простым из них является измерение несущей
частоты сигнала с помощью частотомера. Спектральный анализ позволяет более
точно выделить доплеровскую составляющую в энергетическом спектре сигнала с
фотоприемника. Чтобы повысить точность этих измерений, обычно,
предварительно сигнал пропускают через полосовой фильтр, отрезая
низкочастотную компоненту с фликкер-шумом и высокочастотные шумы.
Существуют и более сложные специальные методы, осуществляющие, например,
аналоговую иди цифровую демодуляцию сигнала. В этом случае говорят о
применении специализированных процессоров доплеровского сигнала.
Другой альтернативой является компьютерная обработка сигнала с
ислользованием современных алгоритмов, например, быстрого преобразования
Фурье (БПФ) и метода периодограмм. Использование ЭВМ существенно расширяет
возможности проведения экспериментов с помощью ЛДА. Заметим, однако, что
этот метод используется, как правило, лишь в случае измерения сравнительно
невысоких скоростей.

§2.2. Описание установки и программного интерфейса
В настоящей работе для измерения скорости вращения матового диска
используется ЛДА, сопряжённый с персональным компьютером и построенный по
дифференциальной схеме без использования модулятора и линии задержки. Для
обработки сигнала применяется спектральный анализ с использованием цифровых
фильтров и различных методов усреднения спектров.
Блок-схема установки, на которой выполняется задача, изображёна на
рис. 8. Излучение одномодового He-Ne (типа ЛГН-207А) с длиной волны
[pic]=0,53 мкм или полупроводникового лазера ([pic]=0,56 мкм) (1) и
мощностью около 2 мВт направляется на телескопический расширитель пучка
(2). Телескоп расширяет лазерный пучок от исходного диаметра 1 мм до 22 мм
(по уровню интенсивности е-1). На телескоп надевается вращающаяся вдоль
горизонтальной оси маска (3) с двумя симметрично расположенными
отверстиями, вырезающими из фронта расширенного пучка два узких
параллельных пучка. Линза (4) фокусирует пучки и сводит их на поверхности
светорассеивающего диска (5). 3-х координатный столик с микрометрическими
подачами, на котором крепится диск вместе с вращающим его мотором и
редуктором, обеспечивает перемещение диска в трех взаимно ортогональных
плоскостях [pic]. Ось вращения диска параллельна оптической оси системы ,
включающей элементы (1 - 4). Пусть система координат, приведенная на рис. 8
связана с центром диска, ось [pic]направлена вдоль оптической оси, а ось
[pic] - вертикально вверх. Пусть также угол поворота маски отсчитывается от
оси [pic].
[pic]
Рассеянное диском излучение собирается линзой (6) на фотоприёмник (7)
(ФД-11А). Выходной электрический сигнал с ФД поступает на усилитель (8),
собранный на микросхемах 140УД8, с коэффициентом усиления порядка 1000 и
полосой пропускания от 200 Гц до 8 КГц на уровне 0,7. Сигнал с выхода
усилителя подаётся параллельно на контрольный осциллограф (9) и линейный
вход звуковой карты, расположенной в компьютере (10). Звуковая карта
настроена на максимальную из возможных частот оцифровки 44100 Гц и 16-
разрядный режим работы. Коэффициент усиления входного тракта
устанавливается в окне настройки звуковой карты (см. ниже).
Конструктивно элементы (6 - 8) укреплены на поворотном устройстве с
вертикальной осью вращения (О), которая обозначена на рисунке 8. Это
устройство позволяет исследовать зависимость доплеровского сигнала от угла
рассеяния. Выведение измерительного объема на ось вращения осуществляется
настройкой телескопа. Рассеивающая поверхность стеклянного диска выводится
на ось вращения путем перемещения диска вдоль оси [pic] с помощью
координатного столика. Перемещением диска вдоль осей [pic] и [pic] можно
менять величину и направление скорости рассеивающих частиц, пересекающих
измерительный объем.
Программный интерфейс к задаче реализован в среде LabView и
предоставляет следующие возможности:
1. Получение выборки с помощью звуковой карты и контроль оцифрованного
сигнала и его спектральной плотности мощности (СПМ) с помощью
соответствующих графиков (пункт "Получение сигнала").
2. Пункт "Получение сигнала" включает также обработку сигнала и
спектра: цифровую фильтрацию и усреднение выборочных спектров оцифрованного
сигнала различными методами.
3. Снятие зависимости проекции скорости рассеивающих частиц от
следующих параметров: ориентации маски (3), а также координат [pic] и [pic]
столика (пункт "Зависимость от параметра").
4. Поскольку анализируемый сигнал является случайным, то погрешность
оценки проекции скорости рассеивающих частиц из анализа выборочных спектров
зависит от числа их усреднений. Эта зависимость и исследуется в этом пункте
("Оценка точности измерений в зависимости от числа усреднений").
Для получения оценки проекции скорости на направление, определяемое
положением маски (3), используется формула (4). Доплеровская частота
находится из анализа выборочных спектров сигнала, регистрируемого
фотоприемником. Сигнал перед спектральной обработкой может быть подвергнут
цифровой фильтрации с выбором типа фильтра и его полосы пропускания.
Выборочные спектры вычисляются с использованием алгоритма БПФ.
Алгоритм определения значения доплеровской частоты следующий. Сначала
в усредненном спектре ищется гармоника максимальной амплитуды, которая, как
правило, принадлежит доплеровской компоненте. Затем относительно нее по
уровню 0,5 вычисляется эффективная ширина доплеровского контура [pic] (см.
рис.5). Далее в интервале частот [pic] находится частота [pic],
соответствующая положению центра тяжести той части доплеровского спектра,
которая ограничена частотами [pic] и [pic]. Найденное значение [pic] и
принимается равным доплеровской частоте [pic]. В общем случае [pic], хотя
эти значения могут быть близки. Изложенный алгоритм определения
доплеровской частоты эффективен для усредненного спектра, если доплеровский
контур достаточно хорошо выделяется на фоне шумовых спектральных компонент,
при этом алгоритм оказывается слабо чувствительным к случайным выбросам
спектральной плотности. Заметим, что чувствительность ЛДА, т.е. минимально
возможная проекция скорости, которая может быть измерена, определяется
фликкер-шумом, поскольку при уменьшении проекции скорости доплеровский
контур оказывается «погруженным» в спектр фликкер-шума (см. Рис. 5) и
алгоритм будет «цепляться» не за доплеровский контур, а случайную
спектральную компоненту, принимающую максимальное значение в выборочном
спектре. Кроме того, алгоритм оказывается чувствительным к импульсным
помехам в сигнале, которые, например, могут иметь сейсмическое
происхождение (удар по столу, на которой расположена установка, резкое
закрывание двери в помещении и пр.). Импульсная помеха приводит к появлению
в высокочастотной части выборочного спектра случайных интенсивных
спектральных компонент, превышающих спектральные компоненты доплеровского
контура. При обработке такого выборочного спектра описанный выше алгоритм
будет давать значение скорости не соответствующее доплеровской.
Программа построена по принципу меню (см. рис. 9) и состоит из трех
подпрограмм: «получение сигнала», «зависимость от параметра» и "точность
измерений в зависимости от числа усреднений".
[pic]
Рис. 9. Лицевая панель программы «Меню»

Принцип работы с этими подпрограммами один и тот же - по команде
пользователя со звуковой карты оцифровывается сигнал заданной длины («Длина
выборки»), а затем подбираются оптимальные параметры цифровой фильтрации и
обработки его спектра. В программе «получение сигнала» цифровая фильтрация
и обработка спектра проводятся независимо. В остальных двух программах
исходная выборка предварительно пропускается через цифровой фильтр, а затем
ее спектр обрабатывается одним из методов. После этого фиксируются значения
амплитуды максимальной гармоники («max амплитуда»), положение центра
тяжести доплеровского контура (доплеровская частота) и проекция скорости
рассеивающих частиц на направление маски («скорость диска»).
В программе «получение сигнала» (рис. 10) демонстрируются возможности
обработки спектра и влияние обработки на результат измерения центра тяжести
доплеровского контура и оценку проекции скорости. Регистрация выборки
заданной длины производится по нажатию клавиши «получить сигнал», процедуры
фильтрации сигнала и усреднения его спектра производятся независимо. По
нажатию клавиши «выход» происходит возврат в программу «Меню».
Наверху отображаются оцифрованный доплеровский сигнал и под ним его
спектр. Двух ползунковый регулятор служит для изменения пределов по оси
частот графика спектра. Слева от графика «Спектр сигнала» выводятся
результаты обработки спектра: максимальная амплитуда гармоники,
доплеровская частота и проекция скорости.
Полученная выборка может быть подвергнута цифровой фильтрации с
выбранным типом фильтра и заданием полосы пропускания. Верхняя и нижняя
частоты пропускания фильтра указываются положением курсоров на графике
«Спектр сигнала» (с помощью кнопки, расположенной справа под графиком
«Спектр сигнала» курсоры приводятся на середину графика). Спектр сигнала
после фильтрации отображается на нижнем левом графике. Изменяя тип фильтра
и полосу его пропускания можно проследить за изменением спектра и оценок
доплеровской частоты и скорости. Слева от графика выводятся результаты
обработки спектра сигнала, подвергнутого цифровой фильтрации: максимальная
амплитуда гармоники, доплеровская частота и проекция скорости.
Обработка спектров проводится для уменьшения дисперсии и представляет
собой усреднение спектра выборки одним из двух методов: периодограмм и
скользящего среднего. Метод периодограмм состоит в том, что исходный сигнал
(выборка из [pic] отсчетов) разбивается на [pic] секций и для каждой i-ой
секции вычисляется выборочный спектр [pic]. Результирующий спектр [pic]
получается усреднением выборочных спектров
[pic].
Число секций, на которые будет разбиваться исходный сигнал, устанавливается
с помощью движкового регулятора, при этом длина секции высвечивается в
соответствующем окошке. Отметим, что число секций влияет на дисперсию [pic]
и спектральное разрешение [pic] результирующего спектра [pic].
В методе скользящего среднего результирующий спектр получается путем
свертки спектра выборки с прямоугольным спектральным окном, ширина которого
составляет m гармоник. Ширина спектрального окна также влияет на дисперсию
[pic] и спектральное разрешение [pic] полученного спектра.
[pic]

Программа «зависимость от параметра» (рис. 11).
Конкретный смысл параметра, от которого ищется зависимость
доплеровской частоты, может быть разным и зависит от упражнения. Этим
параметром может быть значение координат Х или Y точки измерения
относительно центра вращающеroся диска или значение угла поворота маски
[pic], определяющего положение вектора чувствительности анемометра
относительно вектора линейной скорости движения рассеивающих частиц в точке
измерения.
Для снятия той или иной зависимости необходимо выполнить сдедующие
действия:
1) Предварительно в программе «Получение сигнала» для данных условий
измерений следует оптимальным образом выбрать: сначала параметры цифрового
фильтра (тип и полосу пропускания) а затем метод усреднения спектра и его
параметры, а также параметры методов усреднения (метода периодограмм и
скользящего среднего).
2) Войти в программу "Зависимость от параметра"; при этом производится
выборка и ее цифровая фильтрация, а затем спектр фильтрованной выборки
усредняется одним из методов (см. соответствующие закладки окна - «фильтр /
усреднение» на рис.11).
3) В окне «параметр» ввести значения параметра, и после нажатия
клавиши "Выборка" получить очередную дискретную выборку сигнала и ее
усредненный спектр, а также амплитуду максимальной гармоники и скорость;
4) Полученное значение скорости вместе со значением параметра можно:
вывести на график (клавиша «на график»), запомнить в файл (клавиша
«запомнить»). Клавиша «почистить» очищает график.
5) После этого следует перейти к пункту 3 для измерения скорости при
другом значении параметра.

[pic]

Программа «точность измерений в зависимости от числа усреднений» (рис. 12,
13)
Доплеровский сигнал, полученный от рассеяния на вращающемся
шероховатом диске, является, как отмечалось выше, случайным по своей
природе. В нем также отражается временная нестабильность интенсивности
излучения лазера и шумы входного тракта (фотодиода и усилителя). Как
следствие, выборочный спектр также является случайным, а среднеквадратичное
отклонение или дисперсия являются мерой этой случайности. Поэтому
полученная, описанным выше способом, оценка проекции скорости также
является случайной величиной со своей дисперсией. Используемые в задаче
методы усреднения выборочного спектра (метод периодограмм и скользящего
среднего по спектру) позволяют уменьшить дисперсию спектральной оценки и,
как следствие, уменьшить дисперсию оценки проекции скорости. В методе
периодограмм параметром, определяющим дисперсию, является число секций
[pic], на которые разбивается выборка. В методе скользящего среднего
дисперсия зависит от ширины скользящего спектрального окна, которое можно
характеризовать числом гармоник [pic]. Данная программ предоставляет
возможность исследования зависимости дисперсии скорости от параметров [pic]
и [pic].
В программе используется два режима: настройка и измерение, выбор
которых производится переключателем «настроить/померить» (см. Рис. 12). В
режиме настройка (активна верхняя часть лицевой панели на Рис. 12)
производится выбор типа фильтра и его полосы пропускания (закладка
«фильтр»), а также выбор метода усреднения спектра и параметры усреднения
(закладка «усреднение»). После нажатия клавиши «Сигнал», расположенной в
закладке «усреднение», производится выборка и ее обработка. Результаты
настройки отражаются на соответствующих графиках: спектры исходной и
фильтрованной выборки (закладка «фильтр») и усредненный спектр (закладка
«усреднение»).
В режиме измерения (активна нижняя часть лицевой панели на рис. 13)
сначала надо задать число выборок («число измерений») [pic], по которым
будет вычисляться среднее значение и дисперсия [pic] оценки проекции
скорости. После нажатия клавиши «Пуск» компьютер выполнит набор заданного
числа выборок и их обработку. Результаты обработки отображаются в
соответствующих окнах («среднее» и «дисперсия»). На левый график выводятся
значения скорости, полученные при обработке отдельных выборок. Значение
дисперсии скорости [pic] и число усреднений [pic] можно сохранить в файле
нажатием клавиши «сохранить». На правый график, отражающий зависимость
[pic] ([pic] в методе периодограмм и [pic] в методе скользящего среднего),
можно отобразить (нажатием клавиши «на график») результат обработки для
заданного числа измерений [pic]. На этот же график выводится доверительный
интервал для измеренного значения.
Из-за импульсных помех отдельные выборки могут быть «плохими», т.е.
значения, полученные при их обработке, могут очень сильно отличаться от
истинного значения доплеровской скорости. Очевидно, что результат обработки
«плохой» выборки должны быть исключен из процедуры оценки дисперсии. Такая
возможность предусмотрена в программе. «Плохая» выборка может быть легко
идентифицирована на левом графике, на котором отображены значения скорости
для последовательности выборок, и исключена из обработки для получения
оценки дисперсии. Для этого результат обработки «плохой» выборки следует
выделить курсором и нажать клавишу «выкинуть точку». При этом выделенное
значение будет исключено из процедуры оценки дисперсии и заново пересчитаны
значения средней скорости и дисперсии.


[pic]
Рис. 12. Лицевая панель программы «точность измерений в зависимости от
числа усреднений» в режиме настройки


[pic]

Рис. 13. Лицевая панель программы «точность измерений в зависимости от
числа усреднений» в режиме измерение

§3. Экспериментальная часть

Упражнение 1. Юстировка оптической схемы ЛДА. Получение доплеровского
сигнала.
Ознакомившись с расположением основных блоков установки в соответствии
со схемой, изображённой на рис. 8, включите тумблеры «сеть» блоков питания
лазера, усилителя, осциллографа и компьютера.
ВНИМАНИЕ! Выполняя задачу, соблюдайте правила техники безопасности.
Следите за тем, чтобы лазерное излучение не попало в глаза Вам и Вашим
коллегам. Поверхность рассеивающего диска руками не трогать.
Помещая рассеивающий экран, в качестве которого можно использовать
лист белой бумаги, между маской (3) и линзой (4) (см. рис. 8), убедитесь,
что пучки после маски имеют одинаковую интенсивность при любом угле
поворота маски вокруг оси телескопа. Если интенсивности пучков различны, то
с помощью преподавателя подъюстируйте телескоп. Убедитесь, что после линзы
(4) пучки пересекаются на поверхности матового диска. Если матовая
поверхность диска не проходит через область пересечения пучков, то с
помощью микрометрической подачи сдвиньте диск вдоль оси Z в нужное
положение. С помощью листа белой бумаги пронаблюдайте спекл картину и
диаграмму рассеяния за матовым диском. Включите мотор тумблером,
установленным на кожухе мотора, и убедитесь, что диск равномерно вращается.
Поворотное устройство, на котором укреплён блок приёмной оптики и
фотоприёмник, установите в положение, соответствующее регистрации света,
рассеянного вперёд. Маску установите так, чтобы лучи, вышедшие из нее
лежали в горизонтальной плоскости. Убедитесь в том, что изображение
измерительного объёма попадает на ФП. На экране осциллографа должна
появиться временная развёртка доплеровского сигнала, примерный вид которого
показан на рис. 6.
Если на экране осциллографа развёртка луча есть, но доплеровского
сигнала нет, проверьте, не реализуется ли, по крайней мере, одна из
следующих возможностей:
- интенсивности пучков сильно отличаются;
- зондирующие пучки не пересекаются на поверхности диска (подобрать
положение диска вдоль оси Z);
- изображение измерительного объёма не попадает на фоточувствительную
поверхность ФД;
- диск не вращается;
- неправильно установлены чувствительность и/или длительность
развёртки осциллографа.
Перемещением диска вдоль горизонтальной оси получите доплеровский
сигнал максимальной величины. Перед запуском задачи на компьютере (иконка
расположена на рабочем столе) следует убедиться в правильности настройки
звуковой карты. Откройте окно настройки звуковой карты (двойной «клик» на
иконку динамика в нижнем правом углу рабочего стола) и убедитесь, что
активны только две опции: Volume control и Line in (остальные должны
находиться в состоянии Mute), причем в опции Line in движок должен быть
установлен на середину шкалы. Запустите задачу и перейдите в пункт
«получение сигнала». Пронаблюдайте доплеровский сигнал и его спектр.
Вертикальная ось на графике «Доплеровский сигнал» (см. рис. 10)
откалибрована в единицах АЦП звуковой карты, динамический диапазон которого
составляет -32700 - + 32700. Объясните, чем определяется ширина
доплеровского спектра и от чего зависит дисперсия его оценки.
Проанализируйте изменение спектра выборки при использовании фильтров
различных типов, меняя их полосу пропускания. Используя методы усреднения
спектра, проанализируйте влияние параметров этих методов (числа секций и
ширины спектрального окна) на усредненный спектр и значение оценки проекции
скорости. На основе этого анализа выберите наилучшие на Ваш взгляд тип
фильтра и значения параметров обработки спектра. Объясните свой выбор.
Упражнения 2 и 4 следует проводить с выбранными Вами типом фильтра и
параметрами обработки.

Упражнение 2. Анализ зависимости параметров доплеровского сигнала от
качетва юстировки и параметров оптической схемы ЛДА.
В этом упражнении используется программа «получение сигнала».Выводя
доплеровские спектры, компьютер указывает максимальное значение
спектральной плотности, на которое осуществляется нормировка. Это значение
непосредственно связано с величиной отношения сигнал/шум, т.е. с качеством
сигнала. Отсюда следует возможность анализа качества сигнала как по выборке
(или осциллографу), так и по выборочному спектру.
Получите ряд доплеровских спектров, соответствующих разному качеству
юстировки системы, а также при различных положениях маски и фотоприемника.
Проанализируйте полученные результаты.
Установите оптимальную юстировку. Перемещая диск с помощью
микрометрических винтов в плоскости [pic], добейтесь, чтобы область
пересечения пучков находилась на максимальном расстоянии от центра вращения
и располагалась на вертикальный диаметр диска [pic]. Установите положение
маску (3) так, чтобы отверстия в ней располагались по горизонтали. При
этом, как легко убедиться, измеряемая проекция скорости принимает
максимальное значение. Снимите зависимость интенсивности максимальной
гармoники и проекции скорости от углового положения фотоприемника.
Объясните полученные результаты.

Упражнение 3. Оценка зависимости точности измерения скорости от числа
усреднений
Цель упражнения состоит в исследовании зависимости дисперсии
[pic]оценки скорости от усреднения выборочного спектра доплеровского
сигнала. В методе периодограмм дисперсия оценки скорости определяется
числом секций [pic], на которые разбивается исходная выборка, а в методе
скользящего среднего - шириной спектрального окна, выраженного в числе
гармоник [pic], с которым выполняется процедура скользящего среднего; при
этом [pic].
Для получения экспериментальной оценки дисперсии скорости следует
получить серию из [pic] однотипных опытов (получение выборки, ее обработка
и вычисление скорости), по результатам обработки которых может быть
вычислена средняя скорость и дисперсия скорости. Эти возможности
предоставляет программа "точность измерений в зависимости от числа
усреднений".
Установите положение маски (3) и диска (5) в такое положение, чтобы
измеряемая проекция скорости принимала максимальное значение. Установите
фотоприемник на регистрацию рассеяния вперед (при этом отношение сигнал/шум
будет максимальным). Для серии [pic] и ряда значений параметра [pic]
измерить среднюю скорость и дисперсию скорости [pic]для обоих методов
усреднения (метода периодограмм [pic] и скользящего среднего [pic]).
Результаты представить в виде графиков с координатами: [pic] - ось ординат,
[pic] - ось абсцисс, на которых должен быть отражен доверительный интервал,
определяемый значением [pic]. Повторите измерения, увеличив угол
регистрации рассеянного излучения.

Упражнение 4. Снятие зависимостей проекции скорости от положения
измерительного объема на диске и положения маски.
При подготовке к задаче получите зависимость [pic] проекции скорости
[pic] точки нa диске с координатами [pic] на направление ([pic]),
соединяющее центры отверстий в маске для заданной скорости вращения диска
([pic]об/сек) (см. Рис. 14). Постройте графики зависимостей [pic].
Цель упражнения состоит в последовательном получении экспериментальный
зависимостей [pic] с использованием пункта меню «зависимость от параметра».
Изменять координаты при снятии зависимостей следует так, чтобы на графике
получилось 10-15 точек. Пo полученным зависимостям следует определить
чувствительность метода - минимальное значение проекции скорости, которое
может быть измерено и объяснить чем определяется эта чувствительность.
Поставьте личный рекорд измерения минимальной скорости.

[pic]

Упражнение 5. Однолучевая схема.
Классическая однолучевая схема ЛДА при своей практической реализации
требует опорного пучка для выделения доплеровской частоты фотоприемником.
Поскольку часть света проходит через матовый диск без рассеяния, то в
данной установке можно реализовать однолучевую схему, если закрыть одно из
отверстий в маске. Убедитесь, что наблюдается доплеровский спектр, причем
его положение не зависит от положения матового диска относительно линзы (4)
(см. рис. 8). Пересечение падающего пучка с диском выведите на вертикальный
диаметр диска. Снимите зависимость доплеровской скорости от координаты
вертикального положения пучка [pic]относительно диска. На экспериментальную
зависимость нанесите теоретическую

Литература
1. Б.С. Ринкевичус, Лазерная анемометрия, М., «Энергия», 1978
2. Ю.Д. Дубнищев, Б.С.Ринкевичус, Методы лазерной доплеровской
анемометрии, М., «Наука»
3. F. Durst, A. Melling, J.H. Whitelaw Principles and practice of laser-
doppler anemometry, Academic Press, 1976

-----------------------
[pic]

[pic]

[pic]

[pic]

[pic]

О

О

[pic]


9


8


7


6




5

3

4




2

1

Рис. 3. Блок-схема дифференциального ЛДА с частотной модуляцией пучка. 1 -
лазер, 2 - блок формирования зондирующих пучков, 3 - однополосный
модулятор, 4 - линия задержки, 5 - блок передающей оптики, 6 - исследуемый
объект, 7 - блок приемной оптики, 8 - блок выделения д.с.ч., 9 - устройство
обработки сигнала.

б)

а)

Рис. 1. Взаимное расположение векторов рассеяния в схеме с одним (а) и
двумя (б) зондирующими пучками. [pic] - вектор скорости рассеивающей
частицы объекта. На рис.1а [pic] - частота падающего излучения, [pic] -
волновые вектора падающего и рассеянного излучений, [pic] - разностный
волновой вектор. На рис.1б [pic] - частоты пучков падающего излучения,
[pic] - волновые вектора излучения рассеянного первым и вторым пучком,
[pic] - разностный волновой вектор.

[pic]

[pic]

[pic]

Рис. 14. К выводу зависимости [pic]. [pic] - ось диска, вращающегося с
угловой скоростью [pic]об/сек; [pic] - скорость точки[pic], в которой
расположен рассеивающий объем; [pic] - след плоскости, в которой лежат
интерферирующие пучки и осью [pic]

U

?

x

y

X

Y

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

Рис. 4. Импульс фототока от одной частицы

Рис. 5. Энергетический спектр импульса фототока от одной частицы

[pic]

t

Рис. 6. Пример реализации реального сигнала ЛДА

Z

Y

X

7

6

5

4

3

2

1

8

10

9

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

k1

k2

x

U

K

Рис. 2. Интерференционная модель ЛДА (И - индикатриса рассеяния, О -
приемный объектив).

О

И

?

Т

L

[pic]

Рис. 8 Блок-схема экспериментальной установки: 1) лазер; 2) расширитель
пучка; 3) маска; 4) линза; 5) рассеивающий диск, установленный на 3-х
координатном столике; 6) линза; 7) фотодиод; 8) усилитель; 9) контрольный
осциллограф; 10) звуковая карта, встроенная в компьютер.

??

Рис. 10. Лицевая панель программы «получение сигнала».

Рис. 11. Лицевая панель программы «зависимость от параметров».


А