考虑鲁棒性的超临界翼型激波控制鼓包减阻研究

本文主要采用CFD方法研究超临界翼型的激波控制鼓包减阻技术,分析了鼓包形状参数对减阻效果和鲁棒性的影响规律.计算结果显示,鼓包的减阻效果受位置影响较大,当鼓包最高点与干净翼型的激波位置相同时减阻效果较好,鼓包高度过高对减阻效果不利,而较长的鼓包可在更大的高度范围内实现减阻.鼓包还可以通过弱化激波,抑制附面层分离,延缓超临界翼型抖振现象的发生.计算结果显示,鼓包减阻技术整体而言工作范围较窄.但经过设计,较长且较低的鼓包可以在较大的升力范围内具有减阻效果,并且减阻效果对形状变化及雷诺数变化不敏感,还能有效提高阻力发散马赫数,鲁棒性要明显优于较短较高的鼓包,具有工程应用的潜力.     

跨音速透平扭叶片的气动优化设计研究

以并行自适应差分进化算法为核心，耦合曲面造型方法以及计算流体动力学求解技术，发展了一种适用于叶轮机械三维气动优化设计的全局自动气动优化算法．利用该算法，以等熵效率最高为目标，在满足流量约束的条件下对跨音速扭叶片进行了气动优化设计．对优化结果的详细分析表明，最优叶栅的等熵效率比原始叶栅提高了1．1％，气动性能有显著的改善，算法具有良好的优化性能．在跨音速条件下，载荷分布对叶栅的气动性能有着巨大的影响，采用前加载设计可有效地减弱斜激波的强度，减少激波损失，提高流动效率．因此，通过优化叶栅型线来改变叶栅的载荷分布可有效地提高叶栅的气动性能。     

电弧放电诱导喷气流偏转数值模拟研究

利用电弧放电，在拉瓦尔喷管的扩散段内诱导出斜激波，使得扩散段内的超声速气流经过该斜激波后发生偏转喷出喷口，从而产生推力矢量。为获得不同条件对电弧放电诱导喷气流偏转效果的影响规律，在电弧放电与超声速流场耦合的数学模型基础上，开展喷气流马赫数1.5的条件下，电弧放电区域大小、放电温度、放电位置及喷管落压比对电弧放电诱导喷气流偏转效果影响规律的数值模拟研究，并通过初步的实验验证。实验结果表明电弧等离子体诱导激波产生喷气流偏转的方案可行，数学模型建立正确。     

运动激波和气泡(串)相互作用的三维数值模拟

对激波和气泡（串）相互作用诱导的大变形界面演化进行了三维数值研究.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level Set方程追踪多流体界面.在三维情况下,采用GFM确定界面边界条件.将速度分量分别外推,再沿法线分解,得到伪流体切向速度分量,给出了不同时刻运动激波和气泡（串）作用后的压力、密度云图分布图像.计算结果表明：采用速度分量外推方法适合三维界面边界条件处理,并与二维情况下切向速度直接外推的计算结果兼容,验证了切向速度外推方法的正确性.该方法模拟多个流体界面也是有效的,得到了高分辨率的三维变形界面.     

修正变光滑长度SPH方法及其应用

提出了修正的变光滑长度SPH方法及其算法实现.与传统变光滑长度SPH法相比,新方程组基于对称形式核函数近似,并从根本上考虑了变光滑长度影响.方程组中,密度演化方程与变光滑长度方程隐式关联;而SPH动量方程和能量方程在Springel提出的全守恒SPH方程组基础上,改进其分散核近似形式为对称核近似形式得到.由于新方程组密度演化方程和变光滑长度方程隐式耦合,采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程的办法,迭代过程增加的计算量相对很少.给出了2个一维激波管算例和二维 Sedov算例来验证方法的有效性.结果表明,新方法很好地解决了传统SPH法中变光滑长度影响,特别是在模拟二维 Sedov问题时新方法能得到比Springel方法更准确的压强峰值位置,中心压强也更接近理论值.特别适合模拟爆炸与冲击、大变形大扭曲等密度和光滑长度变化剧烈的问题.     

双S弯进气道锤激波动态特性研究

发动机喘振产生的瞬时高压引起的进气道内锤激波载荷是飞机进气道结构强度设计的关键因素之一。通过对飞行器进行内外流一体化非定常三维数值计算,分析双S弯进气道内锤激波三维流场非定常特性的演化过程,研究了不同马赫数,不同超压比对双S弯进气道锤激波载荷的影响。研究表明：锤激波经过进气道弯道时,弯道外侧流场压力远大于内侧;锤激波离开进气道入口后,进气道内流场经过数个周期逐渐恢复至初始状态;相同进气道反压时,来流马赫数越小,锤激波在进气道内部传播速度越快,且进气道内部压力系数峰值越大;相同来流马赫数下,随着超压比的增大,锤激波在进气道内部传播速度加快,进气道内部压力系数峰值增大。     

激波诱导甲烷点火化学反应区的可视化实验研究

研究激波诱导甲烷/空气混合物点燃的反应区特征, 对于认识其气相爆轰规律有着重要意义. 在激波管平台上, 进行了当量比甲烷/空气混合物的激波点火实验, 并采用化学发光成像和平面激光诱导荧光(PLIF, planar laser induced fluorescence)等方法得到了不同工况下的反应区特征. 结果表明, 由于在弱点燃工况下诱导区的着火延迟期沿激波传播方向的梯度较小, 给非线性化学反应过程提供了更多的时间, 反应特征出现显著的不均匀性; 而随着激波的增强, 这种不均匀性显著减小, 在强点燃工况反应区呈规则的爆轰波结构. 单脉冲PLIF方法测量的OH基相对分布揭示的反应区特征与文中其他方法的测量结果相一致, 并和前人的相关研究结论相符, 为激波管内激波和气相物质相互作用条件下进行高时空分辨率测量提供了新的方法和思路.     

耀斑激波的各向异性传播

对飞船观测和行星际闪烁记录的分析表明,耀斑激波的最快传播方向偏离耀斑法线:在赤道面内朝东偏转,在子午面内偏向日球电流片.即便扣除背景太阳风的对流效应,上述非对称传播特征仍然存在.可是,迄今对日球赤道面内耀斑激波的数值模拟给出的最快激波传播方向几乎与耀斑法线平行,与观测结果相左.本文将证明,适当修改下边界的激波引入方式,可以获得最快传播方向东偏之结果. 既然耀斑源区随太阳自转,采用共转坐标系比较方便.在该坐标系中,日球赤道面内的理想 MHD二维平面流动满足如下方程: