Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Елютин Експерименталная физика.Лабораторный практикум 2011

.pdf
Скачиваний:
11
Добавлен:
12.11.2022
Размер:
5.19 Mб
Скачать

Схема установки изображена на рис. 34.8. Излучение источника 1 – галогенной лампы накаливания, проходит через систему нейтральных фильтров и пленку, укрепленных в держателях, и через входную щель попадает в монохроматор 2. Из входной щели монохроматора излучение направляется в блок приемника 3, основой которого служит фотоэлемент. Выходной сигнал с фотоэлемента после усиления через кабель и разъем подается на гнездо блока питания 4. Источник излучения 1 также подключается к блоку питания 4. Регистрация интенсивности излучения в относительных единицах осуществляется на шкале цифрового вольтметра 5. Вход цифрового вольтметра соединен с гнездами блока питания 4. Необходимые длины волн устанавливаются рукояткой 2а монохроматора и отсчитываются по цифровому механическому счетчику с точностью 0,2 нм.

Рис. 34.8

Чувствительность установки регулируется диафрагмой 3а и ступенчатым переключателем 3б. Рукоятка 3в устанавливает «нуль» вольтметра. На лицевой панели блока питания 4 имеется тумблер 4а – «сеть» и сигнальная лампа, тумблер питания фотоэлемента 4б «накал-недокал», обеспечивающий работу источника излучения в

31

двух режимах, и гнезда для подключения вольтметра 4в. Рабочий предел вольтметра – 20 В.

Снятие спектров пропускания осуществляется в режиме «недокал» источника излучения 1 в области 450-630 нм, которая является областью наибольшей чувствительности установки. В этом диапазоне длин волн используется входная щель монохроматора шириной 0,25 мм, а выходная щель шириной 0,05 мм. Во избежание перегрузки вольтметра, световой поток ослабляется с помощью двух шестикратных нейтральных фильтров. В этом случае положение ступенчатого переключателя чувствительности 3б – «3».

До начала измерений необходимо установить «0» на табло вольтметра с помощью рукоятки 3в, предварительно перекрыв доступ света в монохроматор. Следует помнить, что с вольтметра производится отсчет установившегося значения напряжения (для чего требуется выждать 10 – 15 секунд).

Измерения производят вначале без пленки с тем, чтобы снять аппаратную функцию спектральной чувствительности фотоэлемента U0 ( ). При перегрузке вольтметра из-за большой величины све-

тового потока необходимо ввести один, а при необходимости и два светофильтра. Коэффициент ослабления света k указан на самом светофильтре. Светофильтры не убираются до конца измерения.

После этого производят измерение спектра пленки As2S3 толщиной 1,8 мкм U ( ) в том же диапазоне длин волн. На коротковолновом участке диапазона показания вольтметра могут стать малыми, что потребует снятия фильтра, с тем, чтобы восстановить значащие цифры в разрядах дисплея вольтметра.

Коэффициент пропускания вычисляется по формуле

 

D( ) U ( ) U0 ( ).

(34.6)

Если в ходе выполнения работы снимался фильтр, то необходимо вводить поправочный коэффициент, т.е. использовать формулу

D( ) kU ( ) U0 ( ).

(34.7)

Порядок проведения эксперимента

1. Включите блок питания монохроматора, вольтметр и источник света. Дайте приборам прогреться не менее 15 минут.

32

2.Установите на счетчике монохроматора длину волны 600 нм. Произведите установку «0» вольтметра, перекрыв предварительно доступ света в монохроматор.

3.Измерьте без пленки аппаратную функцию U0 ( ) в диапазоне

длин волн 630-450 нм через каждые 10 нм. Результат занесите в таблицу 34.2.

4. Произведите с пленкой измерения спектральной функции U ( ) в том же диапазоне длин волн. Отметьте измерения, где был снят светофильтр. Результат занесите в таблицу 34.2.

Таблица 34.2

, , нм

U0 ( )

U ( )

D( )

 

 

 

 

630

 

 

 

 

 

 

 

620

 

 

 

 

 

 

 

 

 

 

 

460

 

 

 

 

 

 

 

450

 

 

 

 

 

 

 

5.Вычислите D( ) по формулам (34.6) и (34.7).

6.Постройте график D( ) и по графику определите ГР и .

7.Вычислите n по формуле (34.5) и Eg формуле (34.2). Срав-

ните полученный результат с табличным.

Контрольные вопросы

1.Как образуются энергетические зоны в кристалле?

2.Объясните «прозрачность» веществ с точки зрения зонной теории.

3.Объясните физический смысл ГР.

4.Что происходит с запрещенной зоной кристаллического полупроводника при появлении дефектов в его структуре?

33

Литература

1. Яворский Б.М., Пинский А.А. Основы физики в 2-х томах.

Т2. М.: Физматлит. 2003.

2.Савельев И.В. Курс общей физики. Т 3. М.: Наука. 1989. С. 197-210.

3.Физические величины: Справочник. / Под ред. И. С. Григорьева, Е.З. Мейлихова. М.: Энергоатомиздат. 1991.

РАБОТА 35

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ТРАЕКТОРИЙ ЗАРЯЖЕННЫХ ЧАСТИЦ В ОДНОРОДНЫХ ПОЛЯХ

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

Введение

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

Движение частицы в однородном электрическом поле с учетом сопротивления среды

При моделировании движения положительно заряженного тела (частицы) в однородном электрическом поле (для определенности поле направлено вертикально вниз) сила сопротивления среды, в которой движется частица, была выбрана пропорциональной квадрату скорости частицы. Расчеты показали сильное влияние сопротивления на форму траектории. Модель позволяет определить время движения частицы в заданной области пространства, мгновен-

34

ную скорость, геометрические параметры траектории, величину диссипации энергии при движении в среде.

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

L

v2 sin2

,

 

 

0

 

 

 

 

qE m

 

 

 

 

 

с максимальным значением

 

L

 

 

v2

 

,

(35.1)

 

0

 

 

 

 

МАКС

 

qE m

 

 

 

 

 

где v0

– скорость влета частицы,

– угол влета по отношению к

эквипотенциальной поверхности. Максимальное смещение частицы против направления поля составит

H v02 sin2 . 2 qEm

Время движения

T

 

 

2v0 sin

,

 

 

 

 

 

ПОЛЕТА

 

 

 

qE m

 

 

 

 

 

 

 

время движения при максимальной дальности:

 

 

 

 

 

 

v0

.

 

T

 

 

 

 

2

(35.2)

 

 

 

 

 

МАКС

 

 

qE m

 

 

 

 

 

Траектория движения частицы без учета сил сопротивления выглядит следующим образом:

y

y

0

x tg

x2 (qE m)

.

(35.3)

2v2 cos2

 

 

 

 

 

 

 

 

 

 

0

 

 

 

Уравнение второго закона Ньютона для движения с учетом силы сопротивления задается векторным равенством

ma qE F

.

(35.4)

СOПР

 

 

Предполагается, что сила сопротивления FСОПР направлена по касательной к траектории против движения тела и может быть выбрана либо пропорциональной вектору скорости тела с обратным знаком, т.е. в виде

35

FСОПР rЛИНv,

либо пропорциональной квадрату модуля скорости тела в виде

 

 

v

 

F

r

v2.

 

СОПР

КВ v

В зависимости от выбора, вектор силы сопротивления имеет следующие проекции:

FСОПР rЛИНvx , rЛИНvy ,

или

FСОПР rКВvx v2x v2y , rКВvy vx2 vy2 , ,

где rЛИН, КВ – соответствующий коэффициент сопротивления среды,

в которой движется тело.

Уравнение (35.4) в проекциях на оси координат приводит к системе дифференциальных уравнений по времени для координат x и y

d 2 x

 

 

 

r

 

 

 

 

 

 

 

 

ЛИН

 

v

x

 

 

dt2

 

 

 

 

 

 

 

m

 

 

d 2 y

 

 

 

r

 

 

 

qE

 

 

 

 

ЛИН

v

y

 

 

dt2

 

 

m

 

 

 

m

 

 

 

 

r

 

 

 

 

 

 

 

vx2 vy2

 

 

 

 

 

КВ

vx

 

 

,

 

 

 

 

 

 

 

m

 

 

 

 

(35.5)

 

 

r

 

 

 

 

qE

 

vx2 v2y

 

 

 

 

 

КВ

vy

 

 

 

.

 

 

 

 

m

 

 

 

 

m

Кинетическая энергия частицы равна

K

mv2

(35.6)

 

.

2

 

 

 

Потенциальная энергия заряженной частицы:

 

qE( y y0 ).

(35.7)

Для удобства дальнейшего рассмотрения перейдем к безразмерным переменным:

Y

( y y0 )

, X

(x x0 )

,

t

,

 

 

 

 

 

 

 

 

 

LМАКС

TМАКС

 

 

 

LМАКС

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

v

 

 

,V

 

 

 

 

 

 

 

vy

 

 

 

 

 

 

 

 

(35.8)

2

x

 

y

 

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

v0

 

 

 

 

 

 

 

 

 

 

v0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

v2

 

 

 

H

1

 

b

r

 

 

 

 

 

2

0

, b

r

, h

 

sin2 .

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

qE

 

 

L

 

2

ЛИН

 

ЛИН

 

 

 

 

 

КВ

 

КВ qE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

МАКС

 

 

36

В обозначениях (35.8) уравнение (35.3) для движения без учета силы сопротивления среды примет следующий вид:

Y Y X tg

X 2

 

.

(35.9)

2cos2

 

0

 

 

Уравнения Ньютона (35.5) при учете силы сопротивления среды можно записать в безразмерных переменных (35.8) следующим образом:

dVx

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

bЛИНVx

bКВVx Vx

Vy

,

d

dVy

 

2 bКВVy

 

 

,

2 bЛИНVy

Vx2 Vy2

 

 

d

 

 

 

 

 

 

 

 

dX

Vx ,

 

 

 

 

(35.10)

 

 

 

 

 

d

 

 

 

 

 

 

 

 

dY

d

Vy .

Выражения для потенциальной и кинетической энергии частицы в обозначениях (35.8) будут выглядеть как

Wp

 

 

2Y ,

 

 

(35.11)

mv2

2

 

 

 

0

 

 

Vx2

Vy2

 

 

W

 

K

 

 

.

(35.12)

mv2

2

 

2

k

 

 

 

 

 

0

 

 

 

 

 

 

Потери энергии на работу против силы сопротивления среды составят

W 1 2Y

Vx2 Vy2

.

(35.13)

 

2

 

 

Систему дифференциальных уравнений (35.10) представим в виде уравнений в конечных разностях, когда дифференциал функции заменяется разностью ее последовательных значений в узлах (n 1) и (n) расчетной сетки, разделенных временным шагом h. Таким образом, система уравнений в конечных разностях для движения заряженной частицы в однородном электрическом поле будет иметь следующий вид:

37

 

Vx n 1

Vx n

 

 

 

V n

 

b V n

 

 

 

 

 

 

 

 

 

b

 

V n 2 V n 2

,

 

 

h

 

 

 

 

ЛИН x

 

КВ x

 

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vy n 1

Vy n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

n

2 b V

V n 2 V

n 2

,

 

 

 

 

2 bЛИНVy

 

КВ

y

x

y

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X n 1

X n

 

Vx n ,

Y n 1

Y n

Vy n ,

 

n 1 n h,

 

(35.14)

 

 

h

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2 n 2

W n 1 2Y n Vx Vy . 2

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

X (0) 0, Y (0) 0, (0) 0,Vx(0)

2 cos ,Vy(0)

2 sin .

Точность такой численной процедуры (схемы Эйлера) составляет o(h).

Сравнение численного решения с табуляцией классической формулы (35.9) показывает влияние силы сопротивления среды. Заметим, что координата X ( ) в уравнении (35.10) получается как решение соответствующего дифференциального уравнения. Эта величина не может быть заменена известным выражением

X ( ) 2 cos из-за наличия ускорения вдоль оси X , как результата действия силы сопротивления. Результат одного из расчетов (b 1, 45 ) представлен на рисунке 35.1. График 1 – без учета сопротивления среды, график 2 – с учетом сопротивления среды, график 3 – потери энергии на сопротивлении среды.

Моделирование пузырьковой камеры

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

38

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

Рис. 35.1

Прохождение заряженной частицы через перегретую жидкость приводит к образованию вдоль траектории частицы «зародышевых» центров кипения. За время порядка 0,5-3 мс пузырьки, образующиеся на зародышах, достигают размеров 50-300 мкм и могут быть сфотографированы при освещении их импульсным источником света.

Во время движения в пузырьковой камере на частицу действуют сила Лоренца q[v B] и сила сопротивления среды, выбранная в на-

шей задаче в виде r v v2 (здесь v – вектор скорости частицы,

КВ v

q – ее заряд, B – индукция магнитного поля, rКВ – коэффициент сопротивления среды). Следовательно, второй закон Ньютона для частицы можно записать в виде

 

 

v

 

 

ma

q[v B] r

v2.

(35.15)

 

 

КВ v

 

Будем рассматривать случай, когда частица влетает в камеру со скоростью v0 перпендикулярно к направлению вектора индукции

магнитного поля B. В этом случае движение частицы будет происходить в плоскости, перпендикулярной к векторуB. Движение частицы удобно рассматривать в системе координат, у которой ось Z

39

совпадает с направлением вектора B, а ось Y составляет с начальным направлением вектора скорости v0 угол (рис. 35.2).

Рис. 35.2

Запишем уравнение (35.13) в проекции на оси такой системы координат:

 

d 2 x

 

 

dy

 

 

 

 

 

 

 

 

 

m

 

qB

r v

 

 

 

v2 v2 ,

 

 

dt2

 

 

 

 

 

 

 

 

 

 

dt

КВ

x

 

x

y

 

 

 

d 2 y

 

 

dx

 

 

 

 

 

 

 

 

 

m

 

qB

r

 

v

 

 

v2

v2

,

(35.16)

dt2

 

dt

 

 

 

 

 

 

КВ

 

y

 

x

 

y

 

 

m

d 2 z

 

0.

 

 

 

 

 

 

 

 

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Начальные условия для задачи (35.16) выбраны следующими:

x(0) x0 , y(0) y0 , vx

v0 sin , vy

v0 cos .

Перейдя к безразмерным переменным в первых двух уравнениях системы (35.16), можно получить следующую систему:

 

dX

 

 

dVx

 

2 Vy bКВVx

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

Vx ,

 

 

Vx

Vy

 

 

 

 

 

 

 

 

 

d

d

 

 

 

(35.17)

 

 

dY

 

 

dVy

 

2 Vx bКВVy

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

Vy ,

 

 

 

 

Vx

Vy

 

 

 

 

 

 

 

 

 

d

 

 

d

 

 

 

 

 

 

 

vx, y

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

r v2

 

r v2

 

 

 

где введены переменные

 

 

, b

 

 

 

КВ 0

 

КВ 0

, V

x,y

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

КВ

 

2 mv0 T

 

qBv0

 

v0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]