凸拐角附近激波与湍流边界层干扰的数值模拟研究

采用大涡模拟方法数值研究了来流马赫数Ma∞=3、雷诺数Reθ=2070的凸拐角附近激波与湍流边界层干扰问题。对于尖拐角和钝拐角两种情况,计算考察了激波入射到拐角下游时的流场结构、非定常分离流以及湍流统计特性。研究发现:尖拐角能够使流动产生局部集中的顺压梯度,导致壁面摩擦系数在拐角处出现峰值,因此分离点难以越过拐角向上游运动,而再附点发生前后运动。与尖拐角情况不同的是,钝拐角产生的局部顺压梯度不足以阻止分离点越过拐角,因此分离点和再附点都会发生运动。此外,尽管凸拐角产生的顺压梯度能抑制分离点附近的近壁小尺度涡结构,但在经过入射激波作用之后,近壁小尺度涡结构在再附点附近显著增强。     

急流与低层大气重力波能量的相关性研究

通过分析武汉(30.5oN, 114oE)上空2000至2004年的Radiosonde常规观测数据, 对当地对流层(1~10 km)与低平流层(18~25km)的某些惯性重力波特征进行了研究, 发现重力波能量与急流强度变化之间存在显著的相关. 同时分析了海口(20oN, 114oE)(2000至2004年)、北京(40oN, 116oE)(2001年12月至2003年2 月)的Radiosonde观测数据并进行了比较, 发现较高纬度地区的相关性明显大于低纬度地区. 通过对武汉地区2006年1月5天Radiosonde加密观测数据的分析, 讨论了重力波动能与势能的高度变化, 进一步发现纬向风的垂直剪切与急流中心的分布分别与重力波势能和动能的强弱分布相对应.      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

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

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

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

多循环脉冲爆轰发动机外流场的实验研究

利用内径为30mm,长度为1000mm的爆轰管,对多循环脉冲爆轰发动机外流场进行研究.爆轰管内充满氢气-氧气-氮气预混气,采用底端中心点火.实验中采用YA-16高速阴影系统,拍摄了爆轰外流场的阴影照片.根据实验结果,分析了外流场的动力学结构和悬吊激波产生的动力学机理,讨论了两次循环管外流场结构间的差异.第一次循环中的爆轰波溢出爆轰管后,斜压效应和Helmholtz不稳定导致涡环的产生,涡环形状的变化又引起悬吊激波的角度变化,在此过程中引导激波始终为球面形状,在第二个循环中的相应时刻,激波前的气体是向外流动的,导致引导激波逐渐成为葫芦形状.     

吸气式脉冲爆震发动机原理性试验研究

设计了内径分别为60、70mm的两组吸气式两相脉冲爆震发动机,在略高于常压状态下成功进行了以汽油为燃料、以空气为氧化剂的吸气式两相脉冲爆震发动机原理性试验.进气道内采用无阀的进气结构,试验中在爆震管内部安装Shchelkin螺旋来促进爆震波的生成,所测量的爆震波接近充分发展的C-J爆震波.内径60、70mm的发动机最高工作频率分别为15、20Hz.进气道内压强变化幅度在0.1MPa左右,说明所设计的进气道与爆震管匹配良好,这将为脉冲爆震发动机用于工程实际提供技术储备.     

气动阀型式对脉冲爆震发动机爆震特性的影响

通过改变气动阀结构、堵塞比及进气阻力系数,研究了其对脉冲爆震发动机(PDE)爆震波压力特性的影响.为了改善燃油雾化、蒸发和掺混,低充填速度PDE宜采用双旋流式气动阀.同一类型气动阀堵塞比大,爆震效果好,但会使PDE充填速度和工作频率降低,其最佳堵塞比为60%～70%.不同型式气动阀爆震燃烧效果除和堵塞比有关外,还和正反向流动阻力系数有关,阻力系数越大,爆震燃烧效果越好.结果表明,结构简单的气动阀具有单向阀和气动雾化喷嘴的功能,能够用于PDE.     

等离子体覆盖目标散射特性的RKETD-FDTD分析

采用一种既保证计算的高效率,又有较高的计算精度的龙格库塔指数时程差分时域有限差分法(RKETD-FDTD)研究了等离子体的散射特性.该算法解决了电磁波在色散介质中传播的计算问题,导出了在等离子体介质中RKETD-FDTD迭代公式.文中分别计算了等离子体平板的反射系数和非均匀等离子体覆盖导体柱的散射特性,所得结果与解析结果相符合,并且表明等离子体涂层选择合适的碰撞频率,能有效地减小目标的雷达散射截面(RCS).     

共振相位调制型光纤SPR传感器理论分析与数值模拟

提出了一种基于表面等离子体渡共振(Surface plasmon resonance,SPR)相位调制的新型光纤传感模型.根据光波导理论和多层反射理论,构建了涉及多个共振激励参量的理论分析模型,并研究了共振相位调制型SPR效应激励机理和特性.模拟结果表明,入射光波长、纤芯折射率和金属膜介电常数实部绝对值增大,均会导致相位跃变所对应的液体折射率检测范围向高折射率方向发生偏移;而入射光角度和金属膜虚部逐渐增大,将会引起所对应的相位调制灵敏度逐渐退化.因此通过选择共振激励参量,可以实现对折射率检测范围和灵敏度的有效调节.