密切曲面内锥乘波前体进气道设计和试验研究

介绍了密切曲面内锥乘波前体进气道（Osculating Inward turning Cone Waverider Inlet,OICWI）的一体化设计方法,对该型乘波前体进气道的性能进行了数值分析,针对该型一体化乘波前体进气道完成了风洞试验研究。理论设计结果和设计状态无粘模拟结果一致,设计状态下的计算结果表明,前体进气道具有较高的总压恢复、较好出口流场均匀度及较高的流量捕获率。试验研究结果表明,改型一体化前体进气道在马赫数5～7条件下顺利启动,流场波系及压力分布同数值分析结果吻合。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

基于损伤权重的混合多钉连接件疲劳寿命预测方法

复合材料与金属材料混合多钉连接形式是当代飞机结构中最常见的连接形式,因此对于混合多钉连接件疲劳性能的研究有助于提高对飞机结构疲劳损伤的认知。针对以Ti-6Al-4V钛合金为螺栓的ZT7H/QY9611碳纤维增强树脂基复合材料连接板与30CrMnSiNi2A金属连接板混合多钉连接结构进行数值分析和试验研究。利用有限元方法,分别对复合材料和紧固件进行疲劳损伤预测,估算复合材料连接板上螺栓孔附近的损伤,依据复合材料和紧固件的损伤量所占权重,提出了以紧固件分布和复合材料连接板铺层层数为参数的经验公式,进而有效提高混合多钉连接结构疲劳寿命的预测精度。利用试验结果将螺栓孔损伤形式进行分类讨论,探索混合多钉连接件的损伤演化方式;利用C扫描技术得到复合材料分层损伤结果,与模拟结果进行对比分析,进一步解释了模型的合理性。与试验结果对比可以看到,考虑损伤权重的寿命预测值与试验值的对数误差仅为1.1%,相对于不考虑损伤权重方法的8.4%的对数误差,该模型寿命预测精度显著提高。     

航空发动机帽罩热气膜防冰的加热特性

针对复合材料帽罩采用的一种热气膜防冰系统,采用计算流体力学方法研究了其内部冲击换热和外部气膜加热效果。发展了二维轴对称计算方法并与全三维计算结果进行了对比,搭建了热气吹风实验台,利用红外测温仪和热电偶测量帽罩温度,验证了数值计算方法。全面研究了不同射流雷诺数、相对冲击距和气膜缝结构参数对热气加热特性的影响,结果表明:增大射流雷诺数有利于提高内部冲击换热效果和外部气膜加热效果;最佳相对冲击距随射流雷诺数的增大而增大,射流雷诺数为10000~40000范围内最佳相对冲击距在5~8内取得;气膜射流能大幅提高缝后外壁面温度,加热效率与气膜缝宽和位置有关,前缘开设的气膜缝还会提升前缘冲击换热效果。      

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

基于混合方法的膨胀腔非稳态流动及气动噪声分析

采用混合方法对气流作用下的膨胀腔的非稳态流动和气动噪声进了数值仿真分析。建立基于有限体积法的流体仿真模型,采用大涡模拟方法对膨胀腔内部非稳态流动进行计算。数值仿真结果捕捉到了腔体内部非稳态流动以及涡量的发展变化过程,并与实验观测结果基本吻合。采用积分插值的方法将流场结果插入声学网格上,基于气动声学理论计算腔体内部的流噪声源。建立声学仿真模型,将流噪声源导入声场网格中计算远场响应点的声压,并与实验结果进行对比,声压级以及共振频率都吻合良好。结果表明:混合仿真方法准确直观揭示出管内流噪声产生的机理和过程,尾管辐射噪声的仿真结果与实验结果及理论计算结果的偏差都在5%以内。      

细长圆锥前体非对称涡流场的等离子体控制

应用一对单介质阻挡放电(SDBD)等离子体激励器对顶角为20°的圆锥-圆柱组合体圆锥段分离涡流场进行了主动控制试验研究。试验在3.0m×1.6m低速低湍流度风洞中进行,迎角为45°,基于圆锥段底面直径的雷诺数为5×104。流动控制分为等离子体激励器关闭,左舷或右舷等离子体激励器分别开启,左右舷等离子体激励器占空循环3种模式。试验结果包括7个测量截面上的周向压力分布以及积分得到的截面当地力和力矩以及圆锥段力和力矩。研究结果表明,在圆锥头部尖端处迎风面两侧对称放置一对SDBD等离子体激励器,采用合适的激励器形式,并通过适当的电学参数,可以实现对细长旋成体侧向力和力矩的比例控制。通过对模型及等离子体激励器制作的改进,相对于前人相应的研究结果,本文中侧向力和力矩随占空比变化的线性度得到了改善。     

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

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

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

马赫数为6的高超声速钝锥湍流边界层空间演化的直接数值模拟

首先求解了来流马赫数为6的零攻角高超声速钝锥边界层的层流基本流场,在选定的计算域入口引入一组有限幅值的T-S波扰动,用高精度差分格式对流动进行了直接数值模拟.引入的扰动触发了转捩,从而得到了空间模式下的湍流边界层.研究了湍流平均场与脉动场的统计特性,给出了相干结构的流动显示图,并对强雷诺比拟的结论进行了检验.     

Large-eddy simulation of shock-wave/turbulent boundary layer interaction with and without SparkJet control

The efficiency and mechanism of an active control device ‘‘Spark Jet" and its application in shock-induced separation control are studied using large-eddy simulation in this paper.The base flow is the interaction of an oblique shock-wave generated by 8° wedge and a spatially-developing Ma = 2.3 turbulent boundary layer.The Reynolds number based on the incoming flow property and the boundary layer displacement thickness at the impinging point without shock-wave is20000.The detailed numerical approaches were presented.The inflow turbulence was generated using the digital filter method to avoid artificial temporal or streamwise periodicity.The numerical results including velocity profile,Reynolds stress profile,skin friction,and wall pressure were systematically validated against the available wind tunnel particle image velocimetry(PIV) measurements of the same flow condition.Further study on the control of flow separation due to the strong shock-viscous interaction using an active control actuator ‘‘Spark Jet" was conducted.The single-pulsed characteristic of the device was obtained and compared with the experiment.Both instantaneous and time-averaged flow fields have shown that the jet flow issuing from the actuator cavity enhances the flow mixing inside the boundary layer,making the boundary layer more resistant to flow separation.Skin friction coefficient distribution shows that the separation bubble length is reduced by about 35% with control exerted.     

高超声速飞行器跨流域气动力／热预测技术研究

分析了跨流域流动特点、工程背景需求、国内外研究现状,针对高超声速飞行器跨流域气动力／热预测存在的主要问题,从粘性干扰参数与克努森数搭接的跨流域模拟准则、微量天平结构优化及测力技术、红外大面积中低量值热流测量技术、流场显示与测量技术和 N-S／DSMC 自适应紧耦合数值模拟方法等方面,总结了近年来取得的研究进展,并讨论了下一步研究方向。     

中心进气转静系转盘风阻扭矩数值模拟

采用计算流体力学(CFD)商业软件CFX从流动参数及几何结构两个方面对中心进气转静系盘腔内转盘风阻扭矩特性进行了数值模拟研究,同时搭建一个转盘风阻扭矩实验台测量自由盘风阻扭矩.中心进气转静系转盘风阻扭矩与经典经验公式进行了对比,结果表明:利用数值模拟方法,得到的无量纲扭矩系数与经典经验关系式计算得到的数据总体相差小于10%,利用湍流参数来表征转盘转速以及叠加流流量,对于分析转盘风阻扭矩非常便利.利用湍流参数定量分析时需注意叠加流以及转速的量级同样影响着风阻扭矩的数值.几何结构对于风阻扭矩也有影响,但相比流动参数影响较小.      

基于IPDG方法的超声速混合层流动数值模拟

为了满足超声速混合层高精度模拟需求,实现了基于内罚方法的间断伽辽金(IPDG)方法数值模拟.通过将黏性通量作为辅助变量使得Navier-Stokes方程降阶,并利用间断伽辽金方法进行空间离散,最后采用Newton-Krylov隐式方法对空间半离散方程进行时间推进.相对于有限体积法数值精度提高到了 3阶.将该方法应用于对...     

无网格法耦合RNGk-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

 基于移动最小二乘无网格方法,耦合RNG (Re-Normalisation Group) k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM (Advection Upstream Splitting Method)⁺-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP (Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNG k-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

压缩性对细长体涡流非对称发展的影响

通过数值方法对大迎角细长体湍流流场的模拟,探讨压缩性对细长体非对称绕流发展的影响。结果表明细长体顶端的极小扰动诱发显著的非对称绕流,非对称的涡系结构沿轴向是逐步演化的;在亚临界横流马赫数区间,马赫数越高非对称越显著;在超临界横流马赫数区间,细长体两侧出现横流激波,非对称的发展被抑制,马赫数越高非对称越弱。     

基于非Favre过滤的动力学非线性亚格子模型

采用直接过滤的Navier-Stocks（NS）方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,并仅保留级数的一阶导数项。因为截断级数产生的误差相当于丢失了模型的部分耗散作用,采用一种动力学非线性亚格子模型来补偿丢失的耗散影响。用模型模拟了超音速湍流平板边界层及绕双椭球高超声速湍流,并将所得结果与经验公式和计算结果进行了对比。对比表明模型可以有效的模拟可压缩湍流流场。     

基于混合算法的涡轴发动机稳态性能仿真模型

以自由涡轮式涡轴发动机为研究对象,建立了涡轴发动机的稳态性能仿真模型,提出了基于天牛须算法和N+1点残量法的求解发动机模型的混合算法(BAS-N+1混合算法),利用发动机台架试车数据对仿真计算结果进行了验证。结果表明,该稳态性能仿真模型各参数的求解误差在3%以内。与PSO-N+1混合算法相比,BASN+1混合算法求解精度更高,收敛更快。BAS-N+1混合算法既保留了智能算法对初猜值误差的包容性,也拥有接近经典迭代算法的收敛速度和精度,能够实现涡轴发动机稳态仿真模型的高精度大范围快速收敛。     

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

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

Hybrid method based on flame volume concept for lean blowout limits prediction of aero engine combustors

The Lean Blowout (LBO) limit is crucial for the aircraft engines. The semi-empirical (such as Lefebvre’s LBO model and Flame Volume (FV) model), numerical and hybrid methods are widely utilized for the LBO limit quick prediction. An innovative hybrid method based on the FV concept is proposed. This method can be classified as a semi-empirical/physical based hybrid prediction method. In this hybrid method, it is assumed that the flame volume varies nearly linearly with the fuel/air ratio near the LBO. The flame volume is obtained directly by the numerical simulation using the threshold value of the visible flame boundary as 900 K. Then the final LBO limits is determined by the FV model. On the premise of keeping the good generality of prediction, the hybrid method based on the FV concept can further improve the prediction accuracy. The comparison with the prediction of the existing available methods on fifteen combustors shows that the hybrid method based on the FV concept achieves better prediction accuracy. The prediction uncertainties between the experimental results and the predicted values by the hybrid method based on the FV concept are within about ±10%.