Быстрые и медленные МГД-волны в термически активном плазменном слое

Обложка


Цитировать

Полный текст

Аннотация

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

Полный текст

Introduction

The solar atmosphere is the rarified and highly-magnetized plasma. Due to this fact, there are numerous magnetic structures, which exist in the solar corona including coronal loops, prominences, polar plumes and etc. In fact, these magnetic structures play the role of the wave-guides for the compressional and compressional MHD waves, which routinely observed in the solar corona [1; 2]. The geometry of these wave-guides is prescribed by the balance between the gas-dynamic and magnetic forces and lead to the dispersion of the compressional modes.

However, the magnetic structures owe their existence not only to mechanical but also to the thermal balance. The coronal cooling/heating rates are temperature and density dependent. As a consequence, some compressional perturbation can disturb balance between these processes, and can be amplified or damped by the additional impact from the non-adiabatic processes. Thus, there are some feedback between waves and plasma. In other word, the coronal plasma is the thermally active one [3; 4]. Furthermore, the non-adiabatic processes affect propagation speed of the compressional perturbation [5].

Further, we analyze the combined influence of the thermal activity and the magnetic structuring on properties of the compressional magnetohydrodynamic (MHD) waves propagating in the highly magnetized plasma. Most similar studies use strong limitations on the wavelength and the magnitude of the external magnetic field [6; 7]. In our study, we will not use such strong assumption and apply the most general approach.

The paper is organized as follows. In the Section 2, we introduce the initial equations and obtained dispersion relation. In the Section 3, we show calculated dispersion curves. In the Section 4, we discuss the obtained results.

1. Model and dispersion relation

We analyze the waves in the fully ionized optically thin coronal plasma. We neglect the gravitational forces and the effects of viscosity, thermal conductivity and finite magnetic conductivity. The difference of basic equations from equations describing the ideal plasma is in the energy transport equation:

ργγ1DDtPργ=ρW(ρ,T)=Lρ,TΓρ,T, (1.1)

Here, ρ, T, and P are the plasma density, temperature and pressure, respectively. The parameter γ, is for the adiabatic index. In the equation (1) the function Wρ,T describes the net heat-loss function which is the difference between radiative cooling Lρ,T and some heating Γρ,T processes. These rates balance each other in the case of the stationary conditions Lρ0,T0=Γρ0,T0, or Wρ0,T0=0.

To analyze the properties of the MHD waves we apply the standard methods of the perturbation theory. In other words, we linearize basic equations using substitution of the following form: a=a0+a1,a1/a0~ε<<1, where a is some plasma parameter. We model the coronal wave-guide by the magnetic slab with the magnetic field directed along z-axis as follows:

B0x=Bi,xx0,Be,xx0. (1.2)

Here, B0 and x0 are the stationary value of magnetic field strength and the slab half-width, respectively. Searching the solution of the linearized equation in the form a1=a~1xeiωt+kzz, one may obtain dispersion relation for the set (fast/slow and body/surface) of the sausage/kink magnetoacoustic waves in the magnetic slab as follows:

(kz2cA2iω2)kxekxi=ρ0eρ0ikz2cA2eω2​​coth(kxix0)tanh(kxix0). (1.3)

The complete derivations steps can be found in our previous study (see [8]). We use the and hyperbolic functions for the kink and sausage modes, respectively. We introduce the following characteristic temporal scale to describe the thermal activity influence:

τV=CV/W0T, τP=CPT0/W0TT0W0ρρ0, (1.4)

where CV and CP are the specific heat capacities at constant volume and constant pressure, respectively. We also use the derivatives W0T=Q/T|ρ0,T0,W0ρ=Q/ρ|ρ0,T0

Some characteristic speeds of fast and slow MHD waves are shown below:

cA=B024πρ0, cS=γkBT0m, cSQ=τVτPγkBT0m, cT=cS2cA2cA2+cS2, cTQ=cA2cSQ2cA2+cSQ2. (1.5)

Here, speed cA is for the Alfven speed. The speed cS is the usual speed of sound for ideal plasma. The so-called tube speed cT is a result of the pure wave-guide dispersion effect. The speed cSQ is long-wavelength limit value of the slow-wave or acoustic perturbation in the case the thermally active uniform plasma. And the last but not the least is the speed cTQ, which is a result of the both wave-guide and thermal activity dispersion effects.

In the dispersion relation (3) we also use following notations:

kxi,e2=AQi,e2mQi,e2+iωτVi,eAi,e2mi,e2AQi,e2+iωτVi,eAi,e2,

m2=kz2cA2ω2kz2cS2ω2cA2+cS2kz2cT2ω2, mQ2=kz2cA2ω2kz2cSQ2ω2cA2+cSQ2kz2cTQ2ω2,

A2=cA2+cS2kz2cT2ω2, AQ2=cA2+cSQ2kz2cTQ2ω2.

Further, we will use dispersion relation (3) to calculate the dependencies of phase velocities and increment/decrement of MA waves on the wavenumbers in the solar atmosphere conditions.

2. Results

Let us apply the obtained dispersion relation to the coronal conditions. In the current study we will focus on properties of waves in the highly magnetized coronal loop. The magnetic slab parameters corresponding to the mentioned loop are shown in Table.

 

Table

Slab parameters used for calculations

Таблица

Параметры слоя, используемые для расчетов

 Parameter

 Value

Magnetic field strength inside the slab B0i

 100 G

Temperature T0i

 6 MK

Number density inside the slab n0i

 1011 cm

Density contrast n0i/n0e

 10

Slab width 2x0

 2 Mm

 

The parameterization of the radiative cooling function L(ρ,T)=χρTα, has been calculated using the CHIANTI Version 10.0 database [9]. The heating rate has been modeled as Γ(ρ,T)=hρ0.5T3.5. Such heating scenario has been seismologically proposed by [10] using the observations of the damped slow magnetoacoustic waves in the long-lived coronal plasma structures. The dispersion curves for body fast/slow sausage/kink magnetoacoustic waves are shown in Figures 1, 2.

 

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

 

Рис. 2.2. Декремент волн перетяжек и изгибных волн в сильно замагниченной корональной плазме (см. таблицу). Верхняя панель соответствует волнам перетяжек. Нижняя панель предназначена для изгибных волн. Разные цвета соответствуют различным модам

 

It can be easily seen that the slow modes in the thermally active plasma can be found between sound speed cSi and the modified tube speed cTQi. In the ideal plasma case, the long-wavelength limit is cTi. The fast modes in the plasma with the thermal misbalance range between cAe and cAi. The slow waves are highly affected by both thermal activity and wave-guide dispersion. The impact of the thermal activity on the fast-wave dispersion is negligible.

One may notice that decrement of both fast sausage and kink modes are lower than decrement of slow waves. This effect is in agreement with result obtained for the uniform plasma [4]. However, the magnetic structuring leads to non-monotonic behavior of wave decrement, with maximum in the long-wavelength part of the spectrum. In uniform plasma behavior was shown to be monotonic [4]. The calculated decrement of rather weak and cannot explain observed fast wave damping. The slow-wave decrement is comparable with the observed decay time. Moreover, in highly magnetic plasma decrement of slow-waves become greater (compare with results obtained for the hot loop in [8]. This is due to the fact that with decrease of plasma beta/increase of magnetic field the slow wave becomes more acoustic than magnetic mode.

3. Discussion and conclusions

In the presented study we analyzed the combined influence of the thermal activity and the magnetic structuring on properties of the compressional magnetohydrodynamic (MHD) waves. Using perturbation theory, assumption of strong magnetic structuring and the slab geometry, we obtain the dispersion relation for the set (fast/slow and body/surface) of the sausage/kink magnetoacoustic waves. We solve the obtained dispersion relation numerically and use to the higly-magnetized solar corona conditions. We showed that slow-waves are higly affected by both thermal activity and wave-guide dispersion. In particular the long-wavelength becomes equal to cTQ, which is defined not only by geometry of the wave-guide but also by the acting non-adiabatic processes. As a result, the usage of the value cT for the helioseismological needs may cause mistakes. On the contrary, the oscillation of the fundamental modes may be used for phenomenological determination of unknown coronal heating function. We also showed the phase velocity dispersion for the fast MHD waves remains the wave-guide geometry. In the magnetically structured plasma the wave decrement becomes non-monotonic with maximum in the long-wavelength part of the spectrum. The calculated slow wave decrements are comparable with the observed decay times.

The study was supported in part by the Ministry of Education and Science of Russia by State assignment to educational and research institutions under Project No. FSSS-2020-0014 and No. 0023-2019-0003, and by RFBR, project number 20-32-90018. CHIANTI is a collaborative project involving George Mason University, the University of Michigan (USA), University of Cambridge (UK), and NASA Goddard Space Flight Center (USA).

×

Об авторах

Дарья Вадимовна Агапова

Самарский национальный исследовательский университет имени академика С.П. Королева; СФ ФИАН

Автор, ответственный за переписку.
Email: agapovadaria2019@gmail.com
ORCID iD: 0000-0002-3957-7339

магистрант физического факультета; инженер, теоретический сектор

Россия, г. Самара; г. Самара

Сергей Александрович Белов

Самарский национальный исследовательский университет имени академика С.П. Королева; СФ ФИАН

Email: mr_beloff@mail.ru
ORCID iD: 0000-0002-3505-9542

аспирант кафедры физики; научный сотрудник, теоретический сектор

Россия, г. Самара; г. Самара

Нонна Евгеньевна Молевич

Самарский национальный исследовательский университет имени академика С.П. Королева; СФ ФИАН

Email: nonna.molevich@mail.ru
ORCID iD: 0000-0001-5950-5394

доктор физико-математических наук, профессор кафедры физики; главный научный сотрудник, и.о. зав. теоретическим сектором

Россия, г. Самара; г. Самара

Дмитрий Игоревич Завершинский

Самарский национальный исследовательский университет имени академика С.П. Королева; СФ ФИАН

Email: d.zavershinskii@gmail.com
ORCID iD: 0000-0002-3746-7064

кандидат физико-математических наук, доцент кафедры физики; научный сотрудник, теоретический сектор

Россия, г. Самара; г. Самара

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

  1. Nakariakov V.M. [et al.] Kink Oscillations of Coronal Loops. Space Science Reviews, 2021, vol. 217, issue 6, article number: 73. DOI: http://doi.org/10.1007/s11214-021-00847-2.
  2. Wang T.J. Waves in Solar Coronal Loops. In: Low-Frequency Waves in Space Plasmas, 2016, pp. 395–-418. DOI: http://dx.doi.org/10.1002/9781119055006.ch23.
  3. Zavershinskii D., Kolotkov D., Riashchikov D., Molevich N. Mixed Properties of Slow Magnetoacoustic and Entropy Waves in a Plasma with Heating/Cooling Misbalance. Solar Physics, 2021, vol. 296, issue 6, no. 96. DOI: https://doi.org/10.1007/s11207-021-01841-1.
  4. Zavershinskii D.I., Molevich N.E., Riashchikov D.S., Belov S.A. Nonlinear magnetoacoustic waves in plasma with isentropic thermal instability. Physical Review E, 2020, vol. 101, issue 4, p. 43204. DOI: https://doi.org/10.1103/PhysRevE.101.043204.
  5. Belov S.A., Molevich N.E., Zavershinskii D.I. Dispersion of Slow Magnetoacoustic Waves in the Active Region Fan Loops Introduced by Thermal Misbalance. Solar Physics, 2021, vol. 296, issue 8, no. 122. DOI: http://doi.org/10.1007/s11207-021-01868-4.
  6. Zhugzhda Y.D. Force-free thin flux tubes: Basic equations and stability. Physics of Plasmas, 1996, vol. 3, issue 1, p. 10–21. DOI: http://dx.doi.org/10.1063/1.871836.
  7. Zavershinskii D.I., Kolotkov D.Y., Nakariakov V.M., Molevich N.E., Ryashchikov D.S. Formation of quasi-periodic slow magnetoacoustic wave trains by the heating/cooling misbalance. Physics of Plasmas, 2019, vol. 26, issue 8, p. 82113. DOI: https://doi.org/10.1063/1.5115224.
  8. Agapova D.V., Belov, S.A., Molevich, N.E., Zavershinskii D.I. Dynamics of fast and slow magnetoacoustic waves in plasma slabs with thermal misbalance. Monthly Notices of the Royal Astronomical Society, 2022, vol. 514, issue 4, pp. 5941–5951. DOI: http://doi.org/10.1093/mnras/stac1612.
  9. Del Zanna G., Dere K.P., Young P.R., Landi E. CHIANTI - An Atomic Database for Emission Lines. XVI. Version 10, Further Extensions. The Astrophysical Journal, 2021, vol. 909, no. 1, p. 38. DOI: http://doi.org/10.3847/1538-4357/abd8ce
  10. Kolotkov D.Y., Duckenfield T.J., Nakariakov V.M. Seismological constraints on the solar coronal heating function. Astronomy & Astrophysics, 2020, vol. 644, issue 1, p. A33. DOI: http://doi.org/10.1051/0004-6361/202039095.

Дополнительные файлы

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

Скачать (593KB)
3. Рис. 2.2. Декремент волн перетяжек и изгибных волн в сильно замагниченной корональной плазме (см. таблицу). Верхняя панель соответствует волнам перетяжек. Нижняя панель предназначена для изгибных волн. Разные цвета соответствуют различным модам

Скачать (272KB)

© Агапова Д.В., Белов С.А., Молевич Н.Е., Завершинский Д.И., 2022

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах