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

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

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

Рейнольдса

в

такой задаче

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

формуле

Re =

= Sl\<i(b -

a)u-1

при

Qg = 0 (v

- кинематическая

вязкость).

 

При

небольших

числах Рейнольдса движение жидкости ла-

минарно,

его

можно

описать

явной формулой.

Далее при

уве­

личении частоты Qj возникают диссипативные структуры - так называемые вихри Тейлора. При увеличении частоты вращения

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

(аналог инвариантного тора) и далее при

Re = Rec

движение

становится турбулентным. Таким образом,

переход

к хаосу,

наблюдаемый в эксперименте при определенных значениях Slya,b,L,v, происходит в соответствии со сценарием Рюэля-

Такенса [265].

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

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

рения

этой

компоненты

дают

массив

.... а^.

В

работе

[225]

проводилась статистическая обработка более 32 тыс.

точек.

Рассматривая переменную

а,(0» можно

выделить

сред­

нее

время Т,

через которое эта переменная принимает

одно и

то

же

значение а(^у > 0) (в динамических системах это

среднее время

возвращения

на плоскость Пуанкаре).

В

экспе­

рименте выбирался временной интервал ~ 300 Т (300 орбит), примерно по 100 точек на орбите (Д/~10 Т). По этим дан­

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

гическая энтропия. По ним же оценивалась размерность ат­

трактора изучаемой системы. Оказалось, что при увеличении числа Рейнольдса (Re > Rec) размерность аттрактора увели­

чивается,

однако в интервале Rec s

Re s 1, 3

Rec оиа не

превышает

5,4. Таким образом, было

убедительно

показано,

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

253

небольшой размерности, описывающего течение Куэтта - Тей­

лора, оказывается

справедливым.

Этот важный

эксперимен­

тальный результат

указывает на

принципиальную

возможность

описания таких явлений с помощью сравнительно

небольшого

числа обыкновенных дифференциальных уравнений. Близкие за­

кономерности имеют место и для течения жидкости между вра­ щающимися сферами [205, 206]. Такая система представляет

интерес в связи с задачами физики атмосферы.

Другой пример перехода от диссипативных структур к турбулентным режимам дает конвекция Рэлея - Бенара, проис­

ходящая в подогреваемом снизу слое жидкости [322]. Обра­ ботка соответствующего эксперимента проводилась по 15 ты­

сячам

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

At (Т/Ы ~ 2 +

* 10).

По ним определялся корреляционно

показатель. Вна­

чале задавалась размерность р ^-векторов и рассчитывался

корреляционный показатель V , далее

величина

р

увеличива­

лась до тех пор, пока значение

v

не

 

переставало

меняться.

Отметим, что в работе расчет проводился

не по всем N2 рас­

стояниям, а по расстояниям от

л1 точек до остальных точек

множества (в работе [322] п1 =

100,

в

работе

[275]

=

= 15). Этот прием, позволяющий намного сократить объем вы­ числений, но уменьшающий точность определения величины V,

широко используется в настоящее время

в физических

рабо­

тах.

Экспериментальное

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

рэлей-бенаровской

неустойчивости показало,

что во многил

конвективных

тече­

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

небольшой размерности (2,5 £ р s 6). Увеличение числа Рейнольдса или переход к системам, где нужен учет других

усложняющих факторов, обычно приводит к росту размерности [322].

В упоминавшихся работах рассматривалось движение жид­

кости в

ограниченном

объеме, поэтому

естественно

ожидать,

что из-за вязкости движение на

малых

масштабах

будет по­

давляться

и аттрактор

окажется

конечномерным.

В

связи с

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

254

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

реактивной струи [223] и течение типа «затопленной струи»

[73]. В этих работах были вычислены корреляционные показа­

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

нии пространственной координаты. Оказалось, что и здесь движение определяется странным аттрактором небольшой раз­ мерности.

Г=ЩЮС(&)

 

df

 

 

 

V

I

/

 

 

 

♦♦

 

 

 

 

 

 

 

 

 

 

 

 

-*<**W>

-1,0 - 0 ,5

0,5

-1 ,0 - 0 ,5

0,5

 

е=1о9еде

 

®=г°9«,е

 

Рис.

6.20

 

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

странных аттракторов в естествознании. Успех анализа хао­

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

зависит от

того,

насколько

успешно решен

ряд методических

вопросов.

Обратим внимание

на

некоторые

из них.

 

 

 

 

 

В большинстве упоминавшихся

выше работ

по

временному

ряду измерений вычисляется корреляционный показатель. Ти­

пичная зависимость корреляционного

интеграла и тангенса

наклона кривой f = log1QC от величины

е =

log1Qe,

приве­

денная в работе [226], представлена

на

рис.

6.20.

В этой

255

статье проводилась обработка данных, полученных в натурном эксперименте при анализе течения Куэтта - Тейлора. На кривой, определяющей зависимость s = df/de можно выделить четыре области, онн типичны для гидродинамических и ряда других задач.

Область А В этом интервале масштабов размеры ячеек

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

сколько точек, онн не позволяют оценить вероятности р.

 

Область В. В этой области

выборка недостаточна,

чтобы

передать канторову структуру аттрактора. Кроме того,

в ней

обычно существенна погрешность

эксперимента или точность,

с которой известны точки аттрактора.

 

Область С. Точки кривой в

этом интервале характеризу­

ют фрактальную размерность аттрактора. Чтобы расширить об­ ласть С, можно увеличить размер выборки, повысить точность

эксперимента,

наилучшим образом выбирать переменные, кото­

рые

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

ды

обработки

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

 

Область D. Размеры ячеек разбиения стремятся к разме­

рам

аттрактора и не характеризуют его канторову структуру.

При исследовании странных аттракторов более высокой размерности картина оказывается близкой, однако с. увеличе­ нием размерности уменьшается длина участка С. Для того,

256

чтобы увеличить ее, часто используют прием, предложенный в

работе

[380].

 

 

 

 

Пусть выборка невелика: N < 2т

, где т

-

время

убы -

вания

автокорреляционной функции (в

единицах

At),

т -

раз­

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

ячейке разбиения будут преобладать точки, лежащие на том

же участке траектории (см. рис. 6.21). На этих масштабах

аттрактор будет выглядеть как одномерная кривая. Зависи­

мость log2C = /(log2e) в этом случае окажется искаженной.

Соседние точки на траектории, попадающие в одну ячейку по­

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

отражают свойства аттрактора, а характеризуют только спо­ соб выбора At. Чтобы восстановить зависимость С(е), отра­

жающую свойства аттрактора, нужно просто не учитывать бли­

жайших соседей, лежащих на одном отрезке траектории. Это можно сделать, если модифицировать соотношение (6.63) для

корреляционного

интеграла

следующим образом:

 

 

C(r,N,W)

= £ 3 Е

Ni n e(r -\x1+п- х

(\).

(6.83)

При W = 1

формула

(6.83) переходит в

стандартное вы­

ражение для корреляционного интеграла. Значение W > 1 со­ ответствует отбрасыванию того или иного числа ближайших

соседей

на

траектории.

Естественно

отбрасывать

число

точек

W, пропорциональное т (именно они и лежат на линейном

участке).

Тестовые

примеры,

разобранные

в

статье

[380],

показывают,

что

в

качестве

W

лучше

брать величину

т 1п ( т / 2 ).

Вычисляя С по этой формуле, можно «убрать одномерный участок» и получить зависимость, близкую к той, которая

характеризует канторову структуру изучаемого множества. Во

многих случаях,

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

ких систем, этот

прием оказывается

очень полезным.

 

В работе [140] он используется при исследовании тур­

булентных

режимов течения

вязкой

жидкости в

круглой трубе

при

Re =

4000.

Зависимости

x(t) в

этом случае

были получе-

9 Т.С.

Ахромеева и

др.

 

 

 

257

ны

в результате вычислительного эксперимента. В качестве

x(t)

используются различные компоненты скорости в несколь­

ких

точках.

. 6

2

Рис. 6.22

На рис. 6.22 показана типичная картина течения в бес­

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

трехмерная краевая задача для уравнения Навье - Стокса. Показаны мгновенные отклонения от средних значений скорос­

ти в каждой точке на один момент времени в двух перпенди­

кулярных плоскостях, проходящих через ось трубы. Результа­

ты этого расчета были любезно предоставлены авторам В.Г.Приймаком. Обширная библиография, касающаяся трехмер­ ных гидродинамических расчетов, приведена в работе [Д20].

Интервал дискретизации At (x(nkt) = x j выбирался с

помощью критерия, в котором используется взаимная информа­ ция (6.82). При расчете корреляционного интеграла исполь­

зовалась

формула

(6.74).

Обрабатывались выборки длиной

2 * 104 +

6 - 104 .

 

 

Зависимость

In С =

Д1п е) для различных компонент

скорости и тангенса наклона от In е в одной точке трубы показана на рис. 6.23. Видно, что v = 15 ± 2. В силу того,

что

размерность достаточно велика, длина линейного участка

мала.

Увеличение

размерности

вложения

(до тех

пор, пока

она

не слишком велика) практически не

меняет

наблюдаемой

картины.

 

 

 

 

 

*р=20

о р = 3 0 орш 40

+р-Ю *р-*20 ор =30 ар=>40

Рис. 6.23. Зависимость корреляционного показателя

и

локального

наклона,

рассчитанного по г и ^-компонентам скорости, в точке с

радиальной коор­

динатной

Г = 0,596 (в единицах D /2).

Видно, что

значение V,

определенное

тангенсом

наклон" линейного участка,

практически

не

зависит

от

размер—

 

ности вложения. Длина выборки

N

-

4

 

 

 

10

 

 

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

ристики течения, если

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

точках трубы.

Очевидно,

эти скорости будут меняться в раз­

ных пределах,

поэтому

в

каждом случае надо заново выбирать

значения At.

На рис.

6.24 показано,

как меняется наклон

259

кривой In С = Д1п е) в разных точках трубы. В пределах точности расчетов корреляционный показатель, вычисленный по одной из компонент скорости в различных точках, оказы­ вается одним и тем же.

+р<*10 *p =ZO op =30 o p -4 0 +p=10 *p -ZO o p -3 0 ap=40

Рис. 6.24. Зависимость корреляционного

показателя

и

локального

наклона,

рассчитанная по г-компонеитам скорости

в

точках

с

радиальными

координа­

тами г = 0,290 и г

=

0,989

 

 

 

В настоящее время очень мало известно о геометрии

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

Вместе с тем во многих случаях геометрические характерис­

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

Приведем пример, связанный с обсуждавшейся выше зада­ чей. При изучении геометрии аттракторов в теории динами­

ческих систем полезную информацию дает использование так

называемого базиса Карунена -

Лоева {£ } [140], состоящего

из собственных векторов матрицы

ковариаций

260

 

 

А

г

п+Г

Ц

=

J>k<i* , = I и^к+

 

 

 

 

п

 

 

*./

 

 

 

+

( t -

1)Д < ) U(tk +

(j

-

1 )Д 0 ;

i. i = 1.

P-

 

Здесь {jtn}

-

анализируемая

выборка,

представляющая

набор

^-векторов, построенных по

значениям

величины U,

измеря­

емой в моменты tk, tk + bt,

.... Обычно tk = kbt; A

-

мат­

рица pxp,

xk ( -

i—я

координата

вектора

 

 

 

Рис. 6.25

На рисунке 6.25 показана проекция отрезка тректории аттрактора на плоскость £2). Траектория построена по значениям z-компоненты скорости в точке с радиальной коор­ динатой г = 0,596. Наблюдаемая картина ясно указывает на существование быстрого и медленного временных масштабов. Быстрым движениям соответствует мелкая рябь на отвечающих медленным движениям орбитах большого диаметра. Это под­ тверждают и результаты частотного анализа [140].

261

Рис. 6.26

 

 

 

Анализ проекций аттрактора

на

различные

плоскости

(£t„ £.) позволяет обратить внимание

на

следующую

интерес­

ную закономерность. Будем следить (рис. 6.26) в течение

небольших, следующих друг за другом, интервалов времени за

проекцией той же траектории (£3, £4). На каждом интервале

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

как если бы она лежала в окрестности сложного цикла или

инвариантного

тора. В

течение следующего промежутка

време­

ни она оказывается

в

окрестности другого неустойчивого

цикла (тора).

В итоге

на

больших временах проекция

выгля­

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

262

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