高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。     

带凹窗斜劈高速湍流气动光学效应研究

采用数值方法分析研究了带凹窗斜劈高速湍流气动光学特性受窗口深度的影响。湍流数值模拟采用耦合J-B模型的RANS/LES混合方法即DES方法,依靠双时间步LU-SGS方法开展非定常迭代求解,并且利用物理光学方法和波前重构技术等计算气动光学效应,用于分析近场波前畸变和远场光斑分布情况。模拟结果表明,在窗口较浅时,窗口前缘处产生较强激波,带来更大的波前倾斜影响,而窗口较深时,分离剪切层、旋涡等流场结构导致气动光学效应的非定常特性更为显著;高速湍流导致的气动光学效应很强,在所模拟条件下,波前倾斜、均方根光程差、峰谷值以及Strehl比分别达到69μrad、0.22λ、1.2λ和0.31,其对应的跳动值38μrad、0.04λ、0.8λ和0.43。     

二次流动对气冷涡轮叶栅流场的影响

针对典型跨声速高压涡轮叶型平面叶栅吸力面单排孔气膜冷却,采用数值模拟方法,比较分析了加入气膜冷却前后流场变化。结果表明,由于二次流动的影响,加入气膜冷却以后吸力面后部接近下壁面处没有受到冷气保护而直接暴露于主流高温燃气,在实际高压涡轮中将极大的降低叶片寿命。没有气膜冷却情况下,吸力面接近下壁面处边界层仍有可能因受到二次流动的影响发生转捩;加入气膜冷却情况下,气膜孔中心位置下游边界层由于射流和主流的相互作用将转变为湍流边界层,而由于孔间距的影响,只有射流和主流充分掺混以后才能影响到整个叶片的范围。     

试验数据中不确定参数的修正方法研究

采用Kriging方法和CFD技术对风洞试验中的不确定参数——马赫数和迎角进行修正研究。以RAE2822翼型的试验数据作为研究对象,在不确定参数的变化区域内,分别建立以升阻力和压力分布为目标的函数,并采用Kriging方法对函数值进行插值,得到CFD结果和试验值的吻合较好的不确定参数值,以达到对试验参数修正的目的。计算结果说明方法对数值模拟结果精度的提高有一定的研究和应用价值。     

一个新的可压缩性修正的k-ε模型

考虑结构可压缩性修正的影响,发展了一个同时考虑结构可压缩性修正和膨胀可压缩性修正的k-ε湍流模型,新模型包括Chang可实现性、Heinz湍流动能产生项以及Sarkar可压缩性三部分修正.新模型扩宽了以往发展的可压缩性修正模型的适用范围,适用于高超声速(M>5)复杂湍流流动中.通过对多个复杂超声速横侧射流工况的计算,验证了新模型的预测效果.与实验结果相比表明,几个工况下新模型的预测精度都显著高于标准k-ε模型.流体分离强度越大,新模型的修正效果越显著.与标准k-ε模型相比,新模型计算结果与实验更加接近.     

超声速混合层转捩过程实验研究的进展

超声速湍流机理的实验研究是一件十分困难的工作.在2000年以来,本研究小组在低噪声超声速混合层风洞研究、超声速流动精细结构测量技术研究方面取得了重要进展,这给超声速混合层湍流精细结构的研究奠定了基础.为了研究超声速混合层及其气动光学问题,在研制的超声速混合层风洞中,主要以基于纳米技术的平面激光散射技术(Nano-trace Planar Laser Scattering,简称NPLS)为基础,研究了几种对流马赫数的超声速混合层从层流到湍流转捩过程K-H不稳定涡的空间结构,以及K-H不稳定涡的空间结构随着时间的发展过程.实验结果清晰地反映了湍流混合的不稳定性与转捩的精细结构,以及转捩过程的展向精细结构.     

圆转矩形喷管射流掺混特性试验研究

本文采用热线风速仪对圆转矩形收敛喷管的射流流场与掺混特性进行了试验研究,测量了喷口下游不同截面上,射流宽、窄对称面上射流的速度、雷诺剪应力、湍流强度,得到其沿径向和轴向的分布规律。发现在径向方向上,雷诺剪应力先增大后减小,最大值出现在射流与外流的交界面附近,湍流强度则逐渐减小;在轴向上,雷诺剪应力在近喷口区存在小幅波动,而后逐渐减小,湍流度则先略有增大而后变化逐渐减慢。射流特性在宽对称面与窄对称面上的分布规律相同。     

环截面管中的爆燃转爆震实验

为了实现在环截面爆震管封闭端同时点火,采用了先在预燃室中点火,再通过八路支管将火焰引入到爆震管封闭端,各个支管长度一样,因此可看成在环管封闭端多点同时点火.实验结果显示,在爆震管封闭端火焰汇聚形成激波点火,可在环管封闭端附近直接产生爆震波.在爆震管中安放扰流片时,流动阻力使爆震速度偏低,无扰流片时,爆震速度接近理论值.不加扰流片当量比在1.0左右,火焰速度会出现爆燃或爆震两种随机状态,当量比到1.1时可避免这种状况.      

巡飞器无舵偏情况下的气动数值分析

采用SSTk-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程实现对巡飞器无舵偏情况下湍流流场的数值模拟。对巡飞器低雷诺数绕流流场进行了计算,分析了其在无舵偏情况下气动特性随迎角、侧滑角的变化情况。仿真结果表明,巡飞器采用充气式机翼后,具有很好的失速性能与较好的大迎角稳定性。     

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

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