БЕССЕТОЧНЫЙ ПЛАНАРНЫЙ МЕТОД PARTICLE IMAGE VELOCIMETRY (2022)
На сегодняшний день многопроходный метод PIV (Particle Image Velocimetry) широко используется в области экспериментальной механики жидкости и газа из-за его высокой надежности при решении практических задач. Однако он имеет известное ограничение, связанное с ошибками, возникающими при вычислении производных скорости, необходимых для деформации обрабатываемых PIV-изображений при повышении производительности метода. Поскольку количество ошибок увеличивается с применением схем более высокого порядка, на практике чаще всего ограничиваются первым порядком, что в свою очередь приводит к снижению пространственного разрешения. В данной работе предлагается метод, допускающий применение схем более чем второго порядка, что позволяет заметно повысить точность измерения скорости и ее производных и тем самым увеличить пространственное разрешение. Метод не требует восстановления ошибочных векторов скорости, позволяет избежать численного расчета производных скорости и легко применим на практике.
Идентификаторы и классификаторы
- eLIBRARY ID
- 50104877
За последние три десятилетия метод цифровой трассерной визуализации (PIV — Particle Image Velocimetry) стал надежным инструментом измерения различных кинематических характеристик потока жидкости или газа. За это время разработано множество алгоритмов, расширяющих возможности метода PIV и повышающих достоверность получаемых данных [1, 2]. Идея метода заключается в получении и дальнейшей обработке цифровых изображений частиц-трассеров, предварительно засеиваемых в исследуемый поток жидкости. При отслеживании трассеров полагается, что они движутся с локальной скоростью потока. Визуализация осуществляется за счет света, рассеиваемого этими частицами (рис. 1). Для подсветки требуемой области потока чаще всего используется импульсный лазер или лазер непрерывного действия, а также светодиодная подсветка совместно с необходимой оптической системой. Для получения последовательности изображений используются фото- или видеокамеры, установленные перпендикулярно плоскости светового ножа.
Список литературы
- M. P. Tokarev, D. M. Markovich, and A. V. Bil’skii, “Adaptive Particle Image Processing Algorithms for Calculating Instantaneous Velocity Fields”, Vychisl. Tekhnol. 12 (3), 109-131 (2007). EDN: KVVOZX
- C. J. Kähler, T. Astarita, P. P. Vlachos, et al., “Main Results of the 4th International PIV Challenge”, Exp. Fluids 57 (6), 1-71 (2016). DOI: 10.1007/s00348-016-2173-1 EDN: YCPPPJ
- M. Raffel, C. E. Willert, F. Scarano, et al., Particle Image Velocimetry: A Practical Guide (Springer, Cham, 2018).
- B. J. Kim and H. J. Sung, “A Further Assessment of Interpolation Schemes for Window Deformation in PIV”, Exp. Fluids 41 (3), 499-511 (2006). DOI: 10.1007/s00348-006-0177-y
- M. P. Tokarev, Development of Algorithms and Software for Processing Images in Digital Tracer Visualization Methods, Candidate’s Dissertation in Technical Sciences (Institute of Thermophysics, Novosibirsk, 2010). EDN: QEUCUH
- H. T. Huang, H. E. Fiedler, and J. J. Wang, “Limitation and Improvement of PIV. Part I: Limitation of Conventional Techniques due to Deformation of Particle Image Patterns”, Exp. Fluids 15 (3), 168-174 (1993). DOI: 10.1007/BF00189883
- H. T. Huang, H. E. Fiedler, and J. J. Wang, “Limitation and Improvement of PIV. Part II: Particle Image Distortion, a Novel Technique”, Exp. Fluids 15 (4), 263-273 (1993). DOI: 10.1007/BF00223404
- F. Scarano, “Iterative Image Deformation Methods in PIV”, Meas. Sci. Technol. 13 (1), R1-R19 (2001). DOI: 10.1088/0957-0233/13/1/201
- F. Scarano and M. L. Riethmuller, “Iterative Multigrid Approach in PIV Image Processing with Discrete Window Offset”, Exp. Fluids 26 (6), 513-523 (1999). DOI: 10.1007/s003480050318 EDN: AWCHLH
-
F. Scarano and M. L. Riethmuller, "Advances in Iterative Multigrid PIV Image Processing", Exp. Fluids 29 (Suppl 1), S051-S060 (2000). DOI: 10.1007/s003480070007 EDN: ATDKWZ
-
S. T. Wereley and C. D. Meinhart, "Second-Order Accurate Particle Image Velocimetry", Exp. Fluids 31 (3), 258-268 (2001). DOI: 10.1007/s003480100281 EDN: ARULPL
-
L. Gui and J. M. Seiner, "An Improvement in the Nine-Point Central Difference Image Correction Method for Digital Particle Image Velocimetry Recording Evaluation", Meas. Sci. Technol. 15 (10), 1958-1964 (2004). DOI: 10.1088/0957-0233/15/10/002 EDN: HBJASJ
-
J. Westerweel, D. Dabiri, and M. Gharib, "The Effect of a Discrete Window Offset on the Accuracy of Cross-Correlation Analysis of Digital PIV Recordings", Exp. Fluids 23 (1), 20-28 (1997). DOI: 10.1007/s003480050082 EDN: AWRXAH
-
J. Chen and J. Katz, "Elimination of Peak-Locking Error in PIV Analysis Using the Correlation Mapping Method", Meas. Sci. Technol. 16 (8), 1605-1618 (2005). DOI: 10.1088/0957-0233/16/8/010
-
J. Westerweel and F. Scarano, "Universal Outlier Detection for PIV Data", Exp. Fluids 39 (6), 1096-1100 (2005). DOI: 10.1007/s00348-005-0016-6
-
S. Tiwari and J. Kuhnert, Grid Free Method for Solving the Poisson Equation, Volume 25 of Berichte (Fraunhofer Institut Techno- und Wirtschaftsmathematik, Kaiserslautern, 2001).
-
M. Stanislas, K. Okamoto, C. J. Kähler, et al., "Main Results of the Third International PIV Challenge", Exp. Fluids 45 (1), 27-71 (2008). DOI: 10.1007/s00348-008-0462-z
-
B. Lecordier and J. Westerweel, "The EUROPIV Synthetic Image Generator (S.I.G.)", in Particle Image Velocimetry: Recent Improvements (Springer, Berlin, 2004), pp. 145-161. DOI: 10.1007/978-3-642-18795-7_11
-
M. Frigo and S. G. Johnson, "The Design and Implementation of FFTW3", Proc. IEEE 93 (2), 216-231 (2005). DOI: 10.1109/JPROC.2004.840301
Выпуск
Методы и алгоритмы вычислительной математики и их приложения.
Параллельные программные средства и технологии.
Другие статьи выпуска
В статье рассмотрен подход к автоматизированному распараллеливанию программ для кластеров с помощью системы SAPFOR (System FOR Automated Parallelization). Главной целью системы SAPFOR является автоматизация процесса отображения последовательных программ на параллельные архитектуры в модели DVMH, которая является моделью программирования, основанной на директивах. Помимо этого, система SAPFOR позволяет выполнять автоматически некоторый класс преобразований над исходным кодом программы по запросу пользователя через графический интерфейс. На определенных классах задач пользователь системы SAPFOR может рассчитывать на полностью автоматическое распараллеливание, если программа была написана или приведена к потенциально параллельному виду. Также в статье описаны подходы к построению схем распределения данных и вычислений на распределенную память в модели DVMH. Эффективность полученных алгоритмов построения схем аспределения данных и вычислений продемонстрирована на примере некоторых приложений из пакета NAS Parallel Benchmarks.
В данной статье рассматривается балансно-характеристический численный метод решения гиперболических систем уравнений на треугольных расчетных сетках. Описываются основные шаги алгоритма на примере решения двумерных уравнений мелкой воды. Метод верифицирован и проведено его сравнение с методами, разработанными другими авторами, на основных тестах для уравнений мелкой воды над ровным дном.
In this paper, we consider the Fredholm integral equations of the second kind and construct a new iterative scheme associated to the Nyström method, which was elaborated by Atkinson to approximate the solution over a large interval. Primarily, we demonstrate the inability to generalize the Atkinson iterative methods. Then, we describe our modified generalization in detail and discuss its advantages such as convergence of the iterative solution to the exact solution in the sense norm of the Banach space С0[a,b]. Finally, we give a numerical examples to illustrate the accuracy and reliability of our generalization.
В статье рассматривается параллельный алгоритм решения задач глобальной оптимизации и обсуждается его реализация с использованием набора инструментов Intel oneAPI. Предполагается, что целевая функция задачи задана как “черный ящик” и удовлетворяет условию Липшица. Изложенный в статье параллельный алгоритм использует схему редукции размерности на основе кривых Пеано, которые непрерывно и однозначно отображают отрезок вещественной оси на гиперкуб. В качестве средства для реализации параллельного алгоритма использован инструментарий Intel oneAPI, который позволяет писать один код как для центрального процессора, так и для графических ускорителей. Приведены результаты вычислительных экспериментов, полученные при решении серии сложных задач многоэкстремальной оптимизации.
Исследованы закономерности процессов испарения и конденсации чистого пара в методе решеточных уравнений Больцмана. Выполнено моделирование этих процессов при постоянных во времени потоках пара на границе расчетной области. Показано, что в этом случае осуществляются квазистационарные режимы испарения и конденсации. Предложен простой численно эффективный метод задания потока пара на плоской границе расчетной области путем вычисления функций распределения на входящих характеристиках метода решеточных уравнений Больцмана. В расчетах показано, что поток массы при испарении плоской поверхности пропорционален разности плотностей насыщенного и окружающего пара при данной температуре поверхности, что хорошо согласуется с законом Герца-Кнудсена. Результаты трехмерного и одномерного моделирования методом решеточных уравнений Больцмана совпадают с высокой точностью. Показано, что отношение разности плотностей к потоку вещества на границе фаз при заданной температуре линейно зависит от времени релаксации как для испарения, так и для конденсации. Исследовано влияние температуры на интенсивность потоков испарения и конденсации чистого пара. Обнаружена зависимость процессов испарения и конденсации от времени релаксации, которое определяет кинематическую вязкость флюида.
Процесс распараллеливания программ может быть затруднён ввиду их оптимизации под последовательное выполнение. Из-за этого полученная параллельная версия может быть неэффективной, а в некоторых случаях распараллеливание оказывается невозможным. Решить указанные проблемы помогают преобразования исходного кода программ. В данной статье рассматривается реализации в системе автоматизированного распараллеливания SAPFOR (System FOR Automated Parallelization) преобразований последовательных Фортран-программ, позволяющих облегчить работу пользователя в системе и существенно снизить трудоемкость распараллеливания программ. Применение реализованных преобразований в системе SAPFOR продемонстрировано на прикладной программе, решающей систему нелинейных дифференциальных уравнений в частных производных. Также было произведено сравнение производительности полученной параллельной версией с версиями, распараллелеными вручную с использованием DVM и MPI технологий.
Работа посвящена построению параллельных алгоритмов решения прямой начально-краевой задачи и обратной задачи о восстановлении правой части для уравнения диффузии с дробной производной по времени. При использовании дополнительной информации о решении в некоторой внутренней точке обратная задача сводится к прямой задаче для вспомогательного уравнения. После применения конечно-разностных схем задачи сводятся к решению систем линейных алгебраических уравнений. Разработанные алгоритмы основаны на методе параллельной прогонки и реализованы для многоядерных процессоров с использованием технологии OpenMP. Проведены численные эксперименты для исследования производительности разработанных алгоритмов.
Издательство
- Издательство
- МГУ
- Регион
- Россия, Москва
- Почтовый адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- Юр. адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- ФИО
- Садовничий Виктор Антонович (РЕКТОР)
- E-mail адрес
- info@rector.msu.ru
- Контактный телефон
- +7 (495) 9391000
- Сайт
- https://msu.ru/