跨声速层流翼型的混合反设计/优化设计方法

跨声速层流翼型设计须兼顾优良的超临界特性和自然层流特性,因而对设计方法提出了更高的要求。针对现有反设计方法和直接优化设计方法的不足,发展了一种适用于跨声速层流翼型的混合反设计/优化设计方法。该方法引入了基于经验的局部流场特征作为反设计目标,翼型性能指标作为直接优化设计目标,然后加权形成了混合反设计/优化设计总目标,并同时考虑了气动和几何约束。优化算法采用基于自适应并行加点技术的代理优化,流动数值模拟采用耦合基于线性稳定性理论的e^N转捩自动判定的雷诺平均Navier-Stokes（RANS）方程求解器。针对现代中短程民用客机需求,以NPU-LSC-72613翼型为基准,开展了层流翼型减阻的混合反设计/优化设计。分别将局部目标压力分布、总阻力作为反设计和直接优化设计目标,得到了较好的优化结果,验证了方法的有效性。经过2轮优化结果显示混合反设计/优化设计总目标显著下降。所设计翼型吸力面局部压力分布与目标压力分布基本一致,总阻力下降15.5%;吸力面和压力面层流范围均大于55%倍弦长,激波强度显著减弱,说明所设计翼型同时具有优良的超临界和层流特性。将所设计翼型配置到机翼上,通过三维数值模拟进行校验,结果显示所设计跨声速层流机翼升阻比提高了6.64%,在一定升力系数范围内,气动性能均有显著提高,验证了所设计跨声速层流翼型在机翼设计中的适用性。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

地面燃气轮机单管燃烧室流量分配试验

对设计的100kW地面燃气轮机单管燃烧室空气流量分配进行了试验研究。在常压环境下采用堵孔法分别得到了旋流器、主燃孔和掺混孔的流量特性曲线,分析得到燃烧室不同结构的流量系数;试验研究了不同进口流量条件下燃烧室的流量分配,测量得到了燃烧室总压损失。研究发现:随着进口空气流量的增加,燃烧室的流量分配比例基本保持稳定,并且与燃烧室的设计空气比例基本吻合;随着燃烧室进口流量（雷诺数）的增加,旋流器、主燃孔和掺混孔的流量系数呈线性降低;随着进口流量（雷诺数）的增加,燃烧室总压损失逐渐增大;对主燃孔和掺混孔的流量特性测量中,两种测量方法得到的试验结果稍有差别,总体上看两种测量方法的试验结果较为接近。通过对比分析证明两种试验测量方法真实可靠,该研究结果可为100kW地面燃气轮机燃烧室的设计与优化提供依据。      

宽高比及表面粗糙度对矩形微尺度通道流动特性的影响

实验研究了不同宽高比矩形金属微尺度通道的流动特性,并与硅通道对比来探究表面粗糙度的影响。金属微尺度通道宽为0.4mm,通道宽高比分别为0.50、0.67、和1.00,相对表面粗糙度范围为0.1609%~0.2145%。硅通道宽度及高度都为0.4mm,相对表面粗糙度为0.00325%。实验工质为空气,实验的雷诺数范围是250~3000。实验结果表明:在粗糙微尺度通道会发生转捩提前,并且转捩雷诺数随着宽高比增大而减小。在层流区,微尺度矩形通道摩擦因数随着雷诺数增大而减小,在转捩区出现明显增大后再减小。表面粗糙度很小时,未发现转捩提前现象。      

改进的离心压气机滑移因子模型

建立一种基于Qiu模型的改进的离心压气机滑移因子模型。该模型在Qiu模型的基础上,增添二次流项来修正滑移因子,同时考虑了出口叶片角与最大叶片角的差值及相应位置。采用亚声速叶轮Krain4、跨声速叶轮Krain6和IET200叶轮对该改进模型进行CFD验证,并对模型二次流项中常系数c的取值进行讨论。结果表明:与Qiu模型及其他传统滑移因子模型相比,该改进模型的精度较高,平均相对误差在2%以下,适用于亚声速及跨声速离心压气机。      

跨声速轴流压气机的失速发展机理

为了研究跨声速轴流压气机由近失速工况发展为失速工况的动态演化机理,对跨声速轴流压气机转子NASA Rotor 37进行了多通道全三维数值模拟,着重分析了激波及前缘溢流对失速发展的影响。结果显示:在峰值效率工况下,叶片通道内存在一道斜激波;在失速工况下,斜激波演化为脱体激波。泄漏涡通过脱体激波后发生破碎,破碎的泄漏涡在向通道下游发展的过程中,受逆压梯度的影响,在通道中部形成一个明显的涡结构。在失速工况初期,由于泄漏涡的自维持现象,前缘溢流现象随着叶顶阻塞区的周期性发展而间歇性出现;随着流量的降低,通道阻塞程度逐渐增加,会出现前缘溢流一直存在的现象,这一特性可以作为流场开始急剧恶化的标志。     

R141b在矩形微尺度通道中的两相流传热特性

设计搭建水力直径分别为1mm和0.5mm的矩形微尺度通道实验台,研究了以R141b型制冷剂作为工质的两相流沸腾传热特性。实验取热流密度为1～16kW/m2、质量流速为111.1～333.3kg/(m2·s)和质量干度为0～1,分析了三者对平均传热系数的影响,探究影响换热的主导因素。结果表明:热流密度较高时,平均传热系数随热流密度增加而减小,流动换热主要受到沸腾传热的影响;当质量流速较大且热流密度较低时,平均传热系数随热流密度增加而有所增长;热流密度较低时,平均传热系数随质量流速变化明显,热流密度升高到一定值后,质量流速对平均传热系数的影响很小;当质量流速处于111.1～333.3kg/(m2·s)时,平均传热系数随质量干度的增加而减小。      

不平衡激励作用下周向加肋机匣振动能量传递机理

为了研究周向加肋机匣对整机振动的抑制机理,采用有限元法建立了加肋机匣-支承-转子耦合结构,结合振动功率流法分析了在转子不平衡激励作用下周向加肋机匣中的振动能量传递机理。 结果表明:机匣周向加肋筋诱导出的能量涡流场能够分流和耗散掉部分由转子传递至机匣的振动能量;机匣周向加肋筋改变了振动能量的传递路径,减小了机匣与支板连接处振动能量的回流现象;振动能量传递到机匣周向加肋筋后发生传递波波形的转换,携带大部分振动能量的弯曲波转换为纵波并沿加肋筋周向传递,阻断振动能量沿轴向传递至整机其他部位。 研究结果可为航空发动机结构设计以及整机减振提供参考。      

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。