Динамический хаос (ИПИС, ФКС)
.pdfТаким образом, модель Лоренца описывает динамическую систему с трёхмерным фазовым пространством, мгновенное состояние которой определяется набором трёх переменных (x, y, z) , а оператор эволюции задан уравнениями (110).
Физический смысл величин, фигурирующих в уравнениях Лоренца, можно установить из соотношений (84) и (86):
x – характеризует скорость вращения конвекционных валов; y и z – отвечают за распределение температуры, соответственно, по горизонтали и по вертикали;
b – определяется геометрией конвекционной ячейки, а
именно, отношением её вертикального и горизонтального |
|
размеров a h l ; |
|
|
– число Прандтля; |
r R R |
– отношение числа Рэлея R к его критическому |
значениюc |
R . |
|
c |
При численном решении системы (110) на компьютере Лоренц обнаружил, что при определённых значениях параметров b , и r в исследуемой системе устанавливался хаотический режим движения. Эта хаотичность в движении являлась следствием нелинейности системы уравнений Лоренца. Возможность существования хаотических траекторий у детерминированной системы стало важным открытием, повлиявшим впоследствии на дальнейшее развитие теории динамического хаоса.
Рис. 62. Зависимости динамических переменных x , y и z от
времени, полученные численным интегрированием уравнений Лоренца при 10 , b 83, r 28 .
Рис. 107. Чувствительность решений уравнения Лоренца к
изменениям начальных условий в хаотическом режиме: r 28, |
||||
10 и b 8 3 . |
Начальные |
условия: |
x0 |
10, 485 , |
y0 10, 485 , z0 |
27 , где |
, и |
|
– случайные |
добавки порядка 0,01. |
|
|
|
Физические задачи, сводящиеся к системе Лоренца
Оказывается, помимо задачи о конвекции в подогреваемом снизу слое жидкости, существует ещё много других физических систем, в которых происходящие в них процессы описываются теми же самыми уравнениями Лоренца (110). К таким задачам, например, относят задачу о конвекции в замкнутой петле, задачу о водяном колесе, задачу о динамике одномодового лазера и т.п.. Рассмотрим их более подробно.
Конвекция в замкнутой петле
Рассмотрим замкнутую в кольцо трубку, заполненную жидкостью (Rubenfeld, Siegman, 1977). Трубка подогревается снизу и охлаждается сверху (см. рис. 58). При достаточно большой интенсивности подогрева в такой трубке возможно возникновение конвекционного течения жидкости.
Рис. 58. Задача о конвекции в замкнутой кольцевой трубке. X характеризует скорость течения, Y и Z – отклонение температуры от среднего значения в указанных точках.
Пусть – угловая координата, отсчитываемая от нижней точки трубки против часовой стрелки. Тогда зависимость температуры от угла T T ( ) должна быть периодической функцией с периодом 2 , поэтому её можно представить в виде ряда Фурье. Ограничиваясь учётом только первых членов разложения (первой гармоники), можем записать
T ( ) T0 (1 Y sin Z cos ) , |
(111) |
|||
где величины Z и Y |
характеризуют отклонение темпе- |
|||
ратуры от среднего значения T0 в |
нижней точке трубки |
|||
(при 0 ) и |
в крайней правой |
точке (при 2 ), |
||
соответственно. |
Через |
X обозначим скорость |
течения |
жидкости. Поскольку жидкость практически несжимаема, то эту скорость можно считать одной и той же в любом сечении трубки, т.е. не зависящей от .
Сконструируем уравнения, описывающие изменения во времени динамических переменных X , Y , Z , из простых качественных соображений.
Если провести тщательный анализ движения жидкости в трубке, то можно заключить, что изменение её скорости течения обуславливается архимедовой силой, которая
оказывается пропорциональной величине Y , |
а также си- |
лой вязкого трения, пропорциональной X . Поэтому первое |
|
уравнение будет иметь вид |
|
|
|
X cY X , |
(112) |
где c и – некоторые постоянные.
Теперь составим уравнения для Y и Z . Предположим
сначала, что имеет место течение с постоянной скоростью
X Const . Тогда можно записать
T T ( X t) T0 1 Y sin( X t) Z cos( X t) .
Вводя обозначение X t , имеем
|
|
|
T T (Y cos Z sin ) T ( XY cos X Z sin ) . |
||
|
0 |
0 |
Отсюда видно, что перенос температуры потоком жидкости учитывается членом вида X Z в уравнении для Y и членом ( XY ) в уравнении для Z . В оба уравнения следует добавить ещё члены, учитывающие релаксацию, которые должны иметь вид, соответственно, ( DY ) и ( DZ) , где D – константа, а также необходимо учесть подогрев добавлением постоянного члена A в уравнение для Z .