超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

二维空腔超声速流动特性及形状影响数值研究

内埋式武器装载是实现战斗机隐身、超声速巡航和高机动性能的综合选择,已逐步成为新型战斗机武器装载的发展趋势。采用雷诺平均N-S方程和S-A湍流模型,模拟了矩形空腔和非矩形空腔的流动特性。通过分析空腔的流线图、压力云图和底部压力分布曲线,可知长深比是影响空腔流动类型和流动特性的关键因素,改变空腔后壁倾角、底面斜坡角度和前壁倾角对回流区结构及压力分布特性有明显影响。     

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。     

基于DES的二维和三维空腔流动特性研究

空腔流动广泛存在于航空航天工程中,对其流动特性的研究具有十分重要的工程意义。采用数值方法对比模拟了二维和三维空腔流动,控制方程采用N—S方程,空间离散采用有限体积方法,对流通量计算采用Roe格式,非定常时间离散采用双时间步长方法,湍流粘性计算采用基于SA模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致,所得的二维和三维空腔内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较,发现空腔底部压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同,得出在空腔宽度有限时,横向流动简化了流动结构．削弱了振荡幅值．但使振荡过程复杂化．其三维效应不能忽略．     

超声速湍流边界层中横向声速喷流的混合LES/RANS模拟

通过对超声速湍流边界层中横向声速喷流的计算,对一种混合大涡/雷诺平均Navier-Stokes(LES/RANS)模拟方法进行了测试,该方法采用一个依赖于到壁面的距离及当地湍流参数的混合函数结合两方程k-ω SST(shear stress transport)湍流模型和混合尺度亚格子模型来封闭湍流项.计算结果表明:混合模拟方法能够捕捉到喷流/湍流边界层相互干扰的非定常大尺度结构,且对分离区长度、壁面静压峰值和膨胀区静压分布的计算精度要高于RANS(Reynolds-averaged Navier-Stokes)方法.      

超声速空腔流动波系演化及噪声控制研究进展

空腔结构在先进飞行器和动力装置中应用十分广泛,如飞机内埋武器舱、起落架舱、导弹光学头罩、超燃冲压发动机燃烧室、推力矢量发动机喷口等。这些空腔结构在实际应用中普遍面临严峻的气动噪声问题。研究表明,空腔噪声的产生与腔内非定常流动和复杂波系结构关系密切。对超声速空腔流致噪声已开展的研究进行综述,从而为系统地认识波系结构的演化规律、噪声产生机理和噪声控制机理提供参考。首先,介绍了闭式、过渡式和开式3种不同类型超声速空腔流动的主要特点,针对气动噪声最为严重的开式空腔流动,分析了7种典型波系结构形成、传播和相互作用等规律。然后,从空腔压力振荡反馈回路的4个关键环节出发,分析了当前广泛应用的3种压力反馈模型的异同。最后,总结梳理了超声速空腔噪声控制的5条路径,指出在第5条控制路径"干扰反馈激波的传播"研究方面存在的不足。并指出了当前研究存在的部分问题,就未来研究方向提出了建议。     

三维效应对超声速湍流燃烧流动大涡模拟的影响研究

本文对空间发展Mac=0.61超声速湍流燃烧混合层流动开展了二维/三维大涡模拟（LES）研究,探讨了三维效应在超声速湍流燃烧流动的大涡模拟中的重要性。研究发现,相较于三维LES计算结果,二维LES大幅低估了超声速湍流燃烧混合层的增长率,湍动能的二维LES结果显著偏大而雷诺剪切应力显著偏小。因而,二维LES难以准确模拟超声速湍流燃烧混合层的增长率和湍流脉动特性。二维流动中漩涡拉伸、压缩机制失效,致使二维LES的漩涡结构与三维LES结果显著不同：二维LES计算所得流场呈现清晰的旋涡生成、发展、并对和破碎等演化现象,而三维LES的湍流结构呈不规则状。湍流与燃烧间存在强烈的耦合相互作用,由于二维LES难以准确模拟湍流特性,燃烧场进而在二维LES下呈现一系列特殊的、不同于真实燃烧流动情况的数值现象。漩涡运动带动燃料和氧化剂掺混,掺混界面为火焰前锋。由于漩涡结构不同,致使二维LES的火焰前锋面面积不同于三维LES结果,导致其难以获得准确的燃料效率,最终导致二维LES计算获得的统计平均火焰温度显著偏低。     

空腔流动的拉格朗日涡动力学分析

空腔流动中包含着丰富的涡结构相关物理现象,开展空腔流动中的涡动力学分析能够为相应的流动控制和降噪提供理论依据。采用直接数值模拟获得开式空腔流动的高精度流场数据,用拉格朗日拟序结构识别流场中的涡结构,并结合Q 判据和涡量通量分析,对空腔中涡的生成、脱落、对流以及撞击破裂过程展开研究,分析其动力学特性。结果表明:在一个流动周期内,拉格朗日拟序结构的卷起与脱离伴随着前缘涡的生成与脱落,向下游发展;拉格朗日拟序结构在后缘点的破裂伴随着Q 与涡的破裂;空腔内部的主涡则在整个流动周期内都维持相对稳定状态。     

脱体涡模拟方法研究及其应用

基于现有的非定常雷诺平均N-S方程的求解方法,研究与发展了脱体涡模拟方法（DES）。在紊流附面层内用RANS方法,湍流模型采用S-A模型,在其他区域结合S-A模型运用Smagorinski的LES模型,并将S-A模型中的物面耗散项进行适当的改进,使得改进的模型既充当了RANS中的S-A紊流模型,又充当了Smagorinski大涡模拟模型。为验证所发展方法的有效性,将研究的DES方法应用到一些典型流场计算中,并与已有结果或实验结果进行比较。     

空间推进方法求解抛物化Navier-Stokes方程及其验证

采用伪时间空间推进方法求解完全气体PNS方程。PNS方程采用有限体积方法离散,在推进方向上采用一阶精度迎风格式,在推进面的展向和法向采用三阶精度的MUSCL插值和AUSM类通量构造格式;在推进面上的伪时间推进采用二维LU-SGS迭代方法。通过算例证明,空间推进方法能得到正确的压力、摩阻和热流分布,且计算时间比时间迭代方法快一个数量级以上。     

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景,采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型,采用高精度WENO格式进行空间离散,并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致,来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构,并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析,揭示了流场结构特征及其时空演化规律,并获得了流场密度脉动特性。      

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

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

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。     

圆截面U形管流有限元模拟

弯曲圆管内流动具有很强的工业应用背景,管内二次涡的大小、位置及发生的频率会对流动系统的阻力及热质交换产生重大作用,研究弯曲圆管内的流动具有重要的理论意义.本文应用湍流大涡模拟方法对弯道流动进行了数值模拟.用时间推进法,给出了隐含流线迎风耗散作用的二阶完全展开Taylor-Galerkin有限元离散格式,数值计算的结果与实验结果符合较好.表明大涡模拟方法适用于边界形状复杂,存在各向异性的大尺度涡的流动仿真.     

后向台阶超声速流动的动力混合RANS／LES模拟

针对可压缩湍流流动,构建了1种动力混合RANS/LES模式。运用该模式模拟了超声速后向台阶流动,结果再现了二维后向台阶剪切层的拟序结构,计算得到了统计平均的压力分布,回流区长度与实验结果符合较好。与定系数的一般混合RNAS/LES模拟结果的对比表明,动力混合RANS/LES模式在对流场拟序结构的捕捉、流动时均压力分布方面表现较优。     

超声速细长飞行器前置小翼展开过程的非定常特性

针对大长细比飞行器外形超声速大迎角条件下前置小翼展开过程引起的非定常问题,应用结构动网格技术和基于脱体涡模拟（DES）的非定常数值模拟技术进行了研究,获得了详细的小翼展开过程流场结构非定常变化特性,并分析了小翼展开引起的法向力、纵向压心系数等气动特性随展开角度的变化规律。研究结果表明：超声速大迎角条件下前置小翼展开过程对小翼附近区域以及尾舵区域产生了强烈的干扰影响。小翼完全展开后,压心前移4.1%,降低了飞行器静稳定性,法向力系数增加15.6%,提高了气动过载,对飞行器机动性能均产生有利影响。     

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

碳氢燃气超声速剪切流动燃烧数值模拟

数值模拟碳氢富燃燃气与空气二维超声速剪切流动的湍流燃烧.NND对流项差分格式高精度捕捉激波, B-L和普朗特混合长度模型模拟湍流粘性系数, Arrhenius反应动力学公式计算化学反应速率, 模拟非平衡态燃烧.计算的壁面压力值较好地符合了实验值.分析了剪切流驻焰稳定的波系结构, 从预测的流动参数和组分浓度分布判断燃烧状况的好坏, 从而为改善优化燃烧室结构提供工程参考.      

高温高压超声速气液两相流场数值模拟

采用欧拉-拉格朗日方法对燃气-蒸汽发射动力装置内高温高压超声速燃气中横向喷雾的气液两相流进行了数值模拟研究,气相采用RK(Runge-Kutta)-AUSM+(advection upstream splitting method+)格式求解Navier-Stokes方程,液相应用颗粒轨道模型,两相之间的耦合通过在气相...