各向异性磁流体激波跃变关系的无量纲化求解

用无量纲化方法，推出了各向压力异性等离子体中磁流体激波跃变条件及其解析解的一种普遍而简单的形式，并讨论了它在各向同性等离子体中各种特殊情况下的结果．     

行星际激波传播至1AU的渡越时间预测

考虑了激波爆发源角宽度、能量、驱动时间、激波速度及其与背景太阳风之间的相互作用,利用流体力学扰动方程建立起一个激波扰动传播模型,用于研究激波从太阳传播到地球轨道附近(1 AU处)所需要的时间(渡越时间)问题.为印证扰动传播模型的适用性,利用1979-1989年间的27个激波事件,以及1997年2月到2000年1月间的68个激波事件,对激波到达地球轨道附近的渡越时间进行了预测,并将结果与STOA和ISPM预报模型结果进行了比较.实验表明,该模型在所有95个事件中,渡越时间相对误差小于10%的事件数占总事件数的25.26%;相对误差小于20%的占总事件数的50.53%;相对误差小于30%的占总事件的65.26%.      

行星际结构与垂直无碰撞激波的相互作用

应用一维混合模拟方法数值研究了两种行星际结构──反向磁场和高密度等离子团与垂直无碰撞激波的相互作用．结果表明，随着激波上游区磁场的反向，下游区磁场将逐渐改变符号，且等离子体密度和速度分别呈现较强的湍动．激波上游和下游的物理量依然满足Rankine-Hugonoit关系．当高密度的等离子体团通过垂直无碰撞激波时，部分质子被激波反射，部分质子被加速并进入下游区域．由于质子速度分布为非Maxwell分布，在激波下游也激发出较强的湍动．     

运用一维太阳风MHD CE/SE模型预报行星际激波到达时间

在一维太阳风磁流体(MHD)数值模拟中, 应用时空守恒元和解元(Conservation Element/Solution Element, CE/SE)方法, 建立了一个行星际激波扰动传播模型(1D-MHD (CE/SE)模型), 用来预报行星际激波到达时间. 收集了1997年2月至2002年8月间的137个激波事件, 对激波到达地球轨道附近的传播时间进行了预测, 并将预报结果与STOA, ISPM, HAFv.2以及SPM模型所得结果进行比较. 对于相同的样本事件, 1D-MHD (CE/SE)模型给出的渡越时间平均绝对值误差并不大于其他4个模型, 且该模型预报的相对误差小于10 %的事件占25.6 %, 小于30 %的事件占69.3 %, 小于50 %的事件占87.6 %, 其预报精度与其他模型相比基本相当. 这表明该模型在空间天气的激波到达时间预报方面有潜在的应用价值.      

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

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

无碰撞激波上游低频波的激发及其对激波耗散机制的影响

由于很大一部分来自激波上游的粒子被激波面所反射,因此在准平行无碰撞激波的上游存在着等离子体束流．通过一维混合模拟方法,计算了束流密度较小(nb／n0=0．02)和较大(nb／n0=0．2)两种不同情形下的等离子体束流不稳定性．结果表明,在束流密度较小时,束流激发的主要是平行于背景磁场方向传播的右旋波,此波动只能对束流粒子产生影响,包括减速和加热．在束流密度较大时,束流可同时激发平行和反平行于背景磁场方向传播的右旋波,除能对束流粒子产生影响外,还可通过非共振作用加速和加热背景粒子．文中对准平行无碰撞激波耗散机制的影响也进行了讨论。     

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

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

磁云边界层中的重联慢激波观测分析

在Petschek模型中，排空区边界处的一对慢激波是能量耗散的重要机制.已有大量行星际空间的Petschek型磁场重联排空区观测事件被报道，但是只有少量的排空区边界处观测到了慢激波.针对一例位于磁云边界层中的Petschek型磁场重联排空区观测事件，在排空区靠近磁云一侧边界处证认了一例慢激波.激波跃变层两侧的磁场和等离子体参数满足Rankine-Hugoniot关系，且激波上下游的中间马赫数均小于1，上游的慢马赫数为2.94（>1），下游的慢马赫数为0.65（<1），符合慢激波的观测特征.磁云内部的等离子体β值很低，局地阿尔芬速度高，同时磁云边界层中可能发生丰富的磁场重联活动，这可能是磁云前边界处慢激波形成的原因.      

r型网格自适应在间断Galerkin有限元激波捕捉中的应用

间断Galerkin（DG）有限元方法因计算精度高、适用于非结构网格等特点得到广泛研究和应用，其在数值模拟包含强间断流场时存在残差收敛性和计算鲁棒性差问题，均匀分布的网格加剧这一问题并影响激波分辨率。针对该问题，发展了r型网格自适应方法，实现间断Galerkin有限元数值模拟过程中网格自适应加密。基于网格点归一化的压力值作为r型网格自适应中网格点移动驱动力的重要权值，并将网格自适应后的网格点位移变化量与网格点之间的初始位移之比作为驱动力的另一重要权值，实现网格沿激波方向各向异性自适应加密，并且激波附近网格点的相邻网格点同步向激波方向移动。发展了适合间断Galerkin有限元方法的Venkatakrishnan限制器。并列NACA0012翼型超声速算例及三维并列圆柱相互干扰算例结果表明:基于r型网格自适应的间断Galerkin有限元方法能够清晰锐利捕捉激波，提高模拟精度，具有良好的收敛性和鲁棒性。      

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

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

磁流体斜激波的碰撞

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

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

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

磁流体斜激波的汇合

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

耀斑激波传播的数值研究

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

激波在行星际介质中的能量耗散

导出激波下游介质相对上游介质能流通量增量公式,并由HeliosA,B飞船太阳风观测资料得出不同流速太阳风流中各参量随日心距的幂律变化。以此作为背景值分别计算出磁能、内能、动能和总能在不同日心距离处的能量耗散率。结果指出激波后介质以动能增加为主,内能次之,磁能最少;总能耗率在近日处较大,但下降较快。从0.3-1.0AU,不同强度激波总能耗随初始Alfvén激波数A₁₀增大而增大,对A₁₀从2.0-10.0的计算结果与观测值一致。      

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

为了模拟太阳色球层、过渡层和日冕中扰动的传播特征,本文在一自洽的非等温、非均匀等离子体初态下,应用二维时变可压缩磁流体动力学模拟,数值研究了过渡层附近局地加热引起的扰动传播过程.结果表明,扰动分别以局地快磁流波速度向四周传播,数值模拟结果与局地快磁流波推算结果符合得很好,其特征可以解释SOHO/EIT观测到的波动事件.      

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

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