引力场中磁场重联的数值研究──（Ⅰ）磁岛的形成和合并

数值研究了引力场中电阻撕裂模不稳定性所引起的磁场重联，结果表明，在电流片长度Ｌ＝１Ｈ、半宽度δ＝Ｈ／２４的情况下，电流片两端附近将出现两个Ｘ线的磁场重联，并形成磁岛和高温高密度的等离子体团．磁岛宽度随着时间而增长，在ｔ≈５５τＡ时达到饱和，最大饱和岛宽约为４δ．同时从ｔ≈３７τＡ开始，磁岛中心位置逐渐下降；在ｔ＝５７τＡ时发生结合不稳定性，磁岛与底部附近的闭合磁场区合并，导致磁场湮灭和磁能的快速释放．另一方面，由于引力场和等离子体非均匀性的影响，顶部附近随时间而增长的等离子体外流速度达１．１４ＶＡ∞；磁岛中等离子体向下运动的最大速度达１．４１ＶＡ∞，且大于局地声速；在磁岛前方可形成快激波．这些结果可用于解释双带耀斑中后随耀斑环的形成、磁场湮灭和磁能释放、以及Ｄｏｐｐｌｅｒ速度图上观测到的红移现象．     

引力场中磁场重联的数值研究（Ⅱ）磁力线足点的剪切运动

本文基于二维三分量可压缩磁流体动力学模拟，数值研究由于磁力线足点在光球层的剪切运动引起日冕电流片中的磁场重联过程。结果表明，磁力线足点的剪切运动作为引起强迫磁场重联的一种触发机制，将加速磁场重联的发展和磁岛的合并过程。结合不稳定性导致等离子体急剧加速，在β_∞＝０．１的情况下其加速度达到０．３４ν_Ａ∞／τＡ，等离子体的最大下落速度可达１．９０ν_Ａ∞，大于纯电阻撕裂模情况。还讨论了β_∞值对这种磁场重联过程的影响。β_∞值越小，磁场重联和磁岛合并过程发展得越快。     

引力场中磁场重联的数值研究（Ⅲ）等离子体团的形成和运动

应用解析和数值相结合的方法求得了含有闭场区、中性片和开场区的二维局域磁静力平衡日冕基态，进而在此基态下数值研究了由电阻撕裂模不稳定性引起的磁场重联过程．结果表明，在电流片长度远大于其半宽度的情况下，仍将发生具有两个Ｘ线的磁场重联，形成磁岛和高温高密度的等离子体团．等离子体团运动的方向取决于闭合区底部等离子体压力和磁压的比值β_０．当β_０较大时，等离子体团向下运动，引起结合不稳定性和磁岛的合并，等离子体团中的等离子体不断落入闭合磁场区两侧．当β_０较小时，等离子体团向上运动，导致电流片中等离子体的喷发．结果还表明，磁张力的分布是决定等离子体团运动方向的主要因子．     

非稳态多重X线重联的类型与特征

为了进一步认识间歇性多重Ｘ线重联的特征，本文分析计算结果，发现入流自边界进入后向两侧发散，沿着ｘ轴（ｚ＝０）和出流边界（ｚ＝２），压强Ｐ和磁场Ｂ的分布形态表明这种非稳态重联属快模式扩张型．数值结果还表明，进入扩散区的入流马赫数Ｍ_ｉ与磁雷诺数Ｒｍ之间基本满足关系式Ｍ_ｉ≈１．５Ｒ_ｍ^－１／２─１．７５Ｒ_ｍ^－１／２，而且等离子体的加速与压强Ｐ_ｔ(Ｐ_ｔ     

地球远磁尾中的磁场重联（Ⅰ）流动撕裂模不稳定性

本文数值研究了地球远磁尾中流动撕裂模不稳定性所引起的磁场重联过程.结果表明,在短暂的指数增长之后,当磁岛宽度接近等离子体片厚度时,流动撕裂模不稳定性的增长率大大降低,最终磁岛宽度趋于饱和.磁岛的饱和宽度随着磁Reynolds数S的增大而减小,随着剪切流动层宽度δV的增大而增加.在S→∞的情况下,流动撕裂模将退化为流动颈缩模,磁场重联不再发生.当飞船通过由流动撕裂模不稳定性所形成的磁岛时,即可观测到磁场B_z分量由北转南或由南转北的现象.     

多层电流片中撕裂模不稳定性的数值研究

应用二维磁流体动力学模拟方法数值研究了三层电流片中电阻撕裂模不稳定性的特征及磁场重联过程．结果表明，这是一种复杂的非稳态磁场重联．在初期阶段，三个电流片中分别由撕裂模不稳定性引起磁场重联，形成薄而长的磁岛．随着撕裂模不稳定性的非线性发展，每个磁岛的宽度都逐步增大，以至导致新的磁场重联发生．同时，三个电流片的强度都逐渐减弱，且原中心反向电流区最终消失．部分磁能不断地转化为等离子体的热能和动能，引起等离子体的加热和加速．多层电流片中撕裂模不稳定性引起的自发重联，可能对太阳耀斑、日冕加热、太阳风与磁层耦合等有重要影响．     

引力场中磁场重联的数值研究(IV)多层电流片

应用数值方法研究了日冕多层电流片中电阻撕裂模不稳定性的非线性演化和磁场重联过程,结果表明,计算区域顶部附近两侧电流片中的磁岛和等离子体团向上抛射,并携出大量的磁能和热能;中心电流片中的磁岛向下运动,逐渐演变成底部含有３个磁闭合区的冕流结构。进而在中心电流片中再次发生磁场重联,多次形成向下运动的小型磁岛,并与底部磁闭合区发生结合不稳定性。同时在磁闭合区中也有磁场重联发生,导致中心小磁闭合区的湮灭。初始电流片之间的距离趋近,上述演化过程越快。日冕多层电流片中的磁场重联过程可能对日冕物质抛射和日冕加热有着重要影响。     

引力场中磁场重联的数值研究（IV）多层电流片

应用数值方法研究了日冕多层电流片中电阻撕裂模不稳定性的非线性演化和磁场重联过程,结果表明,计算区域顶部附近两侧电流片中的磁岛和等离子体团向上抛射,并携出大量的磁能和热能;中心电流片中的磁岛向下运动,逐渐演变成底部含有３个磁闭合区的冕流结构。进而在中心电流片中再次发生磁场重联,多次形成向下运动的小型磁岛,并与底部磁闭合区发生结合不稳定性。同时在磁闭合区中也有磁场重联发生,导致中心小磁闭合区的湮灭。初始电流片之间的距离趋近,上述演化过程越快。日冕多层电流片中的磁场重联过程可能对日冕物质抛射和日冕加热有着重要影响。     

各向异性等离子体中撕裂模不稳定性和磁场重联的数值研究

应用二维三分量时变可压缩磁流体动力学模拟方法,数值研究了各向异性等离子体(P_⊥≠P_∥,T_⊥≠T_∥)中的撕裂模不稳定性和磁场重联过程.计算结果表明,在短暂的线性增长之后,不稳定性将趋于非线性饱和.线性增长率随着各向异性程度|P_⊥/P_∥|增强而增大,随着等离子体β值减小及磁场y分量By增大而降低.在强垂直各向异性(P_⊥>P_∥)的情况下,电流片中磁场重联形成的磁岛达到非线性饱和后,在X型点附近形成空腔结构;随着空腔的增大,磁岛逐渐变小,并最终消失.在P_⊥>P_∥情况下,仅在电流片中心区域可以激发火蛇管不稳定性,电流片中不能形成大型磁岛.     

磁通量浮现驱动磁场重联的数值研究

应用二维时变可压缩磁流体动力学数值模拟,研究了双极－单极磁场中由于磁通量浮现驱动的磁场重联过程．结果表明,双极场与单极场间磁力线的重联形成上升的冷而密的等离子体团,磁场演变成鞭状结构．向上运动的等离子体团到达其最大高度后将回落和弥散.等离子体团最大上升速度达０．１４ＶＡ,等离子体最大上升速度达０．２７ＶＡ,ＶＡ为下边界处的Ａｌｆｖｅｎ速度．随着磁通量浮现幅度的增大,等离子体上升速度增加,重联过程发展得较快．背景等离子体β１值（β１为等离子体压力与磁压之比）越小,等离子体团中密度增量越大．磁Ｌｕｎｄｑｕｉｓｔ数Ｓ在１０３－１０６之间的改变对等离子体的速度和密度增量影响并不明显．与电阻撕裂模不稳定性引起的自发重联相比,磁通量浮现会更有效地驱动双极－单极场中的磁场重联过程,寻致日冕Ｈα冲浪和Ｘ射线喷流的形成．     

北京重力地球潮参数和重力仪力学特性的测定

根据北京白家疃测站８台ＧＥＯ潮汐重力仪对比观测的资料，采用标准线性体流变模型模拟重力仪及记录系统的力学特性，计算各台仪器的格值改正系数、主潮汐波的振幅衰减系数和相位滞后；确定北京重力地球潮的振幅因子δ和相让滞后为δ（Ｍ_２）＝１．１７８４，δ（Ｏ_１）＝１．１８０４，δ（Ｍ_２）／δ（Ｏ_１）＝０．９９８３，（Ｍ_２）＝０．０７°，（Ｏ_１）＝０．１０°．     

MHD激波解的讨论

从Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系出发，以激波切向磁场ξ、上游激波角θ_１和等离子体β_１值为参数研究各类激波解的特性及相互关系，阐明各类激波强度随介质β值的变化规律．结果指出：（１）在ξ＜－１的Ⅰ_ｂ型中间激波区和ξ＞４的快激波区存在双解；（２）慢激波可以直接和中间激波连接，但不能和快激波直接连接；（３）各类激波强度（用激波密度比衡量）随β_１值均有变化：Ⅰ_ｂ型中间激波上支随β_１值增加而下降，下支则上升；Ⅰ_ａ型中间激波和慢激波都随β_１值增加而下降；快激波双解区上支随β_１值增加而下降。下支上升；快激波在１＜ξ＜４区间的解随β_１值增加而上升。     

地球远磁尾中的磁场重联——（Ⅱ）局部驱动力的影响

基于可压缩磁流体动力学模型,数值研究了由太阳风引起的局部驱动力对地球远磁尾中磁场重联的影响.结果表明,在远磁尾等离子体片中将发生强迫磁场重联,并形成磁岛和等离子体团.形成磁岛的特征时间很大于流动撕裂模不稳定性引起磁岛非线性饱和的特征时间.磁岛宽度随着磁Reynolds数S的增大而减小,随着尾瓣中等离子体压力与磁压之比值β_∞的降低而减小.认为太阳风引起的局部驱动力对地球远磁尾等离子体片中磁场重联的影响,可能不如流动撕裂模不稳定性那样显著.     

岩石磁化率各向异性的频散现象

在磁化磁场为恒定磁场（无定向磁力仪─-Ａ）与１０ｋＨｚ交变磁场（ＭｉｎｉｓｅＰ旋转磁力仪─-Ｍ）的两种情况下，作岩样比测时发现，后者测得的主磁化率值к^Ｍ较小，而且沿磁性叶理面方向ｆ的主磁化率值к_ｍａｘ与к_ｉｎｔ比沿其法线方向ｎ上的к_ｉｎｔ减少得更多．作者以涡流磁滞机构对磁性载体材料的“磁屏蔽”作用发生了各向异性的影响来解释此现象；并推测磁化率各向异性存在频散现象．     

引力场中磁场重联的数值研究（Ⅲ）等离子体团的形成和运动

应用解析和数值相结合的方法求得了含有闭场区、中性片和开场区的二维局域磁静力平衡日冕基态，进而在此基态下数值研究了由电阻撕裂模不稳定性引起的磁场重联过程．结果表明，在电流片长度远大于其半宽度的情况下，仍将发生具有两个Ｘ线的磁场重联，形成磁岛和高温高密度的等离子体团．等离子体团运动的方向取决于闭合区底部等离子体压力和磁压的比值β０．当β０较大时，等离子体团向下运动，引起结合不稳定性和磁岛的合并，等离子体团中的等离子体不断落入闭合磁场区两侧．当β０较小时，等离子体团向上运动，导致电流片中等离子体的喷发．结果还表明，磁张力的分布是决定等离子体团运动方向的主要因子．     

2008年1月5日磁层亚暴期间磁尾近地重联事件的THEMIS观测与分析

本文主要应用THEMIS卫星的磁场和等离子体流观测数据,分析了2008年1月5日08∶51~08∶57 UT亚暴膨胀相期间磁尾的一个近地重联事件.在亚暴膨胀相期间,地面的全天空成像仪清楚地记录到了极光的极向扩展,THEMIS的P5卫星在地球同步轨道附近观测到了磁场的偶极化现象.在亚暴膨胀相末期的08∶51~08∶57 UT期间,P3(X_GSM~-9.12R_E) 和P4 (X_GSM ~-9.40R_E) 同时观测到了一对方向相反的高速等离子体流.这对方向相反的高速等离子体流是由磁尾的重联现象所引起.重联的位置被估计位于X_GSM ~-9.12R_E 和X_GSM~-9.40R_E之间较小的空间范围内.并且,在重联位置的两侧,重联的Hall效应被P3和P4两颗卫星观测到.因此,这一磁尾重联事件发生在距离地球非常近的空间范围内.     

地磁场Sq的经度效应和UT变化

使用国际地球物理年（ＩＧＹ）和国际地球物理合作年（ＩＧＣ）期间全球地磁台网的资料和中国台站的同期地磁记录，对地磁场太阳日变化Ｓ_ｑ的经度效应和ＵＴ变化进行了分析研究．将适用于瞬时全球磁场的球谐分析法与适用于平均场的双调和分析法相结合，提出了一种分析Ｓ_ｑ磁场的新方案，得到了（θ，Ｔ，ｔ）坐标系中组成Ｓ_ｑ场的三个部分，即仅随地方时ＬＴ变化的部分Ｓ_ＬＴ，仅随世界时ＵＴ变化的部分Ｓ_ＵＴ以及既随ＬＴ变化又随ＵＴ变化的部分Ｓ_ＬＵＴ．从全球来看，在ＩＧＹ／ＩＧＣ期间，这三部分的强度之比分别为：Ｘ分量１．０：０．２：０．３，ｙ分量１．０：０．１：０．６，Ｚ分量１．０：０．４：１．０．Ｓ_ｑ的经度效应和ＵＴ变化虽然表现了同一物理过程，但它们有不同的表达形式，其决定因素是地磁轴对地理轴的倾斜，地磁场的非偶极子成分（即区域性异常）和地球内部电性的横向不均匀性（包括海陆分布、地壳上地幔电导率的区域差异）     

低轨道磁化等离子体中运动航天器等离子体鞘层特性

用二维静电粒子模拟程序研究了在低轨道磁化等离子体中航天器充电过程和等离子体鞘层结构．采用磁化等离子体流经航天器的物理模型,并且考虑零级对流电场的影响．结果表明航天器在低轨道空间飞行时,如果不考虑光电子流,航天器会在短时间内(ｔ＝５４ω_ｐｅ^－１,ω_ｐｅ是电子等离子体频率)快速充电到平衡电位。在不考虑航天器运动的情况下其周围等离子体鞘层结构是围绕航天器对称的,等离子体鞘层宽度约２．４个德拜（Ｄｅｂｙｅ）长度．当考虑航天器运动时,航天器的平衡负电位绝对值减小．且在背离航天器运动方向上等离子体鞘层宽度为８．３个德拜长度．磁场的增大使负的平衡电位绝对值变小,当航天器运动从平行于磁场方向转向垂直于磁场方向时,航天器负平衡电位的绝对值减小．     

准周期性磁结构的形成过程及其特征的数值分析

在入流马赫数ＭＡ为０．１－０．４范围内，选取不同的入流，逐一考察在可压缩导电介质中磁场重联的类型随磁雷诺数Ｒｍ的变化．大量的数值结果表明：在不同入流驱动下，磁雷诺数（Ｒｍ）均有临界值存在，如果等离子体系统逐渐趋于稳定的单Ｘ线重联；若间歇性的次级撕裂产生了重复出现的磁结构，而且Ｍａ越大产生间歇性次级撕裂的临界值越高，与ＭＡ之间基本上符合线性关系．此外还发现，当多重Ｘ线重联间歇性地出现时，有关物理量发生准周期性地振荡．     

驱动重联过程中等离子体团的聚合

本文对驱动重联过程作了二维可压缩磁流体动力学模拟．对于长的计算域（１：４），在持续入流的作用下形成了双重磁岛结构，这些磁岛在向外运动的过程中不断靠近，逐渐聚合成为一个大的等离子体团，并继续向外运动．这说明在地球磁层顶和磁尾所观测到的大尺度磁结构，很可能是在驱动重联过程中，由较小的等离子体团聚合后形成的．并初步揭示了在空间等离子体中，磁岛的聚合过程难以直接观测的动力学原因．