运用一维太阳风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 %, 其预报精度与其他模型相比基本相当. 这表明该模型在空间天气的激波到达时间预报方面有潜在的应用价值.      

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

色散和耗散对Alfven孤立波的影响

利用变数变换和解析积分的方法,分别找到了带色散项的DNLS方程和带耗散的MKdV方程的一类解析解,给出了存在孤立波的条件。      

两冕流间CME事件数值模拟的改进

通过偶极子场和六极子场适当叠加，改进猜解磁场，使猜解磁场在太阳南北极符号相反，然后采用理想磁流体力学方程组(MHD)，由猜解磁场与太阳风流动相互作用计算出稳态自洽解，得到定性上与观测比较接近的具有两个冕流的背景结构．在两个冕流间采用具有同心圆磁场位形的触发模型触发CME事件，研究CME的日冕传播特征．模拟结果表明，CME被约束在两冕流间传播，CME闭磁场位形和磁云横截面磁场位形相似，可以解释1AU处观测磁云的部分特征；在CME附近，存在压力和Lorentz力起主要作用的区域，这可以为分析1AU处CME事件的观测数据提供帮助．     

行星际磁云边界层中的电子分布函数特征

基于WIND 飞船观测的1995---2006 年间的磁云事件, 研究了磁云边界层中电子的流动图样, 以及电子速度分布函数的特点与电子加热和加速的关系, 得出以下结果. ①磁云边界层中存在的电子流动图样, 包括各向同性、双向流动和单向流动等. ② 相比于背景太阳风和磁云本体, 通常情况下磁云边界层中电子分布函数的核心热电子成分 (E<60eV) 增多, 超热电子成分 (E>60eV) 在沿磁场垂直方向上增多, 而在沿磁场平行或反平行方向以单方向增加为主, 此外, 还在近1/10 的磁云边界层中观测到了高能电子的明显增多. ③对比研究了磁云边界层与磁云驱动的激波对电子速度分布函数的调制作用, 经过激波, 电子分布函数的超热电子成分在各方向上都有增加, 不同于磁云边界层中在沿磁场平行或反平行方向上超热电子成分以单方向增加为主, 表明二者有不同的形成机制. ④考察了磁云边界层中的波活动增强和电子分布函数及离子流量增加的对应关系. 上述观测和对比分析进一步表明了磁云边界层是一种重要的动力学结构, 磁重联是一种可能的形成机制.      

嵌套闭磁场结构内CME产生和传播的数值模拟

给出了特殊类型的日冕物质抛射(CME)数值模拟定性结果,这种CME核心闭磁场结构前半部分磁力线的方向与太阳整体偶极场磁力线的方向相反．计算结果表明,这种CME核心闭磁场结构磁力线与太阳整体偶极场反向磁力线之间存在过渡磁场结构,在向外传播时过渡磁场结构所占的面积逐渐增大．这一结果可以用来解释飞船为什么能够观测到一类双极磁云,这类磁云前半部分磁场方向与太阳整体偶极场方向相反．为了模拟这一数值结果,强调需要采用包含嵌套闭磁场的冕流背景结构,并在合适的位置触发CME．     

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

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

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

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

无振荡、无自由参数格式在数值模拟盔型磁场中太阳风流动的应用

提出了一种离散二维三分量理想磁流体力学守恒型方程组（MHD）的算法（NNDMHD），它有效地控制了磁场散度不为零误差对动量方程组的影响，把气体动力学中计算跨音速流动问题的有效算法——无振荡、无自由参数（NND）格式推广应用到MHD方程组中。利用该算法首先对常见一维和二维算例进行数值试验，得到比较好的结果，消除了间断处的非物理振荡。然后对太阳风在子午面轴对称盔形磁场位形中流动进行数值试验，在这个算例中，物理量沿径向变化大，NNDMHD格式仍然能够有效地控制磁场散度离散不为零误差导致的非物理流动。这个算例的计算结果表明：在网格划分比通常情况和稀4倍时，该算法仍保持很好的计算稳定性。     

－类非线性系统反馈镇定的充要条件

讨论－类非线性系统全局渐近反馈问题．文内通过构造显示状态反馈和使相应闭环系统全局渐近稳定的显示李氏函数，而得到了使系统全局渐近反馈镇定的充分必要条件．