Котротких Андрей Сергеевич. Динамика концентраций, определяемая нелинейным уравнением "реакция-диффузия" и его обобщениями : Котроткіх Андрій Сергійович. Динаміка концентрацій, що визначається нелінійним рівнянням "реакція-дифузія" і його узагальненнями

Краткое описание

Котротких Андрей Сергеевич. Динамика концентраций, определяемая нелинейным уравнением "реакция-диффузия" и его обобщениями: диссертация ... кандидата Физико-математических наук: 01.01.01 / Котротких Андрей Сергеевич;[Место защиты: ФГБОУ ВО «Воронежский государственный университет»], 2018

Введение к работе
Актуальность темы.Тема диссертации находится на стыке двух направлений из списка «Основные научные направления ВГУ» (раздел «Наука» в портале ВГУ): 1. Аналитические, геометрические и численные методы исследования дифференциальных уравнений; 2. Теория функций и функциональный анализ.
Анализом бифуркационных эффектов начали заниматься еще в XIX веке, и к настоящему времени накопилось большое количество методик по их прогнозированию и «полезному использованию», появились многочисленные публикации и монографии. Однако потребность в развитии новых методов бифуркационного анализа, соответствующих новым запросам практики и современным достижениям вычислительных технологий, сохраняется до сих пор.
Сопровождающее бифуркацию изменение параметров внешнего воздействия (температуры, электромагнитного поля, механического сжатия и пр.) на сложную физическую систему (раствор, смесь, сплав и т.п.) в некоторых случаях приводит к потере устойчивости исходной фазы и, как следствие (как отклик системы), к се переходу в новое состояние (с новыми структурными свойствами). Такой переход сопровождается спинодальным расслоением (распадом), выраженным в изменении локальных концентраций компонентов, в образовании сначала зернистой структуры, а затем кластеров и доменов новой фазы. Структурную перестройку физической среды часто объясняют на основе нелинейных диффузионных уравнений Кана-Хилларда и Свифта-Хойенберга. Близким, но более простым уравнением, также способным моделировать структурные перестройки, является широко известное уравнение «реакция-диффузия» с кубической нелинейностью
w = A(w)+ Aw+w3-С, w = w(x), х eU сМ2,
рассмотренное при краевых условиях Неймана. Исследование посткритических структурных перестроек физических систем, моделируемых данным уравнением и его обобщениями, является весьма актуальной задачей, требующей для своего решения разнообразных методов современного математического анализа и новых вычислительных средств.
Степень разработанности темы.Бифуркационный анализ краевых и начально-краевых задач развивался в Воронежской математической школе, начиная с трудов М.А. Красносельского и его учеников — П.П. Забрейко, В.В. Стрыгина, Ю.Г. Борисовича, Ю.С. Колесова, Э.М.
Мухамадиева, Н.А. Бобылева и др.
Условия зарождения и развития пространственно однородных периодических режимов, описываемых начально-краевыми задачами для квазилинейных параболических уравнений изучались в ярославской школе динамических систем (в многочисленных трудах Ю.С. Колесова, А.С. Кащенко, С.Д. Глызина и других представителей этой школы). Для изучения условий зарождения периодических режимов и построения асимптотических представлений ветвей периодических решений были созданы специальнык процедуры нормализации уравненений, посредством которых определялись основные динамические характеристики бифурциру-ющих колебательных режимов. Фактически были разрабатаны методы инвариантных интегральных подмногообразий и обобщенных нормальных форм, с помощью которых анализ исходного уравнения сводится к изучению конечномерных нелинейных систем обыкновенных дифференциальных уравнений специального вида. Развитие предложенных конструкций опиралось на более ранние идеи, изложенные в известных трудах Хартмана, Митропольского, Лыкова, Бибикова, Брюно, Хэссарда, Казаринова, Вэна, Гукснхсймсра, Холмса и др. С помощью новых методов были получены новые результаты о существовании, устойчивости и асимптотических представлениях колебательных режимов в ситуациях с достаточно сложными вырождениями динамических систем.
В недавно опубликованной работе А.В. Казарникова и СВ. Ревиной1получены формулы асимптотических приближений к бифурцирующему из нуля периодическому решению обобщенной системы Релея с диффузией. Получние закритической ветви автоколебаний проведено на основе (невариационной) схемы Ляпунова-Шмидта, ранее предложенной В.И. Юдовичем.
Анализ многомодовых посткритических состояний включает, как известно, задачу вычисления значений собственных функций возмущенных самосопряженных операторов. В многочисленных трудах известных российских и зарубежных ученых созданы для решения этой задачи как общие, так и специальные методы. Важное место в арсенале таких средств занимает идея использования регуляризованных следов (В.А. Садовничий, В.В. Дубровский, СИ. Кадченко, СН. Какушкин и др.2). В боль-
1Возникновение автоколебаний в системе Рэлея с диффузией / А.В. Казарников, СВ. Ревина //Вестник Южно-Уральского государственного университета, 2016, т. 9, №2, с.16-28.
2Нахождение собственных значений и собственных функций методом регуляризованных следов/ СИ. Кадченко, СН. Какушкин // Челябинск: Изд. ЮУрГУ, 2015. 246 с.
шом цикле работ А.Г. Баскакова и его учеников для аналогичных задач был разработан метод подобных операторов3,4. Эти разработки имеют хорошую перспективу применения в многомодовом посткритическом анализе.
Значительные результаты были достигнуты школой Ю.И. Сапронова, усилиями которой построены теоретические и конструктивные схемы анализа многомодовых и нелокальных бифуркаций. Были рассмотрены также важные примеры использования новых исследовательских схем в теории упругости, теории фазовых переходов и гидродинамике.
Известно, что один из базовых принципов исследования бифуркаций решений начально краевых задач для параболических и более общих уравнений основан на том, что уравнение
— + Av = f(t}v)}0<t}</tv(0) = v0}
гдеf(t, x)при каждомtЄ [0, а]— нелинейный оператор (при условии, что операторАпорождает сильно непрерывную полугруппуT(t)),сводится к интегральному уравнению
v(t) = T(t)v0+ /T(t- s)f(s,v(s))ds J0
(метод Дюамеля).
В настоящей диссертации рассмотрен более простой подход, основанный на том, что рассмотренные бесконечномерные динамические системы являются градиентными. Это обстоятельство позволяет использовать прямой подход к построению траекторий спуска в точки минимума функционала энергии. Такой подход требует предварительного изучения бифуркаций стационарных точек функционала энергии в условиях многомодового вырождения (в порождающей точке минимума). Основы локального анализа в такой ситуации были заложены в работе М.А. Красносельского, Н.А. Бобылева, Э.М. Мухамадиева5и в работах Ю.И. Сапронова, Б.М. Даринского, С.Л. Царева (локальные и нелокальные
3Гармонический анализ линейных операторов / А. Г. Баскаков. - Воронеж: Изд. ВГУ, 1987. -165 с.
4Метод подобных операторов в спектральном анализе несамосопряженного оператора Дирака с негладким потенциалом / А. Г. Баскаков, А. В. Дербушев, А. О. Щербаков // Известия РАН. Сер. матем. 2011. - Т. 75, №3. - С. 3-28.
5Красносельский М.А., Бобылев Н.А., Мухамадиев Э.М. Об одной схеме исследования вырожденных экстремалей функционалов классического вариационного исчисления / ДАН СССР. - 1978. - Т. 240, № 3. - С. 530-533.
бифуркационные задачи)6,7,8,9,10и др.
В диссертации рассмотрены начально краевые задачи для уравнения «реакция-диффузия» с кубической нелинейностью, уравнения Кана-Хилларда, нелинейного обобщения уравнения Фусса-Винклера-Циммер-мана и для уравнения Свифта-Хойенберга — при обычных и обобщенных краевых условиях Дирихле и Неймана. Модельное уравнение «реакция-диффузия» с кубической нелинейностью используется, например, при изучении формирования раскраса шерсти животныхп, а более сложные уравнения Кана-Хилларда и Свифта-Хойенберга — при изучении посткритических фазовых переходов12,13,14,15.
Цель работы.Развитие и применение новых методов бифуркационного анализа актуальных нелинейных начально-краевых задач, соответствующих новым запросам практики и современным достижениям вычислительных технологий. В частности, развитие методов анализа мно-гомодовых и нелокальных бифуркаций.
Методы исследования. Вдиссертации использованы методы функционального анализа, теории нелинейных фредгольмовых операторов, вариационного исчисления, теории особенностей гладких функций и фредгольмовых функционалов, теории приближенных вычислений.
Научная новизна. 1.В диссертационной работе изложена новая (авторская) версия нелокальной редуцирующей схемы Ляпунова-Шмидта (применительно к рассмотренным бесконечномерным динамическим системам).
6Конечномерные редукции в гладких экстремальных задачах / Ю.И. Сапронов // Успехи матем. наук. Т. 51, №1, 101-132 (1996).
7Бифуркации экстремалей фредгольмовых функционалов / Б.М. Даринский, Ю.И. Сапронов, С.Л. Царев // Современная математика. Фундаментальные направления. Том 12 (2004) - С. 3-140.8Гнездилов А.В. Бифуркации критических торов для функционалов с 3-круговой симметрией / Функц. анализ, 2000. Т. 34, вып. 1. - С. 83-86.
9Костин Д.В. Об одной схеме анализа двухмодовых прогибов слабо неоднородной упругой балки / Доклады Академии наук. 2008, Т. 418, № 4, - С. 295-299
10Костин Д.В. Функциональный анализ и многомодовые прогибы упругих систем / Д.В. Костин, Ю.И. Сапронов. Воронеж: ИПЦ ВГУ, 2012, 207 с.
пМарри Дж. Нелинейные дифференциальные уравнения в биологии. Лекции о моделях / М.: Мир. 1983. 399 с.
12Cahn J.W., Hilliard J.E. Free energy of a nonuniform system. I. Interfacial free energy / J. Chem. Phys. 1958. Vol. 28. - P. 258-267.
13Скрипов В.П., Скрипов А.В. Спинодальный распад (Фазовый переход с участием неустойчивых состояний) / УФИ. Т.123, вып.2. 1979. - С.93-231.
14Swift J., Hohenberg P.S. Hydrodynamic fluctuations at the convective instability / Phys. Rev. 1977. V. A15. - P.319-328.
15Кулагин H.E., Лерман Л.М., Шмакова Т.Г. Фронты, бегущие волны и их устойчивость в обобщенном уравнении Свифта-Хойенберга / ЖВМ, 2008, том 48, № 4, с. 693-712

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

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

Впервые построены траектории прямого спуска в точки минимума функционала энергиию из случайно заданных начальных точек (для рассмотренных начально-краевых задач).

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

Полученные общие результаты:
исследованы бифуркации стационарных состояний и траектории спуска бесконечномерных динамических систем типа уравнение «реакция-диффузия» с кубической нелинейностью, уравнение Кана-Хилларда, обобщенного уравнение Фусса-Винклера-Циммермана и уравнениу Свифта-Хойенберга (при обычных и обобщенных краевых условиях Дирихле и Неймана);
предложена новая методика приближенного вычисления ветвей бифурцирующих решений (рассмотренных уравнений) при малых и конечных значениях закритического приращения параметра, созданная на основе вариационной версии процедуры Ляпунова-Шмидта и на использовании ритцевских аппроксимаций ключевой функции по заранее заданному набору собственных функций (мод бифуркаций) главной линейной части градиента функционала энергии;
приведены оценки размера области функционального пространства состояний, на которой допускается нелокальная конечномерная редукция;
в случае локальной редукции найдены главные части ключевых функций и вычислены асимптотические представления ветвей экстремалей по малому закритическому приращению (векторного) параметра;
дано описание алгоритмов и программ соответствующих вычислений (вMaple);
представлены графические изображения линий уровня ключевой функции и функций концентрации вещества, полученных в результате вычисления.
Полученные конкретные результаты. 1.Обоснование примени-
мости методов «фредгольмова анализа» в бифуркационном анализе рассмотренных бесконечномерных динамических систем.

Описание отдельных типовых многомодовых бифуркаций стационарных состояний в случаях рассмотренных уравнений — «реакция-диффузия» с кубической нелинейностью, Кана-Хилларда, обобщенного уравнения Фусса-Винклера-Циммермана и уравнения Свифта-Хойенберга (при обычных и обобщенных краевых условиях Дирихле и Неймана).

Построение и анализ трасс спуска уравнения «реакция-диффузия», редуцированного в подпространство функций с нулевым средним.

Теоремы о главных частях локальных ключевых функций.

Асимптотические представления ветвей бифурцирующих решений.

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

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

Построение компьютерных графических иллюстраций.

Практическая и теоретическая значимость.Работа носит теоретический характер. Представленные в ней научные результаты могут быть использованы в анализе зарождений и развитии посткритических состояний сложных систем.
Апробация результатов диссертации.Результаты диссертации докладывались на ВЗМШ-14, ВЗМШ-15, ВЗМШ-16, ВЗМШ-17, ВВМШ-13 , а также на семинаре по математическому моделерованию (руководитель - проф. В.А. Костин), семинаре проф. Б.М. Даринского по фазовым переходам в кристаллах и семинаре по нелинейному стохастическому анализу (руководитель - проф. Ю.Е. Гликлих).
Публикации.Основные результаты диссертации опубликованы в работах [1 -12]. Работы [2],[6],[10 - 12] опубликованы в журналах из перечня рецензируемых научных журналов и изданий, рекомендованных ВАК Минобрнауки РФ. Из совместных работ [1 - 5] в диссертацию вошли результаты, полученные диссертантом лично.
Структура и объем диссертации.Диссертация состоит из введения, 4 глав и 25 параграфов. Объем работы — 91 страницу. Библиогра-