球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

日冕三维太阳风模拟

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

Effect of current sheets on the power spectrum of the solar wind magnetic field using a cell model

A puzzling observation of solar wind MHD turbulence is the often seen Kolmogorov scaling of k^-5/3$k^{-5/3}$, even though the solar wind MHD turbulence is dominated by Alfvénic fluctuations. Recently Li et al. (2011) proposed that the presence of current sheets may be the cause of the Kolmogorov scaling. Here, using a cell model of the solar wind we examine the effect of current sheets on the power spectrum of the solar wind magnetic field. We model the solar wind as multiple cells separated by current sheets. We prescribe the spectra of turbulent magnetic field in individual cells as IK-like and examine the spectra along trajectories that cross multiple boundaries. We find that these spectra become softer and are consistent with the Kolmogorov-scaling. Our finding supports our recent proposal of Li et al. (2011).     

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

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

基于CPU/GPU的日冕偏振亮度并行计算模型

三维磁流体力学(MHD)数值模拟是用来研究日冕和太阳风最常用的方法之一, 其中将计算得到的日冕电子数密度转化为日冕偏振亮度(Polarization Brightness, PB)是与观测对比的重要方法. 由于待转换电子数据网格密度、PB数据网格密度和计算模型的复杂度, 使得日冕偏振亮度的计算比较耗时, 利用单CPU计算无法达到近实时转换日冕偏振亮度的要求, 从而影响了数值模拟的验证效率. 本文在CPU/GPU环境下, 利用CUDA编程技术, 提出了一个日冕偏振亮度并行计算模型. 实验结果表明, 该模型比CPU上的串行模型计算速度提高了31.86倍, 达到了近实时模拟与观测数据比对的计算要求.      

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

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

The effects of the solar magnetic polarity and the solar wind velocity on Bz-component of the interplanetary magnetic field

We have studied the effect of both solar magnetic polarity and the solar wind velocity on the Bz-component of the interplanetary magnetic field, IMFBz, for the minimum activity of the solar cycles 21, 22, 23 and 24. We made a statistical study of IMFBz in the first section which is considered as an extension of Lyatsky et al. (2003). They made a statistical study of IMFBz for two periods of minimum solar activity 22 and 23 related to 1985–1987 and 1995–1997 when the solar magnetic field had opposite polarity. Our results seem to be consistent with the results obtained by Lyatsky et al. (2003). We found that there is a dependence of IMFBz on the IMFBx and the solar magnetic polarity for the minimum periods of the selected four solar cycles. In addition, we found that there is a dependence of IMFBz on the solar wind velocity.     

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.     

How to address the accuracy of empirical magnetic field models?

Empirical magnetic field models are compared with high-altitude magnetic field measurements and results from an MHD simulation. Comparison of the T96 model and observations from GOES-8 and GOES-9 shows that if the observed solar wind and IMF parameters are used to compute the model field, the model field is more stretched than the observed field. On the other hand, if measurements made by one spacecraft are used to find the model parameters that give a best-fit field at that location, the RMS error can be reduced also at the other spacecraft four hours away in local time. Comparison of T96 model and MHD simulation results shows that the empirical models have a thinner current sheet than the MHD simulation, but that the lobe field values are quite similar to each other. Furthermore, a comparison of an event-oriented, modified T89 model and MHD simulation by Pulkkinen et al. [2000] reveals that if the empirical model is constructed by fitting to in-situ measurements, the resulting model is very similar to the MHD simulation magnetic field. These results indicate that an efficient method of utilizing the present-day empirical models is to select model parameters based on measurements from a few individual points.     

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

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

三维MHD背景场下太阳高能粒子平均自由程模拟

太阳高能粒子（SEP）的平均自由程是研究SEP传播的重要参数，由SEP的物理性质和太阳风物理性质决定.使用MHD-SEP模型对三维MHD背景场下的平均自由程进行了探讨，利用该模型具有可提供接近物理真实的太阳风背景场的优势，对SEP的平均自由程进行了定性分析.分别对太阳活动高年和低年选取2个卡林顿周进行模拟，定性分析其空间变化，并研究平均自由程与径向太阳风速度的相关性.结果表明:该方法得到的平均自由程空间分布与以往研究得到的关于平均自由程的结论相吻合，可以用来定性确立平行平均自由程；该模型可以反映不同事件中平行平均自由程分布的不同特征；表现了平均自由程与径向太阳风速度有很好的负相关关系.结果可为未来缓变SEP平均自由程研究作参考.      

Invariant Modulation of IMF Clock Angle on the Solar Wind Energy Input into the Magnetosphere

By use of the global PPMLR Magnetohydrodynamics (MHD) model, a serial of quasisteady- state numerical simulations were conducted to examine the modulation property of the interplanetary magnetic field clock angle θ on the solar wind energy input into the magnetosphere. All the simulations can be divided into seven groups according to different criteria of solar wind conditions. For each group, 37 numerical examples are analyzed, with the clock angle varying from 0° to 360° with an interval of 10°, keeping the other solar wind parameters (such as the solar wind number density, velocity, and the magnetic field magnitude) unchanged. As expected, the solar wind energy input into the magnetosphere is modulated by the IMF clock angle. The axisymmetrical bell-shaped curve peaks at the clock angle of 180°. However, the modulation effect remains invariant with varying other solar wind conditions. The function form of such an invariant modulation is found to be sin(θ/2)2.70 + 0.25.      

Numeric simulations of in situ observations of interacting ejecta near 1 AU

We use a simple numerical model (González-Esparza, J.A., Santillán, A., Ferrer, J. A numerical study of the interaction between two ejecta in the interplanetary medium: one and two dimensional hydrodynamic simulations, Ann. Geophys. 22, 3741–3749, 2004) to study the evolution of three events in the solar wind reported by Wang et al. (Wang, Y.M., Ye, P.Z., Wang, S. Multiple magnetic clouds: several examples during March–April 2001. J. Geophys. Res. 108, 1370, 2003, doi:10.1029/2003JA009850), where two interacting ejecta detected in situ by ACE near 1 AU were related to CMEs observed previously by SOHO-LASCO. The study is based on a 1-D hydrodynamic model using the ZEUS code (Stone, J.M., Norman, M. ZEUS 2-D: A radiation magnetohydrodynamics code for astrophysical flows in two dimensions, I, the hydrodynamics algorithms and tests, Astrophys. J. 80, 753, 1992). Although this model cannot address either the magnetic field dynamics or the complex geometrical effects intrinsic in the three-dimensional nature of the phenomena, it illuminates the transferring of momentum and evolution of interacting large-scale solar wind disturbances in those cases where there is no merging (magnetic reconnection) between the two ejecta. This model can reproduce, in some cases, characteristics of the events such as transit times and flow signatures as inferred from the two-point measurements from spacecraft.     

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

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

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.     

较小尺度磁螺旋线管相互作用的数值研究

在子午面内,偶极子场和六极子场适当叠加得到势磁场,势磁场与太阳风长时间相瓦作用得到特殊的冕流背景结构.在这种背景结构下,两个较小尺度的磁螺旋线管模型能够连续浮入到计算域,在计算域内相互作用,触发了日冕物质抛射(CME).在数值模拟这一过程时,较小尺度的磁螺旋线管模型具有同心圆形磁场结构,模型中心等离子压强与边界压强之比m=2,模型的半径分别取为a=0.07 R.和a=0.1 Rs(Rs为太阳半径).在这两种情况下,得到了两种典型的计算结果.当a=0.07 Rs时,两个磁螺旋线管模型相瓦作用,在7 Rs内融合成一个磁螺旋线管模型,向外传播;当a=0.1 Rs时,两个磁螺旋线管模型相互作用,作为一个整体向外传播,在计算域内没有融合到一起,基本上保持各自的磁场结构.      

Interplanetary Physics in Mainland China

During the past two years (2016-2018), great achievements have been made in the Chinese research of interplanetary physics, with nearly 100 papers published in the academic journals. The achievements are including but not limited to the following topics:solar corona; solar wind and turbulence; filament/prominence and jets; solar flare; radio bursts; particle acceleration at coronal shocks; magnetic flux ropes; instability; instrument; Coronal Mass Ejections (CMEs) and their interplanetary counterparts; Magnetohydrodynamic (MHD) numerical modeling; solar energetic particles and cosmic rays. The progress further improves our understanding of the eruptions of solar activities, their evolutions and propagations in the heliosphere, and final geoeffects on our Earth. These results were achieved by the Chinese solar and space scientists independently or via international collaborations. This paper will give a brief review of these achievements.      

太阳风向磁层输入电磁能量研究

利用全球磁流体力学（MHD）模拟结果,通过确立包含磁层顶的太阳风流线内边界来识别三维磁层顶位形,并以极尖区位置作为磁层顶日侧与夜侧的分界线,在此基础上定量研究了不同条件下穿过磁层顶向磁层内输入的电磁能量. 研究发现,磁层顶的能量传输与太阳风条件密切相关,磁重联是控制电磁能量传输的重要机制. 结果表明,当IMF（行星际磁场）南向时,极尖区后方的磁尾附近存在电磁能输入最大值,当IMF北向时,电磁能输入最大值发生在极尖区附近;南向IMF条件下,在IMF强度增大或太阳风密度增大时,磁层顶电磁能传输的电磁能量比北向IMF条件时增加更显著. 太阳风通过调节磁层顶面积间接影响到磁层顶能量传输大小. 研究还发现,北向IMF与南向IMF条件下穿过磁层顶的电磁能输入的比值范围约为10%～30%,此比值一定程度上反映了北、南方向IMF与地磁场磁重联效率的比值.      

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.