Главная > Регулярная и стохастическая динамика (Лихтенберг А., Либерман М.)
<< Предыдущий параграф Следующий параграф >>
Пред.
След.
Макеты страниц

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

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

Рассмотрим теперь связь исследованных выше конечномерных хаотических потоков с движением сплошной среды («жидкости») с бесконечным числом степеней свободы. Сначала мы обсудим вопрос, в каком смысле конечномерные модели аппроксимируют жидкость, а затем опишем некоторые механизмы перехода от регулярного движения к турбулентному.
7.4а. Представление Фурье

Будем описывать поведение жидкости с помощью дифференциальных уравнений в частных производных вида
\[
\frac{\partial \boldsymbol{Q}(x, t)}{\partial t}=\hat{\mathscr{L}}(\boldsymbol{x}) \boldsymbol{Q}(\boldsymbol{x}, t),
\]

где $\boldsymbol{Q}-M$-мерный вектор состояния жидкости [компонентами $\boldsymbol{Q}$ могут быть, например, давление $p(\boldsymbol{x}, t)$, скорость жидкости $v(\boldsymbol{x}, t)$, плотность $\rho(\boldsymbol{x}, t)$ и т. д. ], $\boldsymbol{x}$ – обычный радиус-вектор с компонентами $x, y, z$, а $\hat{\mathscr{L}}(\boldsymbol{x})$ – независящий от времени нелинейный дифференциальный оператор. Обычный метод изучения системы (7.4.1)_состоит в переходе к представлению Фурье для вектора $Q$ :
\[
\boldsymbol{Q}(\boldsymbol{x}, t)=\sum_{k} q_{k}(t) e^{i k \cdot \boldsymbol{x}},
\]

где
\[
q_{k}(t)=\frac{1}{(2 \pi)^{3}} \int d^{3} x \boldsymbol{Q}(\boldsymbol{x}, t) e^{-i k \cdot x} .
\]

Подставляя (7.4.2) в (7.4.1) и используя ортогональность функций $e^{i k \cdot x}$, получаем уравнения движения в виде ${ }^{1}$ )
\[
\dot{q}_{k}=V_{k}\left(q_{1} . . q_{k}\right) .
\]

Если оставить в сумме (7.4.2) $N$ «наиболее существенных» мод, то задача сведется к конечному числу ( $M N$ ) обыкновенных диффе-
1) Такое представление годится, конечно, только для модельных задач с простейшими граничными условиями. В более общем случае обычно используется разложение по собственным функциям соответствующего системе (7.4.1) линейного уравнения (см., например, книгу [534]). – Прим. ред.

ренциальных уравнений первого порядка, описывающих временну́ю эволюцию этих мод. Такой метод называется приближением Галё́ркина.

Рассмотрим в качестве примера задачу Рэлея – Бенара о тепловой конвекции (рис. $7.31, a$ ). Слой жидкости толщиной $h$ в поле тяжести подогревается снизу при постоянной разности температур $\Delta T=T_{1}-T_{0}$. Движение жидкости описывается уравнениями

Рис. 7.31. Конвекция Рэлея-Бенара.
$a-$ схема олыта; $\Delta T=T_{1}-T_{0}>0 ; \sigma-$ стационарная конвекция при $\Delta T>\Delta T_{c}$.
Навье–Стокса. Ограничиваясь двумерным движением ( $\partial / \partial z \equiv 0$ ), введем функцию потока $\psi(x, y, t)$, которая связана со скоростью жидкости $v(x, y, t)$ посредством формулы
\[
v=
abla \times(\check{z} \psi) .
\]

Введем также функцию $\Theta(x, y, t)$, описывающую отклонение температуры $T(x, y, t)$ от линейной зависимости по $y$ :
\[
\Theta=T-T_{1}+\frac{\Delta T}{h} y .
\]

В отсутствие конвекции $\Theta=0$. Для выбранных переменных задача сводится к двум уравнениям в частных производных [283]:
\[
\frac{\partial}{\partial t}(\Delta \psi)=\left(\frac{\partial \psi}{\partial x} \frac{\partial \Delta \psi}{\partial y}-\frac{\partial \psi}{\partial y} \frac{\partial \Delta \psi}{\partial x}\right)+v \Delta^{2} \psi-g \alpha \frac{\partial \Theta}{\partial x},
\]

\[
\frac{\partial \Theta}{\partial t}=\left(\frac{\partial \psi}{\partial x} \frac{\partial \Theta}{\partial y}-\frac{\partial \psi}{\partial y} \frac{\partial \Theta}{\partial x}\right)-\frac{\Delta T}{h} \frac{\partial \psi}{\partial x}+x \Delta \Theta .
\]

Здесь $v$ – кинематическая вязкость, $g$ – ускорение силы тяжести, $\alpha$ – коэффициент теплового расширения и $x$ – температуропроводность ${ }^{1}$ ). Примем граничные условия в виде $\Theta=\psi=\Delta \psi=0$ при $y=0$ и $y=h$, что соответствует фиксированным $T_{0}$ и $T_{1}$ и свободной поверхности жидкости. При малых $\Delta T$ имеется устойчивое равновесное состояние $\psi=\Theta=0$, соответствующее покоящейся жидкости и молекулярной теплопередаче. Еще лорд Рэлей изучал линейную устойчивость этого состояния и показал, что выше некоторого критического значения $\Delta T_{c}$ оно становится неустойчивым и в жидкости возникают циркулирующие потоки, (рис. 7.31, б):
\[
\begin{array}{l}
\psi=\psi_{0} \sin \left(\frac{\pi a x}{h}\right) \sin \left(\frac{\pi y}{h}\right), \\
\Theta=\Theta_{0} \cos \left(\frac{\pi a x}{h}\right) \sin \left(\frac{\pi y}{h}\right),
\end{array}
\]

где параметр $a$ характеризует периодичность движения по $x$. Введем безразмерное число Рэлея, характеризующее разность температур:
\[
R_{a}=\frac{g \alpha h^{3} \Delta T}{v x} .
\]

Критическое значение числа Рэлея, определяющее возникновение устойчивой конвекции, равно
\[
R_{c}=\frac{\pi^{4}\left(1-+a^{2}\right)^{3}}{a^{2}}
\]

и принимает минимальное значение $27 \pi^{4} / 4$ для $a=1 / \sqrt{2}$.
При дальнейшем увеличении $R_{a}$ выше $R_{c}$ регулярная конвекция становится линейно неустойчивой. Эксперимент показывает, что конвекция становится при этом нестационарной и нерегулярной. Для анализа этого случая, следуя Зальцману [359], разложим $\psi$ и $\Theta$ в двойной ряд Фурье по $x$ и $y$, так что коэффициенты разложения будут зависеть только от $t$. Оставляя конечное число членов, получаем представление движения в конечномерном фазовом пространстве фурье-амплитуд. Зальцман численно нашел случаи хаотического движения ${ }^{2}$ ). Лоренц [283] исследовал упрощен- $\qquad$
1) Величина $\omega=-\Delta \psi$ характеризует вращение элемента жидкости вокруг оси $z$ и называется завихренностью ( $\check{z} \omega=
abla \times \boldsymbol{v})$. – Прим. ред.
\”) См. благодарности в статье Лоренца [283].- Прим. ред.

ную систему, в которой оставлены только три фурье-амплитуды:
\[
\begin{array}{l}
\frac{a}{x\left(1+a^{2}\right)} \psi=\sqrt{2} X(t) \sin \left(\frac{\pi a x}{h}\right) \sin \left(\frac{\pi y}{h}\right), \\
\frac{\pi R_{a}}{R_{c}(\Delta T)} \Theta=\sqrt{2} Y(t) \cos \left(\frac{\pi a x}{h}\right) \sin \left(\frac{\pi y}{h}\right)-Z(t) \sin \left(\frac{2 \pi y}{h}\right) .
\end{array}
\]

Здесь $X$ – амплитуда конвективного движения, $Y$ – разность температур между восходящими и нисходящими потоками, а $Z$ отклонение вертикального профиля температуры от линейного. Подставляя (7.4.11) в уравнение (7.4.7), приходим к модели Лоренца:
\[
\begin{array}{l}
\dot{X}=-\sigma X+\sigma Y, \\
\dot{Y}=-X Z+r X-Y, \\
\dot{Z}=X Y-b Z,
\end{array}
\]

где $\sigma=v / x-$ число Прандтля, $r=R_{a} / R_{c}$ – приведенное число Рэлея, $b=4\left(1+a^{2}\right)^{-1}$, а точка означает производную по безразмерному времени $\tau=\pi^{2} h^{-2}\left(1+a^{2}\right) x t$.

Модель Лоренца и ее странный аттрактор уже рассматривались в $\S 1.5$ и выше в этой главе. Здесь же нас интересует вопрос: в какой мере эта модель представляет поведение жидкости в задаче Рэлея–Бенара? На первый взгляд обе системы очень далеки друг от друга, поскольку модель Лоренца является чрезвычайно упрощенной с ее всего лишь тремя модами для двух функций состояния жидкости $\psi$ и $\Theta$. Увеличение числа мод до пяти, семи и даже четырнадцати сохраняет некоторые черты поведения модели, включая и образование странного аттрактора. Однако переход к хаотическому движению может происходить при этом через разные последовательности бифуркаций [98] (дополнительную библиографию см. в работе [180]). Более того, численное моделирование двумерной конвекции, согласно (7.4.7), показывает отсутствие турбулентного движения ${ }^{1}$ ). В этом состоит существенное отличие от трехмерной конвекции Рэлея-Бенара, в которой турбулентность наблюдается экспериментально.

В работе [131] было показано, что двумерные решения уравнения Навье-Стокса можно описать асимптотически (при $t \rightarrow \infty$ ) конечным числом мод $N_{2}$. Этот результат был обобщен на трехмерные решения, причем $N_{3} \propto N_{2}^{3}$ [132]. Поэтому при фиксированных параметрах системы и начальных условиях конечномерные
1) Это замечание непонятно. По ловоду особенностей двумерной турбу. лентности см., например, работы [542, 543].- Прим. ред.

модели с некоторым минимальным значением $N_{2}$ (или $N_{3}$ ) представляют все физически существенные свойства реального течения ${ }^{1}$ ). Треве [412] предложил численный метод проверки anocmeриори, было ли выбранное чис.то мод достаточным. Однако пока не существует никакого метода определения необходимого числа мод априори до численного моделирования²).
7.4б. Переход к турбулентности

Обсудим теперь различные гипотезы относительно механизма перехода от регулярного течения к гидродинамической турбулентности. За исключением самой ранней гипотезы Ландау [251], все предложенные механизмы связаны с конечномерными моделями. Имеющиеся экспериментальные данные не позволяют сделать определенного выбора между моделями, поскольку в опытах часто присутствуют черты разных механизмов. Наше описание возникновения турбулентности базируется на обзорах Хеллемана [180], Отта [324] и Экмана [112 ].

Хотя явление турбулентности известно уже сотни лет, детальные измерения характеристик потока вблизи перехода были выполнены только в последнее десятилетие ${ }^{3}$ ). Возникновению турбулентного режима движения обычно предшествует возбуждение
1) Точнее конечное число мод движения полностью определяет и все остальные моды (см. также работу [545]). Физически это связано с сильным затуханием высоких мод, которые, таким образом, просто повторяют с уменьшенной амплитудой колебания основных моди имеют, в частности, тот же тип спектра (дискретный или непрерывный). Иначе можно сказать, что движение на аттракторе вообще имеет конечное число мод, однако эти моды не совпадают с невозмущенными модами приближения Галёркина. Отсюда следует также, что движение имеет только конечное число положительных показателей Ляпунова. Отметим, что конечномерность аттрактора (гипотеза Хопфа [546]) доказана только для двумерного течения [131,545], тогда как в трехмерном случае это принимается в работе [132] в качестве гипотезы. -. Прим. ред.
2) Грубые априорные оценки получены в упомянутой работе [132] и их чожно представить в виде $N_{2} \sim \mathrm{Re}^{3} \sim N_{3}^{1 / 3}$, где $\mathrm{Re}$ – число Рейнольдса. Јднако эти оценки (как и критерий в работе [412]) являются лишь достагочными и, вообще говоря, весьма неэффективными. Например, в случае днородной турбулентности они существенно превосходят естественную ценку числа степеней свободы по Ландау и Лифшицу [539]: $N_{2} \sim \mathrm{Re}^{3 / 2} \sim N_{3}^{2 / 3}$. Это хорошо согласуется также с полученной недавно строгой оценкой снизу необходимой): $N_{2} \sim \operatorname{Re}[560]$. В любом случае такие оценки указывают ға существенный дефект используемых конечномерных моделей турбулентюсти, число степеней свободы которых фиксировано, а параметры измеяяются. – Прим. ред.
3) Это краткое замечание представляет обширные многолетние исследоьания турбулентности весьма однобоко (см., например, работы $[535,536]$ ). 3 последнее время, особенно после нашумевшей работы Рюэля и Тэкенса 355 ], действительно были проведены некоторые специальные эксперименты, іногда, кстати сказать, весьма поспешно (см. примечание редактора на $\therefore 483$ ), по проверке ряда новых гипотез (см. ниже). – Прим. ред.

колебаний одной или нескольких независимых частот и их гармоник, а иногда и субгармоник. Основные эксперименты по детальному исследованию перехода к турбулентности относятся к различным вариантам конвекции Рэлея – Бенара, а также вихрей Тейлора, возникающих в слое жидкости между вращающимися цилиндрами (так называемое круговое течение Куэтта). Эти примеры мы и используем ниже при сравнении с различными гипотезами о механизме возникновения турбулентности.

Первоначальная картина возникновения турбулентности, предложенная Ландау, была основана на представлении об иерархии неустойчивостей. При увеличении некоторого параметра, например числа Рейнольдса или числа Рэлея, нелинейные колебания жидкости теряют устойчивость и появляются все новые и новые независимые частоты движения $\omega_{1}, \omega_{2}, \omega_{3} \ldots$.. При этом должно наблюдаться квазипериодическое движение с одной, двумя, тремя и т. д. основными частотами. Таким образом, мы приходим к последовательности бифуркаций Хопфа, т. е. к движению по поверхности некоторого тора возрастающей размерности. Движение выглядит все более и более сложным, однако непрерывный спектр и хаотическое движение возникают лишь при бесконечном числе бифуркаций. Модель Ландау представлена схематически в табл. 7.2.

Хотя в экспериментах и наблюдается до четырех [158] независимых частот, резкий переход к непрерывному спектру не согласуется с моделью Ландау ${ }^{1}$ ). Помимо этого, теоретически было показано (см. работу [I12]), что последовательность бифуркаций Хопфа, как и само квазипериодическое движение, не являются типичными. Судя по рассмотренным выше примерам резкого ${ }^{2}$ ) возникновения непрерывного спектра, связанного с образованием странного аттрактора, можно ожидать, что именно такой механизм $\qquad$
1) Поучительно отметить, что эта «квазипериодическая» модель турбулентности была очевидно навеяна аналогиями с квантовой механикой. Вот соответствующая цитата из оригинальной работы Ландау [251]:
«Задача об определении возможных значений «частот» $\Omega$ при данных граничных условиях движения представляет собой Eigenwert-проблему. В результате ее решения должен получиться спектр «собственных частот». .. Можно предполагать, что частоты непрерывного спектра соответствуют движениям $\boldsymbol{v}_{1}$ [возмущение], не затухающим на бесконечности [в пространстве], а частоты дискретного спектра – движениям, достаточно быстро затухающим на бесконечности (подобно тому, как это имеет место во многих других Eigenwert-проблемах) … Соответственно этому ниже мы можем рассматривать только «частоты» $\Omega$ дискретного спектра».
В квантовой механике это действительно так.- Прим. ред.
2) Представление о «резком» переходе к хаотическому движению является условным. Так, например, в механизме бифуркаций удвоения хаотический аттрактор появляется в критической точке ( $C=C_{\infty}$ ) лишь номинально, поскольку мощность хаотической компоненты движения в этой точке равна нулю и плавно возрастает в результате обратных бифуркаций (см., например, рис. 7.22).- Прим. ред.

Таблица 7.2. Модели и механизмы перехода к турбулентности

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

Рюэль и Тэкенс [355] предложили другой механизм возникновения турбулентности, согласно которому сначала происходят две последовательные бифуркации Хопфа, как и в модели Ландау, однако затем нелинейность разрушает трехчастотное движение и образуется «странный» аттрактор (табл. 7.2). По первоначальной гипотезе требовалась размерность потока не менее четырех. Предположение о неустойчивости трехчастотного аттрактора в типичном случае было позднее доказано, а минимальная размерность сокращена до трех [317] ${ }^{1}$ ).

Модель Рюэля-Тэкенса исследовалась численно на примере простого двумерного отображения [100]. Были обнаружены переходы от устойчивого фокуса к предельному циклу, затем к двухчастотному движению и, наконец, к странному аттрактору. В этой связи важно отметить, что в отличие от модели Лоренца с тремя модами в модели конвекции Рэлея-Бенара, использующей 14 мод, также обнаружен квазипериодический аттрактор на некоторой двумерной поверхности в 14-мерном фазовом пространстве [98].
1) Хотя работа [355] оказала фактически заметное влияние на развитие как теоретических, так и экспериментальных исследований по турбулентности, предложенный ее авторами конкретный механизм хаоса, связанный с резонансами высоких гармоник, представляется слишком «нежным» для такого «грубого» явления, как гидродинамическая турбулентность. Что же касается минимальной размерности потока, то этот вопрос был выяснен еще в классической работе Лоренца [283] (см. также работу [533]).- Прим. ред.

Некоторые экспериментальные данные, по-видимому, подтверждают модель Рюэля-Тэкенса. Так, в спектрах мощности появляется сначала одна, затем вторая и, возможно, третья независимая частота. На пороге появления третьей частоты внезапно возникает широкополосный шум, который свидетельствует о переходе к хаотическому движению. Экспериментально исследовались как вихри Тейлора в жидкости между вращающимися цилиндрами [125], так и конвекция Рэлея-Бенара [5]. На рис. 7.32 из популярной статьи Суинни и Голуба [396] показаны спектры скорости жидкости для течения Куэтта (слева) и для конвекции Рэлея-Бенара (справа). В обоих случаях перед переходом к непрерывному спектру наблюдается сначала одна, а затем две независимые частоты $f_{1}$ и $f_{2}$. Однако это зависит, вообще говоря, от начальных условий и иногда частоты $f_{1}$ и $f_{2}$ оказываются синхронизованными ${ }^{1}$ ). В другом эксперименте по течению Куэтта [158] наблюдались по крайней мере четыре независимые частоты. Это указывает на то, что переход к турбулентности происходит не всегда после двух бифуркаций Хопфа, как в модели Рюэля-Тэкенса.

Третья модель перехода к турбулентности, предложенная Фейгенбаумом, связана с последовательностью бифуркаций удвоения периода [122]. Переход начинается с бифуркации Хопфа из устойчивого фокуса в предельный цикл с частотой $f_{1}$. При дальнейшем
1) То есть отношение частот $f_{1} / f_{2}=m / n$ рационально (в пределах точности измерений $\sim 10^{-4}$, см., например, работу [541]). Экспериментальная независимость (несоизмеримость) частот определяется с той же точностью для $|m|,|n|<20 .-$ Пим. ред.

увеличении параметра происходят последовательные бифуркации удвоения, приводящие к периодическому движению с частотами $f_{1} / 2, f_{1} / 4, f_{1} / 8$ и т. д. Эта последовательность сходится при некотором критическом значении параметра, при котором возникает странный аттрактор (табл. 7.2). Механизм удвоения периода был описан в $\$ 7.2$ и п. 7.3а, где показан его универсальный характер. Существенно, что вблизи критической точки движение является близким к одномерному (см. п. 7.3a).

Рис. 7.32. Спектры скорости жидкости между вращающимися цилиндрами при трех значения частоты внутреннего цилиндра $f_{0}$ (слева), то же для трех значений $\Delta T$ в задаче Рэлея-Бенара (справа) (по данным работы [396]).

Модель Фейгенбаума хорошо подтверждается численными экспериментами на простых моделях. Как мы видели выше, бифуркации удвоения периода найдены и во многих динамических системах с малой размерностью, таких, как аттрактор Рё́слера, отображение Хенона, уравнение Дюффинга и др. Некоторые эксперименты по конвекции Рэлея-Бенара обнаруживают эти бифуркации, а также некоторые признаки их универсальности. Спектры скорости высокого разрешения в эксперименте Рэлея-Бенара с водой, показанные на рис. $7.33[155,157]$, демонстрируют некоторые из бифуркаций удвоения. Было проведено также сравнение экспериментальных значений амплитуд субгармоник с предсказаниями модели Фейгенбаума. Қак показано в п. 7.2б, отношение амплитуд развитых субгармоник должно быть равным $\gamma \approx 6,6$. На рис. $7.33,2$

расстояние между верхней прямой, проведенной по пикам гармоник частоты $f_{2} / 2$, и нижней прямой соответствует множителю 6,6 . Видно, что большинство пиков гармоник частоты $f_{2} / 4$ примерно ложатся на нижнюю прямую в соответствии с теорией ${ }^{1}$ ). В другом

Рис. 7.33. Бифуркации удвоения периода в эксперименте Рэлея-Бенара (по данным работы [157]).
a-в – при увеличении числа Рэлея возникают субгармоники; ¿ – хаотические полосы; $\partial$ – непрерывный спектр; $e$ – при уменьшении $R_{a}$ появляется вторая частота $f^{*}$.

эксперименте по конвекции Рэлея-Бенара в жидком гелии [263] наблюдались последовательные субгармоники $f_{1}, f_{1} / 2, f_{1} / 4, f_{1} / 8$ и $f_{1} / 16$ в спектре температуры. Отношение амплитуд последовательных субгармоник также равно примерно 6,6 (см. работу [123]).
1) Сравнение экспериментальных данных на рис. 7.33 , а для фиксированного значения параметра $C$ с законом подобия (7.2.44) для разных $C$ возможно вследствие того, что при бифуркации удвоения амплитуды прежних частот (четные гармоники на рис. $7.33,2$ ) остаются неизменными согласно (7.2.40). При сравнении необходимо, во-первых, учитывать, что масштабный фактор $\gamma$ относится к амплитудам гармоник, а во-вторых, брать только нечетные до бифуркации гармоники, т. е. только второй и шестой пики на рис. 7.33 , 2. В результате среднее $\gamma \approx 4,4$ (см. примечание редактора на с. 440).- Прим. ред.

Однако, как мы видели, другие эксперименты не подтверждают эти свойства ${ }^{1}$ ). Можно предположить, что причиной этого является внешний шум, который нарушает тонкую структуру спектра, и в зависимости от условий эксперимента субгармоники то появляются, то пропадают. Такое поведение действительно наблюдалось в численном моделировании Кратчфилдом и Хьюберманом [96].

Рассмотрим, наконец, четвертый механизм возникновения турбулентности, лежащий в основе модели Помо и Манневиля [293] и связанный с переходом к хаотическому движению с перемежаемостью. В этой модели при увеличении некоторого параметра периодическая траектория непосредственно превращается в хаотическую с перемежаемостью в результате обратной тангенциальной бифуркации.

Следуя Экману [112], покажем, каким образом это происходит на примере потери устойчивости предельного цикла периода 3 для одномерного квадратичного отображения (7.2.4). В п. 7.2 (см. рис. 7.15 и пояснения к нему) было показано, что при $C=C_{0}^{(3)}=$ $=(1-\sqrt{8}) / 2$ в результате тангенциальной бифуркации рождаются устойчивая и неустойчивая траектории периода 3. Они сохраняются в небольшом интервале значений $C<C_{0}^{(3)}$. При $C=C_{0}^{(3)}$ они сливаются и исчезают, рождая так называемое хаотическое движение с перемежаемостью. Под этим понимается следующее: можно показать, что хотя движение является хаотическим, однако при малой разности $C-C_{0}^{(3)} \sim \varepsilon$ типичная траектория будет находиться в окрестности периодической траектории в течение $\sim \varepsilon^{-1 / 2}$ итераций. При этом периодическая траектория
\[
x_{3 i}=f\left(f\left(f\left(x_{3}\right)\right)\right), \quad i=1,2,3,
\]

соответствует бифуркационному значению $C=C_{0}^{(3)}$. В течение этого времени траектория выглядит как периодическая, а затем, после выхода из этой области,- как хаотическая. Поэтому в модели Помо и Манневиля нерегулярное движение продолжается в течение некоторого случайного промежутка времени, сменяясь периодическим со средней длительностью $\sim\left|\mu-\mu_{c}\right|^{1 / 2}$, где $\mu_{c}-$ бифуркационное значение некоторого параметра $\mu$.

Численное моделирование квадратичного отображения подтверждает такое поведение [112]. Оно существует и для модели Лоренца в некотором интервале параметров [293]. Подобный переход к турбулентности наблюдался во многих экспериментах, включая конвекцию Рэлея-Бенара [273] и так называемую химическую турбулентность [352] (см. также дополнение А). Однако в этих случаях перемежаемость связана с переходом между
1) Речь идет, по-видимому, о данных на рис. 7.32. – Прим. ред.

двумя квазиустойчивыми состояниями системы и остается неясным, описывается ли такое поведение моделью Помо – Манневиля ${ }^{1}$ ).

Следует подчеркнуть, что единого механизма перехода к турбулентности не существует ${ }^{2}$ ). Если, например, система обнаруживает бифуркации удвоения периода, то можно предсказать зависимость этих бифуркаций от параметра. Однако пока неизвестно, каким образом можно узнать заранее, в какой системе это будет происходить. Отметим также, что все эти модели описывают только возникновение турбулентности и ничего не говорят о свойствах развитой турбулентности (см., например, [295]).
1) Хотя термин «перемежаемость» появился недавно, подобные процессы рассматривались уже довольно давно под более удачным, на наш взгляд, названием – «структурная турбулентность» (см., например, работы $[537,538]$, где имеется подробная библиография). В частности, появление структур в хаотическом режиме простой диссипативной модели описано в работе [530]. Такие флуктуирующие структуры часто встречаются и в гамильтоновых системах. Типичный пример – движение в узком стохастическом слое сепаратрисы маятника (резонанса) (§3.5). Здесь имеются три структуры (вращение в двух направлениях и колебания), между которыми происходят случайные переходы. – Прим. ред.
${ }^{2}$ ) См., например, работу [541].- Прим. ред.

Categories

1
Оглавление
email@scask.ru