基于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翼型气动性能极大改善,流动分离得到有效控制。     

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

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

基于代理模型的大型民机机翼气动优化设计

先进的气动优化设计思想与方法,对于提升大型民机气动与综合性能具有至关重要的意义。探讨了大型民机超临界机翼气动优化设计的基本准则和要点,并结合代理优化算法,提出了一种面向工程应用的多轮次高效全局气动优化设计方法。首先,通过一系列解析函数/翼型/机翼优化测试算例进行了验证。其次,将所提出的方法与人工修型相结合,开展了针对宽体客机超临界机翼的两轮气动优化设计,使其气动性能得到显著改善。最后,采用不同的雷诺平均Navier-Stokes（RANS）方程求解器对安装优化机翼的全机巡航构型进行了典型状态的气动性能综合评估。研究结果表明,所提出的代理优化算法具有很高的优化效率、较强的约束处理和全局优化能力;将所发展的基于代理模型的多轮次气动优化设计方法与人工修型相结合,能够获得满足设计要求的气动外形,验证了该方法在大型民机超临界机翼气动设计中的有效性和工程实用性。     

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

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构,层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（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以下,叶片表面波纹度和前后缘曲率良好.通过进口级的数值模拟表明:该方法生成的加工坐标叶片与冷态叶型的误差对气动性能影响小,流量-总压比-等熵效率特性变化微小,叶片进出口气流速度、方向及流场细节等变化较小,满足工程应用的需求.      

超临界压力下碳氢燃料裂解与流动传热模拟的快速算法

为了快速预测发动机冷却通道内碳氢燃料在考虑热裂解时的流动传热特性,基于正癸烷热裂解反应机理,建立了一套模拟正癸烷裂解吸热和超临界压力传热现象的快速算法。采用三维物性库查表算法计算裂解反应混合物的热物性,同时对组分输运方程进行简化,简化后只需求解1个组分输运方程。通过与现有的实验和数值结果进行比较,检验了快速算法的计算效率和可靠性。结果表明,快速算法与求解全组分输运方程的算法精度相当,但计算效率提升了约20倍。最后采用该算法对三维矩形截面冷却通道内的超临界压力正癸烷裂解与传热过程进行数值模拟,进一步考察了本文快速算法的计算精度及其工程应用价值。