采用γ-Re_(θt)模型的转捩流动计算分析

为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Reθt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序(HUNS3D)。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Reθt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明:γ-Reθt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

高超声速计算中的气体动理学格式

回顾了高超声速连续流部分的计算流体力学(CFD)方法,总结了近些年兴起的气体动理学格式。阐述了该格式的构造机制,强调了将物理规律直接用于构造数值方法的思路。结合一些应用实例,例如激波相互作用、激波边界层相互作用以及边界层分离等高超声速问题,说明了这种构造思路给数值模拟带来的优点。从高超声速的发展历程来看,气体动理学格式的构造过程包含了更基础的物理规律,而且具有多尺度的特性。这些特性有助于研究复杂的高超声速问题。介观或者微观角度直接构造数值方法的发展趋势为高超声速计算工具指出了可能的发展方向。     

后掠翼边界层横流不稳定转捩预测模型

通过对经典Falkner-Skan-Cooke三维边界层相似解的理论分析和数值求解,结合二维边界层转捩判据的思想,采用由试验数据标定的C1准则关系式求解横流不稳定转捩位移厚度雷诺数,建立了针对固定前缘后掠角机翼的横流转捩判据,并且通过方程求解和数据拟合得到了该转捩判据的数学结果.应用该模型对30°前缘后掠角的ONERA-M6机翼和45°前缘后掠角的NLF(2)-0415无限展长机翼进行了横流不稳定转捩数值模拟.模拟结果显示:改进后的转捩模型预测所得到的转捩位置精度较高,均与后掠翼横流试验数据吻合较好,从而证明了构建的横流不稳定转捩判据的合理性和实用性.      

电弧放电等离子体激励控制超声速压气机叶栅激波/边界层干扰仿真研究

为研究电弧放电等离子体激励对超声速压气机叶栅激波/边界层干扰的控制作用,建立了模拟等离子体激励作用效果的唯象学模型,进一步以ARL-SL19超声速叶栅为对象,通过数值仿真研究了电弧放电等离子体与叶栅通道内部流动的相互作用及其对叶栅流动损失的影响。结果表明:等离子体唯象学模型能够较好模拟电弧放电等离子体诱导产生冲击波的气动特性。电弧放电等离子体激励对叶栅通道内部流动主要具有三种作用效果:在放电区,注入的热量会产生阻塞效应,增加近壁面气流的流动损失;在激波/边界层相互作用区,能够改变激波系结构,减小激波损失;在尾迹区,冲击波会诱导产生脱落涡。     

凸拐角附近激波与湍流边界层干扰的数值模拟研究

采用大涡模拟方法数值研究了来流马赫数Ma∞=3、雷诺数Reθ=2070的凸拐角附近激波与湍流边界层干扰问题。对于尖拐角和钝拐角两种情况,计算考察了激波入射到拐角下游时的流场结构、非定常分离流以及湍流统计特性。研究发现:尖拐角能够使流动产生局部集中的顺压梯度,导致壁面摩擦系数在拐角处出现峰值,因此分离点难以越过拐角向上游运动,而再附点发生前后运动。与尖拐角情况不同的是,钝拐角产生的局部顺压梯度不足以阻止分离点越过拐角,因此分离点和再附点都会发生运动。此外,尽管凸拐角产生的顺压梯度能抑制分离点附近的近壁小尺度涡结构,但在经过入射激波作用之后,近壁小尺度涡结构在再附点附近显著增强。     

MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

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

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

雷诺数对增升装置流动特性影响的计算研究Ⅱ——缝道流动特性

在1×106~30×106的雷诺数范围内,马赫数为0.197的情况下,计算并分析了NHLP-2D多段翼型的缝道流动规律,提出了依据缝道出口速度分布定义名义边界层厚度δ/c的方法。研究发现,δ/c随雷诺数的增大而单调减小,且减小速度随雷诺数的增大而明显减缓,符合雷诺数对边界层的影响规律,说明本文定义的δ/c可用于多段翼型边界层厚度的定量研究。对于缝翼缝道,主翼壁面处的δ/c高于缝翼尾缘处的,且二者随迎角的变化规律相反。襟翼缝道处于主翼的强下洗流场中,襟翼上的δ/c随迎角几乎不变。当1×106≤Re≤2×106时,缝道边界层总厚度δT/c随雷诺数有明显的变化,当Re≥3×106时,δT/c随雷诺数的变化率减小。本文研究范围内,δT/c都没有严格地进入雷诺数自准区。当Re≥15×106时,δT/c随雷诺数接近线性变化趋势,为雷诺数规律的外推提供了参考。     

不同波系配置的鼓包压缩面流动特性实验研究

为探究波系配置对鼓包压缩面流动特性的影响,通过纹影和基于纳米粒子示踪的平面激光散射系统(NPLS)观测以及压力测量等实验手段,研究了三种不同波系配置的鼓包压缩面诱导的复杂流动,获得了鼓包压缩面上的压力分布、边界层分布特性以及鼓包压缩面上的精细流场结构。结果表明:三种不同波系配置的鼓包压缩面均具有边界层排移能力,在鼓包对称面上,沿流向边界层不断降低,边界层不断向两侧排移,在鼓包侧面形成堆积。尽管如此,压缩过程的不同,边界层的排移过程存在一定差异,对于单锥鼓包,其最先将低能边界层排移掉,而双锥和等熵锥鼓包略为迟缓,但三种鼓包末端的边界层厚度大体一致。此外,三种鼓包压缩面对称面上前中部压力梯度存在较大差异,但鼓包末端具有相同的压力分布。因此,在保证总偏转角一致的情况下,具体选用哪种鼓包配波方式对于边界层排移、增压能力而言没有显著的差异。     

局部凹槽对高超声速边界层中第二模态扰动演化的影响

采用直接数值模拟方法研究了局部矩形凹槽对来流马赫数为6.0的高超声速平板边界层中第二模态扰动演化的影响。引入了透射系数(定义为凹槽下游与上游扰动渐近幅值的比率)来量化凹槽的局部散射效应。数值结果表明:对于较浅的凹槽,频率较低的第二模态扰动被促进,而高频扰动的规律相反。对于大多数情况,凹槽深度的增加使得透射系数减小,这表明凹槽对扰动的促进作用减弱或抑制作用增强。当凹槽的深度超过某一临界值时,透射系数与凹槽深度的关系变为正相关,这表明了另一机制的出现。临界深度受扰动的频率影响:频率越低,临界深度越大。