具有无源控制空腔时正激波/湍流附面层干扰的数值模拟

采用雷诺平均N-S方程和B/L代数湍流模型计算了具有无源控制空腔时正激波/湍流附面层干扰流场。计算与实验结果的比较表明,本文方法可较准确地预测激波结构、激波与附面层干扰区流动基本特征及波后流动分离状态、激波位置、波前马赫数等参数。      

30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

超临界翼型跨声速抖振CFD计算和POD分析

使用二阶迎风Roe格式、隐式时间推进和S-A(Spalart-Allmaras)一方程湍流模型,通过求解基于格心格式有限体积法的RANS方程模拟了OAT15A超临界翼型的跨声速抖振流场。在模拟出激波/附面层相互干扰诱发的抖振现象的基础上,对翼面激波运动过程中的气流分离泡变化规律开展研究,揭示出激波变化相位略微超前于升力系数相位的非定常现象。为剖析跨声速抖振的物理成因,将抖振计算的非定常解作为快照,应用本征正交分解POD(Proper Orthogonal Composition)方法提取POD模态,从流场相干结构的角度对跨声速抖振现象进行了分析。     

基于代理模型的自适应后缘翼型气动优化设计

针对具有自适应后缘的跨声速翼型,基于代理模型和遗传算法相结合的优化方法,开展考虑自适应后缘结构约束的翼型气动优化设计研究。结果表明,低升力系数下翼面流场没有明显的能量损失,不同设计升力系数得到的自适应后缘翼型阻力相差不大;高升力系数下翼面附近存在激波,翼型阻力主要由激波强度决定。对于以低升力系数为设计点的基本翼型,通过后缘自适应变弯来调整载荷分布,可以降低翼面激波强度,减小翼型阻力;对于以高升力系数为设计点的基本翼型,可以直接通过气动优化来消除翼面激波,使得翼型阻力达到最小。因此对于带有自适应后缘的翼型,为了实现宽升力系数范围内的阻力最小,应首先以高升力系数为设计点完成基本翼型的气动优化设计,然后以低升力系数为设计点完成自适应后缘外形的气动优化设计。     

NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

多段翼型高精度数值模拟技术研究

采用五阶精度的加权紧致格式(WCNS)数值模拟了NLR7301两段翼型、30P-30N三段翼型的复杂流场,主要目的是考核WCNS格式模拟多段翼型复杂流场的能力,研究湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应实验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、转捩位置对该翼型压力分布和典型站位速度型的影响。研究结果表明,基于WCNS格式,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对边界层速度型和阻力系数的模拟精度较差;模拟实验的转捩位置可以改善边界层和尾迹区的模拟精度,提高阻力的数值模拟精度。     

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

马赫数2.048的二元平板紊流附面层及其与激波干扰的实验研究

本文介绍了M=2.048的二元平板紊流附面层和附面层与激波(激波板相对气流为8.9°)干扰实验的设备、测试仪器、数据处理以及主要的实验结果。实验用的风洞试验段截面积为0.6×0.6米~2,每米长雷诺数月Re=2.1×10~7,以激波干扰区附面层厚度为参考长度的雷诺数Reδ_0)≈10~6。 二元平板紊流附面层速度型的实验结果接近于u/ue=(h/δ)~(1/7),而且与二元可压缩附面层方程有限差分解的数值结果符合较好。实验结果与二元平板不可压附面层壁面定律和尾迹定律公式计算值的变化规律接近,但计算结果偏大些,而与速度亏损定律的结论则是一致的。 文中亦给出了由实验结果推算的激波附面层干扰时动涡粘系数ε_m沿流线的定性分布规律:在附面层内层区域,ε_m通过干扰区是逐渐增大的,至干扰区下游时恢复到某一大干上游区的ε_m值。而在外层区域,ε_m在干扰区是先减小后再增大,这时下游区的ε_m略大于上游区的ε_m值。     

粘性流绕俯仰振动翼型流动的数值模拟

本文直接从N-S方程出发,利用LU-ADI格式和Baldwin-Lomax代数湍流模型研究了粘性流绕俯仰振动翼型的流动。数值实验表明,本文的数值结果同实验吻合较好。通过数值模拟手段研究了振动对流场中激波和分离这两大主要特征的影响,结果发现;(1)激波滞后于攻角的变化,如当NACA0012翼型在负攻角状态下上翼面存在激波而下翼面无激波。(2)由于翼型的振动,分离被减轻升力相对提高。(3)随翼型攻角和振幅的变化,翼型振动的升力回线走向可不同。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

基于压力分布的开孔翼型气动特性数值模拟研究

以NACA0012为基准翼型,生成C型结构网格,采用SA湍流模型求解NS方程,获得基准翼型压力分布数值模拟结果。找出翼型上表面前缘低压区域和后缘高压区域,在两个区域进行开孔,让后缘高压区域气流通过导流管引流到前缘区域,加速翼型上表面气流流动,增大翼型速度环量。通过基准翼型与开孔翼型气动特性对比分析可知:小迎角下,基于压力分布的开孔翼型阻力系数略有增大,但能有效提高翼型的升力系数。     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

突跃型与平滑型斜爆震波起爆机制数值模拟

采用带有化学反应的Euler方程,对突跃型与平滑型2种形态的斜爆震波(ODW)在起爆机制和结构特点等方面的差异进行了研究.结果表明:斜激波(OSW)后是否存在亚声速区是判定斜爆震形态的依据.当波后火焰抬升斜激波面使得亚声速区存在时,横向激波与三波点结构就会出现并形成突跃型斜爆震.当波后所有区域皆为超声速区时,波后火焰燃烧无法影响到上游激波面,则会形成平滑型斜爆震.通过斜激波关系式给出了横向激波形成所需的临界斜激波角度,只有当斜激波角度大于此临界角度时才能形成突跃型斜爆震波,反之则形成平滑型斜爆震波.      

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

"X"布局高超声速倒置进气道激波与附面层干扰抑制研究

采用标准k-ωSST湍流模型数值计算方法,针对二元高超声速倒置进气道激波与附面层严重的干扰现象,采用分流楔抑制激波与附面层干扰方法,对有无分流楔的进气道性能及流动机理特征进行了详细的研究。结果表明:采用分流楔的流动控制方法,有效抑制了激波/附面层产生的分离包对进气道内流动的干扰;提高倒置进气道的气动性能,进气道的总压恢复系数和流量捕获系数均有提高,计算模型的壁面总阻力系数得到一定程度的减小。数值计算结果表明,在分流楔尾迹中强剪切流动在一定程度上缓解了激波与附面层干扰的强度;在分流楔后缘存在稳定的横向涡,由于气流进入尾迹驻涡是来自附面层外的总压较高的高能流体,提高了附面层的抗逆压能力;由于尾迹驻涡的存在使得分离涡没有向弹体周向扩散,减小了阻力。该方法实现了对高超声速倒置进气道激波/附面层干扰的抑制,揭示了其抑制的机理。     

"X"布局高超声速倒置进气道激波与附面层干扰

采用标准k-ω SST湍流模型数值计算方法,针对二元高超声速倒置进气道激波与附面层严重的干扰现象,采用分流楔抑制激波与附面层干扰方法,对有无分流楔的进气道性能及流动机理特征进行了详细的研究。结果表明:采用分流楔的流动控制方法,有效抑制了激波/附面层产生的分离包对进气道内流动的干扰;提高倒置进气道的气动性能,进气道的总压恢复系数和流量捕获系数均有提高,计算模型的壁面总阻力系数得到一定程度的减小。数值计算结果表明,在分流楔尾迹中强剪切流动在一定程度上缓解了激波与附面层干扰的强度;在分流楔后缘存在稳定的横向涡,由于气流进入尾迹驻涡是来自附面层外的总压较高的高能流体,提高了附面层的抗逆压能力;由于尾迹驻涡的存在使得分离涡没有向弹体周向扩散,减小了阻力。该方法实现了对高超声速倒置进气道激波/附面层干扰的抑制,揭示了其抑制的机理。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

突扩管中受限射流的计算方法研究

突扩型燃烧室的结构简单可靠,近年来已引起广泛注意。为了计算时均流场和湍动流场,必须考虑迴流区及附面层全方位的压力梯度,给计算方法带来困难。 本文由时均的N—S方程和—ε湍流模型出发,将Patankar解耦技术推广到可压缩流,采用牛顿迭代法和稀疏矩阵技术,由连续和动量方程直接求解压力和速度场,继而获得温     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。