CC BY
8
5. Bel'trami Je. Opyt interpretacii neevklidovoj geometrii [Experience in interpretation of non-Euclidean geometry] / Beltrami E. // Sbornik: «Ob osnovanijah geometrii» [The collection: "On the foundations of geometry"]. - M.: Izd-vo tehniko-teoreticheskoj literatury, 1956. - P. 180-212. [in Russian]
6. Staudt von G. K. Ch. Beitrage zur Geometrie der Lage, Erstes Heftr [Contributions to the Geometry of the situation, the First Heftr] / G. K. Ch. von Staudt - Nurnberg, Verlag von Bauer und Raspe (Julius Merz), 1856. - 129 p.[ in Germany]
7. Glagolev N.A. Proektivnaja geometrija [Projective geometry] / N.A. Glagolev - 2nd ed. - M.: Vysshaja shkola, 1963. -344 p. [in Russian]
8. Staudte von G. K. Ch. Geometrie der Lage [The geometry of the situation] / G. K. Ch. von Staudt. - Nurnberg: Verlag von Bauer und Raspe (Julius Merz), 1847. - 216 p. [in Germany]
9. Klejn F. O tak nazyvaemoj neevklidovoj geometrii [About the so-called non-Euclidean geometry] / F. Klein // Sbornik: «Ob osnovanijah geometrii» [The collection: "On the foundations of geometry"]. - M.: Izd-vo tehniko-teoreticheskoj literatury, 1956. - P. 253-303. [in Russian]
10. Fok V.A. Teorija prostranstva, vremeni i tjagotenija [Theory of space, time and gravity] / V.A. Fok - 2nd ed. - M.: Izd-vo fiziko-matematicheskoj literatury, 1961. - 563 p. [in Russian]
DOI: https://doi.org/10.23670/IRJ.2018.72.6.03
О КВАЗИРЕЛЯТИВИСТСКИХ УРАВНЕНИЯХ ДЛЯ ПЛАЗМЫ СОЛНЕЧНОЙ КОРОНЫ И ВАРИАНТ
ОБЪЯСНЕНИЯ ЕЕ ВЫСОКОЙ ТЕМПЕРАТУРЫ
Научная статья
Трощиев Ю. В.*
ORCID: 0000-0001-6805-7001, Московский государственный университет им. М.В. Ломоносова, Москва, Россия
* Корреспондирующий автор ([email protected])
Аннотация
В статье выводятся квазирелятивистские уравнения из нерялитивистских внесением в них соответствующих изменений. Они включают в себя уравнения непрерывности, законы сохранения импульсов и энергий как для протонов, так и для электронов. Это необходимо, так как гидродинамические и электрические силы, действующие на них, примерно одинаковы, а, следовательно, их импульсы тоже примерно одинаковы. Обычная скорость протонов в солнечной короне порядка 500 км/с, а соответствующая скорость электронов для такого импульса составляет 0.953 с, поэтому электроны необходимо учитывать и они релятивистские. Выражения для столкновительных членов тоже адаптированы к релятивистским скоростям. Тепло, выделяющееся в результате трения между электронными и протонными потоками, оказывается достаточным для нагрева плазмы. Концентрация ne =1010 см-3 и количество выделяющейся при рассеянии энергии компенсируют большое время свободного пробега.
Ключевые слова: плазма, солнечная корона, релятивистские эффекты.
ON THE QUASIRELATIVISTIC EQUATIONS FOR THE SOLAR CORONA PLASMA AND A VARIANT OF
RATIONALE OF ITS HIGH TEMPERATURE
Research article
Troshchiev Yu. V.*
ORCID: 0000-0001-6805-7001, Moscow State University, Moscow, Russia
* Corresponding author ([email protected])
Abstract
The quasi-relativistic equations are derived from nonrelativistic equations bringing amendments in them. They comprise of the continuity equations, and the conservation laws for impulses and energies both for protons and for electrons. It is so because the hydrodynamic and electric forces acting on them are similar, so, their impulses are also similar. The rife proton velocity in solar corona is 500 km/s and the consistent electron velocity for such impulse is 0.953 c, hence, the electrons must be taken into account and they are relativistic. Expressions for the collisional terms are also adapted to relativistic velocities. Heat emitting to warm protons due to friction between electron and proton flows is sufficient to heat plasma. Concentration ne = 1010 cm-3 and amount of scattering energy compensate large free path time.
Keywords: plasma, solar corona, relativistic effects.
Introduction
It is substantial to decide correctly what processes in the solar corona must be investigated and what processes can be neglected. Relativistic equations of solar corona plasma are derived in this article from the classic equations [1]. Influence of the relativistic electrons at plasma heating is investigated. Electrons are often neglected in plasma because of the supposition of plasma neutrality. Moreover, solar corona plasma is often treated as collisionless (collisions in plasma are not bumping but veering over interaction). Nonetheless, electrons suffer collisions and it plays a role in plasma heating. The problems of solar corona plasma heating are being intensively investigated now. But the problem of high temperature 106 K of the solar corona is still considered as unsolved. Some nonstandard source of extremely high temperature is proposed and investigated in this article. The result was reported by me at the conference [12].
Compulsory Consideration of Relativistic Effects while Modeling of the Solar Corona
If plasma is neutral then protons and electrons are settled equally. Thus, gradients of proton and electron pressures are similar and their impulses induced by pressure during a small time period dt are also similar.
Usual velocities of gas flows are 500 km/s in the quiet solar corona. If protons move at Vp=500 km/s then the classic equation for impulses mpVp=meVe gives superlight velocity Ve for electrons. mp, me - masses of proton and electron correspondingly. The relativistic equation for impulses
m-pVp _ rneVe ^
gives electrons velocity
Ve = c = 0.9 5 3 c (2)
and increasing of their orthogonal mass in 3.3 times. Here c is light velocity. If the proton wind is 10 times stronger then electrons velocity is 0.9995 c and orthogonal mass increases in 31.5 times.
This is for gas dynamic forces. Let now plasma is diverted a little from neutrality. Then electric forces will bring the same impulses to protons and electrons during a small time dt because charges of proton and electron are the same. Magnetic forces change only the direction of impulses. So, we can suppose that the normal situation is moving of protons and electrons with approximately same impulses. This abstract reasoning leads to infer: 1) electrons must be taken into account in spite of their little mass because they carry similar to protons impulses, 2) if the devices register proton velocities 500 km/s in the solar corona then the electrons move there at relativistic velocities. The mechanisms of electron acceleration can be diverse. It seems to be an antinomy that protons and electrons move with different velocities in neutral plasma but this antinomy is resolved in the last section of this article that concerns the results of investigation.
The derived hereinafter relativistic equations are more precious than classical ones but they are derived from them, so, it is more natural to call the derived equations quasirelativistic.
Quasirelativistic Equations of Completely Ionized Plasma
It is a way of constructing of the equations in what follows. Some terms in section 4 are given for two-dimensional case without loss of generality. It is possible to write three-dimensional terms and include additional physical processes. The quasirelativistic system of equations includes the terms similar to the classical system in [1] with special corrections for relativistic case and for solving of the problem of high temperature of the solar corona. The corrections are formulation of impulse and energy conservation laws directly for impulses and energies and not for velocities and temperatures, adding one more term and rearranging of one term from the heat flows to the heat sources. The derivation of the friction term and of the corresponding heat source is made in an unusual manner to allow hydrodynamic velocities to be much greater than thermal velocities.
Quasirelativistic equations of completely ionized plasma are the following:
dyie
-^ + div(neVe) = 0, drti
= 0,
deiVe) 9pe
ne d = - — -ene( E + [VeB]) + F/r,
cZ,(p.) dp,
--^ + zteni(E + M) - F/r-
de(Ee)
ne + pediv(Ve) + div(qe) = Qe,
n + pjdiv (V;) + div ( q 0 = Q (3)
where r = (x,y,z) - right-handed coordinates, t- time, indices 'e', '/'denote electrons and ions correspondingly; na -concentrations, - average impulses of the particle, - average heat energies of the pResearch articles, Va - velocities, ma - masses, pa - pressures, e - electron charge, Z{ - ion charge in electron charge units, E - intensity of electric field, B -induction of magnetic field, Ff. - includes friction force between proton and electron flows and thermal force (Ff. is interaction between electrons and protons, so, its directions are opposite), qa - heat flows, Q - heat sources. Viscosity is neglected. Particle derivative is
^ = £+(VaV) , (4)
Symbolic index a can assume values 'e' or '/'. SI system of units is applied in (3).
The equations (3) are three-dimensional. It is more convenient to investigate two-dimensional equations in many cases. Reflecting reality they spend much less computer time and memory for calculations.
Two-dimensional equations look alike (3) with another sense of variables. Namely, Bx = By = 0, Bz = B and no variable depends on z coordinate. Components Ez = 0, Bx and By will remain zeros under these assumptions.
The more easy two-dimensional case is additionally Vez = Viz= 0. z-components of friction and heat flows are linear combinations of the vectors u=Ve-V;, VTe and VTi, so, Fz = 0 and qez = qiz = 0 under zero z-velocities. It is important that they will remain zeros under relativistic summing of the velocities in this case.
Physical Values Transformations for Different Coordinate Systems
It is more convenient to calculate some physical values in the fixed coordinate system and some ones in the systems moving with the flows of protons or electrons. Thus, two relativistic orthogonal masses are introduced: mass in the fixed system for electron-proton collisions
. C2
and mass in the moving system for heat exchange processes inside the flow
m a = . , (6)
J
where R - universal gas constant, f - molar mass. One wave over a symbol denotes the fixed system. It can be omitted if there is no ambiguity. Two waves denote moving system.
Free path times of the pResearch articles (5) and (6) are different:
and
mjc4 1- 1-V//C
T = --—--(7)
e 4тге4п;|7е-7;|Л ' v '
A = lnf^ I Wl ) (8)
~ _ 9(fcBra)2(l + Vl-Qq)2 (9)
a' 1 6 ttZ2Z?e4 n iA c,/^' ( )
Jbl + 4 - ba
3 ki>Tn
mr,c
« (23.4 - 1.15 lgn + 3.45 lg fe at fe < 50 eV
m = • & _ _e (10)
I 25.3 — 1.15 lgn + 2.3 lgfe at fe > 50 eV
Gauss units system is conventional in astrophysics, so, it is used here and in what follows. Additionally, temperature in electron-volts is remained in practical formulas from [1] (e.g., (10)). A - Coulomb logarithm, kB - Boltzmann constant, Ze=1.
The fact that even at relativistic velocities electrons are at least 10 times lighter and move faster than protons is used in (7) and (9).
Let the set of three vectors e1 = V^VJ, e2, e3 in the fixed coordinate system is a right-hand coordinates. They are column-vectors in coordinate notation. Then the matrix A = (e1,e2,e3)r is a matrix of transformation of the fixed coordinate system to coordinates defined by these orts. The matrix AT = A-1 is a reverse transformation. The relativistic transformation of the velocity to the velocity is defined now
m = Bm"-Va (11)
a 1 - (V a,m a) /c2 , (11)
B=A :DA, D = di ag ( 1 , / 1 - Va2/c2 ,/ 1 - Va2/c2 ). diag denotes a diagonal matrix. The matrix B does not depend on the choice of e2 and e3. It can be represented independently from them: B = P + / 1 - Va2/c2 ( E - P) , where E - unit matrix, P -projector to the linear span of ei.
The inverse transformation of the velocity is correspondingly:
Va= B m a: Vf (12)
a 1 + (V a,m a) /c2 ( )
The rule of transformation of the force from the moving system to the fixed one follows from these formulas:
Va , dPa p _ \l_Yl . d£a . dVa . p
dpa ^ c2 dVa dVa dpa
Let now the average energy of the particles in the moving system is < Ea) . Then their root mean square (rms) velocity is defined by the expression
F = I l-^--^F = I l-ia-.^.^.^F. (13)
->1 r2 rfK. л г2 ли JV JS V >
< ^ 2=C2( 1-( ^S) j ■ (14)
Let the square of the rms velocity in the fixed system be an averaging over the sphere of radius of the square of the velocity :
<ia> 2 = ^ J (ma (ma) ) 2 dn= C2 - (c2 - 7a2 ) ( C2 - < £ > 2) (15)
(in (l + Va (Pa)/c2) - In (l - Va (i)/c2fj it) '
where dQ - element of frustum. Then average heat energy is
<£a>=maC2( h i.2/ 2- ) 1 I ■ (16)
Wl-(Va)2/d !-V2/c2
Time and volume also change while transition from one coordinate system to another but they compensate each other because of invariance of spatiotemporal intervals on Lorentz transform. So, the formulas (14)-(16) are transformations of heat energy for transitions from one coordinate system to another.
Physical Terms Included in the Equations
It is possible now to write expressions for friction, heat flows and sources. Friction consists of two components - friction force and thermal force:
F/r - F/ 10w + Fr. (17)
Friction force is
F/;ow--^p. (18)
Thermal force is
Fr = -- II-
3 -, v& dPa dVa dVa ne
[h VfJ , (19)
2 e2 dVa dVa dpa o>ete I
where - cyclotron frequency in the coordinate system moving with the flow, h - unit vector of magnetic
induction direction.
The source being generated by electron heat flow consists of two components:
dTv (q S)+;^cGV (q %) . (20)
Usage of expressions from [1] and bringing relativistic effects to them provides:
3 neT( 2 5>eTe
4u - ^ M, qer - -*-eVfe -[hVq] , (21)
ке = 4 . 6 6 Jill , (22)
mea.i|Tei
where 3e = e B/m ec - cyclotron frequency in the fixed coordinate system.
The source being generated by ion heat flow also consists of two components:
^ (q + (q T) . (23)
Usage of expressions from [1] bringing relativistic effects to them and additionally inclusion of the term q- gives:
q- = -Ш[Ы] qT = (24)
_ ^ = . (25)
where m, - cyclotron frequencies: 3; = eB/m;c, = eB/m;c. Another way to find q- is given hereinafter.
The flows (21) and (24) are written in supposition of great values of cyclotron frequencies. Crux of the matter is that they define heat emitting due to electron and ion collisions. It follows from (18) that force of friction is
F/i0w = ft(V - ve) , ft = ^ > 0 . (26)
The kinetic energy is E = c (/ m2c2 + p2 — m c) . Heat dW emitting in the unit value during dt due to friction is apparently
ned <Ee> + nd <E> + dW = 0 , (27)
so,
d w = — (n e^fdp e + n i^dp ;) * ft(Vi — Ve) 2 d t > 0. (28)
Electrons scatter more than protons and Ve>>Vi (it was assumed in expression for Ffow too). If all electron kinetic energy transforms to heat in collisions then
dWe = ne<£e> P. (d t) = ^d t, (29)
dWj « 0.5 n^p^dt) = 7^-dt,
2re
where p„(dt) - probability of electron-proton collision during dt. Then,
dwe _ {Ee> _ (ve)
dWi (Et) 2 (Vi)
(30)
since
where ra - 1 - V^2/c2. So,
Ea=^|,Pa=!;7^,ra<l,<Pe>«<Pi> (31)
ra+ra ra
dWi*d wg*^<7e><7i> (32)
(iy 2re
These formulas define the minimal estimation for heating of proton gas. To find the upper estimation it is natural to assume dWe = dW, = 0.5dW.
Results of Investigation
Temperature 10.27106 K corresponds to proton velocity 500 km/s and temperature 8.99109 K corresponds to electron velocity 0.953 c. So, the high temperature source for plasma heating exists in spite of solar surface temperature 6000 K. This mechanism is transformation of hydrodynamic velocities into heat velocities in electron-proton collisions. Estimation of energy coming and losing is another physical problem which is out of the framework of this article. But it worth to note that Sun produces much more energy than it is necessary to heat the corona and that solar corona radiation is very faint.
The lower estimations of temperature increasing depending on proton velocity at concentration 1010 are in the Table 1 according to the formula (32).
Table 1 - Rate of temperature increasing (K/s) depending on velocity of protons (km/s)
10 100 500 1000
dT/dt 20000 16000 110 39
Preliminary computer calculations of the system (3) by difference schemes corroborated increasing of temperature. The initial pressure gradients cause then plasma oscillations that seem to be the main mechanism of electron-proton collisions. Langmuir frequency is approximately 106 s-1 here. These oscillations also answer the question how it is possible for protons and electrons to move with different velocities in neutral plasma. And the calculations confirm similarity of their pulses.
Конфликт интересов Conflict of Interest
Не указан. None declared.
Список литературы / References
1. Braginskii S. I. Transport processes in a Plasma / S. I. Braginskii // Reviews of Plasma Physics, ed. M. A. Leontovich, Consultants Bureau - 1965 - V. 1 - P. 205.
2. Загадки солнечной короны // Наука в России - 2011 - № 4 - С. 61-63.
3. Molotkov I.A. Solar corona top heating / I.A. Molotkov, N.A. Ryabova // Geomagn. Aeron - 2016 - V. 56 - № 3 - P. 264-267. doi: 10.1134/S0016793216030130
4. Badalyan O.G. To the problem of solar coronal heating / O.G. Badalyan, V.N. Obridko // Astron. Lett. - 2007 - V. 33 -N 3 - P. 182-191. https://doi.org/10.1134/S106377370703005X
5. Somov B.V. Plasma astrophysics. Part II. Reconnection and flares / B.V. Somov // Springer - 2006 - 417 P.
6. Kovalev V.A. The Role of Collisions in the Particle Acceleration in Solar-Flare Magnetic Traps / V.A. Kovalev, B.V. Somov // Astron. Lett. - 2003 - V. 29 - N 6 - P. 409-415. doi: 10.1134/1.1579790
7. Bilenko I.A. On the regimes of heating during solar flares / I.A. Bilenko, V.A. Kovalev // Astron. Lett. - 2009 - V. 35 - N 11 - P. 791-797. doi: 10.1134/S1063773709110061
8. Kovalev V.A. The acceleration of evolution caused by a decrease of dissipation / V.A. Kovalev // The complex systems
- 2017 - N 1(4) - P. 13-20.
9. Bret A. Collisional behaviors of astrophysical collisionless plasmas / A. Bret // J. Plasma Phys. - 2015 - V. 81. doi:10.1017/S0022377815000173
10. Kurkina E. S. Formation of Thermal Structures with Blowup During Solar Flares / E. S. Kurkina, E. D. Kuretova, V. A. Kovalev // Comput. Math. Model. -2015 - V. 26 - P. 144-155. doi:10.1007/s10598-015-9262-5
11. Kurkina E. S. Small-Scale Heat Localization with Blowup in the Magnetic-Tube Cross Section During a Solar Flare / E. S. Kurkina, Y. V. Troshchiev, V. A. Kovalev, et al // Comput. Math. Model. - 2016 - V. 27 - P. 395-416. doi:10.1007/s10598016-9330-5
12. Troshchiev Y.V. Quasirelativistic equations for the solar corona plasma and a variant of rationale of its high temperature / Y.V. Troshchiev // Theses of the conference High Enegry Phenomena in Relativistic Outflows VI, Moscow, Russia, IKI RAS, 11-15 September - 2017 - P. 22.
Список литературы на английском языке / References in English
1. Braginskii S. I. Transport processes in a Plasma / S. I. Braginskii // Reviews of Plasma Physics, ed. M. A. Leontovich, Consultants Bureau - 1965 - V. 1 - P. 205.
2. Zagadki solnechnoi korony [Mysteries of the solar corona] // Nauka v Rossii [Science in Russia] - 2011 - № 4 - P. 6163. [in Russian]
3. Molotkov I.A. Solar corona top heating / I.A. Molotkov, N.A. Ryabova // Geomagn. Aeron - 2016 - V. 56 - № 3 - P. 264-267. doi: 10.1134/S0016793216030130
4. Badalyan O.G. To the problem of solar coronal heating / O.G. Badalyan, V.N. Obridko // Astron. Lett. - 2007 - V. 33 -N 3 - P. 182-191. https://doi.org/10.1134/S106377370703005X
5. Somov B.V. Plasma astrophysics. Part II. Reconnection and flares / B.V. Somov // Springer - 2006 - 417 P.
6. Kovalev V.A. The Role of Collisions in the Particle Acceleration in Solar-Flare Magnetic Traps / V.A. Kovalev, B.V. Somov // Astron. Lett. - 2003 - V. 29 - N 6 - P. 409-415. doi: 10.1134/1.1579790
7. Bilenko I.A. On the regimes of heating during solar flares / I.A. Bilenko, V.A. Kovalev // Astron. Lett. - 2009 - V. 35 - N 11 - P. 791-797. doi: 10.1134/S1063773709110061
8. Kovalev V.A. The acceleration of evolution caused by a decrease of dissipation / V.A. Kovalev // The complex systems
- 2017 - N 1(4) - P. 13-20.
9. Bret A. Collisional behaviors of astrophysical collisionless plasmas / A. Bret // J. Plasma Phys. - 2015 - V. 81. doi:10.1017/S0022377815000173
10. Kurkina E. S. Formation of Thermal Structures with Blowup During Solar Flares / E. S. Kurkina, E. D. Kuretova, V. A. Kovalev // Comput. Math. Model. -2015 - V. 26 - P. 144-155. doi:10.1007/s10598-015-9262-5
11. Kurkina E. S. Small-Scale Heat Localization with Blowup in the Magnetic-Tube Cross Section During a Solar Flare / E. S. Kurkina, Y. V. Troshchiev, V. A. Kovalev and others // Comput. Math. Model. - 2016 - V. 27 - P. 395-416. doi:10.1007/s10598016-9330-5
12. Troshchiev Y.V. Quasirelativistic equations for the solar corona plasma and a variant of rationale of its high temperature / Y.V. Troshchiev // Theses of the conference High Enegry Phenomena in Relativistic Outflows VI, Moscow, Russia, IKI RAS, 11-15 September - 2017 - P. 22.
