三维舵面绕流的N-S方程计算研究

基于粘性非结构混合网格,采用格心格式的有限体积法,应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟。操纵面铰链力矩的计算结果与风洞试验结果的对比分析表明,计算方法对粘性流动具有良好的适应性。对舵面铰链力矩随迎角、马赫数、雷诺数变化的规律也进行了分析研究。     

湍流燃烧模型对氢燃料超燃室流场模拟的影响

采用化学平衡的假定概率密度函数(PDF)模型和火焰面模型计算了德国宇航研究中心的超燃室反应流,计算结果与有限速率反应模型的和实验的结果进行了对比.使用有限体积法离散Favre平均的N-S方程,湍流模型采用k-ε模型.研究表明:(1)有限速率反应模型在喷氢孔近场,化学平衡的假定PDF模型在喷氢孔远场不能准确捕捉流场的细致结构,而火焰面模型对全流场预测较好,后两种模型的计算时间较有限速率反应模型节省约38%;(2)超燃室内湍流和燃烧相互作用不可忽略,从预测精度和计算效率来看,火焰面模型有较好的工程应用前景.      

雷诺数变化对翼型边界层发展及失速特性的影响

翼型失速及其边界层发展是飞行器设计中的基础科学问题,而雷诺数变化对其影响很大。针对后缘失速翼型,采用Menter k-ω SST模型及耦合扰动放大因子输运方程的转捩模型,进行雷诺数变化对层流-湍流转捩边界层特性和失速特性的影响分析。结果表明:雷诺数增大时,对于转捩边界层,当地涡量雷诺数增大,转捩前移且分离泡减小,流动能量耗散减小,翼型整体表面剪切效应增强,动能更充沛,流动自持能力增强,压力分布可以维持较长距离的梯度抵抗分离能力增强;因此雷诺数增大使翼型失速迎角提高、升力系数增加。     

基于湍动能输运的一方程转捩模型

传统转捩模型的构造十分依赖转捩机理的发现和理解,对机理的模化也存在误差。本文基于尺度自适应的KDO湍流模型,以输运变量r=μt/μ重新标定模型参数,使模型整体达到了流动结构自适应,可实现湍流/转捩一体化计算。根据长度尺度的不同,可细分为KDO-tran和CKDO-tran模型。CFD计算评估了经典的T3A,T3B平板边界层旁路转捩,T3A-平板边界层自然转捩,Aero-A翼型分离泡转捩,DLR-F5,6:1椭球横流转捩及可压缩平板边界层、尖锥转捩。该模型无法精确捕捉层流-湍流的过渡,但能准确捕捉转捩的起始位置(Transition onset),在高Re数、复杂流动的表现尤其突出。该模型捕捉转捩的机制在于流动结构的自适应和湍流输运特性的保存,不引入任何转捩机理却能捕捉多种类型的转捩现象,具有较好的适用性。     

轴流压气机插板式进气畸变数值仿真

以航空发动机稳定性评价体系中常用的插板实验为背景,分别进行低速、高亚声速、跨声速轴流压气机与插板畸变发生器的耦合数值仿真,研究插板畸变条件下压气机的流场特性以及不同类型压气机在插板畸变影响下的失速起始机制。研究发现:转子进口截面的周向流动是影响低速压气机稳定边界的主要因素;插板角涡对转子叶顶泄漏涡的扰动是影响高亚声速压气机稳定边界的重要原因;畸变气流影响下的叶片通道内激波强度和位置的变化是影响跨声速压气机稳定边界的关键因素。      

涡轮导向叶片综合冷却效率实验研究

为了研究流动参数对涡轮导向叶片综合冷却效率的影响,采用红外热像仪对叶片表面的温度分布进行了测量,得到了叶片的综合冷却效率随流量比、温比、主流进口雷诺数和湍流度的变化规律。实验过程中,次流与主流的流量比分别为0.15,0.18,0.20,0.22和0.24;主次流温比分别为1.4,1.7,1.93和2.2;主流进口雷诺数分别为1.0×10~5,1.1×10~5,1.2×10~5,1.3×10~5和1.4×10~5;主流进口湍流度分别为0.506%,8.156%,14.92%。结果表明,综合冷却效率在前缘处最低,沿流向逐渐升高;增大流量比会显著提高叶片的综合冷却效率,在温比为1.93时,流量比由0.15增大至0.24,综合冷却效率平均增加29.3%;温比和主流进口湍流度的增大均不利于综合冷却效率的提升,流量比为0.20时,温比由1.4增大至2.2,综合冷却效率平均下降46.5%,湍流度由0.506%增大至14.92%,综合冷却效率平均降低15.5%;主流进口雷诺数对叶片综合冷却效率的影响很小。     

小型回流燃烧室多工况碳黑分布预测

为了研究小型回流燃烧室不同条件下的工作状态,并分析其产生积碳现象的原因,采用k-ε湍流模型、EDC湍流燃烧等模型对该燃烧室多个工况进行数值计算。结果表明,燃烧室性能参数计算值与台架数据较吻合,各个状态下参数相对误差均小于1.37%。根据各工况的温度分布的变化,发现燃烧室火焰随工况降低逐渐向火焰筒头部收缩。计算了各工况蒸发管内的燃油蒸发率,发现蒸发率在低工况时明显降低。根据出口温度分布情况,计算出口温度分布系数OTDF为0.35。分析碳黑粒子的分布,发现积碳现象主要在火焰筒头部,其原因主要是燃烧不充分和冷却气膜的缺失。     

高湍流度时全气膜涡轮叶片表面冷却和换热特性的实验研究

为获得高主流湍流度时全气膜涡轮叶片表面的冷却和换热特性,在跨声速风洞中实验研究了质量流量比(MFR)和主流雷诺数(Re)对叶片表面气膜冷却效率和换热系数比的影响。在叶片前缘布置了5排圆形孔,在吸力面和压力面分别布置了3排和6排圆形孔,实验结果由嵌入在叶片中截面的热电偶测得。实验中基于弦长的主流雷诺数的范围为3.0×105～9.0×105,叶栅出口马赫数Ma为0.8, MFR的范围是5.5%～12.5%,主流湍流度Tu为14.7%。实验结果表明:主流雷诺数升高显著增强了叶片表面的换热,使层流边界层到湍流边界层的转捩位置提前。对于吸力面S/C0.2的区域(S/C为当地弧长与弦长之比),气膜冷却效率受MFR影响明显,当MFR大于7.7%时提高MFR会导致气膜冷却效率降低;该区域的换热系数比在中低雷诺数时受MFR影响较小,在高雷诺数时随MFR升高而升高。压力面S/C-0.7区域的气膜冷却效率随MFR升高而升高,-0.7S/C-0.4区域的气膜冷却效率受MFR影响较小,对于整个压力面而言,MFR升高提高了叶片表面的换热系数。相对于叶片其它区域,压力面后半段区域和吸力面的气膜冷却效率受雷诺数影响较大。     

三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。     

涡流发生器周向相对位置和高度对高负荷风扇性能的影响

为了研究涡流发生器周向相对位置和高度对高负荷风扇性能的影响,根据风扇的流动特点,设计了在第二级静子叶根入口前加涡流发生器的流动控制方案,并以此为基础提出了多种不同周向位置和高度的涡流发生器方案,通过计算对采取各种方案下的流场进行了分析。研究表明,涡流发生器对风扇第二级静子角区气流分离有较好的控制作用;涡流发生器的周向位置对第二级静子角区气流分离和损失的影响较大,采取方案C时可以更好地抑制角区气流分离,减少局部损失;涡流发生器高度过高会使静子压力面出现不同程度的低速区,同时也会引起静子通道内局部损失增加,在所研究的范围内,当涡流发生器高度降低1%叶高时,其对吸力面角区分离的控制效果更加明显。