Общие численные методы
Общая структура редуцированных базисов для задач агрегационной кинетики различной размерности
С. А. Матвеев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