跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

基于大涡模拟的压缩拐角激波/边界层干扰研究

为了提高对激波/边界层相互干扰的基本理解,采用大涡模拟(LES)对来流马赫数Ma=2.9,转角为24°的压缩拐角激波与湍流边界层相互干扰进行了研究。采用回收/调节方法作为入口湍流生成技术,并在超声速平板湍流边界层上进行了验证。采用涡识别方法和数值纹影图等流动显示方法,研究了干扰区内激波与边界层相互干扰的结构变化特征。通过对比分析湍动能和雷诺正应力在不同流向位置的分布规律,研究表明:经过激波干扰后湍动能主要集中在边界层的外层,并在拐角附近形成一个低湍动能区;雷诺正应力流向分量和法向分量在边界层内的分布呈现为单峰结构,而展向分量呈现为双峰结构。运用间歇因子对分离激波的大尺度流向运动进行研究,发现激波围绕着平均分离点作前后运动,运动的尺度等于进口湍流边界层厚度的72%。证实了拐角下游G?rtler流向涡对的存在,并对其展向分布和空间演化特性进行了详细研究。     

压缩拐角激波与旁路转捩边界层干扰数值研究

为了研究激波与旁路转捩边界层的干扰机理,采用直接数值模拟(DNS)方法对来流马赫数Ma∞=2.9,24°压缩拐角内激波与转捩边界层的相互作用进行了系统的研究。考察了旁路转捩干扰下压缩拐角内分离区形态和激波波系结构的典型特征。比较了转捩干扰与湍流干扰流动结构的差异,并分析了造成差异的原因。研究了拐角内转捩边界层的演化特性,探讨了转捩干扰下脉动峰值压力和峰值摩阻的分布规律及形成机制。研究结果表明:相较于湍流干扰,两侧发卡涡串的展向挤压使得分离区起始点以V字型分布,且分离激波沿展向以破碎状态为主,激波脚呈现多层结构;拐角内的干扰作用急剧加速了边界层的转捩过程;转捩干扰下的拐角内峰值脉动压力以单峰结构出现在分离区的下游,同时干扰区内的强湍动能和高雷诺剪切应力使得其局部峰值摩阻系数要高于湍流干扰。     

高负荷跨声速涡轮激波损失机理及控制技术研究

为了指导高性能、高负荷跨声速涡轮的设计,对其叶栅内的流场结构、尾缘波系结构、减小激波损失的机理及其控制技术进行了分析研究。结果表明:跨声速涡轮尾缘流场结构复杂,存在分离膨胀波、分离激波、基底区、再附激波、尾迹、吸力面反射波甚至激波边界层相互干扰等流动现象。通过采用收缩-扩张通道,喉道后采用直线型吸力面,减小吸力面尾缘弯折角、尾缘厚度和尾缘附近局部修型等措施,从而减弱激波强度,减小激波损失。     

真实气体效应对升力体舵面局部流动分离的影响

数值模拟分析了高马赫数低雷诺数条件下激波边界层干扰、激波与激波相互作用、流动分离再附等流动现象的特点以及高温真实气体效应的影响。分别采用量热完全气体、平衡气体、化学非平衡气体模型对升力体由于舵面偏转引起的局部流动分离情形进行了数值模拟。研究了飞行高度、壁面温度及来流马赫数对流动分离的影响。计算结果表明:真实气体效应使空气在边界层内发生离解反应,边界层内温度降低,粘性减小,动能损失减小,克服逆压梯度的能力更强,从而使分离区明显减小。分离区的减小改变了分离/再附激波的位置和强度,进而对局部压力及热流分布产生重要影响;随高度增加,平衡气体较完全气体分离区相对减小量增大,平衡气体效应对流动分离/再附现象的影响越大;壁温对分离区影响较大,随壁温升高,分离区增大;随马赫数增大,分离区减小,真实气体和完全气体的差异增大,真实气体效应的影响更加显著。     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

等离子体气动激励控制超声速边界层分离的实验研究

等离子体气动激励与超声速气流相互作用已成为高速流动控制领域的研究热点。激波与边界层相互作用现象广泛存在于超声速飞行器之中。本文进行了等离子体气动激励控制压缩角区和激波诱导边界层分离的实验,通过流场纹影显示和壁面静压测量,研究等离子体气动激励如何影响激波、激波如何影响边界层特性的科学问题。实验结果表明:施加毫秒量级表面电弧放电能够前移压缩角区的诱导斜激波,使分离区后移,分离区域增加,但激波强度减弱,流场总压增加;施加微秒量级表面电弧放电能够抑制激波诱导边界层分离,使分离区减小,流场总压减小。基于实验结果,认为毫秒量级表面电弧放电激励控制超声速气流的主要机理为放电过程的焦耳热效应;微秒量级表面电弧放电激励控制超声速气流的主要机理为焦耳热效应和冲击波效应共同作用。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰，平板边界层上形成分离区，又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游，引入小扰动的TS波后，扰动波传播通过带有这些复杂流动现象的流场时，扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析，研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现，激波的出现不同于一般的压缩波，在亚声速区与超声速区对扰动波演化的影响是不同的，此外，分离对扰动有稳定作用。     

马赫6柱-裙构型激波/湍流边界层干扰摩阻统计特性

为探究激波/边界层干扰表面摩阻统计特性,对马赫6柱-裙构型的激波/湍流边界层流场进行了直接数值模拟,推导了柱坐标形式的时均摩阻分解公式,将其与对流效应、流向不均匀性效应、分子黏性耗散效应、曲率效应和湍动能耗散关联起来,在此基础上对摩阻脉动和时均摩阻的统计特性进行了研究。概率密度分布结果表明干扰区内摩阻脉动概率密度明显偏离正态分布,摩阻脉动在该区域的间歇性较强,这与湍流脉动在激波干扰后存在显著内在可压缩性效应有关;谱功率密度结果表明摩阻脉动的能量在干扰前中频部分的能量占据主导,峰值对应频率约为0.14,干扰区后附近区域峰值位置转移到高频区域。采用摩阻分解公式对激波干扰前后的时均摩阻进行分解,结果表明时均摩阻的主要贡献来自分子黏性耗散效应和湍动能耗散效应;在激波干扰前分子黏性应力效应占据主导位置,干扰后湍动能耗散效应占据主导位置;激波干扰后对流效应明显增强,曲率效应的贡献在激波干扰前后基本保持不变;流向不均匀项贡献由正贡献变为负贡献,且所占比例有所上升,这主要是由压力梯度项的变化导致的。对分解后各项积分内函数沿法向分布的分析表明,分解后各项在边界层内分布呈现较大的差异,特别地,湍动能耗散项...