冕环中扭转Alfvén波和快波的耦合共振

推广完备正交函数系(OFSE)展开方法,应用它研究冕环中无耗散扭转Alfven波和快波耦合的时变过程.采用这种方法,从数学上描述了总体波模的固有角频率,从理论和数值计算结果两方面分析了Alfven波和快波耦合时冕环中磁力线Alfven波固有角频率.采用这种方法,分析了耦合驱动项和Alfven波共振的关系.计算结果给出了耦合共振出现位置Alfven波和快波的波幅特征,发现如果足点驱动角频率不等于冕环总体波模的固有角频率,在耦合共振位置当扭转Alfven波幅出现δ断面时,快波振幅沿径向梯度很大,有时近似为一个间断,非常有利于快波耗散;如果足点驱动角频率等于总有体波模的固有角频率,Alfvn波振幅出现δ断面时在径向出现丰富的小尺度结构,共振位置近似为一个间断,非常有利于Alfven波耗散.     

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

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

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

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

无振荡、无自由参数格式在三维太阳风流动数值模拟中的应用

从快速、有效地进行空间天气数值预报的需要出发,针对1998年5月份行星际太阳风暴事件,采用全三维流体力学(HD)模型进行数值试验,考察了无振荡、无自由参数(NND)格式在三维太阳风流动数值模拟中的应用.计算中内边界的密度分布由观测K日冕数据来确定,由此得到的源表面密度分布具有和源表面电流片相似的结构.数值试验表明: 虽然该格式不需要加入人工粘性,但是具有很好的计算稳定性.在扰动计算时激波间断耗散小,在三维任何方向上间断所占的网格数比较少,没有数值色散现象.      

1997年1月空间天气事件的三维MHD模拟

以1997年1月空气天气事件期间的观测为依据，在构造了比较接近真实的背景太阳风基础上，进一步利用三维时变的MHD模式，模拟了CME（日冕物质抛射）激发的扰动在行星际空间的传播过程，对地球空间环境的影响及行星际磁场南向分量Bz在1AU的时间经历。模拟结果与WIND卫星的测量进行了比较。结果表明，模拟与观测得到的扰动得到地球的时间、地球空间环境各量的变化及Bz的时间经历基本一致。     

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

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

嵌套闭磁场日冕背景结构对触发CME特征的影响

在子午面内,研究具有嵌套闭磁场结构冕流背景对触发日冕物质抛射(CME)特征的影响.在冕流背景磁场结构内有三个小尺度的闭磁场结构,其中中间小尺度闭磁场结构的磁场方向和冕流整体偶极磁场方向相反.CME触发模型在这个小尺度闭磁场结构下方浮出,它具有同心圆形磁场结构,半径为a=0.1Rs(Rs为太阳半径);CME触发模型前半部分磁场方向和触发位置处小尺度闭磁场结构的磁场方向相反,与冕流整体偶极磁场方向相同.数值模拟结果表明,当CME触发模型中心等离子压力与边界压力之比m≥2时,上浮CME触发模型可以触发CME;当m＜2时,上浮CME触发模型不能触发CME,计算结论的误差小于1%.      

太阳风中磁场重联离子扩散区的大尺度特征

利用ACE和WIND卫星2007年1月6日的联合探测, 在1AU附近发现了一个等离子体密度极低的Petschek-like重联喷流区. 该喷流区内部出现了非常明显的Hall双极磁场、等离子体密度下降区以及与Hall电流相符的低能段电子投掷角分布. 这些特征与重联离子扩散区的Hall效应非常吻合, 说明很可能在太阳风中观测到了一个离子扩散区. 分析表明, 与之相关的磁场重联为准稳态快速完全反向重联, 其扩散区以一对慢模波为边界, 空间尺度达到80个离子惯性长度, 表现出了大尺度重联的特征.      

行星际日冕物质抛射地磁效应研究的支持向量机方法初步研究

行星际日冕物质抛射（Interplanetary Coronal Mass Ejection,ICME）与地球磁层相互作用并带来地磁暴等地磁扰动.从Richardson和Cane提供的近地球ICME列表中筛选出ICME事件集,基于ICME扰动期间的行星际等离子体与磁场数据提取出特征.通过计算各特征的费舍尔分值（Fisher Score）,对这些特征进行选择,发现行星际磁场南北向分量持续时间小于-10nT且激波等扰动所带来的ICME扰动开始时,太阳风速度的增量等特征与ICME事件的地磁效应密切相关.这与现有的传统统计研究结果一致.以这些特征为基础,训练得到的径向基函数支持向量机能够以0.78±0.08的准确率判断ICME事件是否会产生中等及以上强度的地磁暴（Dst ≤-50nT）.      

基于CPU/GPU的日冕偏振亮度并行计算模型

三维磁流体力学(MHD)数值模拟是用来研究日冕和太阳风最常用的方法之一, 其中将计算得到的日冕电子数密度转化为日冕偏振亮度(Polarization Brightness, PB)是与观测对比的重要方法. 由于待转换电子数据网格密度、PB数据网格密度和计算模型的复杂度, 使得日冕偏振亮度的计算比较耗时, 利用单CPU计算无法达到近实时转换日冕偏振亮度的要求, 从而影响了数值模拟的验证效率. 本文在CPU/GPU环境下, 利用CUDA编程技术, 提出了一个日冕偏振亮度并行计算模型. 实验结果表明, 该模型比CPU上的串行模型计算速度提高了31.86倍, 达到了近实时模拟与观测数据比对的计算要求.