波浪形非均匀间隙封严结构影响涡轮性能的数值模拟

利用数值模拟的方法研究了波浪形非均匀间隙封严结构和均匀轴向间隙封严结构下轮缘封严气流对涡轮性能的影响.研究表明:燃气入侵与出流结构受到静盘、动盘及主流切向速度的影响,以低于动盘转速同向旋转,并改变了转子的进气条件,增强了压力面马蹄涡强度,因此对转子出口流场造成很大影响.封严气流与上游导叶尾迹的相互作用引起转子通道内熵增,造成涡轮效率的下降.与均匀轴向间隙封严结构相比,波浪形非均匀间隙封严结构使大的入侵与出流结构破碎为小的结构,对涡轮性能的负面影响减小,涡轮效率提高了0.9%.结果证明了波浪形非均匀间隙封严结构在具有较好的封严效果的同时提高了涡轮性能.      

转静干涉对跨声速吸附式压气机流场的影响分析

跨声速低反动度吸附式压气机的静叶流道中不但有附面层抽吸,还会存在由于低反动度设计而产生的激波,因而流动现象异于常规压气机,非定常效应明显,故采用非定常数值模拟方法,对跨声速低反动度吸附式压气机的内流场展开研究。通过选取两个最具代表性的时刻,探究了转静子干涉对该压气机气动性能以及内流场的影响。结果表明:非定常不同时刻的效率峰谷值之差达0.937%,流动的非定常性很强。转子尾迹对静子流道中展向抽吸缝前的区域干涉作用明显。在静叶流道中,分离主要存在于三个位置:30%弦长位置、尾缘顶部以及尾缘根部,其中静叶吸力面侧30%弦长位置的分离原因是激波与附面层干扰,且该处部分低能流体的吸除方式为"螺旋路径吸除",并详细分析了其对流场的积极和消极影响。此外还建立了激波尾迹干涉——激波与附面层干扰——非定常效率波动这三者之间的联系,找出了非定常效率波动的原因。     

硅微陀螺平板结构的压膜阻尼研究

压膜阻尼对硅微陀螺的精度具有重要影响。平板结构是硅微陀螺中一种典型结构,本文通过对雷诺方程的推导,得出平板在正弦规律小振幅刚性运动时的阻尼力解析式,并给出低频时阻尼系数的近似算法。利用Matlab和CoventorWare软件对运动平板的阻尼进行仿真计算,得出运动平板的阻尼系数和刚度系数,并分析平板参数对阻尼系数及刚度系数的影响。建立不同类型的带孔平板模型,通过比较可以看出,当阻尼孔按一定方式排列时,可以降低结构的阻尼。     

目标函数误差估算及网格自适应处理

发展了一种基于目标函数误差估算的网格自适应准则,进而通过网格自适应处理提高目标函数计算的效率和准确性。首先描述了目标函数的误差估算及修正的方法,该方法通过伴随方程将原方程的残值误差与目标函数联系起来。然后,基于误差估算建立网格自适应准则,以减少目标函数修正后的剩余误差,提高目标函数计算结果的准确性,并将方法进一步发展至多目标问题。最后将该准则应用于Euler方程控制的NACA0012翼型无粘可压流动的网格自适应模拟。数值计算成功地捕获了与升力、阻力和力矩等目标函数相关的特征流动区域,计算出符合指定精度要求的目标函数值,验证了本文所发展的方法。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

应用协同射流控制的临近空间螺旋桨高增效方法

基于雷诺平均Navier-Stokes方程与多块搭接网格技术,数值模拟了协同射流(CFJ)翼型及桨叶的黏性绕流,分析了CFJ技术的增升减阻效果及工作机理。研究了CFJ的功率及效能比的分析方法,量化分析了CFJ的能量利用率。开展了喷口大小、喷口动量系数等参数对CFJ翼型性能的影响规律研究,并在此基础上开展了应用CFJ的临近空间螺旋桨高增效方法研究。结果表明:数值模拟结果与实验值吻合良好,在不同状态下,CFJ控制技术均能显著改善翼型气动性能。其中,最大升力系数提高了60%~130%;阻力系数降低了100%~440%,部分小迎角工况甚至出现负阻力系数,升阻比显著提高;翼型失速特性明显改善,失速迎角提高了近10°;能量利用率高,效能比可达440%。最终,在最优参数条件设置下,采用基于CFJ控制技术的临近空间螺旋桨可提高效率5%以上。     

基于Kane方程的拖曳式诱饵释放过程动态特性分析

研究了航空拖曳式诱饵释放过程中的动态特性。根据诱饵的运动状态,将释放过程划分为自由状态、放索状态和拖曳状态。利用Kane方程建立了拖曳式诱饵释放的多体系统动力学模型,其中拖索离散为若干段刚性杆,诱饵视为刚体与拖索铰接,其上的作用力包括铰约束力、气动力和重力。针对放索过程中第1个索段质量时变引起的变质量动力学问题,采用Generalized-α算法进行时域求解。在此基础上分析了载机不同飞行高度、飞行速度、拖索释放速度以及拖曳点位置对诱饵释放过程中动态特性的影响。结果表明:在低空高速下释放,诱饵的俯仰角幅值较小,质心相对位置变化较稳定,收敛速度较快,但载机飞行速度过大时,诱饵容易靠近载机尾流区并受其影响;放索速度增大时,诱饵俯仰角幅值增大,质心会出现纵向沉浮运动;拖曳点远离重心时,俯仰角震荡幅值增大,当靠近重心时,收敛性变差,应合理设计并优化拖曳点位置和放索速度。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

基于能量法和特征值法的颤振预测数值方法研究

为研究适用于工程设计的颤振数值预测方法,采用能量法和特征值法对出现颤振现象的真实叶片进行预测。通过将2种数值方法的预测结果与试验结果进行对比,验证数值预测方法的精度和工程适用性。结果表明:2种方法都较准确的判断了某典型风扇转子叶片的气弹稳定性,与试验结果基本吻合,具有较高的工程实用价值。     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。