磁鞘快速流引起的磁层顶凹陷事件

最近研究表叽磁层顶凹陷对磁层-电离层耦合具有重要作用.但是,磁层顶凹陷现象的确认需要多颗卫星的联合观测,目前为止报道的磁层顶凹陷事例非常少.本文利用THEMIS 5颗卫星的联合观测结果,分析了一例由磁鞘快速流引起的磁层顶凹陷事件.2007年7月21日10:00UT-10:45UT,位于日下点磁层顶附近的THEMIS卫星...     

IMF北向时太阳风粒子向磁层输运的试验粒子模拟研究

以ACE卫星实时观测数据驱动的全球磁流体模拟为背景场,选取2003年10月22－24日行星际磁场（IMF）持续北向的事件,使用试验粒子方法,对太阳风粒子向磁层输运的过程进行模拟研究,分析北向IMF下太阳风粒子注入磁层过程中粒子在磁层内的空间分布和时间演化特征。IMF北向期间,进入环电流区域的粒子在晨侧区域的密度大于昏侧,且晨侧的粒子分布范围更广。背阳面磁鞘中的太阳风粒子可以通过低纬边界层进入磁层,但很难通过南北侧磁层顶进入磁层。进入磁尾的太阳风粒子聚集形成冷而密的等离子体片（CDPS）,模拟中CDPS的空间分布和密度大小与观测数据符合。在IMF长时间北向期间,磁尾的粒子数量呈现随时间增长的趋势,并存在约20 min的小幅度准周期变化和约5~6 h的较大幅度的准周期变化。     

关于IMF北向分量特强时的太阳风－磁层－电离层耦合

Relationship between the solar wind-magnetosphere energy coupling function and the energy consumption in the magnetospheric system is tested by the samples during period of northward B_z(>10 nT) of IMF. It is shown that the energy function which can well describe the solar wind-magnetosphere coupling when B_z>0 is no more applicable for B_z>10 nT, but the coupling still exists, the effects of which are mainly found at the polar cap. Dst often turns to positive when B_z> 10 nTand the coupling mechanism is different from that when B_z<0. New coupling function has to be introduced to describe the solar wind-magnetosphere coupling for the case of B_z> 10 nT.      

关于IMF北向分量特强时的太阳风－磁层－电离层耦合

关于ＩＭＦ北向分量特强时的太阳风－磁层－电离层耦合沈长寿，资民筠（北京大学地球物理系，１００８７１）（中国艺术研究院，北京）关键词行星际磁场，太阳风－磁层能量耦合函数为检验Ａｋａｓｏｆｕ提出的能量融合函数“’在ＢＺ强北向时能否反映太阳风与磁层扰动间的...     

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

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

基于能量色散特征计算磁层顶重联位置的方法研究

磁层顶磁场重联是太阳风向磁层输入能量的主要方式.重联如何触发一直是空间物理研究的难点,其机制仍然有待深入研究.由于卫星穿越磁层顶时,很难恰好穿越重联发生的区域,因此难以观测到重联的触发条件.本文利用THEMIS卫星观测,确立了反演磁层顶重联点的方法.当重联刚开始发生时,卫星能够观测到离子的能量色散特征,可利用其计算卫星到重联发生位置的距离.沿着磁力线模型追踪该距离即可反演出磁层顶发生重联的位置.与其他方法进行了对比分析,结果显示本文方法比其他方法具有更高的精度.     

2.5维自适应磁场重联MHD模式

磁雷诺数(R_m)是影响磁场重联的重要因素. 真实的物理环境中R_m往往很高, 例如, 在行星际空间和太阳日冕中R_m通常大于10⁴量级. 高R_m条件下的磁重联表现出很多异常特性, 然而高R_m条件下的磁场重联数值模拟需要很高的时空分辨率, 否则很难分辨出重联过程中形成的薄电流片. 本文基于自适应软件包PARAMESH将并行自适应网格技术引入磁场重联数值模拟, 建立了一个2.5维自适应磁场重联MHD模式, 研究高磁雷诺数条件下重联的动态演化过程, 进而将不同磁雷诺数的参数进行对比研究. 结果表明, 该模式可以自动捕捉到磁场重联产生的奇性电流片, 高磁雷诺数条件下产生的慢激波结构可提供一种快速磁能释放机制.     

行星际磁场与磁尾磁通量绳形成的关系

收集了Cluster卫星在2001-2005年间观测到的磁尾磁通量绳事件,并对磁通量绳(magnetic flux rope)形成及其内部磁场结构与行星际磁场(IMF)的关系作了统计研究.考虑磁通量绳被观测到时行星际磁场的条件,在所有73个磁通量绳事件中,行星际磁场By分量占有主导地位的事件有80%,且78%的事件具有与行星际磁场By分量相同方向的核心场.行星际磁场通过在磁层顶与地球磁场相互作用改变南北等离子体片内磁场相对方向,形成有利于磁通量绳形成的磁场位形,并且行星际磁场By分量的方向对磁通量绳内部核心场的方向具有决定性影响.从统计结果来看,磁通量绳的形成并不会依赖于行星际磁场Bz分量的方向.     

行星际激波对地球磁层的压缩效应分析

2004 年11月9日WIND飞船探测到一个典型的行星际激波. 激波前行星际磁场为持续约50 min的弱南向磁场, 越过激波面, 磁场发生北向偏转且太阳风动压脉冲增强. 在此强动压脉冲增强结构作用下, 磁层被压缩至一个很小的区域. 激波作用于磁层时引起地球同步轨道 各区域高能粒子通量的响应, 但是不同磁地方时的高能粒子通量的响应不同, 表现出双模式扰动, 即在晨昏两侧各能段的电子和质子通量显著增强, 在子夜侧发生类似于亚暴的无色散粒子注入现象. 扰动从向阳面传输到背阳面, 向阳面粒子通量最先增强, 随后背阳面靠近晨昏两侧, 粒子通量开始增强, 最后子夜侧粒子通量表现出无色散高能粒子注入的特点. 另外, 在靠近正午侧, 质子通量先于电子通量发生响应, 在子夜侧电子通量则先于质子通量发生响应. 利用位于向阳面正午两侧的GOES-10 和 GOES-12卫星观测数据发现, 激波作用于磁层时靠近晨侧的磁场变化表现出简单压缩效应, 而靠近昏侧的磁场变化则显然不同, Bx分量减弱, Bz分量几乎减为零, 而By分量则显著增强. 此外, 位于近地磁尾低纬尾瓣区的TC-1卫星观测到激波触发的尾瓣SI现象.     

空间等离子体压力各向异性对磁场重联的影响

基于二维时变可压缩磁流体动力学模拟,数值研究了等离子体压力各向异性对磁场重联的影响,发现一个小的压力各向异性(P_⊥=1.02P_//)即可大大加速磁场重联的发展,这可能是由于磁镜不稳定性与撕裂模不稳定性共同起作用.在P_⊥P_//的情况下,撕裂模受到抑制,电流片中不能形成大型磁岛.     

各向异性等离子体中磁场重联的混合模拟研究

应用二维三分量混合模拟方法数值研究了各向异性等离子体中的磁场重联过程.计算结果表明,当等离子体垂直于磁场方向的压强大于平行方向的压强时（P_⊥/P_∥=1.5）,等离子体不稳定性的增长率会大大增强,重联速度也会加快;当等离子体垂直方向的压强小于平行方向的压强时（P_⊥/P_∥=0.6）,会出现火蛇管不稳定性,将抑制撕裂模不稳定性和磁场重联过程.     

磁场重联的二维粒子模拟研究

利用二维全粒子模拟方法研究了无碰撞等离子体中的磁场重联过程,得到了不同区域的离子和电子速度分布。计算结果表明,电子和离子在扩散区中的不同动力学特性产生的Hall电流使磁场的y分量B_y呈现四极形分布。离子和电子的速度分布偏离了初态时的Maxwell分布,呈现非局域的多重分布。同时由于磁场重联而产生的电场使电子在X点附近得到加速和加热,因而在电子的能谱分布中形成一高能尾。     

三维非线性撕裂模耦合和快速磁场重联

基于三维不可压缩电阻性ＭＨＤ方程，在长柱形位型下，数值研究了电阻撕裂模不稳定性所引起的磁场重联过程。研究结果表明，撕裂模的非线性相互作用和耦合，将产生许多模式的强裂解稳和导致快速的磁场重联。这一过程的进一步发展，导致具有螺旋形结构，不同模式的磁岛磁力线重叠，从而引起较大区域内的磁场各态经历和磁力线随机走向，形成撕裂模端动。     

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

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

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

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

日冕冲浪形成的磁流体动力学模拟

应用二维时变可压缩磁流体动力学模拟,数值研究了双极-单极磁场中电阻撕裂模不稳定性引起的磁场重联过程,用于模拟日冕冲浪的形成.结果表明,在包含有三区——双极场、电流片和单极场的磁静力平衡初态下,双极场和单极场中的磁力线将会直接重联,磁场演变成鞭状(whip)结构.由弯曲磁力线支撑的等离子体团向上运动到最高位置后,逐渐下落和弥散.等离子体团上升速度可达到0.10v_A(v_A为双极场中的Alfv'én速度).模拟结果证实日冕冲浪的形成可能与双极-单极场中的磁场重联密切相关.      

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

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

A CLASS OF TWO-STEP TVD MACCORMACK TYPE NUMERICAL SCHEME FOR MHD EQUATIONS

In this paper, a new numerical scheme of Total Variation Diminishing （TVD） Mac-Cormack type for MagnetoHydroDynamic （MHD） equations is proposed by taking into account of the characteristics such as convergence, stability, resolution. This new scheme is established by solving the MHD equations with a TVD modified MacCormack scheme for the purpose of developing a scheme of quick convergence as well as of TVD property. To show the validation, simplicity and practicability of the scheme for modelling MHD problems, a self-similar Cauchy problem with the discontinuous initial data consisting of constant states, and the collision of two fast MHD shocks, and two-dimensional Orszag and Tang’s MHD vortex problem are discussed with the numerical results conforming to the existing results obtained by the Roe type TVD, the high-order Godunov scheme, and Weighted Essentially Non-Oscillatory （WENO） scheme. The numerical tests show that this two-step TVD MacCormack numerical scheme for MHD system is of robust operation