AMR Simulations of Magnetohydrodynamic Problems by the CESE Method in Curvilinear Coordinates

The objective of this paper is to present new extensions of the space – time conservation element and solution element (CESE) method for simulations of magnetohydrodynamic (MHD) problems in general curvilinear coordinates by using an adaptive mesh refinement (AMR) grid system. By transforming the governing MHD equations from the physical space (x,y,z) to the computational space (ξ,η,ζ) while retaining the form of conservation, the CESE method is established for MHD in the curvilinear coordinates. Utilizing the parallel AMR package PARAMESH, we present the first implementation of applying the AMR CESE method for MHD (AMR-CESE-MHD) in both Cartesian and curvilinear coordinates. To show the validity and capabilities of the AMR-CESE-MHD code, a suite of numerical tests in two and three dimensions including ideal MHD and resistive MHD are carried out, with two of them in both Cartesian and curvilinear coordinates. Numerical tests show that our results are highly consistent with those obtained previously by other authors, and the results under both coordinate systems confirm each other very well.     

Cross Helicity and Turbulent Magnetic Diffusivity in the Solar Convection Zone

In a density-stratified turbulent medium, the cross helicity 〈u′⋅B′〉 is considered as a result of the interaction of the velocity fluctuations and a large-scale magnetic field. By means of a quasilinear theory and by numerical simulations, we find the cross helicity and the mean vertical magnetic field to be anti-correlated. In the high-conductivity limit the ratio of the helicity and the mean magnetic field equals the ratio of the magnetic eddy diffusivity and the (known) density scale height. The result can be used to predict that the cross helicity at the solar surface will exceed the value of 1 gauss km s⁻¹. Its sign is anti-correlated to that of the radial mean magnetic field. Alternatively, we can use our result to determine the value of the turbulent magnetic diffusivity from observations of the cross helicity.     

Numerical Study for the Bursty Nature of Spontaneous Fast Reconnection

In this paper, spontaneous fast reconnection in a neutral current sheet, which is initially perturbed by a localized resistivity, is studied by the newly developed Space-Time Conservation Element and Solution Element (CESE) method. After the initial perturbation is switched off, an anomalous resistivity is allowed to occur if a threshold of the local electron-ion drift velocity is exceeded. For a given threshold value, the amount of the reconnected magnetic flux introduced by the initial perturbation is very crucial for the onset of the anomalous resistivity. The numerical results indicate that fast reconnection can develop self-consistently with slow shocks extending between the diffusion region and a large-scale plasmoid-like structure, which is pushed forward by the reconnection outflow. A Petschek-like configuration is then built up, but it can not be sustained as a quasi-steady state. In fact, during the reconnection evolution, the diffusion region undergoes an elongation process so that after the dynamic process is nonlinearly saturated secondary tearing is subject to occur at the center of the system. This leads to enhanced and time-dependent reconnection. The reconnection evolution is further studied in various physical situations, also confirming the bursty nature of the spontaneous fast reconnection mechanism.     

Transition-Region Explosive Events: Reconnection Modulated by p-Mode Waves

Transition-region explosive events (TREEs) have long been proposed as a consequence of magnetic reconnection. However, several critical issues have not been well addressed, such as the location of the reconnection site, their unusually short lifetime (about one minute), and the recently discovered repetitive behaviour with a period of three to five minutes. In this paper, we perform MHD numerical simulations of magnetic reconnection, where the effect of five-minute solar p-mode oscillations is examined. UV emission lines are synthesised on the basis of numerical results in order to compare with observations directly. It is found that several typical and puzzling features of the TREEs with impulsive bursty behaviour can only be explained if there exist p-mode oscillations and the reconnection site is located in the upper chromosphere at a height range of around 1900 km < h < 2150 km above the solar surface. Furthermore, the lack of proper motions of the high-velocity ejection may be due to a rapid change of temperature along the reconnection ejecta.     

The Correlation between the Magnetic and Velocity Fields on the Full Solar Disk

A possible correlation between the magnetic and velocity fields has been analyzed based on the SOHO/MDI magnetograms and Dopplergrams. It is found that the observed large-scale weak magnetic field (weaker than 50 G (gauss)) is correlated with the velocity statistically. The curves of u⋅b with latitude, where u and b are the velocity and magnetic fields in a rectangular region (±15^○ in longitude, ±45^○ in latitude) on the Sun, show the same patterns in the years 2000, 2004, and 2007. The patterns indicate that u and b are positively correlated near the equator but are anti-correlated at the middle latitudes. For a strong magnetic field between 50 G and 3000 G, the curves of u⋅b with latitude show the same tendencies at the middle latitudes. Near the equator, however, the slope of the curve is positive in 2000 and is negative in 2004 and 2007. In addition, we give an estimation for the amplitude of the cross helicity h _χ (h_c=[`(u·b)]h_{\chi}=\overline{\mathbf{u}\cdot\mathbf{b}}) inferred from the MDI data, which is of the order of 10³ G m s⁻¹ near the center of the solar disk.     

A New Look at Mode Conversion in a Stratified Isothermal Atmosphere

Recent numerical investigations of wave propagation near coronal magnetic null points (McLaughlin and Hood: Astron. Astrophys. 459, 641, 2006) have indicated how a fast MHD wave partially converts into a slow MHD wave as the disturbance passes from a low-β plasma to a high-β plasma. This is a complex process and a clear understanding of the conversion mechanism requires the detailed investigation of a simpler model. An investigation of mode conversion in a stratified, isothermal atmosphere with a uniform, vertical magnetic field is carried out, both numerically and analytically. In contrast to previous investigations of upward-propagating waves (Zhugzhda and Dzhalilov: Astron. Astrophys. 112, 16, 1982a; Cally: Astrophys. J. 548, 473, 2001), this paper studies the downward propagation of waves from a low-β to high-β environment. A simple expression for the amplitude of the transmitted wave is compared with the numerical solution.     

On discontinuous plasma flows in the vicinity of reconnecting current sheets in solar flares

The question about the interpretation of numerical experiments on magnetic reconnection in solar flares is considered. A correspondence between the standard classification of magnetohydrodynamic discontinuities and the parameters characterizing the mass flux through a discontinuity and the magnetic field configuration has been established within a classical formulation of the problem on discontinuous magnetohydrodynamic flows. A pictorial graphical representation of the relationship between the angles of the magnetic field vector relative to the normal to the discontinuity plane on both its sides has also been found. The relations between the parameters of a two-dimensional discontinuous flow have the simplest form in a frame of reference where the magnetic field lines (B) are parallel to the matter velocity (u)—the deHoffmann-Teller frame. The question about the transformation of the magnetic field configuration when passing to a “laboratory” frame of reference where (v · B) ≠ 0, i.e., an electric field is present, is considered in this connection. The result is applied to the analytical solution of the problem on the magnetic field structure in the vicinity of a reconnecting current sheet obtained previously by Bezrodnykh et al. The regions of nonevolutionary shocks are shown to appear near the endpoints of a current sheet with reverse currents.     

Evaluating Mean Magnetic Field in Flare Loops

We analyze multiple-wavelength observations of a two-ribbon flare exhibiting apparent expansion motion of the flare ribbons in the lower atmosphere and rising motion of X-ray emission at the top of newly-formed flare loops. We evaluate magnetic reconnection rate in terms of V _r B _r by measuring the ribbon-expansion velocity (V _r) and the chromospheric magnetic field (B _r) swept by the ribbons. We also measure the velocity (V _t) of the apparent rising motion of the loop-top X-ray source, and estimate the mean magnetic field (B _t) at the top of newly-formed flare loops using the relation 〈V _t B _t〉≈〈V _r B _r〉, namely, conservation of reconnection flux along flare loops. For this flare, B _t is found to be 120 and 60 G, respectively, during two emission peaks five minutes apart in the impulsive phase. An estimate of the magnetic field in flare loops is also achieved by analyzing the microwave and hard X-ray spectral observations, yielding B=250 and 120 G at the two emission peaks, respectively. The measured B from the microwave spectrum is an appropriately-weighted value of magnetic field from the loop top to the loop leg. Therefore, the two methods to evaluate coronal magnetic field in flaring loops produce fully-consistent results in this event.     

Linear Analysis of Axial Sheared Flow in Astrophysical Jets

A linear analysis of axial sheared flow in magnetohydrodynamic (MHD) jets with helical magnetic fields is presented. A linearized set of ideal MHD equations allows the investigation of plasmas with both magnetic shear and flow shear included in the equilibrium profile. These equations are integrated numerically by following the linear development in time of an initial seed perturbation. Global instability growth rates are obtained after the numerical solution converges to the fastest growing mode. It is shown that axial sheared flow reduces the growth of current-driven instabilities in plasma jets with constant magnetic pitch P = rB _z /B _θ.     

Active-Region Magnetic Structure Observed in the Photosphere and Chromosphere

The full magnetic vector has been measured in both the photosphere and chromosphere across sunspots and plage in NOAA Active Region 8299. We investigate the vertical magnetic structure above the umbral, penumbral and plage regions using quantitative statistical comparisons of the photospheric and chromospheric magnetic data. The results include: (1) a general decrease in average magnetic flux density with height; (2) the direct detection of the superpenumbral canopy in the chromosphere; (3) values for dB/dz which are consistent with earlier investigations when derived from a straight difference between the two measurements, but which are somewhat small when derived from the B=0 condition, (4) a monolithic structure in the umbrae which extends well into the upper chromosphere, with a very complex and varied structure in penumbrae and plage, as evidenced by (5) a uniform magnetic scale height in the umbrae with an abrupt jump to widely varying scale heights in penumbral and plage regions. Further, we find (6) evidence that field extrapolations using the photospheric flux as the boundary may not agree with expectations or with observed coronal structures as well as those which use the chromospheric magnetic flux as the extrapolation starting point.     

Sweeping-magnetic-twist mechanism for the acceleration of jets in the solar atmosphere

A magnetodynamic mechanism for the acceleration of jets in the solar atmosphere (surges, Brueckner's EUV jets, and so on) is proposed, and a 2.5-dimensional MHD simulation is performed to show how this mechanism operates in the situation of the chromosphere-corona region of the solar atmosphere. It is seen from the result of simulation that together with the release of the magnetic twist, e.g., into a reconnected open flux tube, the mass in the high density twisted loop is driven out into the open flux tube due both to the pinch effect progressing with the packet of the magnetic twist into the open flux tube, and to the j × B force at the front of the packet of the unwinding twist in the off-axis part of the tube. The former, the progressing pinch, is accompanied by an accelerated hot blob, while the latter, the unwinding front of the magnetic twist, drives a cool cylindrical flow, both with velocities of the order of the local Alfvén velocity. One of the characteristic properties of the jet in our model is that the jet, consisting of hot core and cool sheath, has a helical velocity field in it, explaining the thus-far unexplained observed feature.The sudden release of the magnetic twist into an open flux tube is most likely to be due to the reconnection between a twisted loop and the open flux tube. The mass is driven out in the relaxation process of the magnetic twist from the twisted loop to the open flux tube.     

Evaluation of a Selected Case of the Minimum Dissipative Rate Method for Non-Force-Free Solar Magnetic Field Extrapolations

The minimum dissipative rate (MDR) method for deriving a coronal non-force-free magnetic field solution is partially evaluated. These magnetic field solutions employ a combination of three linear (constant-α) force-free-field solutions with one being a potential field (i.e., α=0). The particular case of the solutions where the other two α’s are of equal magnitude but of opposite sign is examined. This is motivated by studying the SOLIS (Synoptic Optical Long-term Investigation of the Sun (SOLIS), a National Solar Observatory facility) vector magnetograms of AR 10987, which show a global α value consistent with an α=0 value as evaluated by (∇×B) _z /B _z over the region. Typical of the current state of the observing technology, there is no definitive twist for input into the general MDR method. This suggests that the special α case, of two α’s with equal magnitudes and opposite signs, is appropriate given the data. Only for an extensively twisted active region does a dominant, nonzero α normally emerge from a distribution of local values. For a special set of conditions, is it found that (i) the resulting magnetic field is a vertically inflated magnetic field resulting from the electric currents being parallel to the photosphere, similar to the results of Gary and Alexander (Solar Phys. 186:123, 1999), and (ii) for α≈(α _max/2), the Lorentz force per unit volume normalized by the square of the magnetic field is on the order of 1.4×10⁻¹⁰ cm⁻¹. The Lorentz force (F _L) is a factor of ten higher than that of the magnetic force d(B ²/8π)/dz, a component of F _L. The calculated photospheric electric current densities are an order of magnitude smaller than the maximum observed in all active regions. Hence both the Lorentz force density and the generated electric current density seem to be physically consistent with possible solar dynamics. The results imply that the field could be inflated with an overpressure along the neutral line. However, the implementation of this or any other extrapolation method using the electric current density as a lower boundary condition must be done cautiously, with the current magnetography.     

Coronal magnetic field equilibrium with discrete flux sources

A model of the equilibrium structure of the coronal magnetic field is developed, taking account of the fact that field lines are rooted in the photosphere, where field is concentrated into isolated flux tubes. The field is force-free, described by ×B = B, with constant; this field has special physical significance, being the state of mininum energy after small-scale reconnections, and is also mathematically convenient in that the principle of superposition can be used to construct complex geometries. First a model of a single loop is presented, with a flux source and sink pair at the photosphere; both point flux tubes and finite radius flux tubes are considered. Then more complex topologies with multiple sources and sinks are investigated. It is shown that significant topology changes arise for different values of, indicating the possibility that there can be energy changes through magnetic reconnection if the field evolves ideally and then relaxes to a linear state.     

MHD waves and instabilities of a temperature-anisotropic plasma in the solar corona as a source of its heating

The MHD instabilities of a temperature-anisotropic coronal plasma are considered. We show that aperiodic mirror instabilities of slow MHD waves can develop under solar coronal conditions for weak magnetic fields (B < 1 G) and periodic ion-acoustic instabilities can develop for strong magnetic fields (B > 10 G). We have found the instability growth rates and estimated the temporal and spatial scales of development and decay of the periodic instability. We show that the instabilities under consideration can play a prominent role in the energy balance of the corona and may be considered as a large-scale energy source of the wave coronal heating mechanism.     

Nonlinear excitation of small-scale Alfvén waves by fast waves and plasma heating in the solar atmosphere

We study a nonlinear mechanism for the excitation of kinetic Alfvén waves (KAWs) by fast magneto-acoustic waves (FWs) in the solar atmosphere. Our focus is on the excitation of KAWs that have very small wavelengths in the direction perpendicular to the background magnetic field. Because of their small perpendicular length scales, these waves are very efficient in the energy exchange with plasmas and other waves. We show that the nonlinear coupling of the energy of the finite-amplitude FWs to the small-scale KAWs can be much faster than other dissipation mechanisms for fast wave, such as electron viscous damping, Landau damping, and modulational instability. The nonlinear damping of the FWs due to decay FW = KAW + KAW places a limit on the amplitude of the magnetic field in the fast waves in the solar corona and solar-wind at the level B/B ₀∼10⁻². In turn, the nonlinearly excited small-scale KAWs undergo strong dissipation due to resistive or Landau damping and can provide coronal and solar-wind heating. The transient coronal heating observed by Yohkoh and SOHO may be produced by the kinetic Alfvén waves that are excited by parametric decay of fast waves propagating from the reconnection sites.     

The Creation of Outflowing Plasma in the Corona at Emerging Flux Regions: Comparing Observations and Simulations

In this paper we analyse the flux emergence that occurred in the following polarity area of an active region on 1 – 2 December 2006. Observations have revealed the existence of fast outflows at the edge of the emerging flux region. We have performed 3-D numerical simulations to study the mechanisms responsible for these flows. The results indicate that these outflows are reconnection jets or pressure-driven outflows, depending on the relative orientation of the magnetic fields in contact (i.e. the emerging flux and the active region’s field which is favourable for reconnection on the west side and nearly parallel with the pre-existing field on the east side of the emerging flux). In the observations, the flows are larger on the west side until late in the flux emergence, when the reverse is true. The simulations show that the flows are faster on the west side, but do not show the east flows increasing with time. There is an asymmetry in the expansion of the emerging flux region, which is also seen in the observations. The west side of the emerging flux region expands faster into the corona than the other side. In the simulations, efficient magnetic reconnection occurs on the west side, with new loops being created containing strong downflows that are clearly seen in the observations. On the other side, the simulations show strong compression as the dominant mechanism for the generation of flows. There is evidence of these flows in the observations, but the flows are stronger than the simulations predict at the later stages. There could be additional small-angle reconnection that adds to the flows from the compression, as well as reconnection occurring in larger loops that lie across the whole active region.     

Global Simulation of Magnetospheric Space Weather Effects of the Bastille Day Storm

We present results from a global simulation of the interaction of the solar wind with Earth's magnetosphere, ionosphere, and thermosphere for the Bastille Day geomagnetic storm and compare the results with data. We find that during this event the magnetosphere becomes extremely compressed and eroded, causing 3 geosynchronous GOES satellites to enter the magnetosheath for an extended time period. At its extreme, the magnetopause moves at local noon as close as 4.9 R _E to Earth which is interpreted as the consequence of the combined action of enhanced dynamic pressure and strong dayside reconnection due to the strong southward interplanetary magnetic field component B _z, which at one time reaches a value of −60 nT. The lobes bulge sunward and shield the dayside reconnection region, thereby limiting the reconnection rate and thus the cross polar cap potential. Modeled ground magnetic perturbations are compared with data from 37 sub-auroral, auroral, and polar cap magnetometer stations. While the model can not yet predict the perturbations and fluctuations at individual ground stations, its predictions of the fluctuation spectrum in the 0–3 mHz range for the sub-auroral and high-latitude regions are remarkably good. However, at auroral latitudes (63° to 70° magnetic latitude) the predicted fluctuations are slightly too high. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1014228230714     

<Emphasis Type="Italic">In Situ</Emphasis> Signatures of Interchange Reconnection between Magnetic Clouds and Open Magnetic Fields: A Mechanism for the Erosion of Polar Coronal Holes?

We outline a method to determine the direction of solar open flux transport that results from the opening of magnetic clouds (MCs) by interchange reconnection at the Sun based solely on in-situ observations. This method uses established findings about i) the locations and magnetic polarities of emerging MC footpoints, ii) the hemispheric dependence of the helicity of MCs, and iii) the occurrence of interchange reconnection at the Sun being signaled by uni-directional suprathermal electrons inside MCs. Combining those observational facts in a statistical analysis of MCs during solar cycle 23 (period 1995 – 2007), we show that the time of disappearance of the northern polar coronal hole (1998 – 1999), permeated by an outward-pointing magnetic field, is associated with a peak in the number of MCs originating from the northern hemisphere and connected to the Sun by outward-pointing magnetic field lines. A similar peak is observed in the number of MCs originating from the southern hemisphere and connected to the Sun by inward-pointing magnetic field lines. This pattern is interpreted as the result of interchange reconnection occurring between MCs and the open field lines of nearby polar coronal holes. This reconnection process closes down polar coronal hole open field lines and transports these open field lines equatorward, thus contributing to the global coronal magnetic field reversal process. These results will be further constrainable with the rising phase of solar cycle 24.     

Magnetic reconnection in plasmas

A review of the present status of the theory of magnetic reconnection is given. In strongly collisional plasmas reconnection proceeds via resistive current sheets, i.e. quasi-stationary macroscopic Sweet-Parker sheets at intermediate values of the magnetic Reynolds numberR _m , or mirco-current sheets in MHD turbulence, which develops at highR _m . In hot, dilute plasmas the reconnection dynamics is dominated by nondissipative effects, mainly the Hall term and electron inertia. Reconnection rates are found to depend only on the ion mass, being independent of the electron inertia and the residual dissipation coefficients. Small-scale whistler turbulence is readily excited giving rise to an anomalous electron viscosity. Hence reconnection may be much more rapid than predicted by conventional resistive theory.     

Galaxy rotation curves: the effect of $\vec{j} \times\vec{B}$ force

Using the Galaxy as an example, we study the effect of [(j)\vec] ×[(B)\vec]\vec{j} \times \vec{B} force on the rotational curves of gas and plasma in galaxies. Acceptable model for the galactic magnetic field and plausible physical parameters are used to fit the flat rotational curve for gas and plasma based on the observed baryonic (visible) matter distribution and [(j)\vec] ×[(B)\vec]\vec{j} \times\vec{B} force term in the static MHD equation of motion. We also study the effects of varied strength of the magnetic field, its pitch angle and length scale on the rotational curves. We show that [(j)\vec] ×[(B)\vec]\vec{j} \times\vec{B} force does not play an important role on the plasma dynamics in the intermediate range of distances 6–12 kpc from the centre, whilst the effect is sizable for larger r (r≥15 kpc), where it is the most crucial.