Все выпуски

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

В прямоугольной области исследуются нелокальные краевые задачи для одномерного нестационарного уравнения конвекции-диффузии дробного порядка с переменными коэффициентами, описывающие диффузионный перенос той или иной субстанции, а также перенос, обусловленный движением среды. Методом энергетических неравенств выводятся априорные оценки решений нелокальных краевых задач в дифференциальной форме. Построены разностные схемы, и для них доказываются аналоги априорных оценок в разностной форме, приводятся оценки погрешности в предположении достаточной гладкости решений уравнений. Из полученных априорных оценок следуют единственность и устойчивость решения по начальным данным и правой части, а также сходимость решения разностной задачи к решению соответствующей дифференциальной задачи со скоростью $O(h^2+\tau^2)$.

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

В работе рассматривается задача оптимального управления одномерным процессом, заданным стохастическим дифференциальным уравнением, в котором управление воздействует как на коэффициент сноса, так и на коэффициент диффузии, при этом диффузионная составляющая линейна по управлению $$dx(t) = b(t,x(t),u(t))dt +\sigma(t,x(t))u(t)dW(t),\qquad x(0) = x_0.$$ Здесь $x(t)$ - фазовая координата, $u(t)$ - управляющая функция, $W(t)$ - винеровский процесс. Доказана теорема, которая предоставляет структуру решения рассматриваемого уравнения в виде суперпозиции функций $x(t)=Φ(t,u(t)W(t)+y(t))$, в котором $Φ(t,v)$ - известная функция, полностью определяющаяся коэффициентом $σ(t,x)$, и не зависит от управления, а $y(t)$ - решение потраекторно-детерминированного дифференциального уравнения с мерой вида

$$dy(t) = B(t,y(t),u(t))dt - W(t)du(t).$$

Выявленная структура решения позволяет вместо исходной стохастической задачи оптимального управления исследовать новую эквивалентную задачу с фазовой переменной $y(t)$, которая является потраекторно-детерминированной задачей оптимального импульсного управления. При детерминированном рассмотрении новой задачи решения последней могут оказаться упреждающими функциями, поэтому в работе предлагается метод, который позволяет добиться неупреждаемости оптимальных решений. Суть метода заключается в модификации функционала потерь в новой потраекторно-детерминированной задаче специальным образом подобранным интегральным слагаемым, которое позволяет гарантировать неупреждаемость решений.

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

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

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

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

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

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