Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Нестационарные структуры и диффузионный хаос

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
14.86 Mб
Скачать

Стохастические режимы были обнаружены и в средах, которые описываются уравнением

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)

в

двумерном

 

 

 

 

 

 

=

(ц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

 

 

 

 

 

 

Рассматривая

конечномерные динамические

системы, мы

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

размерностью. Эта связь определялась формулой

Каплана-

Йорке

[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 +

Соседние файлы в папке книги