RUS  ENG
Полная версия
ЖУРНАЛЫ // Журнал вычислительной математики и математической физики // Архив

Ж. вычисл. матем. и матем. физ., 2022, том 62, номер 4, страницы 553–563 (Mi zvmmf11381)

Общие численные методы

Общая структура редуцированных базисов для задач агрегационной кинетики различной размерности

С. А. Матвеевab, А. П. Смирновa, И. В. Тимохинabc, Е. Е. Тыртышниковabc

a 119991 Москва, Ленинские горы, стр. 52, МГУ им. М.В. Ломоносова, факультет ВМК, Россия
b 119333 Москва, ул. Губкина, 8, Институт вычислительной математики им. Г.И. Марчука Российской академии наук, Россия
c 119991 Москва, Ленинские горы, МГУ, Московский центр фундаментальной и прикладной математики, Россия

Аннотация: Данная работа посвящена исследованию структуры линейных пространств, позволяющих приближать решения систем дифференциальных уравнений, возникающих в задачах агрегационной кинетики с источником и стоком частиц, на больших интервалах времени с высокой точностью. Основным параметром для семейства рассматриваемых моделей агрегационной кинетики является их размерность $N$, соответствующая максимальному допустимому размеру частицы в системе. При больших значениях предельного допустимого размера частицы $N$ исследование таких задач оказывается затруднительным, так как начинает требовать значительных вычислительных ресурсов. Для ускорения расчетов мы используем идеологию редукции размерности при помощи поиска базиса искомого маломерного пространства методом снимков. В результате использования такого подхода можно установить существование такого пространства, построить его базис и свести исходную нелинейную задачу высокой размерности к задаче существенно меньшей размерности, позволяющей организовать вычисления более эффективно без существенной потери качества численных решений. Главным результатом данной работы является найденная возможность использования базиса задач большей размерности $N$ для решения задач меньшей размерности $M<N$ без подготовки и построения базиса для размерности $M$. В численных экспериментах мы демонстрируем, что достаточно просто взять первые $M<N$ компонент векторов базиса задачи размерности $N$ и использовать полученную систему в качестве базиса для новой задачи размерности $M$. С использованием такого базиса нам удается получать превосходную точность численных решений редуцированной задачи размерности $M$, хотя сами по себе решения задач размерности $M$ и $N$ не обязаны быть близки. Дополнительно в работе предлагается несколько оценок для норм решений уравнений необратимой агрегации с источником. Из оценок следует, что прямая подстановка векторов базиса $N$ в задачи размерности $M>N$ не приводит к успеху, а расширение базиса для задач большей размерности требует дополнительных усилий, что согласуется с результатами экспериментов.
Библ. 28. Фиг. 3.

Ключевые слова: уравнение Смолуховского, редукция модели, агрегация, метод снимков.

УДК: 517.983

Поступила в редакцию: 04.10.2021
Исправленный вариант: 04.10.2021
Принята в печать: 16.12.2021

DOI: 10.31857/S004446692204010X


 Англоязычная версия: Computational Mathematics and Mathematical Physics, 2022, 62:4, 538–547

Реферативные базы данных:


© МИАН, 2024