多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

SIP-AMR-CESE模式对CR2055背景太阳风的模拟

使用三维太阳行星际自适应网格守恒元解元太阳风模型（SIP-AMR-CESE MHD）,模拟从太阳表面到地球轨道附近的太阳风.该模型使用六片网格技术,同时利用PARAMESH软件包实现网格自适应.在该模型的基础上,通过增加广义拉格朗日乘子（GLM）磁场散度误差消去方法,完善网格加密放粗判据,微调加速加热形式等方法,使模拟结果与观测更好地符合.另外,通过控制不同时刻的计算区域,显著提高了模型的计算效率.在此基础上,给出了模型改进后模拟得到的CR2055太阳风稳态解与观测的对比分析.      

上游太阳风中高能电子向同步轨道区的传输

通常认为，同步轨道区的电子通量增加是由于磁暴或者上游太阳风高速流的扰动所引起．近来的观测表明，起源于太阳活动的行星际高能电子也是引起同步轨道电子通量增加的重要原因之一．Zhao等在研究2000年7月14日太阳剧烈活动时发现，同步轨道区相对论电子通量巨幅增加时没有观察到上游太阳风高速流的扰动，并且磁暴发生在电子通量事件之后．采用解析磁场模型和实际磁场模型(T96模型)模拟来自太阳的相对论电子在磁尾中的运动特性．计算结果表明，当行星际磁场南向时，进入到磁尾的行星际相对论电子可以从较远的磁尾区域运动到同步轨道区域．这一研究结果从理论上论证了起源于太阳活动的高能电子可以对同步轨道区相对论电子通量的增加产生重要的作用．     

1998年10月18——-19日磁暴主相的行星际源分析

采用高时间分辨率的地磁指数SYM-H, 同时考虑日地连线引力平衡点(L1点)太阳风地磁效应的滞后性, 精确分析了1998年10月18---19日大磁暴主相的行星际源. 分析结果表明, 磁暴主相的行星际源仅为行星际激波和行星际日冕物质抛射之间的太阳风(Sheath), 磁云对磁暴主相没有贡献. 这个磁暴事例的研究表明, 行星际磁场南向分量与太阳风动压的乘积是影响磁暴主相发展的关键参数.      

三维试验粒子轨道法在磁层粒子全球输运中的应用

根据磁层粒子动力学理论, 通过偶极磁场模型验证利用三维试验粒子轨道方法模拟近地球区(r < 8R_e)带电粒子运动特征的可靠性. 在此基础上, 以太阳风和磁层相互作用的全球MHD模拟结果为背景, 利用三维试验粒子轨道方法, 对非磁暴期间南向行星际磁场背景下太阳风离子注入磁层的情形进行数值模拟, 并对北向行星际磁场背景下太阳风离子注入极尖区以及内磁层的几种不同情形进行了单粒子模拟. 模拟结果反映了南向和北向行星际磁场离子向磁层的几种典型输入过程, 揭示出行星际磁场南向时太阳风粒子在磁层内密度分布的晨昏不对称性以及其在磁鞘和磁层内的大致分布, 并得出统计规律. 模拟结果与理论预测和观测结论相一致, 且通过数值模拟发现, 行星际磁场北向时靠近极尖区附近形成的非典型磁镜结构对于能量粒子经由极尖区注入环电流区域过程有重要的影响和作用.      

一个基于观测的太阳风背景模型

建立由太阳光球磁场和日冕偏振亮度等观测约束的单流体太阳风模型,包括日冕和太阳风的等离子体密度、速度和磁场,温度还有待于以后处理.这里采用高山观测台(HAO)MKⅢ的日冕偏振亮度(pB)在1.36Rs上的观测概图,根据Guhathakurta在1996年发展的日冕电子密度反演模型确定日冕的电子密度分布.同时采用Wilcox太阳观测台(WSO)的光球磁场视向分量的观测概图作为底部边界,根据Zhao等在1994年发展的水平电流-电流片(HCCS)模型得到全球磁场.Phillips在1995年及McComas在2003年分别用Ulysses第一次和第二次跨极飞行的观测发现,归一化到1 AU的太阳风动量流密度除了在10°～30°的纬度范围内略低以外几乎不变.根据这一结论,结合已经得到的密度数据,就可以得到日冕和太阳风的速度.将上面的模型应用于1918卡林顿自转周稳态太阳风的研究,结果与太阳活动极小期的观测基本相符,但是与观测相比较低速高密度区偏大,因此密度模型还有待改进.      

日冕三维太阳风模拟

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

背景太阳风数值模拟的熵守恒格式

背景太阳风研究是根据行星际扰动传播情况预测空间天气状况的基础,磁流体（MHD）模拟是背景太阳风研究的重要手段.采用一种新的数值计算方式,利用Ideal GLM-MHD将计算过程中产生的磁场散度以c_h的速度向计算区域外传播,从而消去磁场散度;重构部分使用受约束的最小二乘法,将磁场散度作为约束条件添加到重构中,进一步对重构后的磁场梯度进行修正;通量计算采用满足热力学第二定律的熵守恒格式,该格式能够确保在计算过程中熵不增,保证数值稳定.研究结果表明,该方法应用于太阳风数值模拟的求解得到了更加稳定的结果.      

嫦娥一号卫星太阳风离子探测器离子流量反演太阳风参数与初步结果分析

嫦娥一号卫星(Chang'E-1)上搭载的两台太阳风离子探测器(SWID-A/B)是国际上首次在200 km极月轨道观测等离子体环境的探测仪器.SWID-A/B的科学目标是探测月球附近等离子体与月球的相互作用,获得月球附近的太阳风速度、密度和温度.太阳风离子探测器的观测数据是各能量成分离子流量的直接反映,包含了太阳风离子的速度、密度和温度信息.本文设计了一种利用离子流量数据反演太阳风速度、密度和温度的算法,并通过模拟太阳风离子注入探测器的过程,验证了算法的可行性.对月球附近太阳风离子基本特征的分析研究表明,在太阳活动低年,空间环境扰动水平相对较低时,行星际太阳风运动到月球附近后依然保持着相同的变化趋势;太阳风离子的速度和密度与在上游行星际空间时相近;太阳风离子的温度则比在上游行星际空间时高103 K.     

月球尾迹的二维MHD模拟

采用理想的二维单流体MHD方程,对太阳风通过月球时所形成的尾迹结构进行数值模拟,得到了太阳风尾迹的粒子分布及磁场分布.模拟结果表明,在月球背阳面的本影区,太阳风粒子密度急剧下降,行星际磁场增强.当行星际磁场与太阳风流动方向平行时,尾迹被拖得很长,而磁场与太阳风流动方向垂直时,尾迹较短.      

Radio sounding of the solar wind acceleration region with spacecraft signals

Data from coronal radio-sounding experiments carried out on various interplanetary spacecraft are used to derive the empirical radial dependence of solar wind velocity and density at heliocentric distances from 3 to 60 solar radii for heliolatitudes below 60° and for low solar activity. The radial dependencies of solar wind power and acceleration are derived from these results. Summaries of the radial behavior of characteristic parameters of the solar wind turbulence (e.g., the spectral index and the inner and outer turbulence scales), as well as the fractional density fluctuation, are also presented. These radio-sounding results provide a benchmark for models of the solar wind in its acceleration region.     

Some characteristics of propagation of flare- or CME-associated interplanetary shock waves

Many interplanetary shock waves have a fast mode MHD wave Mach number between one and two and the ambient solar wind plasma and magnetic field are known to fluctuate. Therefore a weak, fast, MHD interplanetary shock wave propagating into a fluctuating solar wind region or into a solar wind stream will be expected to vary its strength.It is possible that an interplanetary shock wave, upon entering such a region will weaken its strength and degenerate into a fast-mode MHD wave. It is even possible that the shock may dissipate and disappear.A model for the propagation of a solar flare - or CME (Coronal Mass Ejections) - associated interplanetary shock wave is given. A physical mechanism is described to calculate the probability that a weak shock which enters a turbulent solar wind region will degenerate into a MHD wave. That is, the shock would disappear as an entropy-generate entity. This model also suggests that most interplanetary shock waves cannot propagate continuously with a smooth shock surface. It is suggested that the surface of an interplanetary shock will be highly distorted and that parts of the shock surface can degenerate into MHD waves or even disappear during its global propagation through interplanetary space. A few observations to support this model will be briefly described.Finally, this model of shock propagation also applies to corotating shocks. As corotating shocks propagate into fluctuating ambient solar wind regions, shocks may degenerate into waves or disappear.     

The accuracy of present models of the high altitude polar magnetosphere

The Polar satellite has explored the high-latitude, high-latitude magnetosphere out to 9 Earth radii (Re). The magnetic field data returned from this mission can be used both to provide data for new empirical models and to test existing models. Tests include comparing the observed location of the polar cusp with its position in the empirical models and comparing the strength of the magnetic field in the surrounding region. Near the cusp the magnetosphere is quite sensitive to solar wind conditions. In particular the energy density of the cusp plasma depends on the pressure of the solar wind applied to the interface of the cusp and the sheath. The applied pressure in turn depends on the shape of the magnetopause and the orientation of that interface, both controlled by the direction of the interplanetary magnetic field. Magnetohydrodynamic (MHD) models provide a coarse picture of the magnetosphere at high latitudes. While generally quite realistic, these too require testing against observations because even the MHD models must make some simplifying assumptions.     

A Brief Review of Interplanetary Investigations in China from 2014 to 2016

Great progress has been made in the research of solar corona and interplanetary physics by the Chinese scientists during the past two years (2014-2016). Nearly 100 papers were published in this area. In this report, we will give a brief review to these progresses. The investigations include:solar corona, solar wind and turbulence, superhalo electron and energetic particle in the inner heliosphere, solar flares and radio bursts, Coronal Mass Ejections (CMEs) and their interplanetary counterparts, Magnetohydrodynamic (MHD) numerical modeling, CME/shock arrival time prediction, magnetic reconnection, solar variability and its impact on climate. These achievements help us to better understand the evolution of solar activities, solar eruptions, their propagations in the heliosphere, and potential geoeffectiveness. They were achieved by the Chinese solar and space scientists independently or via international collaborations.      

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

A solar wind driven model of geosynchronous plasma moments

We describe a tabular specification model of the density and temperature of ions and electrons at geosynchronous orbit as a function of magnetic local time and solar wind parameters. This model can be used to provide boundary conditions for numerical ring current models. Unlike previous specification models of geosynchronous plasma moments, this model is parameterized by upstream solar wind conditions. We find that solar wind parameters are a better predictor of geosynchronous ion density than magnetospheric indices, and as upstream parameters they are often more appropriate as model inputs since they causally precede the model outputs. Of the upstream parameters that were tested, the best predictors of geosynchronous conditions were the solar wind flow pressure and the magnitude and Z-component of the interplanetary magnetic field.