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

книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы

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

404 Рл. 6. Взаимодействие электромагн. полей с сист. многих ЧастпйЦ

электронным флуктуациям, ответственным за пик вблизи сор_.

5. Боковой максимум смещается на крыло линии спектра и сужается при возрастании электронной темпе­ ратуры выше температуры ионов. Это можно понять, если вспомнить, что в плазме возбуждаются коллективные ион­

ные колебания с фазовой скоростью (pf/k = ]/ (Э_/тп+. 6. При наличии относительного дрейфового движения электронов и ионов профиль линии спектра становится асимметричным и один из пиков сильно возрастает. Этот эффект тоже связан с возбуждением коллективных ион­ ных волн. Как уже указывалось в гл. 3, относительный дрейф электронов и ионов может приводить даже к неустой­

чивости.

4.4. Применение флуктуационно-диссипационной теоремы

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

Весьма последовательным шагом было бы установле­ ние связи между результатами теории рассеяния и вычис­ лениями проводимости с помощью флуктуационно-дисси­ пационной теоремы (см. § 5 гл. 1):

0 £ ■ Щ ег$ ] - Um

<I » (к, ю) Р>„. (4.56)

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

§5. ПРОЦЕССЫ ПЕРЕНОСА ИЗЛУЧЕНИЯ

5.1.Уравнение переноса излучения

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

§ 5. Процессы переноса излучения

405

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

4 - < i -

(5.1)

Здесь L — характерная длина изменения величины е<о (г). Из уравнений Максвелла для сплошной среды

V .Hw= 0,

~

V.(euEu) = 0,

~

(5.2)

У X Н(о —

~

^ь)8об(!)Е(,),

У х Ещ =

 

 

получим волновое

уравнение

 

 

 

ДЁИ+

еа (г) Ёи-I- V (15». Ё „) = 0.

(5.3)

В рамках приближения геометрической оптики мы можем (локально) задать решение в виде плоской волны, распространяющейся в направлении s, и, учитывая усло­ вие (5.1), опустить третий член в выражении (5.3). Тогда

дЩа

|

со2

„ ^

п

(5.4)

ds2

+

С2

8“Е“ — 0.

 

С помощью метода Вентцеля — Крамерса — Бриллюэна (ВКБ) найдем приближенное решение в виде

Еш(s) = Cwехр (* “г j V г<*ds )

(5.5)

для электрического поля и соответственно

Нш(s) = Da exp ( i j У ea ds}

(5.6)

для магнитного поля. Из этих решений получаем абсо­ лютную величину вектора Пойнтинга

S = (2л)-1 j

ЁШ1ЯЙ2ехр [— i (coj—со2) *]«&!)!dco2- (5.7)

Усреднение по

периоду времени, достаточно большому

по сравнению с периодами всех колебаний, дающих вклад

406 Гл. 6. Взаимодействие электромагн. полей с сист. многих частиц

в выражение (5.7), приводит к формуле

(S) = я-1Re |

j

Еа1Н%2exp [—i (coi—<o2) £] dcoi cfto2}. (5.8)

 

(01, G>2>0

 

Подставив сюда

решения в виде (5.5) и (5.6), найдем

{S) = я"1 Re |

j

C0lD l2exp £ —

j (xffll+ хИ2) dsj x

 

0)1,(O2>0

 

X exp | i [c-1 j ((О1ПМ1—Юг^ма) ds—

 

 

(«! — (o2)

| cftoi dw21 . (5.9)

Здесь мы ввели показатель преломления па и коэффи­ циент пространственного поглощения ха:

/гт = Re еш и хт = — 1т ]/е ш«

l / " у-

• (5.10)

С

Г Б о L

ГС<|) J

Рассмотрим теперь волновой пакет, содержащий ча­ стоты (Dj и со2, близкие к со. Тогда можно использовать

приближенное равенство

xra]

+

хШ2 «

2хи. В результа­

те мы имеем

 

 

 

 

 

 

 

П(02®2

 

 

 

d ИщО) _ (01 — й>2

(5.11)

 

(0)i—со2) da

Ug (со)

 

где ug (со) — локальная

групповая скорость

для

волны

с частотой со. Таким образом, получим

 

 

(S)a = я-1 exp £ — j хшdsj Re

j

CaiDZ2X

 

 

 

 

(Hi, G )2 > 0

 

 

 

X exp [i (coj

co2) ( j

 

—«)] dtoidaz.

(5.12)

Следовательно, величина

(S)a удовлетворяет

дифферен­

циальному уравнению

 

 

 

 

 

 

 

d(S)a

-*ш (S)a

 

1

d ( S ) a

 

(5.13)

d s

 

u g (о)

d t

 

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

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

§ 5. Процессы переноса излучения

407

этом интенсивность излучения равна

г __ (S)a

ю Дй Дсо '

Используя закон Спелля, записанный в виде

«aAQ = const,

а также уравнение (5.13) и формулу (5.14),

уравнение переноса излучения х)

(5.14)

(5.15)

получим

t

^

I и, _i

д I и,

7а

 

(5.16)

7,

Я?

2 Г

а

2

Пи,

.

 

 

«ш

os

Пи,

 

 

Очевидно, что в правой части уравнения (5.16) представ­ лены эффекты поглощения волн и отсутствуют эффекты излучения и рассеяния. Это не удивительно, поскольку в уравнениях Максвелла (5.2) не учитываются ни тепло­ вое излучение, ни флуктуации. Если же учесть эффекты излучения, то вместо (5.16) мы получим уравнение

 

t

d Iи, .

д Iи,

/<а ^о)7а

(5.17)

 

и8

п%,

ds на

гса

 

 

в

котором через

/а обозначен

коэффициент излучения

(т.

е. мощность излучения на частоте со единицей объема

среды в единичный телесный угол в направлении луча s и в единичном интервале времени и частоты). В это урав­ нение мы можем подставить выражение для коэффициента излучения, полученное в гл. 5. Для системы, близкой к равновесной, можно использовать закон Кирхгофа, со­

гласно которому

НаВи,-

(5.18)

U =

Здесь Ba — спектральная плотность излучения черного тела, определяемая формулой Планка:

Ви

fid)3

4я3с2 0 .

(5.19)

4я3с2 е>м>/@_1

х) Обратим внимание на то обстоятельство, что

используемое

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

величину ге“* под знак дифференциального оператора d/ds или оставить вне его. При более строгом выводе надо учитывать член

уравнения (5.3), содержащий величину Ей) • V еш. Для этого необхо­ димо учесть различие направлений поляризации волн и использо­ вать ВКБ-метод в приближении более высокого порядка по )7Ь.

408 Гл. 6. Взаимодействие электромазн. полей с сист. многих частиц

[Соотношение (5.19) учитывает два возможных направле­ ния поляризации.]

Если в уравнение переноса также включить эффекты рассеяния, то — по аналогии с уравнением Больцмана — можно представить вклад этих процессов рассеяния в виде дополнительного «члена столкновений» в правой части

уравнения (5.17)

 

 

 

/

4

6/ю\

Л-

[ [Ia' ( Q ' ) - I a (Q)]x

 

'

па‘

/расс

Им

J

 

 

 

 

 

 

X

(со, й' —Q) da' dQ’.

(5.20)

Здесь аШ' (со, Й' — й) — дифференциальное сечение рас­ сеяния, определяемое формулой (4.27).

5.2. Замечания по поводу граничных условий

Уравнение переноса, рассмотренное выше, следует дополнить граничными условиями. Получение этих гра­ ничных условий — довольно сложная задача. Предполо­ жим, что плазма ограничена объемом V. Тогда распре­ деление полей излучения вне этого объема дает граничные условия для решения уравнения переноса.

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

при

наличии пограничного слоя конечной толщины

(0 <

х •< а), когда поле в нем может быть описано с помо­

щью комплексной диэлектрической проницаемости, урав­

 

§ 5. Процессы, переноса излучения

409

нение (5.3)

принимает вид

 

 

 

- ^ - + кг (х)Е = 0,

где

Л2 = -§-е(х).

(5.21)

Предположим далее, что к (х)

= к0 при х ^ 0,

а при

х ^ а к (х)

== Art. Тогда решения уравнения (5.21)

можно

записать как

 

 

 

и

Е = Е0 (eiho*+ ре-«о*)

для xs^

(5.22)

E = E0xeihlx

для

х > а .

 

 

Здесь р и т — коэффициенты отражения и прохождения для амплитуд (но не интенсивностей!) полей.

Чтобы найти решение уравнения (5.21) в интервале значений х [0, а], используем несколько модифицирован­ ный по форме обычный метод медленно меняющейся фазы (см. [11]), в соответствии с которым будем искать решение в виде

Е (х) = А ехр | i j [& (х) — ф(х)] dx | .

(5.23)

Это приводит к дифференциальному уравнению Риккати для ф (х):

ф’ + Икф = к' + 1ф2.

(5.24)

Можно ожидать, что величина ф (х) в большинстве случаев мала, поскольку при ф (х) = 0 решение (5.23) переходит в решение, совпадающее с приближением ВКБ. Логарифмическое дифференцирование выражений (5.22) и (5.23) с последующим их приравниванием при х — а дает граничное условие

ф (а) = 0,

(5.25)

а при х = 0 — коэффициент отражения

_ < Н 0 ) _ _

(5.26)

н 2к0-ф(0)

Таким образом, рассматриваемая задача будет решена, если с помощью уравнения (5.24) и условия (5.25) мы смо­ жем определить ф (0).

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

410 Гл. 6. Взаимодействие электромагн. полей с сист. многих частиц

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

Используя метод вариации произвольной постоянной и условие (5.25), найдем

X

 

X

Z

 

ф (х) = exp [ —2i

к (£)

j

exp |^2i ^ к (£) dgj X

 

b

o

 

b

 

 

 

 

X [k' (z) -f 1фг (z)] dz.

(5.27)

Если ф мало, то можно рассматривать лишь линейное приближение. Тогда

а2

ф(0) = — j exp £2i j к (|)

к' (z) dz

(5.28)

о0

икоэффициент отражения определяется простым соотно­ шением р = ф (0)/2к0. Поскольку интенсивность рассма­

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

ах

j exp ^2i j к (g) dlj к' (x) dx

о0

(5.29)

4

В пределе, когда a-> 0, из (5.29) вытекает обычная фор­ мула Френеля

lim R —

У е ( a ) - V e (0) 2

(5.30)

а-*0

2 Уе (0)

 

Входящую в (5.29) величину

Р(*) =

1

2

УГ{0)

 

d У е (х)

(5.31)

dx

 

согласно формуле Френеля (5.30), можно трактовать как «дифференциальный коэффициент отражения» Этот диф­

§ 5. Процессы переноса излучения

411

ференциальный коэффициент отражения умножается на

 

 

X

 

 

e x p [2 i^ j / F ( |j d g ] ,

(5.32)

 

 

о

 

которая

учитывает

поглощение (через мнимую

часть

Im |/*ё)

и явления

интерференции в выделенном

нами

слое плазмы.

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

Vr

Ф и г . 34. Коэффициент отражения R как функция толщины погра­

ничного слоя плазмы а, измеренной в длинах волны X; co„i =

0,2 со,

A l =

О О .

 

 

1 — для распределения (5.33) плотности

электронов в

слое; 2 — для

линей­

ного распределения плотности электронов

в слое.

 

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

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

412 Гл. 6. Взаимодействие электромаен. полей с сист. многих частиц

Риккати. Нами были проведены такие вычисления для пограничного слоя плазмы с изменением плотности по

Ф и г. 35. Коэффициент отражения R как функция толщины по­ граничного слоя плазмы а, измеренной в длинах волны А,; Шр1 =

= 0,8 о.

кривая 1 — A i =

оо; кривая 2 — Л , =

10.

 

закону

 

 

 

(х) =

щ ^1 —cos-^-j

.

(5.33)

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

Re е

р

Im е

1

/

2 \ 1/и сор-

In Л

1

<0р_

0)2 >

8

\

~п }

о 3

Л

Л

о)3

(5.34)

На фиг. 34 и 35 приведены зависимости для коэффициента отражения соответственно при о)р1/(о = 0,2 и сол1/(о == 0,8 для величины A j, равной 10 и оо.

 

ЦИТИРОВАННАЯ ЛИТЕРАТУРА

1.

Випетапп О., Phys. Rev., 115,

503

(1959).

2.

Kahn F. D., Journ. Fluid Mech.,

14,

321 (1962),

3. Fried В■D., Phys. Fluids, 2, 337 (1959),

 

Jlumepamypd

413

4.

Furth H. P., Phys. Fluids, 6, 48 (1963).

 

5.

Weibel E. S., Phys. Rev. Letters, 2, 83 (1959).

 

6.

Dawson J ., ОЪегтап C., Phys. Fluids, 5, 517 (1962).

7.

Dawson J., Oberman C., Phys. Fluids, 6, 394 (1963).

8.

Rosenbluth M. N., Rostoker N., Phys. Fluids, 5, 776 (1962).

9.

Rernstein I. B., Trehan S. K., Weenink M. P. H.,

Nucl. Fusion,

 

4, 61 (1964).

 

10.Rostoker N ., Rosenbluth M. N., Phys. Fluids, 3, 1 (1960).

11.Geffcken W., Ann. Phys. (Leipzig), 40, 835 (1941).

Дополнительная литература

Bekefi G., Radiation Processes in Plasmas, W iley, New York (1966) (см. перевод: Дж. Бекефи, Радиационные процессы в плазме, изд-во

«Мир», 1971).

Dawson / . , Lecture Notes, Princeton University, Princeton, N. J.,

1964.

Yol. I

(eds.

Dawson J., в книге Advances in Plasma Physics,

A. Simon, W. B. Thompson), W iley (Interscience), New

York,

1968

(см. перевод в сборнике «Физика высокотемпературной плазмы»,

изд-во «Мир», 1972).

Dupree Т. Н„ Phys. Fluids, 6, 1714 (1963).

Dupree Т. Н., Phys. Fluids, 7, 923 (1964).

Birmingham T., Dawson J., Oberman C., Phys. Fluids, 8, 297 (1965).

Bernstein I. B., Trehan S. K., Nucl. Fision, 1, 3 (1961).

Ramsden S. А., в книге Physics of Hot Plasmas (eds. B. J. Rye, J. C. Taylor), Oliver and Boyd, Edinburgh, 1970.

Kofink W„ Ann. Phys. (Leipzig), 1, 119 (1947).

Rawer K., Ann. Phys. (Leipzig), 35, 385 (1939). Margenau H ., Phys. Rev., 109, 6, 1958.

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