Численное моделирование внезапного изменения топологии магнитного поля

Вычислительные технологии

Том 2, № 6, 1997

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВНЕЗАПНОГО ИЗМЕНЕНИЯ ТОПОЛОГИИ МАГНИТНОГО ПОЛЯ

Г. И. ДудниковА, В. П. Жуков, Т. В. ЛисЕйкинА Институт вычислительных технологий СО РАН, Новосибирск, Россия

e-mail: dudn@adm.ict.nsc.ru

The influence of constant background magnetic field directed along hiperbolic null line on the reconnection process initiated by magnetosonic wave in the Earth's magnetosphere is discussed.

1. Введение

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

В работе исследуется влияние однородного постоянного магнитного поля, направленного вдоль гиперболической нулевой линии, на течение плазмы в ее окрестности. Распространение магнитозвуковых волн в отсутствие тороидального магнитного поля подробно описано в [2]. Сильное тороидальное магнитное поле позволяет разложить исходные уравнения по малому параметру: отношению полоидального магнитного поля к тороидальному. В результате получается хорошо известная в теории пилообразных колебаний в токамаке модель Кадомцева [3]. Согласно этой модели различные величины в задаче определенным образом зависят от упомянутого параметра. Основной целью данной работы является проверка этих зависимостей для задачи о распространении магнитозвуковых волн в окрестности -линии магнитного поля. Несмотря на то, что данная задача имеет существенное значение для понимания пределов применимости широко используемой модели Кадомцева, она ранее не исследовалась систематически.

2. Постановка задачи

Математическая модель, приведенная в [1, 4], представляет собой МГД уравнения с учетом эффекта Холла в двумерном (d/dz = 0) плоском случае. Задача решается в квадратной

© Г. И. Дудникова, В. П. Жуков, Т. В. Лисейкина, 1997.

32

Г. И. Дудникова, В. П. Жуков, Т. В. Лисейкина

области: О = {—х0 < х < х0, —у0 < у < , х0 = у0 = 1. Начальные условия имеют вид Их = И* =const — г-компонента магнитного поля, А = 2 (х2 — у2) — г-компонента векторного потенциала, р = в =const — давление, р =1 =const — плотность, V = 0 — скорость плазмы,

И± = (0Аг/ду, —ЗАХ/дх) — поперечные нулевой линии компоненты магнитного поля (х = 0, у = 0 — нулевая линия). Плазма выводится из равновесия при помощи возмущения вектор-потенциала Аг, задаваемого на границе области в виде цилиндрической волны

ное значение электрического поля на границе. Граничные условия допускают протекание плазмы через границы расчетной области [2].

3. Результаты численного моделирования

При И* = 0 течение хорошо изучено [2]. Рассмотрим влияние величины И* на картину течения в отсутствие эффекта Холла. Прежде всего отметим, что если при И* = 0 плотность плазмы изменяется на величину 5р порядка 1, то при увеличении И* изменение плотности стремится к 0 (приближение несжимаемой жидкости). Это обстоятельство широко используется при создании упрощенных моделей [3]. Расчеты показывают, что зависимость 5р от И* выглядит следующим образом:

откуда ёр ~ И—1, что соответствует [3].

Распределение -компоненты электрического тока на квазистационарной стадии процесса (I > 6) (рис. 1) мало меняется по сравнению со случаем И* = 0. Как и в отсутствие продольного магнитного поля (И*), образующийся токовый слой располагается вблизи се-паратрисных поверхностей.

характер-

И* = 2, ёр = 0.236, И* = 5, ёр = 0.12, И* = 10, ёр = 0.077,

\

Рис. 1. Распределение тока (Н* = 10.0).

Рис. 2. Распределение (Н* = 2.0).

Распределение квадрата поперечной компоненты электрического тока ^ = (дНг/ду, —дНг/дх) показано на рис. 2 (при Н* = 0 эта компонента отсутствует). Областью локализации токового слоя является окрестность сепаратрисных поверхностей. С увеличением Н* уменьшается максимальное значение величины ^: Н* = 2, тах(?| ) = 0.178, Н* = 5, тах(?1) = 0.042, Н* = 10, тах(^1) = 0.012. То есть тах(^) ~ Н—2 в соответствии с [3].

Наличие продольного магнитного поля (Н*) существенно влияет на распределение давления плазмы (рис. 3). Видно, что имеет максимум в месте расположения токового слоя (рис. 3, б, в), а не минимум, как в случае отсутствия Н* (рис. 3, а). При увеличении Н* этот максимум тоже увеличивается, но незначительно.

На рис. 4 изображена поверхность продольной компоненты магнитного поля Нг. Видно, что с ростом Н* уменьшается максимальное значение величины 5 = тах(Нг — Н*)/Н*: Н* = 2, 5 = 0.0375, Н* = 5, 5 = 0.0096, Н* = 10, 5 = 0.0026.

Величина Р^ = Р + Н /2 — Н*2 / 2 мало зависит от значения Н*, а ее распределение (рис. 5) сходно с распределением давления плазмы в случае Н* = 0 (см. рис. 3, а).

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

В заключение отметим, что наличие продольного магнитного поля существенно влияет на распределение плотности плазмы, ее давления, а также на распределение квадрата поперечной компоненты тока и оказывает малое влияние на распределение р + Н| / 2, тока и на поле скоростей плазмы. Образующиеся токовые слои располагаются в окрестности сепаратрисных поверхностей. Условия применимости упрощенной модели [3] весьма хорошо выполняются уже при достаточно небольших значениях однородного магнитного поля Н*.

Рис. 3. Распределение давления плазмы: Н* =0 (а), 2.0 (б), 10.0 (в).

Рис. 4. Распределение величины Нг — Н* (Н* = 2.0).

Рис. 5. Распределение Р^ = р + Н^ — Н*2 (Н* = 5.0).

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

[1] Буланов С. В., ДудниковА Г. И., ИнновЕнков И.Н., Пичушкин В. В. Токовые слои в окрестности критических точек магнитного поля. Влияние эффекта Холла. Численное моделирование формирования токовых слоев. Тр. ИОФАН, 51, 1995, 101-123.

[2] Брушлинский К. В., Заборов А. М., СыровАтский С. И. Численный анализ токового слоя в окрестности магнитной нулевой линии. Физика плазмы, 6, №8, 1986, 297-311.

[3] Кадомцев Б. Б., Погуце О. П. Нелинейные спиральные возмущения плазмы в то-камаке. ЖЭТФ, 65, вып. 2(8), 1973, 575-589.

[4] Буланов С. В., ДудниковА Г. И., Жуков В. П. и др. Пересоединение магнитных силовых линий в окрестности критических точек. Физика плазмы, 22, №10, 1996, 867-895.

Поступила в редакцию 8 ноября 1995 г., в переработанном виде 11 июня 1997 г.