Способ измерения поля скоростей в газовых и конденсированных средах

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

 

Изобретение относится к измерительной технике и может найти применение в экспериментальной гидро- и аэродинамике, в промышленных технологиях, связанных с необходимостью исследования полей скоростей в потоках газовых и конденсированных сред, а также скоростей движения поверхностей.

Известны способы, предназначенные для измерения поля скоростей в газовых и конденсированных средах методами когерентной и некогерентной оптики. К ним относятся лазерные доплеровские измерения поля скоростей (в англоязычной литературе - Doppler Global Velocimetry - DGV) [Dubnishchev Yu., Chugui Yu. V., Kompenhanc Jurgen, Belousov P.P., Belousov P.Ya. Method and Device for Laser Doppler Measurement of a Velocity Field of a Fluid Flow. EP 1978369 B1, G01P 3/36. Publication 03.03.2010 Bulletin 2010/09], измерения по трекам изображений частиц (Particle Image Velocimetry - PIV) [Raffel M., Willert CE., Kompenhans J. (1998). Particle Image Velocimetry. Springer, Berlin]. Измерение поля скоростей методом DGV заключается в последовательном выполнении следующих операций: освещают сечение исследуемой среды лазерной плоскостью, конвертируют доплеровский сдвиг частоты рассеянного света в интенсивность и регистрируют изображения лазерной плоскости в конвертированном свете. Распределение интенсивности конвертированного света в изображении лазерной плоскости на фотоматрице отображает поле скоростей в исследуемом сечении.

Недостатком способа DGV является зависимость числа каналов регистрации частотно-конвертированных изображений лазерной плоскости от размерности координатно-измерительного базиса (числа проекций) вектора скорости. К другим недостаткам следует отнести требование тонкой температурной стабилизации конвертора частота-интенсивность, обычно выполняемого на основе молекулярной ячейки поглощения или интерферометра, необходимость автоматической привязки когерентно-оптической передаточной функции конвертора к частоте излучения лазера и использование только когерентного излучения. В методах PIV недостатком являются ограничения на пространственно-временное разрешение исследуемых процессов из-за проблем с идентификацией треков при высокой концентрации частиц и многократной экспозиции.

Кроме того, известен способ [Shigeo Hasokawa, Akio Tomiyama. Spatial Filter Velocimetry Based on Time-Series Particles Images. Exp. Fluid (2012). 52: p.1361-1372], являющийся прототипом предложенного изобретения, основанный на эмуляции динамической пространственной фильтрации изображений движущихся частиц, регистрируемых на фотоматрице при видеосъемке. Частицы движутся в выделенном световой плоскостью сечении исследуемой среды. При считывании изображений частицы в процессе воспроизведения фильма на программном уровне в пиксельной структуре формируется периодическая дискретизирующая последовательность с заданным пространственным периодом и фазовой скоростью, которая эмулирует динамический пространственный фильтр. Этот фильтр трансформирует скорость пространственного перемещения изображения частицы вдоль координатной оси, заданной дискретизирующей последовательностью, в частоту. Измеряя эту частоту и учитывая коэффициент линейного увеличения оптической системы, формирующей изображение, получают величину скорости движения частицы в заданной световой плоскости.

Недостатком этого способа является невозможность получения информации о пространственном поле 3D векторов скоростей в исследуемых газовых и конденсированных средах.

Задачей (техническим результатом) предлагаемого изобретения является расширение функциональных возможностей на измерение 3D поля скоростей в газовых и конденсированных средах.

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

На Фиг. 1 показана структурная схема устройства, реализующего предложенный способ.

На Фиг. 2 показаны спектры эмулированных фильтрованных сигналов на пикселах.

На Фиг. 3 показаны r, g и b - сечения световых плоскостей.

На Фиг. 4 показана корреляционная функция Ksrg (τ) фильтрованного сигнала.

Устройство (Фиг. 1) содержит: конфокальные объективы 1 и 2 с пространственным фильтром 3, расположенным в совмещенной фокальной плоскости, которые образуют телецентрическую систему, проектирующую на фотоматрицу видеокамеры 4 изображение структурированного зондирующего поля, сформированного в исследуемой среде. Видеокамера подключена к компьютеру 5. Телецентрическая система [Бычков P.M., Чугуй Ю.В. Беседы о геометрической оптике. Изд-во Сибирского отделения российской академии наук, 2011, 476 с.] обеспечивает формирование изображения структурированного светового поля на фотоматрице с одинаковым коэффициентом увеличения в широких пределах, независимо от положения на оптической оси. Зондирующее поле состоит из параллельных световых плоскостей. Длины волн λr, λg и λb соответствуют световым чувствительностям пикселей в структуре фотоматрицы. В rgb-кодировке это красный, зеленый и синий цвета.

Способ осуществляется следующим образом. Оптические неоднородности, присутствующие в среде, отслеживают ее движение и, попадая в световые плоскости, рассеивают свет. Изображения оптических неоднородностей в рассеянном свете, сформированные в телецентрической проекции с известным коэффициентом передачи, фиксируются пиксельной структурой соответствующей цветовой чувствительности на фотоматрице и записываются в последовательности кадров цветного видеофильма. Реструктуризация полей скоростей движения оптических неоднородностей в световых плоскостях выполняется путем динамической пространственной фильтрации на длинах волн λr, λg и λb при компьютерном воспроизведении фильма. Периодическая структура динамического пространственного фильтра эмулируется в процессе считывания пиксельной структуры изображений частиц с известным межкадровым интервалом воспроизводимого фильма. Считывание выполняется по строкам и столбцам пиксельной структуры кадра, которым можно сопоставить ортогональные оси x и у в декартовой системе.

Рассмотрим динамическую пространственную фильтрацию изображения частицы, сформированного в световом поле с длиной волны λr. Изображение частицы или оптической неоднородности рассматривается как пространственный сигнал, фиксируемый пиксельной структурой фотоматрицы. Дискретизирующая последовательность, эмулирующая одномерный динамический пространственный фильтр при считывании одной строки по оси x, может быть описана выражением:

где Δх - пространственный полупериод дискретизирующей последовательности, состоящий из заданного числа пикселей; V0xх0; τ0 - межкадровый интервал;

V0x - фазовая скорость по оси х эмулированного пространственного фильтра с периодом 2Δх; δ(х-V0xt-mΔх) - дельта функция Дирака;

t - время;

m - порядковый номер пространственного полупериода. Суммирование идет по числу полупериодов в дискретизирующей последовательности, которая эмулирует динамический пространственный фильтр.

Пусть одномерное сечение пространственного сигнала, изображение которого сформировано в свете с длиной волны λr (красный цвет), описывается функцией s(x-Vxt), где Vx - скорость движения изображения по оси х. Тогда для дискретизированного пространственного сигнала имеем:

Разложим дискретизирующую последовательность ηr(x) в РВД Фурье:

где

Учитывая фильтрующие свойства δ -функции, получаем из (4)

Здесь учтено: einπ=(-1)n; ei2nπ=1. Подставляя (5) в (3) получаем

Выражение для дискретизированного сигнала находится подстановкой (6) в (2):

Пренебрегая в (7) высшими гармониками (|n|>1), получаем с точностью до знака:

Здесь - пространственная частота дискретизирующей последовательности. В (8) К является волновым числом эмулированного гармонического пространственного фильтра cos[K0x(x-V0xt)], фазовая скорость которого V0x;

Vx проекция скорости движения изображения частицы по оси x.

Переходя в (8) к лабораторной системе отсчета, получаем:

где ωx=K0xVx; ω=K0xV0x.

Отсюда

Эмуляция пространственного фильтра по оси у (по столбцам матрицы) осуществляется по аналогии с описанной выше процедурой. Выполняя считывание столбца по оси у, получаем для дискретизированного сигнала с точностью до знака:

где ωy=K0yVy;

ω=K0yV0y;

Отсюда для у-проекции скорости движения изображения частицы находим:

Частоты ω и ω возникают при динамической фильтрации изображений неподвижных частиц и неоднородностей фона. Выражения (9) и (11) описывают эмуляцию динамической пространственной фильтрации изображения в плоскости х, у. Частоты ωх, ω, ωу и ω определяются, например, посредством преобразования Фурье эмулированных сигналов. Переход от частот ωх, ω, ωу и ω к скоростям Vx и Vy на выбранных пространственных масштабах выполняется по формулам (10) и (12).

Пространственная реконструкция поля векторов скоростей реализуется при формировании структурированного зондирующего поля в виде параллельных световых плоскостей на длинах волн, соответствующих цветовой чувствительности пикселей в структуре фотоматрицы (например, λr, λg и λb). Действительные величины компонент скорости движения частицы определяются как Vx/Г и Vy/T, где Г - коэффициент увеличения оптической системы, формирующей изображения. Совершенно аналогично выполняется измерение поля х- и у-компонент вектора скорости движения частиц в других световых плоскостях, цвет которых, например, зеленый (g) и синий (b), согласован с цветовой чувствительностью пикселей фотоматрицы.

Измерение векторов скоростей Vx и Vy в параллельных лазерных плоскостях, сформированных излучением с длинами волн λr, λg и λb, осуществляется эмулированием пространственной фильтрации. Пример хроматической селекции сигналов, полученных в эксперименте, показан на Фиг. 2, где компонента спектра на частоте 10,4 Гц соответствует движущейся частице, а компонента спектра на частоте 6,3 Гц - неподвижному фону. Эмуляция пространственной фильтрации движения частицы по оси х получена на спеклах с цифровой чувствительностью λr (Фиг. 2, а), λg (Фиг. 2, b) и λb (Фиг. 2, с).

Формирование зондирующих параллельных r, g и b лазерных плоскостей на заданных расстояниях друг от друга в сочетании с хроматической селекцией изображений на фотоматрице позволяет измерять 3D поле скоростей. Измерение z-компоненты векторов скоростей, ориентированной по нормали к световым плоскостям, выполняется путем корреляционной обработки хроматически селектированных сигналов.

Обратимся к Фиг. 3. На этом чертеже показаны сечения r, g и b - световых плоскостей в исследуемой среде. Ось z является нормалью к световым плоскостям. Поскольку для формирования световых плоскостей, зондирующих исследуемую среду, обычно используются лазерные пучки, сечения световых плоскостей описываются гауссовыми функциями:,и, где wr, wg и wb - радиусы перетяжек лазерных пучков по оси z [Ю.Н. Дубнищев. Теория и преобразование сигналов в оптических системах. С.Петербург, Изд-во «Лань», 2011, 368 с.]. Выберем эти перетяжки одинаковыми: wr=wg=wb. Пусть в направлении z движется рассеивающая частица, которую, не нарушая общности рассуждений, будем считать точечной, моделируемой дельта-функцией δ(z-Vzt). Расстояния между r и g световыми плоскостями, как это показано на Фиг. 1 и Фиг. 3, равно d1, а расстояния между g и b световыми плоскостями равно d2. Световые поля, рассеянные частицей при пересечении пространственного интервала d1 между плоскостями, сформированными на длинах волн λr, и λg, описываются выражением:

где

При интегрировании в (13) мы воспользовались фильтрующим свойством δ-функции. Найдем фурье-спектр сигнала (13):

Переходя от фурье-спектра к энергетическому спектру сигнала, получаем:

Корреляционную функцию сигнала srg(t) находим, выполняя фурье-преобразование энергетического спектра Wrg(ω):

Структура полученной корреляционной функции показана на Фиг. 4. Корреляционная функция (16) имеет максимумы в точках τ=0, τ=τ1 и τ=-τ1. Поскольку , скорость движения частицы на пространственном интервале d1 определяется как

Описанный выше способ определения z-компоненты скорости заключается в определении через корреляционную функцию сигнала временного интервала τ1, в течение которого частица проходит расстояние d1 между r и g лазерными плоскостями. Совершенно аналогично определяется скорость частицы между g и b лазерными плоскостями:

Равенство скоростей Vz1=Vz2 означает, что на пространственном интервале между r и g плоскостями градиент z-компоненты скорости отсутствует. Если скорости vz1 и Vz2 не равны, градиент скорости Vz на этом пространственном интервале d1+d2 определяется как

Следовательно, предложенный способ позволяет определять и поле компонент скоростей, ориентированных нормально к r, g и b световым плоскостям, сформированным в исследуемой среде. Таким образом, задача расширения функциональных возможностей на измерение 3D поля скоростей в газовых и конденсированных средах решается предлагаемым способом за счет того, что структурированное зондирующее световое поле формируется в исследуемой среде в виде параллельных световых плоскостей на длинах волн, соответствующих цветовой чувствительности пикселей в структуре фотоматрицы, регистрирующей изображения световых плоскостей. Изображения световых плоскостей одновременно формируются на фотоматрице в телецентрической проекции. Эмуляция динамической пространственной фильтрации выполняется на пиксельных структурах, цветовая чувствительность которых согласована с длиной волны соответствующей световой плоскости. Получают корреляционные функции попарных изображений световых плоскостей и определяют пространственное распределение отношения попарного расстояния между световыми плоскостями к временному интервалу между экстремумами в соответствующих корреляционных функциях.

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



 

Похожие патенты:

Изобретение относится к измерительной технике и касается способа измерения скорости течения жидкости с рассеивающими свет частицами. Способ включает в себя освещение потока жидкости одновременно двумя пучками лазерного излучения и определение спектра мощности P12(f) отраженного сигнала.

Изобретение относится к контрольно-измерительной технике и позволяет исследовать кинематические характеристики гидропотоков. Способ, основанный на совместном использовании лазерной доплеровской анемометрии (ЛДА) и цифровой трассерной визуализации (PIV), включает установку CCD камер под углом, вычисленным с помощью корректирующего модуля пробоотбора взвеси калибровочных частиц, определение временного интервала между сериями изображений, фиксирование и запись изображений засеянных частиц и статистическое условное осреднение мгновенных полей скорости, при этом внесение корректировок в параметры пороговой чувствительности CCD камер осуществляют в продолжение исследований при уменьшении регистрируемых событий на 10% или более, либо через каждые 3 часа.

Изобретение относится к контрольно-измерительной технике и позволяет исследовать потоки жидкости и газа. Изобретение основано на совместном использовании ЛДА и PIV.

Изобретение относится к измерительной технике и может быть использовано, в частности, в прикладной метеорологии для оперативного дистанционного измерения скорости и направления ветра.

Изобретение относится к измерениям турбулентностей атмосферы с помощью лидарной системы, в частности на борту летательных аппаратов. .

Изобретение относится к измерительной технике и может быть использовано, в частности, в прикладной метеорологии для оперативного дистанционного определения скорости и направления ветра.

Изобретение относится к лазерным двухточечным оптическим расходомерам и предназначено для использования преимущественно при транспортировке природного газа. .

Изобретение относится к контрольно-измерительной технике и может быть использовано в стабилизаторах скорости движения различных подводных объектов (ПО). .

Группа изобретений относится к медицинской технике, а именно к средствам определения характеристик потока крови. Устройство содержит светоизлучающий блок, выполненный с возможностью излучения света в направлении элемента, блок регистрации света, выполненный с возможностью регистрации света, рассеянного обратно на элементе, оптический блок, выполненный с возможностью пространственного разделения участка элемента падения света элемента и участка элемента регистрации света элемента друг от друга, при этом оптический блок содержит элемент разделения светового пути, выполненный с возможностью разделения пути излучаемого света и пути обратно рассеянного света, и блок определения, выполненный с возможностью определения характеристики потока объекта на основе света, указывающего на излучаемый свет, и регистрируемого обратно рассеянного света. Способ осуществляется посредством работы устройства. Использование изобретений позволяет повысить чувствительность при измерении за счет улучшения отношения сигнал/шум. 2 н. и 11 з.п. ф-лы, 4 ил.

Изобретение относится к контрольно-измерительной технике и позволяет исследовать газожидкостные вихревые течения с любым соотношением жидкости и газа. Способ основан на совместном использовании ЛДА и PIV, включающем пропускание через измерительный объем лазерного излучения, проведение измерений с получением полного периода пульсаций, определение на основе полученной информации временного интервала между сериями изображений, по которым вычисляют мгновенные PIV поля скорости, освещение исследуемого вихревого течения когерентным лазерным светом, фиксирование изображений двумя CCD камерами, принимающими отраженный свет, и запись информации в заданном интервале времени. При этом при диагностике вихревого течения, индуцированного вращающимся прецессирующим ядром (ПВЯ), одновременно формируют сигнала скорости и опорные сигналы с помощью пьезокерамических гидрофонов или с помощью прецизионных конденсаторных микрофонов, расположенных в устройстве для реализации способа, непосредственно внутри исследовательского контейнера после завихрителя вдоль по потоку в виде отдельных сопряженных пар, диаметрально расположенных в горизонтальной и вертикальной плоскостях. Технический результат заключается в расширении технических возможностей и уменьшении ошибок, связанных с резким изменением соотношения жидкости и газа в вихревом потоке. 2 н. и 7 з.п. ф-лы, 6 ил.

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

Наверх