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.     

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.     

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.     

Modelling of Solar Wind, CME Initiation and CME Propagation

Simulations of coronal mass ejections (CMEs) evolving in the interplanetary (IP) space from the Sun up to 1 AU are performed in the framework of ideal magnetohydrodynamics (MHD) by the means of a finite-volume, explicit solver. The aim is to quantify the effect of the background solar wind and of the CME initiation parameters, such as the initial magnetic polarity, on the evolution and on the geo-effectiveness of CMEs. First, three different solar wind models are reconstructed using the same numerical grid and the same numerical scheme. Then, different CME initiation models are considered: Magnetic foot point shearing and magnetic flux emergence. For the fast CME evolution studies, a very simple CME model is considered: A high-density and high-pressure magnetized plasma blob is superposed on a background steady state solar wind model with an initial velocity and launch direction. The simulations show that the initial magnetic polarity substantially affects the IP evolution of the CMEs influencing the propagation velocity, the shape, the trajectory (and thus, the geo-effectiveness).     

The Solar Wind - Inner Heliosphere

The solar wind in the inner heliosphere, inside ~ 5 AU, has been almost fully characterized by the addition of the high heliographic latitude Ulysses mission to the many low latitude inner heliosphere missions that preceded it. The two major omissions are the high latitude solar wind at solar maximum, which will be measured during the second Ulysses polar passages, and the solar wind near the Sun, which could be analyzed by a Solar Probe mission. Here, existing knowledge of the global solar wind in the inner heliosphere is summarized in the context of the new results from Ulysses.     

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.     

MKDVB and CKB Shock Waves

This review presents results from MKDVB and CKB systems to illustrate the properties of magnetohydrodynamics intermediate shock waves. The intermediate shocks can have a family of shock structures; the shock structure and the evolution are related. Furthermore, there exists a class of time-dependent intermediate shocks. This revised version was published online in August 2006 with corrections to the Cover Date.