地球弓激波的三维模拟

利用磁流体动力学(MHD)全球模拟结果,根据弓激波的跃变特性确定出弓激波位置,建立了一个新的综合考虑了快磁声马赫数、太阳风动压、行星际磁场强度以及磁层顶曲率半径的弓激波三维位型模型.将新模型与以往模型的模拟结果进行比较发现,新的弓激波全球模型结果可靠,解决了部分现有模型不能描述弓激波三维位型的问题.研究结果表明,在行星际磁场北向时,随着快磁声马赫数的增大,弓激波日下点距离减小,但是在行星际磁场南向时,快磁声马赫数的变化对弓激波日下点距离影响不大;弓激波位型在赤道面与子午面上存在明显的不对称性,而且随着行星际磁场的转向,这种非对称性也会发生相应改变;行星际磁场南向,B_z值较小时,子午面内弓激波位型已经不是简单的抛物线,出现了明显的类似于极尖区磁层顶的凹陷变化区.      

一种确定模拟的行星际MHD激波局地参数的新方法

简要阐述了分析模拟的行星际磁流体力学（MHD）激波的局部性质时，采用无厚度局部平面激波这一假设的合理性，说明了在激波未扰动区域（激波上游），物理量在几个小时内的变化很小这一事实，利用平面激波的分析方法，提出了分析模拟的行星际MHD激波的新方法，包括激波位置的确定，上下游状态参数的选择，激波局部参数的计算以及激波的分类，最后应用这种方法对一个二维的MHD模拟结果进行了分析。结果证实了过去文献关于磁流体力学混合激波空间连接和时间演化的链式规则，而且说明位于太阳赤道附近的慢激波和中间激波最终会发展为快激波。     

行星际空间MHD模拟激波跃变关系分析

分析了模拟得到的可以传播到1AU以远的日地空间磁流体力学激波与Rankine-Hugnoniot跃变关系的符合程度．通过对模拟激波的结构及其在传播过程中的演化进行的分析，提出了模拟激波的定位方法；基于所提出的定位方法，利用向长青提出的确定MHD激波局地参数的方法计算了模拟得到的激波与Rankine-Hugnoniot跃变关系的偏差．结果表明在激波传播到100Rs以后，激波对中前向快激波与Rankine-Hugnoniot关系的符合达到很高的程度，相对误差在10^-2数量级以内；并且在激波传播到150 Rs以后，相对误差在10^-3数量级以内．这个结果说明文中所使用的有限差分数值格式能较好地模拟激波．     

黄道面内行星际激波相互碰撞作用过程研究—能量效应

本文用两维半MHD数值模型,数值模拟研究了两邻近扰动源所产生的激波在行星际空间黄道面内不同能量时的相互碰撞过程。在内边界(18R_s)两扰动中心的间距取为36°。结果表明:两弱激波(速度在500km/s左右以下)不会产生汇合,而是各自独立地传播;两中等强度激波(速度在1000km/s左右)将发生汇合,但在IAU尚可分辨;两强激波(速度在2000km/s以上)则在1AU以内已发生汇合,汇合后形成一个新激波,其磁场结构与单激波明显不同。激波能量越大,两激波汇合的越快。     

2 AU以内的“渐进型”太阳高能粒子事件模拟

太阳高能粒子（Solar Energetic Particle,SEP）事件是影响地球空间以及深空辐射环境的主要因素之一。“渐进型”太阳高能粒子事件中的高能粒子主要来自于日冕物质抛射（Coronal Mass Ejection,CME）所驱动的激波扩散加速（Diffusive Shock Acceleration,DSA）过程。CME驱动的激波在行星际的传播过程中,其结构不断演化,进而影响到高能粒子的加速过程。本文利用二维太阳高能粒子加速和传播模型,对发生于2014年4月18日的太阳高能粒子事件实例进行了数值模拟。模型考察了黄道面上2 AU的距离以内包含地球所在位置的4个不同点,分别计算了每个点上高能粒子的通量。数值模拟的结果表明：黄道面内不同位置的观察点,与激波波前的磁力线连接不同,从而导致观察点处高能粒子的通量有着显著的差异。该模型的计算结果可以为深空探测计划开展辐射环境研究提供必要的输入。     

黄道面内行星际激波相互碰撞作用过程研究—间距效应

本文用两维半MHD数值模型,模拟研究了不同间距的两扰动源所产生的激波在行星际空间黄道面内的相互碰撞过程。在内边界(18R_s),两扰动中心的间距取为24°、36°、48°、54°、72°、120°,模拟了两激波速度均为1000km/s和2000km/s两种情形。结果表明:1.速度为1000km/s的两中等激波,相距24°时,在1AU以内强烈相互作用而汇合成一新激波;相距36°时,在1AU附近发生汇合,但尚可分辨;相距48°时,在1AU以内有相互作用,但可近似认为是独立传播的;相距72°时,将独立向前传播。2.速眨为2000km/s的两强激波,相距≤48°时,将强烈汇合成一新激波;相距54°时,在1AU附近产生汇合,但尚可分辨;相距72°时,相互作用不足以汇合成新的激波;相距120°时,其传播相互独立地进行。并给出了上述情形下两激波汇合的时间和位置,详细比较了汇合激波与相同能量的单激波的等离子体特性。     

太阳色球层-日冕中扰动的传播(Ⅱ)

首先讨论了磁流波传播的线性特征，然后构造了球坐标中-自治的非等温，非均匀等离子体初态，应用二维时变可压缩磁流体动力学模拟，数值研究了色球层底部压力脉冲所引起扰动的全球传播过程，结果表明，在极区，压力脉冲导致的扰动传播可以区分成两类不同模式的波动，快磁声波与慢磁声波，而在赤道附近，传播扰动是快模磁声波，在源区附近还存在一非传播的扰动，模拟结果的特征有助于解释SOHO/EIT观测到的波动事件。     

全球电离层发电机方程的一种数值解法

为了将建立在赤道为对称面假设基础上的二维半球域对称）电离层发电机理论扩展到全球,本文在解发电机方程时,去掉赤道对称面这一假设条件和两极附加的边界条件,使通常的边值问题变成一种“无边界’问题．它更合理地描述了电离层发电机过程,而且可以给出半球发电机理论所不能给出的特征,例如,可以对Ｓｑ电流体系的ＵＴ变化、季节变化、非偶极子地磁场部分的发电机效应等南北半球不对称的情况进行数值计算．计算实例表示,在简单的赤道对称情况下,它可以复现二维半球（对称）发电机模型的结果．在不对称情况下,电流体系的强度、焦点位置、电流方向等显示出较为复杂的变化,而半球发电机理论则无法解释这些特征．     

1998年5月2日日冕亮度观测图的数值研究

采用数值手段模拟了1998年5月2日日冕亮度观测图，计算模式改进为球坐标下系下特殊的二维理想磁流体（MHD）模型，即把（r,φ）坐标建立在SOHO观测日冕亮度的子午面上，消除了子午面极区的几何奇异性，根据SOHD日冕观测布置磁极子得到初始磁场位形，内边界条件采用自治的投影特征边界条件，计算迭代出稳态的多 磁场结构，得得到了与观测基本一致的亮度图，计算结果表明在太阳表面附近磁场位形对太阳凤等离子参数分布起控制作用。     

Helios飞船观测到的行星际日冕物质抛射事件

根据异常低的质子温度判据，从Heliosl和2飞船等离子体观测数据中(0．3-1AU)识别出大约160个行星际日冕物质抛射事件(ICME)，并在统计意义上分析了ICME在内日球空间的传播和演化规律．     

磁流体斜激波的碰撞

讨论了磁流体斜激波之间的碰撞及其与接触间断的相互作用规律,主要结论如下:(1)两个快激波碰撞后交换位置,同时出现一接触间断和一慢稀疏波对。(2)两个慢激波碰撞后交换位置且强度减弱,同时出现一接触间断和一块激波对。(3)一前向快激波与一后向慢激波碰撞后交换位置,快激波强度增加,慢激波强度减弱,同时出现一后向快激波、一负接触间断和一前向慢稀疏波。(4)一前向快激波与一正(负)接触间断相互作用后交换位置,快激波减弱,同时出现一后向快稀疏波(快激波)、一后向慢激波和一前向慢激波(慢稀疏波).(5)一前向慢激波与一正(负)接触间断相互作用后交换位置,慢激波减弱,同时出现一后向慢稀疏波(慢激波)和一快稀疏波(快激波)对。      

磁流体斜激波的汇合

本文讨论磁流体快、慢激波的汇合作用规律,主要结论如下:(1)两个前向快激波汇合之后,形成一更强的前向快激波,尾随一前向慢稀疏波、一正接触间断(后侧密度大于前侧)、一后向慢激波和一后向快稀疏浚。(2)两个前向慢激波汇合之后,形成一更强的前向慢激波,尾随一正接触间断、一后向慢稀疏波和一后向快激波2在前向慢激波前方出现一前向快波,它或为稀疏浚(中、小激波角情况),或为激波(大激波角情况).(3)前向快激波会追上前向慢激波而发生汇合,之后互换位置且强度减弱,尾随一正接触间断和一后向稀疏波对。      

磁流体中间激波的相互作用

本文讨论磁流体中间激波的相互作用规律．主要结论：中间激波汇合的产物为后向快简单波、后向慢简单波或慢激波、接触间断、前向慢激波和前向快激波，其中后向波成份和接触间断很弱．当左（右）激波较强时，中国激波碰撞产物为后（前）向快激波、后（前）向慢简单波或慢激波、负（正）切向间断、前（后）向慢简单波和前（后）向快激波．     

耀斑激波传播的数值研究

采用一维磁流体力学模型和激波装配法,分析耀斑激波在行星际空间的传播特性和激波下游区的波动结构,并就激波装配法和激波捕捉法的模拟结果的精度和可靠性进行比较.      

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.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Cone model-based SEP event calculations for applications to multipoint observations

The problem of modeling solar energetic particle (SEP) events is important to both space weather research and forecasting, and yet it has seen relatively little progress. Most important SEP events are associated with coronal mass ejections (CMEs) that drive coronal and interplanetary shocks. These shocks can continuously produce accelerated particles from the ambient medium to well beyond 1 AU. This paper describes an effort to model real SEP events using a Center for Integrated Space weather Modeling (CISM) MHD solar wind simulation including a cone model of CMEs to initiate the related shocks. In addition to providing observation-inspired shock geometry and characteristics, this MHD simulation describes the time-dependent observer field line connections to the shock source. As a first approximation, we assume a shock jump-parameterized source strength and spectrum, and that scatter-free transport occurs outside of the shock source, thus emphasizing the role the shock evolution plays in determining the modeled SEP event profile. Three halo CME events on May 12, 1997, November 4, 1997 and December 13, 2006 are used to test the modeling approach. While challenges arise in the identification and characterization of the shocks in the MHD model results, this approach illustrates the importance to SEP event modeling of globally simulating the underlying heliospheric event. The results also suggest the potential utility of such a model for forcasting and for interpretation of separated multipoint measurements such as those expected from the STEREO mission.     

太阳耀斑行星际激波传播中的追赶效应

本文采用二维MHD模型对具有不同间隔时间的2个耀斑先后爆发,模拟研究它们所对应的行星际激波间的追赶效应,并和单个耀斑所产生的行星际激波相比较。研究结果表明,间隔时间一天以内的2个耀斑激波在行星际空间向外传播时,激波之间有明显的相互作用发生,间隔时间的长短决定了激波传播过程中追赶效应的强弱。根据数值试验结果,追赶效应可归纳为4类,(1)强追赶效应,(2)中等追赶效应,(3)弱追赶效应,(4)无追赶效应。属于强追赶效应的2个耀斑激波传播至1AU处,产生的行星际扰动非常相似于单个耀斑激波的扰动。