基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

热/力耦合作用下基于应力分析的冰破坏准则

根据飞机热除冰的物理过程,考虑外部空气动力和蒙皮表面加热的作用,建立了NACA 0012翼型前缘冰层应力计算模型。采用有限元方法和平面三角形单元对控制方程组进行了求解,获得了外部空气动力和蒙皮表面加热对冰层黏附界面应力的影响规律。研究表明:蒙皮表面不加热时,来流速度影响了黏附界面应力的强度,来流攻角影响了黏附界面应力的分布,冰-蒙皮间黏附界面切应力最大值随来流速度呈近似线性增大趋势,但外部空气动力很难造成冰层破坏。蒙皮表面加热时,冰-蒙皮间黏附界面的耦合应力和冰层内部的主应力随着热流密度的增大而增大,很容易超过剪切强度,这是造成冰破坏的关键因素。耦合冰-蒙皮剪切强度随界面温度的变化关系,初步建立了基于应力分析和热/力耦合作用的冰破坏判断准则。外部空气动力产生的界面应力和蒙皮表面加热产生的界面热应力之和,必须大于与蒙皮表面温度相关的剪切强度,则冰层发生破坏,破坏位置是耦合应力超过剪切强度的区域。      

射频介质阻挡放电改善NACA 0015翼型气动性能的实验

介质阻挡放电(DBD)均匀稳定、易于敷设,是机翼/翼型等离子体流动控制(PFC)中最常用的激励方式。射频介质阻挡放电激励频率高、放电功率大,且能在流场中产生明显的加热,应用潜力大。采用射频电源驱动DBD激励器产生等离子体,分析放电的体积力、热特性和诱导流场特性,开展了射频介质阻挡放电改善NACA 0015翼型气动性能的实验,研究了占空比、调制频率、载波频率和电源功率等参数对流动控制效果的影响规律。结果表明:射频等离子体激励的体积力效应随激励电压的增大而增加;射频等离子体激励产生的热量在诱导的流场中进行传导,加速流场;当来流速度为20m/s,Re=3.36×10~5时,在翼型前缘施加激励,使翼型临界失速迎角推迟1°,最大升力系数增大6.43%,且在过失速迎角下仍具有流动控制效果,使升力下降变缓;调制频率越大,控制效果越好;存在最佳占空比、载波频率和功率,占空比对流场控制效果的影响最显著,最佳占空比、载波频率和功率分别为20%,460kHz和50W。射频等离子体激励以体积力效应、热效应和诱导壁面射流改善失速流场,使得NACA0015翼型气动性能极大改善,流动分离得到有效控制。     

三维叶轮机叶片黏性反问题设计改进方法

提出了一种改进的适用于三维黏性流场的叶轮机叶片反问题设计方法。该方法假设叶片的中弧线具有虚拟移动速度,其位移量由目标载荷与实际载荷的差值计算,并利用黏性底层厚度进行限制。采用三次B样条曲线插值方法对叶片中弧线进行光顺,新叶型通过更新后的中弧线和给定的叶片厚度得到。算例验算结果表明:该方法能够根据设计意图对叶型进行修改,鲁棒性强,收敛速度快,叶片的可变自由度高,不依赖于特定的网格和求解器,具备一定的通用性。      

基于深度学习的混合翼型前缘压力分布预测

提出了一种基于深度学习的混合翼型前缘压力分布预测方法,通过对翼型几何特征提取、压力分布曲线的参数化,建立了卷积神经网络模型（CNN）,并利用计算流体力学（CFD）的计算结果作为其训练样本,实现对混合翼型前缘压力分布的预测。结果表明:两种方法计算结果的拟合优度大于0.98,基于深度学习的计算方法耗时1.7 s,CFD方法耗时大于50 s,计算时间大大缩短。该方法能够在满足计算精度的条件下提高计算效率并可应用于其他的翼型设计过程。      

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构,层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟,研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型,基于结构化拼接网格,对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性,雷诺数对翼型气动特性具有显著影响。随雷诺数降低,时均分离泡外形增大、位置后移,平均阻力系数增大,特别是在较低雷诺数下,翼型升阻力系数随时间出现振荡现象。进一步研究表明,造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低,流动黏性增大,导致分离剪切层速度梯度减小,流动发生转捩及再附位置后移,直至翼型表面不再发生转捩和再附。     

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

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

基于任意多项式中弧线的单级高效率风扇设计

为了在压气机通流设计阶段考虑叶片弯掠效应,开发了基于流线曲率法的通流设计程序,提出一种基于四次多项式的任意中弧线叶片造型方法,并推导了任意回转面上的中弧线表达式。以此方法为基础,采用通流设计与叶片造型相互迭代的方式开展大流量跨声速风扇设计研究。此风扇级的设计点为巡航状态,设计流量为155kg/s、压比为1.54。研究结果表明:在设计状态,此风扇级的总压比为1.545,转子和级效率分别为0.939、0.916;在设计转速下,失速裕度为17%,转子和级最高效率分别为0.945、0.923;在起飞状态,流量接近440kg/s,效率与巡航状态相当,压比高于巡航状态。     

基于冷态叶型生成叶片加工坐标的方法

建立了基于冷态叶型生成叶片加工坐标的方法,包括排序、插值加密、叠加位移量、去倒圆及端区叶型拟合等步骤.以某多级轴流压气机为例,应用该方法生成的加工坐标叶片除端区外主流通道叶高内误差在0.01mm以下,叶片表面波纹度和前后缘曲率良好.通过进口级的数值模拟表明:该方法生成的加工坐标叶片与冷态叶型的误差对气动性能影响小,流量-总压比-等熵效率特性变化微小,叶片进出口气流速度、方向及流场细节等变化较小,满足工程应用的需求.      

高超声速飞行器宽速域翼型/机翼设计与分析

对于现代高超声速飞行器的设计而言,除了需要保证高超声速的性能外,还必须兼顾满足工程需求的亚跨超声速特性。首先,采用雷诺平均Navier-Stokes（RANS）方程流动求解器,结合基于Kriging模型的代理优化算法,开展了高超声速飞行器宽速域翼型的优化设计研究,设计出了一种下表面具有双"S"形特征的新翼型。综合性能评估结果表明,该翼型相比于常规的高超声速翼型,在跨声速和高超声速下具有更加优良的气动特性;其跨声速状态下的升阻比达到78.9,高超声速状态下的升阻比达到5.94,能够实现宽速域内良好的综合气动性能。其次,开展了仿德国"桑格尔号"（SANGER）空天飞机运载机机翼的气动特性研究,对配置宽速域翼型与常规高超声速翼型的机翼进行了气动力特性综合对比分析。结果表明,配置新翼型的机翼在宽速域范围内整体气动性能更优,说明所设计的宽速域翼型在三维机翼上也具有一定的实用价值。