Версия 2.0 -- cashocs: Вычислительное программное обеспечение для оптимизации формы и оптимального управления на основе сопряженных методов

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

• ID статьи: 2306.09828
• Название: Version 2.0 -- cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software
• Автор: Sebastian Blauth (Fraunhofer ITWM, Кайзерслаутерн, Германия)
• Классификация: math.OC (Оптимизация и управление)
• Время публикации: Подано на arXiv 13 ноября 2023 г., опубликовано в SoftwareX 24
• Ссылка на статью: https://arxiv.org/abs/2306.09828
• Версия программного обеспечения: v2.0.12
• Лицензия: GNU GPL v3.0

Аннотация

В данной статье представлено обновление версии 2.0 программного обеспечения cashocs. cashocs — это вычислительное программное обеспечение на основе сопряженных методов, специально разработанное для автоматизированного решения задач оптимизации, ограниченных уравнениями в частных производных (УЧП), включая оптимизацию формы и оптимальное управление. С момента первого выпуска программное обеспечение получило множество новых функций и практических инструментов, что сделало его более гибким и эффективным. Наиболее важные новые функции включают: фреймворк пространственного отображения, возможность решения задач топологической оптимизации на основе метода уровней, поддержку параллельных вычислений через MPI и возможность обработки дополнительных (государственных) ограничений.

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

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

Задачи оптимизации, ограниченные УЧП, широко распространены в инженерии и научных вычислениях и включают:

1. Задачи оптимизации формы: оптимизация производительности системы путем изменения геометрической формы
2. Задачи оптимального управления: оптимизация поведения системы путем управления входными сигналами
3. Задачи топологической оптимизации: оптимизация производительности конструкции путем изменения распределения материала

Значимость исследования

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

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

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

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

1. Фреймворк пространственного отображения: впервые реализована полная технология пространственного отображения в программном обеспечении для оптимизации, ограниченной УЧП, поддерживающая задачи оптимизации формы и оптимального управления
2. Поддержка топологической оптимизации: интегрирован решатель топологической оптимизации на основе метода уровней, включая новейшие квазиньютоновские методы
3. Возможности параллельных вычислений: реализована поддержка параллельных вычислений через MPI, позволяющая программному обеспечению работать на системах высокопроизводительных вычислений
4. Механизм обработки ограничений: добавлены методы квадратичного штрафа и расширенный метод множителей Лагранжа для обработки дополнительных ограничений на состояние и управление
5. Улучшение удобства использования: включает автоматическое масштабирование, пользовательские скалярные произведения, улучшенный рабочий процесс перестроения сетки и другие функции

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

Архитектура программного обеспечения

cashocs 2.0 основан на модульной архитектуре и поддерживает три основных класса задач оптимизации:

Задачи оптимизации, ограниченные УЧП
├── cashocs.ShapeOptimizationProblem
├── cashocs.OptimalControlProblem  
├── cashocs.TopologyOptimizationProblem
├── Дополнительные ограничения
│   ├── cashocs.ConstrainedShapeOptimizationProblem
│   └── cashocs.ConstrainedOptimalControlProblem
└── Задачи пространственного отображения
    ├── cashocs.space_mapping.shape_optimization.SpaceMappingProblem
    └── cashocs.space_mapping.optimal_control.SpaceMappingProblem

Технология пространственного отображения

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

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

Метод топологической оптимизации

Использует метод уровней для решения задач топологической оптимизации:

1. Поддерживаемые алгоритмы: включает традиционные алгоритмы топологической оптимизации и новый квазиньютоновский метод
2. Автоматическая система сопряженных уравнений: хотя топологическая производная должна быть предоставлена пользователем, автоматическое вывод системы сопряженных уравнений уже реализован
3. Расширение области применения: расширяет применимость cashocs на область топологической оптимизации с информацией о топологической чувствительности

Механизм обработки ограничений

Реализованы два основных метода обработки ограничений:

1. Метод квадратичного штрафа: обработка ограничений путем добавления штрафных членов в целевую функцию
2. Расширенный метод множителей Лагранжа: гибридный метод, сочетающий множители Лагранжа и штрафные функции
3. Широкая применимость: поддерживает различные типы ограничений, включая ограничения на состояние и управление

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

Технологический стек

• Язык программирования: Python
• Основные зависимости: FEniCS, NumPy, PETSc, MPI
• Обработка сетки: meshio, Gmsh
• Контроль версий: git
• Документация: https://cashocs.readthedocs.io/

Примеры постановки задач

Случай оптимизации формы с пространственным отображением:

• Цель: достижение равномерного распределения потока на трех выходных трубопроводах
• Детальная модель: уравнения Навье-Стокса при числе Рейнольдса 1000, решаемые с помощью Ansys Fluent
• Грубая модель: линейная система Стокса, решаемая в FEniCS на грубой сетке

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

Производительность метода пространственного отображения

Из результатов на рисунке 1 видно:

• Скорость сходимости: сходимость достигается примерно за 5 итераций
• Эффект оптимизации: окончательная геометрическая форма обеспечивает почти идеальное равномерное распределение потока
• Вычислительная эффективность: значительное повышение эффективности по сравнению с прямой оптимизацией детальной модели

Производительность параллельных вычислений

• Без изменения кода: большинство последовательных задач оптимизации, ограниченных УЧП, могут быть непосредственно распараллелены
• Высокопроизводительные вычисления: поддержка решения крупномасштабных задач оптимизации на системах HPC
• На основе FEniCS: полное использование поддержки MPI в FEniCS с минимальными изменениями кода на стороне пользователя

Улучшение перестроения сетки

• Расширенная применимость: новый рабочий процесс перестроения сетки поддерживает задачи пространственного отображения и обработки ограничений
• Обновление синтаксиса: пользователи должны следовать немного отличающемуся, но все еще интуитивному синтаксису
• Узкое место производительности: перестроение сетки само по себе не распараллелено из-за ограничений Gmsh и meshio

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

Развитие технологии пространственного отображения

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

Алгоритмы топологической оптимизации

• Интегрированы новейшие алгоритмы топологической оптимизации, предложенные Amstutz и Andrä
• Реализован квазиньютоновский метод топологической оптимизации, предложенный Blauth и Sturm
• Остается место для улучшения в автоматическом вычислении информации о топологической чувствительности

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

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

1. Полнота функциональности: cashocs 2.0 значительно расширил функциональность программного обеспечения, став комплексным решением для оптимизации, ограниченной УЧП
2. Ценность промышленного применения: фреймворк пространственного отображения особенно подходит для промышленных приложений и может эффективно использовать существующие коммерческие решатели
3. Повышение вычислительных возможностей: поддержка параллельных вычислений MPI позволяет программному обеспечению решать крупномасштабные практические задачи
4. Улучшение удобства использования: множество улучшений повышают удобство и гибкость использования программного обеспечения

Ограничения

1. Вычисление топологической производной: топологическая производная в топологической оптимизации по-прежнему требует ручного ввода пользователем
2. Узкое место перестроения сетки: процесс перестроения сетки не распараллелен и может стать узким местом производительности для некоторых приложений
3. Зависимость от инструментов: зависимость от ограничений внешних инструментов, таких как Gmsh и meshio

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

1. Автоматическое вычисление топологической производной: разработка методов автоматического вычисления топологической производной
2. Параллельное перестроение сетки: улучшение параллельной производительности рабочего процесса перестроения сетки
3. Дополнительные алгоритмы оптимизации: интеграция дополнительных передовых алгоритмов оптимизации и методов обработки ограничений

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

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

1. Высокая инновационность: впервые в одном программном обеспечении интегрированы пространственное отображение, топологическая оптимизация и параллельные вычисления
2. Высокая практическая ценность: особенно подходит для промышленных приложений и может интегрироваться с существующим коммерческим программным обеспечением
3. Передовые технологии: использует новейшие методы численной оптимизации и алгоритмы
4. Полная документация: предоставляет подробные учебные материалы и документацию
5. Открытый исходный код: лицензия GNU GPL обеспечивает доступность программного обеспечения

Недостатки

1. Зависимость от платформы: высокая зависимость от экосистемы FEniCS
2. Кривая обучения: для непрофессиональных пользователей овладение концепциями оптимизации, ограниченной УЧП, требует определенных усилий
3. Узкие места производительности: сериализация перестроения сетки ограничивает масштабируемость некоторых приложений

Влияние

1. Академический вклад: предоставляет важный инструмент для исследований в области оптимизации, ограниченной УЧП
2. Промышленная ценность: может значительно повысить эффективность решения промышленных задач оптимизации
3. Развитие сообщества: открытый характер способствует формированию активного сообщества пользователей и разработчиков

Сценарии применения

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

Библиография

Ключевые ссылки включают:

1. Blauth, S. (2021). cashocs: A Computational, Adjoint-Based Shape Optimization and Optimal Control Software. SoftwareX, 13, 100646.
2. Blauth, S. (2023). Space Mapping for PDE Constrained Shape Optimization. SIAM J. Optim., 33(3), 1707–1733.
3. Amstutz, S., & Andrä, H. (2006). A new algorithm for topology optimization using a level-set method. J. Comput. Phys., 216(2), 573–588.

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