基于修正边界层转捩的翼型气动特性研究

边界层转捩是决定翼型流场特性的重要因素,因此在航空工程中开展边界层转捩研究具有重要意义.从转捩流动的物理特征出发,引人间歇函数对k-ω SST两方程湍流模型的Wilcox转捩模式进行了修正,进而对传统的NACA0012翼型的流场特性和气动性能进行了研究,并与原始边界层转捩的计算结果以及试验结果进行了比较.研究表明,改进后的模型对转捩位置具有较好的预测能力;在采用修正边界层转捩模型的情况下,翼型的阻力预测精度有了一定程度的提高.     

超燃进气道分离流动的凸起物控制机理研究

双模态超燃冲压发动机由于压力扰动可能发生不起动现象,造成推力严重下降,对飞行稳定性与飞行安全具有很强的破坏性。不起动初始阶段主要受到激波与边界层相互作用引起的流动分离影响,采用5阶特征型WE-NO格式与3阶TVD型Runge-Kutta格式的高精度数值方法,求解三维Navier-Stokes方程,研究了利用凸起物作为涡发生器的被动控制方法,及其对高超声速流动分离现象的控制效果。结果表明高精度数值格式能够捕捉到复杂精细的流动分离结构,总结了摩阻、压力等在分离再附位置的变化规律;发现凸起物通过诱导形成局部流向涡进行流动控制,能够改变压力分布,减弱分离强度,影响分离结构。     

振动尖塔对风洞模拟大气湍流边界层的作用

针对大气边界层风洞模拟中常常出现的湍流度随高度衰减太快的问题，本文尝试了一种新的被动模拟方法——振动尖塔法。该方法主要采用了具有弹性底座的尖塔型旋涡发生器。尖塔群受风洞自由来流驱动，发生随机振动，向下游流场注入了强烈的低频扰动，致使试验区的湍流度、积分尺度和风谱惯性子区的宽度都明显增加，改善了风洞的模拟性能。     

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

氢气泡法图像的数字化处理分析

对氢气泡法产生的时间线进行数字化，把图像的灰度分布存入数组进行运算，识别氢气泡时间线的前缘和后缘，根据距氢气泡发生线最近的两条时间线的前缘或后缘间距以及氢气泡时间曲线的时间间隔，计算流速分布。本文还提出一种测量湍流结构的速度的方法，简称结构相关法；对两幅一定时间间隔的氢气泡图像进行数字化得到灰度数组，圈定拟序结构的一部分，利用相关原理比较两图像的数组，得到所圈定拟序结构的位移，进而求得其运动速度。     

两方程湍流模型的可压缩性修正及其应用

本文在不可压两方程湍流模型的基础上进行可压缩性修正,给出了ｋ－ε、ｋ－ε两方程模型的可压缩性修正方程,并对亚声平板湍流流动和跨声喷管流动进行了数值模拟,得到了较好的结果。从而证实了两方程湍流模型经过可压缩性修正,能较好地模拟可压缩流动。     

SFD方法模拟高超声速进气道中激波与边界层干扰

采用ＳＦＤ方法,在ＥＸＣＥＬ软件上利用ＮＮＤ格式数值模拟了激波与平板边界层之间的相互作用以及激波与管道壁面边界层的相互作用,即高超声速进气道隔离段内激波串流动现象。该方法的特点是简单、直观而且不用常规编程,计算结果基本反映了流场的基本结构,表明ＳＦＤ方法是一种具有潜力的计算方法。     

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

二维过渡区喷流干扰流场的DSMC/EPS仿真

应用DSMC／EPS混合算法，数值模拟过渡区有RCS横向喷流干扰的二维有限平板高超声速流场，分析了流场结构、分离形态及物面气动系数分布等特性。结果表明，不同：Knudsen数下，喷流干扰流场存在差异，过渡区中的稀薄气体效应得以体现，证明了DSMC／EPS混合算法是研究稀薄，连续混合区域高超声速流动的有效方法。     

高超声速圆锥边界层转捩数值研究

以类似于处理湍流脉动的方式来模型化非湍流脉动,考虑湍流的间歇性特征及引起超声速流动失稳机制,研究包括转捩起始位置及转捩区范围的高超声速流动转捩问题,发展可反映扰动模态和可压缩性影响的两方程SST转捩/湍流模式.针对尖锥模型,在马赫数6、攻角O°～2°时,采用第二模态(Second-mode)机制,预测并比较完全层流态、完全湍流态,以及自然转捩发展过程中的表面热流、表面摩阻等气动特性.