可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

立楔诱导的高超声速流动分离数值分析

高超声速飞行器进气道等关键部件引起的激波与边界层相互作用将导致流动分离,从而改变当地压力分布与局部受热情况,影响飞行稳定性与飞行安全,因此需要对高超声速流动的分离现象进行细致研究。采用高精度5阶特征型WENO格式与3阶TVD型Runge-Kutta方法,求解三维Navier-Stokes方程,对立楔诱导的高超声速激波与边界层相互作用引起的分离流动流场结构进行了细致的数值模拟与分析。结果表明,5阶特征型WENO格式分辨率远高于类TVD格式;Ma=6时得到清晰的激波结构、分离涡结构及其演化过程和壁面极限流线的拓扑结构,证明了WENO格式应用于高超声速分离流动的可行性与高分辨率;对不同来流Mach数的对比证明Mach数的增大抑制流动分离,导致分离涡减小。     

洛仑兹力控制高超声速进气道边界层分离的数值模拟

采用空间HUE格式、时间LU-SGS推进、sst-kw湍流模型、多块结构网格程序,对磁流体动力学(Magnetohydrodynamic:MHD)控制高超声速二维进气道边界层分离进行了数值研究.研究发现,不施加控制时,数值模拟得到的壁面静压和实验结果符合良好,进气道喉道处分离区占据喉道高度的1/3左右.通过施加MHD控制,消除了进气道内部的边界层分离,总压恢复系数从0.502提高到0.56,喉道处流场畸变系数减小了18.6%.     

翼型表面结冰准定常数值模拟

对翼型表面的结冰过程进行了准定常数值模拟。在每个时间步长内完成网格随着壁面边界的移动而更新、周围流场和水滴撞击特性重新计算、冰形计算及壁面边界的重构工作,如此循环直至所需的结冰计算时间。采用拉格朗日轨迹追踪法获得水滴撞击特性,应用考虑壁面粗糙度影响的边界层积分法计算壁面的对流换热系数,在此基础上求解Messinger结冰热力学模型,冰层始终沿壁面外法线方向生长。对NACA0012翼型在不同环境下生成的三种典型冰形进行了预测,并与实验结果进行了比较,表明本文所述方法的有效性。     

"回收/调节"方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离...     

等离子体控制边界层法流动仿真研究

使用泊松方程、漂移-扩散模型和N-S方程对介质阻隔面放电等离子体控制边界层流动进行了一体化数值模拟.仿真结果表明介质阻隔面放电等离子体通过推-拉空气可以对低速气流进行控制,正向推力和反向拉力均可以在距壁面一定高度处形成射流,正向射流自身增加边界层动量,更重要的是通过引射主流将主流中的能量输运到边界层中,从而减小边界层厚度;反向射流与来流相撞后形成回流区将边界层向外推,能够增加边界层厚度.     

γ-Reθ转捩模型在二维低速问题中的应用

介绍了基于当地关联的γ-Reθ转捩模型.通过Spalart建议的环境源项法对自由流湍流度衰减的有效控制,实现了该转捩模型的实用化.利用中国空气动力研究与发展中心的TRIP软件平台上实现的该转捩模型,对S809翼型、NLR7301两段翼型和30P30N三段翼型3个低速翼型绕流问题进行了计算,将不同迎角下的计算结果与试验结...     

平流层飞艇艇身外形研究

为了完成特定的任务,平流层飞艇需要克服风场保持长期定点,因此要求其阻力最小。飞艇总阻力中艇身阻力占60-70%,因此对飞艇而言,针对艇身外形进行研究得到阻力小而且实际可用的外形是非常重要。本文采用势流-边界层耦合方法与混合遗传算法对平流层飞艇艇身的外形进行了优化。外部势流采用在艇身表面分布点源的Hess\|Smith面元法求解,边界层计算采用积分边界层方法,阻力系数采用
Squire\|Young方法计算得到。最优外形通过由遗传算法和Nelder\|Mead单纯形法组成的混合遗传算法优化得到。通过优化分析得到了一种实际可用的优化外形,具有在湍流和层流两种流态下阻力系数都比较小的优点。     

超声速边界层转捩拟序结构大涡模拟

针对超声速边界层转捩问题，以五阶迎风和六阶对称紧致格式数值求解三维可压缩滤波Navier-Stokes方程，对马赫数4．5，雷诺数10000的空间发展超声速平板边界层的谐波扰动转捩问题进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法，亚格子尺度模型为修正Smagorinsky涡粘性模型。通过入口叠加一对线性最不稳定第一模态斜波扰动的方法，得到了从线性，弱非线性扰动增长、交替A涡结构出现到演化为发卡涡的转捩过程；针对剪切层结构等现象，给出了该转捩拟序运动的详细讨论。比较显示，转捩结构及摩擦系数曲线等同理论分析吻合。     

壁面温度分布的高超声速平板边界层稳定性分析及转捩预测

针对来流马赫数为4.5、6.0和7.0的高超声速平板边界层,取30km高空处的气体参数,壁面为等温、绝热和温度分布等3种不同条件,采用eN方法进行转捩预测。其中,壁面温度分布条件下,在等温壁（冷壁）和绝热壁之间,给出4种流向分布进行分析。取扰动的初始幅值为0.3‰,以幅值达到1.5%作为转捩判断依据。结果表明:当温度为来流温度时,等温壁面条件的转捩位置最靠前,并随马赫数增大更加靠前;绝热壁面条件的转捩位置随马赫数增大而推后;壁面温度分布条件下,在相同时刻,马赫数越大,转捩位置越靠前。相同马赫数下,壁面温度较高者转捩位置较靠后（马赫数为7.0时,不完全满足此规律）。在马赫数为4.5和6.0时,绝热壁面条件转捩由第一模态主导,其余情况主导转捩的都是第二模态。