![](/user_photo/_userpic.png)
книги / Нестационарные структуры и диффузионный хаос
..pdfСтохастические режимы были обнаружены и в средах, которые описываются уравнением
qt + iqxx - 2i\q\2q = (е ,- €2\q\2)q +
+ c3qxx - |
en |
00 |
7^ |
E exp(2rcmf/T). |
|
|
|
n = — CO |
Здесь также удается построить динамические системы небольшой размерности, позволяющие описывать хаотические режимы [347].
§ 9.6. Априорные оценки размерности аттракторов
Обычно при исследовании процессов в нелинейных сре
дах, в которых наблюдаются стохастические режимы, размер
ность странного аттрактора оценивается в ходе вычислитель
ного или натурного эксперимента. Однако в последние годы появились математические теории, позволяющие получить
априорные оценки размерности аттракторов уравнений в частных производных.
Идеально было бы по данному уравнению в частных производных (или бесконечномерному отображению) строить
некоторую конечномерную систему, |
аттрактор |
которой |
был |
бы |
в определенном смысле близок к |
аттрактору |
исходной |
задачи. |
|
По-видимому, таких результатов |
в настоящее время |
нет. |
Однако на этом пути были получены принципиальные результа ты, развиты новые представления, которые позволяют глубже
понять особенности хаотических режимов в |
нелинейных |
зада |
|
чах. Обратим внимание на некоторые из |
них. |
|
|
1. Первоначально большинство |
результатов в |
этой |
области касались гидродинамической турбулентности и урав
нения |
Навье |
- |
Стокса. Однако развитые подходы оказались в |
||||
достаточной |
мере универсальны |
и |
применимы к |
большому |
|||
классу |
систем |
реакция - |
диффузия |
[28], уравнению Кура- |
|||
мото - |
Цузуки |
[232,253] |
и ряду |
других уравнений |
в частных |
производных.
391
Будем рассматривать движение вязкой несжимаемой жидкости в двумерной ограниченной области Q, которое описывается уравнением Навье - Стокса
gy + v*Vv = -Vp + uV2v + f,
|
|
V-v = О, |
(9.21) |
|
|
|
|
где v |
= |
характеризует скорость жидкости, |
р - дав |
ление, |
v - |
вязкость, f - объемная сила, которая |
считается |
не зависящей от времени. Обычно рассматривают двумерную область Я, на границе 8S которой v = 0, либо прямоугольник
Я с периодическими граничными условиями. Нас будут интерес совать установившиеся режимы течения жидкости.
|
Из-за наличия вязкости система является диссипатив |
|
ной, |
поэтому естественно ожидать, что характеристики тече |
|
ния |
будут определяться ее аттрактором. |
По-видимому, |
впервые внимание на тесную связь между теорией турбулент ности и исследованием аттракторов динамических систем было обращено в работе О.М.Ладыженской [130]. Ряд замечаний, касающихся истории этого направления, и обширная библио графия содержатся в обзорах [29, 131].
Рассматривая другие нелинейные среды, мы обращали внимание на существование конечного числа параметров
порядка, эффективно описывающих динамику изучаемого про
цесса. |
Этот вопрос подробно исследовался для уравнения |
Навье - |
Стокса. |
При традиционном гидродинамическом подходе число
степеней свободы, эффективно описывающих турбулентное те
чение, определяется |
как |
|
|
|
|
|
||
|
N ~ |
|
( / / / / , |
|
|
|
|
(9.22) |
где d = 2 или d - |
|
3, в зависимости |
от |
размерности |
прост |
|||
ранства, |
в котором |
|
рассматривается задача, lQ - |
характер |
||||
ный линейный размер |
области, |
занятой |
жидкостью; |
- |
малый |
|||
масштаб, |
начиная |
с |
которого |
эффекты, |
связанные |
с |
вяз- |
392
костью, полностью определяют |
движение. |
Было |
показано, что |
|||
в |
трехмерном |
случае |
|
|
|
|
|
/с = |
v3/4/cV4, |
где |
е = у< Vu |
>2, |
(9.23) |
в |
двумерном |
|
|
|
|
|
|
1с = |
(ц3/% )1/6, |
X = |
v< Ды >2, |
|
(9.24) |
где прямые скобки < > соответствуют усреднению по ансамблю [133, 238]. В двумерном случае в настоящее время удалось получить строгие результаты, близкие к оценке (9.22).
При получении таких оценок уравнение Навье — Стокса записывают в следующем функциональном виде:
|
|
|
|
|
+ vAu + В(и,и) = Pf, |
|
|
|
(9.25) |
||||
где |
v - |
безразмерный |
коэффициент |
кинематической |
вязкости. |
||||||||
При |
этом |
рассматривают |
функциональные |
пространства векто |
|||||||||
ров |
Н и V, Н - |
замыкание множества |
|
|
|
|
|
||||||
|
|
|
V = |
{ |
w е |
IR2 , |
(о е |
С " | Я |
, d i v |
о> = |
0 |
} |
|
в пространстве |
о |
функций, |
интегрируемых |
с |
квадратом, со |
||||||||
L |
|||||||||||||
скалярным |
произведением |
|
|
|
|
|
|
|
|||||
|
|
|
2 |
|
|
|
dx' |
и = (« г ы2); |
v |
= |
(or |
v2). |
|
|
(«. |
») |
= Е._1 |
£ Y / |
V - замыкание множества V в пространстве Соболева со скалярным произведением
|
|
|
|
2 |
ди. dv. |
|
|
<u’ v> = |
Е. , I |
ЯхГЯГ. dx’ |
|
Р в |
формуле |
(9.25) |
- |
оператор |
ортогойального проектирова |
ния |
из Z.2 на |
Н\ А = |
-ЯД, а В(и,и) = Р(иЧ)и. |
393
Наряду с уравнением Навье - Стокса, обычно рассматри вают задачу на собственные значения, называемую задачей
Стокса:
V2w + 4g = -Aw
|
|
|
|
|
|
div w = 0, |
|
|
(9.26) |
|
|
|
|
|
|
|
|
|
|
||
где |
на |
функцию |
w |
накладываются те |
же |
граничные |
условия, |
|||
что |
и |
на |
функцию |
v в уравнении Навье - |
Стокса. Так |
же как |
||||
градиент |
давления |
определяется |
уравнением (9.21), |
градиент |
||||||
величины |
g определяется задачей (9.26). |
|
|
|
||||||
|
Домножим |
уравнение |
(9.21) на |
функцию w с div w = 0 |
||||||
(это |
может быть |
одна |
из |
собственных |
функций |
задачи |
||||
(9.26)). |
Тогда |
|
|
|
|
|
|
|
||
|
|
|
(Vp*w) |
= Г div(pw) |
dx = J р w п ds, |
|
||||
|
|
|
|
|
Q |
|
ds |
|
|
|
что в силу краевых условий равно нулю. Поэтому скалярное умножение на такую функцию (проектирование с помощью оператора П) позволяет исключить давление (см. (9.25)).
Задача Стокса также может быть представлена в виде
Лю = Лю.
Она имеет набор собственных значений
0 < Л^ - Л2 - Л3 s ...—» со
и собственных векторов
wr w2............
Будем разлагать решение уравнения Навье - Стокса по этим собственным векторам wfr):
v(r,/) = Е eft) v/fr).
394
В работе [266] приведен следующий результат. Пусть п достаточно велико - настолько, что выполнено любое из неравенств
А. > 2у-2 lim supllvll2,
|
|
|
|
|
/->ю |
|
|
|
|
|
|
Лп+1 |
> у-2 |
lim |
sup Af2, |
(9.26) |
|||
|
|
|
|
/->00 |
|
|
|
|
|
|
|
An+1 |
> у-1 |
lim |
sup М2, |
|
|||
|
|
|
|
/->00 |
|
|
|
|
|
где |
и Af2 - |
максимальные |
амплитуды v и rot v во всей |
||||||
области |
Q в данный момент |
времени. Тогда |
|||||||
|
|
|
V ~ v(n) |
= Е сIt) |
w /r) |
||||
|
|
|
|
|
|
/= 1 |
‘ |
|
* |
в том |
смысле, |
что |
если |
при |
t |
—> |
оо |
v является постоянной, |
периодической или квазипериодической функцией, то такой же функцией является и у^п\ Можно ожидать, что такое соотно шение будет иметь место и для хаотических’ режимов. Однако
этот результат, вероятно, пока не получен.
Утверждения такого типа оказываются очень важными.
Они являются теоретической основой для экспериментального
исследования вязкой |
жидкости. Из |
них следует, что |
измере |
||
ние конечного набора |
величин, |
зависящего |
от |
числа |
|
Рейнольдса, позволяет получить информацию о типе |
течения. |
||||
2. Следующий |
вопрос |
связан с |
оценкой сложности |
наблю |
даемых хаотических режимов или с оценкой размерности ат тракторов уравнений в частных производных. При его иссле
довании вводится важное понятие максимального |
аттрактора |
||||||
[29]. |
Полугруппой {S<t t £ 0}, отвечающей уравнению (9.21), |
||||||
|
|||||||
называется семейство |
операторов S(, |
действующих |
в банахо |
||||
вом |
пространстве |
Е, |
S(: Е —> |
Е, |
причем S(uQ = |
u(t), |
где |
u(t) |
- решение |
задачи (9.21). |
При |
этом предполагается, |
что |
задача (9.21) однозначно разрешима. Максимальным аттракто ром полугруппы {S(} или уравнения (9.21) называется такое замкнутое ограниченное множество 11 с Е, которое обладает
395
1 ) свойством |
инвариантности: |
S(U = И для / ^ |
0; |
|
|||
2) |
свойством |
притяжения: |
для |
любого |
ограниченного |
||
множества В с |
Е |
расстояние |
dist£ |
(S((B), |
11) —» 0 |
при |
|
t —* «, |
где |
|
|
|
|
|
|
|
dist |
(X, Y) = sup |
inf IIx |
- y\\ . |
|
|
|
|
|
|
*€* |
yGY |
E |
|
|
Максимальный аттрактор существует для большого класса диф ференциальных уравнений и отображений. Достаточные условия его существования определяются следующей теоремой [29].
|
Т |
е о р е м |
а |
9.1. |
Пусть |
полугруппа {S(} |
удовлетво |
|
ряет следующим условиям: |
|
|
|
|
||||
|
1 ) полугруппа {S^ равномерно ограничена, то есть для |
|||||||
любого |
R > 0 существует C(R) такое, что IIS^II < |
C(R), если |
||||||
Hull |
£ |
R 0 £ / £ со; |
|
|
|
|
|
|
|
2) существует компактное в Е поглощающее множество |
|||||||
В^, |
т. |
е. для любого |
ограниченного множества В с |
Е |
сущест |
|||
вует такое число Т, |
что при t г Т |
SfB с |
|
|
||||
|
3) операторы |
S(E —* |
Е непрерывны при t 2: |
0. |
Тогда у |
полугруппы {S,} имеется компактный максимальный аттрактор.
Максимальный аттрактор определяет асимптотическое поведение решений изучаемых уравнений при всех возможных начальных данных. Поэтому он может описывать множество предельных режимов более широкое, чем множество режимов, наблюдаемых в натурном или вычислительном эксперименте.
Это можно проиллюстрировать на примере так называемых градиентных систем
ЭТ = - «(“ ) Ш’ “ <“ ) > °- |
(9-28) |
где дФ/ди - вариационная производная. Уравнение (9.28) обладает глобальной функцией Ляпунова, значения которой убывают вдоль траекторий уравнения всюду, кроме стационар ных точек. При выполнении естественных условий максималь ные аттракторы таких уравнений имеют вид
U = U М(г.),
1 =1 1
396
где |
{г^} |
- |
множество |
стационарных |
решений |
уравнения |
||||||
(9.28), s |
- |
число |
таких |
решений, М(г!) |
- |
максимальное |
||||||
неустойчивое |
инвариантное |
многообразие, |
|
проходящее |
через |
|||||||
точку |
г.. |
Многообразие |
М(г^ |
состоит |
из |
всех |
траекторий |
|||||
ц(т), выходящих при возрастании т из точки г . |
|
|
|
|||||||||
|
Таким |
|
образом, |
наряду |
со |
стационарными |
точками |
макси |
мальный аттрактор содержит и их неустойчивые многообразия.
Это делает его даже в простейших случаях достаточно слож
ным объектом. Например, в уравнении (3.8) при любых значе ниях параметров существует бесконечное количество симмет ричных решений, которые могут определять неустойчивые многообразия. Они также должны входить в максимальный ат
трактор, |
как |
и |
те |
множества, |
которые характеризуют |
поведе |
||||||||
ние |
решений |
при |
t —* со в случае типичных |
начальных данных. |
||||||||||
|
Условия |
теоремы 9.1 |
могут быть |
проверены |
для |
уравне |
||||||||
ния |
Навье |
- |
Стокса |
[27,29]. |
Например, |
для |
функции u(x,t) |
|||||||
выполнено |
неравенство |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
ехр(-Л vt) + |
« г » , |
|
|
|
||
|
llu(r,0H |
о - |
IIып(г)Н „ |
— — , |
|
|
(9.29) |
|||||||
|
|
|
L |
|
и |
LT |
|
1 |
v\x |
|
|
|
||
где |
Aj— |
первое |
собственное значение задачи (9.26). Следо |
|||||||||||
вательно, |
при |
t —» со траектории |
будут |
лежать |
внутри шара |
|||||||||
|
|
|
Sf : |
|
Hull |
9 s |
||/|| |
/A.W. |
|
|
|
|
(9.30) |
|
|
|
|
|
' |
|
|
1Г |
Г |
1 |
|
|
|
|
|
|
Рассматривая |
конечномерные динамические |
системы, мы |
обращали внимание на связь между динамическими характе ристиками системы, которые определяются ее ляпуновскими показателями, и геометрией аттрактора, его фрактальной
размерностью. Эта связь определялась формулой |
Каплана- |
|||
Йорке |
[259]. Оказывается, эта |
связь существует и для урав |
||
нений |
в |
частных производных. |
Именно она и дает ключ к |
|
оценке хаусдорфовой размерности максимальных аттракторов. |
||||
3. В конечномерной системе сумма первых k показателей |
||||
Ляпунова |
определяла скорость |
уменьшения A-мерных |
фазовых |
397
объемов. Обобщая это свойство, |
Ю.С.Ильяшенко |
ввел понятие |
Л-сжимающих систем. |
|
|
О п р е д е л е н и е . |
Отображение области евклидова |
|
или гильбертова пространства в |
себя называется |
k-сжимающим |
с константой q < 1 , если оно уменьшает A-мерные объемы, умножая их на коэффициент, не превосходящий q.
Оценку хаусдорфовой размерности максимального аттрак тора в конечномерном случае дает следующая теорема [105].
Т е о р е м а 9.2. Пусть g - дважды гладкое k-сжимаю- щее отображение замыкания В ограниченной области евклидова
пространства в |
себя. |
Тогда отображение g |
имеет притягива |
||
ющее множество (аттрактор) А = |
gnB, |
хаусдорфова |
раз |
||
мерность которого не превосходит k. |
|
|
|
||
Оказалось, |
что |
утверждение |
теоремы |
естественно |
обоб |
щается и на бесконечномерные системы. При этом могут быть получены результаты, применимые к широкому классу различ
ных |
уравнений |
[27, |
29]. |
|
|
|
|
|
|
|
|
|
|
|
||||
|
Например, в работе [27] приведен следующий результат. |
|||||||||||||||||
|
Пусть |
полупоток |
S( таков, |
что |
|
|
|
|
|
|
||||||||
|
IIS/и ) |
- |
S/»)ll |
2 — Ии-оИ |
2 |
exp(Bt) |
|
|
|
|
(9.31) |
|||||||
(в |
случае |
|
|
уравнения |
|
Навье |
- |
Стокса |
В |
= |
A^/4v, где |
|||||||
Ап= |
шах |
шах |
|ы(х)|),и выполнено |
неравенство |
|
|
|
|
||||||||||
|
иеи |
х бП |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
II (Е - |
Pn)(Sfu - |
S(t>)ll |
2 |
£ |
и |
llu |
— v\\ 2, |
w < |
1 , |
|
(9.32) |
||||||
где |
Р - |
оператор |
ортогонального |
проектирования |
на |
конечно |
||||||||||||
мерное |
подпространство, |
натянутое |
на |
векторы |
|
|
w1......w^. |
|||||||||||
(В случае уравнения |
Навье |
- |
Стокса |
w - |
собственные |
векторы |
||||||||||||
задачи |
(9.26), |
а |
неравенство |
(9.32) |
было |
|
установлено |
|||||||||||
О. А.Ладыженской). |
Тогда |
хаусдорфова |
размерность |
максималь |
||||||||||||||
ного аттрактора |
полугруппы |
S( конечна. |
|
|
|
|
|
|||||||||||
|
Более |
того, |
П.Константином |
и И.Фояшем было показано, |
||||||||||||||
что |
предельную |
емкость |
и |
|
хаусдорфову |
размерность |
макси |
398
мального аттрактора можно оценить сверху с помощью формулы Каплана - Йорке, куда вместо обычных ляпуновских показате лей должны входить так называемые глобальные ляпуновские показатели [237].
Обозначим через Р ^ 0) произведение
* У 0 “ 0) = mx(t, uQ) ... mN(t, «0),
где m.(t, uQ) связан с ляпуновскими показателями А.(«0) соотношением
|
|
\(и |
) = |
lim |
(1 / 0 |
log |
m(t,uQ). |
|
|
|
||||
|
|
1 |
|
|
t-*a> |
|
|
|
1 |
|
|
|
|
|
Пусть аттрактор лежит внутри шара |
радиуса |
р |
В"(0). |
Введем |
||||||||||
величины |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ у о |
= |
sup |
PN(t,u0), |
N i l , |
|
|
||||||
|
|
|
|
|
v V 0) |
|
|
|
|
|
|
|||
|
|
mlt) |
= |
sup |
mU,uQ). |
|
|
|
|
|||||
|
|
|
|
|
V |
V |
° |
> |
|
|
|
|
|
|
Глобальные ляпуновские показатели Д. и fi. определя |
||||||||||||||
ются следующими |
формулами. Для / а 1 и л/ 2: 1 |
|
|
|
||||||||||
|
|
Д. = |
lim |
sup |
т- |
log |
1 |
|
|
|
|
|||
|
|
1 |
|
t-* oo |
|
1 |
|
|
|
|
|
|
||
|
|
Яд, = lim |
} |
log |
/у * ) . |
|
|
|
|
|
||||
|
|
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
|
P] = nv |
|
PN = |
|
|
|
|
|
|
|
|||
В работе |
[237] |
были |
доказаны следующие утверждения. |
|||||||||||
I е |
о р |
5 м |
а |
9.3. |
Пусть |
X - |
инвариантное |
множество, |
||||||
которое |
является |
ограниченным |
в |
V, например, |
X |
- |
макси |
мальный аттрактор. Тогда хаусдорфова размерность определя ется формулой
р, + --+ М /о
£ /о |
(9.33) |
|
399
Подчеркнем, что из большой размерности максимального аттрактора не следует, что будут наблюдаться хаотические режимы. При типичных начальных данных может происходить
выход на простейшее периодическое или стационарное ,реше-
ние. Вместе с теи оценки хаусдорфовой размерности - гаранти
руют, |
что фрактальная размерность странного аттрактора, |
если он |
наблюдается, не превышает данной величины. |
Было бы интересно выяснить, как в натурных или вычис лительных экспериментах меняется размерность странных аттракторов, которые описываются уравнениями Навье - Сток са при изменении обобщенного числа Грасгофа, когда оно достаточно велико. Однако, по-видимому, имеющаяся экспери ментальная техника и возможности компьютеров не позволяют проводить такие исследования.
5.Методы и представления, связанные с оценками
максимального |
аттрактора |
уравнения Навье |
- |
Стокса, |
оказа |
|||||||
лись в полной мере применимы к уравнению Курамото - |
Цузуки |
|||||||||||
[232, |
253]. |
|
|
|
|
|
|
|
|
|
|
|
|
Рассмотрим |
|
эти |
результаты |
более |
подробно, |
|
следуя |
||||
работе [232]. Перепишем уравнение Курамото |
- |
Цузуки |
в |
виде |
||||||||
А( = |
RA + |
(1 |
+ |
i»)Axx - |
(1+/д)4|Л|2, х € |
[0,1], |
|
|
||||
4(*,0) = AQ(X), |
4(0,/) |
= 4(1,/), |
|
|
|
|
(9.41) |
|||||
|
|
|
|
|
||||||||
считая, |
что оно |
задано |
в |
области |
единичной |
длины |
с |
перио |
||||
дическими краевыми условиями. |
|
|
|
|
|
В обсуждаемой работе приводится априорная оценка функции A(t) в норме La (то есть во множестве ограниченных
функций). Аттрактор X уравнения (9.41) лежит |
в шаре |
радиу |
|||
са R в пространстве L2> |
|
|
X удовлетворя |
||
Норма функции 114(х,/)11г |
на |
аттракторе |
|||
|
ли |
|
|
|
|
ет неравенству |
|
|
|
|
|
lim |
supll4112 £ R + 28R2{\ + [1 |
+ |
(1 + 6 )/S 2fl]1/2}, |
(9.42) |
|
t->a |
X |
|
|
|
|
|
00 |
|
|
|
|
где б = max (0 ,-2 + |1 +