A stellar wind passing through the reverse shock is deflected into the astrospheric tail and leaves the stellar system either as a sub-Alfvenic or as a super-Alfvenic tail flow. An example is our own heliosphere and its heliotail. We present an analytical method of calculating stationary, incompressible, and field-aligned plasma flows in the astrotail of a star. We present a recipe for constructing an astrosphere with the help of only a few parameters, like the inner Alfven Mach number and the outer Alfven Mach number, the magnetic field strength within and outside the stellar wind cavity, and the distribution of singular points of the magnetic field within these flows. Within the framework of a one-fluid approximation, it is possible to obtain solutions of the MHD equations for stationary flows from corresponding static MHD equilibria, by using noncanonical mappings of the canonical variables. The canonical variables are the Euler potentials of the magnetic field of magnetohydrostatic equilibria. Thus we start from static equilibria determined by the distribution of magnetic neutral points, and assume that the Alfven Mach number for the corresponding stationary equilibria is finite. The topological structure determines the geometrical structure of the interstellar gas - stellar wind interface. Additional boundary conditions like the outer magnetic field and the jump of the magnetic field across the astropause allow determination of the noncanonical transformations. This delivers the strength of the magnetic field at every point in the astrotail region beyond the reverse shock. The mathematical technique for describing such a scenario is applied to astrospheres in general, but is also relevant for the heliosphere. It shows the restrictions of the outer and the inner magnetic field strength in comparison with the corresponding Alfven Mach numbers in the case of subalfvenic flows.

We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized accretion disk equilibrium and linear wave modes about this equilibrium. We show the existence of novel MHD instabilities in these two--dimensional equilibria which do not occur in an accretion disk in the cylindrical limit. The disk equilibria are numerically computed by the FINESSE code. The stability of accretion disks is investigated analytically as well as numerically. We use the PHOENIX code to compute all the waves and instabilities accessible to the computed disk equilibrium. We concentrate on strongly magnetized disks and sub--Keplerian rotation in a large part of the disk. These disk equilibria show that the thermal pressure of the disk can only decrease outwards if there is a strong gravitational potential. Our theoretical stability analysis shows that convective continuum instabilities can only appear if the density contours coincide with the poloidal magnetic flux contours. Our numerical results confirm and complement this theoretical analysis. Furthermore, these results show that the influence of gravity can either be stabilizing or destabilizing on this new kind of MHD instability. In the likely case of a non--constant density, the height of the disk should exceed a threshold before this type of instability can play a role. This localised MHD instability provides an ideal, linear route to MHD turbulence in strongly magnetized accretion disk tori.

We present analytical and numerical studies of magnetorotational instabilities occuring in magnetized accretion disks. In these studies we make use of the linearised compressible MHD equations. These calculations are performed for general radially stratified disks in the cylindrical limit. In particular, we investigate the influence of nonvanishing toroidal magnetic field component on the growth rate and oscillation frequency of magnetorotational instabilities in Keplerian disks. We find the persistence of these instabilities in accretion disks close to equipartition. Our calculations show that these eigenmodes become overstable (complex eigenvalue), due to the presence of a toroidal magnetic field component, while their growth rate reduces slightly. Furthermore, we demonstrate the presence of magneto-rotational overstabilities in weakly magnetized sub-Keplerian rotating disks. We show that the growth rate scales with the rotation frequency of the disk. These eigenmodes also have a nonzero oscillation frequency, due to the presence of the dominant toroidal magnetic field component. The overstable character of the MRI increases as the rotation frequency of the disk decreases.

We present a numerical study of an idealized magnetohydrodynamic (MHD) configuration consisting of a planar wake flow embedded into a three-dimensional (3D) sheared magnetic field. Our simulations investigate the possibility for in-situ development of large-scale compressive disturbances at cospatial current sheet -- velocity shear regions in the heliosphere. Using a linear MHD solver, we first systematically chart the destabilized wavenumbers, corresponding growth rates, and physical parameter ranges for dominant 3D sinuous-type instabilities in an equilibrium wake--current sheet system. Wakes bounded by sufficiently supersonic (Mach number $M_s > 2.6$) flow streams are found to support dominant fully 3D sinuous instabilities when the plasma beta is of order unity. Fully nonlinear, compressible 2.5D and 3D MHD simulations show the self-consistent formation of shock fronts of fast magnetosonic type. They carry density perturbations far away from the wake's center. Shock formation conditions are identified in sonic and Alfv\'enic Mach number parameter space. Depending on the wake velocity contrast and magnetic field magnitude, as well as on the initial perturbation, the emerging shock patterns can be plane-parallel as well as fully three-dimensionally structured. Similar large-scale transients could therefore originate at distances far above coronal helmet streamers or at the location of the ecliptic current sheet.

Obtainable computational efficiency is evaluated when using an Adaptive Mesh Refinement (AMR) strategy in time accurate simulations governed by sets of conservation laws. For a variety of 1D, 2D, and 3D hydro- and magnetohydrodynamic simulations, AMR is used in combination with several shock-capturing, conservative discretization schemes. Solution accuracy and execution times are compared with static grid simulations at the corresponding high resolution and time spent on AMR overhead is reported. Our examples reach corresponding efficiencies of 5 to 20 in multidimensional calculations and only 1.5 -- 8 % overhead is observed. For AMR calculations of multi-dimensional magnetohydrodynamic problems, several strategies for controlling the $\nabla \cdot \BB =0$ constraint are examined. Three source term approaches suitable for cell-centered $\BB$ representations are shown to be effective. For 2D and 3D calculations where a transition to a more globally turbulent state takes place, it is advocated to use an approximate Riemann solver based discretization at the highest allowed level(s), in combination with the robust Total Variation Diminishing Lax-Friedrichs method on the coarser levels. This level-dependent use of the spatial discretization acts as a computationally efficient, hybrid scheme.

A complete analytical and numerical treatment of all magnetohydrodynamic waves and instabilities for radially stratified, magnetized accretion disks is presented. The instabilities are a possible source of anomalous transport. While recovering results on known hydrodynamicand both weak- and strong-field magnetohydrodynamic perturbations, the full magnetohydrodynamic spectra for a realistic accretion disk model demonstrates a much richer variety of instabilities accessible to the plasma than previously realized. We show that both weakly and strongly magnetized accretion disks are prone to strong non-axisymmetric instabilities.The ability to characterize all waves arising in accretion disks holds great promise for magnetohydrodynamic spectroscopic analysis.

Axisymmetric stellar wind solutions are presented, obtained by numerically solving the ideal magnetohydrodynamic (MHD) equations. Stationary solutions are critically analysed using the knowledge of the flux functions. These flux functions enter in the general variational principle governing all axisymmetric stationary ideal MHD equilibria. The magnetized wind solutions for (differentially) rotating stars contain both a `wind' and a `dead' zone. We illustrate the influence of the magnetic field topology on the wind acceleration pattern, by varying the coronal field strength and the extent of the dead zone. This is evident from the resulting variations in the location and appearance of the critical curves where the wind speed equals the slow, Alfven, and fast speed. Larger dead zones cause effective, fairly isotropic acceleration to super-Alfvenic velocities as the polar, open field lines are forced to fan out rapidly with radial distance. A higher field strength moves the Alfven transition outwards. In the ecliptic, the wind outflow is clearly modulated by the extent of the dead zone. The combined effect of a fast stellar rotation and an equatorial `dead' zone in a bipolar field configuration can lead to efficient thermo-centrifugal equatorial winds. Such winds show both a strong poleward collimation and some equatorward streamline bending due to significant toroidal field pressure at mid-latitudes. We discuss how coronal mass ejections are then simulated on top of the transonic outflows.

We discuss steady-state transonic outflows obtained by direct numerical solution of the hydrodynamic and magnetohydrodynamic equations. We make use of the Versatile Advection Code, a software package for solving systems of (hyperbolic) partial differential equations. We proceed stepwise from a spherically symmetric, isothermal, unmagnetized, non-rotating Parker wind to arrive at axisymmetric, polytropic, magnetized, rotating models. These represent 2D generalisations of the analytical 1D Weber-Davis wind solution, which we obtain in the process. Axisymmetric wind solutions containing both a `wind' and a `dead' zone are presented. Since we are solving for steady-state solutions, we efficiently exploit fully implicit time stepping. The method allows us to model thermally and/or magneto-centrifugally driven stellar outflows. We particularly emphasize the boundary conditions imposed at the stellar surface. For these axisymmetric, steady-state solutions, we can use the knowledge of the flux functions to verify the physical correctness of the numerical solutions.

We investigate the Kelvin-Helmholtz instability occuring at the interface of a shear flow configuration in 2D compressible magnetohydrodynamics (MHD). The linear growth and the subsequent non-linear saturation of the instability are studied numerically. We consider an initial magnetic field aligned with the shear flow, and analyze the differences between cases where the initial field is unidirectional everywhere (uniform case), and where the field changes sign at the interface (reversed case). We recover and extend known results for pure hydrodynamic and MHD cases with a discussion of the dependence of the non-linear saturation on the wavenumber, the sound Mach number, and the Alfvenic Mach number for the MHD case. A reversed field acts to destabilize the linear phase of the Kelvin-Helmholtz instability compared to the pure hydrodynamic case, while a uniform field suppresses its growth. In resistive MHD, reconnection events almost instantly accelerate the buildup of a global plasma circulation. They play an important role throughout the further non-linear evolution as well, since the initial current sheet gets amplified by the vortex flow and can become unstable to tearing instabilities forming magnetic islands. As a result, the saturation behaviour and the overall evolution of the density and the magnetic field is markedly different for the uniform versus the reversed field case.