使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

高超声速边界层转捩的若干问题及工程应用研究进展综述

高超声速飞行器绕流存在着激波、边界层、流动分离、稀薄气体效应和高温气体效应等多种复杂流动现象的空气动力学问题,其中高超声速边界层转捩既是空气动力学的基础问题,也是高超声速流动研究的热点和难点。若能对边界层转捩进行准确预示及有效控制,则可以实现对飞行器气动力热特性的精细设计,改进飞行器性能,提高任务执行能力。文章针对工程中具有复杂外形飞行器存在的典型失稳特征进行了研究进展回顾,提出了工程实际中亟需解决的复杂边界层转捩问题,明确了高超声速边界层转捩研究的工程应用方向。文章最后还对高超声速边界层的流动控制进行了回顾,以期在今后高超声速飞行器设计中实现对边界层的流动控制,提高飞行器的飞行性能。     

“流动稳定性与转捩”专栏简介

正转捩是指流动由层流状态发展为湍流状态的失稳过程,是经典物理学遗留的少数基础科学问题之一,与湍流问题一起被称为"世纪难题"。高超声速边界层转捩会导致摩阻和热流增大3～5倍,严重影响飞行器的性能与飞行安全。"十三五"期间,国家科技部实施了"大科学装置前沿研究"重点研     

超声速压缩拐角激波/边界层干扰动力学模态分解

压缩拐角激波与边界层干扰问题广泛存在于高速飞行器的外部和内部流动中,其非定常复杂流场结构对飞行器气动性能影响显著。动力学模态分析将有助于进一步加深理解激波与边界层干扰流场不同特征频率对应的流动结构及动力学特性,为揭示其复杂流动机理提供参考。本文采用动态模态分解(DMD)方法对来流马赫数为2.9、24°压缩拐角内激波与超声速边界层干扰下的非定常流动进行了模态分析。评估了稀疏改进动态模态分解方法在压缩拐角流动中的适用性,研究了湍流干扰和转捩干扰下典型特征频率对应的动力学模态空间结构差异及其原因,分析了转捩边界层展向非均匀性对低频/高频模态动力学机制的影响规律。研究发现,湍流干扰与转捩干扰下拐角干扰区内均存在两类截然不同的动力学模态:低频模态和高频模态。低频模态结构集中在分离激波及分离泡剪切层的根部,表征为分离泡的大尺度膨胀和收缩运动;高频模态空间分布则以平均声速线附近正负交替结构为主,对应为边界层内不稳定波沿剪切层往下游的传播。转捩边界层的展向结构对低频模态运动特性影响明显,而对高频模态的影响则相对较小。     

基于湍流模型的转捩预测方法研究

阻力的精确计算,一直是低雷诺数飞行器设计的关键.而与附面层分离和转捩相关的摩擦阻力的精确计算则是此类飞行器阻力精确计算的关键.因此,开展了基于湍流模型的转捩预测方法的研究,应用几种湍流模型对平板边界层进行计算分析,研究各湍流模型预测流动转捩的能力.由于各模型对扰动的过于敏感,引起转捩提前发生,使预测与实验差距较大,由于k-ωSST两方程湍流模型预测能力优于其它模型,在此模型基础上引入了间隙函数进行粘性系数修正,对平板和NACA 0012翼型的转捩重新进行了预测,计算结果表明改进后的模型对转捩位置具有较好预测能力.     

高超声速边界层转捩实验综述

高超声速边界层转捩直接影响飞行器表面的摩擦系数与热流分布,对于高超声速飞行器的气动布局以及热防护设计至关重要。尽管高超声速边界层层/湍流转捩的相关研究已经开展长达半个多世纪,但是由于高超声速流动的复杂性以及触发转捩的因素繁多,研究人员对于转捩过程的认识并不透彻,阻碍了先进高超声速飞行器的设计。地面风洞实验作为高超声速空气动力学设计的重要手段之一,在可预见的将来仍是研究高超声速边界层转捩不可或缺的方法。本文以高超声速边界层稳定性与转捩的风洞实验为重点,按照边界层自然转捩的发展过程,分别回顾了国内外在边界层感受性问题以及线性化阶段风洞实验研究的现状,文章最后总结了风洞实验在未来高超声速边界层转捩研究中的工作与意义,并针对未来的实验研究给出了几点建议。     

高超声速三维边界层转捩数值研究进展及预测软件

高超声速边界层转捩问题已经成为飞行器设计需要考虑的一个关键因素，也是高超声速流动研究的热点和难点。三维边界层中存在横流、流向涡、接触线三种典型的流动结构和失稳特征，他们诱导转捩的机理不同。本文总结了低雷诺数湍流、γ-Re_θt、k-ω-γ、层流动能4种常用转捩模型的构造、适用性及应用情况，介绍了适用于三维边界层的BiGlobal、 PSE3D、 TriGlobal三种全局稳定性分析方法，分析了基于转捩预测e^N方法的e Malik、 LASTRAC、 LILO、HyTEN等转捩软件的优缺点，展望了高超声速三维边界层转捩研究及软件开发的未来。     

雷诺数变化对翼型边界层发展及失速特性的影响

翼型失速及其边界层发展是飞行器设计中的基础科学问题,而雷诺数变化对其影响很大。针对后缘失速翼型,采用Menter k-ω SST模型及耦合扰动放大因子输运方程的转捩模型,进行雷诺数变化对层流-湍流转捩边界层特性和失速特性的影响分析。结果表明:雷诺数增大时,对于转捩边界层,当地涡量雷诺数增大,转捩前移且分离泡减小,流动能量耗散减小,翼型整体表面剪切效应增强,动能更充沛,流动自持能力增强,压力分布可以维持较长距离的梯度抵抗分离能力增强;因此雷诺数增大使翼型失速迎角提高、升力系数增加。     

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

高超声速边界层转捩机理及应用的若干进展回顾

高超声速边界层转捩对飞行器的热传递、表面摩阻和流动分离等有重要影响,尤其是再入飞行器和吸气式巡航飞行器。然而,人们对边界层转捩机理中的很多问题认识还不清楚,或存在争议。本文从扰动波演化的角度回顾了高超声速边界层感受性、线性稳定性和非线性作用的国内若干研究进展,并以基于谐波共振的人工转捩技术为例示范了这些机理认识在转捩控制上的应用。扰动的产生和发展是认识边界层转捩机理的核心。通过研究扰动波来认识边界层转捩机理,开展应用创新研究对提升飞行器性能具有重要意义。     

基于改进的k-ω-γ转捩模式预测高超声速飞行器气动特性

为验证边界层转捩对高超声速飞行器气动特性的影响,采用改进的k-ω-γ转捩模式对类X-51A高超声速飞行器进行了全机边界层转捩预测.不仅系统分析了飞行攻角和雷诺数对边界层转捩的影响规律,同时研究了边界层转捩对飞行器气动力和进气道性能的影响.发现边界层转捩对飞行器升力系数和俯仰力矩系数影响较小,对阻力系数影响较大.针对本文...     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

高超声速进气道强制转捩流动的大涡模拟

吸气式高超声速飞行器常在进气道边界层内布置粗糙颗粒或涡流发生器强制流动转捩为湍流以确保发动机正常启动。为了清晰认识强制转捩过程,采用隐式大涡模拟方法,对强制转捩问题开展了数值模拟研究。针对钻石形和斜坡形涡流发生器,计算得到涡流发生器诱导的流动结构,显示出强制转捩流动由涡流发生器产生的反向旋转流向涡对主导。扰动沿流向增长和发展,导致流向涡对以偶模式或奇模式失稳,偶模式失稳产生对称形式的涡对破碎,而奇模式失稳则导致非对称（弯曲）形式的涡对破碎。流向涡对破碎后产生一系列发卡涡并最终促使边界层转捩为湍流。最后就计算网格和数值耗散对隐式大涡模拟结果的影响以及计算的收敛性进行了讨论。     

高超声速边界层转捩研究现状与发展趋势

高超声速飞行器边界层容易经历层流/湍流转捩,层流流动和湍流流动在摩擦阻力、热交换、噪声和掺混等方面有巨大差别,转捩问题已成为制约高超声速技术突破的基础科学问题之一,是当前国际学术研究的热点与难点。本文详细分析了国内外高超声速边界层转捩研究现状,并将其归为三类:已知主要原因的现象与规律、已知部分原因的现象与规律、未知或矛盾的现象。其中已知主要原因的现象与规律包括壁温、马赫数和噪声影响;已知部分原因的现象与规律主要有头部钝度、熵层和攻角影响;未知或矛盾的现象主要有单位雷诺数影响、转捩区长度、转捩区摩阻和热流分布等。同时介绍了高超声速边界层转捩影响因素研究、转捩机理研究、转捩预测方法及模型研究、促进/推迟转捩的控制方法研究、以及一些公开的飞行试验等方面的进展。最后指出,在今后的高超声速边界层转捩研究中,建议把单个影响因素独立出来研究,尽量避免多因素相互干扰;高超声速边界层失稳研究需要特别关注横流失稳、熵层和模态相互作用;转捩预测需考虑三维边界层和来流扰动的影响;转捩控制研究应重点关注高效、低阻、低热的控制方法;转捩飞行试验十分重要,飞行试验和静音风洞发挥的作用会越来越明显。过去60多年的研究经验表明在未来的研究中应该注重多种手段相结合。     

应用一个新的具有转捩敏感性两方程低雷诺数湍流模型预测转捩流动

本文基于Launder-Sharma(LS)低雷诺数的 两方程湍流模型发展了一种新的具有转捩敏感性的低雷诺数湍流模型。针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,新模型进行了改进。新模型只使用当地物理量,不需要求解壁面距离、Y+和边界层积分参数。新模型能够适用于广泛的流动,且容易应用到通用的CFD程序中。对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,新模型能够预测转捩流动,并能对自由湍流变化给出合理的响应。 关键词：航空、航天推进系统;转捩模型;转捩流动;数值模拟     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

一种基于湍流模式的转捩预测方法

通过对平板边界层进行计算分析,基于湍流模式来预测流动转捩是可行的,但预测出的转捩位置不准,比实际情况靠前。因此,本文从转捩流动的物理特征出发对湍流模式中的混合尺度作出了修正,提出了一种考虑间歇函数影响的k-ωSST两方程湍流模式。采用该模式来对平板湍流边界层和不同攻角下的NACA0012翼型绕流进行了转捩预测,计算结果表明,改进后的模式对转捩位置具有较好的预测能力。     

平板边界层中展向波包型扰动引起的转捩

过去在平板边界层转捩及湍流的研究中,主要考虑的扰动在展向是均匀分布的,这样有利于研究,但在实际问题中扰动形式是多样的,边界层可能是三维的,扰动在展向是不均匀的。对于以往研究的扰动来说,三维平板边界层中的展向非均匀扰动是比较复杂的扰动形式,更接近自然转捩,因此研究这种扰动引起的转捩和湍流具有重要的实际意义。基于此,针对平板边界层,控制方程为无量纲化的Navier-Stokes扰动方程,时间上采用三阶精度的差分格式,空间上展向采用伪谱方法,流向和法向采用高阶精度差分格式,应用数值模拟的方法研究了小幅值和有限幅值展向波包两种情况。通过数值模拟和线性稳定性理论分析小幅值波包的演化,得到了小幅值波包的演化符合线性稳定性理论(LST);分析了有限幅值的展向波包型扰动引起的转捩和湍流,描述了物理空间和谱空间中波包型扰动的演化特征;同时针对不同展向位置进行分析,结果表明不同展向位置的转捩位置不同,但转捩过程和特征是类似的。