Detection of flow separation on the NACA 0012 profile using persistent homology methods.
Detection of flow separation on the NACA 0012 profile using persistent homology methods.
Abstract
The paper investigates the application of topological data analysis (TDA) methods, in particular persistent homology, for the quantitative assessment of structural changes in the flow field during the flow around an airfoil. Based on steady RANS calculations from the AirfRANS dataset (NACA 0012 profile), an analysis of the topological invariants of the signed vorticity field in the wake zone was performed. It is shown that in the range of stall‑free flow, the number of topological cycles (H1‑homologies) remains stable. When transitioning to the stalled mode (angles of attack α > 11.5°), a sharp increase in the topological complexity of the wake is observed: the number of identified cycles nearly doubles. The results obtained demonstrate that persistent homology methods can serve as an effective indicator for automatically detecting flow separation and analyzing the topology of vortex structures in computational hydro-aerodynamics.
1. Введение
Современная вычислительная гидроаэродинамика (CFD) генерирует огромные объемы данных, описывающих пространственно-временную структуру течений жидкости и газа. Традиционные методы оценки аэродинамических характеристик часто опираются на интегральные коэффициенты (коэффициент подъемной силы CL, коэффициент лобового сопротивления CD) или точечные измерения. Однако такие метрики не всегда позволяют количественно описать сложные изменения в топологии самого течения, возникающие при переходных режимах, таких как отрыв пограничного слоя, образование вихревых зон и турбулизация следа за профилем.
Методы топологического анализа данных (TDA), в частности алгоритмы персистентной гомологии, предлагают математический аппарат для изучения формы и структуры многомерных объектов без потери их ключевых свойств. Персистентная гомология позволяет идентифицировать инвариантные топологические характеристики — компоненты связности, циклы (дырки) и полости — и оценивать их устойчивость к шуму при изменении масштаба рассмотрения . Данный подход особенно эффективен для выделения глобальной структуры потока, такой как зоны рециркуляции и когерентные вихревые структуры, которые трудно отделить от локального шума численного решения стандартными методами пороговой фильтрации.
Теоретические основы персистентной гомологии были заложены Г. Эдельсбруннером и соавторами , а её систематизация в контексте анализа данных подробно изложена в работах Г. Карлссона и в современном обзоре Chazal и Michel . Для практической реализации метода часто используются эффективные алгоритмические решения, такие как библиотека Ripser, которая позволяет быстро вычислять диаграммы персистентности .
Применение TDA к задачам механики жидкости и газа демонстрирует высокую эффективность. В частности, методы были успешно использованы для анализа структур высокоскоростных турбулентных потоков , а также для сравнительного анализа различных численных методов моделирования турбулентности через призму топологических признаков полей скорости . Более современные исследования показывают применимость топологических дескрипторов для ансамблей гидродинамических течений , а также для анализа конвективных циклов в сетевых потоках . Отдельные работы посвящены применению персистентной гомологии непосредственно к изучению структуры течений жидкости, включая визуализацию и классификацию вихревых паттернов .
Важным направлением является интеграция TDA с методами машинного обучения для восстановления полей течения по экспериментальным данным. Например, Suzuki и др. продемонстрировали возможность оценки параметров потока исключительно на основе анализа изображений с использованием персистентной гомологии, что открывает перспективы для гибридных подходов в аэродинамическом мониторинге.
Для проверки гипотез о применении TDA в инженерной аэродинамике актуально использование открытых верификационных наборов данных. Одним из таковых является датасет AirfRANS, содержащий результаты высокодетализированных расчетов уравнений Навье-Стокса в осредненном виде (RANS) для профиля NACA 0012 при широком диапазоне углов атаки . Использование этих данных позволяет исследовать эволюцию поля завихренности от ламинарного обтекания до режимов глубокого срыва потока.
Цель работы заключается в разработке методики применения персистентной гомологии для количественной оценки изменений топологии поля знаковой завихренности ωz при обтекании аэродинамического профиля, с акцентом на детектирование зоны отрыва потока (срыва потока) через анализ роста числа топологических циклов в области высокой завихренности.
2. Методы и принципы исследования
2.1. Исходные данные
В работе использовались результаты стационарных расчётов из открытого датасета AirfRANS
. Расчёты выполнены методом RANS (Reynolds-Averaged Navier-Stokes) в пакете OpenFOAM с использованием модели турбулентности k-ω SST. Число Рейнольдса, построенное по хорде профиля и параметрам набегающего потока, варьировалось в диапазоне Re = 2·106 – 6·106, что соответствует дозвуковому несжимаемому режиму течения (Ma < 0,3). Пограничный слой разрешался полностью (low-Re подход) без применения пристенных функций, что обеспечивает точное восстановление касательных напряжений на поверхности профиля. Геометрия профиля задавалась параметрически для семейств NACA 4-digit и 5-digit, включая профиль NACA 0012, при углах атаки α ∈ [-5°; +15°].2.2. Расчёт поля завихренности
Для выделения областей интенсивного вращения жидкости использовалась проекция вектора завихренности на ось z, перпендикулярную плоскости течения:
где u, v — компоненты вектора скорости.
Градиенты скорости вычислялись на основе конечно-элементной аппроксимации на треугольной неструктурированной сетке. Для каждого треугольника с вершинами (xi, yi) и скоростями (ui, vi), i = 1, 2, 3 производные определялись аналитически через решение системы линейных уравнений. Полученные на элементах значения интерполировались в узлы сетки методом взвешенного усреднения с сохранением знака.
2.3. Персистентная гомология
Для анализа топологической структуры поля ωz применялся метод персистентной гомологии , , , . Поскольку на этапе стохастического прореживания (субсэмплинга) исходная сеточная связность утрачивается, для анализа облака координат применялась фильтрация Вьеториса-Рипса (Vietoris-Rips filtration), реализуемая библиотекой Ripser
, . Данная фильтрация строит симплициальный комплекс на основе попарных евклидовых расстояний между точками облака. Значением функции в каждой точке служила нормированная завихренность:В работе анализировались одномерные гомологии (H1), соответствующие циклам (замкнутым структурам) в поле завихренности.
2.4. Параметры фильтрации
Для исключения топологического шума применялся порог по персистентности:
где b — время рождения (birth) и d — время смерти (death) топологической особенности.
В расчётах использовалось pmin = 0,05 (5% от нормированного диапазона), что обеспечивает баланс между сохранением физически значимых структур и подавлением численного шума.
Границы окна анализа (x > 0,8, |y| < 0,8) выбраны так, чтобы охватывать зону развития следа для всех исследованных углов атаки, исключая при этом пограничный слой на поверхности профиля.
Порог 90-го перцентиля обеспечивает фокусировку на ядрах вихревых структур, исключая фоновые флуктуации.
2.5. Статистическая устойчивость
Для оценки влияния стохастического прореживания каждый режим расчитывался 11 раз с различными начальными значениями генератора случайных чисел (seed = 40...400). В качестве финальной метрики использовались среднее число циклов H1 и его стандартное отклонение, что позволяет количественно оценить устойчивость индикатора к вариациям выборки.
2.6. Анализ зоны следа
Для детектирования срыва потока TDA применялся не ко всему полю, а к выделенной области следа (x > 0,8 хорды, |y| < 0,8), где отбирались точки с завихренностью выше 90-го перцентиля. Это позволило сфокусировать анализ на топологии вихревых структур, исключив влияние пограничного слоя на профиле.
3. Основные результаты
3.1. Визуализация полей течения
Рисунок 1 - Знаковая завихренность (линейная шкала)
Рисунок 2 - Знаковая завихренность (логарифмическая шкала)
Рисунок 3 - Область обратного течения (ux <0)
3.2. Топологический анализ
Для количественной оценки изменений структуры течения применялся метод персистентной гомологии к областям высокой завихренности в зоне следа (x > 0,8 хорды, |y| <0,8).
В табл. 1 приведены результаты расчёта числа топологических циклов (H₁-гомологий) для 21-го исследованного режима.
Таблица 1 - Топологические характеристики поля завихренности при различных углах атаки
Угол атаки α, ° | Коэффициент CL | Число циклов H₁ |
-4,2 | -0,321 | 26 ± 2,7 |
-3,34 | -0,310 | 24,0 ± 1,7 |
-2,65 | -0,030 | 27,6 ± 2,7 |
-1,91 | -0,061 | 25,3 ± 4,9 |
-1,58 | 0,051 | 26,0 ± 4,0 |
-0,55 | -0,053 | 25,5 ± 3,1 |
-0,47 | 0,506 | 24,6 ± 3,5 |
0,42 | 0,067 | 28,4 ± 4,1 |
1,33 | 0,149 | 22,4 ± 3,2 |
5,28 | 0,682 | 26,8 ± 3,1 |
8,2 | 1,079 | 25,5 ± 2,8 |
9,42 | 1,040 | 39,5 ± 3,1 |
9,92 | 1,093 | 42,6 ± 5,5 |
10,15 | 1,322 | 41,6 ± 5,2 |
10,42 | 1,358 | 36,1 ± 2,6 |
10.54 | 1,185 | 37,2 ± 3,8 |
11,32 | 1,388 | 36,9 ± 2,5 |
11,51 | 1,766 | 44,3 ± 4,6 |
12,12 | 1,675 | 43,5 ± 2,7 |
13,71 | 1,620 | 45,4 ± 3,1 |
14,04 | 1,847 | 51,9 ± 6,4 |
Статистический анализ зависимости числа топологических циклов H1 от угла атаки α (рис. 4) позволил выделить три характерных режима обтекания.
Рисунок 4 - Зависимость числа топологических циклов H1 от угла атаки α
В диапазоне отрицательных и малых положительных углов атаки наблюдается плато топологической стабильности. Число выявленных циклов варьируется в узком диапазоне 22-28 при среднем значении около 25±3. Это соответствует устойчивой конфигурации узкого турбулентного следа и пристенного пограничного слоя, где топологическая сложность определяется лишь тонкими структурами сдвигового слоя.
Переходный режим локального отрыва (9,4° ≤ α <11,5°).
Начиная с угла атаки α ≈ 9,4°, фиксируется статистически значимый рост топологической сложности. Число циклов увеличивается до значений 36-42. Этот рост обусловлен зарождением локальной зоны отрыва у задней кромки профиля. Формирующийся небольшой рециркуляционный пузырь усложняет геометрию поля завихренности, создавая новые замкнутые контуры, фиксируемые алгоритмом TDA.
Режим полного срыва (α ≥ 11,5°).
При достижении критического угла α ≈ 11,5° происходит фазовый переход топологии течения. Наблюдается резкий скачок числа циклов до значений 44–52 (рост более чем на 80% по сравнению с безотрывным режимом). Это соответствует формированию развитой зоны отрыва, занимающей значительную часть верхней поверхности профиля. Высокое значение H1 в этой зоне отражает сложную структуру крупных вихревых образований и границ обширной зоны обратного течения.
Таким образом, предложенный топологический индикатор H1 демонстрирует высокую чувствительность к качественной перестройке структуры следа. Резкий рост числа циклов при α ≥ 9,4° и стабилизация на высоком уровне при α ≥ 11,5° позволяют однозначно детектировать переход потока в срывной режим без прямого анализа силовых коэффициентов.
4. Обсуждение
Физическая интерпретация выявленных топологических инвариантов требует учёта специфики осреднённых по Рейнольдсу уравнений Навье–Стокса (RANS). В нестационарных расчётах (LES/DNS) одномерные гомологии (H1-циклы) в поле завихренности традиционно соответствуют дискретным вихревым ядрам или замкнутым линиям тока. Однако в стационарных RANS-расчётах пульсационная составляющая течения моделируется через замыкающие соотношения турбулентной вязкости, а поле скоростей представляет собой усреднённую по времени картину. В данном контексте выявленные топологические циклы отражают не мгновенные вихри, а сложную топологию средней зоны отрыва и сдвиговых слоёв, окаймляющих рециркуляционный пузырь.
Расширенный статистический анализ позволил выявить трёхстадийную эволюцию структуры течения, что уточняет понимание процесса срыва:
Режим безотрывного обтекания (α < 9,4°).
Наблюдается стабильное «плато» топологической сложности (среднее значение ≈ 25 циклов). Это соответствует узкому турбулентному следу и пристенному пограничному слою, где топологические особенности определяются лишь тонкими структурами сдвига.
Переходная зона локального отрыва (9,4° ≤ α < 11,5°).
При достижении угла α ≈ 9,4° начинается плавный рост числа циклов. Физически это соответствует зарождению локальной зоны обратного течения у задней кромки профиля. Формирующийся небольшой рециркуляционный пузырь усложняет геометрию поля завихренности, создавая новые замкнутые контуры.
Режим полного срыва (α ≥ 11,5°).
При критическом угле α ≈ 11,5° происходит топологический фазовый переход. Наблюдается резкий скачок числа циклов до значений ≈ 45 и более (рост более чем на 80%). Это количественно подтверждает формирование развитой зоны отрыва, занимающей значительную часть верхней поверхности профиля. Высокое значение H1 в этой зоне отражает сложную структуру крупных вихревых образований и границ обширной области пониженного давления.
Важно отметить ограничения, накладываемые RANS-подходом на топологический анализ. Поскольку метод не разрешает нестационарное вихреобразование (например, вихревую дорожку Кармана), TDA-метрики отражают лишь усреднённое состояние потока. Это делает предложенный подход оптимальным для задач стационарного проектирования, верификации расчётных сеток и анализа интегральной структуры отрыва, но требует осторожности при интерпретации в контексте нестационарной аэродинамики. Для полного описания эволюции вихревых структур в отрывной зоне целесообразно в дальнейшем применять методику к нестационарным данным (LES/DES), где диаграммы персистентности смогут отслеживать динамику рождения, слияния и разрушения дискретных вихрей во времени. Перспективность такого подхода подтверждается недавними работами по анализу вихревых следов методами персистентной гомологии и графовыми методами TDA для транзиентных режимов течения .
С практической точки зрения, предложенный топологический индикатор демонстрирует высокий потенциал для автоматизированного детектирования срыва потока. В отличие от традиционных методов, опирающихся на анализ коэффициента давления Cp или касательного напряжения на стенке τw, которые требуют привязки к геометрии профиля и могут быть чувствительны к численному шуму вблизи критических точек, TDA-метрика является инвариантной к системе координат и масштабу расчётной области. Резкий рост числа H1-циклов может быть использован как триггер для адаптивного сгущения сетки в зоне отрыва или как критерий остановки итерационного процесса при оптимизации формы профиля. Кроме того, наблюдаемая связь между топологической сложностью следа и аэродинамическими коэффициентами (CL, CD) открывает путь к созданию суррогатных моделей, связывающих топологические инварианты полей течения с интегральными характеристиками летательного аппарата без необходимости проведения полного CFD-анализа на каждом шаге оптимизации. Подобные подходы уже успешно применяются в смежных областях, таких как химическая инженерия , а интеграция TDA с методами машинного обучения может существенно ускорить процесс поиска аэродинамических конфигураций с заданными характеристиками устойчивости к срыву.
5. Заключение
В работе разработана и апробирована методика применения персистентной гомологии для топологического анализа полей завихренности в задачах вычислительной гидроаэродинамики. На основе открытого датасета AirfRANS (профиль NACA 0012) проведён количественный анализ эволюции топологических инвариантов поля знаковой завихренности ωz в зоне следа при изменении угла атаки.
Основные результаты исследования можно сформулировать следующим образом:
1. Предложенный TDA-подход, основанный на фильтрации Вьеториса-Рипса, позволяет выделять инвариантные топологические характеристики (одномерные гомологии H1) структуры течения, устойчивые к численному шуму и локальным флуктуациям полей скорости.
2. Анализ зависимости числа топологических циклов от угла атаки выявил чёткое разделение режимов обтекания:
– безотрывный режим (α < 9,4°): характеризуется стабильным «плато» топологической сложности со средним значением 25±3 цикла, что соответствует узкому турбулентному следу и пристенному пограничному слою;
– переходная зона (9,4° ≤ α < 11,5°): наблюдается статистически значимый рост числа циклов до 40 - 43, обусловленный зарождением локальной зоны отрыва у задней кромки профиля;
– режим полного срыва (α ≥ 11,5°): фиксируется выход на новое плато с высокой топологической сложностью (45 - 52 цикла), что количественно подтверждает формирование развитой рециркуляционной области на всей верхней поверхности профиля.
3. Резкий скачок топологической сложности при α ≈ 9,4° и стабилизация индикатора при α ≥ 11,5° полностью согласуются с классическими представлениями о срыве потока на профиле NACA 0012. Топологический индикатор реагирует на качественную перестройку структуры следа (появление замкнутых контуров вихревых зон), а не на интегральные силовые характеристики.
4. Топологический индикатор демонстрирует высокую чувствительность к изменению режима обтекания и может рассматриваться как координатно-независимый, масштабно-инвариантный критерий для автоматического детектирования срыва потока в стационарных CFD-расчётах.