基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

高马赫数层流摩阻数值计算精度

飞行器在高空高速飞行时黏性效应显著,摩阻预示精度对飞行器的关键气动性能意义重大。目前,摩阻预示主要依赖数值计算,但高马赫数层流摩擦阻力计算与试验测量结果仍存在差距。以具有高速飞行器典型部件特征的球锥、三角翼为对象,结合风洞试验摩阻测量结果,使用国家数值风洞（NNW）数值计算软件和自研CFD程序研究了数值计算中影响摩阻计算精度的格式数值耗散及壁面温度边界条件等重要因素。研究结果表明：数值耗散越小,表面摩阻的计算精度越高;在速度较低的边界层近壁区内关闭熵修正,将有助于提高表面摩阻的预示精度。此外,在高马赫数流动问题的数值模拟中,壁面温度条件对表面摩阻计算同样具有重要影响。最后,基于分析结果和工程需要提出了对高精度摩阻数值预示的研究需求。     

超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

热力学非平衡状态下DSMC绝热壁边界的抽样方法

从直接模拟蒙特卡罗仿真DSMC方法的实现途径出发,经过分析推导,给出了一种在考虑分子转动能、振动能等热力学非平衡状态下统计求解绝热壁温的DSMC抽样方法。然后采用锥体外形的算例验证,分析了该方法在实现绝热壁条件时的统计误差。计算结果表明,此种方法能够有效地实现DSMC方法壁面处的绝热壁条件,随着采样步数的增加,其边界条件实现时的统计误差不断降低,该方法收敛后的统计误差可降至0.5%,能够满足工程要求。     

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。     

立楔诱导高超声速分离流动的被动控制研究

双模态超燃冲压发动机由于压力扰动可能发生不起动现象,造成推力严重下降,对飞行稳定性与飞行安全具有很强的破坏性.不起动初始阶段主要受到激波与边界层相互作用引起的流动分离影响,希望通过控制分离达到改善流动的目的.采用5阶特征型WENO(weighted essentially non-oscillator)格式与3阶TVD(total variation diminishing)型Runge-Kutta(R-K)格式的高精度数值方法,求解三维Navier-Stokes(N-S)方程,研究与分析了凸起物和被动吹吸两种被动控制方法对激波与边界层相互作用导致的高超声速流动分离现象的控制效果.结果表明:凸起物通过诱导流向涡形成,改变空间压力分布,减弱二次分离,影响分离结构;吹吸方式的被动控制技术通过平衡分离区与再附区之间的高低压差,降低逆压梯度,使压力分布与分离区域发生改变.      

翼型湍流尾缘噪声半经验预测公式改进

改进了传统的翼型湍流边界层尾缘噪声BPM半经验预测公式。传统的BPM半经验湍流边界层尾缘噪声预测公式对高攻角和厚翼型在高频范围的预测结果大于实验结果,通过分析比较传统BPM半经验预测公式和Howe翼型尾缘噪声理论模型发现:这主要是由于传统BPM半经验预测公式对压力面声源噪声辐射高估引起的。因此将压力面声源噪声辐射与吸力面声源噪声辐射的幅值比由原来的边界层位移厚度一次方比值改进为二次方比值,进而得到了改进后的BPM半经验预测公式;使用改进后的BPM半经验预测公式对NACA0012翼型在不同来流不同攻角下的噪声辐射进行了预测比较,发现对于NACA0012翼型,改进后的BPM半经验预测公式具有较高精度;另外也预测了较厚的风力机翼型DU-96-W-180,预测结果明显改善。     

上游尾迹入射角对压气机叶片负荷的影响

压气机叶片负荷的提高使得叶片表面边界层更容易分离,利用上游叶排产生的非定常尾迹能够抑制边界层分离.运动圆柱代替上游转子,在保证下游叶片进气速度大小及攻角不变的情况下,改变圆柱运动速度以获得上游尾迹与下游静叶吸力面不同的夹角,发现在低负荷小分离情况下静叶损失系数与上游尾迹入射角无关;高负荷大分离情况下静叶损失系数随上游尾迹入射角的增加而降低.分析尾迹作用下高负荷静叶通道内流场,当进入静叶通道的尾迹与叶片吸力面近似平行时,尾迹诱导边界层增厚.使得叶片表面分离泡随时间大幅值脉动,损失增加;当尾迹与叶片吸力面的角度逐渐垂直时,尾迹抑制了边界层分离,同时叶片表面分离泡位置近似不变,其原因是尾迹以负射流形式进入边界层内部,补充了边界层内部低能流体,使得损失减小.      

CFD中数值耗散和延伸边界对声传播的影响

对已知色散关系的正弦波在二维管道中的传播使用有限体积法进行计算,结果表明二阶精度下每波长40个左右的网格点能够有效地保留声波.在计算域外增加稀疏网格延伸段能够对边界反射起到抑制作用,为以较小的计算量增加实现较好的反射抑制效果,研究给出了延伸段长度和网格轴向尺寸的推荐公式.      

进口导流叶片对转子叶片颤振特性的影响

上游静子叶片的尾迹扰动是下游转子叶片发生强迫响应问题的主要原因,但通常认为不会对颤振特性产生影响。运用CFD技术求解非定常流场,用能量法对典型涡轮风扇NASA Rotor 67进行了气动弹性分析。结果表明:不加进口导流叶片时其叶片最有可能在第二阶模态、60°叶间振动相角的情况下发生颤振;加上IGV(进口导流叶片)后,IGV的尾迹会显著改变转子叶片的颤振特性,当转子叶片的最危险颤振模态频率接近IGV尾迹的扰动频率时,即使单独转子叶片十分稳定,但在尾迹的激励下颤振仍然可能发生;当该频率远离IGV尾迹扰动频率时,尾迹扰动同样会使转子叶片气动阻尼降低。该结果意味着尽管IGV的引入不会明显改变压气机性能,但有可能会对其颤振特性带来明显的负面影响,需要在涡轮机设计中予以考虑。