
книги / Нестационарные структуры и диффузионный хаос
..pdfРейнольдса |
в |
такой задаче |
определяется по |
формуле |
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