翼吊发动机短舱对增升装置气动特性影响研究

增升装置的气动效率对于全机整体性能有重要影响,而当代大多数跨音速运输类飞机均采用翼吊发动机布局。通过对有、无翼吊短舱的两个三维增升装置的数值模拟,研究了翼吊发动机短舱对于增升装置气动性能的影响。结果表明,翼吊发动机短舱会严重恶化增升装置的气动性能。流动机理分析表明,短舱挂架与机翼前缘结合处的缝翼缺口及大迎角时绕过短舱的分离气流是造成增升装置气动性能恶化的两个主要原因。     

基于深度学习驱动的L型定向热疏导机理

碳/碳（C/C）复合材料具有热导率大、比强度高、耐烧蚀和耐冲刷等优异特性,被广泛应用于飞行器的热防护系统中,其有效导热系数对于实际应用而言是重要的热物理性质,尽管可以通过有效介质理论、对热扩散方程直接求解和玻尔兹曼输运方程等传统方法计算C/C复合材料有效导热系数,但这些数值方法通常十分耗时。本文引入深度学习方法,将格子玻尔兹曼（LBM）的三维格子模型作为三维卷积神经网络（3D-CNN）微观结构,不仅解决了三维微观结构模型难以捕获的问题,还便于实现数值计算模型和CNN模型的同步简化,利用3D-CNN快速精准地预测三维三相C/C复合结构的有效导热系数,基于此对内置L型高导热碳纤维丝的定向热疏C/C复合结构的有效导热系数进行快速预测和研究。研究表明,CNN模型在LBM传热计算上表现出强大的学习能力,但在测试样本结构孔隙率过分超出训练集时预测误差将大幅增加,且当孔隙率变化范围从30%~35%变化到55%~60%时,CNN模型"内插"预测的相对误差较模型"外推"降低了0.93%~30.72%。在C/C复合结构中内置L型高导热碳纤维丝可以将高温区域的热量沿纤维方向定向疏导至低温区域。     

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。     

机翼后缘连续变弯度对客机气动特性影响

后缘连续变弯度机翼在提高民用客机气动特性方面有较大的潜力,近年来被广泛关注。基于建立的全局优化设计系统,研究了机翼后缘连续变弯度对宽体客机翼身组合体气动特性的影响。首先,采用自由型面变形(FFD)技术建立了后缘连续变弯度的参数化方法。然后,采用RANS方程作为流场评估方法,针对翼身组合体构型设计点附近升力系数开展了机翼后缘连续变弯度气动减阻优化设计。最后,探索了仅外翼段后缘连续变弯度和内外翼后缘均连续变弯度优化设计结果的异同。优化结果表明,升力系数小于设计升力系数时,在只考虑外翼段后缘连续变弯度的设计中,不易实现激波阻力和诱导阻力同时降低,考虑内翼段后缘连续变弯度后,减阻量较前者更为明显;升力系数大于设计升力系数时,外翼段和内外翼的后缘偏转均可实现诱导阻力和激波阻力的同时降低,且减阻量相差不大。     

一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。      

适用于非周期流固耦合问题的时间谱方法

针对非定常流固耦合设计、分析问题所面临的计算效率与精度、鲁棒性之间的矛盾,采取将Chebyshev谱方法与基于非定常面元和几何非线性梁有限元模型的流固耦合分析方法相结合的方式,建立了针对大展弦比机翼非定常流固耦合优化设计问题的,可与伴随方法相结合的时间谱方法。Chebyshev谱方法直接对流固耦合的控制方程进行处理,利用Chebyshev算子替换系统状态变量,将非定常问题转化为Chebyshev控制点处耦合的定常问题。通过这种方式建立的流固耦合时间谱方法具有较高的计算精度、效率和足够的鲁棒性。验证算例及Goland机翼颤振速度计算实例表明,Chebyshev谱方法的计算精度随着Chebyshev控制点个数的增加而不断增大。只需选取较少的控制点,Chebyshev谱方法便可以达到满足精度要求的计算结果。与此同时,建立的流固耦合时间谱方法不仅适用于周期性非定常问题还适用于非周期性非定常问题。     

复杂几何细节对增升装置气动性能影响研究

 采用数值模拟的方法研究了主翼翼根几何形状、翼吊发动机短舱、缝翼滑轨及襟翼滑轨舱等几何细节对增升装置气动性能的影响。研究结果表明:切割前缘缝翼时,将大部分翼根整流包留在主翼上会在大迎角下产生低能量的分离涡,造成增升装置气动性能显著恶化,而将大部分翼根整流包切割到前缘缝翼上,能破坏低能量分离涡的产生;大迎角下,短舱上表面、挂架表面及缝翼与挂架之间的间隙产生的分离气流会直接流到主翼上表面,形成大范围的死水区,因此,大尺寸的翼吊发动机短舱会造成增升装置失速迎角及最大升力系数的大幅减小,但安装在短舱适当位置、适当形状的涡流片产生的强漩涡能消除大部分的死水区,挽回部分气动性能损失;缝翼滑轨产生的低能量尾迹会混入主翼附面层,使其能量降低造成升力系数减小,极端情况下缝翼滑轨会直接诱发大范围的流动分离,造成增升装置气动性能的显著恶化;襟翼滑轨舱因其较大的几何尺寸会减小襟翼缝道的面积使得襟翼缝道射流加速,有利于吹走襟翼表面的物面分离。     

绕任意刚性翼型,机翼非定常流动的跨声速欧拉方程解

本文在有限体积方法的基础上引入变系数隐式残值光顺及固接网格法求解了绕二维翼型及三维机翼非定常流动的跨声速欧拉方程,计算结果与同类结果及实验结果吻合。     

一种适用于浸入有限元方法的网格自适应方法

针对动边界流固耦合的数值模拟问题，基于浸入有限元方法提出了一种耦合流场特征和几何特征的笛卡儿网格局部加密自适应方法，克服了单个自适应指示因子无法精确捕捉固体运动的特征的不足。在耦合自适应策略中，分别以流场涡量和固体位置作为流场和几何信息指示因子来驱动网格自适应。通过方腔顶盖驱动圆盘流动算例，以圆盘体积守恒和特征点的运动轨迹验证耦合自适应方法的优势。计算结果表明:仅基于流动特征的自适应不能很好地保证圆盘的体积守恒；仅基于几何特征的自适应无法有效追踪圆盘的轨迹；而耦合自适应策略能同时较好地保证两项指标的计算精度，在保证总体计算自由度不变的情况下，圆盘区域速度散度2-范数降低了一个数量级，圆盘的轨迹误差2-范数降低了2个数量级。      

一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。