УДК 550.3 ГАМИЛЬТОНОВА МОДЕЛЬ КОЛЕБАНИЯ МАГНИТНОЙ

advertisement
УДК 550.3
ГАМИЛЬТОНОВА МОДЕЛЬ КОЛЕБАНИЯ
МАГНИТНОЙ СТРУНЫ
Лопатина А. А.
научный руководитель д-р физ.-мат. наук Еркаев Н. В.
Сибирский федеральный университет
Равновесные конфигурации магнитного поля в электропроводящей среде
характеризуются балансом сил Лоренца и градиента давления. В двумерном случае
известны точные аналитические решения, описывающие распределения равновесных
компонент магнитного поля и давления среды. На рис. 1 показаны распределения
магнитных силовых линий, соответствующих одному из наиболее известных
равновесных решений такого рода, полученному Каном.
A
Рис. 1. Равновесное состояние силовых линий
магнитного поля
Это решение, в частности, применяют для моделирования магнитного поля в
хвосте магнитосферы Земли. В случае достаточно высокой электропроводности
магнитные силовые линии непосредственно связаны с частицами среды, смещения
которых приводят к деформациям силовых линий. Показано, что перемещения тонкого
пучка магнитных силовых линий, называемого магнитной трубкой или магнитным
филаментом, могут описываться уравнениями магнитной струны:
𝜕 2 𝑧 𝜕 𝜌𝜕𝑧
𝜕П
−
(
)+
= 0,
2
𝜕𝑡
𝜕𝛼 𝜕𝛼
𝜕𝑧
𝜕 2𝑥
𝜕 𝜌𝜕𝑥
𝜕П
−
(
)+
=0
2
𝜕𝑡
𝜕𝛼 𝜕𝛼
𝜕𝑥
(1)
Уравнения (1) можно записать в форме Гамильтона:
𝜕𝑟̇⃗
𝛿𝐻
𝜕𝑟⃗ 𝛿𝐻
=−
,
=
,
(2)
𝜕𝑡
𝛿𝑟⃗
𝜕𝑡 𝛿𝑟̇⃗
где функция Гамильтона выражается интегралом вдоль магнитной силовой линии:
𝛼0
𝜕𝑟 2 1
𝜕𝑟 2 1 𝐾𝜌(𝑘−1)
𝐻 = ∫ (( ) ∙ + 𝜌 ∙ ( ) ∙ +
+ П/𝜌)𝑑𝛼 .
𝜕𝑡
2
𝜕𝛼
2
𝑘−1
−𝛼0
(3)
Здесь 𝛼 - переменная Лагранжа, меняющая вдоль струны, 𝜌 – плотность вещества,
𝜕𝑟
𝜕𝑡
-
скорость, K – постоянная, зависящая от температуры,  - адиабатическая постоянная
адиабаты,
- полное фоновое давление, определяемое по известному равновесному
состоянию, r(𝛼) - вектор функция, определяющая пространственную конфигурацию
магнитной силовой линии.
Перемещения вдоль оси 𝑥 и 𝑧 представляем виде суммы:
𝑥 = 𝑥0 (𝛼, 𝑡) + 𝑥1 (𝛼, 𝑡), 𝑧 = 𝑧0 (𝛼, 𝑡) + 𝑧1 (𝛼, 𝑡),
где: 𝑥0 (𝛼, 𝑡) и 𝑧0 (𝛼, 𝑡) начальные невозмущенные координаты, а 𝑥1 (𝛼, 𝑡) и 𝑧1 (𝛼, 𝑡) возмущения, которые находим с помощью разложения в ряды Фурье:
∞
𝜋 ∙ (2(𝑘 − 1) + 1) ∙ 𝛼
𝑧1 (𝛼, 𝑡) = ∑ 𝑎𝑘 cos [
],
2 ∙ 𝛼0
𝑘=1
∞
𝜋 ∙ (2(𝑘 − 1) + 1) ∙ 𝛼
𝑥1 (𝛼, 𝑡) = ∑ 𝑏𝑘 cos [
],
2 ∙ 𝛼0
𝑘=1
где: 𝛼0 - заданный параметр, связанный с длиной рассматриваемого отрезка магнитной
силовой линии. Рассматривая режим достаточно медленных колебаний, пренебрегаем
изменениями плотности и полагаем  = 1 в безразмерных переменных.
Для коэффициентов Фурье получаем систему обыкновенных дифференциальных
уравнений:
1
𝑑𝑎𝑘̇
𝜋2
𝜋
2
(𝑥
− (2𝑘 + 1)
𝑎
+
∫
[𝑓
,
𝑧
,
𝛼)
∙
cos
(
∙ 𝛼 ∙ [2𝑘 + 1])] 𝑑𝛼 ,
𝑘
1
1
1
𝑑𝑡
2
4 ∙ 𝛼02
(4)
−1
1
𝑑𝑏𝑘̇
𝜋2
𝜋
2
(𝑥
− (2𝑘 + 1)
𝑏
+
∫
[𝑓
,
𝑧
,
𝛼)
∙
cos
(
∙ 𝛼 ∙ [2𝑘 + 1])] 𝑑𝛼 ,
𝑘
2
1
1
𝑑𝑡
2
4 ∙ 𝛼02
(5)
−1
где:
𝑓1 (𝑥1 , 𝑧1 , 𝛼) =
𝜕П(𝑥0 + 𝑥1 , 𝑧0 + 𝑧1 ) 𝜕П(𝑥0 , 𝑧0 )
−
,
𝜕𝑧
𝜕𝑧0
𝑓2 (𝑥1 , 𝑧1 , 𝛼) =
𝜕П(𝑥0 + 𝑥1 , 𝑧0 + 𝑧1 ) 𝜕П(𝑥0 , 𝑧0 )
−
.
𝜕𝑥
𝜕𝑥0
Система уравнений ( 4, 5 ) интегрировалась численно с использованием пакета
MATHCAD для N = 10 при начальных условиях:
1
𝜋
𝜋 ∙ (2𝑘 + 1) ∙ 𝛼
𝑧1 (𝛼, 0) = 𝐴 ∙ cos (𝛼 ∙ ) , 𝑎𝑘 = ∫ 𝑧1 (𝛼, 0) ∙ cos [
] 𝑑𝛼, 𝑏𝑘 = 0. (6)
2
2
−1
Результаты расчетов представлены на рис. 2, 3.
Рис. 2. Зависимость координаты (z) и вертикальной
скорости (Vz) от времени
Рис. 3. Конфигурации магнитной струны в различных фазах колебаний:
1– при максимальном смещение вверх, 2 – в среднем положение, 3- при
максимальном смещении вниз.
300
250
T
200
0
0.2
0.4
0.6
0.8
1
А
Рис. 4. Зависимость периода колебаний T от
амплитуды A.
Минимальное значение периода зависит от второй производной полного
давления при z = 0 и хорошо аппроксимируется простой аналитической формулой:
1 𝜕 2П
𝜏 = 2𝜋√
= 174
𝜌 𝜕𝑧 2
(7)
Найденное по графику (рис. 4) значение периода при амплитуде 0,5 равно 177,
что довольно близко значению, определяемому формулой (7).
Выводы. На основе Фурье аппроксимации исследованы колебания магнитной
струны конечной амплитуды. Построены конфигурации магнитной струны для
различных фаз колебаний. Найдена зависимость периода колебаний от амплитуды.
Результаты расчетов показывают, что колебания изохронны до амплитуд 0.2.
Предложена аналитическая формула для оценки минимального периода. При больших
амплитудах начинается заметный рост периода колебаний. Объектом приложения
разработанной модели являются так называемые флэппинг колебания, наблюдаемые в
хвосте магнитосферы Земли. Модель позволяет интерпретировать имеющиеся
наблюдения и дает адекватную оценку периода флэппинг колебаний.
Download