ЭЛЕКТРОНИКА И ПРИБОРОСТРОЕНИЕ
УДК 621.385.01
В.Б. Байбурин, С. А. Трачук
ХАОТИЧЕСКАЯ ДИНАМИКА ВЗАИМОДЕЙСТВУЮЩИХ ЗАРЯДОВ
В СКРЕЩЕННЫХ ПОЛЯХ
Проведено численное исследование модели поведения двух зарядов в скрещенных полях. На основе расчета показателей Ляпунова определены области существования регулярного и хаотического поведения систем, получены карты показателей Ляпунова.
V.B. Bayburin, S.A. Trachuk
CHAOTIC DYNAMICS ASSOCIATED CHARGES IN CROSSED FIELDS
Computational solution behavior of model of two charges in crossed fields implemented in the article. On the basis of calculation of Lyapunov exponents it was determined the region of existing regular and chaotic behavior of systems, were received dynamic conditions map.
Особенности движения одиночного заряда в скрещенных электрическом и магнитном полях (при неоднородном магнитном поле) рассматривались в работе [1]. Представляет интерес рассмотреть движение заряда ^ в скрещенных полях в присутствии второго заряда д2, взаимодействующего с первым. Предположим, что оба заряда отрицательны и одинаковы по величине. Следовательно, они будут отталкиваться. Рассмотрим схему в скрещенных полях, аналогичную диоду с продольным (вдоль оси катода) магнитным полем и радиальным электрическим полем Ег (рис. 1).
В прямоугольной системе координат уравнения движения двух зарядов применительно к схеме на рис. 1 можно записать в следующем виде:
> к У гЛ \
\ а / х еёу
Рис. 1
“2 У2
Х1 _ У'; Х2 _ V,
х2
Vy
-^ + Еу, + Ру12; V У2 = -а V + Еу2 + Р
х2
У2
У 21
У1 = V,; у 2 = V
У2
где х1, у1 - координаты первого заряда в плоскости (х,у); х2, у2 - координаты второго заряда; уч,ул - скорость движения первого заряда; УХ2,УУ2 - скорость движения второго заряда;
О - циклотронная частота; ЕХ1, Е - составляющие напряженности радиального электрического поля Ег, действующие на первый заряд; ЕХ2, Еу2 - на второй заряд.
При этом:
ЕУ1 =-
и 0
( г л
а
1п
У1
(х1 )2 +(У1 )'
ЕХ1 =-
и0
Х1
1п
ЕУ2 =■
V гк У
и
( г Л
а
V гк У
(Х1 )2 + (У1 )
1п
Г г Л
а
V Г У
У 2
(х2 )2 +(У2 )
ЕХ2 =■
и
Г г Л
а
1п
(Х2 )2 + (У2 )
V Г У
где Г£ - радиус катода; га - радиус анода; и0 - величина, соответствующая потенциалу анода.
Силы Р , ^у12 и Р, ^у21 - силы электростатического взаимодействия зарядов.
Принимая, что заряды д1 и д2 равны 1 для сил электростатического взаимодействия, можно записать:
РХ12 =
Х2 - Х1. р _
; У12
У2 - У1
,3
г
Р _ Х2 - Х1- р
Х21 ,„3 ’ * У21
У2 - У1
г
Для численного решения уравнения движения и определения величин фазовых переменных использовался неявный метод, изложенный в работе [2]. При этом принималось: а=1, и0=1, гк=0,5, га=10, нулевые начальные скорости зарядов. Кроме того, сам катод считался «прозрачным» для движущихся зарядов, а радиальное электрическое поле внутри катода принималось равным полю на катоде. Как показали расчеты, это допущение не влияет существенно на характер карты динамических режимов в пространстве «катод-анод».
На рис. 2, 3 и 5 представлены области существования хаотических и регулярных траекторий зарядов в координатах (х,у). Эти области («ляпуновское пространство») получены на основании расчета показателей Ляпунова (А) и окрашены в различные оттенки серого цвета. При этом белый цвет соответствует регулярным траекториям (показатель Ляпунова меньше нуля). Градации темного цвета соответствуют различным положительным значениям показателя по мере роста его величины.
Случай, когда присутствует только один заряд, представлен на рис. 2. Видно, что практически вся область в этом случае имеет белый или светло-серый цвет, что соответствует регулярным траекториям.
Х
2
3
г
3
г
Рис. 2. Карта показателей Ляпунова в плоскости координат (х,у) в случае присутствия одного заряда. Максимальное значение показателя Ляпунова Атах=4.030364, минимальное значение показателя Ляпунова Атіп=-0.568831
Карты динамических режимов для случаев, когда учтена сила взаимодействия двух зарядов, представлены на рис. 3, 4:
Рис. 3. Карта показателей Ляпунова в плоскости координат (х,у).
Случай, когда из каждой точки рассматриваемой области оба заряда начинают двигаться с начальным расстоянием между ними г3=0.3. Максимальное значение показателя Ляпунова Атах=3.850257, минимальное значение показателя Ляпунова Атіп=0.905390
Для случая, когда один заряд начинает свое движение из одной и той же фиксированной точки (например, с координатами (1;1)), а другой имеет в качестве начальных координат последовательно координаты каждой точки исследуемой области, карта динамических режимов представлена на рис. 4:
-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5
2.5
2.0
1.5
1.0
0.5
-0.0
-0.5
-1.0
-1.5
-2.0
-2.5
Рис. 4. Карта показателей Ляпунова в плоскости координат (х,у). Максимальное значение показателя Ляпунова Атах=11.170605, минимальное значение показателя Ляпунова Атіп=-0.760053
а) б)
Рис. 5. Типичные траектории зарядов: а - для области хаотических режимов У0Х1 = У0Х2 = У0у = У0у2 = 0, х01=0,7, у01=0,4;
б - для области устойчивого поведения, где У0Х1 = У0Х2 = У0 у1 = У0 у2 = 0, Х01=0,6, у01=-0,1.
Светлая линия соответствует траектории движения первого заряда (первоначальная и смещенная), темная - второго
Типичный вид хаотических (рис. 5, а) и регулярных (рис. 5, б) траекторий в плоскости координат для случая взаимодействующих зарядов представлен на рис. 5. На рис. 5, а видно, что в хаотическом режиме траектории зарядов, стартующих при практически совпадающих начальных условиях, с течением времени быстро расходятся, для них показатель Ляпунова положителен.
Пример регулярных траекторий представлен на рис. 5, б, где две траектории, имеющие ту же незначительную разницу в начальных условиях, по существу, совпадают. Для них показатель Ляпунова меньше 0 или близок к 0.
Таким образом, появление второго заряда, взаимодействующего с первым, существенно меняет карту динамических режимов. Практически вся исследуемая область (рис. 3) соответствует хаотическим траекториям. Вместе с тем, если начальное положение одного из зарядов фиксировано в одной и той же точке (рис. 4), то появляется значительная часть области, соответствующая регулярным траекториям (белый или светло-серый цвета на рис. 4). Физически это понятно, так как заряд, стартующий в этой части области, еще слабо взаимодействует с другим, ввиду относительно большого расстояния между ними.
ЛИТЕРАТУРА
1. Байбурин В.Б., Мантуров А.О., Юдин А.В. Хаотическое поведение зарядов в скрещенных полях // Известия вузов. Прикладная нелинейная динамика. 2002. Т.10. № 6. С.62-70.
2. Поттер Д. Вычислительные методы в физике. М.: Мир, 1975. 320 с.
Байбурин Вил Бариевич -
доктор физико-математических наук, профессор,
Заслуженный деятель науки РФ, действительный член РАЕН, заведующий кафедрой «Программное обеспечение вычислительной техники и автоматизированных систем» Саратовского государственного технического университета
Трачук Светлана Анатольевна -
студентка Саратовского государственного технического университета