与耀斑双带重合的一对共轭电流带

利用SDO (Solar Dynamics Observatory)/HMI (Helioseismic and Magnetic Imager)观测到的矢量磁图,研究了与活动区AR12673上爆发的一个X9.3级耀斑(2017年9月6日)的相关电流分布和演化.结果显示,在该活动区的磁中性线两边存在一对方向相反的电流密度约为0.4 A/m~2的长电流带,可称其为一对共轭电流带.这对共轭电流带在耀斑发生之前、期间以及之后一直存在;并且观测到,该耀斑的两个亮带的位置几乎刚好与两个电流带重叠,它们的形状也极其相似. 9月6日电流总强度演化曲线表明,电流强度在X9.3级强耀斑爆发期间出现快速增加的现象,这种现象持续了几个小时.这一研究结果有力支持了磁准分界面(Quasi-Separatrix Layer, QSL) 3维重联模型.     

太阳爆发过程中的大尺度磁重联电流片:理论和观测

从理论和观测两个方面来介绍和讨论出现在太阳爆发过程中的磁重联电流片及其物理本质和动力学特征。首先介绍在理论研究和理论模型中,磁重联电流片是如何在爆发磁结构当中形成并发展的,对观测研究有什么指导意义。然后介绍观测工作是从哪几个方面对理论模型预测的电流片进行证认和研究的。第三,将介绍观测研究给出了哪些过去所没有能够预期的结果,这些结果对深入研究耀斑一CME电流片以及其中的磁重联过程的理论工作有什么重要的、挑战性的意义。第四,讨论最新的与此有关的理论研究和数值实验。最后,对未来的研究方向和重要课题进行综述和展望。     

日冕中无衰减振荡的研究进展

日冕中冕环的无衰减横向振荡(简称无衰减振荡)自2012年被发现以来,受到了广泛的关注.无衰减振荡具有在日冕中广泛存在以及振幅无明显变小的特性,使之在解释日冕加热和冕震学诊断上都具有相当的潜力.总结了日冕中无衰减振荡的研究进展,包括观测研究得到的一系列结果、提出的理论和数值模型以及基于无衰减振荡进行冕震学诊断的一些尝试,并且展望了未来可以进一步开展的研究.     

耀斑环动力学的数值研究

数值模拟了太阳耀斑环动力学演化的二维磁重联过程。结果表明,在重联Ｘ 点比较高时,演化过程能再现双带耀斑中的耀斑环活动;而重联Ｘ 点比较低时,演化情况能解释致密耀斑的观测特征     

卫星黑子衰减触发的喷流事件

详细分析了一次太阳低层大气磁场重联触发的喷流事件.这次喷流发生在2014年8月1日,爆发自美国国家海洋和大气管理局(National Oceanic and Atmospheric Administration, NOAA)活动区12127边缘的一个卫星黑子处.该喷流爆发包括日浪、紫外喷流、极紫外高温和低温喷流.大熊湖太阳天文台(Big Bear Solar Observatory,BBSO)的Goode Solar Telescope (GST)高分辨率氧化钛(TiO)谱线的光球观测显示,喷流爆发过程中,卫星黑子一直衰减.到喷流结束,卫星黑子面积共减少了80%.在此过程中,太阳动力学天文台(Solar Dynamics Observatory, SDO)日球磁场成像仪(Helioseismic and Magnetic Imager, HMI)的视向磁场观测表明,该卫星黑子对应的负极磁场与相邻的正极磁场发生明显对消,产生喷流足部亮点.根据SDO卫星太阳大气成像仪(Atmospheric Imaging Assembly, AIA)的多波段观测,该足部亮点首先出现在紫外1600?波段.待紫外(1600?)喷流从紫外足部亮点顶部向上喷发,在极紫外波段也观测到相应的亮源.随着足点源亮度突然增强,有明显的极紫外低温喷流和日浪从足部亮点侧面喷发.从GST的高分辨率Hα图像上,可见日浪由许多精细纤维组成,这些纤维扎根在足点源的东南侧.根据从光球层过色球层再到日冕层的多波段高分辨率观测,色球中下层的磁场对消触发了这次喷流事件.向上喷发的物质流可以携带能量进入上层大气,并加热上层大气.研究结果表明,低层大气磁重联可能对解决日冕加热问题起重要作用.     

太阳耀斑磁重联的数值模拟

数值模拟了太阳耀斑中二维磁重联过程。结果表明,当重联Ｘ 点比较高时,演化过程能再现双带耀斑中耀斑环的运动等主要特征;当重联Ｘ 点比较低时,可解释致密耀斑的观测特征。结果还表明,耀斑环上升和重联点上升之间没有直接的联系。     

太阳光球层磁重联的直接证据

本文首次给出了发生在太阳光球磁重联的一个直接的观测证据。 这一磁重联的观测特征是:(1)重联发生在一新浮现磁通量区的一极与极性相反的老磁通量之间;(2)重联前中性线附近磁剪切明显;(3)被重联两极为一对消磁结构,重联发生在稳定的磁通量损失数小时之后;(4)一个级别为C2.9的亚耀斑发生在重联之前。该耀斑以重联区为中心,双带离重联位置2～3万公里,直到耀斑极大相后14分钟,重联仍未发生;(5)重联后,磁对消速率呈增大趋势。     

太阳耀斑区磁场和电流研究的一些新进展

在太阳耀斑区磁场和电流研究方面，文中将着重介绍太阳横向磁场方位的确定，太阳活动区磁场的非热性表示、太阳耀斑前后的活动区磁场变化、以及耀斑核块与活动区纵向电流密度极大点位置的关系等几个重要问题的研究现状。     

日冕三重无力场电流片的磁场重联

采用二维三分量磁流体力学模型,对日冕三重无力场电流片的磁场重联进行了数值研究,揭示了重联过程的基本物理特征．这类重联过程将加热和加速日冕等离子体,并导致多个高温、高密度、高磁螺度的磁岛的形成和向上喷发．这表明,多重无力场电流片的重联可能在日冕磁能释放、上行等离子体团的形成和太阳磁场螺度向行星际空间的逃逸方面起重要的作用．     

日冕三重无力场电流片的磁场重联

采用二维三分量磁流体力学模型,对日冕三重无力场电流片的磁场重磁联进行了数值研究,揭示了重联过程的基本物理特征,这类重联过程将加热和加速日冕等离子体,并导致多个高温、高密度、高磁螺度的磁岛的形成和向上喷发,这表明,多重无力场电流片的重联可能在日冕磁能释放、上行等离子体的形成和太阳磁场螺度向行星际空间的逃逸方面起重要的作用。     

Turbulence Diffusion and Energy Spectrum Analysis of Large-scale Current Sheets in Flares

Magnetic reconnection (MR) is one of the most important physical processes for many dynamical phenomena in the universe. Magnetohydrodynamical (MHD) simulation is an effective way to study the MR process and the physical pictures related to the MR. With different parameter setups, we investigate the influences of the Magnetic Reynolds number and spatial resolution on the reconnection rate, numerical dissipation, and energy spectrum distribution in the MHD simulation. We have found that the magnetic Reynolds number $R_m$ has definite impact on the reconnection rate and energy spectrum distribution. The characteristic time for entering into the non-linear phase will be earlier as the Reynolds number increases. When it comes to the tearing phase, the reconnection rate will increase rapidly. On the other hand, the magnetic Reynolds number affects significantly the Kolmogorov microscopic scale $l_{\mathit{ko}}$, which becomes smaller as $R_m$ increases. An extra dissipation is defined as the combined effect of the numerical diffusion and turbulence dissipation. It is shown that the extra dissipation is dominated by the numerical diffusion before the tearing mode instability takes place. After the instability develops, the extra dissipation rises vastly, which indicates that turbulence caused by the instability can enhance the diffusion obviously. Furthermore, the energy spectrum analysis indicates that $l_{\mathit{ko}}$ of the large-scale current sheet may appear at a macroscopic MHD scale very possibly.     

Alleviation of catastrophic quenching in solar dynamo model with nonlocal alpha‐effect

The nonlocal alpha‐effect of Babcock‐Leighton type is not prone to the catastrophic quenching due to conservation of magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha‐effect, the diamagnetic pumping, and dynamical equation for the magnetic alpha‐effect. The same model shows catastrophic quenching when the alpha‐effect is changed to its local formulation. The nonlocal model shows a preferred excitation of magnetic fields of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal‐to‐polar fields ratio of about a thousand (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The origin of cutoff frequencies for torsional tube waves propagating in the solar atmosphere

Torsional waves supported by magnetic flux tubes have long been thought to bear a high potential for supplying energy and momentum to the upper solar atmosphere, thereby contributing to its heating and to the driving of dynamic events like spicules. This hope rested on the belief that their propagation is not impeded by cutoff restrictions, unlike longitudinal and kink waves. We point out that this applies only to thin, isothermal tubes. When they widen in the chromosphere, and as a result of temperature gradients, cutoff restrictions arise. We compare them to recent observational reports of such waves and of vortex motions and find that their long period components are already affected by cutoff restrictions. An observational strategy is proposed that should permit the derivation of better information on vortex flows from off‐center observations with next generation telescopes (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convective dynamos in spherical wedge geometry

Self‐consistent convective dynamo simulations in wedge‐shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large‐scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solar dynamo models with α ‐effect and turbulent pumping from local 3D convection calculations

Results from kinematic solar dynamo models employing α ‐effect and turbulent pumping from local convection calculations are presented. We estimate the magnitude of these effects to be around 2–3 m s^–1, having scaled the local quantities with the convective velocity at the bottom of the convection zone from a solar mixing‐length model. Rotation profile of the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface shear layer on the dynamo solutions. With these choices of the small‐ and large‐scale velocity fields, we obtain estimate of the ratio of the two induction effects, C _α /C _Ω ≈ 10^–3, which we keep fixed in all models. We also include a one‐cell meridional circulation pattern having a magnitude of 10–20 m s^–1 near the surface and 1–2 m s^–1 at the bottom of the convection zone. The model essentially represents a distributed turbulent dynamo, as the α ‐effect is nonzero throughout the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of latitude. Turbulent pumping of the mean fields is predominantly down‐ and equatorward. The anisotropies in the turbulent diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed negative sign of B _r B _ϕ . Although the activity clearly shifts towards the equator in comparison to previous models due to the combined action of the α ‐effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of the activity still occurs at too high latitudes (between 5° … 60°). Other problems include the relatively narrow parameter space within which the preferred solution is dipolar (A0), and the somewhat too short cycle lengths of the solar‐type solutions. The role of the surface shear layer is found to be important only in the case where the α ‐effect has an appreciable magnitude near the surface. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age‐old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucial role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large‐scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean‐field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some current issues of solar simulations are put forward (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling a Maunder minimum

We introduce on/off intermittency into a mean field dynamo model by imposing stochastic fluctuations in either the alpha effect or through the inclusion of a fluctuating electromotive force. Sufficiently strong small scale fluctuations with time scales of the order of 0.3–3 years can produce long term variations in the system on time scales of the order of hundreds of years. However, global suppression of magnetic activity in both hemispheres at once was not observed. The variation of the magnetic field does not resemble that of the sunspot number, but is more reminiscent of the ¹⁰Be record. The interpretation of our results focuses attention on the connection between the level of magnetic activity and the sunspot number, an issue that must be elucidated if long term solar effects are to be well understood. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

3-D non-linear evolution of a magnetic flux tube in a spherical shell: The isentropic case

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code.We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)