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

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

中高层大气对低层大气脉冲扰动的响应

从基本的大气运动方程出发,采用二维全隐欧拉格式(FICE)建立极坐标下二维声重波传播的数值模式.对小振幅高斯形态声重波传播过程的模拟结果表明,该数值模式能够再现声重波波包向上的稳定传播过程.应用此数值模式模拟了在没有背景风场和存在背景风场条件下中高层大气对10km高空处近似脉冲点源扰动的响应过程,获得了中高层大气对低层大气脉冲扰动响应的图像.结果表明,电离层高度对重力波扰动的响应不但出现在激发源激发之后,还与激发源具有较大距离,扰动的幅度随高度的增大而增大;有背景风时,顺风情况下重力波波动传播路径较逆风情况下平缓,且传播的水平距离比逆风情况下远,波动的振幅较逆风情况下弱,逆向背景风场能够加速重力波向上传播.同时,模拟还表明脉冲扰动正上方可以激发产生一种以6 min为主要周期,具有声学分支特性的快速波动.本文模拟的扰动结果与2004年12月26日苏门答腊-安达曼地震引起的电离层TEC扰动具有类似的图像.      

多离子空间等离子体中快磁声波与粒子的回旋共振特性

快磁声波是空间等离子体中一种接近垂直传播的右旋极化电磁波,能够在等离子体层内外传播.快磁声波与带电粒子的回旋共振相互作用能够导致高能电子随机加速和投掷角扩散、能量质子投掷角扩散等,从而影响辐射带高能带电粒子的动态过程.分别基于完整的色散关系和高密度近似的色散关系,在不同空间等离子体条件下研究多离子空间等离子体中不同传播角的快磁声波色散曲线,并计算了快磁声波与H+,He+和O+离子的最小共振能量.结果表明,当传播角较小时,采用高密度近似与采用完整色散关系计算的离子最小共振能量没有太大差别.在中低密度中强磁场空间等离子体中,传播角≥ 88°时高密度近似色散关系会带来很大的误差,因此应利用完整色散关系计算最小共振能量.      

温度扰动激发的重力波波包非线性传播过程的数值研究

采用谱配置方法建立了一个中层大气扰动传播的全非线性动力学数值模型，利用该模型对具有高期分布的波状温度扰动在可压大气中的演变和传播过程进行数值模拟，模拟结果表明：温度扰动很快就会激发出两支重力波波包，其中一支向上传播，而另一支则向下传播；这两支重力波波包的传播路径与线性理论给出的射线路径有些差异，进一步的数值分析则表明，给定的温度扰动在经过约3h的演变后，完全转化为波动，在这个转化过程中，79％的扰动能量转化为波动能量。     

波状扰动激发的重力波波包在中层大气中能量传播和转换特性的数值研究

采用谱配置方法对具有高斯分布的波状扰动激发重力波过程中的能量传播和转换特性进行数值研究．模拟结果表明：初始时刻单一的水平风、垂直风、密度或温度的波状扰动都可以激发出两支重力波波包，其中一文向上传播，而另一支则向下传播；在相同的波参数下，4种波状扰动所激发的重力波波包的能量传播路径几乎完全一致，与线性理论给出的射线路径却有些差异．进一步的数值分析则表明：4种波状扰动转化为波动的特征时间基本相等，约为2个波动周期，但是在这个转换过程中，不同扰动的能量转换效率却差别很大．     

太阳大气中有限振幅磁场扰动的动力学过程

本文从完整的磁流体动力学方程组出发,通过太阳大气中磁力线管根部有限振幅磁场的扰动,研究了非线性磁场的动力学演化。假设初始磁场位形足β<<1的势场,根部磁力线管磁场扰动,驱动等离子体运动,一部分磁能转换为等离子体动能。等离子体压缩运动具有快磁声波的特性。计算结果给出非线性磁场演化的定量关系,可以解释太阳大气中日冕活动过程。也可用于模拟实验室里高β实验装置中的等离子体的持性。      

有限振幅扰动下高超声速钝楔绕流感受性

为了研究自由流脉冲波作用下高超音速流动感受性特征,采用高阶有限差分方法模拟钝楔超高音速绕流中来流慢声脉冲波作用下的非定常流场,并应用傅里叶频谱分析了边界层扰动模态的演化.结果显示,自由流脉冲扰动下,边界层不同模态波沿流线呈现不同变化规律.球头区边界层感受到的主要是基频波,在球头向非球头区过渡段,主导模态迅速向高频迁移;在非球头区,主导模态频率缓慢向高频迁移.总的来说,沿流向总扰动模态中的低频模态成分变小,高频,特别是二阶谐频附近模态成分显著地增大.沿流线发展,边界层不同扰动模态之间存在伴随能量迁移的模态竞争物理现象,主导模态的发展对其他模态有很强的抑制作用.      

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

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

近地微电离大气中阿尔芬波的截止

利用全电流双流体模型分析了地面磁扰动信号在近地层大气中的传播.结果表明磁扰动信号的传播模式受中性粒子elax-elax离子碰撞频率、磁扰动频率、磁粘滞系数以及大气离化系数等参数的影响.中性粒子碰撞效应导致近地层大气中不存在阿尔芬波传播通道.当磁扰动频率远小于0.001Hz时,磁扰动传播存在趋肤效应,且在近地大气中从高往下逐渐增强;当磁扰动频率远大于0.001Hz时,磁扰动以类似于真空电磁波的模式在近地大气中传播.      

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

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

Alfven波在微电离大气中的衰减特性研究

利用简单的偶极子地磁场模型以及大气电子密度和电导率模式, 分析地面产生的磁扰动以Alfven波的模式传播到近地空间区域. 这种地面的磁扰动可能干扰近地空间卫星对空间磁扰动的观测. 通过对地面磁扰动Alfven波模式1000 km高度内的衰减情况进行模拟, 认为在近地空间采用地磁偶极子模型是合理的. 由于衰减随扰动频率的增大而急剧增强, 分析还得到了近地卫星能够探测到地面磁扰动的最大频率. 计算结果表明, Alfven波的衰减主要集中在高度50 km以下, 这个区域内的大气电导率极其微弱, 使Alfven波的传播受到极大衰减. 0.4 Hz以下的Alfven波沿磁力线传播到1000 km高度后衰减结为原来扰动幅度的千分之一, 因此频率在0.4 Hz以下的Alfven波可能会干扰低轨卫星探测磁场脉动.      

我国中部地区电离层扰动时空尺度的统计分析

将最大熵动态互普方法用于高频多普勒台阵数据处理，从覆盖太阳活动高、低年份的高频多普勒台阵观测数据中，得出了中国中部地区电离层声重波扰动的传播方向、时空尺度等参量之间的统计分布特性，进而探讨了有关扰动源的情况，结果表明，极区源和局地源两种不同的激发源所激发的声波，分别反映了电离层扰动的全球传播和地区特性两种不同的特征。     

双变量波状扰动激发重力波波包的数值研究

采用谱配置方法对双变量波状扰动激发重力波的过程及其能量传播和转换特性进行数值研究．模拟结果表明，当初始给定的两个变量的波状扰动的相对振幅和相位严格满足重力波的偏振关系时，只有水平风场扰动(u′)与温度扰动(7′)的共同作用才能抑制能量向下传输．通常情况下，两个变量的波状扰动都会激发出两文传播方向相反的重力波波包，这两支波包的相对能量明显取决于两个波状扰动之间的相对振幅关系．当两个扰动之间的相对振幅满足特定条件时，向下传播的能量几乎为零，只会有一支向上传播的重力波波包被激发．进一步的数值分析则表明，各种情况下两个变量的波状扰动所激发的重力波波包的能量传播路径和速度基本相同；与单变量波状扰动激发重力波波包的过程相比，能量传播路径几乎完全一致，扰动转化为波动的特征时间也基本相同，但能量转换效率更高．     

2009年2月13日EUV波现象数值研究

结合STEREO卫星的观测和三维磁流体力学数值模拟方法, 采用WSO (Wilcox Solar Observatory)磁场数据和势场源表面模型建立日冕初始磁场, 并在日面活动区加上时变的压强扰动, 对2009年2月13日05:35UT爆发的CME-EUV波(Coronal Mass Ejections-Extreme Ultraviolet wave, 日冕物质抛射-远紫外波)事件进行研究. 从COR1/STEREO-A图像判断, 此次CME前沿速度约340km·s-1, 角宽度约60°; 分析EUVI/STEREO-B 195 ?的差分图像, 可以看到, 环形亮环波前从活动区向四周传播, 亮环波前后面是日冕暗化区, 取四个方向的波前位置进行线性拟合可知, 该EUV波速度为247km·s-1, 数值模拟得到的EUV波速度为245km·s-1, 将计算结果采用IDL可视化后可以看到明显的亮环和暗区结构, 数值模拟结果与卫星观测相一致, 表明该EUV波现象是快磁声波.      

重力波波包在三维非等温大气和水平风场中非线性传播的数值研究

采用全隐欧拉格式（FICE）对重力波波包在三维非等温大气、均匀和剪切风场中的非线性传播进行了数值模拟，给出了重力波波包三维非线性传播的全过程，分析了重力波的传播特性及背景温度场、风场对重力波传播的影响。结果表明：波包扰动速度振幅的增长比在WKB条件下振幅的增长要慢；波包非线性传播的路径、能量传输速度不同于WKB近似下的结果，非线性效应导致了重力波的传播特性的改变；温度场的非均匀性会改变重力波传播的路径和速度；剪切风场使扰动速度振幅的增加变得缓慢，使垂直波长减小。     

空间尘埃等离子体中尘粒电荷的相关涨落对尘埃磁声波的影响

本文运用自洽的三流体方程组, 考虑了尘埃的充电过程, 得到均匀磁化尘埃等离子体中垂直于磁场传播的尘埃磁声波的色散关系, 结合空间环境讨论了尘埃电荷的相关涨落对尘埃磁声波的影响.      

太阳风中低频离子静电孤波

本文给出在太阳风超声速流动条件下,离子静电孤波的传播特性,结果与Helios1,2卫星观测的静电离子噪声做了比较。离子声波扰动的非线性发展使太阳风等离子体呈规则的小尺度起伏,离子声波在马赫锥外传播,因此理论预言密度起伏不沿着太阳风速度方向,而是在横向方向.      

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

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

行星际空间MHD模拟激波跃变关系分析

分析了模拟得到的可以传播到1AU以远的日地空间磁流体力学激波与Rankine-Hugnoniot跃变关系的符合程度．通过对模拟激波的结构及其在传播过程中的演化进行的分析，提出了模拟激波的定位方法；基于所提出的定位方法，利用向长青提出的确定MHD激波局地参数的方法计算了模拟得到的激波与Rankine-Hugnoniot跃变关系的偏差．结果表明在激波传播到100Rs以后，激波对中前向快激波与Rankine-Hugnoniot关系的符合达到很高的程度，相对误差在10^-2数量级以内；并且在激波传播到150 Rs以后，相对误差在10^-3数量级以内．这个结果说明文中所使用的有限差分数值格式能较好地模拟激波．     

低电离层中的行星波扰动及其垂直传播

根据武昌电离层观象台接收的罗兰－Ｃ ＬＦ无线电波的传播资料，利用小波分析方法，研究了１９８５年和１９８６年冬季的低电离层中的行星波扰动及其垂直传播。结果表明：（１）两个冬季中，在７５ｋｍ和９５ｋｍ高度上，均存在明显的行星波扰动，主要成分是５－１０天和１０－２０天周期的波动，扰动幅度随高度增长，但增长值小于理论增长值；（２）在７５ｋｍ高度上的扰动和９５ｋｍ高度上的扰动，具有十分相似的扰动形态。根据行