Все выпуски

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

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

В работе рассмотрен локально-неравновесный процесс затвердевания переохлажденного бинарного расплава. В целях простоты предполагается, что затвердевающая бинарная система находится при постоянных температуре и давлении и имеет две фазы, соответствующие твердому и жидкому состояниям. Математическое описание процесса затвердевания основано на модели фазового поля, обобщающей подход Плаппа (M. Plapp, Phys. Rev. E 84, 031601 (2011)) на случай локально-неравновесных процессов. Для вывода термодинамически согласованных уравнений модели использован метод расширенной необратимой термодинамики в отличие от феноменологического подхода Плаппа. Другое различие с моделью Плаппа состоит в использовании в качестве динамической переменной концентрации, а не химпотенциала примеси. В рамках полученной модели показана эквивалентность описания процесса затвердевания через концентрационное поле и через химпотенциал системы. В силу малости времен релаксации представленная модель сводится к сингулярно-возмущенной системе уравнений в частных производных параболического типа, описывающих динамику фазового и концентрационного полей. В работе предполагается известным описание термодинамических равновесных состояний на основе экспериментально полученных потенциалов Гиббса.

Для проверки полученной модели проведено численное моделирование одномерной задачи затвердевания в приближении разбавленного расплава Si-As, ранее неоднократно исследовавшегося экспериментально. Чтобы численно решить систему сингулярно-возмущенных уравнений, в работе предложен градиентно-устойчивый явный метод интегрирования уравнений второго порядка точности по времени. Для сведения бесконечного пространственного интервала к конечному использован метод «периодического сдвига». Оценка устойчивости получена из численных экспериментов.

Из численного моделирования процесса затвердевания разбавленного расплава Si-As получены профили концентрации и фазового поля, а также коэффициент распределения примеси на фронте затвердевания в зависимости от величины переохлаждения. Для проверки адекватности результатов численных экспериментов использовано аналитическое выражение для коэффициента распределения как функции переохлаждения, полученное из точного решения локально-неравновесной модели с резкой границей. Исследовано влияние параметров модели на процесс затвердевания и поведение численных решений вблизи диффузной границы.

Для нахождения решений гиперболического уравнения Аллена-Кана использован метод первого интеграла, который следует из известной теоремы Гильберта о нулях. Получены точные аналитические решения в виде бегущей волны, определяющие полный класс решений гиперболического уравнения Аллена-Кана. Показано, что в этом классе существует два подкласса решений: подкласс непрерывных решений и подкласс разрывных решений с сингулярностью в начале координат. Такая неединственность решений ставит вопрос об устойчивом аттракторе, то есть о решении бегущей волны, к которому будут стремиться нестационарные состояния, определяемые гиперболическим уравнением Аллена-Кана. Найденные решения включают в себя как частный случай полученные ранее решения для параболического уравнения Аллена-Кана в виде конечного числа $\tanh$-функций.

В работе исследуется алгоритм решения уравнения кристаллического фазового поля (Phase Field Crystal - PFC) в гиперболической постановке. Уравнение описывает фазовые превращения из метастабильного или неустойчивого состояния на масштабе атомной плотности и является дифференциальным уравнением шестого порядка по пространству и второго порядка по времени. Алгоритм основан на методе изогеометрического анализа (IGA) и реализован посредством библиотеки PetIGA. Полученный программный код допускает распараллеливание расчетов, что существенно ускоряет процесс решения задачи. Дана оценка эффективности используемых инструментов при проведении расчетов на высокопроизводительных вычислительных кластерах. Проведен анализ эффективности исследуемого алгоритма при работе с гетерогенными вычислительными системами.

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

При помощи модифицированного метода функции Ляпунова и построения омега-предельного множества соответствующей динамической системы сдвигов сформулированы утверждения о существовании ограниченных на положительной полуоси допустимых управляемых процессов и утверждение о равномерной локальной управляемости соответствующего магистрального процесса.