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

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

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

В метеорологии, астрофизике, ряде других областей ха­

рактерные числа Рейнольдса, при которых турбулентные реше-

»ия представляют наибольший интерес, превышают 10 * 10 .

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

ваемой

асимптотической

неединственностью

[145].

Остановим­

ся на

ней подробнее.

 

 

 

Рассмотрим задачу

Коши для волнового уравнения

 

ии = C2(X)LU,

 

 

 

и(0) = <pQ(x)

exp(tA$0(x)), А » 1,

-со

< х < со.

Начальные данные здесь представляют собой быстро осцил­ лирующую функцию (множитель exp(/AS0(x)) с медленно меняю­ щейся амплитудой осцилляций (множитель V>0(*)). Тогда, применяя асимптотические методы, можно убедиться, что решение такой задачи имеет следующий вид:

и = <p(x,t) exp(/AS(x,/)) + 0(1/А),

где функция S определяется уравнением Гамильтона - Якоби

S2 = c2(VS)2.

То есть решение принадлежит тому же классу быстро осцил­ лирующих функций, что и начальные данные. Естественно ожидать, что это свойство будет характерно и для уравнения Навье - Стокса

д у +

(и,Ч)и + ЧР = ^

Аи,

(V,u)

= 0

(9.65)

 

и = (u^x,t), u2(x,t), ы3(х,/)),

х е Q, Q с IR3, Re » 1

для некоторых классов начальных

данных.

413

В качестве

одного

из

ких

классов

можно

рассматри­

вать течение со

сложным

профилем

скорости

- так

называемую

«сильно когерентную структуру». Начальные данные при этом

имеют вид

 

«|,_0

=

 

+

Г

 

exp{imv(x)/h} aQ(m,x),

 

(9.65)

 

 

 

 

 

m€£

 

 

 

 

 

 

 

 

где А = (Re)_1/3,

x

Й

c

IR3,

Vn,a

C00, 1

-

множество

це-

лых

чисел, Г

ш

0(ш,х)|

<

const

для всех

N £ 0 ,

 

 

 

т€Д

 

 

 

 

 

 

 

 

 

 

 

 

 

а0(-т,х) = Од(/п,х),

 

0^(0,х) = 0,

 

OQ — (орц

2

^03^'

 

 

Функция Уу*) определяет гладкое течение,

члены

в

скобках описывают быстро осциллирующую компоненту.

 

 

Было показано, что в этом случае асимптотическое ре­

шение задачи (9.63) имеет вид

 

 

 

 

 

 

 

и =

V0 + hV^+ ...

+

Г

exp{imv(x)/h}

(a°(m,x,t) +

 

 

 

 

 

 

m€Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ha^(m,x,t) + ...).

(9.67)

 

Возникающие

 

при

построении

такого

ряда

трудности

характеризуются в работе [145] следующим образом: «Мы по­

лучили

цепочку

зацепляющихся соотношений,

в которой N -й

член

влияет

на

первый,

поскольку

уравнение для а0 зависит

от

У 1,

в

свою очередь

V1 зависит

от а1, которое зависит от

V

и

т.

д.

Получается,

что

изменение

члена

при

А

влияет

на

член

А0.

При

подстановке

же

в систему

Навье

-

Стокса

мы

получим

невязку zh:

 

j =

0(h).

Значит, гладкой по­

рядка А невязке отвечают величины, различие между которыми порядка единицы. А это и есть асимптотическая неединствен­

ность. Чтобы аккуратно в этом убедиться, возьмем разложе­ ние в ряд Тейлора по степеням t нашего асимптотического решения. Мы можем из наших соотношений последовательно найти все члены этого ряда и решить задачу, подставив по­

лученный ряд

в систему

Навье -

Стокса.

Тогда

невязка

бу­

дет иметь вид zh + 0(tM).

Выбрав

М

= [jV/e],

где е сколь

угодно мало, мы

получим,

что для

t ~

Ае

невязка имеет

вид

414

гл, т. е. ее N-\ производная по

х, t будет иметь порядок

0(A). В то же

время решения будут

отличаться на о(/)».

Исходя

из этих результатов,

можно предположить, что

при определенных начальных данных и больших числах

Рейнольдса решения уравнения Навье - Стокса

также

неедин-

Ственны

Аналогичная

ситуация,

по-видимому,

возникает в

системах

типа

реакция -

диффузия

в случае, когда в

качест­

ве большого

параметра

выступает

длина области,

и могут

быть использованы асимптотические

разложения.

 

 

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

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

довать вопросы разрешимости этой цепочки уравнений, выяс­

нить число ее решений. Открытым остается и вопрос о том, существует ли в окрестности асимптотического решения,

построенного как некоторый формальный ряд, решение исход­

ной задачи. Новый подход, связанный с представлениями об

асимптотической неединственности, требует дальнейшего исследования.

Переходная турбулентность. Рассматривая хаотические

режимы в нелинейных средах, описываемых детерминированными

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

связано с наличием странного аттрактора в некоторой дина­

мической системе. При этом неявно предполагается, что вре­

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

закончился,

и поведение* системы определяется ее

аттракто­

ром. *

 

 

Поэтому

особый интерес представляет работа

[243], в

которой это

предположение ставится под сомнение

для боль­

шого класса нелинейных сред. Чтобы представить ситуацию, в которой основной интерес представляют не установившиеся, а переходные режимы, рассмотрим, следуя этой работе, отобра­ жение такого вида: .

•Cl =Z7TT.£ К х п‘ ' ) '

(9-68)

1 = - г

 

415

где

 

 

 

 

 

 

f(x)

— sx + u(modl),

x‘n € [0,1 ),

i =

0.....

JV -1 .

Здесь n характеризует временную координату, i

-

простран­

ственную.

Отображение xn+1

= f(xn) устроено таким

образом,

что в течение определенного числа итераций х

возрастает, и

достигая порогового значения х = s-1(l

-

а>),

резко умень­

шается (рис. 9.19).

 

 

 

 

 

Такие отображения ранее использовались в модели ней­

ронов

и

в теории колебательных химических реакций.

Пусть

5 < 1.

В

этом

случае система

не имеет хаотических

режимов

и ведет себя регулярным упорядоченным образом.

В част­

ности,

при

s =

0,91 и и> = 0,1

в ней существует устойчивый

цикл S25

Сумма в правой части формулы (9.68) просто определяет

усреднение по ближайшим соседям. Пусть г = 1. Чтобы оце­

нить, насколько

быстро

происходит

выход на аттрактор,

О

случайных

начальных

конфигураций, рассчи­

задавались 10

танных с помощью метода Монте-Карло. Результат усреднения времен выхода на пространственно однородный или более

416

рожный цикл и давал оценку времени выхода на аттрактор

Рис. 9.20. Типичная картина переходного пространственного хаоса. Заштри­

хованным квадратам соответствуют значения x(i,n ) — 0,5

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

417

пространственной неоднородностью, такое поведение было названо переходным пространственным хаосом. При указанных значениях параметров зависимость TN от длины области N с

хорошей

точностью определяется

формулой (см. рис. 9.21)

 

 

Тн =

Ту'

2

,

9

.

6

9

)

где

 

 

 

 

 

 

 

 

 

 

N = 21,5

± 0,5,

Г,

=

149,5

± 0 , 5

и

<х = 3,0 ± 0,2.

Значение Т^ хорошо согласуется

со

средним

временем

выхода

на

аттрактор

отображения

хп+^

=

/(•*„).

равным

Т ^ -

149,8.

Для других

значений параметров

также

наблю­

дается сверхэкспоненциальный рост. При этом показатель а меняется в интервале 1 < а < 3,2 [243].

Рис. 9.21. Зависимость среднего времени переходного процесса Т от числа элементов N [243]

418

Оказалось, что в этом и ряде других случаев традици­ онный анализ количественных характеристик хаотических

режимов не позволяет отличать переходный пространственный Хаос от хаотических режимов, существующих при t —* ю. На временах t < ТN будет происходить с высокой точностью установление пространственно-временного спектра мощности, функции взаимной информации, распределения вероятностей в разных точках фазового пространст л.

Если переходный хаос имеет отношение к гидродинами­

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

аттракторы таких систем, даже если числа Рейнольдса и соответственно длина области невелики. Авторы работы [243]

оценили по формуле (9.69) время

выхода

на

аттрактор в

модели (9.68), состоящей из 128

ячеек,

при

параметрах,

указанных выше. Если считать, что изучаемый физический

процесс таков, чтО каждая итерация занимает 10-15 сек., то

чтобы дождаться выхода на аттрактор, потребуется Ю40 лет.

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

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

ческих режимов, которые до

настоящего

времени практически

не исследовались.

 

 

Однако сейчас неясно,

насколько

велик класс моделей,

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

выхода на аттрактор растет сверхэкспоненциально. В самом

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

Пространственно - временной хаос в неограниченных об­

ластях. До сих пор мы рассматривали в качестве прообразов хаотических режимов в системах реакция - диффузия диффе­

*

14

419

ренциальные уравнения или отображения, заданные в ограни­ ченных областях. Мы считали также, что в сосредоточенной

системе

аттрактор определяется

циклом.

Однако

недавно в

работе

Л. А. Бунимовича

и Я-Г. Синая

был

предложен

другой

прообраз хаотических

режимов -

так называемый

пространст­

венно -

временной

хаос

[229].

В ней

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

пове­

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

нием, итерации которого ведут себя хаотическим образом.

Этот подход представляет особый интерес для систем реакция - диффузия по двум причинам. В настоящее время

активно изучается большой класс хаотических реакций при

условиях, когда можно не рассматривать пространственные

распределения концентраций. И важно было бы выяснить, что

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

зионные эффекты. Во многих системах основной интерес пред­

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

пространстве в фиксированный момент времени. Такое поведе­

ние трудно исследовать, если двигаться

от

простейших

режи­

мов к

более

сложным, например

увеличивая

длину

области.

Поэтому

особенно

важным оказывается

предельный

случай,

когда

длина области

бесконечна.

 

 

 

 

 

 

 

В работе [229] была рассмотрена

 

следующая

модель.

Пусть

каждой

точке

одномерной

цепочки

i

Z1 (Z1

-

мно­

жество целых чисел) сопоставляется переменная х. € [0,1]. Фазовое пространство изучаемой динамической системы М сос­

тоит из бесконечных в оба конца последовательностей {д^.}, -«о < i < ю. М имеет естественную топологию прямого произ­

ведения

отрезков

и соответствующую <г

- алгебру

1L

Рассмот­

рим функцию /: [0,1] —» R1 такую, что

 

 

 

 

1)

/(0)

= 0;

/(1) =

d, где d г 2 -

целое

число;

 

2)

f

е

С1+аг([0,1])

для некоторого

у >

0

и

/ ' г А =

= const

>

1.

 

 

 

 

 

 

 

420

 

Будем обозначать через f

также

одномерное растягиваю­

щее

отображение /:

х

—» {/(*)},

где

{ }

обозначает дробную

часть

действительного

числа.

 

 

 

 

 

В этом случае f имеет единственную

инвариантную

меру,

которая абсолютно

непрерывна

относительно лебеговой

меры,

ееплотность непрерывна и строго положительна. Обозначим

через

F

отображение

М

на

себя, такое,

что

(Fx).

=

f(x.),

- С О <

{ <

СО.

 

 

 

 

 

 

 

 

 

 

 

Пусть функция а, заданная на отрезке I = [0,1], опре­

деляется

формулами

 

 

 

 

 

 

 

 

 

 

1 )

а(у) = е

для

б ^ г / 5 1 - 6 ,

 

 

 

 

 

 

2 ) а(0) = а(1 ) = 0;

 

 

 

 

 

 

 

 

 

3)

а '(у) -

0 для

0

i

I/

i 6 ; а'(у)

£ 0

для

1 -

у

^ 1 ;

 

4) а е С2([0,1]).

 

 

 

 

 

 

 

 

 

 

Обозначим

через

 

 

величину

max

|a'(i/)|.

Рассмот-

рим оператор усреднения

А

 

» 6 [0 , 1 ]

 

 

 

 

 

 

 

 

 

 

 

 

{Ах). = (1

-

 

 

“ (*■)

 

+• хм ).

 

 

 

 

а(х.))х. + - г ; -

 

 

 

Изучаются только такие операторы А, что AM = М. (В частности, этому условию не удовлетворяет функция а(х) = = е.)

Поведение цепочки осцилляторов определяется отображе­ нием Ф = А о F.

(♦*),

=

 

-

аЩ хЩ х.)

+

Ш (х ))

+

/Ц .+1)).

(1

2 1 (/(*■_,)

 

 

 

 

 

 

 

 

(9.70)

Была доказана следующая

 

 

 

 

Т

е

о р

е м а

9.5. Для

достаточно малых е

и

су­

ществует мера ц, определенная на измеримом пространстве

(М,Ч)

такая,

что

 

 

 

1) д инвариантна относительно отображения Ф;

 

 

2 ) для

любых целых N^ и

индуцированная мера

на

пространстве

конечных последовательностей {*.}, -Л^ <

j <

s JV2,

абсолютно непрерывна относительно меры Лебега;

 

421

3)динамическая система (М, % р, Ф) является пере­ мешивающей;

4)динамическая система (М, % р, S) также является перемешивающей, где S - оператор сдвига в М.

Из этой теоремы следует, что отображение (9.70) обла­

дает сильными эргодическими свойствами: в нем существует и

временной, и пространственный хаос.

В нескольких работах численно исследовались цепочки осцилляторов, динамика которых определяется одномерными

отображениями [64]. Однако, чтобы использовать представле­ ния о пространственно-временном хаосе при решении многих

конкретных задач, нужен дальнейший анализ этого класса ма­ тематических моделей. Важно было бы оценить, насколько ве­

лика

должна быть

длина области, чтобы наблюдалось поведе­

ние,

близкое

к

пространственно -

временному

хаосу.

Интересно было бы выяснить, какие

качественные

эффекты

будут

наблюдаться

при изменении функции f и вида связи а.

По-видимому,

эти

вопросы будут

активно изучаться в

ближайшие годы.

 

 

 

 

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