宽速域翼尖涡及其与斜激波相互作用

针对宽速域飞行器产生的翼尖涡及其与斜激波相互作用问题,采用数值模拟方法,在来流马赫数0.2～6.0范围内,研究了翼尖涡的演化特性,并探讨了斜激波诱导翼尖涡破碎的机理。结果表明,在宽速域范围,翼尖涡的切向速度和环量分布仍然符合低速流动中的自相似关系式;在超声速至高超声速流动中,翼尖涡强度沿流向更快地减弱;当翼尖涡与斜激波相互作用时,涡核处压力越低的翼尖涡更容易破碎,而经典的涡破碎理论难以在宽速域范围内准确预测涡破碎现象;通过引入涡核处的压力亏损效应修正经典理论,大幅提升了预测斜激波诱导翼尖涡破碎的能力。     

流向涡与斜激波相互作用的试验研究和数值模拟

为研究流向涡与斜激波相互作用在超声速燃烧中的应用,进行了由翼产生的流向涡与楔块产生的二维斜激波相互作用的燃烧室冷流试验研究.在不同激波强度下,纹影仪捕捉到了强、中、弱不同的涡/波作用现象.仿真与试验结果符合得较好.试验结果表明:对马赫数为2.3流场,在翼攻角10°时,能产生强流向涡,此工况下,锲角越大,涡/波作用越强.仿真结果表明:马赫数对涡/波作用影响较大,总压影响不明显,总温可影响亚声速回流区的尺寸.      

斜激波与内凹半圆柱面无粘相互作用

针对三维内转式进气道流动中激波入射内凹壁面主导的复杂波系干扰问题，采用斜激波入射内凹半圆柱面的简化构型，在来流马赫数为6的条件下，通过无粘数值模拟结合理论分析，研究了激波角βi=14°～29°的斜激波在内凹半圆柱面反射形成的三维流场。结果表明，流场对称面均会出现显著高于二维情况的逆压梯度。当βi≤25°时，从侧壁到对称面，斜激波经历了从马赫反射（MR）到规则反射（RR）的转变，形成了MR-RR型流场，转变点处产生的桥激波向对称面延伸，桥激波在对称面反射后产生的压力峰值高于二维斜激波反射；当βi≥25°时，斜激波在侧壁和对称面均发生马赫反射，形成了MR-MR型流场，两种马赫反射分界点处产生的桥激波向侧壁发展，侧壁气流在对称面相撞后产生的压力峰值高于正激波后的压力；当βi=25°时，流场存在MR-RR型和MR-MR型双解现象。通过降维分析理论，揭示了两类流场中转变点和分界点的形成机制，并厘清了桥激波的产生原因和初期演化特征。当βi≥18°时，无粘激波干扰所主导的侧壁气流加剧向对称面汇聚，并在对称面附近产生流向涡对。无粘分析获得的认识，有助于揭示内转式进气道中流动汇聚和流向涡对等现象的形成机...     

圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。     

跨声速风扇转子间隙流动结构分析

为探索跨声速转子间隙流动结构,归纳间隙泄漏涡（TLV）和激波相互作用机理,以跨声速转子为研究对象,数值模拟不同间隙下流动特性.此外,着重探究TLV与激波的相互作用,斜激波受到间隙泄漏流的干扰,被削弱打断向上游凹曲;建立三维模型,加深对跨声速叶尖区域流动的全面认识.研究表明:随间隙增大,相同流量下,效率越小,压比也越小;TLV强度更大,偏离吸力面程度增大,沿周向和展向影响范围都越大.压力差提供泄漏流迁移动力,间隙提供泄漏流形成通道.选取h/c=1.0%间隙,随着间隙高度的增大,泄漏流周向运动趋势更明显且二次泄漏现象更剧烈;沿泄漏流方向无量纲流向涡量有少量减少,无量纲螺旋度较高,集中涡特性明显.      

圆形出口内转式进气道流动特征

采用数值仿真的方法研究了内转式进气道的流动特征。研究表明:设计状态在近壁面唇罩激波诱发了二次流,进而发展形成流向涡,造成低能流堆积,隔离段出口流场分布不均,消弱了进气道的抗反压能力。有攻角条件下,口面激波偏离唇罩前缘,激波形态发生改变,激波波面中部展向具有准二维特性,压缩面两侧气流压缩变弱,激波层变薄,出现局部膨胀区;有攻角条件下的无黏流场,在进气道压缩段形成三维流向涡,该流向涡促进高能高速气流向壁面迁移,改善了黏性条件下隔离段出口流场的均匀度。     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

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

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

高负荷涡轮端区非定常流动相互作用研究

采用三维黏性非定常数值模拟方法研究了某型高压涡轮端区非定常流动相互作用,着重研究了上游静子尾迹与转子二次流的非定常作用机制,同时还分析了负荷分布、激波等对端区非定常流动的影响。结果表明,静子尾迹的非定常作用一定程度减小了转子轮毂二次流的径向涡量;尾迹对流向涡量的影响取决于尾迹沿叶高的分布,当吸力面一侧的尾迹具有与二次流方向相反的流向涡量时,二次流的流向涡量减小;非定常效应还使得转子叶片根部负荷略为减小,也一定程度上抑制了转子轮毂二次流的发展。此外,受静子尾缘激波的影响,转子叶片表面负荷分布发生明显的周期性变化,导致叶片表面承受较强的非定常力,在涡轮设计中必须考虑。另外,通过计算涡轮级中的熵增和熵产,定量地分析了端区非定常相互作用产生的损失,并得到了一些有意义的结论。     

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

跨音风扇激波/漏流涡干扰模型探讨

激波/漏流涡干扰系旋转运动与突然的逆压力梯度的强相互作用,引起各自形态和强度发生变化。涡破裂对风扇/压气机效率和稳定性十分不利。本文探讨了激波/漏流涡干扰的物理模型,通过求解涡核运动方程,分析了涡核特征参数对涡破裂的影响,为控制涡破裂提供了根据。      

绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.     

Shock wave configurations and reflection hysteresis outside a planar Laval nozzle

When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection(RR) pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection(MR) pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited condition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these conditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RRMMR transition in the dual-solution domain.     

跨声速压气机转子叶尖非定常流场数值研究

采用时间精确求解方法对某高负荷跨声速轴流压气机转子在98%设计转速下的叶尖非定常流场进行了数值研究.结果显示,激波自身振荡不明显,叶尖区域流动的非定常性主要来源于叶尖泄漏涡的破碎及其与激波之间的相互干涉.对比设计状态与失速状态下叶尖泄漏涡的特点发现:在近失速点时,叶尖区域间歇性出现前缘溢流.分析表明:叶尖泄漏涡在激波后破碎是造成堵塞的主要原因,也是造成spike型失速初始扰动的原因.      

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

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

激波/边界层干扰及微型涡流发生器控制研究进展

激波/边界层干扰是一种发生在超声速/高超声速流动中的普遍现象。该现象将引起分离、流场结构振荡、局部高热通量和压力载荷。主要总结了近十年来激波/边界层干扰特性与微型涡流发生器及其组合体在流动控制中的最新进展。微型涡流发生器是目前研究最多、应用最广泛的控制方法,其流动机理和控制特性被大量挖掘。为了适应来流条件的变化、满足实际工况的需要,应开发定量评估和参数化设计方法。同时,应探索微型涡流发生器与其他控制方法的组合,实现更大程度、更广范围流场的控制。     

可压缩球涡和平面激波相互作用的数值研究

行进了可压缩球涡和平面激波相互作用的数值研究。一方面，用非定常Ｅｕｌｅｒ方程数值模拟了一种可压缩球涡的自身演化，同时，利用Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系，在流场中嵌入运动激波，求解了激波－球涡相干的流动过程。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...