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

首先讨论了磁流波传播的线性特征，然后构造了球坐标中-自治的非等温，非均匀等离子体初态，应用二维时变可压缩磁流体动力学模拟，数值研究了色球层底部压力脉冲所引起扰动的全球传播过程，结果表明，在极区，压力脉冲导致的扰动传播可以区分成两类不同模式的波动，快磁声波与慢磁声波，而在赤道附近，传播扰动是快模磁声波，在源区附近还存在一非传播的扰动，模拟结果的特征有助于解释SOHO/EIT观测到的波动事件。     

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

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

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

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

非等间距显式算法在重力波非线性传播数值模拟中的初步应用

采用3阶TVD型显式时间积分方法,建立了2维可压缩大气中重力波非线性传播的非等间距显式数值模式.对小振幅重力波传播过程的模拟结果表明,该模式能够很好地再现小振幅重力波的传播过程,并能保持能量守恒关系和各扰动量之间的相位关系,与线性重力波理论预测的结果吻合很好.对有限振幅重力波的模拟结果表明,与格点数相同的等间距网格中的模拟结果相比,采用非等间距网格能够以更高的分辨率模拟重力波的不稳定传播直到破碎的整个过程;在非等间距网格中的模拟结果与在加密一倍的等间距网格中的模拟结果相同,说明了用非等间距网格模拟重力波的饱和与破碎过程是正确的,并且采用非等间距网格能够大大减少计算量.      

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

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

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

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

日冕扰动传播的数值模拟——(Ⅱ)热传导和辐射的影响

利用包含完全能量方程的磁流体力学方程组,数值研究了磁静力平衡基态下由于日面冷物质径向喷射所引起的日冕动力学响应。这些结果可以解释日珥的形成及伴随的日冕扰动传播。结果还表明,热传导和辐射损失对扰动传播的基本性质影响很小;对于日珥的形成,也未显出重要作用。      

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

日冕扰动传播的数值模拟(Ⅲ)太阳风的影响

本文数值研究了二维磁流体动力学平衡基态下开场区日面冷物质径向喷射所引起的日冕动力学响应。结果表明:(1)在高密度环前方有一弱扰动区近似以Alfvén速度向外传播;(2)高密度环前缘移动速度随着径向距离而增加, 其增加值近似为局地太阳风速度;(3)高密度环中等离子体的最大径向速度约在4个太阳半径处趋于局地逃逸速度;(4)对于强开放场, 环形结构在θ方向上没有明显的扩张。这些结果可以更好地解释伴随有日珥的日冕物质抛射事件。      

Some characteristics of propagation of flare- or CME-associated interplanetary shock waves

Many interplanetary shock waves have a fast mode MHD wave Mach number between one and two and the ambient solar wind plasma and magnetic field are known to fluctuate. Therefore a weak, fast, MHD interplanetary shock wave propagating into a fluctuating solar wind region or into a solar wind stream will be expected to vary its strength.It is possible that an interplanetary shock wave, upon entering such a region will weaken its strength and degenerate into a fast-mode MHD wave. It is even possible that the shock may dissipate and disappear.A model for the propagation of a solar flare - or CME (Coronal Mass Ejections) - associated interplanetary shock wave is given. A physical mechanism is described to calculate the probability that a weak shock which enters a turbulent solar wind region will degenerate into a MHD wave. That is, the shock would disappear as an entropy-generate entity. This model also suggests that most interplanetary shock waves cannot propagate continuously with a smooth shock surface. It is suggested that the surface of an interplanetary shock will be highly distorted and that parts of the shock surface can degenerate into MHD waves or even disappear during its global propagation through interplanetary space. A few observations to support this model will be briefly described.Finally, this model of shock propagation also applies to corotating shocks. As corotating shocks propagate into fluctuating ambient solar wind regions, shocks may degenerate into waves or disappear.     

数字开关液压系统管路压力波传播建模与分析

为准确描述数字开关液压系统管路压力波传播特性，应用基于传递函数与时间延迟的管路压力波传播分析模型，通过耦合高速开关阀输出流量特性与管路压力波传播分析模型，实现两位两通高速开关阀控缸系统管路压力波传播特性的建模与分析，讨论了不同参数下的管路压力波传播特性，并通过相应的实验加以验证。结果表明，实验结果与分析模型结果一致性较好，高速开关阀动态特性与系统管路压力波传播特性密切相关，管路属性对系统管路压力波传播过程影响显著，将为后续数字开关惯性液压系统的研究奠定基础。      

二维有限长度柔性壁面上T-S波演化的数值研究

受自然界鸟类羽毛的柔性特征启发,利用数值模拟的方式开展了柔性壁面对亚声速边界层中T-S波演化的影响研究。刚性壁面上的数值结果与线性理论吻合得很好,验证了数值方法的可靠性。在此基础上,将部分刚性壁面替换为柔性壁面,结果表明,柔性壁面可以抑制T-S波在空间上的增长,从而推迟边界层流动转捩。壁面的变形不只跟随T-S波的波形,还因为柔性段与刚性段相接的前缘和后缘引起与扰动源频率相同的更大尺度的壁面波动,壁面的实际变形由这几种波叠加而成。开展的参数研究结果表明,增大表面的质量密度对于柔性壁面衰减扰动的效果几乎没有影响;增大表面张力和增加底部支撑的弹性系数可以增加壁面的刚性,减小壁面变形的幅度;增加阻尼可以抑制柔性段前后缘产生的大尺度壁面波动的传播,而对跟随T-S波的变形影响不大。总体上,柔性壁面的变形程度越大,其扰动的抑制效果越强。      

耀斑激波传播的数值研究

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

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

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

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

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

激波风洞7°尖锥边界层转捩实验研究

高超声速边界层转捩对摩阻、传热等有重要影响，飞行器的研制迫切希望能精确预测和控制边界层转捩。在JF8A激波风洞中开展了7°半锥角的高超声速尖锥边界层转捩实验研究，利用响应频率达到1 MHz量级的高频压力传感器对尖锥壁面脉动压力进行了测量，并结合热流测量结果，研究了高超声速尖锥边界层中扰动波的发展过程。实验结果表明:JF8A激波风洞在雷诺数为6.4×106/m状态下核心流的自由流噪声为2.8%；高频脉动压力测量技术能清晰地捕捉转捩过程中的第二模态波及其发展历程，试验状态下模型的第二模态波频率范围为165~206 kHz。当前研究结果能够为高超声速数值方法验证提供数据支撑。      

具有不同等离子体β值的Hall MHD重联中的低频波研究

数值研究具有不同等离子体β值(等离子体热压与磁压之比)的Hall MHD重联过程.对于β值在0.5≤β≤6.5的算例,稳态重联率(6)A/(6)t|st为0.15≥(6)A/(6)t|st≥0.095.本文研究β值在上述范围内的Hall MHD重联中的波动.通过快速傅里叶变换(FFT),将平面外磁场By分量和xz平面内的速度分量vx,vz位于给定点的时间序列转换为功率谱.结果表明,By,vx,vz功率谱的频率范围为Ωci＜ω＜8Ωci,其中Ωci为离子回旋频率.随着β值增大,功率谱峰值能量降低,并且朝着低频端移动.对于电场E采用最小方差分析法(MVA)可以确定波的传播方向k,波矢k与局地磁场B之间的夹角α随β值增大而增大.对于β值最大的算例3(β=6.5),α＜28°.上述结果表明,k准平行于B.为了研究波的偏振特征,在垂直于k的平面内画出E的矢端图.三个算例的矢端图均表明,波动具有右旋偏振特征,在算例1(β=0.5)中矢端图为右旋椭圆偏振,而在算例2(β=2.5)和算例3(β=6.5)中矢端图变成复杂的右旋偏振图形,这些均为哨声波的典型特征.本文的研究为快速磁重联与哨声波之间提供了一种可能的联系.      

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

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