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

Жданов Основы физических процессов 2007

.pdf
Скачиваний:
233
Добавлен:
16.08.2013
Размер:
9.04 Mб
Скачать

D ~ ( x)2 τ ,

(1.62)

где x - среднее смещение частицы при хаотических блужданиях, а τ -время между столкновениями. Для диффузии частиц в газе или плазме x соответствует средней длине свободного пробега λ . Учитывая, что τ = λvT , получаем

D ~ λvT ,

где vT = 3T m -средняя тепловая скорость частиц.

Элементарная кинетическая теория газов приводит к следующим выражениям для коэффициентов диффузии, вязкости и

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

 

 

1

 

 

 

 

 

 

D =

λ vT

 

,

 

 

(1.63)

3

 

 

 

 

 

1

 

 

 

 

 

 

η =

 

ρλ vT

= ρD ,

 

 

(1.64)

3

 

 

 

κ = 1

 

 

 

 

 

 

 

 

ρ c λ v

T

= c η

~ nD

,

(1.65)

3

 

 

 

 

V

v

 

 

 

 

 

 

 

 

 

 

 

 

 

где ρ = mn -массовая плотность, cV - удельная теплоемкость газа при постоянном объеме, Для одноатомных газов cV = (32)(km)

(если определять теплоемкость в системе СГС) Напомним, что температуру мы измеряем в энергетических единицах. Это означает, что во всех выражениях можно формально полагать постоянную Больцмана k равной единице.

Следует подчеркнуть, что выражения (1.63)- (1.65) получены для простого газа, образованного из одинаковых частиц массой m . В этом случае коэффициент D является по существу коэффициентом самодиффузии D11 , т.е. диффузией выделенных («ме-

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

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

61

ты x ). Рассмотрим плоскость, перпендикулярную оси х, и выделим на ней еди-

ничную площадку. Будем предполагать, что молекулы, летящие по направлению к этой площадке, претерпевают последние столкновения на расстоянии

средней длины свободного пробега частиц λ и далее движутся к ней без

столкновений. Учтем также, что при изотропном распределении молекул по направлению скоростей в заданном направлении движется в среднем 1/6

часть молекул. Тогда их поток определяется выражением N = (16)nvT . Результирующая плотность потока частиц, пересекающих выделенную площадку в направлении оси x , определяется как разность потоков молекул, летящих с левой и правой стороны относительно площадки,

J x

= N + N =

1

n(x λ)vT

1

n(x + λ)vT ,

(1.66)

 

6

 

6

 

 

 

где координата х соответствует положению площадки.

 

При оценке средней длины свободного пробега λ 1 nσ можно вос-

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

моделью молекул - твердых

упругих шаров. В

этом случае

σ = π a2 , где a - диаметр шара. Для рассматриваемого слабо неоднородного состояния газа плотность n в выражении для средней длины свободного

пробега, а следовательно и длинуλ , можно приближенно полагать постоянными. Воспользуемся также малостью этой величины. Например, при нормальных условиях ( T = 273K и p =1 атм) плотность газа равна числу Лош-

мидта n 2.7 1019 см-3, , а λ имеет порядок 10-6 см., если a = 5Ao . Это позволяет использовать приближенные разложения плотности в ряды вида

n(x λ)n(x)λ

n

, n(x + λ)n(x)+ λ

n ,

и представить J x ( 1.66) как

x

 

x

 

 

 

 

 

1

n

 

J x ≈ −

 

λ vT x .

(1.67)

3

Это выражение по своей структуре аналогично определению диффузионного потока (1.59), направленного вдоль оси x Коэффициент пропорциональности

в нем и есть коэффициент самодиффузииD11 (1.63).

Аналогичным образом можно получить выражение и для потока тепла qx . Учитывая, что средняя энергия теплового движения частиц равна cV mT ,

где T -температура газа, которая теперь меняется вдоль оси x (плотность n предполагается постоянной), находим

 

1

ρvT cV [T (x λ)T (x + λ)]≈ −

1

T

qx =

 

 

ρcV λ vT x .

6

3

Сравнение с (1.60) показывает, что смысл коэффициента теплопроводности имеет величина

62

κ =

1

ρ c

λ v

T

~ nD ,

(1.68)

 

 

3

V

 

 

 

 

 

 

 

 

 

которая совпадает с определением (1.65).

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

рый движется со скоростью uz (вдоль оси z ). Площадка, через которую рас-

сматривается перенос импульса, располагается в некотором положении x , перпендикулярно оси x . Рассматривая разность значений импульсов, перено-

симых из областей, отстоящих от площадки на расстояниях x λ и

x + λ ,

для результирующего потока импульса получаем

 

1

uz

 

τ xz = − 3 ρλ vT

 

.

(1.69)

x

Коэффициент пропорциональности в этом выражении соответствует коэффициенту вязкости η , определяемому формулой (1.64)

Получаемые в элементарной кинетической теории выражения для коэффициентов переноса дают правильные зависимости от основных параметров газа, отличаясь от точных значений лишь численными коэффициентами. При этом в соответствии с результатами строгой кинетической теории газа коэффициенты переноса для модели молекул – твердых упругих шаров – связаны соотношениями [50 ]:

ρD =

6

η

, κ =

5

c η ,

 

(1.70)

5

2

 

где

 

 

 

 

V

 

 

 

 

 

1

 

 

1

 

 

η = 0.354ρλ v ,

 

λ =

=

.

(1.71)

 

nσ

nπa2

 

 

 

 

 

 

 

 

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

ционально T 1 2 .

Для

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

D ~ (1 p)T 3 2 , где

p -давление газа. Характер температурной зави-

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

63

Вязкость плазмы

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

Выражение для коэффициента вязкости (1.64), даваемое элементарной теорией, можно представить также в виде:

η =

1

mnλ vT

=

1

mnτvT2

= pτ

, p = nT , (1.72)

3

3

 

 

 

 

 

 

где τ - характерное время свободного пробега между столкновениями (время релаксации).

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

Принимая для оценок, что сечение взаимодействий молекулы σ = πa2 = const , получаем

η ~ mn

1

 

T 1 2

~ m

1 2

T

1 2

.

(1.73)

 

 

 

 

 

 

 

 

 

 

 

nσ m

 

 

 

 

 

 

Существенно, что вязкость не зависит от давления газа.

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

делены [51]. Вклад ионов в вязкость определяется

при этом време-

нем релаксации для ионов τii (1.45), так что

 

η

i

= p τ

~ M 1 2T 5 2 Λ .

(1.74)

 

i ii

i i

 

64

Электронная вязкость определяется как электрон-электронными, так и электрон-ионными столкновениями. Поскольку, однако

τee = 2τei , то с учетом выражения для τee (1.44) имеем

ηe ~ me1 2Te5 2 Λ .

(1.75)

Очевидно, что при Te = Ti

ηe ηi = (me Mi )12 <<1,

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

вязкостью ионов для неизотермической плазмы при Te >> Ti , если

выполнено условие Te Ti ~ (Mi me )15 .

Существенно, что в соответствии с (1.74) ионная вязкость гораздо сильнее зависит от температуры, чем вязкость нейтральной компоненты. Отметим, что коэффициенты вязкости ионов и электронов для полностью ионизованной плазмы, получаемые метода-

ми строгой кинетической теории, даются выражениями [51]:

ηi = 0.96 niTiτii , ηe = 0.73 neTeτei .

(1.76)

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

Теплопроводность слабоионизованной плазмы определяется переносом тепла только нейтральными частицами. В соответствии с (1.65)

κ ~η m ~ ( p m)τ ~ m 1 2T 1 2 .

(1.77)

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

ки, основанные на (1.77) и

выражениях для τee (1.44) и τii

(1.45),

дают

 

 

 

 

 

 

 

 

 

 

κ

e

~ m1 2T 5 2

/ Λ ,

κ

i

~ M 1 2T 5 2

/ Λ .

(1.78)

 

e

e

 

 

i

i

 

 

При Te = Ti имеем

65

κi κe = (me Mi )12 <<1.

Это означает, что основной вклад в теплопроводность полностью ионизованной плазмы вносит электронная теплопроводность. Вклад теплопроводности ионов может оказаться сравнимым с электронной теплопроводностью для неизотермической плазмы с

Ti > Te при выполнении условия Ti Te ~ (Mi me )15 .

Коэффициенты теплопроводности, получаемые в строгой кинетической теории, имеют вид [51,52]

κe = 3.16neTeτei / me , κi = 3.9niTiτii / Mi .

(1.79)

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

Диффузия в плазме

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

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

ur =

1

nk mk urk , ρ = nk mk .

(1.80)

ρ

 

k

k

 

 

v

-соответственно масса, плотность (концентра-

Здесь mk , nk и uk

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

66

( ne , ni <<1) скорость u фактически совпадает со скоростью ней-

трального компонента ua . Для полностью ионизованной плазмы из-за условия me << Mi эта скорость практически равна скорости

ионов ui .

Плотность диффузионного потока частиц k -го компонента (будем называть ее просто диффузионным потоком ) определяется как

Jk = nk (urk ur).

(1.81)

Пусть концентрации компонентов nk меняются вдоль оси x . Если

к плазме приложено еще электрическое поле E , то диффузионный поток k -го компонента вдоль направления x может быть представлен в виде:

J

 

= −D

dnk

+ Z

 

μ

 

n

E .

(1.82)

 

dx

 

 

 

k

k

 

k

 

k

k

 

 

Здесь Dk и μk - соответственно коэффициенты диффузии и подвижности k -го компонента в плазме, Zk -кратность заряда (зарядовое число), для электронов Ze = −1 . Из общих термо-

динамических соображений можно показать , что коэффициенты диффузии и подвижности связаны соотношением Эйнштейна:

μk =

e

Dk .

(1.83)

T

 

 

 

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

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

67

ром силы, действующие в выделенном объеме плазмы на частицы k -го компонента за счет неоднородности концентрации (парциального давления) и благодаря наличию электрического поля, уравновешиваются силами трения, действующими в среднем на частицы данного компонента со стороны частиц других компонентов плазмы. Последняя величина носит название диффузионной силы трения и равна, очевидно, усредненной величине суммарного импульса, передаваемого в столкновениях частиц данного сорта с частицами других сортов. Уравнение баланса импульса для k -го компонента плазмы можно пред-ставить в виде [51,53]:

dpk

 

+ nk Zk eE = Rk ,

(1.84)

dx

 

 

 

 

Rk

 

= −nk μkjνkj (uk u j ).

(1.85)

 

 

 

 

j

 

Поясним физический

смысл

отдельных членов, входящих

в уравнение

(1.84). Выделим в плазме малый цилиндрический объем толщиной dx с поперечным сечением S . Сила, действующая на выделенный элемент за счет градиента парциального давления k -го компонента , равна, очевидно,

[ p

 

(p

 

 

 

) ]S = −Sdp

 

 

 

dp

k

 

k

k

+ dp

k

k

=

 

Sdx .

 

 

 

 

 

 

 

 

dx

В отсутствие электрического поля эта величина уравновешивается силой диффузионного трения, равной Rk Sdx , где Rk равна силе диффузионного тре-

ния, возникающей в единице объема плазмы. В результате (dpk dx)= Rk ,

где левая часть этого соотношения есть действующая на частицы k -го компо-

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

(1.84 ).

Вывод выражения для диффузионной силы трения Rk (1.85)

рассматривался в [51,53]. Из общих соображений ясно, что средняя передача импульса при столкновениях частиц должна быть пропорциональна разности направленных скоростей компонентов. Для компонентов, движущихся с одинаковыми скоростями, она отсутствует (при столкновениях частиц одного сорта импульс в среднем не передается). Кроме того, очевидно, что эта величина должна

68

быть пропорциональна плотности частиц k-го компонента и средней частоте столкновений νkj с частицами других сортов. Приве-

денная массаμkj появляется в выражении, поскольку столкновения частиц при вычислении Rk рассматриваются в системе отсчета,

связанной с их центром масс.

Обратимся к случаю слабоионизованной плазмы. Уравнения баланса импульса для электронов и ионов (с кратностью заряда Z )

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

dpdxe neeE = −ne meνea (ue ua ),

 

dpi

+ n ZeE = −n μ ν

(u

i

u

a

).

(1.86)

 

 

 

 

dx

i

i ia ia

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

учитывается, что μea me , заметим,

что для ионов и

атомов одинаковой

массы

μia = ma

2 . Мы полагаем, что

градиент

полного давления

dp dx = 0 .

Это позволяет

нам не

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

Диффузионные потоки в

рассматриваемом нами случае

определены как

 

Je = ne (ue ua ),

Ji = ni (ui ua ).

Тогда с использованием (1.86) получаем

J

 

= −D

dne

μ

n

E

,

 

dx

 

e

e

e

e

 

 

J

 

= −D

 

dni

 

+ Z μ n E ,

 

 

 

dx

 

 

i

i

 

 

i

i

 

где использованы

 

соотношения pe

= neTe

и

коэффициентов диффузии получаем

 

 

 

 

De =

Te

 

,

Di =

Ti

 

 

m ν

 

μ ν

 

 

 

 

e

ea

 

ia ia

(1.87)

(1.88)

pi = niTi . Для

. (1.89)

69

Коэффициенты подвижности μe и μi выражаются через De и Di с помощью соотношения (1.83).

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

νea = na σea (g)g . , νia = na σia (g)g ,

(1.90)

где g -модуль относительной скорости частиц, а означает усре-

днение по максвелловскому распределению частиц по скоростям.

Для модели твердых упругих шаров νea, naσeaveT , νia, naσiaviT и мы приходим к оценкам

De ~

T 3 2

, Di ~

T 3 2

 

e

i

 

 

 

 

.

(1.91)

me1 2σea

Mi1 2σia

Как видно, при Te = Ti и

σea ~ σia имеем

Di De ~ (me Mi )1 2 ,

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

В случае полностью ионизованной плазмы достаточно записать одно уравнение баланса вида

 

dpe

 

n

eE = −n m ν

ei

(u

e

u

).

(1.92)

 

 

 

 

 

 

dx

 

 

e

 

 

 

 

 

e e

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

2 )τee , а τee определено выражением

Здесьνei =1 τei

,

 

где τei = (1

 

(1.53). В

результате

 

 

диффузионный

поток

электронов

Je = ne (ue

ui

)

определяется как

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

= −D

 

dne

 

μ

n

E ,

 

 

 

(1.93)

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

e

 

 

 

ei dx

e

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T 5 2

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

e

 

 

Dei =

 

 

 

 

 

 

 

 

e

 

 

 

 

 

, μe =

 

Dei .

(1.94)

 

 

 

π

 

n

e4 m1 2 Λ

 

T

 

 

4 2

 

 

e

 

 

e

 

 

 

 

 

e

 

 

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

70