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

采用了直接过滤的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-ε模型相比,新模型计算结果与实验更加接近.     

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

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

压气机转子叶片的气动弹性数值模拟

采用了交替迭代算法,对某压气机转子叶片进行了气动弹性数值模拟.自行开发了基于有限元的结构求解器用于结构动力学求解,引用他人开发的非定常流体求解器用于气动力的求解.结构求解器提供叶片的表面位移给流体求解器以改变流场,流体求解器提供气动载荷给结构求解器来计算叶片的变形,界面处理系统在这两个求解器之间进行信息传递.算例表明,这种交替迭代算法在气动弹性数值模拟中是可行的,可以得到叶片的瞬态响应,从而判断叶片是否发生颤振.      

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

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

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

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

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

高超声速飞行器瞬态热试验

为了进行高超声速飞行器热防护系统的初步设计和数值计算的验证,设计开发了高超声速飞行器瞬态气动加热地面试验系统及其控制软件.试验系统能够根据飞行器的飞行轨迹和外形参数加载瞬态热流,实时测出结构表面的热流值和温度,得到飞行器的表面试验热流曲线和温度曲线.试验系统采用真空舱模拟飞行环境,并为此设计了冷却床,在真空环境下能比较真实的模拟热防护系统的下表面热环境,使瞬态热试验的原理更加合理,精度进一步提高.