修正的粘滞律及有磁吸积盘的稳定性研究

本文采用修正的粘滞定律及磁流体力学研究了薄吸积盘内区及外区的稳定性问题。运用微扰方法导出了色散方程，分析了四种情况下吸积盘的不稳定性，结果表明：在同时考虑磁场和修正的粘滞律时，吸积盘中存在着三种振荡模式，其中粘滞模式总是稳定的，磁声速模式（包括向里、向外传播两种模式）通常是不稳定的。这些结果为解释ＢＬ Ｌａｃ天体、Ｓｅｙｆｅｒｔ星系、类星体等活动星系核的光变现象提供了理论依据。     

Current‐driven instabilities in weakly ionized disks

Cool weakly ionized gaseous rotating disk form the basis for many models in astrophysics objects. Instabilities against perturbations in such disks play an important role in the theory of the formation of stars and planets. Traditionally, axisymmetric magnetohydrodynamic (MHD) and recently Hall‐MHD instabilities have been thoroughly studied as providers of an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, the Hall instability against axisymmetric perturbations in incompressible rotating fluid in external poloidal and toroidal magnetic field is considered. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Disk and secondary irradiation in dwarf novae

After a brief review of the disk instability model for dwarf novae, I describe the effects of irradiation of the disc and the secondary (truncation of the inner disk and changes in the S-curve), and present numerical results showing that, as proposed by Brian Warner at the Wyoming conference (Warner, 1998) “The spectrum of post-superoutburst behaviour is generated by the interplay of reactions of the disk and secondary to irradiation”. Interestingly, we are also able to reproduce long lasting outbursts, which are similar to superoutbursts, without assuming the existence of a tidal instability.     

Energy Source for Shingle-like Structures in Spiral-arms

The filaments or shingles observed in neighbouring spiral-arms were found in a previous paper to have a wavelength of ˜ 1.3 kpc. The life-time of the shingles was estimated at ˜ 10⁸ years. The energy source to regenerate the shingles is sought in terms of a JEANS -type instability modified by a two stream instability between population I and disk material. From the dispersion relation it is found that an instability of this type occurs at a critical-wavelength of 1.9 kpc.     

Towards a new standard model for black hole accretion

We briefly review recent developments in black hole accretion disk theory, emphasizing the vital role played by magnetohydrodynamic (MHD) stresses in transporting angular momentum. The apparent universality of accretion-related outflow phenomena is a strong indicator that large-scale MHD torques facilitate vertical transport of angular momentum. This leads to an enhanced overall rate of angular momentum transport and allows accretion of matter to proceed at an interesting rate. Furthermore, we argue that when vertical transport is important, the radial structure of the accretion disk is modified at small radii and this affects the disk emission spectrum. We present a simple model demonstrating how energetic, magnetically-driven outflows modify the emergent disk emission spectrum with respect to that predicted by standard accretion disk theory. A comparison of the predicted spectra against observations of quasar spectral energy distributions suggests that mass accretion rates inferred using the standard disk model may be severely underestimated.     

Self-Collimated Stationary Disk Winds

In this contribution, we first review the theory of self-collimated jets launched from magnetized accretion disks (disk-winds originating from the first AUs). We show why it is crucial to solve in a self-consistent way the interplay between the resistive accretion disk and the ideal MHD jets. Indeed, this is the only way to get exact values for the disk ejection efficiency ξ (the jet mass load issue). Then, we show self-similar calculations of such accretion-ejection structures: first cold jets, then warm jets obtained in the presence of a hot disk chromosphere. Finally, we present for the first time an accretion-ejection flow crossing all three critical points.     

The study of magnetic reconnection in solar spicules

This work is devoted to study the magnetic reconnection instability under solar spicule conditions. Numerical study of the resistive tearing instability in a current sheet is presented by considering the magnetohydrodynamic (MHD) framework. To investigate the effect of this instability in a stratified atmosphere of solar spicules, we solve linear and non-ideal MHD equations in the x?z plane. In the linear analysis it is assumed that resistivity is only important within the current sheet, and the exponential growth of energies takes place faster as plasma resistivity increases. We are interested to see the occurrence of magnetic reconnection during the lifetime of a typical solar spicule.     

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 ionization-cone structures in active galactic nuclei: II. Nonlinear hydrodynamic modelling

In Part I of this paper (Paper I) we described an equilibrium model of a jet in the gravitational field corresponding to the rigid-rotation region of the galactic disk. We used linear stability analysis to find the waveguide-resonance instability of internal gravity waves due to the superreflection of these waves from the jet boundary. In this part of the paper, we perform nonlinear numerical 2D and 3D simulations of the development of this instability. We show that the shocks produced by this instability in the ambient medium of the jet are localized inside a cone with a large opening angle and are capable of producing features that are morphologically similar to those observed in galaxies with active nuclei (NGC 5252 for example).     

Loop models of solar flares: Revisions and comparisons

Due to developments in solar flare observations which appear to show that a particular class of solar flares result from instabilities occurring in magnetic loops we re-examine the Alfvén-Carlqvist flare model to show that it is workable and we update the Spicer loop model of a flare. It is noted that the Alfvén-Carlqvist model of necessity requires an external current driver which must maintain the current driven instability at marginal stability during the duration of the flare. In addition, it is argued that if the Alfvén-Carlqvist model is to work the current density must rise in a time shorter than an MHD or resistive tearing mode time scale. Otherwise, the dominant flare mechanism must be an ideal MHD or tearing type instability. Further, the distinctions between the two models are highlighted and a new hybrid model of the Alfvén-Carlqvist and Spicer models is introduced.     

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.     

Gravitational instability in isotropic MHD plasma waves

The effect of compressive viscosity, thermal conductivity and radiative heat-loss functions on the gravitational instability of infinitely extended homogeneous MHD plasma has been investigated. By taking in account these parameters we developed the six-order dispersion relation for magnetohydrodynamic (MHD) waves propagating in a homogeneous and isotropic plasma. The general dispersion relation has been developed from set of linearized basic equations and solved analytically to analyse the conditions of instability and instability of self-gravitating plasma embedded in a constant magnetic field. Our result shows that the presence of viscosity and thermal conductivity in a strong magnetic field substantially modifies the fundamental Jeans criterion of gravitational instability.     

X射线脉冲星周期的变化与吸积盘边界的K-H不稳定性

本文在中子星磁层与吸积盘之间引入了一个速度、密度、压强和磁场都连续变化的有限厚度的剪切层,以代替Anzer理论中的切向间断面,用磁流体力学方法讨论了中子星磁层与吸积盘交界处等离子体可压缩情况下平面波扰动的K-H不稳定性。结果表明,K-H不稳定性依然存在,径向波矢扰动成为不稳定的主要模式。文中特别讨论了剪切层厚度取值对中子星自转的影响,表明适当调节剪切层厚度就可解释X射线脉冲星周期的变化。将此模型应用到脉冲X射线源Her X-1上,得到较好的结果。     

Near-IR [Fe ii] Emission Diagnostics: Disk Winds in Young Stars

We investigate in this contribution emissivity properties of the [Fe ii] near-IR lines and their application to the study of jets in young stars. We report a model-independent method to determine the gas phase abundance in jets. Then, we present synthetic maps deduced from a recent self-consistent MHD cold disk wind model from which the thermal solution has been solved a posteriori. Kinematics is compared with observations in the [Fe ii] near-IR lines of L1551-IRS5. General kinematical behaviour is well reproduced, the existence of two velocity components in particular. However, we note a relative deficit in emission at intermediate velocity as observed in the jets from optically visible T Tauri stars. We examine predicted rotation signatures for various types of MHD disk winds. Compared to observations of DG Tau, warm disk winds reproduce quite well observed velocity shifts, whereas cold solutions predict too high toroidal velocities by a factor 2–3.     

eta Carinae: Binarity Confirmed

Gas giant planets have been detected in orbit around an increasing number of nearby stars. Two theories have been advanced for the formation of such planets: core accretion and disk instability. Core accretion, the generally accepted mechanism, requires several million years or more to form a gas giant planet in a protoplanetary disk like the solar nebula. Disk instability, on the other hand, can form a gas giant protoplanet in a few hundred years. However, disk instability has previously been thought to be important only in relatively massive disks. New three-dimensional, "locally isothermal," hydrodynamical models without velocity damping show that a disk instability can form Jupiter-mass clumps, even in a disk with a mass (0.091 M middle dot in circle within 20 AU) low enough to be in the range inferred for the solar nebula. The clumps form with initially eccentric orbits, and their survival will depend on their ability to contract to higher densities before they can be tidally disrupted at successive periastrons. Because the disk mass in these models is comparable to that apparently required for the core accretion mechanism to operate, the models imply that disk instability could obviate the core accretion mechanism in the solar nebula and elsewhere.     

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.     

耀斑大尺度电流片的湍流耗散和能谱分析

磁重联在宇宙的许多动力学现象中都是非常核心的过程.磁流体动力学(MHD)数值模拟是研究磁重联过程以及相应物理图像的一种很有效的手段.通过不同的参数组合,来研究MHD数值模拟中磁雷诺数和空间分辨率对磁重联率、数值耗散和能谱分布的影响.对得到的数据进行分析后,发现磁雷诺数对磁重联率和能谱分布有一定的影响.磁雷诺数越大,磁重联过程进入非线性阶段所需的特征时间越短,磁重联率就越早发生跃升.磁雷诺数R_m对耗散开始发挥作用的Kolmogorov微观尺度l_ko有明显影响:R_m越大,l_ko就越小.研究了磁重联过程中包括数值耗散在内的额外耗散对重联过程的影响.结果表明,撕裂模不稳定性开始之前的额外耗散以纯数值耗散为主,撕裂模不稳定性出现之后,额外耗散出现同步跃升,说明不稳定性导致的湍流明显增强了耗散的效果,相当于在局部湍流区引入了超电阻.能谱分析进一步表明,大尺度电流片的l_ko完全可能出现在宏观的MHD尺度上.     

Kink unstable coronal loops: current sheets,current saturation and magnetic reconnection

The kink instability in a coronal loop is a possible explanation of a compact loop flare as it may cause a current sheet to form allowing reconnection to take place and release the free magnetic energy stored in the loop. However, current sheets do not form in all cases. Ali and Sneyd (2001) investigated three different classes of equilibrium (determined by the form of the twist) using a magneto-frictional code. They searched for the equilibria to which the loop might evolve once it had become unstable to the kink instability. They found indications of current-sheet formation for only one class of equilibrium studied. However, as they pointed out, since their code searched for equilibria they were unable to say for certain that the loop would evolve in this way. In this paper we have considered the same three classes of equilibria but have used a code which follows the non-linear 3D MHD (magnetohydrodynamic) evolution of the loop. We have investigated whether or not there are indications of current-sheet formation. In the cases where there is evidence of this we have found that reconnection does occur and releases sufficient magnetic energy to explain a compact loop flare.     

Cosmic ray scattering in compressible turbulence

We study the scattering of low-energy cosmic rays (CRs) in a turbulent, compressive magnetohydrodynamic (MHD) fluid. We show that compressible MHD modes – fast or slow waves with wavelengths smaller than CR mean free paths induce cyclotron instability in CRs. The instability feeds the new small-scale Alfvénic wave component with wavevectors mostly along magnetic field, which is not a part of the MHD turbulence cascade. This new component gives feedback on the instability through decreasing the CR mean free path. We show that the ambient turbulence fully suppresses the instability at large scales, while wave steepening constrains the amplitude of the waves at small scales. We provide the energy spectrum of the plane-parallel Alfvénic component and calculate mean free paths of CRs as a function of their energy. We find that for the typical parameters of turbulence in the interstellar medium and in the intercluster medium the new Alfvénic component provides the scattering of the low-energy CRs that exceeds the direct resonance scattering by MHD modes. This solves the problem of insufficient scattering of low-energy CRs in the turbulent interstellar or intracluster medium that was reported in the literature.