Viscous computations of resonant absorption of MHD waves in flux tubes by fem

A numerical code is presented for computing the stationary state of resonant absorption of MHD waves in cylindrical flux tubes in linear, compressible, and viscous MHD. The full viscosity stress tensor is included in the code with the five viscosity coefficients as given by Braginskii (1965). Also non-zero plasma pressure effects are taken into account, and the finite elements discretization with the Galerkin method has been used. The implementation of the stress tensor and the numerical accuracy of the tensorial viscous MHD code are scrutinized in test case. The test case involves the absorption of waves in cylindrical flux tubes considered by Lou (1990) and Goossens and Poedts (1992) in the context of absorption of acoustic oscillations. The results for the absorption rates obtained with the tensorial viscous code agree completely with the results obtained by Lou in a scalar viscous MHD and by Goossens and Poedts in resistive MHD. This verifies not only the complicated tensor viscous code but again proves that the absorption rate is independent of the actual dissipation mechanism.     

Steady and Time-Dependent MHD Modelling of Jets

A brief review is given of some results of our work on the construction of (I) steady and (II) time-dependent MHD models for nonrelativistic and relativistic astrophysical outflows and jets, analytically and numerically. The only available exact solutions for MHD outflows are those in separable coordinates, i.e., with the symmetry of radial or meridional self-similarity. Physically accepted solutions pass from the fast magnetosonic separatrix surface in order to satisfy MHD causality. An energetic criterion is outlined for selecting radially expanding winds from cylindrically expanding jets. Numerical simulations of magnetic self-collimation verify the conclusions of analytical steady solutions. We also propose a two-component model consisting of a wind outflow from a central object and a faster rotating outflow launched from a surrounding accretion disk which plays the role of the flow collimator. We also discuss the problem of shock formation during the magnetic collimation of wind-type outflows into jets.     

Simulating Magnetized Jets

A suitable model for the macroscopic behavior of accretion disk-jet systems is provided by the equations of MagnetoHydroDynamics (MHD). These equations allow us to perform scale-encompassing numerical simulations of multidimensional nonlinear magnetized plasma flows. For that purpose, we continue the development and exploitation of the Versatile Advection Code (VAC) along with its recent extension which employs dynamically controlled grid adaptation. In the adaptive mesh refinement AMRVAC code, modules for simulating any-dimensional special relativistic hydro- and magnetohydrodynamic problems are currently operational. Here, we review recent 3D MHD simulations of fundamental plasma instabilities, relevant when dealing with cospatial shear flow and twisted magnetic fields. Such magnetized jet flows can be susceptible to a wide variety of hydro (e.g. Kelvin-Helmholtz) or magnetohydrodynamic (e.g. current driven kink) instabilities. Recent MHD computations of 3D jet flows have revealed how such mutually interacting instabilities can in fact aid in maintaining jet coherency. Another breakthrough from computational magnetofluid modeling is the demonstration of continuous, collimated, transmagnetosonic jet launching from magnetized accretion disks. Summarizing, MHD simulations are rapidly gaining realism and significantly advance our understanding of nonlinear astrophysical magnetofluid dynamics.     

Magnetic-Tower Jet Solution for Launching Astrophysical Jets

In spite of the large number of global three-dimensional (3-D) magnetohydrodynamic (MHD) simulations of accretion disks and astrophysical jets, which have been developed since 2000, the launching mechanisms of jets is somewhat controversial. Previous studies of jets have concentrated on the effect of the large-scale magnetic fields permeating accretion disks. However, the existence of such global magnetic fields is not evident in various astrophysical objects, and their origin is not well understood. Thus, we study the effect of small-scale magnetic fields confined within the accretion disk. We review our recent findings on the formation of jets in dynamo-active accretion disks by using 3-D MHD simulations. In our simulations, we found the emergence of accumulated azimuthal magnetic fields from the inner region of the disk (the so-called magnetic tower) and also the formation of a jet accelerated by the magnetic pressure of the tower. Our results indicate that the magnetic tower jet is one of the most promising mechanisms for launching jets from the magnetized accretion disk in various astrophysical objects. We will discuss the formation of cosmic jets in the context of the magnetic tower model.     

Formation of astrophysical jets by a contracting magnetic accretion disk

In the present paper, we discuss an MHD model for the formation of astrophysical jets, in which the directed flows are ejected along the rotation axis of an accretion disk formed from a cloud having a large scale magnetic field parallel to the angular momentum axis of the disk. The acceleration of jets is due to thej×B force in the relaxing magnetic twist which is produced by the rotation of the disk. The characteristic features of the jets, predicted by our mechanism and hopefully to be proven by observations, are the helical velocity and the hollow cylindrical shape of the jet, with a diameter of roughly the size of the region from which the acceretion disk collected its mass. Justification for the assumption of the perpendicular orientation of the disk, or the parallelism of the jets, to the external magnetic field may be provided by the fact that the component of rotation whose axis is perpendicular to the field may have been damped in the earlier phase of the cloud contraction.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 Septemper–6 October, 1984.     

On the MHD Acceleration of Astrophysical Jets

We present a 2.5D magnetohydrodynamic (MHD) simulation of the acceleration of a collimated jet from a magnetized accretion disk. We employ a MHD Adaptive Mesh Refinement (AMR) code (FLASH—University of Chicago). Thanks to this tool we can follow the evolution of the system for many dynamical timescales with a high-spatial resolution. Assuming an initial condition in which a Keplerian disk, thus with no accretion motions, is threaded by a uniform poloidal magnetic field, we show how both the accretion flow and the acceleration of the outflow occur, and we present in detail which are the forces responsible for the jet launching and collimation. Our simulation also shows how the collimating forces due to the self-generated toroidal magnetic field can produce some peculiar knotty features.     

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.     

Local potential analysis of MHD instability

In many astrophysical problems, the study of the stability of an atmosphere in the presence of a magnetic field is of importance. In most cases the MHD instabilities of atmospheres are studied by energy principle of Bernsteinet al. (1958). In this paper, a general method for studying the stability of a system subject to MHD equations of conditions has been proposed. This is based on the local potential concept put forward by Glansdorff and Prigogine (1964). The scheme for securing stability criteria has been demonstrated in two particular cases.     

A multidimensional magnetohydrodynamics code based on semidiscrete central WENO and PPM schemes

We present a multidimensional magnetohydrodynamics code based on semidiscrete central Godunov-type schemes and high order weighted essentially non-oscillatory data reconstruction and incorporate elements of the piecewise parabolic method into the reconstruction schemes. This algorithm avoids spectral decomposition into characteristic waves and is easily implemented. We have implemented this algorithm and carried out one and two dimensional tests, which verifies that this code can achieve the accuracy and efficiency to solve the astrophysical problems.     

Turbulence in Accretion Disks

An accretion disk is an inevitable part of the star forming process. Recent years have witnessed dramatic progress in our understanding of how turbulence arises and transports angular momentum in astrophysical accretion disks. The key conceptual point is that the combination of a subthermal magnetic field and outwardly decreasing differential rotation is subject to the magnetorotational instability. This rapidly generates magnetohydrodynamical (MHD) turbulence, leading to greatly enhanced angular momentum transport. Purely hydrodynamic disks, on the other hand, are stable. Disks that are too cool to couple effectively to the magnetic field will not be turbulent. Fully global three dimensional MHD simulations are now beginning to probe the properties of accretion disks from first principles.     

A Rapid X-ray Flare from Markarian 501

The tree code for the approximate evaluation of gravitational forces is extended and substantially accelerated by including mutual cell-cell interactions. These are computed by a Taylor series in Cartesian coordinates and in a completely symmetric fashion, such that Newton's third law is satisfied by construction and that therefore momentum is exactly conserved. The computational effort is further reduced by exploiting the mutual symmetry of the interactions. For typical astrophysical problems with N=105 and at the same level of accuracy, the new code is about 4 times faster than the tree code. For large N, the computational costs are found to scale almost linearly with N, which can also be supported by a theoretical argument, and the advantage over the tree code increases with ever larger N.     

Steady MHD Solutions for Collimated Winds

A convenient approach to model MHD steady axisymmetric outflows is the so-called self-similar technique wherein the physical variables are factorized and a scaling law is assumed along one of the coordinates. This scaling depends on the astrophysical process under investigation. In this note we summarize all possible self-similar MHD outflow solutions; furthermore, we briefly discuss the main properties of a class of solutions which are self-similar in the meridional direction and allow to analyse in simple terms the dynamical properties of an outflow close to its rotational axis. Special attention is focused on the asymptotic structure of collimated winds. It will be shown that different regimes are possible for jets, in particular they can be either thermally or magnetically confined, depending on the physical conditions of the flow. This analysis is complementary with the well known radial self-similar models which are invoked to study winds from accretion disks. This revised version was published online in July 2006 with corrections to the Cover Date.     

Contraction and Fragmentation of Magnetized Rotating Clouds and Formation of Binary Systems

Using three-dimensional (3D) magnetohydrodynamical (MHD) nested-grid simulations, the fragmentation of a rotating magnetized molecular cloud core is studied. An isothermal rotating magnetized cylindrical cloud in hydrostatic balance is considered. We studied non-axisymmetric evolution of the cloud. It is found that non-axisymmetry hardly evolves in the early phase, but it begins to grow after the gas contracts and forms a thin disk. The disk formation and thus growth of non-axisymmetric perturbation are strongly promoted by rotation and magnetic field strength. We found two types of fragmentations: fragmentation from a ring and that from a bar. These two types of fragmentations occur in thin adiabatic cores with the thickness being smaller than 1/4 of the radial size. For the fragments to survive, they should be formed in a heavily elongated barred core or a flat round disk. In the models showing fragmentation, outflows from respective fragments are found as well as that driven by the rotating bar or the disk.     

How Can Jets Survive MHD Instabilities?

We present the main findings of two recent studies using high-resolution MHD simulations of supersonic magnetized shear flow layers. First, a strong large-scale coalescence effect partially countered by small-scale reconnection events is shown to dominate the dynamics in a two-dimensional layer subject to Kelvin-Helmholtz (KH) instabilities. Second, an interaction mechanism between two different types of instabilities (KH and current-driven modes) is shown to occur in a cylindrical jet configuration embedded in an helical magnetic field. Finally, we discuss the implications of these results for astrophysical jets survival.     

Theoretical and Experimental Studies of Radiative Shocks

This paper deals with the radiative shock from both theoretical and numerical points of view. It is based on the whole experimental results obtained at Laboratoire d'Utilisation des Lasers Intenses (LULI, école Polytechnique). Radiative shocks are high-Mach number shocks with a strong coupling between radiation and hydrodynamics which leads to a structure governed by a radiative precursor. These shocks are involved in various astrophysical systems: stellar accretion shocks, pulsating stars, interaction between supernovae and the interstellar medium. In laboratory, these radiative shocks are generated using high power lasers. New diagnostics have been implemented to study the geometrical shape of the shock and the front shock density. Data were obtained varying initial conditions for different laser intensities and temperature. The modeling of these phenomena is mainly performed through numerical simulations (1D and 2D) and analytical studies. We exhibit results obtained from several radiative hydrodynamics codes. As a result, it is possible to discuss about the influence of the geometry and physical parameters introduced in the 1D and 2D models.     

A new computational fluid dynamics code I: Fyris Alpha

A new hydrodynamics code aimed at astrophysical applications has been developed. The new code and algorithms are presented along with a comprehensive suite of test problems in one, two, and three dimensions. The new code is shown to be robust and accurate, equalling or improving upon a set of comparison codes. Fyris Alpha will be made freely available to the scientific community.     

Dynamic evolution of a quasi-spherical general polytropic magnetofluid with self-gravity

In various astrophysical contexts, we analyze self-similar behaviours of magnetohydrodynamic (MHD) evolution of a quasi-spherical polytropic magnetized gas under self-gravity with the specific entropy conserved along streamlines. In particular, this MHD model analysis frees the scaling parameter n in the conventional polytropic self-similar transformation from the constraint of n+γ=2 with γ being the polytropic index and therefore substantially generalizes earlier analysis results on polytropic gas dynamics that has a constant specific entropy everywhere in space at all time. On the basis of the self-similar nonlinear MHD ordinary differential equations, we examine behaviours of the magnetosonic critical curves, the MHD shock conditions, and various asymptotic solutions. We then construct global semi-complete self-similar MHD solutions using a combination of analytical and numerical means and indicate plausible astrophysical applications of these magnetized flow solutions with or without MHD shocks.     

RAPID: A fast,high resolution,flux-conservative algorithm designed for planet–disk interactions

We describe a newly developed hydrodynamic code for studying accretion disk processes. The numerical method uses a finite volume, non-linear, Total Variation Diminishing (TVD) scheme to capture shocks and control spurious oscillations. It is second-order accurate in time and space and makes use of a FARGO-type algorithm to alleviate Courant–Friedrichs–Lewy time step restrictions imposed by the rapidly rotating inner disk region. OpenMP directives are implemented enabling faster computations on shared-memory, multi-processor machines. The resulting code is simple, fast and memory efficient. We discuss the relevant details of the numerical method and provide results of the code’s performance on standard test problems. We also include a detailed examination of the code’s performance on planetary disk–planet interactions. We show that the results produced on the standard problem setup are consistent with a wide variety of other codes.