Материал: Sb97952

При записи вектора D(t, k) нужно следить за тем, чтобы неравенство v c не нарушалось ни в один момент рассматриваемого времени.

4. Для вычислений используйте функцию rkfixed.

Вычисление решения в 1000 промежуточных точках на отрезке времени

[0, 1·10–13 с]: RK = rkfixed (k, 0, 10–13, 1000, Z).

5. Выведите таблицу с результатами моделирования движения заряженной частицы в электромагнитном поле. Матрица, получаемая в результате решения (табл. 5.1), должна содержать семь столбцов: первый – t, второй – x(t), третий – vx, четвертый – y(t), пятый – vy, шестой – z(t), седьмой – vz.

Таблица 5.1

Результаты моделирования движения заряженной частицы в магнитном поле

26

Траектория частицы изменяется в очень узком пределе, что является следствием малого значения магнитной индукции. Действие электрической составляющей можно наблюдать по изменению составляющих скоростей vy и vz, вследствие чего частица движется по винтовой траектории.

6. Постройте график зависимости z f (x) для всего временного интервала моделирования (рис. 5.1).

z, м

0..0604

RK<5>

00..0602

Рис. 5.1

Для получения доступа к столбцу, задающему координаты, в массиве [RK] необходимо использовать треугольные скобки , в которые вписывается номер столбца.

5.5.Содержание отчета

1.Основные теоретические сведения по моделированию движения заряженной частицы в электромагнитном поле для заданных условий.

2.Привести таблицу с результатами моделирования движения заряженной частицы в электромагнитном поле для заданных условий.

3.Привести график зависимости z f (x) для всего временного интер-

вала моделирования.

4.Привести листинг кода программы, моделирующей движение заряженной частицы в электромагнитном поле для заданных условий.

5.Сформулировать выводы по результатам моделирования движения заряженной частицы в электромагнитном поле для заданных условий.

5.6. Вопросы для самоконтроля

1. По какой траектории в общем случае движется точечный электрический заряд в постоянном (не меняется с течением времени) однородном

27

(одинаковом в любой точке геометрического пространства) электромагнитном поле?

2. Меняется ли кинетическая энергия (энергия движения) заряженной частицы, движущейся в постоянном или однородном электромагнитном поле?

6. МОДЕЛИРОВАНИЕ ПРОЦЕССА РАССЕЯНИЯ ДВИЖУЩИХСЯ ЗАРЯДОВ

Цель работы: ознакомление с методикой моделирования рассеяния движущихся зарядов друг на друге.

6.1. Основные теоретические положения

Движущийся электрический заряд создает вокруг себя электрическое и магнитное поля, зависящие от скорости его движения:

B μ0 q v r , 4π r3

где v/с – безразмерная скорость электрического заряда; – угол между направлением вектора скорости движения заряженной частицы и ее радиу- сом-вектором.

Уравнение движения одной заряженной частицы в электромагнитном поле другой заряженной частицы:

Решив (6.2) относительно радиуса-вектора, можно получить траекторию движения заряженной частицы в электромагнитном поле другой движущейся заряженной частицы.

В случае если скорость движения заряженной частицы много меньше скорости распространения электромагнитной волны в вакууме ( v c ), то релятивистский фактор принимается равным единице ( 1, 0 ). При этом необходимо следить за тем, чтобы неравенство v c не нарушалось ни в один последующий момент рассматриваемого времени.

6.2.Указания к выполнению работы

1.Смоделировать (любым методом, описание выбранного метода привести в основных теоретических положениях) движение заряженной частицы в

28

электромагнитном поле другой движущейся заряженной частицы для заданных условий.

2.Получить таблицу с результатами моделирования движения заряженной частицы в электромагнитном поле другой движущейся заряженной частицы для заданных условий.

3.Построить график зависимости z f (x) для всего временного интер-

вала моделирования.

6.3. Пример типового задания

Смоделировать рассеяние двух движущихся электрических зарядов, изначально находящихся в состоянии покоя. Заряды: q1 = N пКл, q2 = (6 – N) пКл, где N – номер студента в списке группы. Массы зарядов:

y01 = (8 – N) м; z01 = 0; x02 = (13 – N) м; y01 = (–8 + N) м; z02 = 0.

6.4. Пример решения задания

Решение можно выполнить в ПО MS Excel. Время разбивается на отрезки ( dt 2E10 c), в течение которых скорость и ускорение считались постоянными. Исходные координаты заряженных частиц:

x01 = –11 · 10–6 м; y01 = 4 · 10–6 м; z01 = 0; x02 = 11 · 10–6 м; y02 = –4 · 10–6 м; z02 = 0.

Заряды заряженных частиц:

q1 = 2 · 10–31 Кл, q2 = 4 · 10–31 Кл.

Массы заряженных частиц:

m1 = 1.1 · 10–11 кг, m2 = 2· 10–12 кг.

Начальные скорости равны нулю: vx = 0; vy = 0; vz = 0.

Начальные координаты по оси z равны нулю, частицы не перемещаются по данной оси, поэтому задача псевдотрехмерная. Кроме того, при рассеивании частицы не успевают достичь скоростей, при которых необходимо учитывать магнитное поле и релятивистский эффект, поэтому ими можно пренебречь. Тогда для каждой частицы можно записать уравнения:

29

Для получения рассеяния пришлось предположить, что частицы имеют начальные скорости. Сначала вычисляется ускорение первой частицы с учетом скоростей и координат двух частиц в предыдущий отрезок времени (6.7), затем ускорение второй частицы:

v r (v x1 v x2 )2 (x1 x2 ) (v y1 v y2 )2 ( y1 y2 ) , v (v x1 v x2 )2 (v y1 v y2 )2 ,

r (x1 x2 )2 ( y1 y2 )2 ,

Далее последовательно вычисляют скорости (6.5) и координаты (6.4) (табл. 6.1). После этого данные выводятся на график в виде двух зависимостей y f (x) .

30