Исследование аттракторов для некоторых уравнений неньютоновой гидродинамики тема автореферата и диссертации по математике, 01.01.02 ВАК РФ

На правах рукописи

005011726

КОНДРАТЬЕВ СТАНИСЛАВ КОНСТАНТИНОВИЧ

ИССЛЕДОВАНИЕ АТТРАКТОРОВ ДЛЯ НЕКОТОРЫХ УРАВНЕНИЙ НЕНЬЮТОНОВОЙ ГИДРОДИНАМИКИ

01.01.02 — дифференциальные уравнения, динамические системы и оптимальное управление

АВТОРЕФЕРАТ

диссертации на. соискание ученой степени кандидата физико-математических наук

2 4 НОЯ 2011

Воронеж — 2011

005011726

Работа выполнена в Воронежском государственном университете

Научный руководитель: доктор физико-математических наук,

профессор Звягин Виктор Григорьевич. Официальные оппоненты: доктор физико-математических наук,

профессор Орлов Владимир Петрович; кандидат физико-математических наук, доцент Шамин Роман Вячеславович. Ведущая организация: Российский университет дружбы пародов.

Защита состоится 13 декабря 2011 г. в 15.10 на заседании диссерта-ционого совета Д 212.038.22 в Воронежском государственном университете, 394006, г. Воронеж, Университетская пл., 1, ауд. 333.

С диссертацией можно ознакомиться в библиотеке Воронежского государственного университета.

Автореферат разослан «Ш» ноября 2011 г.

Ученый секретарь диссертационного совета Д 212.038.22, доктор физико-математических наук, профессор

ь^^уу^ ГликлихЮ

Общая характеристика работы

Актуальность темы. Изучение траекторпых н глобальных аттракторов пространств траекторий играет важную роль в задачах математической физики и, в частности, гидродинамики. Теория аттракторов динамических систем была перенесена на уравнения математической физики во второй половине прошлого века в работах нескольких математиков; в частности, О. А. Ладыженская, получила классический результат о существовании глобального аттрактора двумерной системы Навье—Стокса. Однако эта теория оказалась неприменима уже к трёхмерной системе Навье—Стокса, для которой не удаётся построить динамической системы в силу отсутствия результатов о глобальном существовании слабых решений начально-краевой задачи или единственности слабых. В обход этих трудностей в конце прошлого века в работах М. И. Вишика, В. В. Чепыжова, а также Дж. Селла была развита теория траскторных и глобальных аттракторов пространств траекторий. В этих теориях пространство траекторий является развитием понятия динамической системы, а глобальный аттрактор является аналогом аттрактора полугруппы. С помощью этой теории удалось доказать существование траекторного и глобального аттракторов трёхмерной системы Навье—Стокса. Приложение теории аттракторов пространств траекторий к другим задачам гидродинамики потребовало обобщения этой теории на случай трансляцнонно неинварнантных пространств траекторий, поскольку такие пространства естественно возникают во многих приложениях. Такое обобщение было сделано в работах В. Г. Звягина и Д. А. Воротникова. Именно этот вариант теории аттракторов используется в настоящей работе.

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

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

Методика исследований. Использовались методы теории аттракторов пространств траекторий, а также топологические методы нелинейного анализа и методы и идеи теории дифференциальных уравнений.

Научная новизна. Все основные результаты диссертации являются новыми. Среди них можно отмстить наиболее важные:

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

2. Доказано существование минимального траекторного и глобального аттракторов слабых решений регуляризованной системы уравнений движения жидких сред с памятью при условии определённого ограничения на параметры системы.

3. Построена визуализация аттракторов возмущений течения Пуа-зёйля в модели Джеффриса.

Достоверность полученных результатов подтверждается математическими методами исследований. Все результаты диссертации доказаны.

Теоретическая и практическая значимость. Работа имеет теоретический характер. Полученные результаты могут применяться при исследовании различных течений вязкоупругих жидкостей и сред.

Результаты диссертации могут найти применение в исследованиях проводимых в научных коллективах Российского университета дружбы народов, Московского, Санкт-петербургского, Воронежского и Ярославского государственных университетов.

Апробация работы. Результаты диссертации докладывались и обсуждались на воронежской зимней математическая школе «Современные методы теории функций и смежные проблемы» (Воронеж, 2011); международной научной конференции, носвящёшюй 70-лстию члена-корреспондента АН Республики Таджикистан Мухамадиева Эр-гашбоя Мирзоевпча (Душанбе, Таджикистан, 2011); семинаре иод руководством профессора В. Г. Звягина (ВГУ, 2011); семинаре иод руководством профессора А. Л. Скубачевского (РУДН, 2011).

Публикации. Основные результаты диссертации опубликованы в работах автора [1]—[6]. Из совместных работ [1), [2] в диссертацию вошли только принадлежащие С.К.Кондратьеву результаты. Работы [2|, [4] соответствуют перечню ВАК.

Структура и объём диссертации. Диссертация состоит из введения, трёх глав, разбитых на 20 параграфов, и списка литературы, включающего 35 источников. Общий объём диссертации 139 страниц.

Краткое содержание работы

Нумерация приводимых ниже теорем совпадает с их нумерацией в диссертации.

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

ствующим аттракторам (минимальному траскторпому и глобальному) исходной системы.

Математической моделью движения слабокоицеитрированпых вод-пых растворов полимеров называется система

1=1 \г=1 /

сНуи = 0, (2)

где ь(х, £) — вектор скорости частицы жидкости, находящейся в точке х в момент времени £; р(х,<) —давление жидкости в точке х в момент времени /(ж, £) — вектор плотности внешних сил; £ = — тензор скоростей деформации,

с с < \ 1 (дг>' , дуз

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

Уравнения (1), (2) рассматриваются в ограниченной области Л С С Мп (п = 2,3) с гладкой границей. На границе области задастся условие прилипания:

V

Также система снабжается начальным условием

V = а, (4) (=0

где а — функция нз некоторого функционального пространства.

Рассматривается автономный случай: / 6 (£2(0))" не зависит от времени Ь.

Пусть {Ка}аек — шкала пространств, порождённая оператором Сток-са А = — 7ГА, где 7г — проектор Лерэ.

В качестве пространства траекторий Н+ задачи (1)—(3) рассматривается множество функций класса

= 0. (3)

да

таких, что тождество

^ ^ (1х + Vu(í) : Чуйх + и J : 4<рйх-

п а п

г,],к-1 п п

выполнено почта всюду на М+ для любой функции ¡р € V3, и имеет место неравенство

»«(ОН! + ИОН-, < До (1 + ||«||1зс(к+,к.,е-а')

для почти всех < > 0 (здесь Д» на- некоторые постоянные, не зависящие от у).

В диссертации показано, что для каждого а е V1 существует траектория у 6 Н+, такая, что г/(0) = а.

Основные результаты о существовании аттракторов формулируются следующим образом.

Пусть 5 6 (0,1] — некоторое число.

Теорема 1.4.1 Пусть / € (Ь2(0,))п. Тогда существует минимальный траекторный аттрактор Ы пространства траекторий /Н+. Аттрактор ограничен в ¿оо(®+; V1), компактен в С(М+; он притягивает в топологии пространства С(Ж+; V1"5) семейства траекторий, ограниченные в £оо(1К+; V1).

Теорема 1.4.2 Пусть / е (2/2(0))"- Существует глобальный аттрактор Л пространства траекторий 'Н+. Аттрактор ограничен в Vх, компактен в V1-"5; он притягивает в топологии пространства V1'6 семейства траекторий, ограниченные в Ьос(М+; V1).

Задача (1)—(3) изучается с помощью аннрокснмацпопного уравне-

пия

п

А

а

йх-дхк

(1 Р

! Vv{t): Чрёх + ее-«1^ ^ У(Д«(0) : У(Ду>)<*г+

я п

+ Vу Уг;^) : У^йж-^.(^(О^г^-

г,],к=1 п п

= 11<рЛх (реУ3), (6) п

где е > 0 — параметр.

В качестве пространства траекторий уравнения (6) рассматривается множество, состоящее из решений этого уравнения, определённых на К+, существенно ограниченных со значениями в V1 и удовлетворяющих оценке

Ф(0)|Ц < 1, (7)

а также из функций вида Тгде и — решение уравнения

(6) на К+,

/г > 0, н Т(/г)и удовлетворяет оценке

Отметим, что неравенство (7) оправдывается тем, что именно такие решения аппроксимационного уравнения используются при предельном переходе к решению исходного уравнения.

Для задачи (1)—(3) рассмотрим новое пространство траекторий добавив к Н+ функции, удовлетворяющие тождеству (5) почти всюду на К+ для любой функции е V3, а также оценке

Показывается, что пространства траектории Н* (г > 0) имеют минимальные траекториыс и глобальные аттракторы. Кроме того, имеет место следующая теорема о сходимости аттракторов:

Теорема 1.4.3. Минимальные траекторные аттракторы Ы; ап-проксимациоиного уравнения (0) сходятся к мишшалыюму траек-торному аттрактору Ы{) пространства траекторий задачи (1) — (3) в смысле полуотклоиения в пространстве С(К+; V1"6), то есть

Глобальные аттракторы Л£ аппрокегшациоииого уравнения (6) сходятся к глобальному аттрактору Ло пространства траекторий Н^ задачи (1)—(3) в смысле полу отклонения в пространстве I/1-'5, то есть

Во второй главе рассматриваются аттракторы слабых решений регуляризованной системы уравнений движения жидких срсд с памятью.

Пусть П — ограниченная область пространства Ж", га = 2,3, с гладкой границей. Рассматривается следующая начально-краевая задача для регуляризованной модели движения жидких срсд с памятью:

вир М ||и - ч-О (е-* 0).

эир Ы ||г/ - ->-0 (г ->■ 0). уеАс г6Л0

- Ц1 БЬ' / е '* £(ю)(э, Zfi(s■,t,x)) дв-

- Ку £(ь) = - grad р + /, (<, х) 6 (0, +оо) х П (8)

сИуг> = 0, {Ь,х)€П

V = 0,

(9) (10)

(И) (12)

г

ь'(0,х) = х € П

п

где v = (v\,... ,vn) — вектор-функция скорости движения частиц среды, р — давление среды, / — плотность внешних сил, не зависящая от t, £ = (¿"у)— тензор скоростей деформации, Л > 0, цо > 0, pi ^ О — постоянные.

Поясним обозначение Zs(s;t,x). Для этого введём стандартные в математических вопросах гидродинамики пространства V н Я.

Обозначим через ("D(íl))n пространство функций класса С°° из П в К", имеющих компактный носитель в fI. Положим

V={ve{V(ü))n: div и = 0} .

Пусть Я —замыкание V по норме пространства (£2(0))", a F —замыкание V по норме пространства (W2 (fi))"; Я — гильбертово пространство со скалярным произведением, индуцированным (L2(^))n, а на V рассматривается норма ||v|]v = ||£(w)||¿2.

Рассмотрим траекторию, определяемую уравнением

т

z(t; t,x) = х + J Ssv{s,z{s;t,x))ds, re [О, Г], (t,x) 6 (0,T) x íl

(13)

В уравнении (13) используется оператор регуляризации Я —> —>• С1(П)п П V для S > 0. Он обладает следующими свойствами: Ss(v) —> v в Я при 5 —J- 0, и порождаемое им отображение Ss: 12(0,Т; Я) L2(0,T; С^П)" П V) непрерывно.

Для каждого v S ¿2(0j Т; V) уравнение (13) имеет единственное решение Zs(v). Из(13) имеемz(t;t,х) = Zs(T\t,x). Если« £ 1-20С(К+; V), то можно считать, что функция Z¿(r;í,x) определена при т > 0 и t > 0. Эта функция и участвует в (8).

Далее предполагаем выполненным условие

Но - > 0. (14)

В качестве пространства траекторий задачи (8)—(12) рассматривается множество функции V б ¿00(К+; Я) П ¿2°П(®+;У) с производной г/ е ¿11ос(К+; К*), таких, что ограничение V па любой отрезок [О, Г] является слабым решением задачи (8)—(12) с некоторым и0 £ Я, зависящим от у, и при любом £ ^ 0 выполняются неравенства

Iе-2*<-">||«/(5)||?, < С? (1 + (16)

«

с некоторыми постоянными С\ > 0 и 7 > 0, не зависящими от V.

Показывается, что для каждого гг° € Я существует траектория V € И+, такая, что г>(0) = г>°.

Пусть в е (0,1]; Уц — замыкание V в пространстве (11-^(0))", Уд* — сопряжённое к этому пространству.

Теорема 2.2.1 Пусть параметры ро, /¿ь Л задачи (8)—(12) связаны соотношением (14) и / 6 V*. Тогда существует минимальный траекториый аттрактор Ы пространства траекторий Н+ задачи (8)—(12). Аттрактор ограничен в ^(Ж-,.; Я), компактен в С(К+; Уд); он притягивает в топологии пространства С(К+; Ув*) ограниченные в норме Ь00(Ш+;Н) множества траекторий.

Теорема 2.2.2 Пусть

параметры /¿о, ^ задачи (8)—(12) связаны соотношением (14), и / € У*. Тогда существует глобальный аттрактор Л пространства траекторий 'Н+ задачи (8)—(12). Аттрактор ограничен в Н, компактен в У^; он притягивает в топологии пространства Уд ограниченные в норме ¿ос(Ж+;Я) множества траекторий.

В третьей главе строится визуализация аттракторов возмущений течения Пуазсйля для среды Джсффриса.

Предположим, что область

О = {{хих2)-. -оо<хх< +оо,-1 < х2 < 1} С Ж2.

заполнена несжимаемой средой единичной плотности, подчиняющейся реологическому соотношению Джеффриса. Пусть V = (г;ь г;2) — ее поле скорости, р — давление, а — девиатор тензора напряжений; v, р и а зависят от координат на плоскости Х\, х2 и времени Предполагается, что иоле внешних сил потенциально. Далее, пусть £ = £(У) — тензор скоростей деформации. Тогда мы имеем следующие уравнения, выполняющиеся при (х1,х2) € £ > 0:

^ О д1 ^ дх{

(18) 1191

Здесь (17) —уравнение Коши движения среды, (18)—условие несжимаемости, (19) — реологическое соотношение Джеффриса с полной производной, г] — вязкость среды, Ах — время релаксации, А2 — время ретардации; г/, Аь А2 — положительные постоянные, и А1 > Аг-

На горизонтальных границах рассматриваем условие прилипания:

V =0. (20)

12=±1

С помощью замены = т]\г!\\ > 0, ц2 = (г/ - > 0, т =

= а - 2//1(£(у) система (17)—(19) записывается в виде

ЪГ + У] - А V + ётайр = г,

<)1 ^ ОХг

г=1

сну v = о,

2

дт 1 ^ 9т

1=1 12

Рассматриваются возмущения течения Пуазёйля в виде бегущей волны:

V = уР + \,р = рР+р,т = Тр + Т,

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

V = 11е (е^'-^У^)) , р = Ке(е^>-^Р(х2)), (21)

г = Ле (е"1^-^^));

здесь

|Ы12!)' Т(х2) = (Т^Ы).

\«Ы )

Возмущения указанного вида делятся на затухающие, периодические и возрастающие. Пространство траекторий строится из затухающих п периодических возмущении; его траскторный аттрактор состоит из пулевого н периодических возмущений, которые в фазовом пространстве изображаются соответственно точкой и эллипсами. Ненулевые возмущения существует не для любых значений параметров а > О н с € С. Для значений

/ц = ИГ4, ц2 = 1(Г5, А! = 0.1.

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

Публикации автора по теме диссертации

[1] Кондратьев С. К. Аттракторы для уравнений моделей движения вязкоупругих сред / В. Г. Звягин, С. К. Кондратьев; Воронежский государственный университет. — Воронеж: Издательско-полиграфический центр Воронежского государственного университета, 2010.-266 с.

[2] Кондратьев С. К. Аттракторы слабых решений регуляризоваи-ной системы уравнений движения жидких сред с памятью / В. Г. Звягин, С. К. Кондратьев. // Известия вузов. Математика. — 2011.-№ 8.-С. 86-89.

[3] Кондратьев С. К. Визуализация аттракторов малых возмущений течения Пуазёйля для системы Джеффриса с полной производной /С.К.Кондратьев // Современные методы теории функций и смежные проблемы : материалы ВЗМШ-2011.—Воронеж, 2011. — 374 с.-С. 176-177.

[4] Кондратьев С. К. Об аттракторах модели движения слабоконцентрированных водных растворов полимеров / С.К.Кондратьев // Вестник В ГУ. Сер. Физика. Математика. — 2010. — № 1, — С. 117— 138.

[5] Кондратьев С. К. О сходимости аттракторов аппроксимаций модели движения слабоконцентрированных водных растворов полимеров / С. К. Кондратьев // Современные проблемы математики и её приложений : материалы международной научной конференции, посвященной 70-летию члена-корреспондента АН Ресиубли-

ки Таджикистан Мухамаднева Эргашбоя Мирзоевича. — Душанбе, 2011.-С. 58-59.

[6] Кондратьев С. К. О сходимости траекторных и глобальных аттракторов аппроксимации автономной трехмерной системы На-вьс—Стокса / С. К. Кондратьев // Вестник ВГУ. Сер. Физика. Математика. - 2009. - № 1. - С. 120—137.

Работы [2), [4] соответствуют перечню ВАК.

Подписано s печать 3.10.11. Формат 60*84 '/i6. Усл. печ. л. 0.9. Тираж 100 экз. Заказ 1365.

Отпечатано с готового оригинал-макета в типографии Издательско-иолиграфического центра Воронежскою государственною университета. 394000, Воронеж, ул. Пушкинская, 3

Введение

1 Аттракторы модели движения слабоконцентрированных водных растворов полимеров

1.1 Уравнения движения

1.2 Элементы теории аттракторов траекторных пространств

1.3 Функциональные пространства и обозначения.

1.4 Постановка задачи и основные результаты.

1.5 Свойства операторов.

1.6 Априорные оценки.

1.7 Существование решений

1.8 Пространство траекторий и аттракторы.

1.9 Сходимость аттракторов

1.10 Сходимость аттракторов аппроксимаций.

2 Аттракторы слабых решений регуляризованной системы уравнений движения жидких сред с памятью

2.1 Система уравнений движения жидких сред с памятью

2.2 Основные результаты.

2.3 Вспомогательные утверждения: свойства используемых операторов и оценка решений.

2.4 Пространство траекторий задачи (2.1.1)—(2.1.2).

3 Визуализация аттракторов возмущений течения Пуазёйля для среды Джеффриса

3.1 Общий подход к визуализации аттракторов.

3.2 Возмущения течения Пуазёйля в модели Джеффриса

3.3 Аттракторы возмущений

3.4 Численные результаты и визуализация

3.5 Вывод уравнения возмущений.

3.6 Вспомогательная спектральная задача

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

Понятие «аттрактор» возникло в теории динамических систем, которая изначально была тесно связана с теорией обыкновенных дифференциальных уравнений. Во второй половине XX в. для уравнений с частными производными в работах О. А. Ладыженской, А. В. Бабина, М. И. Вишика, Р. Темама и других математиков была построена теория бесконечномерных динамических систем. Одним из первых её приложений было доказательство существование глобального аттрактора двумерной системы Навье— Стокса в работах О.А.Ладыженской (см. [18]).

Однако применение теории аттракторов динамических систем в задачах математической физики и, в частности, в гидродинамике ограничено тем, что для построения динамической системы требуется однозначная разрешимость задачи с заданным начальным условием. Однако уже для трёхмерной системы Навье—Стокса не установлено ни существование глобального сильного решения, ни единственность слабого, так что построить динамическую систему не удаётся.

В обход указанных трудностей возникла теория траекторных аттракторов, позволяющая строить глобальные аттракторы для ряда эволюционных уравнений. Эта теория была предложена российскими учёными М. И. Ви-шиком и В. В. Чепыжовым (см., например, статью [29] и монографию [28]) и независимо американским учёным Дж. Шеллом (см. [32]), а впоследствии была усовершенствована с целью применения в задачах неньютоновой гидродинамики в работах В.Г.Звягина и Д.А.Воротникова (см. монографию [36], а также статьи [4,34,35]. Именно этот последний подход к аттракторам используется в настоящей работе.

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

Для пространств траекторий рассматриваются глобальные и траектор-ные аттракторы. Глобальный аттрактор является обобщением аналога аттрактора динамической системы; он представляет собой некоторое множество в фазовом пространстве, притягивающее ограниченные семейства траекторий. Траекторный аттрактор — это определённое множество функций времени со значениями в фазовом пространстве, на которые с течением времени «становятся похожи» траектории. Между понятиями траекторного и глобального аттрактора существует тесная связь: например, если существует минимальный (по включению) траекторный аттрактор, то его сечение в произвольный момент времени (то есть множество значений в некоторый момент времени всех функций, составляющих этот аттрактор), является глобальным аттрактором.

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

Одна из первых задач, в которой нашла приложение теория траектор-ных аттракторов,— трёхмерная система Навье—Стокса (см. [3], [28], [29]). Подход Вишика и Чепыжова накладывал на пространство траекторий ограничительное требование транляционной инвариантности, поэтому его приложение к многим задачам неньютоновой гидродинамики было затруднено, поскольку требовало получения инвариантных относительно сдвигов по времени оценок решений. После того, как в работах В. Г. Звягина и Д. А. Воротникова условие трансляционной инвариантности было снято, удалось построить аттракторы для различных моделей движения вязкоупругих сред, в частности, для системы Джеффриса с полной производной (см. [34]). Изложение теории аттракторов пространств траекторий, от которых не требуется трансляционной инвариантности, можно найти в [34,36].

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

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

Приведём обзор содержания диссертации по главам.

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

Движение несжимаемой жидкости единичной плотности описывается системой в форме Коши где — вектор скорости частицы жидкости, находящейся в точке х в момент времени р(х, — давление жидкости в точке х в момент времени /(х, £) — вектор плотности внешних сил (в рассматриваемом автономном случае / = /(х)); о = (<з^(х)) —девиатор тензора напряжений (симметрическая матрица порядка п)\ символ Бит означает вектор, компоненты которого равны дивергенции соответствующих строк матрицы.

В системе (1), (2) неизвестными являются V, р, а, и общее число неизвестных превосходит число уравнений. Поэтому эта система обычно дополняется некоторыми соотношениями между девиатором тензора напряжений а и тензором скорости пргьппи/гяиыы ? — (р. л

Эти соотношения называются реологическими; различные реологические соотношения описывают различные классы жидкостей. Самое известное реологическое соотношение имеет вид п

1) сЦу у = О,

2) о = 2у£,

3) где V > 0, при котором система (1), (2) превращается в уравнения Навье— Стокса.

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

В работе [23] для учёта релаксационных свойств было предложено реологическое соотношение

4) где V — кинематический коэффициент вязкости, к — время ретардации (положительные величины); й д д г=1 полная производная по времени. После подстановки о в уравнение (1) приходим к следующей системе: ду . ^ ду дАу /^ д£(у)\ , . ,г. а - + Е -К1Г -2к 1I £ ■+ ^ = ■(5) г=1 \г=1 /

Ну V = 0. (6)

Отличие данной системы от системы Навье—Стокса состоит в присутствии члена 2к.ТУ\ч{6£/(1£). Если течение жидкости близко к установившемуся, то этот дополнительный член оказывает малое влияние, и жидкость ведёт себя практически как ньютонова. В случае же турбулетного или ламинарного пульсирующего течения (то есть при наличии местных ускорений) значение дополнительного члена возрастает, и движение значительно отличается от движения ньютоноской жидкости. Экспериментально было установлено (см. [ 1 ]), что именно такое поведение характерно для слабокон 0- (7) ап V центрированных водных растворов полимеров (полиэтиленоксида, полиа-криламида), поэтому система (5), (6) стала называться математической моделью движения слабоконцентрированных водных растворов полимеров.

Уравнения (5), (6) будем рассматривать в ограниченной области с К" (п = 2,3). На границе области рассмотрим условие прилипания: V

Также система снабжается начальным условием V где а — функция из некоторого функционального пространства.

Разрешимость в сильном и слабом смысле различных начально-краевых задач для системы (5), (6) рассматривалась в работах А. П. Осколкова [21], [22] (результаты которого о слабых решениях содержали ошибки, что отмечалось самим автором), а также В. Г. Звягина и М. В. Турбина [9]. Метод доказательства существования слабых решений, предложенный в последней работе, не позволяет получить оценок диссипативного типа, необходимых для доказательства существования аттракторов. В связи с этим оказалось необходимым найти новое доказательство существования, дающее нужную оценку.

Для доказательства существования слабых решений используется ап-проксимационно-топологический метод. Для используемых аппроксимаций также устанавливается существование минимальных траекторных и глобальных аттракторов, а также устанавливается их сходимость к соответствующим аттракторам исходной задачи. Подобные задачи (для га-лёркинских приближений рассматривались в работах М. И. Вишика и В. В. Чепыжова [28,29] и в работе автора [16] в контексте теории аттракторов инвариантных пространств траекторий. Однако установление сходимости аттракторов неинвариантных пространств представляет дополнительные трудности, поскольку минимальный траекторный аттрактор может не содержаться в пространстве траекторий, а в таком случае мы не обладаем непосредственной информацией о составляющих его функциях.

Во второй главе исследуются траекторный и глобальный аттракторы автономной регуляризованной системы уравнения движения жидких сред с памятью: + Е'- ^ Div J* e-{t~s)/xS(v)(S, Z5(s; t, x)) dsjioDiv5(v) = — gradp + /, divv = 0. (9)

Показано, что если параметры системы удовлетворяют неравенству цо — —цД > 0, то система имеет минимальный траекторный и глобальный аттракторы.

В системе (9) участвует Z5(s; t, х) — траектория, определённая по ре-гуляризованному полю скоростей Ss(v) (Sg— оператор регуляризации). Необходимость его введения оператора регуляризации связана с тем, что само слабое решение v € L2(0,T; V) не обладает достаточной гладкостью, чтобы определить траекторию. Введение той или иной регуляризации является стандартной операцией при исследовании задач гидродинамики, в которые входит определение траекторий точки через вектор-функцию скорости (см., например, [24,30]).

Существование слабых решений уравнений движения (9) в более общем неавтономном случае устанавливается в работах В.Г.Звягина, В.Т.Дмит-риенко [10,11]. Она подробно описана в монографии [36].

Во второй главе настоящей работы для слабых решений, построенных в работах [10,11], устанавливаются оценки диссипативного типа. На их основе определяется пространство траекторий, для которого доказывается существование минимального траекторного и глобального аттракторов.

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

Существование траекторного и глобального аттракторов системы Джеффриса (а также равномерных аттракторов в неавтономном случае) рассматривалось в работах В. Г. Звягина и Д. А. Воротникова [4,34,35]. С целью получения визуализации в главе 3 настоящей работы рассматривается более узкое пространство траекторий, состоящее из течений, близких к течению Пуазёйля.

Опишем рассматриваемую модель подробнее.

В плоском канале {(х\,х2) е Е2: — 1 ^ х2 ^ 1} рассматривается движение вязкоупругой среды, удовлетворяющей реологическому соотношению Джеффриса с«)

Здесь г| обозначает вязкость среды, Х\ — её время релаксации, Х2 — её время ретардации; ц, %1 и Х2 — положительные постоянные, и ^ > ?12- Полная система, описывающая движение, состоит из уравнения (10) и уравнений движения несжимаемой среды (1), (2), причём в уравнении (1) полагаем / = 0, что соответствует потенциальному полю внешних сил. На границе канала рассматривается условие несжимаемости.

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

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

В случае системы Навье—Стокса аналогичный подход приводит к уравнению Орра—Зоммерфельда.)

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

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

Изучению течений, близких к течению Пуазёйля для ньютонових жидкостей, посвящён ряд работ, среди которых отметим [26,27,31,33]. В частности, в этих работах используются численные методы. Однако вопросы визуализации аттракторов в этих работах не рассматриваются, как не рассматриваются и неньютоновы жидкости.

Суммируя вышеизложенное, отметим, что в диссертации получены следующие новые результаты:

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

2. Доказано существование минимального траекторного и глобального аттракторов слабых решений регуляризованной системы уравнений движения жидких среде памятью при условии определённого ограничения на параметры системы.

3. Предложена визуализация аттракторов возмущений течения Пуазёйля в модели Джеффриса.

Результаты диссертации докладывались и обсуждались на воронежской зимней математическая школе «Современные методы теории функций и смежные проблемы» (Воронеж, 2011); международной научной конференции, посвященной 70-летию члена-корреспондента АН Республики Таджикистан Мухамадиева Эргашбоя Мирзоевича (Душанбе, Таджикистан, 2011); семинаре под руководством профессора В. Г. Звягина; семинаре под руководством профессора А. Л. Скубачевского (РУДН, 2011).

Исследования, включённые в настоящую диссертацию, поддержаны грантами РФФИ-04-01-00081-а, РФФИ-07-01-00137-а, РФФИ-10-01-00143-а и федеральной целевой программой «Научные и научно-педагогические кадры инновационной России на 2009—2013 годы» (руководитель грантов и программы — профессор В. Г. Звягин).

Основные результаты диссертации опубликованы в работах автора [8, 13— 15], из них [8] и [ 14] соответствуют перечню ВАК для кандидатских диссертаций. Из совместной работы [8] в диссертацию вошли только принадлежащие автору результаты.

1. Амфилохиев В. Б. Течения полимерных растворов при наличии конвективных ускорений / В. Б. Амфилохиев, Я. И. Войткунский, Н.П.Мазаева, Я-С. Ходорковский // Труды Ленинградского ордена Ленина кораблестроительного института. — 1975. — Вып. 96.— С. 3-9.

2. Бетчов Р. Вопросы гидродинамической устойчивости / Р. Бетчов, В. Криминале. — М.: Мир, 1971. — 350 с.

3. Вишик М. И. Траекторный и глобальный аттракторы ЗЭ системы На-вье—Стокса / М. И. Вишик, В.В.Чепыжов // Матем. заметки.— 2002. — Т. 71, вып. 2. — С. 194—213.

4. Воротников Д. А. О траекторных и глобальных аттракторах для уравнений движения вязкоупругой среды / Д. А. Воротников, В. Г. Звягин // УМН. — 2006. — №2. — С. 161 — 162.

5. Дмитриенко В. Т. Конструкции оператора регуляризации в моделях движения вязкоупругих сред// В. Т.Дмитриенко, В. Г. Звягин // Вестник ВГУ. —2004. — № 2. — С. 148—153.

6. Звягин В. Г. Аппроксимационно-топологичёский подход к исследованию задач гидродинамики. Система Навье—Стокса / В. Г. Звягин, В. Т. Дмитриенко. — М.: Едиториал УРСС, 2004. — 112 с.

7. Звягин В. Г. Аттракторы слабых решений регуляризованной системы уравнений движения жидких сред с памятью / В.Г.Звягин, С. К. Кондратьев. // Известия вузов. Математика. — 2011. — №8. — С. 86— 89.

8. Звягин В. Г. Исследование начально-краевых задач для математических моделей движения жидкостей Кельвина—Фойгта / В. Г. Звягин, М.В.Турбин // Гидродинамика. СМФН. — 2009. — Т. 31. —С. 3— 144.

9. Звягин В. Г. О слабых решениях начально-краевой задачи для уравнения движения вязкоупругой жидкости / В. Г. Звягин, В. Т. Дмитриенко //Доклады РАН. —2001. Т. 380, № 3. —С. 308— 311.

10. Звягин В. Г. О слабых решениях регуляризованной модели вязкоупругой жидкости / В. Г. Звягин, В. Т. Дмитриенко // Дифференциальные уравнения. —2002. —Т. 38, № 12, —С. 1633—1645.

11. Картан А. Дифференциальное исчисление. Дифференциальные формы / А. Картан. — М.: Мир, 1971. — 392 с.

12. Кондратьев С. К. Визуализация аттракторов малых возмущений течения Пуазёйля для системы Джеффриса с полной производной /С.К.Кондратьев // Современные методы теории функций и смежные проблемы : материалы ВЗМШ-2011. — Воронеж, 2011.— 374 е. —С. 176—177.

13. Кондратьев С. К. Об аттракторах модели движения слабоконцентрированных водных растворов полимеров / С. К. Кондратьев // Вестник В ГУ. Сер. Физика. Математика.— 2010. — № 1. — С. 117—138.136

14. Кондратьев С. К. О сходимости траекторных и глобальных аттракторов аппроксимаций автономной трехмерной системы Навье—Сток-са / С.К.Кондратьев // Вестник ВГУ. Сер. Физика. Математика.— 2009. —№ 1. —С. 126—137.

15. Крейн С. Г. Шкалы банаховых пространств / С. Г. Крейн, Ю. И. Петунин // УМН. — 1966. —Т. 21, вып. 2(128). —С. 89— 168.

16. Ладыженская O.A. О нахождении минимальных глобальных аттракторов для уравнений Навье—Стокса и других уравнений с частными производными / О.А.Ладыженская // УМН. — 1987. — Т. 42, вып. 6. — С. 25—60.

17. Лионе Ж.-Л. Неоднородные граничные задачи и их приложения /Ж.-Л. Лионе, Э. Мадженес. — М.: «Мир», 1971. —371 с.

18. Наймарк М. А. Линейные дифференциальные операторы / М. А. Наймарк. — М.: «Наука», 1969. — 528 с.

19. Осколков А.П. Начально-краевые задачи с краевым условием проскальзывания для модифицированных уравнений Навье — Стокса / А. П. Осколков // Зап. науч. семин. ПОМИ. — 1994. — Т. 213. — С. 93—115.

20. Осколков А. П. О некоторых нестационарных линейных и квазилинейных системах, встречающихся при изучении движения вязких жидкостей/А. П. Осколков//Зап. науч. семин. ЛОМИ. — 1976. — Т. 59. — С. 133—177.

21. Павловский В. А. К вопросу о теоретическом описании слабых водных растворов полимеров / В.А.Павловский // ДАН СССР. — 1971.— Т. 200, №4. — С. 809—812.

22. Темам Р. Уравнения Навье—Стокса. Теория и численный анализ / Р. Темам. — М.: Мир, 1981, — 288 с.

23. Фурсиков А. В. Оптимальное управление распределёнными системами. Теория и приложения / А. В. Фурсиков. — Новосибирск : Научная книга, 1999, —350 с.

24. Afendikov A. Bifurcations of Poiseuille Flow between Parallel Plates: Three-Dimensional Solutions with Large Spanwise Wavelength / A. Afendikov, A.Mielke // Arch. Rational Mech. Anal.— 1995.— Vol. 129 —P. 101—127.

25. Afendikov A. Multi-pulse solutions to the Navier—Stokes problem between parallel plates / A.Afendikov, A.Mielke // Z. Angew. Math. Phys — 2001 — Vol. 52 — P. 79—100.

26. Chepyzhov V. V. Attractors for equations of mathematical physics / V. V. Chepyzhov, M. I .Vishik.— Providence, RI: AMS Colloquium Publications, 2002.— 363 p.

27. Chepyzhov V. V. Evolution equations and their trajectory attractors / V.V.Chepyzhov, M.I.Vishik // J. Math. Pures Appl— 1997 — Vol. 76 —P. 913—964.

28. Litviniv V. G. Regular model and nonstationary problem for the nonlinear viscoelastic fluid / V. G. Litvinov// Siberial Journal of Diff. Eq.— 1997.— V. 1. No. 4 —P. 351-382.

29. Orszag S. Accurate solution of the Orr—Sommerfeld stability equation / S. Orszag// J. Fluid Mech.— 1971— Vol. 50, Iss. 4 — P. 659—703.

30. Sell G.R. Dynamics of Evolutionary Equations / G. R. Sell, Y. You.— New York : Springer, 1998.—670 p.

31. Skorokhodov S:L. Numerical Analysis of the Spectrum of the Orr— Sommerfeld Problem / S. L. Skorokhodov // Computational Mathematics and Mathematical Physics.— 2007.— Vol. 47, Iss. 10.— P. 1603— 1621.

32. Vorotnikov D. A. Trajectory and global attractors of the boundary value problem for autonomous motion equations of viscoelastic medium / D. A. Vorotnikov, V. G. Zvyagin // J. Math. Fluid Mech — 2008 — Vol. 10.—P. 19—44.

33. Vorotnikov D. A. Uniform attractors for non-automous motion equations of viscoelastic medium / D. A. Vorotnikov, V. G. Zvyagin // J. Math. Anal. Appl.— 2007 — Vol. 325 — P. 438—458.

34. Zvyagin V. G. Topological Approximation Methods for Evolutionary Problems of Nonlinear Hydrodynamics / V. G. Zvyagin, D. A. Vorotnikov.— Berlin, New York : Walter de Gruyter, 2008 — 232 p.