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

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

无碰撞激波的理想MHD模拟

采用了WENO格式数值求解一维理想MHD方程组，模拟了行星际无碰撞激波，研究了垂直无碰撞激波与行星际反向磁场结构和高密度等离子体团的相互作用过程，并与粒子模拟的结果进行比对，两者的结果非常类似．模拟结果表明，对涉及无碰撞激波的大部分现象中理想MHD模拟是准确且可行的，同时相对于粒子模拟又有很好的计算效率，便于扩展至二维或三维的情形．     

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

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

地球磁鞘中磁场起伏特性的分布

本文利用行星际起伏通过激波后的变化的MHD模型, 具体讨论了地球磁鞘中磁场起伏特性在黄道面内的分布.主要结果是:(1)行星际磁场起伏的强度和各向异性在磁鞘中被显著放大;(2)行星际磁场基本位于黄道面内时, 磁鞘中磁场起伏特性(强度、相对起伏和各向异性等)呈现明显的晨一昏不对称性, 早晨侧(准平行激波)显著地高于黄昏侧(准垂直激波);(3)行星际磁场方向对磁场起伏特性在磁鞘中的分布有强烈的控制作用, 早晨侧响应灵敏, 黄昏侧反响不大.相对地讲, 黄昏侧的磁活动较之早晨侧稳定;(4)行星际磁场转南的增强将导致磁鞘中磁场起伏的最大区域自黄道面低纬向北极高纬移动, 南-北不对称性磁活动随之加强, 最强大致出现在磁场与黄道面相交成大约45°时, 而晨-昏不对称性的强弱程度则发生相反变化;(5)行星际磁场的相对起伏增加, 晨-昏不对称性反随之减弱.磁鞘中磁场起伏分布的特性与卫星观测大体符合, 是磁顶、边界层某些晨-管不对称性出现的可能起因之一.      

火星空间磁场结构特征

在火星空间模拟的单流体MHD模型的基础上, 研究了火星空间磁场结构及火星表面局部磁异常对磁场结构的影响. 在太阳风与火星相互作用的过程中, 形成弓激波和磁堆积区, 行星际磁场弯曲并向两极移动且被拖拽变形, 大部分磁力线从火星两极绕过, 通过火星之后在磁尾留下V字形结构. 火星表面附近局部磁异常也对火星磁场结构产生不可忽视的影响. 不同位置和强度的磁异常与太阳风相互作用形成结构及形态各异的微磁层, 如被拖拽的微磁层和存在开磁力线的微磁层等. 局部磁异常改变了近火磁场结构, 并可能改变等离子体的分布.      

电离层等效电流体系对行星际激波的响应研究

行星际激波是导致地球磁层-电离层系统发生扰动的重要原因之一,其可以通过对磁层-电离层系统电流体系的改变来影响地磁变化.本文采用全球三维磁流体力学数值模拟方法,分析了行星际激波作用下电离层等效电流体系的即时响应.模拟结果表明,在激波作用下伴随着异常场向电流对的产生,电离层在午前午后出现一对反向的等效电流涡.这对涡旋一边向极侧和夜侧运动,一边经历强度增强和减弱直至消失的过程.激波过后等效电流体系图像逐渐演化为激波下游行星际条件控制的典型图像.这个响应过程与行星际激波强度有关,激波强度越强,则反向的等效电流涡旋强度越大,寿命也就越短.      

地球弓激波的三维模拟

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

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

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

行星际起伏场通过激波后的变化

本文采用磁流体平面斜激波模型探讨了行星际起伏场(包括磁场和等离子体速度)通过激波后的变化。得到的近似解析解表明:由本模型所得到的结果能较好地解释人们[1,2]对行星际、磁鞘层观测资料进行统计研究所得到的认识。      

磁流体力学激波的特征速度

根据磁流体力学间断面的守恒条件与磁流体力学单波方程的相似性,引入了一个称为磁流体力学激波特征速度的物理量,它是激波在波前后两侧介质中传播速度的几何平均值,当激波很弱时,它趋近于磁流体力学波的速度。本文导出一组以此特征速度为强度参数的激波跃变公式,形式上与单波的公式组非常相像,从而简化了激波跃变量的计算。其中密度跃变公式本身解析地证明:磁流体力学激波与磁流体力学波传播速度之间的关系是由激波是压缩波这一特性直接决定的。      

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.     

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黎曼问题的求解对分析行星际扰动演化趋势和激波相互作用具有重要意义。本文基于MHD共面间断可以分解为前后向快慢激波、中心简单波和接触间断的假定,提出MHD共面黎曼问题的一种三参数迭代解法,运用该解法实现纯法向速度间断的分解,并初步探讨该解法的适用范围。      

行星际激波的多方指数

从Ｈｅｌｉｏｓ等６个飞船１９７２年至１９８７年的太阳风观测资料找出９８６个行星际快激波，利用ＭＨＤ Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ方程组分析了激波上下游参数关系，确定出每一个激波的最佳多方指数。结果指出，多方指数有一分布，其最可几值约在１．７０－１．７５之间。     

Propagation of interplanetary shocks into the Earth’s magnetosphere

This paper is devoted to the study of propagation of disturbances caused by interplanetary shocks (IPS) through the Earth’s magnetosphere. Using simultaneous observations of various fast forward shocks by different satellites in the solar wind, magnetosheath and magnetosphere from 1995 till 2002, we traced the interplanetary shocks into the Earth’s magnetosphere, we calculated the velocity of their propagation into the Earth’s magnetosphere and analyzed fronts of the disturbances. From the onset of disturbances at different satellites in the magnetosphere we obtained speed values ranging from 500 to 1300 km/s in the direction along the IP shock normal, that is in a general agreement with results of previous numerical MHD simulations. The paper discusses in detail a sequence of two events on November 9th, 2002. For the two cases we estimated the propagation speed of the IP shock caused disturbance between the dayside and nightside magnetosphere to be 590 km/s and 714–741 km/s, respectively. We partially attributed this increase to higher Alfven speed in the outer magnetosphere due to the compression of the magnetosphere as a consequence of the first event, and partially to the faster and stronger driving interplanetary shock. High-time resolution GOES magnetic field data revealed a complex structure of the compressional wave fronts at the dayside geosynchronous orbit during these events, with initial very steep parts (10 s). We discuss a few possible mechanisms of such steep front formation in the paper.     

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

行星际慢激波的传播和演化

采用二维理想MHD模型，分别在日球赤道面(二维二分量模型)和日球子午面(二维三分量模型)内研究太阳风中慢激波的传播和演化规律．结果表明，慢激波在向外传播的过程中逐渐演化为由原慢激波和新产生的快激波构成的激波系统，该激波系统在子午面内相对慢激波源中心法线基本对称，而在赤道面内则是不对称的：快激波阵面和慢激波阵面之间存在一个切触点，该处两个激波合并，蜕化为气体激波．上述切触点相对激波源中心法线东偏，且东偏角度在激波系统向外传播过程中不断增加．初步分析表明，行星际磁场的螺旋结构是产生日球赤道面内慢激波传播和演化的东西不对称性的主要原因．