Квантовые компьютеры ускоряют LS-DYNA: первые квантовые вычисления, полезные для прикланого CAE
После прошлой недели я перестал считать, что квантовые компьютеры - это просто академический распил бабла. Итак, мы тут все таки или иначе решаем задачи линейной алгебры. Причем, в механике, когда мы используем неявный решатель, у нас собирается большая и страшная матрица жесткости, работать с которой не очень приятно, и все хотят ускорить данный процесс.
Казалось бы, LS-DYNA знаменита своим явным решателем, где такой проблемы просто нет. Но явный решатель не всегда подходит для задачи медленной динамики. А если мы вспоминаем про связанные задачи, где у нас есть FEM и CFD c FSI (привет медицинские расчеты) или про линейную динамику (привет NVH), то на одном явном решателе не выехать.
О чем же статья от ученых из Ansys и IonQ Так, препритн которой для журнала Quantum Physics был опубликована 17 марта 2025.
Основная проблема: Решение больших систем линейных уравнений Ax=b, возникающих в задачах МКЭ, представляет собой вычислительно трудоёмкий процесс. Для таких систем характерны большие разреженные матрицы. Однако прямые методы решения (например, LU- и LDL^T-разложения) приводят к существенному заполнению («fill-in») матриц, что резко увеличивает затраты времени на решение.
Цель исследования: Снизить заполнение матрицы (минимизация «fill-in») путём оптимального переупорядочивания строк и столбцов матрицы. Для решения задачи переупорядочивания используется квантовый подход с применением алгоритма VarQITE.
Квантовый подход VarQITE: Алгоритм VarQITE основан на концепции эволюции квантового состояния по мнимому времени и применяется для поиска основного состояния гамильтониана, который описывает задачу оптимизации. Задача переупорядочивания матрицы формулируется как задача разделения графа (Graph Partitioning Problem, GPP), которая сводится к задаче бинарной оптимизации (QUBO — Quadratic Unconstrained Binary Optimization). Далее эта задача кодируется в виде гамильтониана, минимальное собственное значение которого соответствует оптимальному решению.
Интеграция с LS-DYNA: Предложен гибридный квантово-классический подход, интегрированный в программный пакет LS-DYNA от Ansys. LS-DYNA применяется для механических задач, задач вибрационного анализа и гидродинамики. В LS-DYNA используется метод вложенного разбиения (Nested Dissection) и многоуровневые алгоритмы для грубого разбиения графа на подграфы меньшего размера. Эти подграфы затем решаются с помощью VarQITE на квантовом симуляторе или квантовом устройстве.
Результаты исследования: Экспериментально показано, что VarQITE может улучшить качество решения задачи разбиения графов в сравнении с классическими эвристиками (например, LS-GPart, встроенным в LS-DYNA). Результаты тестирования включали задачи из области механики (например, моделирование деформации крыши автомобиля), задачи течения жидкости (насосы крови) и вибрационный анализ автомобиля. В реальных тестах (с миллионами узлов и десятками миллионов рёбер) время решения системы при использовании квантового подхода сокращалось до 12% по сравнению с классическими методами.
Реализация на реальном квантовом оборудовании: VarQITE был протестирован на реальных квантовых компьютерах IonQ (Aria и Forte), и результаты были сопоставимы с идеализированными симуляциями. Для улучшения качества полученных решений авторы предлагают гибридный подход: решение задачи на квантовом компьютере с последующим уточнением с помощью классического алгоритма (например, модифицированный алгоритм Фидуччи-Маттиеса).
Заключение: Полученные результаты демонстрируют потенциал использования квантовых вычислений в реальных задачах инженерного анализа. VarQITE, интегрированный в промышленные продукты, такие как LS-DYNA, способен уже в ближайшем будущем предоставить вычислительные преимущества.
Препринт: https://arxiv.org/abs/2503.13128
Accelerating large-scale linear algebra using variational quantum…