高速太阳风中α粒子分布的形成

以多级波粒相互作用模型和1维电磁混合模拟方法研究了高速太阳风中α粒子分布的形成机制．认为α粒子从其产生区到1AU,经历4个特征区域：在冕洞下方的小尺度磁重联区,α粒子产生并得到加速：在约十几个太阳半径日冕区,α粒子形成球壳分布函数：在几十个太阳半径处,在Alfven波和温度各向异性不稳定性的作用下,α粒子分布函数演化为漂移Maxwell分布：当α粒子进入0．3－1AU时,Alfven波的作用使α粒子速度超过太阳风速度约1个局域Alfven速度．     

非线性动力Alfven波的演化

本文讨论了非线性动力Alfven波在1＞＞β＞＞me／mi条件下形成孤子解的限制条件，表明只有亚Alfven速度才能形成压缩型孤子．由动力Alfven波孤子和湍流产生的极大纵向电子速度不会超过本底电子热速度，但由此产生的纵向电场在小尺度范围可以很强，这可能对太阳风的能量传输与粒子加速产生一定影响．     

磁云边界层中的复合重联喷流观测分析

太阳风中的磁场重联通常与行星际日冕物质抛射有关.本文分析了1995年10月18日WIND飞船观测到的一例磁云前边界层中的复合重联喷流事件.该复合排空区由相邻两个不同方向的喷流构成，这两个喷流分别经过Walén关系的证认，符合行星际磁场重联排空区等离子体喷流的特征.结果表明，在磁云前端可能存在众多重联点，从而将磁云本体的磁场剥离，形成比单一重联喷流区更复杂的三维边界层结构.磁云边界层中可能发生多点多次重联，从而不表现出单点重联的排空区特征，这可能是行星际磁场重联排空区较少在ICME前端被观测到的原因之一.      

磁重联区Alfven波及新生离子的加速

利用二维混合数值模拟研究了有速度驱动、低等离子体β值情况下的磁场重联过程，结果表明磁重联过程可以产生Alfven波，该Alfven波动对重新区中的新生离子作用，使得新生离子经历投掷角散射方程，具有球壳分布特征，部分新生离子得到加速，其获得的最大能量约为4（miVA0^2/2)，此加速过程所需的加速时间在100／Ωi量级，是一个极快的加速机制，加速粒子能谱为双幂律谱。     

日冕三维太阳风模拟

COIN-TVD MHD模型是近年发展起来的能有效实现日冕–行星际三维太阳风模拟的模型.本文利用此模型针对日冕区三维太阳风进行研究,为了模拟日冕太阳风的加热加速,对模型中的体积加热项做了调整.在磁流体模拟中,减小磁场散度的误差是关键问题之一,在调整体积加热项后应用扩散法、八波法、扩散八波法,对2199卡林顿周的背景太阳...     

太阳表面水平运动驱动的磁圈与开放磁场重联的模拟研究

太阳大气的诸多观测事件(如耀斑、喷流等)均被归因于磁重联产生的能量转换. 近年来, 关于太阳风起源, 有研究提出了磁重联使闭合磁圈开放为太阳风供应物质的新模式. 在该模式中, 闭合磁圈被光球超米粒组织对流携带, 向超米粒边界运动, 与位于边界的开放磁场相碰撞进而发生磁重联. 该模式中磁重联的驱动及其效应是本文的研究目标. 磁流体力学(MHD)数值模拟是研究太阳大气磁重联物理过程的重要途径. 本文建立了一个二维MHD数值模型, 结合太阳大气温度和密度的分层分布, 在超米粒组织尺度上模拟了水平流动驱动的闭合磁圈与开放磁场的重联过程. 通过对模拟结果的定量分析, 认为磁重联确实能够将闭合磁圈的物质释放, 进而供应给新的开放磁结构并产生向上流动. 该结果为进一步模拟研究太阳风初始外流奠定了基础.      

等离子体彗尾动力学的数值研究

本文基于可压缩磁流体动力学模型,数值研究了尾瓣巾具有超Alfven速流动的等离子体彗尾的动力学特征。结果表明,等离子体片和尾瓣之间的剪切等离子体流动将会激发流动撕裂模不稳定性,引起彗尾等离子体片中发生磁场重联,形成磁岛和高密度的等离子体团。进而模拟了太阳风引起的局部驱动力对等离子体彗尾中磁场重联的影响,其特征时间远大于流动撕裂模。我们认为一些观测到的等离子体彗尾中的四块和彗尾截断事件可能主要与彗尾中剪切等离子体流动所引起的流动撕裂模不稳定性有关。      

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

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

两冕流间CME事件数值模拟的改进

通过偶极子场和六极子场适当叠加，改进猜解磁场，使猜解磁场在太阳南北极符号相反，然后采用理想磁流体力学方程组(MHD)，由猜解磁场与太阳风流动相互作用计算出稳态自洽解，得到定性上与观测比较接近的具有两个冕流的背景结构．在两个冕流间采用具有同心圆磁场位形的触发模型触发CME事件，研究CME的日冕传播特征．模拟结果表明，CME被约束在两冕流间传播，CME闭磁场位形和磁云横截面磁场位形相似，可以解释1AU处观测磁云的部分特征；在CME附近，存在压力和Lorentz力起主要作用的区域，这可以为分析1AU处CME事件的观测数据提供帮助．     

太阳风中小尺度磁通量管边界重联的统计研究

通过对WIND卫星1995—2005年的数据,利用程序筛选及人工识别两种不同方法确定的小尺度磁通量管进行比较,发现程序筛选法中41%的小尺度磁通量管有边界重联现象,与人工识别法确定的小尺度磁通量管的统计结果接近;通过人工识别和程序筛选两种方法确定的小尺度磁通量管的边界重联特征,包括磁场剪切角、磁场强度以及重联耗散区的持续时间等,也具有相同的统计趋势.结果表明,两种方法确定的小尺度磁通量管在重联特性上没有本质区别,因此采用这两种方法得到的数据作为样本来统计小尺度磁通量管前后边界重联事件.本文共确定了71个重联事件,统计结果显示有50个(70%)重联耗散区磁场的减小超过20%,47个(66%)磁场剪切角大于90°;多数重联事件的磁场剪切角大于90°,表明小尺度磁通量管边界中主要发生的是反平行重联.将小尺度磁通量管的前后边界重联分开进行统计,结果显示其前后边界重联的特征是相似的,与磁云前后边界存在差异的性质不同,这意味着太阳风中的小尺度磁通量管并不具有磁云这种大尺度磁通量管的膨胀特征.     

Combining MHD and kinetic modelling of solar flares

Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport in converging magnetic fields.     

Search for a chromospheric resonator above sunspots

We have recently observed thirteen sunspots, both with CDS on SOHO and TRACE. For each sunspot we investigate the oscillations in the chromosphere, transition region and corona. Above the umbra the power spectra show one dominant peak corresponding to a period of 3 minutes. The intensity oscillation amplitude increases with increasing temperature, reaches a maximum for lines emitting close to 1–2 × 10⁵ K, and decreases for higher temperatures. Part of the wave energy penetrates into the corona, channeled into small areas that appear to coincide with the endpoints of sunspot coronal loops. The observations support the hypothesis that the oscillations are caused by upwardly propagating acoustic waves and show no signs of the resonances, equally spaced ≈1 mHz in frequency, predicted by the chromospheric resonator theories.     

Energetic electrons in impulsive solar flares: Radio diagnostics

Radio emissions during and outside solar flares are tracers of energetic electrons from the bottom of the corona to the interplanetary space. This review focusses on impulsive flares, where joint analyses of radio, hard X-ray and γ-ray observations proved to be powerful probes of the properties of accelerated electrons and of the sites in the corona where they are accelerated. Evidence of electron acceleration and transport in the corona from microwave imaging and decimetre wave spectroscopy is reviewed and compared, and recent work on the interpretation of microwave spectra in terms of energetic electron spectra is discussed. The two directions for future instrumentation are the extension to shorter wavelengths, with the aim of probing relativistic electrons, and solar dedicated spectral imaging from centimetric to metric waves to provide a unified view of the acceleration signatures that stem so far from different instruments with either spectroscopic or imaging capabilities.     

Post-flare thermal waves in the solar corona

While imaging giant post-flare arches in the solar corona, the Hard X-Ray Spectrometer aboard the SMM detected thermal disturbances propagating through the corona after two-ribbon flares. The speed of propagation is close to, or below, 10 km s^?1, and no obvious time-variation of the speed is indicated in the HXIS data. For subsequent two-ribbon flares in the same active region, these thermal disturbances (waves) exhibit highly homologous properties; thus the waves appear to propagate through preexisting arches formed after earlier flares. Temperatures of > 20 × 10⁶ K have been detected in these moving phenomena. We suggest that we see here in X-rays upper products of the consecutive reconnections which create the post-flare loops below. Temperature maps in fine field of view of HXIS offer now a new possibility to detect postflare arches in the corona built during two-ribbon flares.     

不同波段太阳射电爆发的准周期振荡

本文对1989年4月9日与4B级大耀斑伴生的两个波段(3750 MHz,2840 MHz)的射电大爆发进行了简单的初步分析。结果表明:这两个波段的爆发不仅各自具有繁多的周期,而且都存在30+3s的准周期振荡。这种现象可能与磁环的准周期振荡有关,或与粒子流有关。      

Nonlinear development of MHD waves in coronal loops

The solar corona, modelled by a low β, resistive plasma slab sustains MHD wave propagations due to footpoint motions in the photosphere. The numerical simulation presents the evolution of MHD waves and the formation of current sheets. Steep gradients at the slab edges, which are signatures of resonance layers are observed. Singularities are removed by the inclusion of finite resistivity. The fast waves develop kink modes. As the plasma evolves the current sheets which provide heating at the edges fragment into two current sheets at each edge which in turn come closer when the twist is enhanced.     

Solar atmosphere wave dynamics generated by solar global oscillating eigenmodes

The solar atmosphere exhibits a diverse range of wave phenomena, where one of the earliest discovered was the five-minute global acoustic oscillation, also referred to as the p-mode. The analysis of wave propagation in the solar atmosphere may be used as a diagnostic tool to estimate accurately the physical characteristics of the Sun’s atmospheric layers.In this paper, we investigate the dynamics and upward propagation of waves which are generated by the solar global eigenmodes. We report on a series of hydrodynamic simulations of a realistically stratified model of the solar atmosphere representing its lower region from the photosphere to low corona. With the objective of modelling atmospheric perturbations, propagating from the photosphere into the chromosphere, transition region and low corona, generated by the photospheric global oscillations the simulations use photospheric drivers mimicking the solar p-modes. The drivers are spatially structured harmonics across the computational box parallel to the solar surface. The drivers perturb the atmosphere at 0.5?Mm above the bottom boundary of the model and are placed coincident with the location of the temperature minimum. A combination of the VALIIIC and McWhirter solar atmospheres are used as the background equilibrium model.We report how synthetic photospheric oscillations may manifest in a magnetic field free model of the quiet Sun. To carry out the simulations, we employed the magnetohydrodynamics code, SMAUG (Sheffield MHD Accelerated Using GPUs).Our results show that the amount of energy propagating into the solar atmosphere is consistent with a model of solar global oscillations described by Taroyan and Erdélyi (2008) using the Klein-Gordon equation. The computed results indicate a power law which is compared to observations reported by Ireland et al. (2015) using data from the Solar Dynamics Observatory/Atmospheric Imaging Assembly.     

太阳大气中非线性波的数字计算

本文用多方过程方程式封闭一维带引力场的流体力学方程组,用不同的γ值计算了不同类型的扰动在HSRA太阳大气模型中向上传播和演化。说明在太阳大气中向上传播的小扰动能够发展成为激波;激波的形成高度基本上不受热力学过程的影响,而与扰动强度关系很大,愈小的扰动传播愈高;但由于太阳大气中湍流场的存在,完全有可能使这样小的扰动在未形成激波前就被破坏掉。这就证实了一些人认为的对流层机械波并不是一个对色球和日冕的有效加热机制的说法。      

内重力波传播的3维传输函数模式研究

在考虑背景风场及大气耗散的条件下,建立了3维内重力波传输函数数值模式.分析了300 km高度3维传输函数在频率波数域的特性,并以近地面单位脉冲点源为激发源,得到了内重力波在3维空间中的时空分布.讨论了不同时空尺度地面方波源激发的内重力波在电离层高度的能量分布特征.结果表明,(1)对内重力波而言,背景大气相当于一个带通滤波器,只有波动周期和波长分别在15～30 min和200～400 km之间的重力波扰动最容易上传到300km高度;(2)在背景风场的作用下水平面上以同心圆扩散的波阵面以及垂直方向上成漏斗状的波阵面发生了变形,并且逆风方向比顺风方向更有利于声重力波由对流层向电离层高度传播;(3)300km高度对时间尺度和空间尺度分别在20～30 min和150～250 km之间的地面方波源响应的总能量最强.      

无碰撞电流片哨声波不稳定性

通过数值求解文献[4]中物理模型A得到的一般形式的色散关系,讨论了无碰撞电流片低频波不稳定性问题.结果表明,哨声波能被无碰撞电流片直接激发.在中性片上(z/di=0),在较宽的波数范围内,斜哨声波是可以传播的,但它基本上是稳定的.在离子惯性区内(z/di<1,电子惯性区外),斜传播的哨声波是不稳定的.在离子惯性区边缘(z/di=1),斜传播的哨声波仍然是不稳定的,增长率更大,不稳定的波频率范围更高.此外,朝向中性片方向传播(kzdi<0)的哨声波比离开中性片方向传播(kzdi>0)的哨声波有更大的增长率.