1997年1月空间天气事件的三维MHD模拟

以1997年1月空气天气事件期间的观测为依据，在构造了比较接近真实的背景太阳风基础上，进一步利用三维时变的MHD模式，模拟了CME（日冕物质抛射）激发的扰动在行星际空间的传播过程，对地球空间环境的影响及行星际磁场南向分量Bz在1AU的时间经历。模拟结果与WIND卫星的测量进行了比较。结果表明，模拟与观测得到的扰动得到地球的时间、地球空间环境各量的变化及Bz的时间经历基本一致。     

1998年5月日冕三维结构的数值模拟

采用三维理想磙流体力学（MHD）模式，内边界条件把二维投影特征线边界方法推广应用到三维计算，有效地稳定了数值计算并保证稳态解的自洽性；初始猜解磁场由1935卡林顿周光球磁场观测数据得到，这样计算得到的1998年5月份期间日冕三维结构比较符合实际，计算结果表明：（1）计算得到的源表面非径向磁场量值在磁中性线附近不超过2μT，表明源表面磁场基本径向。（2）模拟得到的源表面径向磁场量值除了在磁中性线附近的区域外变化不大，这和观测一致。（3）由源表面磁场按平方反比的规律计算出1AU处磁场量值更接近观测值。（4）计算得到的日冕结构和观测定性一致，三维数值模拟结果表明，日冕的三维大尺度背景结构主要是由磁场决定的，在闭磁场处或者电流片附近，太阳风的密度高，速度低；在开场区，太阳风的密度低，速度高。     

嵌套闭磁场日冕背景结构对触发CME特征的影响

在子午面内,研究具有嵌套闭磁场结构冕流背景对触发日冕物质抛射(CME)特征的影响.在冕流背景磁场结构内有三个小尺度的闭磁场结构,其中中间小尺度闭磁场结构的磁场方向和冕流整体偶极磁场方向相反.CME触发模型在这个小尺度闭磁场结构下方浮出,它具有同心圆形磁场结构,半径为a=0.1Rs(Rs为太阳半径);CME触发模型前半部分磁场方向和触发位置处小尺度闭磁场结构的磁场方向相反,与冕流整体偶极磁场方向相同.数值模拟结果表明,当CME触发模型中心等离子压力与边界压力之比m≥2时,上浮CME触发模型可以触发CME;当m＜2时,上浮CME触发模型不能触发CME,计算结论的误差小于1%.      

Superlow-Frequency MHD-Eigenoscillations With a Discrete Spectrum in Parabolic Model of Magnetosphere

A new concept of the global magnetospheric resonator is suggested for fast magnetosonic waves in which the role of the resonator is played by the near part of the plasma sheet. It is shown that the magnetosonic wave is confined in this region of the magnetosphere throughout its boundaries. The representative value of the resonator's eigenfrequency estimated at f ≈ 1 mHz is in good agreement with observational data on ultraThe theory explains the ground-based localization of the oscillations observed in the midnight-morning sector of the high-latitude magnetosphere.      

MHD waves in coronal flux tubes

The basic MHD waves of a coronal flux loop are investigated for the rectangular box model of a plasma with oblique magnetic field and line-tied at the ends. The waves found are completely different from those in a periodic box, representative for tokamaks. They consist of a mixture of Alfvén components with a ballooning factor, favouring minimal field line bending, and fast components without such a factor. Pure Alfvén modes are only found as singular limiting cases of cluster spectra of Alfvén-fast waves, where the fast components are localised in a photospheric boundary layer which is dictated by the requirements of line-tying. This justifies the assumption of continuous spectra in coronal loops, required for the mechanism of resonant Alfvén wave heating. The waves consist of large amplitude Alfvén components in the corona and fast components with a small but rapidly varying amplitude in the boundary layer, so that they appear to have the right signature for effective transfer of energy from the photosphere to the corona.     

A Self-Consistent Determination of the Temperature Profile and The Magnetic Field Geometry in Winds of Late-Type Stars

Cool giant and supergiant stars generally present low velocity winds with high mass-loss rates. Several models have been proposed to explain the acceleration process of these winds. Although dust is known to be present in these objects, the radiation pressure on these particles is uneffective in reproducing the observed physical parameters of the wind. The most promising acceleration mechanism cited in the literature is the transference of momentum and energy from Alfvén waves to the gas. Usually, these models consider the wind to be isothermal. We present a stellar wind model in which the Alfvén waves are used as the main acceleration mechanism, and determine the temperature profile by solving the energy equation taking into account both the radiative losses and the wave heating. We also determine, self-consistently, the magnetic field geometry as the result of the competition between the magnetic field and the thermal pressure gradient. As the main result, we show that the magnetic geometry presents a super-radial index in the region where the gas pressure is increasing. However, this super-radial index is greater than that observed for the solar corona.     

Fluid instabilities around accreting X-ray sources

Some aspects of fluid instabilities occurring in the magnetospheres of accreting neutron stars are discussed. It is pointed out that (i) in the absence of strong differential rotation, the accreting plasma should be drawn out into spiralling, sheet-like structures, resulting in efficient mixing between the two media; (ii) the Rayleigh-Taylor instability also acts to limit the X-ray luminosity in super-critical sources; and (iii) magnetic shear has a strong stabilizing influence on Kelvin-Helmholtz modes, and its presence may allow substantial amounts of material to be supported around the magnetosphere.     

Application of ADER Scheme in MHD Simulation

The Arbitrary accuracy Derivatives Riemann problem method(ADER) scheme is a new high order numerical scheme based on the concept of finite volume integration,and it is very easy to be extended up to any order of space and time accuracy by using a Taylor time expansion at the cell interface position.So far the approach has been applied successfully to flow mechanics problems.Our objective here is to carry out the extension of multidimensional ADER schemes to multidimensional MHD systems of conservation laws by calculating several MHD problems in one and two dimensions: (ⅰ) Brio-Wu shock tube problem,(ⅱ) Dai-Woodward shock tube problem,(ⅲ) Orszag-Tang MHD vortex problem.The numerical results prove that the ADER scheme possesses the ability to solve MHD problem,remains high order accuracy both in space and time,keeps precise in capturing the shock.Meanwhile,the compared tests show that the ADER scheme can restrain the oscillation and obtain the high order non-oscillatory result.     

双楔面诱导的斜爆轰波阵面的磁流体控制

基于高超声速飞行器爆轰推进发展需要,数值研究了由双楔面诱导的斜爆轰波类型及其磁流体(MHD)控制的可行性,探讨了后楔倾角变化对双楔斜爆轰流场与MHD控制的影响.结果表明:对于前、后楔倾角分别为15°与20°的双楔面诱导的斜爆轰,Lorentz力可以将不同来流马赫数条件下的稳定斜爆轰波阵面恢复到设计位置,并可以使不稳定斜爆轰波趋于稳定,但是,无法将不稳定斜爆轰波恢复到设计位置.而当后楔倾角稍微增加时,若爆轰流场稳定,同样可对斜爆轰波阵面进行有效控制.但当后楔倾角大于25°时,稳定斜爆轰波流场失稳,此时MHD无法控制流场稳定性.但若增加气体反应速率,虽然斜爆轰波不稳定,MHD控制却可使其趋于稳定,但无法将爆轰波阵面恢复到设计位置.      

Spectral analysis of 3D MHD models of coronal structures

We study extreme-ultraviolet emission line spectra derived from three-dimensional magnetohydrodynamic models of structures in the corona. In order to investigate the effects of increased magnetic activity at photospheric levels in a numerical experiment, a much higher magnetic flux density is applied at the photosphere as compared to the Sun. Thus, we can expect our results to highlight the differences between the Sun and more active, but still solar-like stars. We discuss signatures seen in extreme-ultraviolet emission lines synthesized from these models and compare them to observed signatures in the spatial distribution and temporal evolution of Doppler shifts in lines formed in the transition region and corona. This is of major interest to test the quality of the underlying magnetohydrodynamic model to heat the corona, i.e. currents in the corona driven by photospheric motions (flux braiding).