История теоретического изучения течения жидкости в картинках и примерах

В реферате, который вы сейчас читаете, предпринята попытка осветить некоторые

этапы развития теории турбулентности. Так как задача описания турбулентности

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

науки в которую развилась эта проблема с течением времени не представляется

возможным. Поэтому автор и не пытался ``объять необъятное'', а поставил себе

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

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

более или менее подрбно.

\beginsection{Как всё начиналось}

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

прусский император Фридрих Великий пригласил Леонарда Эйлера работать в

Потсдам. Согласно популярной

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

одной из задач Эйлера было сооружение фонтана. Как истинный теоретик

он начал свою работу с вывода законов движения жидкости.

В 1755 году в письме Ньютону он

приводит уравнение, описывающее движение жидкости, которое в современных

обозначениях для случая постоянной плотности выглядит так:

$$ {\partial \u\arg \over \partial t} + \u\arg \scal \nabla\, \u\arg

   = -\nabla p\arg .$$

Где $\u\arg$ и $p\arg$ --- скорость жидкости и давление в пространственной

точке~$\r$ в момент времени~$t$, а точка в $\u\arg\scal\nabla$ означает

скалярное произведение. Левая часть уравнения Эйлера --- ускорение

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

этот элемент, которая порождается неоднородностью распределения давления в

жидкости. Таким образом уравнение Эйлера --- это фактически уравнение Ньютона

для элемента жидкости.

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

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

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

Дело в том, что при выводе уравнения Эйлера было упущено немаловажное, как

выяснилось, соображение о наличии

вязкого трения, т.~е. диссипации энергии при трении соседних элементов

жидкости друг о друга. Член учитывающий этот процесс был добавлен в уравнение

Навье в 1827 и Стоксом в 1845. Получившееся уравнение известно как уравнение

Навье-Стокса:

$$ {\partial \u\arg \over \partial t} + \u\arg \scal \nabla\, \u\arg

   = -\nabla p\arg +\nu \nabla^2 \u\arg .$$

Где $\nu$ --- кинематическая вязкость, для воды и воздуха при комнатной

температуре она соответственно равна 0.01 и 0.15 см${}^2$/сек.

Если бы член с вязкостью отсутствовал, то кинетическая энергия $\u^2/2$

сохранялась бы. При наличии этого члена кинетическая энергия рассеивается

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

подкачки энергии извне в конце концов всякое движение в жидкости прекратится.

\beginsection{Простые оценки}

Уже простейшие оценки на решение уравнения Навье-Стокса могут оказаться

очень далеки от реальности. Например, можно оценить скорость течения в великих

реках, таких как Нил и Волга, длина которых достигает тысяч километров, а

перепад высот между истоком и устьем составляет сотни метров. Оценим характерный

угол наклона ложа реки $\alpha$ как $10^{-4}$ радиан, а характерную глубину

$L$ как 10 метров. Приравнивая силу тяжести $\alpha g$ ($g\simeq 10^3$

см/сек${}^2$) к силе вязкого трения $\nu\, d^2u/dz^2 \sim \nu u/L^2$, находим

$u\sim 10^7$ см/сек. Конечно же это абсурд, возможно к большому сожалению

сплавщиков леса. Эта оценка противоречит даже закону сохранения энергии.

Потенциальная энергия накопленная водой на перепаде высот $H$ равна

$\rho g H$. Для таких рек как Нил и Волга $H\sim 5\times 10^4$ см и,

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

$u\sim \sqrt{2\rho g H}\simeq 10^4$ см/сек. Эта оценка расходится с реальностью

на два порядка. Объяснение такого несоответствия было предложено Ренольдсом,

который заметил важность безразмерного отношения нелинейного члена уравнения

Навье-Стокса к вязкому. Если на масштабе $L$ характерное изменение скорости

$U$, то нелинейный член может быть оценен как $U^2/L$, а вязкий --- как

$\nu U/L^2$. Отношение этих членов называется числом Рейнольдса (Re) и равно

$UL/\nu$. Если Re~$\ll 1$, то можно пренебречь нелинейным членом по сравнению

с вязким. В этом случае уравнение Навье-Стокса становится линейным и в

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

без завихрений.

Такое течение называют ламинарным. Если же реализуется обратная ситуация

(Re~$\gg 1$),то решение

уравнения Навье-Стокса становится неустойчивым и приобретает сложный

завихряющийся характер. Такое течение называют турбулентным. Как правило,

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

Re~$\sim 10^7$.

Таким образом мы приходим к существованию трёх областей, где течение жидкости

происходит качественно разным образом. Первая область -- Re~$\ll 1$ --

область ламинарности. Вторая область -- Re~$\sim 1$ -- область зарождения

турбулентности. И третья область -- Re~$\gg 1$ -- область развитой

турбулентности.

Далее мы сосредоточим наше внимание на рассмотрении последних двух областей.

\beginsection{Несколько картинок}

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

энергетическим спектром, дающим представление о распределении кинетической

энергии по частотам движений. Спектральная функция $P(\omega)$ определяется

как квадрат фурье-образа поля скоростей:

$$ \u(\omega)={1\over T} \int\limits_0^T dt\, e^{2\pi i \omega t} \u(t),$$

$$ P(\omega)=|\u(\omega)|^2.$$

Удобно так же пользоваться определёнными геометрическими образами, а именно

пространством состояний жидкости, каждая точка которого отвечает

опре\-делён\-ному распределению скоростей в этой жидкости. Состояниям в близкие

моменты времени соответствуют близкие точки. (Это вообще говоря

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

быть заменено конечномерным --- см. далее).

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

энергетического спектра и в пространстве состояний.

Если поместить какое-либо тело в поток жидкости, например, опору моста в русло

реки, то при очень малых скоростях жидкость течёт ламинарно (Рис.~1).

\Рис{1. Re $=10^{-2}$}{r1.bmp}

Такое течение является стационарным, т.~е. скорость в любой точке

пространства не зависит от времени. Следовательно вся энергия в спектре

сосредоточена на нулевой частоте. В пространстве состояний такое течение

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

т.~е. если начальное течение соответствовало другой точке в пространстве

состояний, то в пределе $t \rightarrow \infty$ любое распределение скоростей

будет стремится к устойчивому. (Строго говоря не любое, а любое из области

притяжения устойчивой траектории).

С ростом скорости в потоке образуются вихри, однако картина продолжает

остоваться стационарной (Рис.~2).

\Рис{2. Re $= 20$}{r2.bmp}

Так как поле скоростей по прежнему стационарно, то никаких изменений

относительно ламинарного течения в спектре не произойдёт. В пространстве

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

точкой, однако, её положение изменится.

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

Возникает нестационарное течение, которое, например, можно наблюдать с моста.

При этом скорость, измеренная в некоторой точке вниз по потоку за мостом,

оказывается периодической функцией времени (Рис.~3).

\Рис{3. Re $=100$}{r3.bmp}

В такой ситуации происходит качественное изменение как энергетического спектра,

так и траектории системы в пространстве состояний. В спектре появляется новая

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

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

турбулентности --- сценарии Ландау --- предполагалось, что по мере увеличения

числа Рейнольдса в системе будет возникать всё больше новых частот. Траектория

системы будет усложнятся:

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

в трёхмерный и далее бесконечный каскад новых торов. Однако, сейчас не

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

другой и данный сценарий крайне маловероятен. Один из сценариев имеющих

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

При ещё большем возрастании числа Рейнольдса (Re) крупные вихри начинают

порождать неупорядоченные внутренние вихри.

В зависимости скорости от времени кроме периодической

компоненты, появляются так же и нерегулярные отклонения.

Спектр представляет собой пики основных частот на фоне сравнительно

малоинтенсивного ``белого шума''. Траектория системы начинает размываться.

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

тора. В этом случае мы наблюдаем следующую картину потока (Рис.~4).

\Рис{4. Re $=10^4$}{r4.bmp}

Если число Рейнольдса возрастёт ещё больше, то возникает чрезвычайно сложное

поле скоростей, и зависимость $\u(t)$, как впрочем и траектория системы,

становится совершенно хаотической.

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

Из спектра исчезают пики частот и возрастает интенсивность шума. Шум практически

равномерно распределён в довольно широком интервале частот.

Картина потока соответствующая такой ситуации изображена на рисунке~5.

\Рис{5. Re $=10^6$}{r5.bmp}

Полной теории, исчерпывающим образом объясняющей возникновение турбулентности

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

Был предложен, однако, ряд возможных сценариев турбулизации течения. Эти

сценарии по большей части основаны на компьютерном исследовании модельных

систем и частично подтверждаются реальными гидродинамическими экспериментами.

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

рассмотрим только один, правда, довольно подробно.

\beginsection{Зарождение турбулентности}

Методы исследования хаотизации гидродинамических систем основаны на том,

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

именно когда в системе появилось всего несколько частот, можно

изучать как устойчивость диссипативной механической системы с конечным

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

например, амплитуды различных фурье-компонент скорости течения.

Математически это сводится к исследовнию эволюции системы, описываемой

эволюционными уравнениями вида

$$\dot {\bf x}(t)={\bf F}({\bf x}),$$

где ${\bf x}(t)$ --- вектор в пространстве $n$ величин

$x^{(1)},x^{(2)},\ldots ,x^{(n)}$, описывающих систему; функция ${\bf F}$

зависит от параметра изменение которого может приводить к изменению характера

движения. Для диссипативной системы дивергенция вектора $\bf \dot x$ в

$\bf x$-пространстве отрицательна, чем выражается сокращение объёмов

$\bf x$-пространства при

движении\footnote{$^{1)}$}{По математической терминологии функцию $\bf F$ называют

векторным полем системы. Если оно не зависит явно от времени (как в

рассматриваемом случае), систему называют автономной.}:

$${\rm div\,} \dot {\bf x}(t)={\rm div\,} {\bf F}=

   \partial F^{(i)}/\partial x^{(i)}<0. $$

Рассмотрим возможные варианты поведения системы описываемой данными

уравнениями. До сих пор молчаливо подразумевалось, что при потере устойчивости

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

даижение. Логически же это вовсе не обязательно. Ограниченность амплитуд

пульсаций скорости обеспечивает лишь ограниченность объёма пространства

состояний, внутри которого располагаются траектории, соответствующие

установившемуся режиму течения вязкой жидкости, но как выглядит картина

траекторий в этом объёме априори ничего сказать нельзя. Исключением

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

В этом случае

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

реализоваться только два типа устойчивых траекторий --- точка и цикл. Это

является следствием собственно двумерия системы и единственности решения

задачи Коши эволюционных уравнений. В случае же большего числа переменных

траектории могут стремиться к предельному циклу или к незамкнутой намотке на

торе (соответственно образам периодического или квазипериодического движений),

но могут вести себя и совершенно по-иному --- сложно и запутанно. Именно эта

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

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

Представить себе сложное и запутанное поведение траектории внутри

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

что все траектории в нём неустойчивы. Среди них могут быть не только неустойчивые

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

не выходя из неё. Неустойчивость означает, что две сколь угодно близкие точки

пространства состояний, передвигаясь в дальнейшем по проходящим через них

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

к одной траектории: ввиду ограниченности области незамкнутая траектория может

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

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

Эта картина имеет ещё и другой аспект --- чувствительная зависимость течения

от малого изменеия начальных условий. Если движение устойчиво, то малая

неточность в задании начальных условий приведёт лишь к аналогичной неточности

в определении конечного состояния. Если же движение неустойчиво, то исходная

неточность со временем нарастает и дальнейшее состояние системы уже

невозможно предвидеть ({\it Н.~С.~Крылов, 1944; M.~Born, 1952\/}).

Притягивающее множество неустойчивых траекторий в пространстве состояний

диссипативной системы действительно может существовать ({\it E. Lorenz,

1963\/}); его принято называть \it стохастическим\/ \rm или \it странным

аттрактором\/\rm\footnote{$^{2)}$}{В отличие от обычных аттракторов (устойчивые

предельные циклы, предельные точки и т. п.); название аттрактора ``странный''

связано со сложностью его структуры. В физической литературе термином

``странный аттрактор'' обозначают и более сложные притягивающие множества,

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

малыми областями притяжения, что ни в физическом, ни в численном

экспериментах их нельзя обнаружить.}.

Он был обнаружен Лоренцем при изучении

системы трёх нелинейных дифференциальных уравнений, приближённо описывающих

некоторые гидродинамические течения.

Странный аттрактор может появиться уже после нескольких бифуркаций возникновения

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

квазипериодический режим (незамкнутая обмотка на торе), создав на торе

странный аттрактор ({\it D.~Ruelle, F.~Takens, 1971}). Это, однако, не может

произойти на второй (начиная с разрушения стационарного режима) бифуркации.

При этой бифуркации появляется незамкнутая обмотка на двумерном торе. Учёт

малой нелинейности не разрушает тора, так чт странный аттрактор должен был бы

быть расположен на нём. Но на двумерной поверхности невозможно существование

притягивающего множества неустойчивых траекторий. Дело в том, что траектории в

пространстве состояний не могут пересекаться друг с другом (или сами с собой);

это противоречило бы причинности поведения классических систем: состояние

системы в каждый момент времени однозначно определяет её поведение в

следующие моменты. На двумерной поверхности невозможность пересечений настолько

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

Но уже на третьей бифуркации возникновение странного аттрактора становится

возможным (хотя и не обязательным!). Такой аттрактор, приходящий на смену

трёхчастотному квазипериодическому режиму, расположен на трёхмерном торе

({\it S.~Newhouse, D.~Ruelle, F.~Takens, 1978\/}).

Принадлежащие странному аттрактору сложные, запутанные траектории расположены

в ограниченном объёме пространства состояний. Классификация возможных типов

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

задачах, в настоящее время неизвестна; неясны даже критерии, на которых

должна была бы основываться такая классификация. Существующие знания о

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

возникающих при компьютерном решении модельных систем обыкновенных

дифференциальных уравнений, довольно далёких от реальных гидродинамических

уравнений. О структуре странного аттрактора можно, однако, высказать некоторые

общие суждения, следующие уже из неустойчивости (седлового типа) траекторий

и диссипативности системы.

Для наглядности будем говорить о трёхмерном пространстве состояний и представлять

себе аттрактор расположенный внутри двумерного тора.

Рассмотрим пучок траекторий на пути к аттрактору (ими описываются переходные

режимы движения жидкости, ведущие к установлению ``стационарной'' турбулентности).

В поперечном сечении пучка траектории (точнее --- их следы) заполняют

определённую площадь; проследим за изменением величины и формы этой площади

вдоль пучка. Учтём, что элемент объёма в окрестности седловой траектории в

одном из (поперечных) направлений растягивается, а в другом --- сжимается;

ввиду диссипативности системы сжатие сильнее, чем растяжение --- объёмы

должны уменьшаться. По ходу траекторий эти направления должны меняться ---

в противном случае траектори ушли бы слишком далеко (что означало бы слишком

большое изменение скорости жидкости). Всё это приведёт к тому, что сечение

пучка уменьшится по площади и приобретёт сплющенную, и в то же время изогнутую

форму. Но этот процесс должен происходить не только с сечением пучка в целом,

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

систему вложенных друг в друга полос, разделённых пустотами. С течением

времени (т.~е. вдоль пучка траекторий) число полос быстро возрастсет, а их

ширины убывают. Возникающий в пределе $t\longrightarrow\infty $ аттрактор

представляет собой несчётное множество бесконечного числа не касающихся друг

друга слоёв --- поверхностей, на которых располагаются седловые траектории

(своими притягивающими направлениями обращённые ``наружу'' аттрактора). Своими

боковыми сторонами и своими концами эти слои сложным образом соединяются друг

с другом; каждая из принадлежащих аттрактору траекторий блуждает по всем слоям

и по прошествии достаточно большого времени пройдёт достаточно близко к любой

точке аттрактора (свойство {\it эргодичности\/}). Общий объём слоёв и общая

площадь их сечений равны нулю.

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

относятся к категории {\it канторовых\/}. Именно канторовость структуры

следует считатть наиболее характерным свойством аттрактора и в более общем

случае $n$-мерного $(n>3)$ пространства состояний.

 Объём странного аттрактора

в своём пространстве состояний всегда равен нулю. Он может, однако, быть

ненулевым в другом пространстве --- меньшей размерности. Последнее определяется

следующим образом. Разобьём всё $n$-мерное пространство на малые кубики с

длиной ребра $\varepsilon$ и объёмом $\varepsilon^n $. Пусть $N(\varepsilon )$ ---

минимальное число кубиков, совокупность которых полностью покрывает аттрактор.

Определим размерность $D$ аттрактора как предел

$$D=\lim_{\varepsilon\to 0}{\ln N(\varepsilon)\over \ln (1/\varepsilon)}.$$

Эта величина известна в математике как

{\it хаусдорфова\/} (или {\it фрактальная\/}) размерность.

Существование этого предела означает конечность объёма аттрактора в $D$-мерном

пространстве: при малом $\varepsilon$ имеем $N(\varepsilon)\approx V\varepsilon^{-D}$

(где $V$ --- постоянная), откуда видно, что $N(\varepsilon )$ можно расматривать

как число $D$-мерных кубиков, покрывающих в $D$-мерном пространстве объём $V$.

Определённая согласно предыдущей формуле размерность не может, очевидно, превышать полную

размерность $n$ пространства состояний, но может быть меньше его и, в отличие

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

множеств\footnote{$^{3)}$} {Покрывающие

множество $n$-мерные кубики могут оказаться ``почти пустыми''; именно поэтому

может быть $D<n$. Для обычных множеств это определение даёт очевидные результаты.

Так, для множества $N$ изолированных точек имеем $N(\varepsilon )= N$  и $D=0$;

для отрезка $L$ линии: $N(\varepsilon)=L/\varepsilon ,\ D=1;$ для площадки $S$

двумерной поверхности: $N(\varepsilon )=S/\varepsilon^2 ,\ D=2 ,$ и т.д.}.

\beginsection{Переход к турбулентности путём удвоения периодов}

При изучении эволюции возущения периодического движения жидкости удобно ввести

{\it мультипликатор\/} $\mu$ периодического движения, как отношение амплитуды

возмущения в момент времени $t$ к амплитуде возмущения в момент времени $t-T$,

где $T$ --- период невозмущённого движения. То есть за период амплитуда

возмущения изменяется в $\mu$ раз. Периодическому движению непрерывной среды

(жидкости) соответствует бесконечное множество мультипликаторов, отвечающих

бесконечному числу возможных независимых возмущений. Когда модули

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

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

равными 1, т.~е. в комплексной плоскости значения $\mu$ пересекают единичную

окружность. Ввиду вещественности уравнений проходить через эту окружность

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

оставаясь вещественными, т.~е. в точках $+1$ или $-1.$ Потеря устойчивости

периодическим движением сопровождается определённой качественной перестройкой

поведения траекторий в пространстве состояний в окрестности ставшего

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

{\it бифуркацией\/}. Характер бифуркации в значительной степени определяется

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

Рассмотрим потерю устойчивости периодическим движением путём похождения

мультипликатора через значение $-1$ или $+1.$

В $n$-мерном пространстве состояний $n-1$ мультипликаторов определяют поведение

траекторий в $n-1$ различных направлений в окрестности рассматриваемой

периодической траектории (отличных от направления касательной в каждой точке

самой этой траектории). Пусть близкий к $\pm 1$ мультипликатор отвечает

некоторому $l$-му направлению. Остальные $n-2$ мультипликаторов малы по модулю;

 поэтому по соответствующим им $n-2$ направлениям все траектории будут со

временем прижиматься к некоторой двумерной поверхности (назовём её $\Sigma$),

которой принадлежит $l$-ое направление и направление указанных касательных.

Можно сказать, что в окрестности предельного цикла пространство состояний при

$t\longrightarrow\infty $ оказывается почти двумерным (строго двумерным оно быть

не может --- траектории могут располагаться по обе стороны $\Sigma$ и переходить

с одной стороны поверхности на другую). Разрежем поток траекторий вблизи

$\Sigma$ некоторой секущей поверхностью $\sigma.$ Каждая траектория, повторно

пересекая $\sigma,$ ставит в соответствие исходной точке пересечения

(назовём её $x_j$) точку пересечения в момент следующего возврата $x_{j+1}.$

Связь $x_{j+1}=f(x_j;R)$ называют {\it отображением Пуанкаре\/} или

{\it отображением последования\/}; она зависит от параметра $R$ (в данном случае

--- числа Рейнольдса\footnote{$^{4)}$}{Или числа Рэлея, если речь идёт о тепловой

конвекции.}), значение которого определяет степень близости к бифуркации ---

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

прижаты к поверхности $\Sigma$, множество точек пересечения поверхности $\sigma$

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

аппроксимировать линией; отображение Пуанкаре станет одномерным преобразованием

$$

x_{j+1}=f(x_j;R),\eqno(1)

$$

причём $x$ будет просто координатой на указанной линии

\footnote{$^{5)}$} {Обозначение $x$ не имеет, разумеется, нияего общего с

координатой в физическом пространстве.}. Дискретная переменная ${j}$ играет

роль времени, измеряемого в единицах периода движения.

Отображение Пуанкаре  даёт альтернативный способ определения характера течения

вблизи бифуркации. Самому периодическому движению отвечает {\it неподвижная

точка\/} преобразования (1) --- значение $x_j=x_*$, не меняющееся при отображении, т.~е.

для которого $x_{j+1}=x_j$. Роль мультипликатора играет производная

$\mu=dx_{j+1}/dx_j$, взятая в точке $x_j=x_*$. Точки $x_j=x_* +\xi$ в

окрестности $x_*$ в результате отображения переходят в

$x_{j+1}\approx x_* +\mu\xi$. Неподвижная точка устойчива (и является

аттрактором отображения), если $|\mu|<1$: повторно применяя ({\it итерируя\/})

отображение и начав с какой-либо точки в окрестности точки $x_*$, мы будем

асимптотически приближаться к последней (по закону $|\mu|^r$, где $r$ ---

число итераций). Напротив при $\mu >1$ неподвижная точка неустойчива.

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

мультипликатора через $-1$. Равенство $\mu=-1$ означает, что начальное возмущение

через интервал времени $T_0$ меняет знак, не изменяясь по абсолютной величине:

ещё через период $T_0$ возмущение перейдёт само в себя. Таким образом, при

переходе $\mu$ через значение $-1$ в окрестности предельного цикла с периодом

$T_0$ возникает новый предельный цикл с периодом $2T_0$ --- {\it бифуркация

удвоения периода\/}.% (Рис.~6)

Если принять условно неподвижную точку отображения Пуанкаре за точку $x=0,$

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

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

$$

x_{j+1}=\, -[1+(R-R_1)]x_j+x_j^2+\beta x_j^3,\eqno(2)

$$

где $\beta >0 $\footnote {$^{6)}$}{Коэффициент при $R-R_1$ может быть превращён

в единицу соответствующим переопределением $R$, а коэффициент при $x_j^2$

обращён в $+1$ переопределением $x_j$ (что и предполагается в (2)).}.

При $R<R_1$ неподвижная точка $x_*=0$ устойчива, а при $R>R_1$ --- неустойчива.

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

отображение (2) дважды, т.~е. рассмотреть его за два шага (две единицы времени) и

определить неподвижные точки вновь полученного отображения; если они существуют

и устойчивы, то они и отвечают циклу удвоенного периода.

Двукратная итерация преобразования (2) приводит (с нужной точностью по

малым величинам $x_j$ и $R-R_1$) к отображению

$$

x_{j+2} = x_j+2(R-R_1)x_j- 2(1+\beta)x_j^3.\eqno(3)

$$

Оно всегда имеет неподвижную точку $x_*=0.$ При $R<R_1$ эта точка единственна

и устойчива (мультипликатор $|dx_{j+2}/dx_j|<1$); для движения с периодом 1

(в единицах $T_0$) интервал времени 2 --- тоже период. При $R=R_1$

мультипликатор обращается в $+1$ и при $R>R_1$ точка $x_*=0$ становится

неустойчивой. Вэтот момент рождается пара устойчивых неподвижных точек

$$

x_*^{(1),(2)} = \pm \left [{R-R_1\over 1+\beta}\right ]^{1/2},\eqno(4)

$$

которые и соответствуют устойчивому предельному циклу удвоенного периода

\footnote {$^{7)}$}{или, как мы будем говорить для краткости,

{\it 2-циклу.\/} Относящиеся к нему неподвижные точки будем называть

{\it элементами цикла \/}. }; преобразование (3) оставляет каждую из этих точек

на месте, а преобразование (2) переводит каждую из них в другую. Подчеркнём,

что цикл единичного периода при описанной бифуркации не исчезает --- он остаётся

решением уравнений движения, но неустойчивым.

Вблизи бифуркации движение остаётся ещё ``почти периодическим'' с периодом 1:

точки последовательных возвратов траектории $x^{(1)}$ и $x^{(2)}$ близки друг к

другу. Интервал $x^{(1)}-x^{(2)}$ между ними является мерой амплитуды колебаний

с периодом 2; она растёт с надкритичностью как $(R-R_1)^{1/2}$.

Многократное повторение бифуркаций удвоения периода открывает один из

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

бесконечно, причём они следуют друг за другом (по мере увеличения $R$) через

все убывающие интервалы; последовательность критических значений

$R_1, R_2,\ldots$ стремится к конечному пределу, за которым периодичность исчезает

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

в этом сценарии с возникновением турбулентности. Мы увидим, что этот сценарий

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

({\it M. J. Feigenbaum, 1978\/})%\footnote {$^{1)}$}{Последовательность бифуркаций

%удвоения периода (нумеруемых далее порядковыми номерами $1, 2, \ldots $)

%не обязательно должна начинаться с первой же бифуркации периодического движения.

%Она может, в принципе, начаться и после нескольких первых бифуркаций с

%возникновением несоизмеримых частот, после их синхронизации за счёт рассмотренного

%ранее механизма.}

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

следуют друг за другом (при увеличении $R$) настолько быстро, что даже в

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

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

быть описана одномерным отображением Пуанкаре, зависящим от одного параметра.

Выбор рассматриваемого ниже отображения естественен в силу следующих

соображений. В значительной части интервала изменения переменной $x$

отображение должно быть ``растягивающим'', $|df(x;\lambda)/dx|>1;$ это даёт

возможность возникновения неустойчивостей. Отображение должно также

возвращать траектории, выходящие за границы некоторого интервала, обратно в

него; противное означало бы неограниченное возрастание амплитуд пульсаций

скорости, что невозможно. Обоим этим требованиям вместе могут удовлетворять

лишь немонотонные функции $f(x;\lambda),$ т.~е. не взамнооднозначные отображения

(1):значение $x_{j+1}$ однозначно определяется предшествующим значением

$x_j,$ но не наоборот. Простейший вид такой функции --- функция с одним

максимумом; в окрестности максимума положим

$$

x_{j+1} = f(x_j;\lambda) = 1-\lambda x_j^2,\eqno(5)

$$

где $\lambda $ --- положительный параметр, который надо рассматривать

(в гидродинамическом аспекте) как возрастающую функцию $R$\footnote {$^{8)}$}

{Подчеркнём, что допустимость не взаимнооднозначных отображений связана с

приближённостью одномерного расмотрения. Если бы все траектории располагались

строго на одной поверхности $\Sigma $ (так что отображение Пуанкаре было бы

строго одномерным), подобная неоднозначность была бы невозможна: она означала

бы пересечение траекторий (две траектории с различными $x_j$ пересекались бы в

точке $x_{j+1}$). В этом же смысле следствием приближённости является

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

отображения расположена в экстремуме отображающей функции (такая точка

может быть названа ``сверхустойчивой'' --- приближение к ней происходит по

закону более быстрому, чем указанный выше).}. Примем условно отрезок

$[-1,+1]$ как интервал изменения величины $x;$ при $\lambda $ между 0~и 2

все итерации отображения (5) оставляют $x$ в том же интервале.

Преобразование (5) имеет неподвижную точку --- корень уравнения

$x_* = 1-\lambda x_*^2$. Эта точка становится неустойчивой при

$\lambda >\Lambda_1,$ где $\Lambda_1$ --- значение параметра $\lambda ,$

для которого мультипликатор $\mu = -2\lambda x_* = -1;$ из двух написанных

уравнений находим $\Lambda_1 =3/4.$ Это --- первое критическое значение

параметра $\lambda ,$ определяющее момент первой бифуркации удвоения периода:

появления 2-цикла. Проследим за появлением последующих бифуркаций с помощью

приближённого приёма, позволяющего выяснить некоторые качественные особенности

процесса, хотя и не дающего точных значений характерных констант; затем будут

сформулированы точные утверждения.

Повторив преобразования (5) дважды, получим

$$

x_{j+2} = 1-\lambda +2\lambda^2 x_j^2-\lambda^3 x_j^4,\eqno(6)

$$

Пренебрежём здесь последним слагаемым --- четвёртой степени по $x_j.$

Оставшееся равенство масштабным преобразованием\footnote {$^{9)}$}{ Это

преобразование невозможно при значении $\lambda =1 $ (при котором неподвижная точка

отображения (6) совпадает с центральным экстремумом: $x_*=0$). Это значение,

однако, заведомо не является интересующим нас следующим критическим значением

$\Lambda_2.$}

$$

x_j\to x_j/\alpha_0, \quad \alpha_0 = 1/(1-\lambda)

$$

приводится к виду

$$

x_{j+2}= 1-\lambda_1 x_j^2,

$$

отличающемуся от (5) лишь заменой параметра $\lambda$ на

$$

\lambda_1 = \phi (\lambda)\equiv 2\lambda^2(\lambda-1).\eqno(7)

$$

Повторяя эту операцию с масштабнымми множителями

$\alpha_1= 1/(1-\lambda_1),\ldots ,$ получим ряд последовательных отображений

того же вида:

$$

x_{j+2^m}= 1-\lambda_m x_j^2,\quad \lambda_m =\phi (\lambda_{m-1}).\eqno(8)

$$

Неподвижные точки отображений (8) отвечают $2^m$-циклам\footnote{$^{10)}$}

{Во избежание недоразумений подчеркнём, что после произведённых масштабных

преобразований отображения (8) должны быть определены теперь на растянутых

интерваах $|x|\le |\alpha_1 \alpha_2\ldots\alpha_{m-1}|$ (а не на $|x|\le 1,$

как в (5),(6)). Однако в силу сделанных пренебрежений выражения (8) могут

фактически описывать лишь область вблизи центральных экстремумов отображающих

функций.}. Поскольку все эти отображения имеют тот же вид, что и (5), то

можно сразу заключить, что $2^m$-циклы ($m=1,2,3,\ldots$)становятся

неустойчивыми при $\lambda_m =\Lambda_1 =3/4.$ Соответствующие же критические

значения $\Lambda_m$ исходного параметра $\lambda$ получаются путём решения

цепочки уравнений

$$

\Lambda_1 =\phi (\Lambda_2),\quad\Lambda_2 =\phi (\Lambda_3),\ldots ,

\Lambda_{m-1} =\phi (\Lambda_m).

$$

%графически они даются построением, показанным на рис.~7.

Очевидно, что при

$m\longrightarrow \infty $ последовательность этих чисел сходится к конечному

пределу $\Lambda_\infty$ --- корню уравнения $\Lambda_\infty =\phi (\Lambda_\infty);$

он равен $\Lambda_\infty =(1+\sqrt 3)/2\approx 1{,}37.$ К конечному пределу стремятся и

масштабные множители: $\alpha_m\to\alpha,$

где $\alpha =1/(1-\Lambda_\infty)\approx -2{,}8.$

Легко найти закон, по которому происходит приближение $\Lambda_m$ к

$\Lambda_\infty$ при больших $m.$ Из уравнения $\Lambda_m =\phi (\Lambda_{m+1}$

при малых разностях $\Lambda_\infty -\Lambda_m $ находим

$$

\Lambda_\infty -\Lambda_{m+1} ={1\over\delta}(\Lambda_\infty -\Lambda_m),\eqno(9)

$$

где $\delta =\phi '(\Lambda_\infty) =4+\sqrt 3 \approx 5{,}73.$ Другими словами,

$\Lambda_\infty -\Lambda_m\sim\delta^m,$ т.~е. значения $\Lambda_m$ приближаются

к пределу по закону геометричекой прогрессии. По такому же закону меняются

интервалы между последовательными критическими числами: (9) можно переписать в

эквивалентном виде

$$

\Lambda_{m+2} -\Lambda_{m+1}={1\over\delta}(\Lambda_{m+1}-\Lambda_m).\eqno(10)

$$

В гидродинамическом аспекте, как уже указывалось, параметр $\lambda $ надо

рассматривать как функцию числа Рейнольдса, соответственно чемупоявляются

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

удвоения периода и стремящиеся к конечному пределу $R_\infty .$ Очевидно,

что для этих значений справедливы те же предельные законы (9),(10) (с той

же постоянной $\delta $), что и для чисел $\Lambda_m .$

Изложенные рассуждения иллюстрируют происхождение основных закономерностей

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

сходятся к пределу $\Lambda_\infty $ по закону (9),(10); появление

масштабного множителя $\alpha .$ Полученные при этом значения харакерных

констант, однако, не точны. Точные значения (полученные путём многократного

компьютерного итерирования отображения (5)) показателя сходимости $\delta$

({\it число Фейгенбаума\/}) и масштабного множителя $\alpha$:

$$

\delta =4{,}6692\ldots ,\quad\alpha =-2{,}5029\ldots\eqno(11)

$$

а предельное значение $\Lambda_\infty =1{,}401$\footnote{$^{11)}$}{ Значение

$\Lambda_\infty $ имеет несколько условный характер, поскольку оно зависит от

способа введения параметра в исходное отображение --- функцию $f(x;\lambda)$

(значения же $\delta$ и $\alpha$ от этого не зависят вовсе).}. Обратим внимание

на сравнительно большое значение $\delta ;$ быстрая сходимость приводит к тому,

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

периода.

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

степенями $x_j^2,$ кроме первой, отображение (8) позволяет установить лишь факт

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

элементы описываемого этим отображением $2^m$-цикла\footnote{$^{12)}$}{ То есть

все $2^m$ точки

$x_*^{(1)},x_*^{(2)},\ldots ,$

переходящие последовательно

друг в друга (периодические) при итерациях отображения (31.5) и неподвижные

(и устойчивые) по отношению к $2^m$-кратно итерированному отображению.

Отметим, что производные $dx_{j+2^m}/dx_j$ во всех точках

$x_*^{(1)},x_*^{(2)},\ldots $ автоматически одинаковы (и потому одновременно

проходят через $-1$ в момент следующей бифуркации).}. В действительности

итерированные отображения (5) представляют собой полиномы по $x_j^2,$

степень которых при каждой итерации возрастает вдвое. Они представляют собой

сложные функции от $x_j$ с быстро возрастающим числом экстремумов, симметрично

расположенных по отношению к точке $x_j=0$ (которая тоже всегда остаётся

экстремумом).

Замечательно, что нетолько значения $\delta$ и $\alpha$, но и предельный вид

самого бесконечно кратно итерированного отображения оказываются в определённом

смысле независящими от вида начального отображения $x_{j+1}=f(x_j;\lambda) :$

достаточно, чтобы зависящая от одного параметра функция $f(x;\lambda)$ была

гладкой функцией с одним квадратичным максимумом (пусть это будет в точке $x=0$);

она не обязана даже быть симметричной относительно этой точки вдали от неё. Это

свойство {\it универсальности\/} существенно увеличивает степень общности

излагаемой теории. Его точная формулировка состоит в следующем.

Рассмотрим отображение, задаваемое функцией $f(x)$ (функция $f(x;\lambda)$ с

определённым выбором $\lambda$ -- см. ниже ),нормированной условтем $f(0)=1.$

Применив его дважды, получим функцию $f(f(x)).$ Изменим масштаб как самой этой

функции, так и переменной $x$ в $\alpha_0 =1/f(1)$ раз; таким образом получим

новую функцию

$$

f_1(x)=\alpha_0f(f(x/\alpha_0)),

$$

для которой снова будет $f_1(0)=1.$ повторяя эту операцию, получим

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

$$

f_{m+1}(x)=\alpha_m f_m(f_m(x/\alpha_m))\equiv\widehat Tf_m,\quad \alpha_m =1/f_m(1).\eqno(12)

$$

Если эта последовательность стремится при $m\longrightarrow\infty$ к некоторой

определённой предельной функции $f_\infty(x)\equiv g(x),$ эта последняя

должна быть ``неподвижной фунуцией'' определённого в (12) оператора $\widehat T,$

т.~е. должна удовлетворять функциональному уравнению

$$

g(x)=\widehat TG\equiv\alpha g(g(x/\alpha)),\quad\alpha =1/g(1),\quad g(0)=1.\eqno(13)

$$

Всилу предположенных свойств допустимых функций $f(x),$ функция $g(x)$ должна

быть гладкой и и иметь квадратичный экстремум в точке $x=0;$ никакого другого

следа от конкретного вида $f(x)$ в уравнении (13) или в налагаемых на его

решение условиях не остаётся. Подчеркнём, что после произведённых при выводе масштабных

преобразований (с $|\alpha_m|>1$) решение уравнения определяется при всех

значениях фигурирующей в нём переменной $x$ от $-\infty$ до $+\infty$ (а не

только на интервале $-1\le x\le +1$). Функция $g(x)$ автоматически является

чётной по $x;$ она должна быть такой, поскольку среди допустимых функций

$f(x)$ имеются чётные, а чётное отображение заведомо остаётся чётным после

любого числа итераций.

Такое решение уравнения (13) действительно существует и единственно (хотя и

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

с бесконечным числом экстремумов, неограниченную по своей величине; постоянная

$\alpha$ определяется вместе с самой функцией $g(x).$ Фактически достаточно

построить эту функцию на интервале $[-1,1]$, после чего она может быть

продолжена за его пределы итерированием операции $\widehat T.$ Обратим внимание

на то, что на каждом шаге итерирования $\widehat T$ в (12) значения функции

$f_{m+1}(x)$ на интервале $[-1,1]$ определяются значениями функции $f_m(x)$

на сокращённой в $|\alpha_m|\approx |\alpha|$ раз части этого отрезка. Это

значит, что в пределе многократных итераций для определения функции $g(x)$

на интервале $[-1,1]$ (а тем самым и на всей оси $x$) существенны всё меньшие и

меньшие части исходной функции вблизи её максимума; в этом и состоит, в

конечном итоге, источник универсальности\footnote{$^{13)}$}{Уверенность в

существовании единственного решения уравнения (13) основана на компьютерном

моделировании. Решение ищется (на интервале $[-1,1]$) в виде полинома высокой

степени по $x^2;$ точность моделирования должна быть тем выше, чем до более

широкой области значений $x$ (вне указаного отрезка) мы хотели бы за тем

продолжить функцию итерированием $\widehat T.$ На интервале $[-1,1] $ функция

$g(x)$ имеет один экстремум, вблизи которого $g(x)=1-1,528x^2$ (если считать

экстремум максимумом; этот выбор условен ввиду инвариантности уравнения (13)

отонсительно изменения знака $g$).}.

Функция $g(x)$ определяет структуру апериодического аттрактора, возникающего

в результате бесконечной последовательности удвоений периода. Но это

происходит при вполне опредеоённом для функции $f(x;\lambda)$ значении параметра

$\lambda =\Lambda_\infty .$ Ясно поэтому, что функции, образованные из

$f(x;\lambda)$ путём многократного итерирования преобразования (12),

действительно сходятся к $g(x)$ лишь при этом изолированном значении $\lambda .$

Отсюда в свою очередь следует, что неподвижная функция оператора $\widehat T$

неустойчива по отношению к её малым изменениям, отвечающим малым отклонениям

параметра $\lambda$ от значения $\Lambda_\infty .$ Исследование этой

неустойчивости даёт возможность определения универсальной постоянной $\delta$ ---

снова без всякой связи с конкретным видом функции $f(x).$

Масштабный множитель $\alpha$ определяет изменение --- уменьшение ---

геометрических (в пространстве состояний)характеристик аттрактора на каждом

шаге удвоений периода; этими характеристиками являются расстояния между

элементами предельных циклов на оси $x.$ Поскольку, однако, каждое удвоение

сопровождается ещё и увеличением числа элементов цикла, это утверждение должно

быть конкретизировано и уточнено. При этом заранее ясно, что закон изменения

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

\footnote{$^{14)}$}{Имеются ввиду расстояния на нерастянутом отрезке $[-1,1],$

условно выбранном с самого начала как интервал изменения $x,$ на котором

расположены все элементы циклов. Отрицательность $\alpha$ означает, что при

бифуркациях происходит также инверсия расположения элементов относительно точки

$x=0.$}. Действительно, если две близкие точки преобразуются через почти

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

$|\alpha|$ раз; если же преобразование происходит через участок функции вблизи

её экстремума --- расстояние сократится в $\alpha^2$ раз.

В момент бифуркации (при $\lambda =\Lambda_m$) каждый элемент (точка) $2^m$-цикла

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

постепенно возрастает, но точки остаются ближайшими друг к другу на всём

протяжении изменения $\lambda$ до следующей бифуркации. Если следить за переходами

элементов цикла друг в друга с течением времени (т.~е. при последовательных

отображениях $x_{j+1}=f(x_j;\lambda)$, то каждая из компонент пары перейдёт в

другую через $2^m$ единиц времени. Это значит, что расстояние между точками

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

этом смысле представляет особый физический интерес.

Расположим все элементы $2^{m+1}$-цикла в том порядке, в котором они обходятся

со временем, и обозначим их как $x_{m+1}(t)$ где время $t$ (измеренное в единицах

основного периода $T_0$) пробегает целочисленные значения $t/T_0 =1,2,\ldots ,2^{m+1}.$

Эти элементы возникают из элементов $2^m$-цикла расщеплением последних на пары.

Интервалы между точками каждой пары даются разностями

$$

\xi_{m+1}(t)=x_{m+1}(t)-x_{m+1}(t+T_m),\eqno(14)

$$

где $T_m=2^mT_0=T_{m+1}/2$ --- период $2^m$-цикла, т.~е. половина периода

$2^{m+1}$-цикла. Введём функцию $\sigma_m (t)$ --- масштабный множитель,

определяющий изменение интервалов (14) при переходе от одного цикла к

следующему\footnote {$^{15)}$}{ Поскольку оба цикла существуют в разных

интервалах значения параметра $\lambda$ (на интервалах $(\Lambda_{m-1},\Lambda_m)$

и $(\Lambda_m,\Lambda_{m+1})),$ и на этих интервалах величины (14) существенно

меняются, то их смысл в определении (15) нуждается в уточнении. Будем понимать

их при тех значениях параметра $\lambda ,$ когда циклы ``сверхустойчивы'';

по одному такому значению имеется в области существования каждого цикла.}:

$$

{\xi_{m+1}(t)\over\xi_m(t)}=\sigma_m(t).\eqno(15)

$$

Очевидно, что

$$

\xi_{m+1}(t+T_m)=-\xi_{m+1}(t),\eqno(16)

$$

и поэтому

$$

\sigma_m(t+T_m)=-\sigma_m(t).\eqno(17)

$$

Функция $\sigma_m(t)$ имеет сложные свойства, но можно показать, что её

предельный (при больших $m$) вид с хорошей точностью аппроксимируется простым

образом:

$$

\sigma_m(t)=\left \{

\eqalign {&1/\alpha  \cr

          &1/\alpha^2\cr}\right.

\eqalign {

&\hbox{\rm\ при\ } \quad 0<t<T_m /2,\cr

&\hbox{\rm\ при\ } \quad T_m /2<t<T_m\cr}\eqno(18)

$$

(при надлежащем выборе начала отсчёта $t$)\footnote{$^{16)}$}{Здесь не будет

проводится в принципе простое, но громоздкое исследование свойств функции

$\sigma_m(t)$. См. {\it Фейгенбаум~М.\/} --- УФН, 1983, т.141, с.343

[Los Alamos Science, 1980, v. 1, p. 4].}.

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

(частотного) движения жидкости, претерпевающей удвоения периода. В

гидродинамическом аспекте величину $x_m(t)$ надо понимать как характеристику

скорости жидкости. Для движения с периодом $T_m$ спектр функции $x_m(t)$

(от непрерывного времени $t$!) содержит частоты $k\omega_m$ ($k=1,2,3,\ldots$) ---

основную частоту $\omega_m =2\pi/T_m$ и её гармоники. После удвоения периода

течение описывается функцией $x_{m+1}(t)$ с периодом $T_{m+1}=2T_m .$ Её

спектральное разложение содержит, наряду с теми же частотами $k\omega_m ,$

ещё и субгармоники частоты $\omega_m$ --- част\'оты $l\omega_m /2, l=1,3,5,\ldots$

Представим $x_{m+1}(t)$ в виде

$$

x_{m+1}(t)={1\over 2}\{ \xi_{m+1}(t)+\eta_{m+1}(t)\},

$$

где $\xi_{m+1}$ --- разность (14), а

$$

\eta_{m+1}(t)=x_{m+1}(t) + x_{m+1}(t+T_m).

$$

Спектральное разложение $\eta_{m+1}(t)$ содержит только частоты $k\omega_m$;

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

$$

{1\over T_{m+1}}\int\limits_0^{T_{m+1}} \eta_{m+1}(t)\, e^{i\pi lt/T} dt=

{1\over 2T_m}\int\limits_0^{T_m} [\eta_{m+1}(t)-\eta_{m+1}(t+T_m)] e^{i\pi lt/T} dt

$$

обращаются в нуль в силу равенства $\eta_{m+1}(t+T_m)=\eta_{m+1}(t)$.

С другой стороны, величины $\eta_m(t)$ в первом приближении не меняются при

бифуркации: $\eta_{m+1}(t)\approx \eta_m(t)$; это значит, что интенсивность

колебаний с частотами $k\omega_m$ тоже остаётся неизменной.

Спектральное же разложение величин $\xi_{m+1}(t)$ содержит, напротив, только

субгармоники $l\omega_m/2$ --- новые частоты, появляющиеся на $(m+1)$-м шаге

удвоений. Суммарная интенсивность этих спектральных компонент определяется

интегралом

$$

I_{m+1}={1\over T_{m+1}}\int\limits_0^{T_{m+1}} \xi_{m+1}^2(t) dt.\eqno(19)

$$

Выразив $\xi_{m+1}(t)$ через $\xi_m(t)$, пишем

$$

I_{m+1}={1\over 2T_m}2\int\limits_0^{T_m} \sigma_m^2(t) \xi_m^2(t) dt.

$$

С учётом (16--18) получим

$$

I_{m+1}={1\over 2}\left( {1\over \alpha^2}+{1\over \alpha^4}\right){1\over T_m}

\int\limits_0^{T_m} \xi_m^2(t) dt =

{1\over 2}\left({1\over \alpha^2}+{1\over \alpha^4}\right) I_m,

$$

и окончательно

$$

I_m/I_{m+1}\approx 10{,}8.\eqno(20)

$$

Таким образом, интенсивность новых спектральных компонент, появляющихся после

бифуркации удвоения периода, превышает таковую для следующей бифуркации в

определённое, не зависящее от номера бифуркации, число раз.

({\it M.~J. Feigenbaum, 1979\/})

\footnote{$^{17)}$}{Это относится не только к суммарной интенсивности появляющихся

субгармоник, но и к интенсивности каждой из них. На каждую субгармонику,

появляющуюся после $m$-ой бифуркации, приходится по две (по одной справа и слева)

субгармоники  после $(m+1)$-ой бифуркации. Поэтому отношение интенсивностей

отдельных новых появляющихся после двух последовательных бифуркаций спектральных

пиков вдвое больше величины (20). Более точное значение этой величины 10{,}48.

Оно получается путём анализа состояния в самой точке $\lambda=\Lambda_\infty$

с помощью универсальной функции $g(x)$; в этой точке присутствуют уже все

частоты и вопрос указанный в сноске 15  не возникает. См. {\it Nanenberg~M.,

Rudnick~J.\/} --- Phys. Rev., 1981, v. 24B, p. 493.}.

\beginsection{Переход к турбулентности через перемежаемость}

Рассмотрим, наконец, разрушение периодического движения при прохождении

мультипликатора через значение $\mu=+1$.

Этот тип бифуркации описывается (в рамках одномерного отображения Пуанкаре)

функцией $x_{j+1}=f(x_j;R)$, которая при определённом значении параметра

(числа Рейнольдса), $R=R_0$, касается прямой $x_{j+1}=x_j$. Выбрав точку

касания в качестве $x_j=0$, напишем вблизи нее разложение функции отображения

в виде

\footnote{$^{18)}$}{Коэффициент при $R-R_0$ и коэффициент (положительный)

при $x_j^2$ можно обратить в единицу соответствующим определением $R$~и $x_j$,

что и предполагается в (22).}.

$$

x_{j+1}=(R-R_0)+x_j+x_j^2.\eqno(21)

$$

При $R<R_0$ %(рис.~23)

существуют две неподвижные точки

$$

x_*^{(1),(2)}=\mp (R_0-R)^{1/2},

$$

из которых одна ($x_*^{(1)}$) отвечает устойчивому, а другая ($x_*^{(2)}$) ---

неустойчивому периодическому движению. При $R=R_0$ мультипликатор в обоих

точках становится равным $+1$, оба периодических движения сливаются и при

$R>R_0$ исчезают (неподвижные точки переходят в комплексную область).

При малой надкритичности расстояние между линией (21) и $x_{j+1}=x_j$ мало

(в области $x_j=0$). На этом интервале значений $x$, следовательно, каждая

итерация отображения (21) лишь незначительно перемещает след траектории, и

для прохождения им всего интервала потребуется много шагов. Другими словами,

на сравнительно большом промежутке времени траектория в пространстве состояний

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

в физическом пространстве регулярное (ламинарное) движение жидкости. Отсюда

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

({\it P.~Manneville, Y.~Pomeau, 1980\/}).

Можно представить себе, что к рассмотренному участку функции отображения

примыкают участки, приводящие к хаотизации траекторий; им отвечает в пространстве

состояний множество локально неустойчивых траекторий. Это множество, однако,

само по себе не является аттрактором и с течением времени точка, изображающая

систему, его покидает. При $R<R_0$  траектория выходит на устойчивый цикл, т.~е.

в физическом пространстве устанавливается  ламинарное периодическое движение.

При $R>R_0$ устойчивый цикл отсутствует и возникает движение, в котором

``турбулентные'' периоды чередуются с ламинарными (отсюда название сценария ---

переход через {\it перемежаемость\/}).

О длительности турбулентных периодов нельзя сделать каких-либо общих заключений.

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

Для этого напишем разностное уравнение (21) в виде дифференциального. Имея в

виду малость изменения $x_j$ на одном шаге отображения, заменим разность

$x_{j+1}-x_j$ производной $dx/dt$ по непрерывной переменной $t$:

$$

dx/dt=(R-R_0)+x^2.

$$

Найдём время $\tau$, необходимое для прохождения отрезка между точками $x_1$~и

$x_2$, лежащими по обе стороны точки $x=0$ на расстояниях, больших по сравнению

с $(R-R_0)^{1/2}$, но ещё в области применимости разложения (21). Имеем

$$

\tau=(R-R_0)^{-1/2}\left.{\rm arctg}[x(R-R_0)^{-1/2}]\right|_{x_1}^{x_2},

$$

откуда

$$

\tau \sim (R-R_0)^{-1/2},

$$

чем и определяется искомая зависимомть;  длительность ламинарных периодов

убывает с ростом надкритичности.

В этом сценарии остаётся открытым как вопрос о пути подхода к его началу, так

и вопрос о природе возникающей турбулентности.

\beginsection{Развитая турбулентность}

\rightline{ Big whirls have little whirls that feed on their velocity,}

\rightline{ and little whirls have lesser whirls and so on to viscosity}

\rightline{ --- in the molecular sense.}

\rightline{\it Richardson L.}

\medskip\noindent

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

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

временем в каждой точке потока ({\it развитая турбулентность\/}); скорость всё

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

нерегулярное изменение скорости имеет место от точки к точке потока,

рассматриваемого в заданный момент времени. В настоящее время полной

количественной теории развитой турбулентности ещё не существует. Известен,

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

Введём понятие о средней скорости движения, получающейся в результате усреднения

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

При таком усреднении нерегулярность изменения скорости сглаживается и средняя

скорость оказывается плавно меняющейся вдоль потока функцией. Будем в

дальнейшем обозначать среднюю скорость посредством~$u.$ Разность $v'=v-u$ между

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

нерегулярное изменение, будем называть {\it пульсационной\/} частью скорости.

Рассмотрим подробнее характер накладывающегося на усреднённый поток нерегулярного,

пульсационного движения. Это движение можно в свою очередь качественно

рассматривать как результат наложения движений ({\it турбулентных пульсаций\/})

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

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

скорость движения). По мере возрастания числа Рейнольдса появляются сначала

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

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

присутствуют пульсации с масштабами от самых больших до очень малых. Основную

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

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

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

величины этого {\it основного\/} (или {\it внешнего\/}) масштаба турбулентного

движения посредством~$l.$ Эти крупномасштабные движения обладают наибольшими

амплитудами. Их скорость по порядку величины сравнима с изменениями $\Delta u$

средней скорости на протяжении расстояний $l$ (имеется ввиду порядок величины не

самой скорости, а её изменения, поскольку именно оно характеризует скорость

турбулентного движения; абсолютная же величина средней скорости может быть

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

движение)\footnote {$^{18)}$}{В действительности, по видимому, масштабы основных

пульсаций в несколько раз меньше, чем характерные размеры $l$, а их скорость ---

в несколько раз меньше, чем $\Delta u.$}. Что же касается частот этих пульсаций,

то они --- порядка отношения $u/l$ средней скорости $u$ (а не её изменения

$\Delta u $) к размерам $l.$ Действительно, частота определяет период

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

отсчёта. Но относительно такой системы вся эта картина движентся вместе со всей

жидкостью со скоростью порядка $u.$

Мелкомасштабные же пульсации, соответствующие большим частотам, участвуют в

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

как мелкую детальную структуру, накладывающуюся на основные крупномасштабные

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

малая часть всей кинетической энергии жидкости.

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

изменения пульсационной скорости вдоль потока (рассматриваемого в заданный момент

времени). На протяжении больших расстояний (сравнимых с $l$) изменение

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

пульсаций и потому сравнимо по величине с $\Delta u.$ На малых же (по сравнению

с $l$) расстояниях оно определяется мелкомасштабными пульсациями и потому мало

по сравнению с $\Delta u $ (но велико по сравнению с изменением средней скорости

на том же малом расстоянии). Такая же картина имеет место, если наблюдать

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

(по сравнению с характеристическим временем~$T\sim l/u$) интервалов времени

скорость испытывает незначительные изменения; в течении же больших промежутков

времени скорость мняется на величины~$\sim\Delta u.$

В число Рейнольдса $R$, определяющее свойства течения жидкости в целом, в

качестве характеристических размеров входит длина $l.$ Наряду с таким числом

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

различных масштабов. Если $\lambda$ --- масштаб пульсаций, а $v_\lambda$ ---

порядок величины их скорости, то $R_\lambda\sim v_\lambda\lambda /\nu .$ Это

число тем меньше, чем меньше масштаб движения.

При больших $R$ велики также и числа Рейнольдса $R_\lambda$ крупномасштабных

пульсаций. Но большие числа Рейнольдса эквивалентны малым вязкостям. Отсюда

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

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

крупномасштабных пульсациях не происходит и заметной диссипации энергии.

Вязкость жидкости становится существенной только для самых мелкомасштабных

пульсаций, для которых $R_\lambda\sim 1$ (масштаб $\lambda_0$ этих пульсаций будет

определён ниже). Именно в этих мелкомасштабных пульсациях, не существенных с точки

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

диссипация энергии.

Мы приходим, таким образом, к следующему представлению о диссипации энергии при

турбулентном движении ({\it L. Richardson, 1922\/}). От пульсаций с б\'ольшими

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

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

энергии от крупно- к мелкомасштабным пульсациям, т.~е. от малых частот к большим.

Этот поток диссипируется, т.~е. кинетическая энергия переходит в тепло, в самых

мелкомасштабных пульсациях. (Именно описание этого процесса, являющееся

перефразированием Свифта, вынесено эпиграфом раздела.)

Разумеется, для поддержания ``стационарного'' состояния

потока необходимо наличие внешних источников энергии, непрерывно передающих её

основному крупномасштабному движению.

Поскольку вязкость жидкости существенна только для самых мелкомасштабных

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

$\lambda\gg\lambda_0,$ не могут зависеть от $\nu$ (более точно, эти величины

не должны меняться при изменении $\nu$ и неизменных остальных условиях, в

которых происходит движение). Это обстоятельство сужает круг величин,

определяющих свойства турбулентного движения, в результате чего для исследования

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

размерностью имеющихся в нашем распоряжении величин.

Применим такие соображения к определению порядка величины диссипации энергии при

турбулентном движении. Пусть $\varepsilon$ есть среднее количество энергии,

диссипируемой в единицу времени в единице массы жидкости\footnote {${19)}$}{Здесь

и далее $\varepsilon$ будет обозначать среднюю диссипацию энергии, а не внутреннюю

энергию жидкости!}. Мы видели, что эта энергия черпается из крупномасштабного

движения, откуда постепенно передаётся во всё меньшие масштабы, пока не

диссипируется в пульсациях масштаба~$\sim\lambda_0.$ Поэтому, несмотря на то,

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

$\varepsilon$ может быть определён с помощью одних только величин, характерных

для крупномасштабных движений. Таковыми являются плотность жидкости $\rho,$

размеры $l$ и скорость $\Delta u.$ Из этих трёх величин можно составить всего

одну комбинацию, обладающую той же размерностью, что и $\varepsilon$, т.~е.

эрг/г$\cdot$c$=$см$^2$/с$^3.$ Таким способом получаем:

$$

\varepsilon\sim {(\Delta u)^3\over l} ,\eqno(1)

$$

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

Турбулентно движущуюся жидкость можно в некоторых отношениях качественно описывать

как жидкость, обладающую некоторой, как говорят, {\it турбулентной вязкостью\/}

$\nu_{turb}$, отличной от истинной кинематической вязкости $\nu .$ Характеризуя

свойства турбулентного движения, $\nu_{turb}$ должно по поряжку величины

определяться величинами $\rho ,\, \Delta u ,\, l .$ Единственной составленной из

них величиной с размерностью кинематической вязкости является $\Delta u\cdot l,$

поэтому

$$

\nu_{turb}\sim\Delta u\cdot l .\eqno(2)

$$

Отношение турбулентной вязкости к обычной

$$

{\nu_{turb}/\nu}\sim R,\eqno(3)

$$

т.~е. растёт с числом Рейнольдса\footnote {$^{20)}$}{В действительности в этом отношении

должен стоять ещё довольно значительный численный коэффициент. Это связано с указанным выше

обстоятельством, что $l$ и $\Delta u$ могут довольно заметно отличаться от истинных

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

$$

{\nu_{turb}/\nu}\sim R/R_0,

$$

где учитывается, что $\nu_{turb}$ и $\nu$ должны в действительности

сравниваться не при $R\sim 1$, а при $R\sim R_0 .$}.

Диссипация энергии выражается через $\nu_{turb}$ формулой

$$

\varepsilon\sim\nu_{turb}(\Delta u /l)^2 ,\eqno(4)

$$

в соответствии с обычным определением вязкости. В то время как $\nu$ определяет

диссипацию энергии по производным от истинной скорости по координатам,

турбулентная вязкость связывает диссипацию с градиентом ($\sim\Delta u /l$) средней

скорости движения.

Наконец, укажем, что порядок величины $\Delta p$ изменения давления на протяжении

области турбулентного движения тоже может быть определён из соображений подобия:

$$

\Delta p\sim\rho(\Delta u)^2 .\eqno(5)

$$

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

можно составить из $\rho ,\ l , $~и $\Delta u.$

Перейдём теперь к изучению свойств развитой турбулентности в масштабах $\lambda,$

малых по сравнению с основным масштабом $l.$ Об этих свойствах говорят как о

локальных свойствах турбулентности. При этом мы будем рассматривать жидкость

вдали от твёрдых стенок,--- точнее, на расстояниях от них, больших по сравнению

с $\lambda.$

О такой мелкомасштабной турбулентности вдали от твёрдых тел можно высказать

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

изотропии. Последнее означает, что в участках, размеры которых малы по сравнению

с $l,$ свойства турбулентного движения одинаковы по всем направлениям; в частности,

они не зависят от направления скорости усреднённого движения. Подчеркнём, что

здесь и далее везде, где говорится о свойствах турбулентного движения в малом

участке жидкости, подразумевается относительное движение жидких частиц в этом

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

целом и которое связано с движением более крупных масштабов.

Оказывается возможным получить ряд существенных результатов о локальных свойствах

турбулентности непосредственно из соображений подобия ({\it А.~Н. Колмогоров, 1941;

А.~М. Обухов, 1941\/}).

Для этого займёмся выяснением, какими параметрами могут вообще определяться

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

по сравнению с расстояниями $\lambda_0,$ на которыхначинает играть роль вязкость

жидкости; ниже будет идти речь именно о таких расстояниях. Этими параметрами

является плотность $\rho$ жидкости и, кроме того, ещё одна характерная для

турбулентного потока величина --- энергия $\varepsilon$, диссипируемая в единицу

времени в единице массы жидкости. Мы видели, что $\varepsilon$ представляет поток

энергии непрерывно передаваемой от пульсаций с б\'ольшими к пульсациям с меньшими

масштабами. Поэтому, хотя диссипация энергии обусловливается в конечном итоге

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

величина $\varepsilon$  определяет свойства движения и в б\'ольших масштабах. Что

касается масштабов $l$ и $\Delta u$ размеров и скорости движения в целом, то

естесственно считать, что (при заданных $\rho$ и $\varepsilon$) локальные свойства

турбулентности от этих величин не зависят. Вязкость жидкости $\nu$ тоже не может

входить ни в какие интересующие нас теперь величины (речь идёт о расстояниях

$\lambda\gg\lambda_0$).

Определим порядок величины $v_\lambda$ изменения скорости турбулентного движения

на протяжении расстояний порядка $\lambda .$ Оно должно определяться только

величиной $\varepsilon$ и, разумеется, самим расстоянием \footnote {$^{21)}$}{

Величина $\varepsilon$ имеет размерность эрг/(г$\cdot$с)$=$см$^2/$с$^3,$ не

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

массы, является плотность $\rho.$ Поэтому последняя вообще не участвует в

составлении величин, размерность которых не содержит размерности массы.}$\lambda .$

Из этих двух величин можно составить всего одну комбинацию с размерностью

скорости:

$$

(\varepsilon\lambda)^{1/3}.\eqno(6)

$$

Таким образом, изменение скорости на протяжении малого расстояния

пропорционально кубическому корню из этого расстояния

({\it закон Колмогорова-Обухова\/}). Величину $v_\lambda$ можно рассматривать

и как скорость турбулентных

движений масштаба~$\lambda$: изменение средней скорости на малых расстояниях

мало по сравнению с изменением пульсационной скорости на этих же расстояниях, и

им можно пренебречь.

К соотношению (6) можно прийти и другим путём, выражая постоянную величину ---

диссипацию $\varepsilon$ --- через величины, характеризующие пульсации

масштаба~$\lambda .$ При этом $\epsilon$ должно быть пропорционально квадрату

градиента скорости $v_\lambda$ и соответствующему коэффициенту турбулентной

вязкости $\nu_{turb\lambda}\sim\lambda v_\lambda :$

$$

\varepsilon\sim\nu_{turb\lambda}\left(v_\lambda\over\lambda\right )^2

\sim {v_\lambda^3\over\lambda} ,

$$

откуда и получается (6).

Поставим теперь вопрос несколько иначе. Определим порядок величины $v_\tau$

изменения скорости в заданной точке пространства, испытываемого ею в течении

промежутка времени $\tau$, малого по сравнению с характеристическим временем

$T\sim l/u$ движения в целом. Для этого замечаем, что благодаря наличию общего

течения каждый данный участок жидкости в продолжение промежутка времени $\tau$

перемещается в пространстве на расстояние порядка произведения $\tau u$ средней

скорости $u$ на время $\tau .$ Поэтому в данной точке пространства по истечении

времени $\tau$ будет находиться участок жидкости, который в начальный момент был

удалён от этой точки на расстояние $u\tau .$ Искомую величину $v_\tau$ можно,

следовательно, получить, подставляя в (6) $\tau u$ вместо $\lambda :$

$$

v_\tau\sim(\varepsilon u\tau)^{1/3}.\eqno(7)

$$

От величины  $v_\tau$ следует отличать изменение $v'_\tau$ скорости данного

перемещающегося в пространстве участка жидкости. Это изменение может, очевидно,

зависеть только от величины $\varepsilon$, и, разумеется, от величины самого

интервала времени $\tau .$ Составляя из $\varepsilon$ и $\tau$ комбинацию

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

$$

v'_\tau\sim(\varepsilon\tau)^{1/2}.\eqno(8)

$$

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

квадратному, а не кубическому корню из $\tau .$ Легко видеть, что при $\tau\ll T$

изменение $v'_\tau$ всегда меньше изменения $v_\tau$\footnote {$^{22)}$}{Неравенство

$v'_\tau\ll v_\tau$ , по существу, уже подразумевалось при выводе (7).}.

С помощью выражения (1) для $\varepsilon$ можно переписать формулы (6),(7) в виде

$$

{v_\lambda\over\Delta u}\sim\left (\lambda\over l\right )^{1/3},\qquad

{v_\tau\over\Delta u}\sim\left (\tau\over T\right )^{1/3}.\eqno(9)

$$

В такой записи ясно видно свойство подобия локальной турбулентности: мелкомасштабные

характеристики различных турбулентных течений отличаются только масштабами измерения

длин и скоростей (или, что то же длин и времён).

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

расстояния $\lambda_0$ определяют собой в то же время порядок величины масштабов

наиболее мелкомасштабных пульсаций в турбулентном потоке (величину $\lambda_0$

называют {\it внутренним масштабом\/} турбулентности в противоположность

внешнему масштабу $l$). Для этого составляем ``локальное число Рейнольдса'':

$$

R_\lambda\sim{v_\lambda\lambda\over\nu}\sim

{\Delta u\cdot\lambda^{4/3}\over\nu l^{1/3}}\sim

R\left (\lambda\over l\right )^{4/3},

$$

где $R\sim\Delta u\cdot l/\nu$--- число Рейнольдса движения в целом. Порядок

величины $\lambda_0$ определяется тем, что должно быть $R_{\lambda_0}\sim 1.$

Отсюда находим

$$

\lambda_0\sim l/R^{3/4}.\eqno(10)

$$

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

величин $\varepsilon$ и $\nu :$

$$

\lambda_0\sim (\nu^3/\varepsilon)^{1/4}.\eqno(11)

$$

Таким образом, внутренний масштаб турбулентности быстро падает при увеличении

числа Рейнольдса. Для соответствующей скорости имеем

$$

v_{\lambda_0}\sim\Delta u/R^{1/4}. \eqno(12)

$$

Она тоже падает с увеличением $R$\footnote {$^{23)}$}{Формулы (10)-(12)определяют

законы изменения соответствующих величин с $R$. Что же касается количественной

стороны дела, то более правильным было бы писать в них отношение $R/R_0$ вместо

$R$.}.

Область масштабов $\lambda\sim l$ называют {\it областью энергии\/}; в ней

сосредоточена основная часть кинетической энергии жидкости. Значения $\lambda

\le\lambda_0$ составляют {\it область диссипации\/} ---

в ней происходит диссипация кинетической энергии. При очень больших значениях

$R$ обе эти области достаточно раздвинуты друг от друга, и между ними расположен

{\it инерционный интервал\/}, в котором

$$

\lambda_0\ll\lambda\ll l;

$$

к нему относятся излагаемые здесь результаты.

Закон Колмогорова--Обухова можно представить в эквивалентной спектральной

(по пространству) форме. Введём вместо масштабов $\lambda$ соответствующие

``волновые числа'' пульсаций $k\sim 1/\lambda ,$ и пусть $E(k)dk$ есть кинетическая

энергия (единицы массы жидкости), заключённая в пульсациях со значениями $k$

в заданном интервале $dk.$ Функция $E(k)$ имеет размерность см$^3/$с$^2$;

составляя комбинацию этой размерности из $\varepsilon$ и $k$, получим

$$

E(k)\sim\varepsilon^{2/3}k^{-5/3}.\eqno(13)

$$

В эквивалентности этой формулы закону (6) легко убедиться, заметив, что квадрат

$v_\lambda$ определяет порядок величины суммарной энергии, заключённой в пульсациях

со всеми масштабами порядка и меньше заданного значения $\lambda .$ К этому же

результату мы придём, интегрируя выражение (13)

$$

\int\limits_k^\infty E(k)dk\sim {\varepsilon^{2/3}\over k^{2/3}}\sim {(\varepsilon\lambda)^{2/3}}\sim v_\lambda^2 .

$$

Наряду с пространственными масштабами турбулентных пульсаций, можно рассматривать

также и их временные характеристики --- частоты. Нижний конец частотного спектра

турбулентного движения лежит при частотах $\sim u/l .$ Верхний же его конец

определяется частотами

$$

\omega_0\sim{u\over\lambda_0}\sim{{u\over l}R^{3/4}},\eqno(14)

$$

отвечающими внутреннему масштабу турбулентности. Инерционной области отвечают

частоты в интервале

$$

{u\over l}\ll\omega\ll {u\over l}R^{3/4}.

$$

Неравенство $\omega\gg u/l$ означает, что по отношению к локальным свойствам

турбулентности основное движение можно считать стационарным. Распределение энергии

по частотному спектру в инерционной области получается из (13) заменой $k\sim\omega /u:$

$$

E(\omega)\sim (u\epsilon)^{2/3}\omega^{-5/3}, \eqno(15)

$$

причём $E(\omega)d\omega$ есть энергия, заключённая в частотном интервале $d\omega$.

Частота $\omega$ определяет период повторяемости во времени движения в данном участке

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

частоты (обозначим её $\omega '$), определяющей период повторяемости движения

в данном перемещающемся в пространстве участке жидкости. Распределение энергии по

спектру этих частот не может зависеть от $u$, и должно определяться только параметром

$\varepsilon$  и самой частотой $\omega ' .$ Снова из соображений размерности

найдём, что

$$

E(\omega ')\sim{\varepsilon/\omega '^2 .}

$$

Эта формула находится в таком же отношении к закону (15), как (8) к (7).

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

находящихся первоначально вблизи друг от друга. Рассмотрим две жидкие частицы на

малом (в инерциальной одласти) расстоянии $\lambda .$ Снова руководствуясь

соображениями размерности, можно заключить, что скорость изменения этого расстояния

со временем

$$

{d\lambda\over dt}\sim (\varepsilon\lambda)^{1/3}.

$$

Интегрируя это соотношение, найдём, что время $\tau$, в течение которого две

частицы, находившиеся первоначально на расстоянии $\lambda_1$ друг от друга,

разойдутся на расстояние $\lambda_2\gg\lambda_1 ,$ равно по порядку величины

$$

\tau\sim {\lambda_2^{4/3}/\varepsilon^{1/3}}.

$$

Обратим внимание на самоускоряющийся характер процесса: скорость расхождения

растёт с увеличением $\lambda .$ Это свойство связано с тем, что к расхождению

частиц, находящихся на расстоянии $\lambda$, приводят только пульсации масштабов

$\le \lambda$; пульсации больших масштабов переносят обе частицы вместе и не

приводят к их расхождению.

Наконец, остановимся на свойствах движения в участках с размерами $\lambda\ll\lambda_0.$

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

плавно. Поэтому можно разложить здесь $v_\lambda$ постепеням $\lambda$ и,

сохранив только первый член, получим $v_\lambda ={\rm const}\cdot\lambda$.

Коэффициент определяется требованием, чтобы при $\lambda\sim\lambda_0$ было

$v_\lambda\sim v_{\lambda_0}$.

Таким образом находим

$$

v_\lambda\sim {{v_{\lambda_0}\over\lambda_0}\lambda}\sim {{\Delta u\over l}\lambda R^{1/2}}.

$$

Этот результат можно получить также и путём приравнивания двух выражений для диссипации

энергии $\varepsilon$: выражения $(\Delta u)^3/l$ (1), определяющего $\varepsilon$

через характеристики

крупномасштабных пульсаций, и выражения $\nu (v_\lambda /\lambda)^2$, определяющего ту же величину через

градиент скорости тех пульсаций, в которых фактически и происходит диссипация.

\beginsection{Заключение}

В заключение остаётся пожалеть о том, что за рамками реферата остались многие

современные достижения теории турбулентности; как то плодотворная идея Каданова

о {\it скейлинге\/}, $\varepsilon$-разложения Вильсона и ренорм-групповые методы,

развитые в квантовой теории поля и с успехом применяемые в теории развитой

турбулентности. Но опять таки ``плюнь в глаза тому кто скажет,

что можно объять необъятное'', засим я завершаю этот реферат, полагая свою

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

%\supereject

\vfill

\break

\beginsection{Список использованной литературы}

\frenchspacing

\item{[1]}{\it Л. П. Каданов\/} ``Пути к хаосу''

Перевод статьи: {\it Kadanoff L. P.} --- Physics Today, December 1983, p. 43.

Физика за рубежом. Серия {\bf А}. Исследования. Сборник статей. Москва,

``Мир'', 1985.

\item{[2]}{\it Л. Д. Ландау, Е. М. Лившиц\/} Теоретическая физика. Том VI.

 ``Гидродинамика'', Москва, ``Наука'', 1988.

\item{[3]}{\it V. S. L'vov, I. Procaccia\/} ``Turbulence: a universal problem.''

 Physics World, 9 (8) 35-40 (1996).

\item{[4]}{\it А. С. Монин, А. М. Яглом\/} ``Статистическая гидромеханика. Том 2'',

Санкт-Петербург, ``Гидрометеоиздат'', 1996.

\bye