Showing posts with label Vlasov-Boltzmann equation. Show all posts
Showing posts with label Vlasov-Boltzmann equation. Show all posts

Friday, August 28, 2020

Why The MHD Regime Is Preferred In Doing Solar Physics



Most solar physicists work in what's called the "MHD" or magneto-hydrodynamic regime, which numerous plasma physicists have criticized as "taking a short cut". For example, treatments based on the analysis of the indigenous magnetic fields currently occupy a predominant role in analysis of energy stored in sunspot regions. The majority of these are based on MHD  theory wherein low frequency, long scale length and quasi-neutral (electric charge) assumptions are introduced into the plasma physics. Such assumptions lead to a basic set of what we call "MHD" equations, which include:

a)( r  / t)  +  Ñ (r v) = 0

b) r (  V / t)  +    r v Ñ ( v) =   -  Ñ p  + 1/c  J X B

c)  E + v X B = 0

More fundamentally, the MHD regime is really an approximation of the one-fluid and two-fluid cases, with many added simplifying assumptions – and these "take a toll"  in ignoring the kinetics and micro-physics.  In truth, MHD does represent the end product of successive information loss in the transition from Vlasov theory to two-fluid theory to one-fluid theory and thence MHD.   But the benefit of using it is self-evident in many flare models, for example, so that one is no longer hampered by being unable to "see the forest for the trees". 

One starts with the Vlasov-Boltzmann equation:

f/ t + v Ñ z f + (F/ m ) Ñ v f = ( f/ t)c

Whereby the lowest moment, for example, is obtained by integrating the preceding equation with F as the Lorentz force:

ò ( f/ t)  dv    + ò  v Ñ z f dv  + (q/ m ) ò (E + v X B) f/ v  =  ò ( f/ t)c   dv


The 2nd moment is obtained by multiplying the original eqn. (Boltzmann) by mv then integrating it over dv.

The progression from using this procedure is that one gets in succession:

Two -fluid theory (e.g. ions and electrons treated as separate fluids), i.e. in fluid eqns.:


1) ( r a  / t)  +  Ñ (r a v a) = 0

where the index  a = e, i.


for electrons, ions.   

2)  r a (  V o /  t)  +    r v a · Ñ ( v a) =   -  Ñ p a +  e a n a  (E + v a x B)

3)  p o  r a  g    =  const.

Plus the relevant form for Maxwell's equations.

More technically as Dwight R.Nicholson notes (Plasma Theory, p. 129) one takes "velocity moments of the Vlasov equation in seven dimensional (x, v, t) space."  From these (ibid.): "an infinite hierarchy of equations in 4-dimensional (x, t) space can be derived."   And: "when an appropriate truncation of this infinite hierarchy is carried out, the standard two fluid theory of plasma physics is obtained."    Thus we can proceed to treat:

 f (x, v, t)  +   v · Ñ x  f o  +   (q/ m) (E + v/ c  x B) Ñ v · f o   = 0

With the normalization: n o  (x, t) =  ò  dv f o  (x, v, t)   from which the 2 fluid velocity equations are obtained.  But given these we go direct to....
!
!
!
!
v
One fluid theory (introducing low frequency, long wave length and quasi -neutral approximations, e.g.   n e=    n i


  ( r  / t)  +  Ñ (r v) = 0

where no separate  index for  e, i, is needed.  Then, more quickly to...
!

!
v
MHD Theory

(This proceeds from 1-fluid theory with a few further assumptions, simplifications)

Apart from this, the most conspicuous agent of solar plasma containment and topological change is the coronal loop –


 Which is fully consistent with a low beta-plasma regime- since otherwise gas kinetic pressure would predominate and the configuration would be ‘open’.  Note: 

beta = (½  r  v 2)/   B 2/2 m 0 

where  r is the fluid density, v the fluid velocity, B the magnetic induction and m 0  the magnetic permeability of free space.) 

With this little exercise, we see that the solar physicist is more than amply justified in employing the MHD approach.   This is directly connected to the gas kinetic pressure (½  r  v 2)  being much less than the magnetic pressure  ( B 2/2 m 0 ), so the density of plasma particles  (n e. ,  n i ) is much less  - reducing the need to overemphasize the micro physics especially as the Vlasov collision term  ( f/  t)c  -> 0. Indeed, for most applications, i.e.  for assessing magnetic fields in  sunspot regions, two fluid theory would be far too complicated while achieving little additional insight.


Monday, March 21, 2016

Coronal Mass Ejections and Collisionless Shocks (1)


The CME captured by the Solar Dynamics and Heliospheric Observatory in February, 2013.

Coronal mass ejections are powerful, energetic events on the Sun which can also have serious repercussions for terrestrial life, see e.g.

http://brane-space.blogspot.com/2012/06/cme-that-you-dont-want-to-see.html

How do coronal mass ejections occur? This is believed to happen via a sudden, mass transient ejection probably via a high Mach number shock wave. It has also been found that collisionless shocks form ahead of CMEs when the latter’s velocity exceeds the Alfven velocity of the ambient plasma in the corona. Recall this Alfven velocity is given by:



VA     =  Ö (Bo / mo   r o)


To fix ideas, let us recall the corona is a collisionless plasma . (Recall in the limit 1/N << 1 and 1/N ® 0  we say a plasma is “collisionless” where:N  =   4p n o  l3 D  / 3   and n o  is the number density with   l D  the Debye length.   (See example problems from Jan, 23rd). Let a coronal plasma have number density n = 10 16 /m3)  and a kinetic temperature T =10 6 K.  The plasma parameter is defined:


L  =  n o  l3 D

Where the Debye length

lD,s =[kT s εo / 4p  Z s2 n s  e2 ] ½

Here:  Z s   = 1 so that:

l D,s   =

[(1.38 x 10-23  (10 6 K) (8.85 x 10 -12 F/m)  / (4p)  (10 16/m3)  e2 ] ½

l D,s   =  1.9 x 10 - 4   m

L  =  7.3 x 10 4   

And N  =   4p n o  l3 D  / 3 

N =  (4p)(7.3 x 10 4   ) / 3  »   3.1 x 10 5   

Therefore:   1/N   =   1/ (3.1 x 10 5 )
 

1/N = 3.2 x 10 - 6   


 1/N ® 0  , hence the  plasma is “collisionless” so that this meets the condition.

Then we may use the collisionless Vlasov-Boltzmann equation for a plasma description, e.g.:


v f x/ x  - q(s)/m(s) [ j/ x  · f s/ v]  = 0

where (f s) is any function of the constant of the motion.

This, along with its collisional form, has been called “the most important equation in plasma physics”.[1] The associated Poisson equation is:

Ñ E = 4p e (ò fe d3 V - òf i d3V) = 0


where the e, i –subscripts denote electrons and ions, respectively and f refers to the appropriate distribution for each.

In the classical theory of shock waves in collisionless plasmas, one would not expect dispersion to be a sufficient condition to incite a shock.  (A dispersion relation implies that a relationship exists between the plasma frequency w and the wave number k. )  What it can do is form a large amplitude fixed wave called a “soliton”. (‘Steepening’ + dispersion yields solitons, while steepening + dissipation yields shocks). These soliton solutions ensure no turbulence (due to shocks), thereby meeting the collisionless condition. To estimate M we use (e.g. Sakai and Ohsawa, 1987)[2]:

B m =  (2M – 1) Bo

Or in terms of the ratio of the magnetic fields:


M =  2 [B m  / Bo] + 1

To reach an M= 80 value, say which might conceivably apply for coronal shocks, one would need:
[B m  / Bo]  »  39.5


I.e. the ratio of the maximum magnetic field intensity to the equilibrium value is at least 39.5.

 Compression is also critical to have an effective shock and we have the condition, based on the plasma number density N:

DN/ N   =  DB/ Bo  [1 - DB/ 2 B] - 1

We already saw for the equilibrium collisionless corona:

N »   3.1 x 10 5   

It is instructional here to determine the trigger or break point for the transient, or CME. This should be reached when the massive influx of dissipating particles is so great that:

DN/ N    ® ¥

 
By inspection we see that if DB = 2B :

¥   =   2B/ Bo  [1 – 2B/ 2 B] - 1

 
To achieve this one would most likely need a strong two-stream instability requiring (cgs units):
M > 1 +  [8p N Te / B2 ]  1/3

This elicits the necessary magnitude for B given our values. A typical soliton –like plot will appear:

No photo description available.
The relevant equation for a soliton embodying a dissipative term to evolve into a shock can be written:

(- v o  + c s   + v)  v / x’  -  m 2 v / x2’  +  a 3 v / x3

Where the underlined term above is the dissipative term. In a collisionless plasma such as the corona under consideration, the dissipation can be attributed to Landau damping and this may be the most common form to excite a shock associated with CMEs. For all such collisionless shocks the dissipative length scale, L << l mfp   , the mean free path.  

The preceding equation can be derived using an analog to Newton’s 2nd law of motion, viz.
a v” = F(v)  = a d2v/ dt2   =    -dF / dv

Where:   F  =  (c s   -  v o) v 2 / 2   +     v3/ 6

Is the Sagdeev Potential.


(To be continued - with problems)


[1] D.R. Nicholson:1983, Introduction to Plasma Theory,  John Wiley, p. 43

[2]J.-I. Sakai and Y. Ohsawa:1987, Space Sci. Rev., 46, 113.