Аттрактор Лоренца
artpgchЧисленное моделирование системы Лоренца
Введение
В данном примере представлено численное моделирование системы Лоренца — классического примера динамической системы, демонстрирующей хаотическое поведение. Система была предложена Э. Лоренцем в 1963 году в рамках задачи о конвекции в подогреваемом снизу слое жидкости.
Математическое описание
Она описывается системой трёх обыкновенных дифференциальных уравнений:
Переменные модели имеют следующую интерпретацию в контексте исходной физической задачи:
-
— интенсивность конвективного движения;
-
— разность температур между восходящим и нисходящим потоками;
-
— отклонение профиля температуры от линейного равновесного распределения.
Система содержит три безразмерных управляющих параметра:
- (число Прандтля) — отношение кинематической вязкости к теплопроводности среды. Определяет соотношение между скоростями диссипации механических и тепловых возмущений.
- (число Рэлея) — параметр, характеризующий интенсивность нагрева. При ρ<1 система имеет единственную устойчивую особую точку; с ростом ρ происходят бифуркации, приводящие к хаотической динамике.
- — геометрический параметр, определяемый отношением характерных размеров конвективной ячейки.
При значениях параметров , , фазовая траектория системы формирует странный аттрактор — притягивающее множество в фазовом пространстве, обладающее фрактальной структурой. Траектория переключается между двумя областями, соответствующими вращению жидкости в противоположных направлениях. Основными свойствами системы в этом режиме являются:
-
экспоненциальная чувствительность к начальным условиям, делающая невозможным долгосрочное прогнозирование;
-
наличие двух центров вращения с нерегулярными переходами между ними.
Определим шаг интегрирования.
Δt = 0.01
Определим исходные параметры, инициализируем массивы и зададим начальные условия.
σ = 10
β = 8/3
ρ = 28
N = 30000
x = zeros(N)
y = zeros(N)
z = zeros(N)
x[1] = 1
y[1] = 1
z[1] = 1
Выполним численное интегрирование явным методом Эйлера.
for i in 1:N-1
dx = σ * (y[i] - x[i])
dy = x[i] * (ρ - z[i]) - y[i]
dz = x[i] * y[i] - β * z[i]
x[i+1] = x[i] + Δt * dx
y[i+1] = y[i] + Δt * dy
z[i+1] = z[i] + Δt * dz
end
Построим аттрактор в трёхмерном евклидовом пространстве.
plot(x, y, z, linewidth = 0.5, legend = false,
aspect_ratio = :equal, grid = true,
title = "Аттрактор Лоренца", xlabel = "x", ylabel = "y", zlabel = "z",
camera = (30, 30), size = (700, 600))
Заключение
В данном примере проведено численное моделирование системы Лоренца при различных режимах управляющих параметров. Построены фазовые траектории, соответствующие классическому странному аттрактору. Несмотря на то что система Лоренца была получена как упрощённая модель атмосферной конвекции, её значение выходит далеко за пределы метеорологии. Основная ценность данной модели состоит в демонстрации принципиальной возможности хаотического поведения в детерминированных системах малой размерности, что имеет применение в следующих областях:
- Прогнозирование погоды и климата. Система Лоренца стала одной из первых математических моделей, показавших ограниченность горизонта предсказуемости атмосферных процессов. На практике это привело к развитию ансамблевых методов прогнозирования, когда вместо одного детерминированного расчёта выполняется серия расчётов с немного различающимися начальными условиями, а результат интерпретируется в вероятностных терминах.
- Гидродинамика и тепломассообмен. Анализ перехода к турбулентности в конвективных течениях используется при проектировании теплообменников, систем охлаждения электронной аппаратуры, а также при изучении мантийной конвекции в геофизике.
- Лазерная физика. Уравнения, математически эквивалентные системе Лоренца, описывают динамику одномодового лазера. Анализ хаотических режимов генерации применяется при разработке источников излучения с заданными характеристиками.
- Химическая кинетика и биология. Хаотические режимы, аналогичные наблюдаемым в системе Лоренца, возникают в автокаталитических реакциях, моделях динамики популяций и нейронной активности.
- Криптография и защита информации. Свойство экспоненциальной чувствительности к начальным условиям используется при построении генераторов псевдослучайных последовательностей и систем шифрования на основе синхронизации хаотических генераторов.
Таким образом, численное исследование системы Лоренца представляет собой не только иллюстрацию фундаментальных свойств нелинейной динамики, но и основу для понимания широкого круга явлений, встречающихся в инженерной практике и естественных науках.