Model of a fluxtube with a twisted magnetic field in the stratified solar atmosphere

We build a single vertical straight magnetic fluxtube spanning the solar photosphere and the transition region which does not expand with height. We assume that the fluxtube containing twisted magnetic fields is in magnetohydrostatic equilibrium within a realistic stratified atmosphere subject to solar gravity. Incorporating specific forms of current density and gas pressure in the Grad–Shafranov equation, we solve the magnetic flux function, and find it to be separable with a Coulomb wave function in radial direction while the vertical part of the solution decreases exponentially. We employ improved fluxtube boundary conditions and take a realistic ambient external pressure for the photosphere to transition region, to derive a family of solutions for reasonable values of the fluxtube radius and magnetic field strength at the base of the axis that are the free parameters in our model. We find that our model estimates are consistent with the magnetic field strength and the radii of Magnetic bright points (MBPs) as estimated from observations. We also derive thermodynamic quantities inside the fluxtube.     

磁流体中间激波的相互作用

本文讨论磁流体中间激波的相互作用规律．主要结论：中间激波汇合的产物为后向快简单波、后向慢简单波或慢激波、接触间断、前向慢激波和前向快激波，其中后向波成份和接触间断很弱．当左（右）激波较强时，中国激波碰撞产物为后（前）向快激波、后（前）向慢简单波或慢激波、负（正）切向间断、前（后）向慢简单波和前（后）向快激波．     

Determining CME parameters by fitting heliospheric observations: Numerical investigation of the accuracy of the methods

Transients in the heliosphere, including coronal mass ejections (CMEs) and corotating interaction regions can be imaged to large heliocentric distances by heliospheric imagers (HIs), such as the HIs onboard STEREO and SMEI onboard Coriolis. These observations can be analyzed using different techniques to derive the CME speed and direction. In this paper, we use a three-dimensional (3-D) magneto-hydrodynamic (MHD) numerical simulation to investigate one of these methods, the fitting method of  and . Because we use a 3-D simulation, we can determine with great accuracy the CME initial speed, its speed at 1 AU and its average transit speed as well as its size and direction of propagation. We are able to compare the results of the fitting method with the values from the simulation for different viewing angles between the CME direction of propagation and the Sun-spacecraft line. We focus on one simulation of a wide (120–140°) CME, whose initial speed is about 800 km s⁻¹. For this case, we find that the best-fit speed is in good agreement with the speed of the CME at 1 AU, and this, independently of the viewing angle. The fitted direction of propagation is not in good agreement with the viewing angle in the simulation, although smaller viewing angles result in smaller fitted directions. This is due to the extremely wide nature of the ejection. A new fitting method, proposed to take into account the CME width, results in better agreement between fitted and actual directions for directions close to the Sun–Earth line. For other directions, it gives results comparable to the fitting method of Sheeley et al. (1999). The CME deceleration has only a small effect on the fitted direction, resulting in fitted values about 1–4° higher than the actual values.     

MHD Spectroscopy of Transonic Flows

In previous publications (Keppens et al.: 2002, Astrophys. J. 569, L121; Goedbloed et al.: 2004a, Phys. Plasmas 11, 28), we have demonstrated that stationary rotation of magnetized plasma about a compact central object permits an enormous number of different MHD instabilities, with the well-known magneto-rotational instability (Velikhov, E. P.: 1959, Soviet Phys.–JETP Lett. 36, 995; Chandrasekhar, S.: 1960, Proc. Natl. Acad. Sci. U.S.A. 46, 253; Balbus, S. A. and Hawley, J. F.: 1991, Astrophys. J. 376, 214) as just one of them. We here concentrate on the new instabilities found that are driven by transonic transitions of the poloidal flow. A particularly promising class of instabilities, from the point of view of MHD turbulence in accretion disks, is the class of trans-slow Alfv’en continuum modes, that occur when the poloidal flow exceeds a critical value of the slow magnetosonic speed. When this happens, virtually every magnetic/flow surface of the disk becomes unstable with respect to highly localized modes of the continuous spectrum. The mode structures rotate, in turn, about the rotating disk. These structures lock and become explosively unstable when the mass of the central object is increased beyond a certain critical value. Their growth rates then become huge, of the order of the Alfv’en transit time. These instabilities appear to have all requisite properties to facilitate accretion flows across magnetic surfaces and jet formation.     

Synthesis of CME-Associated Moreton and EIT Wave Features from MHD Simulations

Soft X-ray (SXR) waves, EIT waves, and Hα Moreton waves are all associated with coronal mass ejections (CMEs). The knowledge of the characteristics about these waves is crucial for the understanding of CMEs, and hence for the space weather researches. MHD numerical simulation is performed, with the consideration of the quiet Sun atmosphere, to investigate the CME/flare processes. On the basis of the numerical results, SXR, EUV, and Hα images of the eruption are synthesized, where SXR waves, EIT waves, and Hα Moreton waves are identified. It confirms that the EIT waves, which border the expanding dimmming region, are produced by the successive opening (or stretching) of the closed magnetic field lines. Hα Moreton waves are found to propagate outward synchronously with the SXR waves, lagging behind the latter spatially by ～27 Mm in the simulated scenario. However, the EIT wave velocity is only a third of the Moreton wave velocity. The synthesized results also suggest that Hα± 0.45Å would be the best off-band for the detection of Hα Moreton waves.     

Solar filament eruptions and their physical role in triggering coronal mass ejections

Solar filament eruptions play a crucial role in triggering coronal mass ejections (CMEs). More than 80% of eruptions lead to a CME. This correlation has been studied extensively during the past solar cycles and the last long solar minimum. The statistics made on events occurring during the rising phase of the new solar cycle 24 is in agreement with this finding. Both filaments and CMEs have been related to twisted magnetic fields. Therefore, nearly all the MHD CME models include a twisted flux tube, called a flux rope. Either the flux rope is present long before the eruption, or it is built up by reconnection of a sheared arcade from the beginning of the eruption.     

Why should the latitude of the observer be considered when modeling gradual proton events? An insight using the concept of cobpoint

The shape of flux profiles of gradual solar energetic particle (SEP) events depends on several not well-understood factors, such as the strength of the associated shock, the relative position of the observer in space with respect to the traveling shock, the existence of a background seed particle population, the interplanetary conditions for particle transport, as well as the particle energy. Here, we focus on two of these factors: the influence of the shock strength and the relative position of the observer. We performed a 3D simulation of the propagation of a coronal/interplanetary CME-driven shock in the framework of ideal MHD modeling. We analyze the passage of this shock by nine spacecraft located at ∼0.4 AU (Mercury’s orbit) and at different longitudes and latitudes. We study the evolution of the plasma conditions in the shock front region magnetically connected to each spacecraft, that is the region of the shock front scanned by the “cobpoint” (Heras et al., 1995), as the shock propagates away from the Sun. Particularly, we discuss the influence of the latitude of the observer on the injection rate of shock-accelerated particles and, hence, on the resulting proton flux profiles to be detected by each spacecraft.     

Three-dimensional path planning for unmanned aerial vehicle based on interfered fluid dynamical system

This paper proposes a method for planning the three-dimensional path for low-flying unmanned aerial vehicle(UAV) in complex terrain based on interfered fluid dynamical system(IFDS) and the theory of obstacle avoidance by the flowing stream. With no requirement of solutions to fluid equations under complex boundary conditions, the proposed method is suitable for situations with complex terrain and different shapes of obstacles. Firstly, by transforming the mountains, radar and anti-aircraft fire in complex terrain into cylindrical, conical, spherical, parallelepiped obstacles and their combinations, the 3D low-flying path planning problem is turned into solving streamlines for obstacle avoidance by fluid flow. Secondly, on the basis of a unified mathematical expression of typical obstacle shapes including sphere, cylinder, cone and parallelepiped, the modulation matrix for interfered fluid dynamical system is constructed and 3D streamlines around a single obstacle are obtained. Solutions to streamlines with multiple obstacles are then derived using weighted average of the velocity field. Thirdly, extra control force method and virtual obstacle method are proposed to deal with the stagnation point and the case of obstacles’ overlapping respectively. Finally, taking path length and flight height as sub-goals, genetic algorithm(GA) is used to obtain optimal 3D path under the maneuverability constraints of the UAV. Simulation results show that the environmental modeling is simple and the path is smooth and suitable for UAV. Theoretical proof is also presented to show that the proposed method has no effect on the characteristics of fluid avoiding obstacles.     

Hall MHD reconnection with an initial guide field By0

An uniform out-of-plane magnetic field component B_y0 is added to the equilibrium Harris sheet with plasma β = 0.5 and L_c = 0.5d_i (where L_c is the half-width of the equilibrium current layer and d_i is the ion inertial length). Driven by the continuous boundary inflows, the magnetic reconnections with the guide field B_y0/B₀ ranging from 0 to 4.0 are investigated using a 2.5D Hall magnetohydro-dynamic (MHD) code developed from a multi-step implicit scheme. The features of the reconnection field are substantially altered in the presence of the guide field. The openness of the magnetic separatrix angle is slightly reduced and the anti-symmetric quadrupolar structure of B_y field and the symmetric distribution of plasma pressure P are replaced by an asymmetric B_y four-wing structure and an asymmetric P plot as a non-zero B_y0 is added. The decoupling of electrons and ions also occurs near the X line in the case with a finite B_y0, but the effect of initial B_y0 on the electron flow is greater than that on the ion flow. The reconnection rates at the X-line drops from 0.151 to 0.06, namely, ∂A/∂t is reduced by a factor of 2.5 as B_y0/B₀ increases from 0 to 4.0. The reduction of reconnection rate might be related to the reducing openness of reconnection layer with the increasing B_y0.     

Estimation of the bias of the Minimum Variance technique in the determination of magnetic clouds global quantities and orientation

Magnetic clouds (MCs) are highly magnetized plasma structures that have a low proton temperature and a magnetic field vector that rotates when seen by a heliospheric observer. More than 25 years of observations of magnetic and plasma properties of MCs at 1 AU have provided significant knowledge of their magnetic structure. However, because in situ observations only give information along the trajectory of the spacecraft, their real 3D magnetic configuration remains still partially unknown. We generate a set of synthetic clouds, exploring the space of parameters that represents the possible orientations and minimum distances of the satellite trajectory to the cloud axis, p. The synthetic clouds have a local cylindrical symmetry and a linear force-free magnetic configuration. From the analysis of synthetic clouds, we quantify the errors introduced in the determination of the orientation/size (and, consequently, of the global magnetohydrodynamic quantities) by the Minimum Variance method when p is not zero.