CC BY
15
92
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2009. № 6
КРАТКИЕ СООБЩЕНИЯ ТЕОРЕТИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ ФИЗИКА
Распределение энергии в сильно возмущенной хаотической системе
Ю. Е. Щадилова
Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра квантовой электроники. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
E-mail: аyes@shg.ru
Статья поступила 04.03.2009, подписана в печать 07.09.2009
Для модели осциллятора Дуффиига, возмущенного сильным гармоническим полем и совершающего хаотическое движение, получено аналитическое выражение для распределения по энергии как функции амплитуды и частоты возмущения.
Ключевые слова: нелинейная динамика, хаос, осциллятор Дуффинга, распределение энергии. УДК: 531.395. PACS: 05.45.-а.
Гамильтоновы нелинейные неавтономные системы общего вида, гамильтонианы которых H(p,q,t) явно зависят от времени, а функция H(p,q,t) сложнее квадратичного полинома по р и q, можно разделить на два класса. Это системы, способные поглотить неограниченное количество энергии от источника внешней силы, хаотическое движение которых может быть инфинитно по действию (или по энергии), а также системы, у которых энергия вечно остается ограниченной. К первому классу относятся, например, ротатор с 5-толчками — механический прототип стандартного отображения [1] и модель стохастической ионизации атома водорода [2]. Ко второму классу относятся модели одномерных нелинейных осцилляторов с жесткой нелинейностью (частота собственных колебаний является возрастающей функцией энергии) во внешнем гармоническом поле [3] и модель Улама для стохастического ускорения Ферми [4]. В системах второго класса существуют стационарные функции распределения динамических величин. Задача аналитического описания таких функций представляет интерес, так как они позволяют определять зависимости моментов функций стохастической компоненты (например, среднее значение энергии) от параметров модели.
Целью настоящей работы является определение функции распределения Р(Е) по энергии для модели осциллятора Дуффинга во внешнем гармоническом поле большой амплитуды, заданной уравнением
х - х + х3 = F sinujt, (1)
где F > 1, ш « 1. Колебания осциллятора Дуффинга были впервые описаны в монографии [5]. Эта модель стала одной из центральных в теории нелинейных колебаний; новую жизнь модели дало открытие в ней хаотического движения [6], исследование которого активно продолжается до настоящего времени [7-10].
Под энергией Е мы понимаем величину энергии невозмущенной системы Е = х2/2-х2/2 + хА/4, взятую в моменты времени tn = тгп/ш (п — целое). В моменты tn возмущение обращается в нуль.
Схожая задача была решена Либерманом и Лих-тенбергом [11]. Эти авторы описали распределение
по скоростям w(u) для модели Улама — одномерного движения точки между двумя упруго отражающими ее стенками, одна из стенок осциллирует по периодическому закону. Распределение было получено, основываясь на малости изменения энергии за период поля, что позволило использовать для описания кинетики энергии модель энергетической диффузии. В моделях одномерных нелинейных осцилляторов с жесткой нелинейностью, например в модели осциллятора Дуффинга (1) в сильном поле, за период внешнего поля энергия может изменяться в десятки раз, что требует иного подхода.
Рассмотрим проекцию траектории частицы на стробоскопическое сечение Пуанкаре {х;х} при фиксированных значениях частоты и амплитуды внешней силы (рис. 1). Точки заполняют большую часть области Е(х,х) <Е+, где Е+ — максимально достижимое значение энергии системы. Поэтому примем для простоты, что они заполняют всю область, и притом равномерно, последнее вытекает из сохранения фазового объема. Это предположение, эквивалентное малости островков устойчивости в сравнении с хаотической компонентой, является главным ограничителем применимости нашего расчета. Теперь отыскание функции распределения по энергии сводится к вычислению зависимости Е+ от параметров внешнего поля F и ш. Заметим, что значение Е+, которое для нашей модели дает критерий перекрытия резонансов Чирикова [12], отличается от найденного при численном интегрировании системы (1) в десятки раз (например, при параметрах внешнего поля F = 2.5, ш= 1 значение Е+, найденное численным интегрированием, равно 10.5, использование критерия перекрытия резонансов Чирикова приводит к значению 1.5).
При больших значениях энергии Е > 1 движение невозмущенного осциллятора Дуффинга слабо ангармо-нично — отношение амплитуд третьей и первой гармоник составляет лишь несколько процентов (например, при начальных условиях х = 0, х = 4, которые соответствуют энергии Е = 8, это отношение равно 0.058). Поэтому в низшем приближении движение может быть
ТЕОРЕТИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ ФИЗИКА
93
-1 0 1 Координата, х
Рис. 1. Стробоскопическое отображение фазовой траектории осциллятора Дуффинга (1) на плоскость {х;х} при параметрах внешней силы Т7 = 2.5, и = 1. Сплошная линия — граница области фазового пространства, доступной при энергиях системы, не превосходящих Е+
описано гармоническим законом х^)=Асо$Ш. В области Е ^ 1 амплитуда А и частота О слабо зависят от энергии: А ос О ос Эта слабая зависимость
позволяет считать амплитуду и частоту в возмущенной системе медленными функциями времени и представить ее движение в виде
х(0 =А(Е(0) соб
П(Е(?))Ж' + <р о =
= Л(0собФ(0. (2)
Скорость изменения энергии дается выражением Е = х^)Е ът ^. Подставляя в него закон движения (2) и пренебрегая членами высших порядков малости, (АО,/А) < 1, получаем
7Т/Ш
Е = Ел-
А£1 эт Ф/7 ът^сИ,
(3)
где £о ~~ начальная энергия.
Максимальную энергию Е+ оценим из формулы (3), используя прямой вариационный метод. Возьмем простейшую пробную функцию вида
Ф(0 = Ш + ¥>о,
(4)
где О — частота усредненного по половине периода движения, — начальная фаза движения.
Подставив (4) в (3), вычислив функционал энергии и продифференцировав его по начальной фазе, получаем
Е(Е0,Е,и)=Е0-2
ЕиАй
?Р
■ ОЛ
7Г О 9 С°Б I о ~
>2 \ 2 и
где А(Ео) — амплитуда усредненного по половине периода движения. Используя асимптотики амплитуды и частоты сильно возмущенного движения и принимая
в качестве начальной энергию движения свободной частицы во внешнем поле Ео = Е/ио2, получаем
' ' — и3
7Г ХЕ1/А\ соз,.-^),
(5)
где ¡1 = л/2 « 1.414, Л = ^Г(1/4)Г-1(3/4) « 1.198. Зависимости энергии Е+ от частоты и и амплитуды Е внешней силы, рассчитанные по формуле (5), представлены на рис. 2. Как видно их них, данные аналитического расчета и численного эксперимента находятся в 10%-м согласии в области параметров внешнего поля Е= 1,...,4 и и = 0.5,..., 1.2, где островки устойчивости не превосходят области, занятой хаотической компонентой.
2 4 6
Амплитуда внешней силы, Т7
15
2 о
Н 10
0.4
0.6 0.8 1.0 Частота внешней силы, со
1.2
Рис. 2. Зависимость максимальном энергии от амплитуды И внешней силы при частоте ш = 1 (а) и от частоты ш внешней силы при амплитуде И = 2.5 (б). Сплошная линия — расчет по формуле (5), точки — результат численного эксперимента
Поскольку в канонических переменных {х; росх} для гамильтоновых систем имеет место сохранение фазового объема, то плотность вероятности в хаотической области на стробоскопическом сечении постоянна [13, с. 466]. Вероятность Р(Е) найти систему в состоянии с энергией, не превосходящей Е,
94
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2009. № 6
§ 0.12 я
6 0.08 и
п
0
и
1 0.04
Я О
5
Э
6
| 0 2 4 6 8 10 ^ Энергия, Е
Рис. 3. Вид функции распределения по энергии при Р = '2, ш=1. Сплошная линия — расчет по формуле (6), точки — результат численного эксперимента
пропорциональна площади 5(£), ограниченной кривой к2/2 - х'2/2 + л,4/4 = Е. Эта площадь для модели двухъямного осциллятора Дуффинга (1) может быть выражена аналитически:
Я(£)= уЧ + 4Е ((VI + 4£ - 1) К(к) + 2Е(к)) ,
где Е(к) и К(к) — полные эллиптические интегралы первого и второго рода соответственно [14],
к2 = (1 + VI + / (2VI +4£) — безразмерный параметр. Тогда функция распределения
Р(Е) = Ы(Е1)5(Е), (6)
где Ы(Е |) — нормировочный коэффициент. Характерный вид функции Р{Е), полученный по формуле (6) с использованием оценки (5) для максимальной энергии, представлен на рис. 3. Из распределения по энергии видно, что при удовлетворительном общем согласии результаты, полученные по формуле (6), дают
завышенную вероятность в области высоких энергий. Это обстоятельство легко понять из графика на рис. 1: области регулярного движения (которые мы в нашем подходе присчитывали к хаотической компоненте) лежат вблизи граничной кривой х2/2 — х'2/2 + х4/А = Е\ .
Автор благодарит П. В. Елютина за интерес к работе и А. А. Никулина за полезные обсуждения.
Список литературы
1. Заславский Г.М., Сагдеев Р.З., Усиков Д.А., Черников A.A. Слабый хаос и квазирегулярные структуры. М., 1991.
2. Делоне Н.Б., Крайнов В.П., Шепелянский Д.Л. // УФН. 1983. 140. С. 355.
3. Chism W., Timberlake T., Reichl L.E. // Phys. Rev. E. 1988. 58, N 2. P. 1713.
4. Заславский Г.M., Сагдеев Р.З. Введение в нелинейную физику: от маятника до турбулентности и хаоса. М., 1988.
5. Duffing G. Erzwungene Schwingungen bei veränderlicher Eigenfrequenz und ihre technische Bedeutung. Berlin, 1918.
6. Holmes P. // Philos. Trans. R. Soc. London, Ser. A. 1979. 292, N 1394. P. 419.
7. Lai Y.-С., Kandangath A., Krishnamoorthy S. et al. // Phys. Rev. Lett. 2005. 94, N 21. P. 214101.
8. Wu B.S., Suna W.P., Limb C.W. // J. of Sound and Vibration. 2007. 307, N 3-5. P. 953.
9. Arie Kapulkinand Arjendu K. Pattanayak // Phys. Rev. Lett. 2008. 101, P. 074101.
10. Buric N. 11 Phys. Rev. A. 2009. 79. P. 022101.
11. Lieberman M.A., Liechtenberg A.J. // Phys. Rev. A. 1972. 5. P. 1852.
12. Chirikov B.V. Ц Phys. Reports. 1979. 52, N 5. P. 263.
13. Лихтенберг А., Либерман M. Регулярная и стохастическая динамика. М., 1984.
14. Корн Г., Корн. Т. Справочник по математике. М., 1973.
Energy distribution in strongly perturbed chaotic system Yu.E. Shchadilova
Department of Quantum Electronics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow, 119991, Russia. E-mail: a yes@shg.ru
The analytical expression for the energy distribution for the Duffing oscillator that moves chaotically under the strong harmonic perturbation is obtained as a function of the perturbation amplitude and frequency.
Keywords: nonlinear dynamics, chaos, Duffing's oscillator, energy distribution. PACS: 05.45.-a. Received 4 March 2009.
English version: Moscow University Physics Bulletin 6(2009).
Сведения об авторе
Щадилова Юлия Евгеньевна — студентка; тел.: (495) 939-36-69, e-mail: yes@shg.ru.
