Showing posts with label coronal loop. Show all posts
Showing posts with label coronal loop. Show all posts

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.

Saturday, March 29, 2014

Solving the Primary Riddle of Magnetic Reconnection in Solar Flares
















A letter in Physics Today (January. p. 8, 'The Search for Magnetic Reconnection in Solar Flares' ) summarized the problem of establishing the theory of  magnetic reconnection associated with solar flares.  Peter Foukal noted that we need a direct observational method, as opposed to mere tantalizing clues. He used as an analogy the nature of sunspots and how a century ago the vortex -like structures (e.g. penumbral whorls)  suggested magnetic fields. However, the actual magnetic field existence wasn't really established until George Ellery Hale used the Zeeman splitting of a spectral line to actually detect spots' magnetic fields. E.g.

Where the left image shows the line-centered sunspot for which the Zeeman effect (right image) is detected and measured. The greater the spectral line splitting the greater the magnitude of the associated magnetic field.

In a similar way, studies of coronal morphology and motions have left solar physicists in a similar position today. There are tantalizing clues for reconnection, but no direct evidence. Foukal cites an earlier article (Physics Today, Sept., 2013) by Johanna Miller that mentions a key signature of reconnection appears to be an intense motional electric field. This would be defined:

E = vB

where v is the plasma velocity, B the magnetic induction, and ℓ  a length element. Foukal notes that "some evidence for such fields has been observed using a relatively simple polarimeter to measure the Stark effect", and that comparison of such observations with the recently developed three dimensional models that Miller describes "might finally enable us to decipher the role of reconnection in solar flares." (The graphic attached shows one such 3-D model for a solar coronal loop with an assumed potential to spawn solar flares via magnetic reconnection).

Key to this is settling, establishing "whether the potential drops expected with reconnection occur across solar structures that produce detectable emission in Stark - affected hydrogen lines."

Make no mistake this would be an enormous advance given that so many current flare parameters, including for energy and power, are more in the way of indirect estimates.

For example, estimation of flare energy from the soft x-ray record is fairly straightforward and entails multiplying the SXR flux (left axis) of the “half-power” points by the time duration (horizontal axis) then by the recorded flare area in square meters.   A segment of such a record is shown below:
No automatic alt text available.
Consider the flare occurring at 04h 30m UT on Nov. 6,  1980 with an estimated half-power point flux of F = 10 -5 W m –2  and a duration t » 3h » 10800s. If the associated Solar Geophysical Data records show it has a flare area of 10 23  m 2    then:

Flare energy » (10 -5 W m –2) (10800s)( 10 23  m 2 )
 »  1.1  x 10 23  J


 The power of the flare is then estimated, by dividing the energy by the time of duration, so:


Power = Energy / time = (1.1  x 10 23  J)/ 10800s


 Flare power » 10 19 W


 Other quantitative estimates are possible but they require that specific models be applied (say for ‘double layers’) or that other ancillary measurements be known (say the magnetic flux, j =  BA, where A is the area of the spot, say, and B the magnetic induction).


  If the electrical resistance R, associated with a pre-flare system is known, then the current I can be estimated and the region assessed for a flare. Thus, if R is known and we know (from basic physics) that:


P =  Io 2 R 


Then the pre-flare current :   Io = [P/ R]1/2


If Io   is then known it is possible to obtain the voltage drop (V(t)) in a double layer (assuming that model applies) since:

 P =  I V(t)

If both V(t) and I are  known, it is feasible to obtain the change in the system’s inductance (dL/dt). Conversely, if we can estimate dL/dt (say using the Spicer (1980) method), and can also estimate I, then we can also estimate the magnitude of the potential drop:

V(t) =  I dL/dt  

For example, say the pre-flare current is estimated to be:  Io »   1.1 x 10 10  A . If  we can also estimate the change in self-inductance, say  dL/dt   »  0.1  H/s then the potential drop would be estimated as:

V(t)  »   1.1 x 10 10  A (0.1  H/s)  »     1.1 x 10 9  V


 Which is a reasonable value though one would like to know its time evolution up the instant of flare eruption.

Clearly, as Foukal indicates, it would be far more satisfactory to be able to obtain a direct measurement of V(t) - say using the Stark effect*.  Fortunately, we may not have too much longer before we can rise above mere estimates of potential drops. Foukal references (ibid.), a "state of the art electrograph installed, for example, on the Advanced Technology Solar Telescope" which could "open the door to a more sensitive study of motional electric fields."  Provided the NASA budget isn't further  cut, this would be a huge boon.

If we can more accurately identify the actual measured (Stark-effect) related locations in coronal loops where flare -associated potential drops occur, we can more confidently predict them - including large, geo-effective flares. Also the type that can cause CMEs or coronal mass ejections. This would have enormous import for our telecommunications, as well as aircraft navigation - which can be disrupted, adversely affected, by large solar flares and particularly those with CMEs.

*See, e.g.