基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

大气紊流环境下的纵向飞行参数响应分析

介绍了大气紊流环境下飞机纵向飞行参数的计算与匹配方法.由于紊流的随机性特征,飞行参数响应分析采用时域与频域相结合的方式,对某型飞机紊流环境下的纵向飞行参数进行了计算分析.采用所介绍的方法,可以获得比较合理的飞行情况来指导型号设计工作.     

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

尾缘锯齿结构的降噪物理机制实验

对比分析了常规尾缘翼型与锯齿尾缘翼型尾缘湍流流场的基本特征,并通过线阵列的方法测量了两种尾缘结构的噪声.结果表明:锯齿尾缘翼型尾缘湍流流场的湍流强度以及3个方向上的湍流强度都相比于常规尾缘翼型有显著减少,声场结果显示锯齿尾缘翼型对尾缘噪声有显著减小,对前缘噪声影响很小.锯齿结构加宽了尾迹区域并加快了大涡的破碎,产生了额外的马蹄涡,湍流脉动衰减率沿着流动方向变大.      

吸气式高超声速飞行器机体推进一体化技术研究进展

吸气式高超声速一体化飞行器最显著的特点是子系统之间的耦合较其他类型飞行器更加强烈,这使得其设计具有挑战性。所有的子系统之间部件相互干涉,包括:气动、推进、控制、结构、装载和热防护等,特别是机体与超燃冲压发动机之间的耦合最为突出。飞行器的前体和后体下壁面既是主要的气动型面,又是超燃冲压发动机进气道外压缩型面和尾喷管的膨胀型面,在产生推力的同时也产生升力和俯仰力矩。机体与发动机的强耦合作用对飞行器的推力、升力、阻力、俯仰力矩、气动加热、机身冷却、稳定性和控制特性有直接的影响。本文介绍了国内外机体推进一体化技术的研究进展,重点介绍了中国空气动力研究与发展中心(CARDC)的相关研究工作,包括:密切曲锥曲面乘波进气道和基于双激波轴对称基准流场内转式进气道设计方法、独创的大尺度脉冲式燃烧加热风洞一体化飞行器带动力试验技术和高超声速内外流耦合数值模拟技术等。对高速飞行中激波边界层相互干扰、流动分离机理、可压缩湍流转捩及其控制、超燃冲压发动机燃烧流动机理等相关基础问题也进行了研究,强调了对高效高精度计算方法的迫切需求。     

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

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

基于自然激励技术的颤振边界预测

为了预测紊流激励条件下机翼的颤振边界,基于自然激励技术提取紊流响应的自由衰减信号,采用矩阵束方法识别模态参数,最后通过Z-W(Zimmerman-Weissenburger)方法计算稳定性判据,拟合判据变化曲线并外推颤振边界.对平板机翼模型进行了数值仿真分析,对单独机翼模型风洞颤振试验数据进行了计算.结果表明:采用自然激励技术与矩阵束方法能够较准确地识别紊流激励响应的模态参数,频率识别误差小于6％,阻尼比识别误差小于30％,结合Z-W方法能够在较低风速较早地预测颤振边界,有助于提高试验的安全性.     

燃油喷嘴雾化特性试验研究

为了检验某型航空发动机燃油喷嘴改进设计效果,利用相位多普勒粒子分析仪对燃油喷嘴的雾化性能参数进行试验研究。得到雾化液滴的索太尔平均直径的空间分布、轴向平均速度、脉动速度及其湍流度的分布情况。结果表明:轴向平均速度呈凹盆状分布,脉动速度呈双峰状分布;喷雾中心湍流度大,喷雾边缘湍流度小。随着供油压力增大,在相同测试截面上,喷雾的范围和中心区域粒径变大,边缘位置粒径变小。在相同供油压力下,随着与喷嘴距离的增加,喷雾范围增大,喷雾的轴向平均速度和脉动速度减小,轴向速度的湍流度波动幅度减小。     

高湍流度时全气膜涡轮叶片表面冷却和换热特性的实验研究

为获得高主流湍流度时全气膜涡轮叶片表面的冷却和换热特性,在跨声速风洞中实验研究了质量流量比(MFR)和主流雷诺数(Re)对叶片表面气膜冷却效率和换热系数比的影响。在叶片前缘布置了5排圆形孔,在吸力面和压力面分别布置了3排和6排圆形孔,实验结果由嵌入在叶片中截面的热电偶测得。实验中基于弦长的主流雷诺数的范围为3.0×105～9.0×105,叶栅出口马赫数Ma为0.8, MFR的范围是5.5%～12.5%,主流湍流度Tu为14.7%。实验结果表明:主流雷诺数升高显著增强了叶片表面的换热,使层流边界层到湍流边界层的转捩位置提前。对于吸力面S/C0.2的区域(S/C为当地弧长与弦长之比),气膜冷却效率受MFR影响明显,当MFR大于7.7%时提高MFR会导致气膜冷却效率降低;该区域的换热系数比在中低雷诺数时受MFR影响较小,在高雷诺数时随MFR升高而升高。压力面S/C-0.7区域的气膜冷却效率随MFR升高而升高,-0.7S/C-0.4区域的气膜冷却效率受MFR影响较小,对于整个压力面而言,MFR升高提高了叶片表面的换热系数。相对于叶片其它区域,压力面后半段区域和吸力面的气膜冷却效率受雷诺数影响较大。