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

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

腔内平板对空腔自激励振荡的影响及预估振荡频率方程的改进

运用全隐式无分裂方法,紊流模型采用修正的B-L代数模型,数值模拟了超音速流通过腔内存在平板的空腔流动.研究了腔内有不同长度平板存在的空腔流动的剪切层自持性振荡,并与无平板存在时的现象进行类比.腔内平板将剪切层与空腔隔开,减弱了剪切层振荡与空腔压强脉动的耦合作用;随着腔内平板长度的增加,剪切层振荡的程度有所减小.基于特征长度可变的思想,本文改进了Heller的频率预估方程,改进后的频率方程能更好地预估空腔的振荡频率.与Heller的方程相比,修改后的方程预估所得到的频率值更接近于闭式空腔声学模型方程的预估结果;各个方程预估空腔振荡频率值的比较进一步验证了改进后的方程的正确性及合理性.     

空腔流动特性数值模拟研究

本文采用数值方法模拟研究了空腔的长深比（L／D）对流动特性的影响。控制方程采用Navier-Stokes方程，计算方法采用有限体积法，湍流模型采用Baldwin-Lomax模型。通过分析三种典型空腔的流线图谱、压力云图、底部压力分布曲线，得到空腔流动特性与其几何参数和来流参数间的相应关系。从数值模拟的结果来看开式空腔流动的空腔内压力分布均匀，流场结构相对稳定，有利于内埋式武器安全平稳地分离。     

背压脉动引起的超声速进气道流固耦合振动分析

冲压发动机燃烧室压力脉动对进气道形成背压脉动,影响了进气道内流动的稳定性。采用数值仿真方法,研究了背压脉动作用下超声速进气道内非定常流动与局部弹性壁板的流固耦合振动问题,分析了壁板阻尼对进气道流固耦合振动问题的影响。结果表明,弹性壁板振动形成了"动态空腔"效应,大幅降低了激波串振荡幅值,增加了非定常流动稳定性。其中,无阻尼弹性壁板模型的激波串振荡幅值下降为刚性壁板模型结果的45.6%。通过增加壁板阻尼的方式,能够增强壁板振动的"动态空腔"效应。当刚度比例阻尼系数为10~(-4)时,激波串振荡幅值降为刚性壁板模型结果的43.7%。     

改进的VOF方法对气液两相流振荡流动和传热计算的影响

在VOF(volume of fluid)模型的基础上,引入level set函数,形成CLSVOF(coupled level set volume of fluid)模型对气液两相流的振荡流动和传热过程进行数值模拟,与只采用VOF模型模拟的结果进行对比分析,并与实验结果相比较.结果表明:无论是针对矩形空腔中两相流的振荡流动和传热,还是活塞冷却油腔中机油的振荡冷却过程,采用CLSVOF模型都可以更为准确地模拟气液两相流的振荡流动规律,传热系数计算值也与实验结果保持了更高的一致性.对于活塞冷却油腔中机油的振荡冷却计算,随着转速即雷诺数的增大,传热系数计算值的误差开始增大.在极限转速3000r/min时, CLSVOF模型的误差为3.8%,VOF模型的误差为7.5%, CLSVOF模型的误差增幅稳定在0.6%左右.CLSVOF模型的误差较小且误差增幅稳定,说明其计算准确性更高.      

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

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

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

小间距梯形扰流柱通道内的流动换热数值计算

为了深入了解涡轮叶片尾缘扰流柱区域的流动结构和换热情况,通过数值计算研究了小间距梯形扰流柱通道内的流动和换热.数值计算结果表明:①数值计算得到的端壁静压分布、弦向出流量分布、总压系数与实验结果符合良好,端壁换热分布规律总体上与实验结果一致.②扰流柱区域的流动方向基本为弦向.沿弦向湍流度是逐渐升高的,靠近弦向出口扰流柱后的湍流度比较低.数值计算结果对涡轮叶片内部冷却计算具有重要的参考价值.      

长深比对空腔流动与声学特性的影响分析

长深比(L/D)是影响空腔流动的最重要外形参数。为了对空腔的流动机理有更深入的认识,采用基于Realizable k-ε湍流模型的脱体涡模拟(DES)数值方法,研究空腔L/D的变化对空腔流动与声学特性的影响。选取的来流马赫数为0.85,基于空腔长度的雷诺数为6.84×106。结果表明:随着L/D的增加,空腔前部的压强系数逐渐降低,后部的压强系数逐渐升高,导致空腔内的静压梯度逐渐增大;整个空腔内的压强脉动强度减弱,噪声环境有所改善。研究结果可为空腔流动类型的划分提供一定的依据。     

用一种非线性k-ε模型计算突扩台阶后的流动

对一种新的非线性k-ε模型的无量纲系数进行了重新拟合.用新的系数计算了槽道流动、方截面管流和等应变率,平均流场无旋流动的正应力分布.计算结果与实验资料符合很好.将该模型并入常规的k-ε模型的计算程序中,并对突扩台阶后的湍流流动进行了数值计算,同线性k-ε模型的计算结果和实验结果进行了比较.结果表明,这种非线性k-ε模型同线性模型相比有较大的改进,同已有的其它非线性模型相比,也有一定的优点.     

内埋式弹舱流场特性及武器分离特性改进措施

为改善内埋弹舱的流场特性以及内埋武器分离特性,采用在弹舱前缘悬细金属条的方法对弹舱流场进行流动控制,并在高速风洞中进行了试验研究。通过分析舱底静态压力试验结果以及脉动压力试验结果,研究了武器模型处于不同分离位置时流动控制对弹舱流场特性的影响;通过测量武器模型力和力矩,研究了流动控制对武器分离特性的影响。研究结果表明：武器模型处于不同分离位置时,该流动控制方法对弹舱底部静态压力分布以及总声压级分布的影响是相似的;当弹舱的流场类型为过渡／闭式穴流动时,采用该流动控制措施能有效降低舱内的静态压力梯度,并能有效改善武器的分离特性;当弹舱的流场类型为开式穴流动时,采用该流动控制措施能有效抑制舱内产生的气动噪声。     

亚声速武器舱空腔流动压力特性及其控制方法

在高速风洞中对某武器舱舱门打开状态下,舱内的静压分布和脉动压力特性进行了试验,并对空腔前缘多孔扰流板和前缘吹气两种流动控制方法分别进行了变参数研究。试验马赫数为0.75,武器舱长度为300 mm,长深比为3,宽深比为1。试验结果表明该武器舱基本构型顶部沿着流向的压力分布比较均匀,脉动压力监测点的频谱特性表现出明显的模态特征,为开式空腔流动类型,而且在所研究的范围内受飞机迎角变化的影响很小。在典型的飞机迎角下,前缘多孔扰流板流动控制方法可以很好地改善武器舱内部流场的稳态特性并降低舱内的脉动压力;多孔扰流板的安装高度、展向宽度为流动控制效果的主要影响参数;在所研究的参数范围内,综合考虑多孔扰流板对脉动压力宽频和单频特性的控制效果,高安装高度、短展向宽度的参数组合形式最优,并且在飞机常用的迎角范围内具有较好的流动控制效果。在典型的飞机迎角下,前缘吹气流动控制方法也可以很好地改善武器舱内部流场的稳态特性并降低舱内的脉动压力;吹气位置、吹气流量为流动控制效果的主要影响参数;在所研究的参数范围内,综合考虑前缘吹气对脉动压力宽频和单频特性的控制效果,大的吹气流量、短的展向宽度的参数组合形式最优,并且在飞机常用的迎角范围内具有较好的流动控制效果。     

Trajectory and attitude deviations for internal store separation due to unsteady and quasi-steady test method

The deviations of trajectory and attitude angle for internal store separation are evaluated by two wind tunnel test methods. One is the Freedrop Test(FDT), which is known as unsteady and time-dependent method of scaled model. The other is the Captive Trajectory System(CTS) test,which is usually regarded as a quasi-steady and time-averaged test technology. The result shows that there is a streamwise adverse pressure gradient on the cavity resulting in a nose-up pitching moment coefficient(>0) ...     

Aerodynamic characteristics of store during lateral jet assisted separation from cavity using free drop technique

Aimed at the problem of store separation from internal cavity, this paper innovatively puts forward a separation scheme of using lateral jet to assist store safe separation. The jet ensures that the store is continuously subjected to down head moment during separation, so as to ensure safe separation. The wind tunnel free drop test technique with lateral jet is established, which can ensure that there is no support interference in the motion process of the store and more truly simulate the motio...     

零质量射流对开式空腔气动噪声抑制效果分析

高速开式空腔流动,腔内存在较复杂的流场结构,在一定条件下腔内存在较为严重的压力、速度等脉动,诱发强烈噪声,声压级(SPL)甚至可高达170 dB,对腔内的储藏物与空腔自身结构安全构成较大威胁,因此开式空腔噪声抑制方法成为争相研究的热点.为此,对跨、超声速流动条件（马赫数Ma=0.9,1.5）下有、无零质量射流时开式空腔...     

篦齿封严泄漏特性的实验

为了研究工程实际尺寸下篦齿封严泄漏特性,设计并搭建了篦齿封严实验台.在进出口压比为1.1~2.0的条件下,探究了节流间隙、齿数、齿腔深度和宽度等对篦齿封严泄漏量和密封腔内压力分布的影响.结果表明:泄漏系数随着压比和节流间隙的提高而增加,但压比达到1.6后,增加趋势逐渐变缓;增加齿数能显著提高篦齿的封严性能,但随着齿数变多,效果逐渐减弱;较浅和较宽的齿腔对于提高篦齿的封严性能更有利;进出口压比、节流间隙和齿腔深度对于篦齿腔内压力系数的分布几乎没有影响;篦齿腔内沿程压力不断降低,但在各齿腔内降低的程度逐渐减小;在篦齿封严中,越接近入口的齿腔对密封性能提高的贡献越大,对封严性能的优劣起着决定性作用.      

Maxwell气固相互作用模型对稀薄高超声速凹腔绕流流动特征和热环境的影响

针对高空航天飞机等再入飞行器表面缺陷或防热瓦缝隙导致的局部压力过高和气动加热问题,采用直接模拟Monte Carlo （DSMC）方法研究了飞行高度为80 km的稀薄流区高超声速凹腔绕流问题,考虑气固相互作用（GSI）模型对凹腔流场特征和表面压力、热流的影响。结果表明：稀薄流条件（80 km）下,GSI为完全漫反射时,在凹腔前缘分离的剪切层再次附着在后缘,在凹腔内部形成一个充满腔体的单涡结构;随着GSI从完全漫反射向镜面反射变化,气体与凹腔表面之间的切向动量交换减弱,即黏性剪切作用减弱,外部气流被卷入凹腔的程度减弱,导致涡结构不断减小直至消失,凹腔底部逐渐出现所谓的"死水区"。与完全漫反射相比,镜面反射或近镜面反射会导致凹腔上游侧面的峰值压力和峰值热流以及下游侧面的峰值压力剧烈增大,在飞行器设计中,应特别留意上述表面的压力载荷和热载荷。     

转静盘腔流动瞬态响应特性研究

以航空发动机转静盘腔内非稳态流动为研究对象,基于用户自定义函数编程设置进口压力按斜坡函数变化时的非稳态边界条件,利用数值模拟方法,在滞后效应、流场演化和物性参数方面,研究了不同无量纲跃升幅值下转静盘腔流动瞬态响应机理。结果表明:无量纲跃升幅值由1.05增至1.2,盘腔内部的当地响应时间延长了81.1%,相对总温的超调量增加了157.0%,其峰值时间增加了31.3%,且进口压力波动的传播方向总是由进口指向出口;瞬态过程中,转静盘腔的容积效应和旋转效应共同作用使出口流量明显滞后进口流量,且其滞留质量流量与相对总温的超调现象密切相关;冲击射流、旋转贴壁射流和静子壁面侧的涡系构成了盘腔流场演化的主体,无量纲跃升幅值越大,流场演化过程越剧烈,响应时间越长;声速和热扩散率与盘腔响应速度密切相关。     

旋转盘腔瞬态响应特性的研究

为保证发动机在飞行包线内正常运转,需研究旋转盘腔的瞬态响应特性。采用1-D模型方法和计算流体力学(CFD方法对旋转盘腔进口压力突升的情况进行非稳态数值计算,所得结果与文献中的结果进行对比,提出1-D模型方法的一些缺点,并证明了CFD计算的正确性。然后用CFD方法并通过用户自定义函数(UDF编程研究了进口压力渐增、正弦变化以及盘腔尺寸对旋转盘腔流动瞬态响应特性的影响。结果表明:1-D模型的计算结果不能显示出CFD模型计算结果的一个高阶震荡;进口压力以不同方式变化,瞬态响应存在不同程度的滞后;进口压力突增和进口压力渐增响应的特征时间比进口压力正弦变化的特征时间分别增加56.0%和106.4%盘腔宽径比由0.2变化到0.39时,腔内均压变高,出口质量流量变低,特征响应时间缩短至40%当宽径比由0.2变化到0.58时,特征响应时间缩短至25%。     

环形中心钝体驻涡燃烧室驻涡腔有无喷射的对比

采用三维雷诺平均Navier-Stokes(N-S)方程、renormalization group (RNG) k-ε湍流模型和标准壁面函数对驻涡腔有无喷射的环形中心钝体驻涡燃烧室的冷态流场进行了数值仿真,分析了驻涡腔有无喷射对环形中心钝体驻涡燃烧室涡系结构、驻涡腔流动参数和燃烧室总体性能的影响.结果表明:相比于无喷射时,驻涡腔添加喷射可以使驻涡腔内形成相对稳定的双驻涡结构;驻涡腔喷射的存在使得环形中心钝体驻涡燃烧室出口截面总压损失系数降低约9.2%,并能提高驻涡腔内的平均气流参数;驻涡腔喷射对环形中心钝体驻涡燃烧室出口截面流动参数沿流道高度方向的变化趋势影响不大.