Table of Contents

Введение: почему велоциметрия изображения частиц имеет значение в динамике жидкости

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

Основные принципы велоциметрии изображения частиц

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

Оригинальное название: The Seeding Particle: The Heart of PIV

Выбор правильных частиц посева имеет решающее значение. Частицы должны быть достаточно малы, чтобы точно следовать движению жидкости (число Стокса намного меньше 1), но достаточно велики, чтобы рассеивать достаточно света для обнаружения. Общие материалы включают:

  • Жидкие потоки: полые стеклянные сферы (10-100 мкм), полиамидные шарики или флуоресцентные частицы для минимизации отражений стенок.
  • Газовые потоки: капли масла (DEHS, оливковое масло), диоксид титана (TiO2) или частицы дыма (1-10 мкм).
  • Двухфазные потоки: частицы, выбранные для соответствия плотности одной фазы или с использованием коэффициента преломления.

Плотность посева должна быть достаточно высокой, чтобы обеспечить надежную корреляцию, но не настолько высокой, чтобы изображения частиц чрезмерно перекрывались.Как правило, 10-25 частиц на окно опроса считаются оптимальными.

Лазерное освещение и визуализация

Системы PIV обычно используют двухимпульсные Nd:YAG-лазеры (532 нм зеленым), способные стрелять двумя импульсами с очень коротким, регулируемым разделением времени (микросекунды до миллисекунд). Лазерный луч формируется в тонкий лист с использованием цилиндрических и сферических линз. Толщина листа - обычно между 0,5 мм и 2 мм - определяет объем измерения в направлении вне плоскости. Высокоскоростные камеры, часто с разрешением 4-16 мегапикселей и частотой кадров от нескольких Гц до нескольких кГц, захватывают изображения частиц. Для 3D-технологий, таких как стерео-PIV, используются две камеры, каждая из которых просматривает лазерный лист под другим углом.

Типы систем и конфигураций PIV

За последние три десятилетия PIV превратился из простого 2D-планарного метода в семейство мощных техник. Выбор зависит от сложности потока и требуемого пространственного/временного разрешения.

2D Планарный PIV

Самая базовая конфигурация: одна камера фиксирует рассеянный частицами свет в одном лазерном листе. Она выдает два скоростных компонента (u, v) в плоскости листа. Этого достаточно для многих устойчивых и слегка трехмерных потоков. Методика проста в реализации и остается рабочей лошадкой многих лабораторий.

Стерео ПИВ (2D-3C)

Используя две камеры, расположенные под наклонными углами к лазерному листу, стерео PIV восстанавливает все три компонента скорости (u, v, w) в освещенной плоскости. Компонент вне плоскости выводится из несоответствия между двумя видами камеры с использованием методов калибровки (например, модель пинхола или полиномиальное картирование третьего порядка). Стерео PIV необходим для потоков со значительным движением вне плоскости, таких как вихри или потоки завихрений.

Томографический PIV (3D-3C)

Наиболее продвинутый планарный подход. Три или более камеры просматривают толстый лазерный объем (например, 5-10 мм). Алгебраические алгоритмы томографической реконструкции (например, MART-Multiplicative Algebraic Reconstruction Technique) реконструируют распределение 3D-частиц в объеме. Затем 3D-перекрестная корреляция дает все три компонента скорости по истинной трехмерной области. Томо-PIV используется для высокотурбулентных потоков и сложных геометрий, таких как струйные потоки или струйные потоки.

Микро-ПИВ (μPIV)

Адаптация ПИВ для микромасштабных потоков (каналы 10–500 мкм). Вместо лазерного листа используется объемная подсветка, а глубина резкости объектива микроскопа определяет плоскость измерения. Флуоресцентные частицы и фильтры дальнего действия подавляют фоновый шум. μPIV широко применяется в микрофлюидиках, биологических потоках и исследованиях теплопередачи в малых масштабах.

Высокоскоростной PIV

Используя лазеры скорости кГц и быстрые камеры, высокоскоростная PIV захватывает разрешённые временем последовательности скоростных полей. Это позволяет изучать переходные явления, неустойчивости потока и спектральный анализ турбулентности. Современные системы могут получать данные на частоте 10–50 кГц, обеспечивая разрешённое во времени представление об эволюции потока.

Пошаговая цепочка обработки данных PIV

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

  1. Предварительная обработка изображения: Фоновое вычитание, нормализация интенсивности и фильтрация (например, минимальное или среднее вычитание скольжения) для удаления лазерных отражений и неравномерного освещения.
  2. Выбор окна допроса: Изображения делятся на квадратные окна (обычно 32×32 или 64×64 пикселей). Выбор балансирует пространственное разрешение (маленькие окна) против корреляционной прочности (большие окна содержат больше частиц).
  3. Анализ перекрестной корреляции: Каждое окно на изображении А коррелирует с областью поиска на изображении В. Обнаружение пика на карте корреляции дает вектор смещения с точностью до субпикселей (например, с использованием гауссовой пиковой установки).
  4. Проверка и удаление выпадающих: Ложные векторы идентифицируются медианной фильтрацией, нормализованными проверками гистограммы или сравнениями RMS. Недействительные векторы заменяются интерполяцией от соседей.
  5. Постобработка: Векторные поля сглаживаются (например, гауссовский фильтр) для уменьшения шума. Производные величины, такие как вихрь, скорость деформации и напряжения Рейнольдса, можно вычислить по пространственным градиентам.
  6. Усреднение сборки: Для статистически неподвижных потоков усредняются сотни или тысячи мгновенных полей для получения профилей средней скорости и турбулентной статистики.

Применение PIV в экспериментальной механике жидкостей

PIV применяется практически во всех отраслях гидродинамики. Ниже приведены подробные примеры, иллюстрирующие его универсальность.

Аэрокосмическая и аэронавтическая

При испытаниях аэродинамической трубы PIV отображает поток вокруг аэродинамических фольг, крыльев и целых моделей самолетов. Он показывает пузыри разделения, передовые вихри на дельта-крыльях и турбулентность пробуждения. Например, исследования на аэродинамической поверхности NACA 4412 при высоком угле атаки используют PIV для захвата механизмов сваливания. Данные PIV непосредственно подтверждают моделирование вычислительной динамики жидкости (CFD), уменьшая потребность в большом количестве кранов давления.

Автомобильные и турбомашинные

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

Биомедицинская инженерия

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

Экологические и геофизические потоки

Полевой PIV (иногда называемый крупномасштабным PIV) использует естественные трассеры, такие как снежинки, пыль или плавучие пузырьки, для измерения полей ветра в атмосферном пограничном слое. В гидравлической технике PIV применяется в открытых каналах для изучения транспорта осадков и турбулентных структур вблизи русла реки.

Промышленные процессы

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

Преимущества PIV перед другими методами измерения потока

Полноценное измерение

В отличие от зондов с горячей проводкой или трубок Пито, которые измеряют в одной точке, PIV обеспечивает полную карту мгновенной скорости. Это имеет решающее значение для захвата структур потока, которые не известны априори.

Ненавязчивый

В поток не вводится зонд, поэтому измерение не нарушает изучаемое явление. Это особенно важно для тонких потоков, таких как биологические жидкости или чувствительная турбулентность.

Высокое пространственное разрешение

Современные ПИВ могут достигать векторного расстояния от 0,1 мм (для микро-ПИВ) до нескольких миллиметров (в больших аэродинамических трубах). Этого достаточно для разрешения вихрей и сдвига слоев.

Гибкость

PIV работает в газах, жидкостях и многофазных потоках. При правильном посеве и оптическом доступе его можно применять в экстремальных условиях (высокая температура, высокое давление или вакуум).

Проблемы и ограничения

Несмотря на свою силу, ПИВ имеет неотъемлемые ограничения, которые должен учитывать каждый экспериментатор.

Стоимость и сложность

Типичная планарная система PIV стоит 50 000-150 000 долларов США для лазера, камеры, оптики и программного обеспечения. Томографические или высокоскоростные установки могут превышать 300 000 долларов США. Кроме того, точное выравнивание лазеров и камер требует квалифицированного персонала.

Оптический доступ

Для установки требуется прозрачное окно (или погружение в жидкость с соответствующим показателем преломления), чтобы впустить лазерный лист и рассеянный свет. Многие промышленные геометрии не обеспечивают такой доступ.

Убыток с внепланового

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

Засев и отрыв частиц

В высокоскоростных потоках газа или ударных волнах частицы могут не следовать за ускорениями верно (конечное число Стокса). Это вводит систематические ошибки. Аналогично, в потоках с сильными градиентами плотности (например, горение) скорость частиц может отклоняться от скорости жидкости.

Вычислительная нагрузка

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

Количественная неопределенность

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

Последние достижения и будущие направления

Shake-the-Box (STB) и Лагранжианское отслеживание частиц

Новая парадигма под названием «Shake-the-Box» (разработанная в DLR Göttingen) сочетает в себе многокамерные изображения с итеративным обнаружением и отслеживанием частиц. Она дает индивидуальные траектории частиц по длинным последовательностям, обеспечивая поля ускорения и градиенты давления от производной материала. STB предлагает значительно более высокое пространственное разрешение, чем традиционные томографические PIV, особенно для турбулентных потоков.

В реальном времени PIV

С появлением программируемых ворот (FPGA) и быстрых алгоритмов корреляции, PIV в реальном времени теперь возможен для управления процессами и мониторинга аэродинамической трубы. Скорости 10-40 векторных полей в секунду позволяют операторам регулировать положение модели или условия потока на лету.

PIV в сочетании с другими методами

Гибридные методы, такие как PIV/PLIF (планарная лазерная флуоресценция) или PIV/schlieren, обеспечивают одновременные измерения скорости и скалярного поля. Это особенно эффективно для исследований горения, смешивания и теплопередачи.

Миниатюризация и портативные системы

Для полевых работ создаются компактные, работающие от батарей системы ПИВ с использованием твердотельных лазеров и небольших камер, которые весят менее 10 кг и которые были развернуты на небольших дронах или в водотоках для измерения потока окружающей среды.

Машинное обучение для PIV

Подходы глубокого обучения (например, сверточные нейронные сети) используются для замены перекрестной корреляции сквозной оптической оценкой потока. Такие методы, как PIV-DCNN, могут обрабатывать большие градиенты смещения и уменьшать шум, хотя они требуют обширных данных обучения.

Лучшие практики для успешного эксперимента PIV

Чтобы получить достоверные результаты, исследователи должны следовать этим рекомендациям:

  • Выполните тщательную калибровку: используйте микрометровую цель или пластину с точечной сеткой для отображения координат изображения на физические координаты.
  • Оптимизируйте время разделения импульсов (Δt), чтобы частицы перемещались на 4-8 пикселей между воздействиями. Слишком большие приводят к потере вне плоскости; слишком маленькие создают плохой динамический диапазон.
  • Проверить качество посева: обеспечить равномерное распределение, адекватную концентрацию и отсутствие частиц, прилипающих к стенкам.
  • Проверяйте оптическую установку: минимизируйте отражения, покрывая поверхности антибликовой краской или используя флуоресцентные частицы с фильтром дальнего действия.
  • Проведите исследование конвергенции: приобретите достаточное количество статистически независимых образцов, чтобы получить стабильные средние и турбулентные величины.
  • Неопределенности документов: сообщайте о систематических и случайных ошибках в соответствии со стандартами, такими как AIAA или DIN для PIV.

Заключение

Велоциметрия изображения частиц превратилась в надежную полнопольную методику измерения, которая лежит в основе бесчисленных открытий в механике жидкости. Ее эволюция от рудиментарных 2D-планарных установок до высокоскоростных, 3D и систем, улучшающих машинное обучение, отражает неустанное стремление поля к большей детализации и точности. В то время как проблемы остаются - особенно в стоимости, оптическом доступе и количественной оценке неопределенности - продолжающиеся достижения в лазерной технологии, датчиках камеры и вычислительных алгоритмах продолжают расширять границы того, что можно измерить. Для любого экспериментального механика жидкости овладение PIV больше не вариант, а необходимость.

Для дальнейшего чтения по теории и практике PIV см. классический учебник Велоциметрия изображения частиц Раффеля, Виллерта, Верли и Компенханса, или посетите технические ресурсы ведущих производителей, таких как Dantec Dynamics и LaVision.