边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。     

动静圆柱坐标系下N-S方程的分区算法及子域干涉处理研究

耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程,统一了动子域和静子域的控制方程;在子域计算中,采用高分辨率的三阶ENN格式及LU-SGS隐式解法;在子域干涉处理中,研究并对比了动静子域干涉面的重合处理与重叠处理;最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气机叶栅内三维紊流流场的数值模拟,验证了该分区算法的有效性,计算结果与实验数据吻合较好。     

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

解决欧拉方程奇异性的方法探讨

提出了一种解决欧拉方程奇异性的工程算法,从飞行器运动方程入手,推导出具体计算公式,并进行了数值计算验证。该方法无论是理论分析还是工程计算都证明是可行的,没有发现原理缺陷和方法误差,不会产生奇异性,在飞行器仿真、模拟和姿态控制中具有普遍的应用和参考价值。     

虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

再入弹抛壳气动特性研究

本文用有限差分方法求解再入母弹头与其所抛壳体之间的无粘干扰流场,获得了体—壳的干扰气动力,并在准定常假设的基础上近似模拟了弹体与壳体的分离过程。数值研究的结果表明,体—壳间存在的干扰会导致壳的压心位置急剧变化,危及母弹安全;抛壳速度是体—壳安全分离的关键参数。     

基于FW-H方程的旋翼气动声学计算研究

由流体力学N S方程导出的非齐次波动方程———FfowcsWilliams Hawkings方程(简称FW H方程),可以精确地描述在静止流体中运动的物体与流体相互作用的发声问题。以FW H方程为理论模型,将旋翼桨叶运动发声问题等效为包含桨叶的任意运动控制面(声源面)的声辐射问题,并在旋翼绕流Euler方程数值模拟的基础上,在时域内计算了悬停旋翼和前飞旋翼的声场。应用于UH 1H和AH 1/OLS两种旋翼模型的气动声学计算表明:计算结果与噪声实验值符合良好;所研制的程序不仅能够较准确地计算单极子噪声和偶极子噪声,而且具有较强的跨音速四极子噪声预测能力。     

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

跨声速粘性流绕振荡翼型的非定常计算

本文采用时间和穴是二阶精度的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式及ｑ－ω二方程微分模式，结合网格自适应技术，数值模拟Ｎ－Ｓ方程，计算了跨声速下的翼型非定常运动，包括俯仰，浮沉和前后平移振荡。结果表明，压力分布和气动系数与实验基本符合，微分模式和自适应网格能够显著提高激波和边界层计算精度。