1 / 72

基于最小能量耗散原则的太阳日冕磁场非无力场外推

基于最小能量耗散原则的太阳日冕磁场非无力场外推. Qiang Hu ( 胡强) CSPAR, University of Alabama in Huntsville, USA qh0001@uah.edu. Acknowledgement. Collaborators :. D. Shaikh C. W. Smith/N.F. Ness W.-L. Teh Bengt U. Ö. Sonnerup A.-H. Wang/S.T. Wu Vasyl Yurchyshyn Gary Zank. Debi P. Choudhary B. Dasgupta

eric-perry
Download Presentation

基于最小能量耗散原则的太阳日冕磁场非无力场外推

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 基于最小能量耗散原则的太阳日冕磁场非无力场外推基于最小能量耗散原则的太阳日冕磁场非无力场外推 Qiang Hu(胡强) CSPAR, University of Alabama in Huntsville, USA qh0001@uah.edu

  2. Acknowledgement Collaborators: D. Shaikh C. W. Smith/N.F. Ness W.-L. Teh Bengt U. Ö. Sonnerup A.-H. Wang/S.T. Wu Vasyl Yurchyshyn Gary Zank Debi P. Choudhary B. Dasgupta Charlie Farrugia G.A. Gary Yang Liu Dana Longcope Jiong Qiu G. M. Webb NASA grants NNG04GF47G, NNG06GD41G, NNX07AO73G, and NNX08AH46G (data provided by various NASA/ESA missions, and ground facilities; images credit: mostly NASA/ESA unless where indicated)

  3. What is extrapolating coronal magnetic field and why? -- obtain three-dimensional magnetic field in a finite volume from observed magnetic field data only available on bottom boundary • Generally hard to directly measure • Critical to energy release of eruptive events • Sophisticated models needed • …

  4. Coronal Magnetic Field Extrapolation (2D3D) • One existing simple model, variational principle of minimum energy (e.g., Taylor, 1974; Freidberg, 1987): Linear force-free field (LFFF,const) Or, Nonlinear FFF ( varies) However, Amari and Luciani (2000), among others, showed by 3D numerical simulation that in certain solar physics situation, …, the final “relaxed state is far from the constant- linear force-free field that would be predicted by Taylor’s conjecture” …, and suggested to derive alternative variational problem.

  5. Gary (2009) An alternative... • Principle of Minimum Dissipation Rate (MDR): the energy dissipation rate is minimum. • (Montgomery and Phillips,1988; Dasgupta et al. 1998; Bhattacharyya and Janaki, 2004) • (Several extended variational principles of minimum energy (Mahajan 2008; Turner 1986) yield solutions that are subsets to the above) • Simple Examples: • Current distribution in a circuit • Total ohmic dissipation is minimum • Velocity profile of a viscous liquid flowing through a duct • Total viscous dissipation is minimum

  6. constraint ( = i, e) variational problem Formulation of the variational problem (for an open system with external drive, or helicity injection) • The generalized helicity dissipation rate is time-invariant. • From the MDR principle, the minimizer is the total energy dissipation rate

  7. Euler-Lagrange equations Eliminating vorticity in favor of the magnetic field

  8. New Approach • For an open system with flow, the MDR theory yields (Bhattacharyya et al. 2007;Hu et al., 2007; Hu and Dasgupta, 2008, Sol. Phys.)   Take an extra curl to eliminate the undetermined potential field , one obtains (a1=-1-3, b1= 13 )

  9. (5) (7) (8) • Equations (5),(7) and (8) form a 3rd order system. It is guaranteed invertible to yield the boundary conditions for each Bi, given measurements of B at bottom boundary, provided the parameters, 1 ,2 and 3 are distinct.

  10. Above equations provide the boundary conditions (normal components only at z=0) for each LFFF Bi, given Bat certain heights, which then can be solved by an LFFF solver based on FFT (e.g, Alissandrakis, 1981). • One parameter, 2 has to be set to 0. The parameters, 1 and 3, are determined by optimizing the agreement between calculated (b) and measured transverse magnetic field at z=0, by minimizing <0.5? Measurement error + Computational error Measurement error

  11. Reduced approach: choose B2=cB’, proportional to a reference field, B’=A’, and B’n=Bn, such that the relative helicity is ABdV-A’B’dV, with A=B1/1+B3/3+cA’. • B’=0, Bz’=Bz, at z=0 Only one layer of vector magnetogram is needed. And the relative helicity of a solar active region can be calculated. (As a special case, B2=0, in Hu and Dasgupta, 2006)

  12. k=0 Reduced approach: Obtain E(k) and N If E(k) <  k=k+1 Y End • Iterative reduced approach: transverse magnetic field vectors at z=0 (En=0.32): n n Y k>kmax

  13. Test Case of Numerical Simulation Data (Hu et al. 2008, ApJ) (b) our extrapolation (12812863) • “exact” solution • (Courtesy of Prof. J. Buechner)

  14. Transverse magnetic field vectors at z=0 (En=0.30): • Figures of merit (Hu et al., 2008, ApJ): Energy ratios

  15. Integrated current densities along field lines: (a) exact (b) extrapolated J_perp  0 J_para

  16. E/Ep_pre=1.26 E/Ep_post=1.30 (Data courtesy of M. DeRosa, Schrijver et al. 2008; Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC (UK) as international partners. It is operated by these agencies in co-operation with ESA and NSC (Norway).

  17. Main Features • More general non-force free (non-vanishing  currents); • Better energy estimate • Fast and easy (FFT-based); • Make it much less demanding for computing resources • Applicable to one single-layer measurement (Hu et al. 2008, 2009) • Applicable to flow

  18. A homotopy formula for vector magnetic potential (based on ):

  19. “In topology, two continuousfunctions … if one can be "continuously deformed" into the other, such a deformation being called a homotopy between the two functions.” -- Wikipedia (Berger, M.)

  20. Relative magnetic helicity via homotopy formula: Bz(x,y) r

  21. Relative magnetic helicity via homotopy formula: B in a 3D volume (e.g., see Longcope & Malanushenko, 2008) r

  22. A simplified vector potential for a potential field? • Multi-pole expansion of a potential field: • For each 2k-th pole, B(k), • Dipole: • Related to spherical harmonic expansion, for example. • For MDR-based extrapolation:

  23. Outlook • Validate and apply the algorithm for one-layer vector magnetograms • Validate the theory – proof of MDR by numerical simulations • Global non-force free extrapolation Stay tuned!

  24. Test with Data-Driven MHD Simulation A. Wang, S. T. Wu et al.

  25. MHD states:

  26. Dec. 12-13 2006 Flare and CME (Schrijver et al. 2008) (SOHO LASCO CME CATALOG http://cdaw.gsfc.nasa.gov/CME_list/) (KOSOVICHEV & SEKII, 2007)

  27. Reduced approach: transverse magnetic field vectors at z=0 (En=0.7-0.9): Measured Computed (Data courtesy of M. DeRosa, Schrijver et al. 2008; Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC (UK) as international partners. It is operated by these agencies in co-operation with ESA and NSC (Norway).

  28. Post-flare case: En=0.28

  29. En(1,3):

  30. A different variational principle Suggests a minimizer for our problem • R. Bhattacharyya and M. S. Janaki,Phys. Plasmas 11, 615 (2004). • D. Montgomery and L. Phillips, Phys. Rev. A 38, 2953, (1988). • B. Dasgupta, P. Dasgupta, M. S. Janaki, T. Watanabe and T. Sato, Phys. Rev. Letts, 81, 3144, (1998) • Principle of minimum dissipation rate (MDR) In an irreversible process a system spontaneously evolves to states in which the energy dissipation rate is minimum.

  31. L Onsager, Phys. Rev, 37, 405 (1931)I. Prigogine, Thermodynamics of Irreversible Processes, Wiley (1955)------------------------------------------------------------------------------------------------------ A theorem from irreversible thermodynamics: Principle of Minimum Entropy Production “The steady state ofan irreversible process, i.e., the state in which thermodynamics variables are independent of time, is characterized by a minimum value of the rate of Entropy Production” Rate of Minimum Entropy Production is equivalent to Rate of Minimum Dissipation of Energy in most cases.

  32. Numerical simulation (Shaikh et al., 2007; NG21A-0206 ) showed the evolutionof the decay rates associated with the turbulent relaxation, viz, Magnetic Helicity KM, Magnetic Energy WM and the Dissipation Rate R. KM = WM = R=

  33. constraint ( = i, e) variational problem Formulation of the variational problem (for an open system with external drive, or helicity injection) • The generalized helicity dissipation rate is time-invariant. • From the MDR principle, the minimizer is the total energy dissipation rate

  34. Euler-Lagrange equations Eliminating vorticity in favor of the magnetic field

  35. Summary of Procedures magnetohydrodynamic

  36. Analytic Test Case: non-force free active region model given by Low (1992) Top View:

  37. A real case: Active Region (AR)8210 (preliminary) Imaging Vector Magnetograph (IVM) at Mees Solar Observatory (courtesy of M. Georgoulis) (Choudhary et al. 2001)

  38. En distribution:

  39. GS Reconstruction • One-fluid Magnetohydrostatic Theory • 2 ½ D: Bz 0 • Co-moving frame: DeHoffmann-Teller (HT) frame • No inertia force Grad-Shafranov (GS) Equation (A=Az): A • Bt = 0 2 Pt(A)=p(A)+Bz(A)/20

  40.  GS Reconstruction Technique • Find z by the requirement that Pt(A) be single-valued • Transform time to spatial dimensions via VHT, and calculate A(x,0), • 3. Calculate Pt(x,0) directly from measurements. • 4. Fit Pt(x,0)/Bz(x,0) vs. A(x,0) by a function, Pt(A)/Bz(A). • A boundary, A=Ab, is chosen. • 5. Computing A(x,y) by • utilizing A(x,0), Bx(x,0), and GS equation. x: projected s/c path o:inbound Am ^ *:outbound

  41.    Pt(A) • Finding z axis by minimizing residue of Pt(A): • Residue=[∑i(Pt,i – Pt,i )2] • /|max(Pt)-min(Pt)| Pt(x,0) 1 2 1st 2nd o: 1st half *: 2nd half i=1…m A(x,0) Enumerating all possible directions in space to find the optimal z axis for which the associated Residue is a minimum. A residue map is constructed to show the uniqueness of the solution with uncertainty estimate.

  42. GS Solver:

  43. Multispacecraft Test of GS Method Cluster FTEs (from Sonnerup et al., 2004; see also Hasegawa et al. 2004, 2005, 2006)

  44. Introduction • Grad-Shafranov (GS) equation: p=jB in 2D • GS technique: solve GS equation using in-situ data, 1D2D • (e.g., Sonnerup and Guo, 1996; Hau and Sonnerup, 1999; Hu and Sonnerup, 2000, 2001, 2002, 2003; Sonnerup et al. 2006)

  45. Small-scale flux ropes in the solar wind (Hu and Sonnerup, 2001)

  46. Features of the GS Reconstruction Technique: • Fully 2 ½ D solution (less fitting) • Self-consistent theoretical modeling; boundary definition (less subjective) • Utilization of simultaneous magnetic and plasma measurements; Non-force free • Adapted to a fully multispacecraft technique (Hasegawa et al. [2004]) • • Limitations (diagnostic measures): • 2D, uncertainty in z (the quality of Pt(A) fitting, Rf) • 2D  Pt(A), Pt(A)  2D • Time stationary (quality of the frame of reference) • Static (evaluating the residual plasma flow) • Numerical errors limit the extent in y direction (rule of thumb: |y||x|, y«x) ?

  47. ACE-Wind comparison : Wind data o  : Predicted

More Related