Модельное предсказательное планирование пути в потоке Навье-Стокса с использованием редуцированных моделей на основе POD

Основная информация

• ID статьи: 2511.22123
• Название: Model Predictive Path Planning in Navier-Stokes Flow with POD-Based Reduced-Order Models
• Авторы: Adam Waterman, Martin Guay (Queen's University, Canada)
• Классификация: math.OC (Оптимизация и управление), cs.SY (Системы и управление), eess.SY (Системы и управление)
• Дата публикации: 27 ноября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2511.22123

Аннотация

В данной работе предложена схема оптимального генерирования траектории в системах, управляемых потоком и описываемых уравнениями Навье-Стокса, путём объединения редуцированных моделей (ROM) на основе собственного ортогонального разложения (POD) и модельного предсказательного управления (MPC). Метод включает: (i) аппроксимацию поля скоростей из данных посредством POD снимков и ортогональной проекции; (ii) вывод динамической модели Галёркина в редуцированных координатах; (iii) применение MPC для планирования управляющих воздействий, направляющих агента через предсказанное поле потока при соблюдении ограничений на состояние и исполнительные механизмы. Благодаря использованию редуцированного моделирования метод обеспечивает управление в реальном времени в высокомерных потоковых средах. Моделирование демонстрирует точную реконструкцию поля потока и эффективное генерирование траектории в реальных ветровых условиях.

Исследовательский контекст и мотивация

Определение проблемы

Данное исследование направлено на решение проблемы планирования траектории в реальном времени в системах, управляемых потоком, в частности для высокомерных потоковых сред, описываемых уравнениями Навье-Стокса. Типичные сценарии применения включают:

• Удержание позиции и навигацию стратосферных аэростатов
• Планирование пути подводных планеров
• Другие мобильные платформы, где движение определяется окружающим потоком, а возможности исполнительных механизмов ограничены

Значимость

1. Вычислительные вызовы: Прямое применение высокоточных моделей Навье-Стокса для оптимизации MPC вычислительно недопустимо и не удовлетворяет требованиям управления в реальном времени
2. Практические потребности: Во многих реальных системах (например, высотные аэростаты) горизонтальное движение в основном определяется окружающим ветровым полем, с ограниченным управлением только в вертикальном направлении
3. Прогностическая ценность: Точное предсказание эволюции потока и использование структуры потока (например, вертикального ветрового сдвига) может значительно повысить эффективность навигации

Ограничения существующих методов

1. Традиционное MPC: При применении к системам, ограниченным PDE, вычислительные затраты чрезмерно высоки
2. Методы глубокого обучения: Такие как Balloon Learning Environment (BLE), не обладают физической согласованностью и интерпретируемостью
3. Существующие ROM-MPC: Сосредоточены главным образом на стабилизации в замкнутом контуре, меньше внимания уделяют генерированию пути в потоках

Исследовательская мотивация

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

Основные вклады

1. Интегрированная схема: Впервые объединены редуцированная модель POD-Галёркина и MPC для задачи планирования пути в потоке Навье-Стокса
2. Проектирование нестационарного наблюдателя: Предложен расширенный фильтр Калмана (EKF), объединяющий измерения от сети неподвижных датчиков и мобильного агента, с реализацией оценки состояния потока через нестационарную матрицу измерений C(t)
3. Генерирование траектории с учётом потока: Разработана задача оптимизации на конечном горизонте, использующая ROM для предсказания поля скоростей, с использованием вертикального управляющего воздействия для эксплуатации слоистой структуры потока
4. Практическая верификация: Использованы данные переанализа ERA5 для моделирования, подтверждающего эффективность и вычислительную производительность метода в реальных метеорологических условиях

Подробное описание метода

Определение задачи

Входные данные:

• Исторические снимки поля потока {uk(x)}^m_
• Измерения от сети неподвижных датчиков y_fixed(t)
• Измерения текущего положения и скорости агента y_mobile(t)
• Целевое положение x_ref

Выходные данные:

• Последовательность управления вертикальной скоростью {u_{z,k+i|k}}^{N-1}_
• Предсказанная траектория {x_{k+i|k}}^N_

Ограничения:

• Пространственные границы x ∈ X ⊂ R³
• Ограничения исполнительных механизмов u_z ∈ U (например, максимальная скорость набора высоты |u_z| ≤ 1 м/с)
• Терминальные ограничения x_{k+N|k} ∈ X_f

Архитектура модели

1. Редуцированное моделирование POD

Аппроксимация поля скоростей:

v(x,t) = Σ(i=1 to n) a_i(t)φ_i(x)

где φ_i(x) — моды POD, a_i(t) — нестационарные коэффициенты мод.

Вычисление мод POD: Получены путём решения задачи на собственные значения:

Uc = λc, U_ij = (1/m)⟨u_i, u_j⟩

Моды упорядочены по энергии λ_ℓ в убывающем порядке; первые n мод захватывают основную кинетическую энергию.

Динамика проекции Галёркина: Разложение поля скоростей подставляется в уравнения Навье-Стокса и проецируется на пространство мод POD:

ȧ_k(t) = L_k a(t) + a(t)^T Q_k a(t)

где:

• L_k: линейный член, элементы L_ = ⟨∇²φ_i(x), φ_k(x)⟩
• Q_k: квадратичный нелинейный член, элементы Q_ = ⟨(φ_j(x)·∇)φ_i(x), φ_k(x)⟩

2. Наблюдатель расширенного фильтра Калмана

Модель в пространстве состояний:

ȧ_k(t) = L_k a(t) + a(t)^T Q_k a(t)  (уравнение состояния)
y(t) = C(t)a(t)                        (уравнение наблюдения)

Нестационарная матрица измерений:

C(t) = [C_fixed; C_mobile(t)]

• C_fixed: оценка мод POD в положениях неподвижных датчиков
• C_mobile(t): оценка мод POD в текущем положении агента

Дискретный EKF:

Этап предсказания:

â_{k|k-1} = f_d(â_{k-1})
P_{k|k-1} = A_{k-1}P_{k-1|k-1}A^T_{k-1} + Q

Этап обновления:

K_k = P_{k|k-1}C^T_k(C_k P_{k|k-1}C^T_k + R)^{-1}
â_{k|k} = â_{k|k-1} + K_k(y_k - C_k â_{k|k-1})
P_{k|k} = (I - K_k C_k)P_{k|k-1}

3. Модельное предсказательное планирование пути

Кинематика агента:

ẋ(t) = v_pred(x(t), t) + u_z(t)e_z

Горизонтальное движение полностью определяется потоком, вертикальное направление управляемо.

Задача оптимизации на конечном горизонте:

min_{u_z} V_f(x_{k+N|k}) + Σ(i=0 to N-1) ℓ(x_{k+i|k}, u_{z,k+i|k})

s.t. x_{k+i+1|k} = x_{k+i|k} + Δt[v_pred(x_{k+i|k}, t_{k+i}) + u_{z,k+i|k}e_z]
     x_{k+i|k} ∈ X, u_{z,k+i|k} ∈ U, x_{k+N|k} ∈ X_f

Функция стоимости:

ℓ(x, u_z) = w_p ||x - x_ref||²_2 + w_u u²_z

Балансирует отслеживание цели и стоимость управления.

Технические инновации

1. Физическая согласованность: Проекция Галёркина сохраняет физическую структуру уравнений Навье-Стокса, в отличие от чистых моделей чёрного ящика на основе данных
2. Мобильное зондирование: Агент сам служит мобильным датчиком, непрерывно улучшая оценку потока через нестационарную матрицу измерений C(t)
3. Использование потока: Оптимизатор активно использует ветровой сдвиг путём регулирования высоты, а не пассивно реагирует на окружающую среду
4. Стратегия скользящего окна: На каждом шаге выполняется только первое управляющее воздействие, затем план пересчитывается на основе обновленного предсказания потока, обеспечивая робастность к неопределённости предсказания
5. Вычислительная эффективность: Редуцированная модель снижает высокомерную PDE к n-мерной системе ODE (в данной работе n значительно меньше размерности исходной сетки), делая возможной оптимизацию в реальном времени

Экспериментальная установка

Набор данных

Данные переанализа ERA5:

• Пространственная область: Центр в 49°N, 81°W, диапазон 2° по долготе и широте (примерно 146 км × 220 км)
• Вертикальный диапазон: Уровни давления от 125 гПа до 10 гПа (примерно 15-30 км высоты)
• Временная дискретизация: 36 дней, случайно выбранных в течение 6 месяцев, равномерная выборка каждый час
• Количество снимков: m снимков поля скоростей для разложения POD
• Предварительная обработка: Удаление среднего смещения, позволяющее ROM захватить динамику ветрового поля

Метрики оценки

1. Финальное расстояние до цели:

   d_f = ||x(T_f) - x_ref||_2
   

2. Средняя степень выравнивания с потоком:

   γ̄ = (1/T_f)∫[ẋ(t)·v_pred(x,t)/(||ẋ(t)|| ||v_pred(x,t)||)]dt
   

   
   Измеряет эффективность использования благоприятных направлений потока планировщиком
3. Стоимость управления:

   u_rms = √[(1/T_f)∫u²_z(t)dt]
   

4. Время выполнения (RT): Измеряется общее время вычисления с помощью MATLAB tic-toc

Методы сравнения

Основное сравнение проводится по различным длинам горизонта предсказания:

• T = 1.5h (короткий горизонт)
• T = 3h (средне-короткий горизонт)
• T = 6h (средне-длинный горизонт)
• T = 12h (длинный горизонт)

Детали реализации

Вычислительная платформа:

• MATLAB R2024a
• AMD Ryzen 2600X (6 ядер, 3.6 ГГц)
• 32 ГБ ОЗУ

Оптимизатор:

• Последовательное квадратичное программирование (SQP)
• Стратегия горячего старта (инициализация входной последовательностью из предыдущего шага)

Параметры моделирования:

• Общее время моделирования: 24 часа
• Интервал дискретизации: Δt = 10 минут
• Целевая область: цилиндрическая область радиусом 50 км
• Максимальная скорость набора высоты: u_max = 1 м/с
• Количество датчиков: 8 (включая агента)

Веса стоимости:

• w_p: вес отклонения положения
• w_u: вес стоимости управления (Подобраны эмпирически для балансирования производительности)

Результаты экспериментов

Основные результаты

Количественное сравнение производительности (Таблица I):

Ключевые находки:

1. Значительное влияние горизонта: Увеличение с 1.5ч до 6ч снижает финальное расстояние с 54.6 км до 13.5 км (улучшение на 75%)
2. Убывающая предельная полезность: Увеличение с 6ч до 12ч даёт улучшение всего на 0.1 км, но время вычисления увеличивается почти в 3 раза
3. Высокая степень выравнивания с потоком: Все конфигурации имеют γ̄ ≥ 0.92, что указывает на эффективное использование потока планировщиком
4. Повышение эффективности управления: Долгосрочное планирование позволяет более плавные регулировки высоты (u_rms снижается с 0.60 до 0.36 м/с)

Точность реконструкции поля потока

На рисунке 3 показано:

• Изменение RMSE предсказания ROM и фактического поля скоростей во времени
• На 15-часовой отметке ошибка предсказания остаётся в приемлемом диапазоне
• Тепловые карты на рисунках 1 и 2 демонстрируют высокое соответствие пространственной структуры предсказанного и реального полей

Анализ траектории

Горизонт предсказания 3 часа (рисунки 4 и 6):

• Агент начинает в области конвекции запад-восток
• Частые регулировки высоты для поиска благоприятных ветровых слоёв
• Финальное расстояние до цели 32.4 км
• Управляющие воздействия относительно активные (u_rms = 0.56 м/с)

Горизонт предсказания 12 часов (рисунки 5 и 7):

• Траектория более плавная и стратегическая
• Избегаются выборы с краткосрочной выгодой, но долгосрочным убытком
• Изменения высоты более постепенные и планомерные
• Финальное расстояние 13.4 км, стоимость управления снижена на 40%

Экспериментальные находки

1. Компромисс горизонта предсказания:
   • Короткий горизонт: локально осуществим, но глобально неэффективен, производит "близорукое" поведение
   • Длинный горизонт: глобально оптимален, но вычислительно дорог, чувствителен к неопределённости долгосрочного предсказания
   • Оптимальный горизонт: 6-12 часов достигает баланса между производительностью и вычислительными затратами
2. Стратегия использования потока:
   • Планировщик учится "ждать": медленное движение в неблагоприятных ветровых слоях в ожидании благоприятных условий
   • Активное использование вертикального ветрового сдвига: регулирование высоты для входа в ветровые слои благоприятного направления
3. Вычислительная осуществимость:
   • Даже при 12-часовом горизонте общее время вычисления 1983 секунды остаётся в приемлемом диапазоне для 24-часового моделирования
   • В практических приложениях можно добиться дальнейшего ускорения параллелизацией и специализированными решателями
4. Робастность:
   • Переплан в скользящем окне обеспечивает естественную робастность к ошибкам предсказания потока
   • EKF непрерывно обновляет оценку потока, снижая дрейф предсказания

Связанные работы

Модельное предсказательное управление

• Управление процессами: Климат-контроль зданий 1,2
• Управление потоком: Управление следом цилиндра 3,4
• Оптимизация, ограниченная PDE: Обзоры 5-7 демонстрируют применение ROM-MPC к крупномасштабным системам

Редуцированное моделирование

• Теория POD: Первоначально введена Lumley 11, формализована Berkooz и др. 12
• Приложения в гидродинамике: Сбалансированное POD Rowley 9, ROM для управления потоком Noack и др. 10
• Комплексные обзоры: Yano 8 и Taira и др. 14 предоставляют всеобъемлющие обзоры анализа мод в аэродинамике и гидродинамике

Методы, управляемые данными

• DMD и оператор Купмана: Proctor и др. 16, Korda и Mezić 15
• Глубокое обучение: Balloon Learning Environment 17,18 для навигации стратосферных аэростатов

Преимущества данной работы

1. Управляемость физической информацией: По сравнению с чистыми методами обучения, сохраняет физическую согласованность уравнений Навье-Стокса
2. Осуществимость в реальном времени: По сравнению с высокоточным CFD-MPC, вычислительная эффективность повышена на несколько порядков величины
3. Фокус на генерирование пути: Отличается от традиционного ROM-MPC, ориентированного на стабилизацию, сосредоточен на навигационной производительности

Выводы и обсуждение

Основные выводы

1. Эффективность схемы: Объединение POD-Галёркина ROM и MPC успешно реализует эффективное планирование пути в потоке Навье-Стокса
2. Возможность управления в реальном времени: Редуцированное моделирование снижает вычислительную сложность с высокомерной PDE до низкомерной ODE, делая оптимизацию в скользящем окне осуществимой в практических временных масштабах
3. Физическая согласованность: Проекция Галёркина сохраняет основные структуры гидродинамики, производя физически обоснованные траектории
4. Практическая верификация: Моделирование с использованием реальных данных ERA5 демонстрирует применимость метода в сложных метеорологических условиях

Ограничения

1. Предположения модели:
   • Предполагается, что агент — материальная точка, игнорируется динамика ориентации и аэродинамические характеристики реальных летательных аппаратов
   • Рассматривается только вертикальное управление, не моделируется горизонтальная тяга
2. Ограничения ROM:
   • Моды POD основаны на исторических данных, могут неточно предсказывать экстремальные потоки вне обучающего набора
   • Предположение о линейном подпространстве может не захватить высоконелинейные или переходные явления
3. Вычислительные затраты:
   • Длинный горизонт (12ч) всё ещё требует значительного времени вычисления (~2000 секунд)
   • Отсутствует прямое сравнение вычислительных затрат с высокоточным CFD или другими методами
4. Гарантии осуществимости:
   • Поскольку горизонтальное движение полностью определяется потоком, невозможно гарантировать строгую рекурсивную осуществимость
   • Отсутствуют теоретические доказательства стабильности или сходимости
5. Объём экспериментов:
   • Тестирование проведено только в одном географическом регионе и сезоне
   • Не рассмотрены шумы измерений, отказы датчиков и другие практические факторы

Направления будущих исследований

1. Интеграция динамики:
   • Включение детальной модели динамики летательного аппарата (ориентация, аэродинамическое сопротивление)
   • Рассмотрение многовходового управления (горизонтальная тяга + вертикальное управление)
2. Гибридная схема обучения:
   • Интеграция с Balloon Learning Environment 17,18
   • Объединение преимуществ физического ROM и глубокого обучения
3. Количественная оценка неопределённости:
   • Интеграция оценок неопределённости предсказания ROM
   • Робастное MPC или стохастическое MPC формулирование
4. Адаптивное редуцирование:
   • Онлайн обновление мод POD для адаптации к новым режимам потока
   • Динамическая регулировка количества сохраняемых мод
5. Многоагентная координация:
   • Расширение на координированное планирование пути для нескольких агентов
   • Использование распределённого зондирования группой агентов

Глубокая оценка

Преимущества

1. Новизна метода:

• Впервые систематически объединены POD-Галёркина ROM и MPC для планирования пути в потоке Навье-Стокса
• Проектирование нестационарной матрицы измерений остроумно, интегрирует мобильный агент как датчик в наблюдатель
• Органичное слияние физической информации и методов, управляемых данными

2. Теоретическая строгость:

• Полная цепь вывода, начиная с уравнений Навье-Стокса, ясна и логична
• Математическое выражение проекции Галёркина строго
• Теоретическая основа проектирования наблюдателя EKF прочна

3. Достаточность экспериментов:

• Использование реальных данных ERA5 повышает достоверность
• Систематическое сравнение по нескольким горизонтам предсказания
• Полный набор количественных метрик (расстояние, выравнивание, стоимость управления, время вычисления)

4. Практическая ценность:

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

5. Ясность изложения:

• Структура ясна, логика последовательна
• Математические обозначения определены чётко
• Графики и таблицы эффективно поддерживают объяснение

Недостатки

1. Ограничения метода:

• Чрезмерные упрощающие предположения: Модель материальной точки игнорирует сложную динамику реальных летательных аппаратов
• Однонаправленное управление: Предположение об управлении только в вертикальном направлении ограничивает применимость
• Отсутствие теоретических гарантий: Нет доказательств рекурсивной осуществимости или стабильности

2. Недостатки экспериментального проектирования:

• Отсутствие базовых методов сравнения: Нет прямого сравнения с другими методами (RL из BLE, традиционное MPC, эвристические методы)
• Единственный сценарий: Тестирование только в одном географическом регионе и сезоне, неизвестна обобщаемость
• Идеализированные условия: Не рассмотрены шумы измерений, задержки связи, динамика исполнительных механизмов и другие практические факторы

3. Недостаточная глубина анализа:

• Влияние точности ROM: Отсутствует систематический анализ влияния количества мод n на производительность планирования
• Отсутствие анализа чувствительности: Выбор весов стоимости w_p, w_u лишён принципиального руководства
• Отсутствие случаев отказа: Не показаны ситуации, в которых метод неэффективен, или граничные условия

4. Оценка вычислительной эффективности:

• Ограниченное значение абсолютного времени: Отсутствует сравнение с высокоточными методами
• Неизвестная масштабируемость: Производительность при большей пространственной области или более высокой разрешающей способности не тестирована
• Сомнительная реальная временность: Время вычисления 2000 секунд для 12-часового горизонта может быть слишком велико для некоторых приложений

5. Проблемы воспроизводимости:

• Код не открыт
• Некоторые гиперпараметры (количество мод POD n, ковариации EKF Q и R) не указаны явно
• Детали предварительной обработки данных (конкретный метод удаления среднего) недостаточно подробны

Влияние

Вклад в область:

• Теоретический вклад: Предоставляет новую теоретическую схему для планирования пути в системах, управляемых потоком
• Методологический вклад: Демонстрирует практическое применение ROM-MPC в сложных PDE-ограниченных системах
• Прикладной вклад: Предоставляет осуществимое решение управления для стратосферных аэростатов и других реальных систем

Практическая ценность:

• Средняя-высокая: Метод демонстрирует хорошую производительность в идеальных условиях, но требует дальнейшей инженерной разработки
• Уровень готовности технологии: Примерно TRL 3-4 (проверка концепции), требуется значительная работа для практического развёртывания

Воспроизводимость:

• Средняя: Теоретическое описание ясно, но отсутствуют код и детальные параметры, полное воспроизведение требует определённых усилий

Ожидаемое влияние:

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

Применимые сценарии

Идеально применимые сценарии:

1. Навигация стратосферных аэростатов: Исходный мотивирующий сценарий, горизонтальное движение определяется ветром
2. Морские планеры: Аналогичные характеристики вертикального управления + потокового привода
3. Беспилотные летательные аппараты для атмосферного зондирования: Требуют энергоэффективного планирования пути с использованием ветра

Сценарии, требующие улучшения:

1. Быстрые маневренные летательные аппараты: Требуют включения полной динамики и многовходового управления
2. Высокопереходные потоки: Такие как турбулентность и конвекция, где линейное подпространство POD может быть неадекватным
3. Задачи с множественными ограничениями: Такие как избежание препятствий и множественные цели, требуют расширения формулирования оптимизации

Неприменимые сценарии:

1. Наземные мобильные роботы: Влияние потока минимально, преимущества метода не проявляются
2. Задачи высокоточного удержания позиции: Модель материальной точки и однонаправленное управление недостаточны
3. Приложения с экстремальными требованиями к реальной временности: Текущее время вычисления может быть слишком велико

Ключевые ссылки

Фундаментальные ссылки:

1. Теория POD:
   • Lumley (1967): Первоначальное введение POD в анализ турбулентности
   • Berkooz, Holmes, Lumley (1993): Формализация POD в гидродинамике
2. Редуцированное моделирование:
   • Rowley (2005): Сбалансированное POD для моделирования потоков
   • Noack et al. (2011): Монография по редуцированному моделированию для управления потоком
3. Приложения ROM-MPC:
   • Benner, Gugercin, Willcox (2015): Обзор методов редуцирования параметризованных динамических систем
   • Yano (2020): Редуцирование моделей в вычислительной аэродинамике
4. Практические приложения:
   • Greaves et al. (2021): Balloon Learning Environment
   • Oldewurtel et al. (2012): Применение MPC в управлении энергией зданий

Общая оценка: Это твёрдая методологическая статья, которая успешно объединяет классическое редуцированное моделирование POD с современным управлением MPC, решая проблему планирования пути в реальном времени в системах, управляемых потоком. Теоретический вывод строг, экспериментальная верификация достаточна, имеется явная практическая ценность. Основные недостатки заключаются в сильных упрощающих предположениях, отсутствии сравнения с другими методами и недостаточной инженерной разработке для практического развёртывания. Статья вносит позитивный вклад в пересечение теории управления и гидродинамики, предоставляя хорошую отправную точку для последующих исследований. Рекомендуется, чтобы будущая работа сосредоточилась на интеграции динамики, количественной оценке неопределённости и верификации на реальных системах.