Showing posts with label Alfven waves. Show all posts
Showing posts with label Alfven waves. Show all posts

Friday, September 14, 2018

Alfven Waves Revisited (2)






















We left off trying to obtain the x, y-components of velocity associated with the wave, by use of Fourier transforms. The x-component is:

w v x   –  c s     k    v x  + B  z    k m  r o   [v  z   B  x  –  v  x  B  z] = 0

The y-component is decoupled from the others (x, z) and can be written:

 w v y   -  B  x 2   k  v y m o   r o = 0

or simply:

w  = [B  x 2   m o  r o] k 2

where the quantity in brackets is the Alfven velocity or alternatively written:

v  A  = [ w/ k] =  B x  / [m o  r o1/ 2

or

v   = Bo  / [m o  r o1/ 2

since Bo    is in the x –direction


For completeness, we should be able to show the z-component equation is:

w v  -  B x  k m o  r o [v  z   B  x –  v  x  B  z] = 0

Here, let me back up and refer readers again to the basic wave equation one can obtain by getting the 2nd derivative of:

 v l / t  = -(c s  2 ) grad  p  l  / r  - 1/ m  r  [B o X Curl B 1]

One can then find the solution in terms of plane waves by assuming:

v  1 = v  1 *[exp ik.x – iwt]


for which taking the second derivative, of  v  1 with respect to t yields the original equation in w we found earlier. All of this the energetic reader should be able to work out, but most of it (after taking derivatives) reduces to brute force algebra!

For completeness, I need to note what happens when you solve the preceding (simultaneous) equations in x, and z.

w 4 +  w2 [- c s  2  k   - B x  k  / m o  r o ] + c s  2  k  4    [B x / m o  r o ] = 0

or:

w 4 -   w(c s    + v  A 2 )k  +  c s  2 v  A cos  2 (Θ) k = 0

and finally,

w = ½[( c s  2 v  A k 2  +   [(c s  2 v  A k  4  –  4 c s  2 v  A 2    cos 2 (Θ) k4 ]1/2 

Now, if one plots the preceding using for the vertical axis (c s    + v  A 2 ) and for the horizontal B o (e.g. x) one will get what is called “Friedrich’s diagram” (see sketch image). It consists of


1)A smaller “dumb bell” or figure-8 shaped graph centered at the origin. This will be for what we call “slow mode” waves

2)A single larger lobe that envelopes the smaller right lobe of the dumb bell. This will be for Alfven waves proper.

3) A circle shaped graph surrounding both 1, 2 above. This will be for what we call the “fast MHD” mode.


The critical thing to note here is that the fast mode is the only MHD wave able to carry energy perpendicular to the magnetic field. This has important ramifications for solar flares, as well as magnetospheric effects (such as the aurora). Meanwhile, the phase velocity (w /k) of the slow mode wave perpendicular to the magnetic field is always zero. In the limit where the sound speed c s  2  < < v  A 2, and the Alfven speed v  A 2<< c s  2, the slow wave disappears. (Which one can easily validate and confirm for the equation in w)

Other properties, points to note:

-the velocity perturbation v 1 is orthogonal to B o

-the wave is incompressible since DIV v 1 = ik.v 1 = 0

-the magnetic field perturbation (B 1 ) is aligned with the velocity perturbation. Since both are perpendicular to k and B o

-the current density perturbation (J 1) exists as a current perturbation perpendicular to k and B o  e.g.

1 = k X  B o

- When c s  2 <<   v  A 2 the fast mode wave becomes a compressional Alfven wave. This has a group velocity equal to its phase velocity w /k

Wednesday, September 12, 2018

Alfven Waves Revisited (1)

Alfven waves are the most important waves propagating in the solar atmosphere, as well as the Earth’s magnetosphere (underpinning the coupling between it and the ionosphere). They are important in that they efficiently carry energy and momentum along the magnetic field.

One way to get a handle (of sorts) on Alfven waves is to look at the analogy with mechanical waves – say propagating along a string put under tension. Consider the reference frame or coordinate system:


^ y
!
!----------------------------------------->x


Say x marks the direction of propagation in the above coordinate system, and y is the direction of transverse (wave) displacement. Then the vertical force component is:

Fy = - T(y/ x)

where T is the tension. Thus, just as the restoring force for a mechanical wave is the string tension T, the restoring force for an Alfven wave is the magnetic tension. This magnetic version of “tension” accelerates the plasma and is opposed by the inertia of the ions (mainly from proton masses m(p))

Now, the wave speed on a string is related to u (mass per unit length), and T such that:


v = (T/ u) 1/2

and as we can see, increasing the string tension increases the wave speed in an analogous way to what magnetic tension does for the Alfven wave. The magnetic tension analog can be expressed (as we shall see) as:


T(M) = B2 / m o


where B is the magnetic induction and mo  is the magnetic permeability for free space:
  4π x 10 -7   H/m.

In what follows we assume a uniform plasma in equilibrium, which will then be subjected to velocity disturbance or perturbation that affects all other key quantities. The treatment is kept as simple as possible (considering the complexity of the subject matter!) , and we don’t veer out of the linear domain. Nevertheless it should be stated at the outset that some details are omitted, or left as work for yourself with hints provided. In this way you will better understand and appreciate the genesis of Alfven waves. In terms of symbols, all have retained their earlier meanings (i.e. from previous plasma and math posts) and this includes the vector operators, DIV, grad, Curl etc.

Examining the origin of these waves always starts with setting out the basic equations for what we call “ideal MHD”:

1) r / = - DIV ( r v)

2)  (r v)/  t  = - DIV (r vv) – grad p + 1/ m o [(Curl B) X B]

3) B / t = Curl (v X B)

4) p/ t = - v. grad p –  g p DIV v

where  v is the fluid velocity, p the pressure, B the magnetic induction, and  g (‘gamma’)                    = - d ln p/ d ln V where V denotes volume

Now, introduce small perturbed quantities (e.g. imagine introducing a small perturbation into the plasma velocity such that vo -> v l , which will also subject the mass density, fluid pressure and magnetic field to perturbation), such that:

r = +   r l

v = v  l

B = B o + B l

p = p o + p l

Now,  substitute these back into the original ideal MHD equations to obtain:


5) r l / t = -  o  DIV v l

6) o  (v l/ t) = - grad p l + 1/ m o  [(Curl B l) X B o ]

7)  l / t  = Curl (l  X B o )

8)  ¶ p  l  / t = - g p o DIV v l


Now, divide through the 2nd equation above by  the mass density r:

¶ v l / t  = -(c s  2) grad  p  l  / r  - 1/ m  r  [B o X Curl B 1]


where ‘c s’ is the sound speed. (Note that the reader should be familiar with a vector identity also used to obtain the preceding!)

Now, using this result and the last two equations of the perturbed set, we apply Fourier transforms such for   ¶  / ¶  t and  ¶ / ¶ k to obtain:

w v l –  c s  2  (kv x )k* +  B o m  r o  [k X k* X ( v l  X B o)] = 0

where w denotes the plasma frequency, k is the wave number vector (k* the vector orientation) and the other quantities are as before. In part 2 we’ll obtain the x and y components of the velocity but readers can try in the meantime to do it themselves!

Friday, August 18, 2017

What Aspects Of The Coming Total Solar Eclipse Are Most Critical For Astrophysics?


Image showing totality in a previous  solar eclipse. The Sun is totally blocked out by the Moon - in the line of sight- and the solar corona is seen expanding from the solar limb


Graphic showing the path for Aug. 21 total eclipse - and the partial zones to either side, with the percentages of the Sun to be covered.

As we approach August 21 and the first total eclipse to appear in in the U.S. in38 years, many are trying to stay ahead of the fervid hype. Here are the only two things you need to know: 1) In most of the U.S. the eclipse will be seen as partial only and you need special glasses in order to observe the transition to maximum coverup, and 2) If you want to observe the total eclipse (see top image) you will have to travel to some place within the band of totality (see graphic).  For example, in our location (Colorado Springs) the Sun will be roughly 85 percent covered.

Here's the skinny on stats for this eclipse: Approximately 12 million Americans will be directly in the totality band shown in the graphic. Another 88 million will be within 200 miles of some place inside the totality band. Naturally then, those living in communities right in the band (e.g. Jackson Hole, WY, St. Joseph, MO) are expecting eclipse watchers to pour in and pump up their economies - buying t-shirts, special glasses and other paraphernalia.

While the basic image of totality will be the appeal for casual observers,  and even partial eclipse will occur for (and excite)  many others, this doesn't hold quite the same spell for solar physicists or astrophysicists. They look to good observations made during totality, namely to be able to record and analyze the solar corona.   This is the outermost region of the solar atmosphere, at millions of kelvins temperature and extending sometimes millions of kilometers into space.

The corona up to now has presented a mystery, especially in terms of its "inverse" increasing temperature profile considered in the context the general temperature profile is steadily decreasing from the Sun's interior to outside.   For example, the Sun's "surface" or photosphere has an effective temperature of 5,777 K .   The corona's temperature by contrast is at least 2 million degrees K.

How did we find this out?  The key breakthrough probably arrived in 1939 when astronomer Walter Grotrian found that a previously discovered spectral line (attributed to "coronium") was actually 13 times ionized iron. Since it takes enormous energy to ionize even one iron atom (meaning stripping one of its outer electrons away) this meant the iron had to be subject to enormously high temperatures.  Indeed, such ions can only exist in plasmas with temperatures between 1 million and 5 million K.

The first clue as the cause of the extraordinary hearing came with high altitude rocket flights in the 1970s, bearing x-ray telescopes. These showed  features down to a resolution of 1 arcsec or 730 km.  Careful observations showed a positive correlation between the x-ray brightness of active regions (i.e. containing complex sunspots) and plasma-filled magnetic loops, e.g.



















The above image was taken by the much more recent Atmospheric Imaging Assembly, built for the Solar Dynamics Observatory . But you easily get the idea noting the brightness (representative image taken at  211Å) and the magnetic loop structures.

Another milestone occurred when Eugene Parker predicted the hot corona must expand into space as a solar "wind". It could not remain stationary or in place with no expansion. The proof of this is left at the end for those interested.

The next question, of course, is where the energy comes from not only to heat the corona but accelerate the solar wind which is an extension of it.  We believe now that MHD or "magnetohydrodynamic" waves of some type provide an answer.   Field lines such as shown in the previous image are rooted at both ends in the Sun's convection zone where they are jostled by the churning of convective cells.  If sufficiently rapid these field line motions could generate magnetic waves capable of carrying energy upward into the loop.   It seems that Alfven waves e.g.

http://brane-space.blogspot.com/2011/04/getting-handle-on-alfven-waves-2.html

are the most effective in reaching coronal heights. (Alas, Alfven waves do not compress the ambient plasma so we still need to identify some other mechanism that can effectively transfer the wave energy to the plasma. A clue may lie with Landau damping, say associated with a beam (or 2-stream)  instability in the plasma. We have a dual Maxwellian profile.


In the region where the slope is positive (f(v)  / v > 0) there is a greater number of faster than slower particles so a greater amount of energy is transferred from particles to associated (e.g. Alfven) waves.  Since f eb contains more fast than slow particles a wave is excited.

But in Landau damping, with the slope negative (f(v)  / v < 0) the number of particles slower than the waves phase velocity exceeds the number of those that are faster. Thus, more particles gain energy from the wave than lose energy to it.

Whatever new data the eclipse generates, and there will be lots of it, we can be confident that the mysteries of the solar corona won't be solved even after it's all analyzed and multiple papers published. But I will keep readers abreast of any new findings say that appear in the next 8-12 months.

Those interested in actual scientific participation are invited to go to the Globe Observer site, where you can collect the app and then be prepared to make temperature measurements.

https://observer.globe.gov/

See also:


----------
Proof that the corona can't be stationary:


A static corona  superficially appears  to be quite reasonable but that’s why we need to test this is so. The first one to do this was Sydney Chapman. He began by first assuming the condition for hydrostatic equilibrium applied:

dp/ dr = -
r {GMs/ r2}

where G is the usual Newtonian gravitational constant, and
r defines the plasma density for the corona, while Ms is the mass of the Sun, and r the distance from the solar center:

r = n(mp)

with n the number density for protons

The coronal pressure (P) is given by:

P = 2 n T

Provided both protons and electrons are assumed to have the same temperature.

The thermal conductivity of the corona is dominated by electron thermal conductivity and takes the form:

k = ko T 5/2

for typical coronal conditions the value of
k is about 20 times the value of copper at room temperature.

The coronal heat flux density is:

q = -
k Ñ T

A static corona means heat inputs cancel heat outputs so that the divergence:

Ñ× q = 0

Assuming a spherical symmetry for the corona one can write:

1/r2 [d/dr (r2
ko T 5/2  dT/dr)] = 0

    Obviously the preceding assumptions mean there must be some distance where the coronal temperature becomes zero.

From the above equation one should be able to show:

d(T 7/2) = 7/2 (F To 5/2)/ 4
p ko d(1/r) = C d(1/r)

where C is a constant.

The integral is:

To 7/2 - T 7/2 = C[ 1/Ro - 1/r]

Now, set the temperature at infinity (T) to zero and obtain:

C = Ro To 7/2 

which fixes the total flux at:

F = 2/7   [4
p Ro ko To ]

After another step, one finds:

T(r) = To (Ro / r) 2/7

    This gives the temperature T at a distance from the Sun= r. This is based on using a defined value (say To = 2 x 106 K) at a defined distance, say Ro = 7 x 108 m.

    For example, at the Earth’s distance (r = 1.5 x 1011 m) one would find: T = 4.3 x 105 K

    This seems fine, until one examines the pressure.

    Analogous to the temperature formalism, we have, the pressure p(r) at some distance r defined by:

p(r) =


p(Ro) exp [7/5 GMs mp/ 2 T(Ro) Ro {( Ro / r)5/7 – 1}]

    Now, if one allows r to approach infinity, e.g. r
®¥ an interesting thing occurs in the equation, as we can see. That is, the denominator of the first term in the end brackets becomes so large (Ro / ¥) that the first term vanishes.

Then we are left with the expression for the pressure:

p(
¥) = p(Ro) [exp – 7k/5 * 1/ T(Ro) Ro]

where ‘k’ denotes a constant composed of all the constant quantities in the previous eqn. (G, M, mp etc)

    Substituting the given values into the above, one finds p(Ro) multiplied by a factor

exp[0] = 1

    The reason is that the exponential of a very small and negative valued magnitude
® 0

Then:

p(
¥)  »  p(Ro)

But this can’t be since the pressure of the coronal base would then be the same as the value at infinity!

    This led astrophysicists to conclude an unphysical result, and that the static coronal model couldn’t be accurate.

    If the static model were accurate, the pressure at infinity should be zero, p(
¥)   = 0, not a small finite pressure that’s effectively equal to the coronal base pressure. This finding led to the further investigations that disclosed a solar “wind” had to flow outwards from the corona.







Sunday, November 8, 2015

Rising Tone Magnetosonic Waves Detected By THEMIS Spacecraft



A previously unknown type of magnetosonic wave was detected by the THEMIS spacecraft as it traveled from the Earth's plasmasphere (green central 'blob') to the outer radiation belt (pink).

Let's accept most people would not know a magnetosonic wave from a magnetometer, or a plasma wave from any other. Despite that the use of plasma waves in actually studying near Earth plasmas has increased in importance given major solar events - such as large flares - can adversely impact the ionosphere as well as atmosphere of Earth.  Among the effects we can include the Ottawa power grid going down in 1989, and also recurrent short wave and other band blackouts - including affecting GPS and satellite TV broadcasts. All the underlying plasma phenomena come under the heading of 'space weather'.

Space weather data are assembled from across a wide spectrum as may be expected when we are trying to ascertain the effects of the Sun and solar wind on our Earth.  One of the more important diagnostics are plasma waves in the near Earth space environment.  Now, the THEMIS ('Time History of Events and Macroscale Interactions During Substorms') spacecraft has evidently detected a new kind of magnetosonic plasma wave which may play an important role in space weather forecasts.

For reference, the THEMIS spacecraft orbits in the magnetosphere near Earth's magnetic equator and collects data from magnetic storms especially near the boundary of the magnetosphere on the dayside as well as from Earth's radiation belts.

Now, two events recorded in June and August, 2010, appear to confirm the existence of a special type of magnetosonic (MS) wave, best described as "rising tone".  All plasma waves are distinguished by their phase velocities, viz. w / k  where w  is the plasma frequency and k the wave number vector. In the case of the standard magnetosonic wave we have:

w2 / k 2    =   2  (A  +  s 2 ) /  (A  +  2 )


Where A  is the Alfven velocity, v s  is the ion sound speed and c the velocity of light. In the limit of low magnetic fields, for which  v A  -> 0  the wave becomes an ordinary ion acoustic wave.  Basically the waves under review get their energy by interacting with protons trapped in Earth's magnetic field spiraling around magnetic field lines.

Historically, the frequencies of MS waves were believed to be temporally continuous, or essentially like the smoothly varying notes from a trombone player - going from one note to the next.  This indicated a simple, linear relation between MS waves and protons.

Now, Fu et al, writing in Geophysical Research Letters(2014) have incorporated the 2010 THEMIS detections to interject a possible complication. That is, a sharp rising tone in their spectrogram (see graphic) indicated a more complex, nonlinear series of interactions between MS waves and protons. This would be more analogous to a flute player performing a series of runs and trills.

To be sure, space and plasma physicists have observed rising tone phenomena in other plasma waves including electromagnetic ion cyclotron waves. In these waves it is the web of forces between  the particles which create currents that enhance the wave's frequency. However, this nonlinear behavior bas never before been observed in MS waves.

It remains to be seen just how useful this new form will be in space weather forecasts but for many of us in the space sciences we look forward to productive research and more refined forecasts

Friday, June 27, 2014

A Cavity Resonator Model Applied to Solar Loops and Flare Triggers (1)



The problem of identifying a unique trigger for solar flares has been pursued for over 4 decades, but with little to show for it. In this post I examine a possible approach that might be productive, especially after higher resolution images become available from a planned solar telescope.


1.     Background to Cavity Resonator Approaches:


On the Sun itself, the 5-minute oscillations, more figuratively described from time to time as a “ringing of the photosphere” were first tied to waves trapped in a resonant cavity by Schatzman, 1956[2]. The gist of the model is that the upper and lower cavity boundaries reflect waves into the cavity and thereby engender a standing wave, which may either be acoustic or gravity.

 Properties of generic coronal cavity resonators were elucidated by Hollweg (1984)[3] whereby a coronal loop is treated as three relatively disjoint regions, separated by discontinuities. The standing waves produced are invoked to account for energy dissipation and heating in the corona. Some aspects of Hollweg’s model (e.g. reflection properties of Alfven waves, quality factor Q, relation to wave number vectors) are also employed in my own resonator model for solar flare inception.

 Where I diverge from both the (e.g. Federov et al) auroral cavity resonator and the coronal one proposed by Hollweg is that in this flare trigger model I adopt a dual resonator for a given compact flare loop. The basic sketch is shown above. Thus, to trigger a specific (e.g. compact) flare the conditions must be such that the Q-(quality) value in each resonator reinforce the other. In my generic model, I include a small coronal arch cavity resonator with some resemblance to Hollweg’s and a large scale loop resonator which depends on the oscillations arising from magnetic (Alfven) waves in combination with the loop’s twist (and associated kink instability)

 The difference is that I go into much more detail to incorporate ex-post facto data into my model to show how the magnitudes of the changing physical quantities vary not only in the corona in its pre-flare state, but during the flare as well. As an application ansatz, dual resonators in the electrical engineering setting often use closed-loop resonators in order to shift down the original resonator and arrive at a very small structure (e.g. Collado et al, 2007)[4]. In the flare trigger model context this would be the kernel or coronal loop apex resonator.  In the engineering context, mirrors are sometimes employed to linear cavities (analogous to the extended coronal loop with its primary cavity at the apex) to obtain simultaneous dual wavelength oscillations. It is precisely within the scope of these dual oscillations that the flare trigger can be conceived – e.g. for specific cases when a dual resonance is achieved and with it the maximum instability.

2.     Motivation:

Having established that a hybrid flare model is the most plausible one to approach the 1B/M4 flare of November 5, 1980, I now single out the key feature for the flare trigger. Before proceeding, let us inspect the gestalt for what this article is all about, as depicted in the schematic below:
No photo description available.

The diagram depicts the generic inputs and processes entering the hybrid flare model (in the main rectangle), which includes components for R (reliability statistics), H (helicity considerations) and Poisson statistics. Thus, the hybrid model seeks to reconcile all of these, as well as recognizing the inputs from two paradigms: the E-J and the B-v.  For example, the B-v paradigm and its assumptions figure more prominently in the reliability analyses as well as the magnetic helicity. (E.g. see: http://brane-space.blogspot.com/2010/10/look-at-magnetic-helicity.html )  The E-J paradigm factors more into the basis for Poisson variations based on the manner in which the current densities (J) arise and how the E-field is generated.  The putative flare trigger attempts to make use of all of these.

What is desired is a model that approximately replicates the event sequence for the region AR2776 such that the flares occurring conform to the average Poisson  activity:  l (av) =  2.6 x 10-5 s-1  and the length, resonance variations described above imply a twisted cavity (dual) resonator over the region defined by (l1 + ℓ1^ +  ℓ2^   +  xi ).  (See e.g. my post of June 22nd: http://brane-space.blogspot.com/2014/06/quantifying-solar-loop-oscillations.html

Where:   0 <   xi  <  1.1 x 106 m

The basic geometry is shown in Fig. 1 (top) in the region of the apex, and primary coronal cavity. It is assumed that with compressional Alfven waves the loop aperture can vary, from a1 to the outer radius taken as r1. This variation could well account for the uncertainty in source-kernel dimensions:


The electric field E(z) shown in Fig. 1 is described according to:


E(z)  =   Eo cos w(t – z/ vp)

Where Eo  denotes the uniform (non-varying field magnitude) and  vp =  c sin(J)

is the phase velocity with J the pitch angle of the twist component for relative helicity (H(R) [T]).

The model works via the basic loop changing its effective resonator length (for which there is an associated resonator angular frequency wo) and twist (F(r)).  Radial surfaces (rs < r)  form in the loop apex (small resonant cavity) for E-field resonating corrections modeled after the J o (ar)   Bessel function.  The J o (ar)   induced field in turn generates azimuthal corrections in the axial B-field that alters the twist of the gross loop.

Thus, the twist dependence is (see also http://brane-space.blogspot.com/2013/04/looking-at-bessel-functions-applications.html):


F(r) =   L J1(ar)/ r J o (ar)   =    L B j (r) / r Bz (r)     = L E z (r)  /  r E j (r)

Such that: E z (r)  ® B j (r) ® E1 z (r) ®   B1 j (r) ® E2 z (r) ®   B2 j (r)


E j (r) ® B z (r) ®  E1 j (r) ® B1 z (r) ®  E2 j (r) ® B2 z (r) . . . . En j (r) ® Bn z (r)

From the secondary, tertiary etc. fields inner nested radii aij are generated which conform to ratios related to the Bessel functions. The key point of the mutually generated fields is that they operate according to a positive feedback which ultimately incepts a resonance condition and explosive release of energy. The radii in turn can be used to obtain wave modes associated with a given oscillation period for the resonator. The relative E-field strengths successively generated for the ideal cavity coronal resonator are defined by the Bessel series:

Jm (x) = (1/ 2m m!) xm [1 -  x 2/ 22 1! (m + 1)  +  x4/ 242! (m + 1) (m + 2) -  ….

.(-1)j x2j / 2 2j j! (m + 1) (m + 2)……(m + j) +  …]

Which may be simplified for the m = 0 case  to:


J0 (x) = 1 -  (x/2)2 + 1/(2!)2 (x/4)4 – 1/ (3!)2 (x/2)6 +  .




Where x for the E-field is defined x =  Ö mo Öo w r  =   2.405



For which a key cut-off radius is defined at the surface rs  = r .  Other surfaces  (si) may be defined for zeros of J0 (x).  Meanwhile, coronal loop oscillation periods and emergence have been well explicated by a number of authors (e.g. Edwin and Roberts, 1983, op. cit., Andries et al, 2005[5])
(More to come)



[1] E.N. Fedorov, V.A. Pilipenko, M.J. Engebretson, and T. J. Rosenber: 2004, ‘Alfven Wave Modulation of the Auroral Acceleration Region’, in Earth Planets Space, 56, p. 649.
[2] E. Schatzman: 1956, Ann. Astrophysics, 19, 45.
[3] Joseph V. Hollweg, Solar Phys., 91, 269, 1984
[4] C. Collado, J. Pozo, J. Mateu and J.M. O’Callaghan: 2007, European Microwave Week.
[5] J. Andries, M. Goosens, J.V. Hollweg, I. Arregui and T. Van Doorselaere,: 2005,  Astronomy & Astrophysics,  430, 1109.