Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO
Рассмотрим теперь связь исследованных выше конечномерных хаотических потоков с движением сплошной среды («жидкости») с бесконечным числом степеней свободы. Сначала мы обсудим вопрос, в каком смысле конечномерные модели аппроксимируют жидкость, а затем опишем некоторые механизмы перехода от регулярного движения к турбулентному. Будем описывать поведение жидкости с помощью дифференциальных уравнений в частных производных вида где $\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$ : где Подставляя (7.4.2) в (7.4.1) и используя ортогональность функций $e^{i k \cdot x}$, получаем уравнения движения в виде ${ }^{1}$ ) Если оставить в сумме (7.4.2) $N$ «наиболее существенных» мод, то задача сведется к конечному числу ( $M N$ ) обыкновенных диффе- ренциальных уравнений первого порядка, описывающих временну́ю эволюцию этих мод. Такой метод называется приближением Галё́ркина. Рассмотрим в качестве примера задачу Рэлея — Бенара о тепловой конвекции (рис. $7.31, a$ ). Слой жидкости толщиной $h$ в поле тяжести подогревается снизу при постоянной разности температур $\Delta T=T_{1}-T_{0}$. Движение жидкости описывается уравнениями Рис. 7.31. Конвекция Рэлея-Бенара. Введем также функцию $\Theta(x, y, t)$, описывающую отклонение температуры $T(x, y, t)$ от линейной зависимости по $y$ : В отсутствие конвекции $\Theta=0$. Для выбранных переменных задача сводится к двум уравнениям в частных производных [283]: \[ Здесь $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, б): где параметр $a$ характеризует периодичность движения по $x$. Введем безразмерное число Рэлея, характеризующее разность температур: Критическое значение числа Рэлея, определяющее возникновение устойчивой конвекции, равно и принимает минимальное значение $27 \pi^{4} / 4$ для $a=1 / \sqrt{2}$. ную систему, в которой оставлены только три фурье-амплитуды: Здесь $X$ — амплитуда конвективного движения, $Y$ — разность температур между восходящими и нисходящими потоками, а $Z$ отклонение вертикального профиля температуры от линейного. Подставляя (7.4.11) в уравнение (7.4.7), приходим к модели Лоренца: где $\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]. Поэтому при фиксированных параметрах системы и начальных условиях конечномерные модели с некоторым минимальным значением $N_{2}$ (или $N_{3}$ ) представляют все физически существенные свойства реального течения ${ }^{1}$ ). Треве [412] предложил численный метод проверки anocmeриори, было ли выбранное чис.то мод достаточным. Однако пока не существует никакого метода определения необходимого числа мод априори до численного моделирования²). Обсудим теперь различные гипотезы относительно механизма перехода от регулярного течения к гидродинамической турбулентности. За исключением самой ранней гипотезы Ландау [251], все предложенные механизмы связаны с конечномерными моделями. Имеющиеся экспериментальные данные не позволяют сделать определенного выбора между моделями, поскольку в опытах часто присутствуют черты разных механизмов. Наше описание возникновения турбулентности базируется на обзорах Хеллемана [180], Отта [324] и Экмана [112 ]. Хотя явление турбулентности известно уже сотни лет, детальные измерения характеристик потока вблизи перехода были выполнены только в последнее десятилетие ${ }^{3}$ ). Возникновению турбулентного режима движения обычно предшествует возбуждение колебаний одной или нескольких независимых частот и их гармоник, а иногда и субгармоник. Основные эксперименты по детальному исследованию перехода к турбулентности относятся к различным вариантам конвекции Рэлея — Бенара, а также вихрей Тейлора, возникающих в слое жидкости между вращающимися цилиндрами (так называемое круговое течение Куэтта). Эти примеры мы и используем ниже при сравнении с различными гипотезами о механизме возникновения турбулентности. Первоначальная картина возникновения турбулентности, предложенная Ландау, была основана на представлении об иерархии неустойчивостей. При увеличении некоторого параметра, например числа Рейнольдса или числа Рэлея, нелинейные колебания жидкости теряют устойчивость и появляются все новые и новые независимые частоты движения $\omega_{1}, \omega_{2}, \omega_{3} \ldots$.. При этом должно наблюдаться квазипериодическое движение с одной, двумя, тремя и т. д. основными частотами. Таким образом, мы приходим к последовательности бифуркаций Хопфа, т. е. к движению по поверхности некоторого тора возрастающей размерности. Движение выглядит все более и более сложным, однако непрерывный спектр и хаотическое движение возникают лишь при бесконечном числе бифуркаций. Модель Ландау представлена схематически в табл. 7.2. Хотя в экспериментах и наблюдается до четырех [158] независимых частот, резкий переход к непрерывному спектру не согласуется с моделью Ландау ${ }^{1}$ ). Помимо этого, теоретически было показано (см. работу [I12]), что последовательность бифуркаций Хопфа, как и само квазипериодическое движение, не являются типичными. Судя по рассмотренным выше примерам резкого ${ }^{2}$ ) возникновения непрерывного спектра, связанного с образованием странного аттрактора, можно ожидать, что именно такой механизм $\qquad$ Таблица 7.2. Модели и механизмы перехода к турбулентности и послужит основой построения современной теории турбулентности, как это и оказывается в действительности. Рюэль и Тэкенс [355] предложили другой механизм возникновения турбулентности, согласно которому сначала происходят две последовательные бифуркации Хопфа, как и в модели Ландау, однако затем нелинейность разрушает трехчастотное движение и образуется «странный» аттрактор (табл. 7.2). По первоначальной гипотезе требовалась размерность потока не менее четырех. Предположение о неустойчивости трехчастотного аттрактора в типичном случае было позднее доказано, а минимальная размерность сокращена до трех [317] ${ }^{1}$ ). Модель Рюэля-Тэкенса исследовалась численно на примере простого двумерного отображения [100]. Были обнаружены переходы от устойчивого фокуса к предельному циклу, затем к двухчастотному движению и, наконец, к странному аттрактору. В этой связи важно отметить, что в отличие от модели Лоренца с тремя модами в модели конвекции Рэлея-Бенара, использующей 14 мод, также обнаружен квазипериодический аттрактор на некоторой двумерной поверхности в 14-мерном фазовом пространстве [98]. Некоторые экспериментальные данные, по-видимому, подтверждают модель Рюэля-Тэкенса. Так, в спектрах мощности появляется сначала одна, затем вторая и, возможно, третья независимая частота. На пороге появления третьей частоты внезапно возникает широкополосный шум, который свидетельствует о переходе к хаотическому движению. Экспериментально исследовались как вихри Тейлора в жидкости между вращающимися цилиндрами [125], так и конвекция Рэлея-Бенара [5]. На рис. 7.32 из популярной статьи Суинни и Голуба [396] показаны спектры скорости жидкости для течения Куэтта (слева) и для конвекции Рэлея-Бенара (справа). В обоих случаях перед переходом к непрерывному спектру наблюдается сначала одна, а затем две независимые частоты $f_{1}$ и $f_{2}$. Однако это зависит, вообще говоря, от начальных условий и иногда частоты $f_{1}$ и $f_{2}$ оказываются синхронизованными ${ }^{1}$ ). В другом эксперименте по течению Куэтта [158] наблюдались по крайней мере четыре независимые частоты. Это указывает на то, что переход к турбулентности происходит не всегда после двух бифуркаций Хопфа, как в модели Рюэля-Тэкенса. Третья модель перехода к турбулентности, предложенная Фейгенбаумом, связана с последовательностью бифуркаций удвоения периода [122]. Переход начинается с бифуркации Хопфа из устойчивого фокуса в предельный цикл с частотой $f_{1}$. При дальнейшем увеличении параметра происходят последовательные бифуркации удвоения, приводящие к периодическому движению с частотами $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]). эксперименте по конвекции Рэлея-Бенара в жидком гелии [263] наблюдались последовательные субгармоники $f_{1}, f_{1} / 2, f_{1} / 4, f_{1} / 8$ и $f_{1} / 16$ в спектре температуры. Отношение амплитуд последовательных субгармоник также равно примерно 6,6 (см. работу [123]). Однако, как мы видели, другие эксперименты не подтверждают эти свойства ${ }^{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}$ итераций. При этом периодическая траектория соответствует бифуркационному значению $C=C_{0}^{(3)}$. В течение этого времени траектория выглядит как периодическая, а затем, после выхода из этой области,- как хаотическая. Поэтому в модели Помо и Манневиля нерегулярное движение продолжается в течение некоторого случайного промежутка времени, сменяясь периодическим со средней длительностью $\sim\left|\mu-\mu_{c}\right|^{1 / 2}$, где $\mu_{c}-$ бифуркационное значение некоторого параметра $\mu$. Численное моделирование квадратичного отображения подтверждает такое поведение [112]. Оно существует и для модели Лоренца в некотором интервале параметров [293]. Подобный переход к турбулентности наблюдался во многих экспериментах, включая конвекцию Рэлея-Бенара [273] и так называемую химическую турбулентность [352] (см. также дополнение А). Однако в этих случаях перемежаемость связана с переходом между двумя квазиустойчивыми состояниями системы и остается неясным, описывается ли такое поведение моделью Помо — Манневиля ${ }^{1}$ ). Следует подчеркнуть, что единого механизма перехода к турбулентности не существует ${ }^{2}$ ). Если, например, система обнаруживает бифуркации удвоения периода, то можно предсказать зависимость этих бифуркаций от параметра. Однако пока неизвестно, каким образом можно узнать заранее, в какой системе это будет происходить. Отметим также, что все эти модели описывают только возникновение турбулентности и ничего не говорят о свойствах развитой турбулентности (см., например, [295]).
|
1 |
Оглавление
|