超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

高超声速飞行器后体/尾喷管优化设计

结合试验设计方法、替代模型技术和遗传算法构建一套改进的优化方法,并将其应用于高超声速飞行器后体/尾喷管的一体化设计.在构建优化方法时,针对替代模型采用渐近全局策略提高精度;针对遗传算法,采用实数编码、多目标定级排序和改进小生境技术.在后体/尾喷管的一体化设计应用中,结合高精度的计算流体力学(CFD,Computational Fluid Dynamics)求解,以两个设计点的推力和升力为目标,以力矩为约束,得到优化问题的Pareto最优前沿面,优化结果在综合性能方面有很大提升.该优化方法可进一步推广应用于更复杂的优化设计当中.      

洞映射方法的研究和改进

针对传统洞映射方法存储大的缺点,对"广义封闭"的概念进行扩展,提出了最小洞映射方法,该方法允许挖洞曲面结束于网格截断面,有效缩小了洞映射区域,节省了存储开销.对适用于广义封闭的洞映射单元识别方法进行了分析,指出了用物面信息判断映射单元属性可靠性较差.发展了一种新的特别适合广义封闭问题的识别方法:"Inverse mark",使用计算网格结点信息自动识别洞外单元,再作为初始点在网格内部推进.研究表明:"Inverse mark"方法计算效率高,可靠性好,自动化程度高,有效提高了重叠网格方法对缝隙等局部复杂结构的适应性     

火箭整流罩外气动噪声环境的大涡模拟研究

基于五阶加权本质无振荡（WENO）格式构造隐式大涡模拟方法（ILES）,对跨声速来流条件下（Ma=0.8）火箭整流罩外噪声环境进行数值模拟。通过与风洞试验结果及国外文献进行对比,ILES方法能够在较粗网格下准确预测壁面湍流脉动特性。跨声速流动在壁面折角处出现分离、再附、激波/边界层干扰现象,均方根脉动压力系数出现峰值,同时该区域噪声能量在全频段都较高,易引起结构抖振效应。最后,根据ILES模拟结果,指出工程常用的外噪声经验公式的不足,并提出改进措施。     

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

吸气式高超声速飞行器多参数灵敏度分析

为实现吸气式高超声速飞行器多参数情况下的灵敏度分析及参数分类,降低设计的复杂度,在吸气式高超声速飞行器参数化建模的基础上,首先采用正交试验设计生成样本,再通过计算流体力学(CFD)进行高精度气动力性能计算,最后运用方差分析法进行气动性能的参数灵敏度分析,在运用较少的试验样本点的情况下,即可完成多参数、复杂构型条件下气动性能的参数灵敏度分析。结果表明,该方法可以正确地分析出参数对飞行器气动性能的敏感程度,得到参数对气动性能的影响规律。同时,通过灵敏度分析的计算样本,还可以初步选出气动性能较优的飞行器构型,为地面试验和优化设计提供参考。     

鸭式布局地效翼数值研究

针对地效飞行器气动特性和纵向静稳定性这两个重要问题,应用数值模拟的方法对鸭式布局地效飞行器进行了研究。在鸭式布局中,主翼和鸭翼互相干扰,流场比较复杂。简化了地效飞行器模型,只考虑了主翼和鸭翼,通过改变鸭式布局地效飞行器中鸭翼相对主翼的位置和角度,分析其对气动特性和纵向静稳定性的影响。     

网格之间流场信息传递的高精度方法和验证

理论分析了插值方法引入误差的原因,提出一种利用时间推进方法和移动网格技术实现网格之间流场信息传递的思想,并构造计算方法。在介绍新的信息传递算法原理的基础上,根据三维情况下传递过程需要满足的时间推进步长稳定性条件,推导出最小推进步数的要求。最后,通过对激波诱导物体运动的非定常流动模拟,比较了信息传递算法和线性插值、加权最小二乘插值方法,计算结果表明,与插值方法相比,信息传递过程没有引入误差,能有效抑制非物理振荡,保持流场计算精度。     

利用CFD技术计算飞行器动导数

利用非定常流场数值计算方法模拟飞行器强迫俯仰振荡仅能得到俯仰力矩系数对迎角变化率和俯仰角速度的动导数之和,而动稳定性分析需要单独的动导数数值.为解决这个问题,利用滑移网格模拟强迫振荡运动,得到俯仰力矩系数对迎角变化率和俯仰角速度的动导数之和.利用旋转参考坐标系模拟定常拉升,得到俯仰力矩系数对俯仰角速度的动导数.利用旋转参考坐标系模拟匀速滚转,得到滚转力矩系数对滚转角速度的动导数.对有翼导弹和水上飞机进行了纵向和横向动导数的计算.计算结果与试验数据、文献数据以及其他方法得到的结果具有较好的一致性,表明提出的方法可用于复杂外形飞行器动导数计算.