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

книги из ГПНТБ / Лавренко, В. А. Рекомбинация атомов водорода на поверхностях твердых тел

.pdf
Скачиваний:
27
Добавлен:
22.10.2023
Размер:
12.65 Mб
Скачать

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

C m s h a J = O m + B V i n

и с помощью (9.21) получаем вместо последнего две однопараметрические системы уравнений

 

 

Ф„ =

-

-

В) рЕ„

 

-

 

X

 

У,

SmnOm;

 

 

 

. .

 

г,-.

 

 

th а.,,1

. th <xnl V

с

 

th ап1

(9.23)

Ф« =

 

 

 

 

 

 

В) рчп

— f

X - —

> ^ 5 „ т Ф

т

 

- 2 - Ф„;

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

.

N

 

 

 

 

 

 

,

/V

 

 

 

 

 

Є

 

V і

 

° т

 

9

• ТІ

 

L V

а , п

ih п /<?

 

 

~

р

2 J cha,,,/0 '""'

 

р

Z J а т ch «„,/

 

a '"t o mn-

 

 

 

 

i7i=I

 

 

 

 

 

 

т = 1

 

 

 

 

 

Далее

воспользуемся

методом

последовательных

приближений

по (Л — В) р =

Я,. При не слишком

больших

р эта величина

оста­

ется меньше единицы даже в

случае

полного

поглощения атомов

образцом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф„ =

2 ^ ' + 1 Ф « ;

Ф« =

2

k + ' q i

 

 

(9.24)

 

 

 

 

 

 

/=0

 

 

 

 

1=0

 

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фп =

-

1

^

- ;

^п

=

^

( -

^ -

Ф п

)

;

(9.25)

В результате функция /г (г/) (9.18) может быть записана в виде

/г (х) = Q (х) + 2

(х) — Bp (х) +

В 2 Я'+'а/ (х) = hx (х) +

 

 

/=0

 

;=0

 

 

 

+ В/г2 (х).

(9.26)

Задаваясь

однажды

выбранной совокупностью параметров в , р

и точек х£ (i

=

1, 2

р), в которые будет помещен зонд, все вхо­

дящие в (9.26)

функции можно вычислить заранее с помощью (9.22—

9.25) на ЭВМ

 

 

 

 

 

 

 

(9.27)

со/ (х) = > , Ф ^ ж ^ Г ; СТ/(х)

= ^ ф / 1 Ж ^ г .

Подлежащая определению величина (А — В) (или X) отыскива-

ется по следующей схеме. Температура зонда является линейной функцией h

 

 

T(x)

=

K[hl(x)

+ Bh2

(х)].

 

Располагая

соответствующим

набором

экспериментальных то­

чек Tj

(х — х{),

применяем

метод

наименьших квадратов,

миними­

зируя

величину

 

 

 

T(xt)

 

 

 

 

У ( К , В А ) = 2

1

(9.28)

 

 

Тс

 

 

 

 

 

 

 

Из ~ШГ ~ ^'

dJ

 

У ч а е м

К = К

(X), В = В (I) и затем

= 0

П 0 Л

путем прямого машинного подбора единственного параметра А,

находим минимум

(9.28). Подбор удобно начинать с Х0,

получаемого

с помощью (9.26) в первом

приближении

по К и для В =

0. В этом

случае целесообразно

использовать

условие

 

=

0

и

в

резуль­

тате

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

ш0

іхї)

 

 

 

 

X

 

 

 

Тс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Li

ТІ

 

 

 

 

 

 

1=1

 

 

 

 

 

 

i'=i

 

 

о /

^

 

 

 

 

 

 

 

со0

(хе)

Щ (х,.) О (*,-,)

р

Р

«о (•*;<)

 

 

х

у

A (*t)

^

 

 

J

r<

Z J

 

ТІ

Li

Tc

 

 

rjf,

 

 

 

 

 

 

 

 

 

 

 

= 1

 

i'=l

 

'

t=i

 

 

i'=I

 

«

 

 

 

Далее для У А) осуществляется перебор

в пределах

от 0,254 %0

до 3,453 Я,0, причем каждое последующее значение X на

10% отли­

чается от предыдущего.Всего26 точек: Xt

l , l ' , %0(t = 0,

1 , 1 3 ) ,

X's 0,9s , X0 (S =

1,2,

 

 

13). При этом число членов в суммах по

/ (9.26) выбирается с запасом, исходя из значения

Я0 . Во всех слу­

чаях оказывается

вполне

достаточным взять

/ м а К с = 1 0 ,

и

точка

минимума

(9.28)

попадает

в указанный

промежуток.

 

 

 

Далее

рассмотрим

разработанный

метод с точки

зрения

приня­

той нами методики исследования температурной зависимости коэф­ фициента гетерогенной рекомбинации атомов. Удобно провести некоторые сравнения с тепловой методикой, описанной в § 7.

Предложенный там метод предусматривал регистрацию

нагре­

вания (за счет рекомбинации) самого исследуемого

образца

— дис­

ка, перекрывающего сечение реакционной трубки и

перемещаемого

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

Изготовление ж е образцов

значительного диаметра может

быть

связано с технологическими

трудностями.

 

Подобный недостаток становится существенным при темпера­

турных измерениях, где приходится измерять малые прибавки

бла­

годаря реакции, идущей на предварительно нагретом до высокой температуры образце.

Не менее серьезными оказываются трудности осуществления эф­ фективного подогрева подвижного образца при условии сохранения

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

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

в отношении рекомбинации тер-

21 22 п /3 /4

 

 

.иозонд,

что

позволяет

получить

 

 

 

 

 

 

информацию о возмущении

рас­

 

 

 

пределения

атомов

газа

образ­

 

 

 

цом .

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате

решения

мате­

 

 

 

матической

задачи

(9.5—9.26) с

 

 

 

помощью ЭВМ

практически

по­

 

 

 

лучена

функция

q

(х,

 

у,

у'),

 

 

 

определяющая

тепловой

поток,

 

 

 

входящий в термозонд, на раз­

 

 

 

личных расстояниях х от источ­

 

 

 

ника постоянной

 

интенсивности

 

 

 

в зависимости

от

коэффициентов

 

 

 

рекомбинации

образца

у

и

шай­

'Рис. 15. Экспериментальная

установка

бы у'- Д л я

определения

 

соответ­

для исследования температурной зави­

ствующих тепловых

потоков за

симости коэффициента

рекомбинации.

счет рекомбинации

q

(х)

 

измеря­

Т (х). Блок-схема

 

 

ют значения

температуры

зонда

установки представлена

на

рис.

15.

 

 

 

 

 

В диффузионной кварцевой трубке 12 диаметром 54 мм на рас­

стоянии / = 40 мм от разрядной трубки /

закреплена

перекрываю­

щая кварцевая шайба 13. Образец — диск

18 диаметром р — 12

мм

(толщина 0,1—0,3 мм)

помещен

во внутреннее

отверстие

 

шайбы.

Нагрев образца до заданных температур в интервале от комнатной до 1200° С осуществляется с помощью печи 20 специально подо­ бранной формы. Электропечь изготовлена из листового вольфрама

толщиной 0,5

мм электроискровой обработкой и крепится к образ­

цу с тыльной

стороны. При исследовании образцов из проводящих

материалов между печью и образцом помещена тонкая изолирую­ щая прокладка из кварца. Температура тыльной стороны образца контролируется вольфрам-рениевой термопарой 21 диаметром 0,1 мм. При необходимости температуру поверхности образца со стороны реакционной камеры проверяют подвижной термопарой-зондом 17;

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

Строго соосно с 12 вварена боковая трубка 7. В этой трубке смонтирована система передвижения 9 зонда с помощью электро­ магнита 8. Зонд — тщательно градуированная термопара, в точке сварки — шарик диаметром ~0,3 мм. Применение столь неболь­ шой подвижной пробы-зонда приводит к значительному сокраще­ нию времени установления теплового равновесия в процессе изме­ рений. Отметим, что в качестве зонда можно использовать не только термопару, но и термосопротивление, проволочную петлю из сплава

с высоким удельным сопротивлением. Проволоки термопар

зонда

находятся внутри

кварцевой

трубки 10 диаметром 3 мм,

проходя­

щей через разрядную трубку

/ . Д л я более строгой центровки

зонда

по оси реакционной

камеры

на входе в последнюю помещена

квар­

цевая центрирующая

шайба

15 с внешним диаметром 6 мм и внут­

ренним — 3 мм,

крепящаяся тремя тонкими кварцевыми

подвеса­

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

При заданной температуре поверхности образца в результате измерений получаем зависимость температуры зонда от х при вклю­ ченном (Т (х)) и выключенном (V (х)) разряде, т. е. при наличии и отсутствии атомов на образце и зонде. Положение зонда относитель­ но образца контролируется с точностью 0,1 мм. Соответствующие значения теплового потока q (х) дают возможность определить па­ раметры у и у'.

Уравнение теплового баланса для зонда можно записать в виде

q (х) +

а,7 І +

о' Т І +

a, (T x -

Tx) +

ах

( Г 2

-

Тх) + а°хТІ =

 

 

 

 

 

 

=

аТх + а°хк0).

 

 

 

 

(9.29)

Здесь

q (х) — тепло,

поступающее

за счет

рекомбинации

атомов на

зонде;

Ти

Т2,

Т0

и

Тх

— температура

образца,

шайбы

(средняя),

стенок

трубки

и зонда

соответственно;

ах,

ах,

а°х — не

зависящие

от температуры параметры, ответственные за теплопроводность

и

конвекцию между зондом и

соответственно

образцом, шайбой

и

стенками трубки; ах,

а[х, а°х, а

— параметры,

характеризующие теп­

лопередачу за счет

излучения образцом, шайбой, стенками трубки

и зондом.

 

 

 

 

Нелинейными по температуре членами, характеризующими теп­ ловой поток за счет конвекции (ограничиваемся линейным прибли­ жением), мы пренебрегаем.

Из

условия теплового равновесия

при 7 \ =

TQ,

q (х)

= 0, Г 2 =

= Тх

— Т0 получаем

ах + ах + о°х = ст.

 

 

 

Подставляем это значение а в (9.29).

 

 

 

Тогда

 

 

 

 

 

 

 

 

q (х) + ах (71 -

Ті) + ох(Т$

~

П) + ахх

-

Тх)

+

 

+ а\. (Тг -

Тх) = ст° х

-

Tt) + <4 х

-

Т0).

Если не пренебрегать излучением шайбы, то для нахождения q (х) требуется решить систему из семи уравнений с семью неизвестны­ ми: q (х),ох, ох, ах, <хх, а° и ах. Но далее будем пренебрегать излу­ чением шайбы ввиду поправочного характера этой величины и малой излучательной способности кварца. Тогда уравнение (9.29) примет вид

q (х) + ихТ\

+

а,

(7\

-

Тх)

+

а,

2

-

Тх)

+

а°хТ*0 =

оТ\

+

 

 

 

 

 

+

ax(Tx-TJ.

 

 

 

 

 

 

 

(9.30)

При

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т1

= Т2

= Тх

= Т0

и

q (х) = 0

 

 

 

 

 

 

ах

,

 

о

 

о

= о — ах.

 

 

 

 

 

 

+ ах = а;

ах

 

 

 

Подставим

значение

а°х в (9.30). Получаем

 

 

 

 

 

q (х) + ох

(71 -

Го) +

ах

(7\ -

Тх)

+

ах

2

-

Тх)

= а (Т\

-

П) +

 

 

 

 

 

+ < 4 ( 7 \ - Т 0 ) .

 

 

 

 

 

(9.31)

Поскольку о не зависит от х (положение зонда), разделим все члены (9.31) на ст, но оставим прежние обозначения — поскольку мультипликативная константа в q (х) для расчета несущественна. Тогда

q (х) + ах (71 -

71)

+

ах

(7\ -

Тх) + ах 2

-

 

Тх)

=

 

 

 

=

Г х

- П

+

а°(Тх0).

 

 

 

 

 

(9.32)

Кроме q (х), теперь осталось четыре

параметра:

 

ах,

ах,

ах,

а°х,

т. е. для вычисления

q (х) надо составить систему из (9.32) и четырех

уравнений при выключенном

разряде для четырех

различных

зна­

чений температуры

образца

7\.

 

 

 

 

Тх — Т0

 

 

 

 

 

 

Пренебрегаем далее

нелинейным

по

членом

в

правой

части (9.32), так как даже

при

Тх

=

600° К он

составляет

всего

около 5% линейного члена.

 

 

 

 

 

 

 

 

 

 

 

 

 

Такое заключение сделано на основании работы [219]. Линейное

приближение выполняется,

 

когда

 

 

 

 

 

 

 

 

 

 

 

1

Т,,)

х

Г 0

) 2

+

4 х

— Тв) Т0

+ &ТІ

1.

(9.33)

- і - л р а а ( Т х -

 

 

 

 

\

 

*

°'

 

 

0

«

В данном случае

радиус

зонда

р =

0,15 см;

а

 

=

0,5—

коэффи­

циент излучательной способности материала зонда (степень черно­

ты),

а = 5,67

• Ю - 1 2

до/с/см2 • сек • град. — универсальная постоян­

ная

Стефана,

X — коэффициент теплопроводности водорода (Я =

= 1,8 • Ю - 3 ) ;

g=

1,111 — значение

коэффициента «перекрыва­

ния»

для шарика в

неограниченном

пространстве.

Тогда

при Тх

=

600° К

величина

(9.33)

составляет

~ 5 % от

линейного

члена.

 

 

 

 

 

Теперь

(9.32) можно записать в виде

 

 

q (х) + их

(71 -

П)

+ а, х

- Тх) + ах

2 -

Тх) = а"х х

- 70),

 

 

 

 

 

 

 

(9.34)

где ах — новый параметр вместо а°х.

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

q(х) + ах(Т\~П)

+

ах(Т, - Тх) = Р, х-

Т0).

(9.35)

Д л я определения q (х)

по

экспериментальным

данным

теперь

должна решаться на ЭВМ система из четырех уравнений с четырь­ мя неизвестными. Ее можно свести к системе из трех уравнений с

тремя неизвестными, так как ах

(излучение

образца)

можно

вычис­

лить — с точностью

до

коэффициента

— заранее, предполагая за­

кон

Ламберта.

 

 

 

 

 

 

 

 

 

 

 

 

 

где

с — константа.

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение

(9.35) можно

 

записать

 

 

 

 

 

 

 

 

q(x)

+ cfl~

/ ( ^ ~ " + р

 

2

(71

-71)

+

ах(7\

- Т х ) =

р , ( Т х - Т 0 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.37)

 

Разделим

(9.37)

на

 

постоянную

с

и

снова оставим прежние

обозначения.

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

<7 (к) +

(1 -

/ ( / ^ ~ * +

р

2 )

(ТІ ~ П)

+

ах

г

-

Тх)

=

р, х-

Тп).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.38)

Уравнение

(9.38) и еще два уравнения

вида

 

 

 

 

 

(1 -

/ ( /

_ ^ ~ *

) V і

-

ЇЇ)

+ ах (Т -

Гх)

=

р., х

-

Т0)

(9.39)

при температурах образца Т = Тг и Т = Т\ (соответствующие из­ мерения Тх при выключенном разряде должны быть проведены)

дают возможность

исключить ах и $х и вычислить

q (х) из (9.38).

Подсчитаем величину коэффициента £ л =

1 —

У ( / _ х ) * + р З }

для опытов, проведенных при диаметре образца

12,0 мм и расстоя­

нии от источника

атомов до образца / =40,0 мм; х — расстояние от

разрядной трубки

до зонда.

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

6

 

 

 

 

 

.V, MM

 

 

 

 

Л", MM

1 X

 

 

Ъ

 

 

 

 

 

V(l- л-)2 -f

P=

 

V U - л*)2 +

р"

ъх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

 

0,976

 

0,024

32

 

0,800

 

0,200

 

 

 

 

 

20

 

0,958

 

0,042

34

 

0,707

 

0,293

 

 

 

 

 

25

 

0,929

 

0,071

35

 

0,640

 

0,360

 

 

 

 

 

28

 

0.894

 

0,106

36

 

0,554

 

0,446

 

 

 

 

 

30

 

0,858

 

0,142

 

 

 

 

 

 

 

 

 

Если

считать,

что р\. =

const

(не зависит от х),

 

то, поделив

в

(9.38)

и (9.39)

все члены

на Р и вычтя

затем (9.39) из (9.38),

можно

получить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д(х)=*(Тхх)(1+ах).

 

 

 

 

 

(9.40)

Д л я

нахождения а,,

надо провести

 

эксперимент

 

(измерить

Тх

и Тх)

при выключенном

разряде

при двух температурах

образца

7\ и Г , (комнатные температуры соответственно Го и Го). Из

двух

уравнений (9.39)

после небольших

преобразований

имеем

 

 

 

 

 

 

 

е

е{Т'х0)-{Тхо)

 

 

 

 

 

(

9 4

] )

 

 

 

 

 

е(Т,-Тх)-(Т2х)

 

'

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

є

= —;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т\-т0*

 

 

 

 

 

 

 

 

Определив

ах

согласно

(9.41) и подставив в (9.40),

находим

зна­

чения

соответствующих

рекомбинационных тепловых

потоков

q (х)

в зависимости

от расстояния от зонда

до источника

атомов вдоль

оси реакционной

трубки.

Последние

подставляются

вместо

Г (х)

в общую

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

рекомбинации.

 

 

 

 

 

 

 

 

 

Т а б л и ц а

7*

 

ю

 

і

 

о

 

 

к

сч

п

С*

§ Ё

а

 

 

3-і

= о 5

о >•

Я о

>,

О

О а>

4-і

О.

„ о> га

5 я

III

*=1

_ Я

S- а

 

S " =

В § 1

 

й ч >,

291

7,5

 

 

 

392

10,1

 

5,0

6,1

486

15,6

 

12,7

13,9

593

12,9

 

9,3

9,8

О

гмпература°К ізца, Н Q.

693

793

894

994

а (х) = 0 (из лучение уч­ тено)

12,3

9,6

13,7

18,4

7-Ю1

х ™f

О 0J

•в-кЗ

х л

а

8,8

5,1

10,4

16,5

а (дг) + 0 (из лучение уч­ тено)

7,7

5,6

11,1

17,6

Р (х) = const.

 

 

В

общем случае,

когда приходится учитывать

теплопередачу

от образца к зонду

как благодаря излучению, так и

теплопровод-

ности, удобно искать численное решение системы (9.38, 9.39) по стандартной программе СП-0043 [228], что и позволяет вычислить экспериментальные значения q (х). Затем с помощью ЭВМ набор теоретических кривых q (х, у, у') совмещается с экспериментальной зависимостью q (х) в результате последовательного перебора с оп­ ределенным шагом значений параметра у. Истинное значение у отвечает минимуму среднеквадратичного расхождения

j = у

Г Их, г у')

і ?

ZJ

q (х)

 

і

 

 

О необходимости строгого учета всех параметров, характеризую­ щих теплообмен зонда, и их влиянии на рассчитываемые значения у свидетельствуют полученные данные по исследованию взаимодей­ ствия атомарного водорода с поверхностью чистого поликристалли­ ческого вольфрама (W — 99,99%) (табл. 7).

§ 10. ОПРЕДЕЛЕНИЕ СТЕПЕНИ АТОМИЗАЦИИ ГАЗА В ПОТОКЕ И КОЭФФИЦИЕНТОВ РЕКОМБИНАЦИИ НА МЕТАЛЛИЧЕСКИХ ПРОВОЛОКАХ

Из исследуемого металла — проволоки

изготовляются

 

две сеточ­

ные

пробы,

помещаемые

в диффузионную

или проточную

труб­

ку [229] с произвольной, но известной скоростью течения

газа.

Решение соответствующей

краевой

задачи

получается

в

виде

трех

различных разложений по функциям Бесселя нулевого

порядка,

«сшиваемых» друг с другом в плоскостях

расположения

проб. Пред­

полагается,

что сеточная

проба не вносит возмущения

в

аэродина­

мику

газового потока. Д л я расчета

входящего потока

атомарного

газа

в пробу-сетку последнюю заменяют

поглощающей

поверхно­

стью, перекрывающей все сечение проточной трубки, с коэффициен­ том рекомбинации, связанным с коэффициентом рекомбинации материала пробы [218].

Настоящая методика дает возможность прямого определения коэффициентов рекомбинации атомов газа на поверхности метал­ лов и в то же время степени атомизации газа в потоке. При этом теп­ ловые эффекты проб могут быть рассчитаны, например, по методике Рогинского и Шехтер [1].

Пусть скорость движения потока V направлена вдоль оси л:. Урав­

нение неразрывности с

учетом гомогенной

рекомбинации

атомов

имеет вид

 

 

 

DAn

| - ( V n ) _ ^ =

- ^ - ,

(ЮЛ)

где п — концентрация; D — коэффициент диффузии атомов; k — константа скорости гомогенной реакции.

Рассмотрим аксиально

симметричную

задачу

для цилиндриче­

ской трубки диаметром 2R

с осью, совпадающей

с осью х.

Считаем,

что концентрация атомов

настолько мала,

что гомогенной

рекомби­

нацией можно пренебречь

и поглощение частиц наблюдается только

в двух сечениях трубки: х = Rd и х =

Rl (d <с I), где

находятся

поглощающие пробы. При этом число частиц,

входящих

в пробу

за единицу времени, равно

 

 

 

 

j>i-ynds.

 

 

(10.2)

Здесь с — средняя тепловая скорость частиц; у — коэффициент гетерогенной рекомбинации атомов; ds — элемент поверхности про­ бы. Интеграл берется по всей поверхности пробы.

Предположим, что проба не возмущает движущийся поток (с точки зрения аэродинамики), и рассмотрим приближение «разма­ зывания» пробы в сечении трубки. Будем считать, что проба запол­ няет все сечение трубки, но в каждом его элементе da она обладает коэффициентом рекомбинации у' = yds/do, где ds — часть поверх­ ности пробы, попадающая в элемент da. Обозначим dslda = k (г). Тогда

 

 

 

 

 

 

y' = k(r)y.

 

 

 

(10.3)

Поток, входящий в элемент сечения трубки da,

равен

 

 

 

 

 

 

 

с1д = у'-^-п(х,

г) da

 

 

(10.4)

при

х =

Rd,

х

= Rl, и уравнение

(10.1) приобретает вид

 

 

А п ~

~w 4г

 

= тг

п

(*•r)[k"

W 6

<* -

R d )

+ Ь W 6

( * - я о ь

 

 

 

 

 

 

 

 

 

 

 

 

(10.5)

Стоящая в правой части уравнения б-функция выражает тот

факт, что мы пренебрегаем размерами

проб

вдоль координаты

х.

Д л я

трубки

постоянного

сечения

можно считать dV/dx = 0.

Далее вводим безразмерные цилиндрические

переменные

| =

xIR,

р =

r/R.

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д°-п

VR

дп

 

дЧ

дп

1

_

 

 

 

 

 

дЪ?

D

dl

+

dp* +

dp

' р

~~

 

 

 

 

=

- g - Rn(l,

Р ) М Р ) « ( Б - < * )

+ М Р ) 8 ( & - 0 ] .

(10-6)

0 •< I < оо (рассматриваем бесконечно длинную трубку).

 

Граничные

условия

 

 

 

 

 

 

 

 

 

 

 

 

 

"(О,

Р) =

п0,

д

п & 0 )

= 0 ,

л(оо,

р) =

0,

 

 

 

 

 

 

 

1 }

+ e o n ( g ,

1) = 0,

 

 

(10.7)

где в 0 =

y0cR/2D,y0

коэффициент

рекомбинации

атомов на стен­

ках

трубки.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение уравнения (10.6) можно получить в виде

 

 

 

 

 

 

n(g,

p) =

S * , ( i ) * t ( p ) .

 

 

 

(Ю.8)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#/(р)

=

- М ^ Р ) .

 

 

 

 

(Ю.9)

 

 

 

 

а Л ( а < )

=

в о Л ( а / ) .

 

 

 

(Ю.І0)

Здесь У т

(а) — функция

Бесселя

порядка

т.

например, в [1011.

 

Корни

уравнения

(10.10) протабулированы,

Принимая

обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р< =

V

W

+

а '

_

"2D" '

 

(

1

0 Л 1 )

 

 

 

 

 

-•/"К2 /?2

,

і ,

v#

 

п

п

( 0 ,

 

 

 

 

=

К W

+ °' + " 2 D -

'

 

( І 0 Л 2 )

где

V = const,

получаем

 

 

 

 

 

 

 

 

 

 

 

f M|e - "« s + М е Ч ^ ( 9 = м Г е - р ' Е Н - ^ " е Ч

0 < | < d ,

 

c t < g < / ,

(10.13)

причем коэффициенты Mt, Nt должны быть определены системой алгебраических уравнений

 

 

М\ +

N] =

 

^

+

a*) '

 

 

 

 

 

 

JB(ai)(Ql

 

 

M\l

-М}

=

- (Nll -

N\)

е( *+'<*

 

 

М}и

-MY

 

= / v ! ' e W ,

 

 

 

- M\ ( - @d + pt) + PlM? + Ml (®d

+ qt) e < ^ + ^ - qtN\l

e ( W = 0,

 

M™ ( в , + pt)

-

ptM?

+ qtN№№

= 0.

 

Здесь Qd =

Rkd;

6, =

 

 

 

Будем считать

^ не завися­

щими от p, что отвечает равномерному

распределению витков про­

бы в сечении

трубки.

 

 

 

 

 

 

 

 

Соседние файлы в папке книги из ГПНТБ