突扩燃烧室流场RANS与LES模拟

以中心突扩燃烧室为模拟对象,分别从旋涡运动、空间速度分布和时间速度脉动等方面对RANS模拟结果、LES模拟结果和试验结果进行了比较分析。表明RANS方法没有模拟出燃烧室中的旋涡运动,与试验结果相差甚远;而LES方法则比较好地模拟出了燃烧室中旋涡运动和速度脉动。得出LES方法是对突扩燃烧室内流场数值研究的有效方法。     

导弹垂直发射燃气排导系统两相流场数值分析

根据燃气排导系统内两相流场的分布特点,利用可压缩雷诺 N-S平均方程对流场进行计算。湍流模型采用标准 k-ε二方程模型,采用颗粒轨道模型对颗粒相进行模拟,利用 FLUENT就发射高度及燃烧室压强对三维流场参数分布的影响进行了重点分析,为优化发射结构设计提供了重要的理论依据。     

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

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

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

湍流涡黏性模型方程中相位差的测量

基于理论上湍流相干结构复涡黏性模型对涡黏性系数的分析,应用热线测速技术测量了平板湍流边界层多尺度相干结构动力学方程中非相干结构成分对相干结构贡献的雷诺应力分量与相干结构流向速度流向变形率之间的相位差.分析了湍流边界层相干结构猝发过程中,非相干结构成分对相干结构贡献的雷诺应力分量与相干结构流向速度流向变形率之间的相位差沿湍流边界层法向的变化规律,为建立更加符合实际的湍流模型提供了实验依据.     

改变结构参数对旋流喷嘴内外两相流场的影响

针对一种超细粉体旋流分散喷嘴,在不同进气角度和出口锥角的条件下,采用雷诺应力湍流模型(RSM)并考虑气体的可压缩性,对喷嘴内部及出口附近的强旋、跨声速流场进行了模拟,同时结合离散相模型研究了1~10μm粒子的运动轨迹和质量浓度分布规律.结果显示:改变旋流喷嘴的进气角度对旋流强度和湍流强度影响不大,而出口锥角改变时,喷嘴内的旋流强度有明显变化;随着颗粒密度或出口锥角增大,或进气角度减小,旋流喷嘴内粒子的临界逃逸粒径相应减小;在喷嘴后方,颗粒相达到较好的扩散均匀性所需的掺混距离约为20倍喷嘴出口直径.      

高速发动机油膜惯性对活塞裙润滑影响

利用一种新的迭代算法 ,分析了高速发动机油膜惯性对活塞裙润滑特性的影响 .该迭代方法是利用有限元法和差分法交替求解Navier Stocks方程和雷诺方程 ,据此导出了适用于高速发动机活塞裙的润滑计算的混和润滑雷诺模型 .新的模型借助惯性系数 ,引入了油膜惯性项 ;同时给出求解含有油膜惯性项的迭代步骤和有限元表达式。计算结果表明 :随着惯性系数和活塞裙的长径比的不同 ,油膜惯性会对油膜的摩擦力、压力和承载力产生不同的影响 ,这种影响对承载力尤为明显 .该润滑模型也可用于中、低速发动机的活塞裙润滑计算以及不计入惯性项 (惯性系数置为零 )的某些润滑问题求解 .     

基于S-A与SSG/LRR-ω两种湍流模型的CHN-T1标模计算与分析

通过求解雷诺平均Navier-Stokes方程,对第一届AeCW-1研讨会上要求的所有工况进行了计算,着重研究S-A湍流模型和雷诺应力模型在CHN-T1标模上的表现,评估两种湍流模型的网格收敛性,对抖振特性的模拟能力以及对雷诺数的敏感性。根据与实验数据的对比分析,得出主要结论如下:两种湍流模型都能趋向于网格收敛,雷诺应力模型的网格收敛更快。在较大迎角的工况,两种湍流模型对机翼上激波诱导的分离流有较大的差别。CHN-T1标模的尾撑,改变了飞机平尾的压力分布,并且在飞机机尾上产生复杂的三维分离流动,对俯仰力矩系数的准确计算产生了较大的影响。考虑静气弹效应,两种湍流模型计算的力矩系数比不考虑静气弹效应更大,雷诺应力模型的偏差更大。在高雷诺数工况,两种湍流模型的计算结果趋于一致。     

一种关于静压气体轴承节流孔系数的计算方法

基于层流边界层方程的分离变量算法和雷诺方程的解析算法,提出了一种关于单节流孔静压气体止推轴承的节流孔系数的计算方法。该方法通过比较层流边界层方程计算获得的气体轴承的质量流量和雷诺方程计算获得的质量流量计算获得了节流孔系数。将计算获得的节流孔系数和节流孔系数为常数0.8代入单节流孔气体止推轴承的雷诺方程中,计算获得的承载力与分离变量算法求解层流边界层方程获得的承载力进行对比,可以发现,相对于采用节流孔系数为0.8来说, 采用该计算的节流孔系数求解雷诺方程的承载力与分离变量算法求解获得的承载力结果精度最大提高了8%。从而验证了该计算节流孔系数方法的正确性。      

基于后缘小翼的旋翼翼型动态失速控制分析

针对后缘小翼(TEF)的典型运动参数对旋翼翼型动态失速特性的控制进行了研究。发展了一套适用于带有后缘小翼控制的旋翼翼型非定常流动特性模拟的高效、高精度CFD方法。通过求解Poisson方程生成围绕旋翼翼型的黏性贴体和正交网格,为保证后缘小翼附近的网格生成质量,建立了基于翼型点重构的方法来描述后缘小翼的偏转运动;为克服变形网格方法可能导致网格畸变的不足,发展了一套适用于带有后缘小翼控制的旋翼翼型运动嵌套网格方法。基于非定常雷诺平均Navier-Stokes(URANS)方程、双时间法、Spalart-Allmaras(S-A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,发展了旋翼翼型非定常气动特性分析的高精度数值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隐式时间推进方法及并行技术提高计算效率。以有试验结果验证的HH-02翼型和SC1095翼型为算例,精确捕捉了动态失速状态下的气动力迟滞效应,验证了本文方法的有效性。着重针对SC1095旋翼翼型的动态失速状态开展后缘小翼的控制分析,提出了可以体现翼型升力、阻力及力矩综合特性的关系式Po和Pc,揭示了后缘小翼振荡频率、相位差和偏转幅值对动态失速特性影响的规律。研究结果表明:当后缘小翼偏转的相对运动频率为1.0,且小翼运动规律与翼型振荡规律之间的相位差为0°时,后缘小翼能够更好地抑制翼型动态失速现象;在此状态下,当偏转幅值为10°时,SC1095翼型最大阻力系数和最大力矩系数可以分别降低19%和27%。