Научная статья на тему 'Нагрев плазмы на разрывных МГД-течениях вблизи области магнитного пересоединения'

Нагрев плазмы на разрывных МГД-течениях вблизи области магнитного пересоединения Текст научной статьи по специальности «Физика»

CC BY
72
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПЛАЗМА / PLASMA / МАГНИТНАЯ ГИДРОДИНАМИКА / MAGNETOHYDRODYNAMICS / РАЗРЫВЫ / DISCONTINUITIES

Аннотация научной статьи по физике, автор научной работы — Леденцов Леонид Сергеевич, Сомов Борис Всеволодович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Леденцов Леонид Сергеевич, Сомов Борис Всеволодович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Нагрев плазмы на разрывных МГД-течениях вблизи области магнитного пересоединения»

Нагрев плазмы на разрывных МГД-течениях вблизи области магнитного пересоединения

Л. С. Леденцовa, Б. В. Сомовb

Государственный астрономический институт имени П. К. Штернберга Московского государственного университета имени М. В. Ломоносова, отдел физики Солнца.

Россия, 119991, Москва, Университетский пр-т, д. 13.

E-mail: a [email protected], b [email protected] Статья поступила 02.12.2012, подписана в печать 21.12.2012.

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

Ключевые слова: плазма, магнитная гидродинамика, разрывы.

УДК: 524.3-78, 523.9. PACS: 96.60.Q-, 96.60.qe, 96.60.Iv, 95.30.Qd, 96.50.Ry.

Введение

Разрывные течения плазмы в магнитном поле реализуются в очень широком классе физических условий. Как следствие они присутствуют в различного рода технических установках и устройствах, имеющих практическое значение [1-3], в лабораторных и численных экспериментах (например, [4-6]), в космических условиях — особенно в связи с эффектом магнитного пересоединения [7-10]. При этом существенным, как правило, является вопрос о нагреве плазмы до самых высоких температур [11].

С теоретической точки зрения при исследовании свойств разрывных течений плазмы принципиальную роль играют следующие факты. Уравнения обычной гидродинамики имеют, как известно, лишь два типа разрывных решений: тангенциальный разрыв и ударная волна [12]. Это связано с тем, что в гидродинамике существует только два типа волн малой амплитуды: энтропийные волны, движущиеся вместе со средой, и звуковые волны, распространяющиеся со скоростью звука. В магнитной гидродинамике (МГД) присутствие магнитного поля в плазме приводит к наличию четырех типов волн малой амплитуды. При этом, в отличие от обычной гидродинамики, нет минимальной скорости типа скорости звука, с которой могут распространяться малые возмущения. Как следствие в МГД картина разрывных течений гораздо более богатая, чем в гидродинамике. Существуют быстрые, медленные, альвенов-ские ударные волны и другие разрывы [13, 14], причем между разрывами различных типов могут существовать непрерывные переходы.

Напомним, что при переходе через поверхность разрыва происходит резкое изменение, скачок параметров плазмы. Тип разрывного МГД-решения, т. е. его характер, определяется изменениями плотности плазмы, скорости ее течения и вмороженного в нее магнитного поля. Кроме того, на поверхности разрыва происходит нагрев плазмы, величина которого, разумеется, тоже зависит от типа МГД-разрыва, но не определяет его классификационные признаки: непрерывность или ска-

чок плотности, наличие или отсутствие перпендикулярных составляющих скорости и магнитного поля В± .

Цель настоящей работы — изучение характера нагрева плазмы разрывными течениями различных типов. На основе полной системы граничных условий на разрыве мы получаем выражение, описывающее в явном виде величину изменения внутренней энергии плазмы. Исходя из установленных ранее соответствий между стандартной классификацией МГД-разрывов и параметрами течения плазмы [15], изучаем свойства нагрева плазмы на различных типах разрывов.

В разд. 1 приведены общая система граничных условий на поверхности МГД-разрыва и соотношения, связывающие конфигурацию магнитного поля вблизи поверхности разрыва с параметрами плазмы (плотность, поток массы и поток магнитного поля через разрыв). Раздел 2 посвящен анализу геометрии магнитного поля для различных типов разрывных течений. В разд. 3 выводится уравнение, описывающее скачок внутренней энергии на поверхности разрыва. Влияние характеристик плазмы и МГД-разрыва на величину этого скачка изучается в разд. 4.

1. Граничные условия на разрыве

В окрестности МГД-разрыва плотность плазмы, ее давление, скорость течения, направление и напряженность магнитного поля могут изменяться скачком на расстоянии, сравнимом с длиной свободного пробега частиц. Физические процессы внутри такого скачка определяются кинетическими явлениями в плазме, как ламинарными, так и турбулентными [16]. В приближении диссипативной МГД внутренняя структура разрывного течения определяется диссипативными коэффициентами переноса (вязкостью и проводимостью), а также теплопроводностью [17, 18]. Однако в приближении идеальной МГД этот скачок имеет нулевую толщину, т. е. происходит на некоторой поверхности разрыва.

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

У

VI /в2

/у2

X

У2

Р1

■Р2

О X

Рис. 1. Изменение магнитного поля В, поля скоростей V и плотности р плазмы на фронте ударной волны X = 0

координат, в которой наблюдатель перемещается вместе с поверхностью разрыва, расположенной в плоскости (у,г). Однородная плазма втекает с постоянной скоростью в разрыв слева и вытекает из него справа. Находясь в рамках идеальной МГД, мы пренебрегаем вязкостью, теплопроводностью и электрическим сопротивлением плазмы. Тогда граничные условия для уравнений МГД на разрыве могут быть записаны в виде следующих законов сохранения [13, 19, 20]:

{Вх} = 0,

|р^х } = 0, {ихВу - иуБх} = 0 {VхВг - игВх} = 0, 1

4п

|р^у - ^ БхБу\ = а

1

PVхVz - 4ПВхВг\ = 0,

}

{

[р + р^2 + ^ =

^х (V2 + е + р) + 4П (В^х - (V ■ В)Вх)} = 0.

(1) (2)

(3)

(4)

(5)

(6)

(7)

(8)

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

{Вх } = Вх2 - Вх1 = 0,

или, что то же самое, сохранение магнитного потока. Величины, отмеченные индексом 1, относятся к левой стороне (рис. 1), соответствующей набегающему потоку плазмы, а индексом 2 — убегающему.

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

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

Для двумерного случая течений плазмы в плоскости (х, у) в [15] из системы уравнений (1)-(7) получена формула, связывающая углы наклона вектора магнитного поля к нормали, опущенной на поверхность разрыва х = 0, с величинами плотности р1 и р2, потоком магнитного поля Вх и потоком массы т = рvх:

^2 =

2 (В2х/4п - т27) + т2{г}

tg 01,

2 (В|/4п - т27) - т2{г}

где tg 0 = Ву/Вх, г = 1/р. Здесь и далее везде тильдой помечены средние значения величин; так, 7 = (г1+г2)/2. Перепишем это уравнение, раскрыв скачки {г} и средние значения 7:

,а т2 ■ 4ш1/В2х - 1

^ 02 = ——т tg 01. т2 ■ 4пг2/В2 - 1

Обозначим т2й = В2х/4пт1 и т2оп = В2х/4пт2. Ранее было показано, что т^ и топ — величины потока массы, протекающего через волну выключения и волну включения соответственно [15]. Отметим, что тоц < топ, так как по теореме Цемплена [12] на разрыве г2 < г1. Уравнение для углов наклона магнитного поля примет простой вид

^ 02=тт-4 ^ 01.

(9)

т2/т2п - 1

Кроме того, в [15] из условия существования нетривиальных решений рассматриваемой системы МГД уравнений найдено ограничение на возможные значения потока массы через разрыв, или

д,

4пг'

или

т2 <

(10)

В2 + В 2

т2 > —т-^. 4пг

(11)

22 т2 < тД,

2 2 2 т2 > тд + т^.

Введем обозначения тД = Вх/4п7 и т^ = Ву/4п7. Величина тд соответствует альвеновскому потоку массы [15]. Так как г2 < 7 < г1, то т^ < тд < топ. Условия (10), (11) примут вид

2 2 (12)

(13)

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

2. Конфигурация поля

Уравнение (9) вместе с условиями (12)-(13) описывает зависимость углов наклона магнитного поля от величины потока массы, протекающего через разрыв. Функция (9) содержит пару особенностей. При т = тоц имеем 02 = 0, а при т = топ вне зависимости от угла 01 получаем 02 ^ ±п/2.

Вид зависимости 02(т2, 01) можно задать либо двумя параметрами тоц и топ, либо, например, величинами р1 и {р} . Нас интересуют классификационные признаки разрывных течений, т. е. качественные изменения соотношения между углами 01 и 02 при варьиро-

вании т2. Поэтому будем рассматривать формулу (9) пока без каких-либо конкретных приложений к определенным физическим условиям в среде. Значения параметров уравнения (9) и условий (12)-(13) выберем из соображений наглядности и простоты. Пусть значения т2((, т|, т^п соотносятся как 3 : 4 : 6. Квадрат потока массы будем измерять в единицах т|/4.

На рис. 2 показаны графики зависимостей в2(т2) для трех значений угла падения магнитного поля на разрыв: в1 = 5°, 25° и 45°. Кривые 02(т2) ведут себя одинаковым образом. Во-первых, они пересекаются в одной точке при т2 = т2й = 3. Во-вторых, для каждой кривой 02 ^ — 01 при т2 ^ т2к = 4. В-третьих, все они имеют область, не удовлетворяющую условиям (12)-(13). Располагается эта область в окрестности точки т2 = т2 = 6.

02 60

40

20

О

-20

-40

-60

I п ! ш IV______

" N Ч \ ~~~~

\\

1 \

Ч \

M i\ i_L

О

10 12 14 16 18 /п

Рис. 2. Зависимости угла наклона магнитного поля в2 за плоскостью разрыва от квадрата потока массы вещества т2 при различных значениях угла в1. Случай в1 = 5° обозначен сплошной линией, в1 = 25° — штриховой, в1 = 45° — пунктирной. Все графики в2(т2) для разных в1 оканчиваются на кривой, нанесенной тонкой линией

Выделим на рис. 2 четыре области, для каждой из которых характерно свое поведение зависимости 02(т2). В области I (0< т2<3) с ростом т2 происходит уменьшение составляющей Ву2 вектора В2 магнитного поля, причем 0 < 02 < 01, т. е. при переходе через поверхность разрыва тангенциальная составляющая магнитного поля ослабевает. В точке т2 = 3 при пересечении плоскости разрыва тангенциальная компонента Ву2 становится равной нулю. В области II (3< т2<4) компонента Ву2 отрицательна, увеличивается по модулю, но теперь — 01 < ^2<0. В области III (4 < т2 < 6), как и в области II, при пересечении плоскости разрыва тангенциальная компонента Ву меняет знак. Однако теперь Ву увеличивается по модулю (62 < — в1). Наконец, в области IV (т2 > 6) происходит усиление магнитного поля (02 > 01) с сохранением знака.

Таким образом, мы видим, как при постепенном увеличении потока плазмы меняется характер соотношения между углами наклона магнитного поля, а следовательно, и тип МГД-разрыва. Области I и II отвечают медленным ударным МГД-волнам, соответственно не обращающим (51) и обращающим (51) тангенциальную компоненту магнитного поля. Область III

соответствует трансальвеновской ударной волне ( Tr ), а область IV — быстрой ударной волне (5+ ).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

На границе областей II и III m2 = m|. Область определения функции 02(m2, 91) слева от величины m| задается условием (12), а справа — условием (13). Так как при m2 ^ m| величина By ^ 0, условие (13) преобразуется к виду m2 > m|. Поэтому функция 02(m2, 91) как в области II, так и в области III определена вблизи значения m| (см. рис. 2). Однако при увеличении m2 растет и правая часть условия (13). При некотором значении m2 в области III эти величины сравниваются. Наблюдается максимально сильная (в наибольшей степени усиливающая энергию магнитного поля) трансальвеновская ударная волна. При дальнейшем увеличении потока массы значение m2 не может удовлетворить условиям (12)-(13) пока m2 вновь не сравняется с m| + m\ . Происходит это в области IV, где наблюдается наиболее сильная быстрая ударная волна.

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

B2 + B2

2 x ~ ^y

m2 =--

4nr

Выразим величину By1 :

By1 = ±2\! 4nrm2 - B2 - By2,

y2,

где символ «+» соответствует области IV, а « — » — области III. Разделим левую и правую части полученного уравнения на Bx :

tg91 = ±2^m2/m2 - 1 - tg

(14)

Подставив выражение (14) в (9), получим уравнение искомой кривой:

tg 92 = ±

2 m2/m2ff - 2

m2/ m2ff + m2/ m2n - 2

/

m2 1.

mA

(15)

График кривой (15) представлен на рис. 2 тонкой линией.

3. Скачок внутренней энергии

Чтобы определить эффективность нагрева плазмы, обратимся к граничному условию (8), представляющему собой закон сохранения энергии. Используя уравнение (2), найдем из (8) скачок внутренней энергии

{е} = — {т}—т{ихР}—4пт

{B2üx - (v ■ B)Bx} . (16)

Правая часть уравнения (16) состоит из трех слагаемых. Используя средние величины скоростей г>х, г>у и г>г, представим первое слагаемое в виде

— { у} = —К }— ^у{иу }— ЪгШ.

Скачки тангенциальных компонент скорости выразим через скачки тангенциальных компонент магнитного

поля при помощи уравнений (5) и (6) вместе с условиями (1)-(2). Получим

= {ВУ}, = {Вг}.

Теперь первое слагаемое в правой части уравнения (16) выглядит следующим образом:

2

\ V21 ~ г 1 ууВх <0 л VzBх

-( 2] = -vх{vх}- 4ут

В}- {Вг}. (17)

4пт

Во второе слагаемое, записанное в виде 1 г 1 Р г 1 V г

--{VхP} =--{Vх }--{р},

т т т

подставим скачок давления из уравнения (7), а именно

{р} = -тК}- В{Ву}- 4П{Вг}.

Здесь так же, как и в (16), использовано условие (2). Второе слагаемое уравнения (16) в итоге принимает вид

1 Г 1 Р Г 1 , ~ Г 1 , ^ВУ ГО 1 , {в л

--{^р} = — {vх}+^ {vх} + ^птт {Ву} + {Вг}.

т

т

(18)

В третьем слагаемом раскроем скалярное произведение (V ■ В):

1

4пт

{B2Vх - (V ■ В)Вх} =

1

4пт

{(VхВу - VyBх)Ву + (VхВг - VzBх)Bz} .

Затем применим к полученному уравнению условия (3) и (4):

1

4пт

{В^х - (V ■ В)Вх} =

щВУл - ^Вх {Ву}- т. (19)

4пт { у' 4пт {

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

г 1 Р г 1 , ^хВу - УуВ^ VхBy - VyBх

{е} = - т + 4-т {Вл} —4- {Вл} +

+ УхВг - УгВх в 1 _ VхBz - VzBх В |

4пт г 4пт г '

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

{е} - -т<*}-Щт-{В-}

Вынося -{Vх}/т = -{г} за скобки, находим окончательно уравнение, выражающее скачок внутренней энергии на разрыве через скачки обратной плотности и тангенциальных компонент магнитного поля:

{е} = -{г} р +

{Ву }2 + {В-}2\ 16п ) .

(20)

Для двумерных разрывов соотношение (20) принимает совсем простой вид

{Ву}2 '

{е} = -{г} [р +

(р+щ

(21)

4. Условия вблизи области пересоединения

Уравнение (20) позволяет сделать определенные выводы относительно изменения внутренней энергии плазмы при пересечении поверхности разрыва. Во-первых, внутренняя энергия плазмы растет, так как -{г} > 0 по теореме Цемплена, а р и {Ву}2 — положительные величины. Во-вторых, изменение внутренней энергии складывается из двух частей: термодинамической, определяемой уравнением состояния вещества, и магнитной, связанной с изменением структуры магнитного поля вблизи поверхности разрыва.

Рассмотрим адиабатическое приближение р ~ р7, где 7 — показатель адиабаты. Термодинамическая часть скачка внутренней энергии увеличивается с ростом скачка плотности на разрыве: -{г}р ~ {р}7. Зависимость от р1 определяется величиной показателя 7, а именно

{е} = -{г }р ~ р1-2.

Для идеального газа 7 = (г + 2)/г, где г — число степеней свободы частиц плазмы (см. [22]). Скачок внутренней энергии увеличивается с ростом р1 при 7 >2 (г <2) и уменьшается при 7 <2 (г >2). Ограничение числа степеней свободы, таким образом, влияет на характер нагрева вещества.

Магнитная часть скачка внутренней энергии зависит от конфигурации магнитного поля, а значит, и от типа разрыва. Для расчетов, представленных на рис. 2, зависимости скачка внутренней энергии от потока массы через разрыв, вычисленные с использованием формулы (21), показаны на рис. 3. Здесь за нуль принято значение термодинамической части скачка, не зависящей от потока массы. Видно, что максимальный скачок внутренней энергии при заданных параметрах плазмы осуществляется наиболее сильной трансальве-новской ударной волной, причем его величина быстро растет с увеличением угла падения магнитного поля 01 . Соотношения между эффективностью нагрева плазмы другими типами разрывов зависят от конкретных условий среды. Так, нагрев медленными ударными волнами

{е} 25

20

15

10

5

I !п 1 ;1П 1 IV

: 1 ; 1 ; 1 / 1 ; /I / ......г—!' 1/ V .......

^.....

0

10 12 14 16 18 ти

Рис. 3. Зависимости скачка внутренней энергии {е} от потока массы вещества при различных значениях угла 01: 5° — сплошная линия, 25° — штриховая, 45° — пунктирная

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

Наибольшего разогрева на разрывах можно ожидать от плазмы, претерпевающей большой скачок плотности в магнитном поле с резко меняющейся геометрией. Именно такие условия выполняются в области магнитного пересоединения. При слиянии двух противоположно направленных магнитных потоков образуется токовый слой, который замедляет процесс пересоединения и накапливает в себе свободную магнитную энергию. В дальнейшем эта энергия реализуется при разрыве токового слоя (см., например, [23]). Магнитное пересоединение осуществляет быстрое преобразование свободной энергии в энергию частиц плазмы и сопровождается образованием сложной картины МГД-раз-рывов в областях с резкими изменениями магнитного поля и поля скоростей. Такие системы разрывов наблюдаются как в лабораторных, так и в численных экспериментах [4, 6].

В то же время подобные разрывные структуры могут быть внесены в качестве самостоятельного элемента в ту или иную модель пересоединения. На рис. 4 представлена токовая система самосогласованной аналитической модели магнитного пересоединения [24, 25]. К концу токового слоя (СЬ) присоединены две ударные волны конечной длины (на рис. 4 показана только правая половина токового слоя, левая достраивается симметрично). Вещество втекает в слой сверху и снизу, вытекает слева и справа. Ранее нами было показано, что вблизи торцов токового слоя, где образуются обратные токи, ударные волны являются трансальве-новскими [15]. По мере удаления от токового слоя скачок напряженности магнитного поля на разрыве уменьшается, как уменьшается до нуля и скачок плотности. Таким образом, наилучшие условия для нагрева плазмы осуществляются вблизи области обратных токов.

Рис. 4. Конфигурация электрических токов (жирные прямые отрезки) состоит из токового слоя (О) и четырех присоединенных к его торцам поверхностей разрыва конечной длины R. L — полуширина токового слоя

Заключение

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

Работа выполнена при финансовой поддержке РФФИ (грант 11-02-00843-а)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Список литературы

1. Sutton G.W., Sherman A. // Engineering Magnetohydro-dynamics. N.Y,; San Francisco; Toronto, 1965.

2. Лукьянов С.Ю. // Горячая плазма и управляемый ядерный синтез. М., 1975.

3. Морозов А.И. // Введение в плазмодинамику. М., 2006.

4. Biskamp D. // Nonlinear Magnetohydrodynamics. Cambridge, UK, 1997.

5. Имшенник В.С., Боброва Н.А. // Динамика столкнови-тельной плазмы. М., 1997.

6. Space Plasma Simulation / Ed. by J. Buchner, C. T. Dum, M. Scholer. Berlin; Heidelberg; N.Y., 2003.

7. Petschek H.E. // AAS-NASA Symposium on the Physics of Solar Flares. NASA SP-50 / Ed. by W.N. Hess. Washington, 1964. P. 425.

8. Брушлинский К.В., Заборов А.М., Сыроватский С.И. // Физика плазмы. 1980. 6. С. 297.

9. Hones E.W.Jr. (Ed.) // Magnetic Reconnection in Space and Laboratory Plasmas. Washington, DC, 1984.

10. Magnetic Reconnection in Space and Laboratory Plasmas / Ed. by M. Hoshino, R. L. Stenzel, K. Shibata. Tokyo, 2001.

11. Orta J.A., Huerta M.A., Boynton G.C. // Astrophys. J. 2003. 596. P. 646.

12. Ландау Л.Д., Лифшиц Е.М. Курс теоретической физики. Т. 6. Гидродинамика. М., 1986.

13. Сыроватский С.И. // УФН. 1957. 62. С. 247.

14. Anderson J.E. // Magnetohydrodynamic Shock Waves. Massachusetts (Cambridge, US), 1963.

15. Леденцов Л.С., Сомов Б.В. // Письма в Астрон. журн. 2011. 37. С. 151.

16. Tideman D.A. and Krall N.A. // Shock Waves in Collision-less Plasma. N.Y.; L.; Sydney, 1971.

17. Сиротина Е.П., Сыроватский С.И. // ЖЭТФ. 1960 39. С. 746.

18. Зельдович Я.Б., Райзер Ю.П. // Физика ударных волн и высокотемпературных гидродинамических явлений. М.,

1966.

19. Шерклиф Дж. // Курс магнитной гидродинамики. М.,

1967.

20. Ландау Л.Д., Лифшиц Е.М. // Курс теоретической физики. Т. 8. Электродинамика сплошных сред. М., 1982.

21. Somov B.V. // Plasma Astrophysics. Pt I. Fundamental and Practice. N.Y., 2012.

22. Ландау Л.Д., Лифшиц Е.М. // Курс теоретической физики. Т. 5. Статистическая физика. Ч. 1. М., 1976.

23. Somov B.V. // Plasma Astrophysics. Pt II. Reconnection and Flares. N.Y., 2012.

24. Безродных С.И., Власов, В.И., Сомов Б.В. // Письма в Астрон. журн. 2007. 33. C. 153.

25. Безродных С.И., Власов, В.И., Сомов Б.В. // Письма в Астрон. журн. 2011. 37. C. 133.

Plasma heating by discontinuous MHD flows near the region of magnetic reconnection L.S. Ledentsova, B.V. Somovb

P. K. Sternberg State Institute of Astronomy, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.

E-mail: a [email protected], b [email protected].

We had used a complete system of boundary conditions for MHD equations on a surface of discontinuity to obtain an expression that describes the change in internal energy of plasma flowing through the discontinuity. The dependence of plasma heating from the density and the configuration of magnetic field near the surface of discontinuity (in other words, from the type of MHD flow) are studied. Conditions of plasma heating at discontinuities in the self-consistent analytical model of magnetic reconnection are discussed.

Keywords: plasma, magnetohydrodynamics, discontinuities. PACS: 96.60.Q-, 96.60.qe, 96.60.Iv, 95.30.Qd, 96.50.Ry.

Received 2 December 2012.

English version: Moscow University Physics Bulletin 2(2013).

Сведения об авторах

1. Леденцов Леонид Сергеевич — инженер I категории; тел.: (495) 939-16-44, e-mail: [email protected].

2. Сомов Борис Всеволодович — докт. физ.-мат. наук, профессор, зав. отделом; тел.: (495) 939-16-44, e-mail: [email protected].

i Надоели баннеры? Вы всегда можете отключить рекламу.