Научная статья на тему 'Влияние эффектов вязкости на поведение изменения энтропии в ударной волне, имевшей место после коронального выброса массы 13 декабря 2006 года'

Влияние эффектов вязкости на поведение изменения энтропии в ударной волне, имевшей место после коронального выброса массы 13 декабря 2006 года Текст научной статьи по специальности «Математика»

CC BY
152
84
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Астрофизический бюллетень
WOS
Scopus
ВАК
Область наук
Ключевые слова
солнечный ветер — Солнце: корональные выбросы массы (CME) — Солнце: активность / solar wind—Sun: coronal mass ejections (CMEs)—Sun: activity

Аннотация научной статьи по математике, автор научной работы — Х. Чавуш, А. Курт

Корональные выбросы массы (CME) и солнечный ветер — два основных проявления солнечной активности. Эти явления могут вызвать ударную волну. Ударная волна возникает там, где солнечный ветер меняется со сверхзвукового (по отношению к окружающей межпланетной среде) на дозвуковой. Главной целью данной работы является применение алгоритма и результатов, приведенных в наших недавних публикациях, к ударной волне, возникшей после CME 13 декабря 2006 г., и изучение поведения энтропии во время этого проявления солнечной активности.

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

Effects of Viscosity on the Behavior of Entropy Change in the ShockWave That Occurred after the December 13, 2006 CoronalMass Ejection

Coronalmass ejections (CMEs) and the solar wind are the two main demonstrations of solar activity. These events can drive a shock wave. The shock wave occurs where the solar wind changes from being supersonic (with respect to the surrounding interplanetary medium) to subsonic. The main purpose of this study is to apply the algorithm and the results given in our recent papers to the shock wave that happened after the December 13, 2006 CME, and evaluate the behavior of entropy during this solar activity.

Текст научной работы на тему «Влияние эффектов вязкости на поведение изменения энтропии в ударной волне, имевшей место после коронального выброса массы 13 декабря 2006 года»

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2015, том 70, № 2, с. 231-236

УДК 523.98-782; 523.62-726

ВЛИЯНИЕ ЭФФЕКТОВ ВЯЗКОСТИ НА ПОВЕДЕНИЕ ИЗМЕНЕНИЯ ЭНТРОПИИ В УДАРНОЙ ВОЛНЕ, ИМЕВШЕЙ МЕСТО ПОСЛЕ КОРОНАЛЬНОГО ВЫБРОСА МАССЫ 13 ДЕКАБРЯ 2006 ГОДА

© 2015 Х. Чавуш*, А. Курт

Университет 18 марта Чанаккале, Чанаккале, 17100 Турция

Поступила в редакцию 15 октября 2014 года; принята в печать 29 декабря 2014 года

Корональные выбросы массы (CME) и солнечный ветер — два основных проявления солнечной активности. Эти явления могут вызвать ударную волну. Ударная волна возникает там, где солнечный ветер меняется со сверхзвукового (по отношению к окружающей межпланетной среде) на дозвуковой. Главной целью данной работы является применение алгоритма и результатов, приведенных в наших недавних публикациях, к ударной волне, возникшей после CME 13 декабря 2006 г., и изучение поведения энтропии во время этого проявления солнечной активности.

Ключевые слова: солнечный ветер — Солнце: корональные выбросы массы (CME) — Солнце: активность

1. ВВЕДЕНИЕ

Корона является внешней частью солнечной атмосферы. Она характеризуется сложными плазменными и магнитными структурами. В короне возникают, главным образом, два структурных элемента: замкнутые и открытые магнитные структуры. Тип активности характеризуется взаимодействием между плазмой и магнитным полем в солнечной короне [1]. Корональные выбросы массы (coronal mass ejections, CME) происходят из-за замкнутой магнитной петли, которая порождает огромное облако плазмы в межпланетном пространстве [2]. Быстрый поток плазмы (т.е. солнечный ветер) вызван незамкнутыми структурами магнитных силовых линий [3]. В качестве источников медленного солнечного ветра недавно также были предложены струи по краям корональных дыр [4]. Из-за взаимодействий с локальной межзвездной средой сверхзвуковые потоки частиц в солнечном ветре могут вызвать ударную волну.

Существует несколько способов вызвать ударную волну в солнечном ветре, например, взрывные волны от Солнца, CME, а также взаимодействия между быстрыми и медленными потоками [5]. Ударная волна возникает вследствие того, что частицы солнечного ветра испускаются Солнцем со скоростями примерно 550—700 км с-1, а скорость звука (в межпланетной среде) составляет приблизительно 100 км с-1 [6, 7]. Однако из-за

E-mail: h_cavus@comu.edu.tr

флуктуаций ее точное значение сильно зависит от плотности.

В своей работе Еселевич и Еселевич [8] показали, что перед фронтальной структурой CME постепенно формируется возмущенная область, это вызвано взаимодействием CME с невозмущенным солнечным ветром. Ее размер увеличивается по мере удаления CME от Солнца. В узкой области фронта возмущенной зоны наблюдается скачок. Характеристики этой возмущенной зоны удовлетворяют критериям поршневой ударной волны. Эта волна является столкновительной при радиальных расстояниях r < 6 R© от центра Солнца (Я© — радиус Солнца) и становится бесстолкновительной при r > 6 R©.

Ударные волны изучаются уже более ста лет. Однако, мало работ посвящено распределению энтропии по ударной волне. Поведение энтропии — сложный для понимания предмет. Мордучов и Либби [9] изучали распределение энтропии по слою ударной волны без учета эффектов вязкости и теплопередачи. Они нашли, что энтропия в отличие от других физических параметров, которые ведут себя монотонно, увеличивается во фронтовой части ударной волны до максимального значения в центре фронта, а затем уменьшается в другой половине фронта. Может показаться, что это противоречит второму закону термодинамики, однако авторы объясняют эту ситуацию следующим образом: второй закон действует для системы в целом, т.к. энтропия увеличивается в прифронтовой области ударной волны.

231

ЧАВУШ, КУРТ

Таблица 1. Значения физических параметров за фронтом ударной волны, возникшей после CME 13 декабря 2006 г. [13]

_я 7li, СМ Pi, пРа Д,К mi, км с 1 /Я

1.8 0.982 1.5 х 102 * * 5 * 573 1.17 х 107

Хамад [10] изучал поведение энтропии в ударных волнах в обычном пылевом газе, используя уравнения Навье—Стокса. Он нашел, что максимум профиля энтропии приходится на фронт ударной волны. Он также показал, что энтропия увеличивается по волне в зависимости от числа Маха за фронтом волны и плотности частиц. Хамад [11] изучал влияние числа Прандтля на поведение энтропии в обычной ударной волне. Он представил решение для простого газа и нашел, что энтропия уменьшается по мере увеличения вязкости и, наоборот, увеличивается с возрастающими значениями теплопроводности и зафронтового числа Маха.

Вандегрифф и др. [12] разработали модель для предсказания прихода ударных волн из межпланетного пространства на Землю. Лиу и др. [13] исследовали ударную волну, возникшую после CME 13 декабря 2006 г., используя данные STEREO, Ulysses и ACE. Чавуш [14, 15] (далее работы I и II) исследовал эффекты вязкости, разработав новый алгоритм. В работе I описываются механизм и формирование ударной волны. В работе II модель была применена к CME 13 декабря 2006 г., чтобы найти значения физических параметров прифронтовой области ударной волны. Основная цель данной работы, в отличие от предыдущих, — исследовать поведение энтропии для ударной волны, возникшей в солнечном ветре после CME 13 декабря 2006 г. Чтобы описать такой тип ударного процесса, нужно решить уравнения Навье—Стокса. Математически этот процесс можно аппроксимировать гидродинамическим случаем. При такой процедуре вязкое поведение жидкости можно описать, используя число Рейнольдса, как в работах I и II. Результаты сравниваются с другими подобными работами в разделе 4, вместе с обсуждением и заключением.

2. ФИЗИЧЕСКАЯ ФОРМУЛИРОВКА ЗАДАЧИ

2.1. Основные физические свойства

Физическая структура солнечной атмосферы

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

магнитное и плазменное давления играют важные роли, попеременно доминируя. Доминирование этих сил определяется параметром плазмы Д

Этот коэффициент принимает разные значения в зависимости от магнитного поля. Гэри [16]

утверждал, что он имеет большие значения (в > 1) в области ускорения солнечного ветра, теоретически равные бесконечности. Маттеус и др. в своей работе [17] приводят его значение от 44 до бесконечности. Следовательно, можно с легкостью утверждать, что давление газа в динамике жидкой или газообразной среды играет решающую роль в динамике солнечного ветра [18].

CME можно идентифицировать по профилям плотности, давления, температуры и скорости. Мы используем ударную волну, возникшую после CME 13 декабря, 2006 г. Соответствующие параметры за фронтом этой волны можно получить из данных миссии ACE, они приведены в таблице 1 [13]. Этот CME — один из самых больших со времен Хэллоуинской бури 2003 г. Цель данной работы — вычислить поведение энтропии, применяя алгоритм и результаты приведенные в работах I и II для ударной волны, вызванной CME 13 декабря 2006 г.

Как было указано выше, локальная скорость звука в межпланетном пространстве равна приблизительно 100 км с-1 [7], следовательно, в этой области может возникнуть ударная волна со значениями, приведенными в таблице 1.

2.2. Формулировка энтропии

Чтобы вычислить поведение энтропии по ударной волне, рассмотрим уравнение энтропии [19],

T5S = 5e + р5(1/р), (1)

где T, 5, e, p и р обозначают соответственно температуру, малые изменения, внутреннюю энергию, давление и плотность. Внутреннюю энергию можно определить как e = h — p/р, используя энтальпию h. Для идеального газа уравнение (1) можно переписать в виде

6S = cv6\n^ - R6\np. (2)

р

cv = R/(y — 1) и R обозначают соответственно удельную теплоемкость при постоянном объеме и газовую постоянную. Уравнение (2) можно переписать, пользуясь определением R = cp — cv, через удельную теплоемкость при постоянном давлении (т.е. cp = y R/(y — 1)) следующим образом:

5S = cv 5 lnp — cP 5 ln р. (3)

Последнее уравнение определяет S с точностью до константы, которая подбирается таким образом, что при cv = const, уравнение (3) становится

S = cv ln(p/pY). (4)

В последнем уравнении y обозначает отношение удельных теплоемкостей (cp/cv), ее значение равно 5/3 для одноатомного газа. Значение энтропии за

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 2 2015

ВЛИЯНИЕ ЭФФЕКТОВ ВЯЗКОСТИ НА ПОВЕДЕНИЕ ИЗМЕНЕНИЯ ЭНТРОПИИ

233

Рис. 1. Зависимости числа Рейнольдса Re2 от М1 (слева) и Rei (справа).

фронтом ударной волны, возникшей после CME 13 декабря 2006 г., можно вычислить, используя физические параметры, приведенные в таблице 1.

2.3. Общее уравнение Значения параметра плазмы в

(3 = ^- (5)

pmag

больше в солнечном ветре [16—18, 20]. Следовательно, газовое давление играет более важную роль, чем магнитное поле. Солнечный ветер может приводиться в движение в основном газовым давлением при высоких температурах. Задача становится схожей с гидродинамическим случаем, приведенным в работах Паркера [21,22] и Холзера и Аксфорда [23].

Физическая формулировка настоящей задачи будет структурирована, как в работе I. Для сжимаемой вязкой ударной волны в равномерном течении, общее уравнение для коэффициента сжатия (т.е. отношение плотностей за и перед фронтом) было получено в виде (работа I)

4 1 3 Rei

(y -1) m2 + i

к

2

4 1 1 —

3 Re

у m2 +1

к

(6)

+

Y + 1 2

4 1 3 Re2

Mf

0.

Индексы 1 и 2 обозначают области за и перед фронтом соответственно, Rei и Re2 — зафронтовое и прифронтовое числа Рейнольдса, к и Mi — коэффициент сжатия и зафронтовое число Маха. Величины Mi, Rei и Re2 имеют важное значение для распределений физических параметров перед фронтом. Общее уравнение можно решить численно, вводя подходящие условия. Согласно работе Еселевич и Еселевич [24], у удовлетворяет соотношению 5/3 < у < 3 для бесстолкновительного

фронта ударной волны. Для нахождения значений физических параметров перед фронтом следует использовать уравнение (6) и соотношения скачка Рэнкина—Гюгонио [25] . Используя уравнение (4), изменение энтропии можно выразить через отношение давлений и коэффициент сжатия следующим образом:

S2 - Si

Су ln

К-к-у

Pi

(7)

2.4. Значение числа Рейнольдса перед фронтом в солнечном ветре

Число Рейнольдса (Re) играет важную роль в динамике этой области. Обычно, оно имеет большие значения (Re > 1) в солнечном ветре. Re принимает значения между 10i2 и 10i4 в области ускорения солнечного ветра [4, 26]. Для упрощения задачи число Рейнольдса перед фронтом Re2 было выражено как функция адиабатического индекса, Rei и Mi, как в работе Бруна и др. [27] и в работе I. Слева на рис. 1 отношение Re2/Re1 приводится как функция Mi. Оно уменьшается по мере увеличения числа Маха за фронтом, а отношение чисел Рейнольдса равно единице (т.е. Rei = Re2) при Mi ~ 2, как показано на рисунке. Это значение Mi считается точкой перехода от слабых ударных волн к сильным [14, 28].

Изменение Re2 относительно Re1 для разных значений M1 для одноатомного газа дано справа на рис. 1. Значения числа Рейнольдса за фронтом взяты из литературы: 10i2—10i4 для области солнечного ветра [4, 26]. Используя эти значения, мы показываем, что Re2 увеличивается с Rei. Еще один важный результат, полученный из этого графика, состоит в том, что Re2 больше для меньших значений M1.

3. ПОВЕДЕНИЕ ЭНТРОПИИ ДЛЯ УДАРНОЙ ВОЛНЫ ПОСЛЕ CME 13 ДЕКАБРЯ 2006 Г.

Мы попытались применить модель для нахождения параметров прифронтовой области ударной

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №2 2015

234

ЧАВУШ, КУРТ

Рис. 2. Распределение разности энтропии S2 — Si в зависимости от Mi.

Рис. 3. Зависимость разности энтропии S2 — Si от Re2/Rei.

9.0

7.5

6.0

4.5 3.0

1.5

■ Cs2 /cSi * К

к Re2/Rei о ц/ц - M2 /Mi

.. ■ *4 ‘

И*;:'

-

44444444

3.4

6.8 10.2 S2- Si

13.6

17.0

Рис. 4. Вариации некоторых параметров по отношению к разности энтропии S2 — Si.

волны, возникшей после CME 13 декабря 2006 г., используя некоторые специальные решения уравнений (6) и (7) и соотношения скачка Рэнкина— Гюгонио [25] и применяя метод, приведенный в работе I, и результаты из работы II. Метод реализован в виде алгоритма для системы Maple 9.5. Результаты для поведения энтропии представлены в таблице 2 и на рис. 2—5. В этих вычислениях для солнечного ветра [4] было использовано значение зафронтового числа Рейнольдса 1013.

В таблице 2 представлены вариации энтропии для разных значений чисел Маха и Рейнольдса за фронтом. Распределение S2 — S1 уменьшается с возрастающими значениями отношения чисел Рейнольдса. Для более высоких значений Re1, S2 — S1, как и ожидалось, больше. Критическое значение разности энтропии было найдено равным 2.25 для точки поворота М1 = 2.045, которая приравнивает отношение чисел Рейнольдса к единице, как показано на рис. 1.

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

На рис. 2 показана зависимость изменения энтропии S2 — S1 от М1. Видно, что изменение энтропии отрицательно для меньших значений М1 (т.е. М1 < 1.2). Этот случай нарушает второй закон термодинамики (т.е. S2 — S1 < 0). Таким образом, слабые ударные волны (т.е. М1 < 1.2) почти изо-энтропийны. С другой стороны, на рис. 3, оно

16

12

со

оУ 8 4 0

7

Рис. 5. Изменение разности энтропии в зависимости от числа Маха за фронтом: (/) Mi = 1.2,

Re2/Rei = 1.704; (2) Mi = 1.6, Re2/Rei = 1.278;

(4) Mi = 2, Re2/Rei = 1.022; (4) Mi = 2.045, Re2/Rei = 1; (5) Mi = 2.5, Re2/Rei = 0.818;

(6) Mi =3, Re2/Rei = 0.682; (7) Mi =4,

Re2/Rei = 0.511; (8) Mi = 5, Re2/Rei = 0.409.

уменьшается с возрастающим отношением чисел Рейнольдса. Ударная волна становится изоэнтропийной для возрастающих значений Re2/Re1.

На рис. 4 показаны вариации отношений некоторых параметров в зависимости от S2 — S1. Заполненные квадраты показывают отношение скоростей звука в областях перед и за фронтом. Оно возрастает с разностью энтропии. Треугольники представляют зависимость коэффициента сжатия. Он возрастает с ростом S2 — S1, как и отношение скоростей звука. Однако отношение М2/М1, показанное штриховой линией, уменьшается. Отношение скоростей перед и за фронтом показано кружками. Оно уменьшается с возрастанием разности энтропии. Как и отношение скоростей, отношение чисел Рейнольдса также обратно пропорционально разности энтропии. Все эти соотношения равны единице для изоэнтропийного случая (т.е. S2 — S1 = 0). Другими словами, при S2 — S1 = 0 ударной волны не возникает, т.к. к = 1 (т.е. отсутствует сжатие).

На рис. 5 представлена зависимость разности

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 2 2015

ВЛИЯНИЕ ЭФФЕКТОВ ВЯЗКОСТИ НА ПОВЕДЕНИЕ ИЗМЕНЕНИЯ ЭНТРОПИИ

235

Таблица 2. Изменение энтропии для различных значений Mi и Re2/Re1

Ml = 1.2 Mi = 1.6 Mi = 2.045 Mi = 2.5 Mi = 4 Mi = 5

Re2/Rei = 1.704 Re2/Rei = 1.278 Re2/Rei = 1 Re2/Rei =0.818 Re2/Rei = 0.511 Re2/Rei =0.409

Со to 1 Со 0.06 0.80 2.25 4.57 11.97 16.37

энтропии относительно М1 (или отношения чисел Рейнольдса). Она возрастает с М1 подобно рис. 2. Очень слабые ударные волны (т.е. Mi < 1.2) и здесь почти изоэнтропийны. S2 очень близка к зафронтовому значению. С другой стороны, для сильных ударных волн (M1 > 2) разница большая. Она уменьшается с возрастающим отношением чисел Рейнольдса, как на рис. 3.

4. ОБСУЖДЕНИЕ И ЗАКЛЮЧЕНИЕ

Основная цель данной работы — применить алгоритм и результаты из работ I и II к ударной волне, возникшей во время CME 13 декабря 2006 г. для оценки поведения энтропии во время данного эпизода солнечной активности. CME — основное проявление солнечной активности. Эти явления могут не только вызывать геомагнитные бури, но и приводить в действие ударные волны в межпланетном пространстве. Изучение этих явлений очень важно для исследования космической погоды. С современными наблюдательными возможностями понимание ключевого вопроса, а именно, эволюции физических параметров, остается очень сложной задачей. Когда происходит корональный выброс массы, помимо взаимодействия с окружающим газом происходит сложный процесс с участием магнитной и тепловой энергии. Около Солнца доминирует магнитное давление. Давление газа становится преобладающим вдали от Солнца. Следовательно, к исследованиям ударных волн, вызванных CME, можно применять гидродинамическую модель [16—18]. Мы проанализировали поведение энтропии для ударной волны, возникшей 13 декабря 2006 г., применяя численную одномерную гидродинамическую модель. Выводы, полученные из модели, представлены по пунктам. •

• Значения числа Рейнольдса перед фронтом больше для меньших M1.

• Для зависимости числа Маха за фронтом (M1) разность энтропии S2 — S1 растет с M1. Для очень слабых ударных волн с M1 < 1.2 волна становится почти изоэнтропийной.

• Для вариаций энтропии в зависимости от Re2/Re1 ударная волна становится изоэнтропийной при возрастающих значениях отношения чисел Рейнольдса.

Таблица 3. Значения физических параметров перед фронтом ударной волны после CME 13 декабря 2006 г.

Mi Re2/Rei _Q 712, CM U2, км с -1 P2, Pa

4.5 0.45 6.27 131 2.4 x 10~8

• Для зависимости коэффициента сжатия к, S2 — S1 увеличивается с возрастанием к.

• Для зависимости отношения скоростей звука разность энтропии быстро растет с увеличением отношения скоростей звука.

• В отличие от вариаций к и отношения скоростей звука, S2 — S1 уменьшается с возрастанием числа Маха и отношения скоростей потока.

В работе II было найдено, что значения физических параметров прифронтовой области ударной волны (например, плотность приблизительно равна 6 см_3) соответствуют варианту М1 ~ 4.5 при сравнении результата с работой Лиу и др. [13]. Используя этот результат, мы привели полный список физических параметров прифронтовой области ударной волны в таблице 3.

При сранении поведение энтропии с работами Хамада [10—11] было найдено хорошее соответствие для зависимости скоростей. Он нашел похожую обратную пропорциональность по отношению к скорости. Масек и Редаелли [29] оценили энтропию для солнечного ветра без ударной волны, проанализировав данные космического аппарата Helios. Они подтверждают увеличение энтропии. В другой своей статье Редаелли и Масек [30] работали над экспонентой Ляпунова и энтропией солнечного ветра и нашли, что обе положительны. Ударное нагревание солнечного ветра и влияние неоднородности энтропии на корреляцию «плотность—температура» в солнечном ветре были изучены Вангом и др. [31] и Маларой и др. [32] соответственно. Они дают значение энтропии того же порядка, что и в настоящей работе. Было продемонстрировано, что модель представленная в данной работе, дает значения, близкие к литературным.

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №2 2015

236

ЧАВУШ, КУРТ

БЛАГОДАРНОСТИ

Автор с благодарностью отмечает техническую

и вычислительную поддержку Фонда научных

исследований Университета 18 марта Чанаккале

(проект № 2011/012).

СПИСОК ЛИТЕРАТУРЫ

1. E. R. Priest, Solar Magnetohydrodynamics, (D. Riedel Publ. Company, Dordrecht, 2000).

2. S. K. Antiochos, C. R. De Vore, and J. A. Klimchuk, Astrophys. J. 510,485(1999).

3. M. Stix, The Sun, (Springer, Berlin, 1991).

4. J. E. Borovsky and H. O. Funsten, J. Geophys. Research Space Physics 108, SMP 13-1 (2003).

5. P. A. Sturrock and J. R. Spreiter, J. Geophys. Research 70, 5345 (1965).

6. T. Suzuki, Monthly Notices Royal Astron. Soc. 349, 1227(2004).

7. V M. Nakariakov, L. Ofman, and T. D. Arber, Astron. and Astrophys. 353, 741 (2000).

8. V Eselevich and M. Eselevich, Astrophys. J. 761, 68 (2012).

9. M. Morduchow and P. A. Libby, J. Aeron. Sci. 16, 674 (1949).

10. H. Hamad, Zeitschrift fur angewandte Mathematik und Physik, 49,827(1998).

11. H. Hamad, Acta Mechanica 140, 65 (2000).

12. J. Vandegriff, K. Wagstaff, G. Ho, and J. Plauger, Advances Space Research 36, 2323 (2003).

13. Y. Liu, J. G. Luhmann, R. Muller-Mellin, et al., Astrophys. J. 689, 563 (2008).

14. H. Cavus, Advances in Astronomy 2013, 1 (2013).

15. H. Cavus, Astrophysical Bulletin 70, 117 (2015).

16. G. A. Gary, Solar Phys. 203, 71 (2001).

17. W. H. Matthaeus, S. Ghosh, S. Oughton, and D. A. Roberts, J. Geophys. Research 101, 7619 (1996).

18. J. A. Gonzalez-Esparza, P. Corona-Romero, and E. Aguilar-Rodriguez, in Proc. XXIX Int. Conf. on Phenomena in Ionized Gases (Cancun, Mexico, 2009).

19. P. Wesseling, Principles of Computational Fluid Dynamics (Springer, Berlin, 2001).

20. D. Tsiklauri, V M. Nakariakov, and T. D. Arber, Astron. and Astrophys. 395, 285 (2002).

21. E. N. Parker, Astrophys. J. 128, 664(1958).

22. E. N. Parker, Astrophys. J. 134, 20 (1961).

23. T. E. Holzer and W. I. Axford, Annual Rev. Astron. Astrophys. 8, 31 (1970).

24. M. V Eselevich and V G. Eselevich, Astronomy Reports 55, 359(2011).

25. Y. B. Zel’dovich and Y. P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena (Dover Publications Inc., New York, 2002).

26. I. Veselovsky, Astrophys. and Space Sci. 277, 219 (2001).

27. F. C. Bruhn, K. Pauly, and V Kaznov, in Artificial Intelligence, Robotics, and Automation in Space, Ed. by B. Battrick, ESA SP-603, id. 46.1 (2005).

28. Z. U. A. Warsi, Fluid Dynamics: Theoretical and Computational Approaches, 2nd ed. (CRC Press, Boca Raton Fla., 1998).

29. W. M. Macek and S. Redaelli, Phys. Rev. E 62, 6496 (2000).

30. S. Redaelli and W. M. Macek, Planetary Space Sci. 49, 1211 (2001).

31. Y. C. Whang, S. Liu, and L. Burlaga, J. Geophys. Research, 95, 18769(1990).

32. F. Malara, L. Primavera, and P. Veltri, Phys. Rev. E, 59 6023 (1999).

Перевод Е. Чмырёвой

Effects of Viscosity on the Behavior of Entropy Change in the ShockWave That Occurred after the December 13, 2006 Coronal Mass Ejection

H. Cavus and A. Kurt

Coronal mass ejections (CMEs) and the solar wind are the two main demonstrations of solar activity. These events can drive a shockwave. The shockwave occurs where the solar wind changes from being supersonic (with respect to the surrounding interplanetary medium) to subsonic. The main purpose of this study is to apply the algorithm and the results given in our recent papers to the shock wave that happened after the December 13, 2006 CME, and evaluate the behavior of entropy during this solar activity.

Keywords: solar wind—Sun: coronal mass ejections (CMEs)—Sun: activity

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 2 2015

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