用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。     

基于CFD降阶模型的阵风减缓主动控制研究

飞行器飞行时会受到大气紊流的影响,降低飞行品质。阵风减缓控制是改善飞行器飞行性能的关键技术。现有的阵风响应分析多以离散阵风为研究对象,对更加真实描述大气紊流的连续型阵风时域分析关注较少。采用成形滤波器方法将频域形式给出的大气紊流信号转换为时域信号。在跨声速区域内,利用系统辨识技术,基于计算流体力学(CFD)方法建立阵风激励下的气动载荷状态空间降阶模型(ROM)。为方便控制器设计,借助平衡模态法进行模型的进一步降阶。使用模型预测控制(MPC)算法通过控制操纵面偏转实现阵风减缓主动控制。以AGARD445.6标模作为仿真算例,验证基于ROM设计的阵风减缓控制律的有效性。仿真结果表明,在跨声速飞行状态下,模型预测控制器能够在满足操纵面偏转范围的约束下,对连续阵风激励下的翼根弯矩输出进行有效抑制。     

发动机进气帽罩防冰热载荷的数值模拟研究

对某发动机进气帽罩防冰热载荷的计算方法进行研究。利用FLUENT软件对帽罩周围的流场进行模拟,同时采用FLUENT中的离散相模型计算帽罩周围的水滴运动轨迹,得到帽罩表面的水滴撞击效率;利用数值计算结果及流动换热经验关系式分析帽罩表面的各项热流,获得该发动机帽罩的防冰热载荷。计算中还比较了不同湍流模型对防冰热载荷的影响。     

湍流度和吹风比对叶片前缘双出口孔射流冷却效率影响

为了探讨湍流度对一种新型气膜孔射流气膜冷却影响,利用商业软件提供的有限体积法求解N-S方程,对湍流度分别为0.4%、10%和20%时的双出口孔射流冷却效率进行数值模拟。吹风比变化范围为0.5到2.0。首先将圆柱孔射流冷却效率计算结果和实验数据进行了对比,二者吻合较好。计算得到了双出口孔射流气膜冷却下的流场、径向平均冷却效率。结果表明,湍流度和吹风比对冷却效率都有较大影响。湍流度为0.4%和10%时,最高冷却效率在吹风比1.0时获得;湍流度30%时,最高冷却效率在吹风比2.0时获得。吹风比0.5和1.0时,冷却效率随着湍流度的增加而降低;吹风比2.0时,冷却效率随着湍流度的增加而增加。     

结构因素对离心通风器性能影响的数值研究

在高转速下,航空发动机离心通风器内腔存在强烈的湍流流动,而稳定工作后通风器内腔中漩涡呈周期性产生和淹灭。小尺寸颗粒受流体微团湍流随机脉动的影响,产生相对于平均流的随机脉动运动。为了研究通风孔偏心距和辐板顶圆半径对通风器性能的影响,采用离散相模型(DPM)模拟通风器内颗粒的运动轨迹。并应用随机游走(DRW)模型模拟连续相湍流瞬时速度脉动对颗粒轨迹的影响,采用随机涡寿命模型确定随机追踪模型的积分时间。结果表明:辐板顶圆半径的增大有助于提高通风器的分离效率,同时也增大了腔内流通阻力;通风孔偏心距对减小通风阻力的作用明显,但降低了离心通风器的分离效率。     

旋转封严篦齿风阻温升的试验研究与数值分析

选取具有典型封严结构的真实尺寸台阶型封严篦齿,通过高速旋转条件下的风阻温升特性试验,进行了进出口压比(1.05~2.8)、封严间隙与齿尖厚度比(相对封严间隙)(2.4~4.0)、雷诺数(1900~28000)及转速(0~12000r/min)等参数对篦齿泄漏流风阻温升的影响研究.并基于试验条件,采用renormalization group(RNG)k-ε湍流模型,对同一篦齿的风阻温升特性进行了数值计算.结果表明:篦齿泄漏流温升随进出口压比、雷诺数、封严间隙与齿尖厚度比的增加而减小,但随转速的增加而增大.数值计算结果与试验数据吻合 良好,台阶型封严篦齿第1道和第2道封严齿内的风阻温升最为显著,前两道封严齿内的温升占总温升的60%. 根据试验结果,提出了计算台阶型篦齿泄漏流温升的经验公式,其计算结果与试验数据具有较好的符合性.      

Gamma-Theta经验转捩模型在DES中的应用

与RANS湍流模型相比,DES湍流模型的边界层外自由流部分存在LES区,该区域的涡粘性会有所降低。为了研究了这一特征对γ-Reθt转捩模型的影响,以平板实验为算例,选取基于k-w SST湍流模型的DES模型和DDES模型,检验γ-Reθt转捩模型在DES模型和DDES模型中的表现。结果表明：DES模型和DDES模型使转捩位置推迟,SST—DES模型由于对网格过度敏感而没有触发转捩,SST—DDES模型使转捩位置后移较少;DDES模型比DES模型更适合γ-Reθt转捩模型;转捩经验关系需要进行修正后才能将γ-Reθt应用于DDES模型。     

考虑模型抖振力跨向不完全相关性效应的气动导纳识别

首先通过考虑节段模型上抖振力跨向不完全相关性效应推导了节段模型风洞试验中作用在模型断面上的分布抖振力谱和由底支式天平测到的模型总抖振力谱之间的关系。然后,以准平板断面为例,进行了格栅湍流场节段模型的测力和同步测压试验,获得了模型总抖振力谱以及模型抖振力跨向相关性函数。接着,采用等效导纳法以及抖振力自谱和抖振力脉动风速交叉谱综合残量最小二乘法分别识别了准平板节段模型等效气动导纳和六分量气动导纳,讨论了模型抖振力跨向不完全相关性效应对气动导纳识别结果的影响,并把识别得到的准平板断面气动导纳与平板断面气动导纳的理论结果——Sears函数进行了比较。结果表明：忽略抖振力跨向不完全相关性效应（即假设节段模型分布抖振力沿跨向完全相关）而直接采用平均抖振力作为断面上分布抖振力的传统方法会导致气动导纳识别结果偏小,并且,其偏小程度会随着频率的增加而增加;此外,相对于抖振升力和扭矩相关的气动导纳分量而言,由于抖振阻力的跨向相关性比抖振升力和扭矩的跨向相关性显得更弱,因此抖振力跨向不完全相关性效应对阻力相关气动导纳分量识别结果的影响更大;利用识别得到的六分量气动导纳反算的作用在模型上的分布抖振力谱与试验中实测结果非常接近,经抖振力跨向不完全相关性效应修正后的竖向脉动风速对应的升力和扭矩气动导纳分量的识别结果与Sears函数也比较接近,从而验证了用于六分量气动导纳识别的自谱-交叉谱综合最小二乘法的可靠性。