扩散作用对y-Re优模型转捩预测的影响

Gamma．theta转捩模型是Menter及其合作者提出的一个基于经验关系的转捩模型,其经验关系式和模型参数是通过来流米雷诺数均在150万以下的仍系列平板实验在CFX程序中标定的。首先通过33A平板修正了转捩经验关系式．然后对S＆K平板实验进行了数值模拟,研究了较高雷诺数下模型控制方程扩散项系数以及与数值扩散相关的流向网格尺度对gamma．theta模型转捩预测结果的影响。模拟的结果表明,这些因素对转捩位置都有明显影响,并且gamma方程的扩散项系数还影响转捩完成后全湍区的摩阻值。三种因素对转捩预测影响的趋势分别为：gamma方程的扩散项系数越大,theta方程的扩散项系数越小,流向网格尺度越大,则转捩位置越靠前。     

超临界翼型稳健型优化设计研究

基于标准遗传算法、RBF神经网络以及类函数/型函数变形技术建立了翼型气动优化设计系统,在选取3次CST方法即分别对翼型上下表面采用4个设计变量进行参数化,对某型客机基本翼型在给定设计指标下进行优化设计,分别研究了巡航状态下的气动优化,以及结合蒙特卡洛分析方法马赫数随机平均分布下的翼型稳健型优化设计。结果显示,优化后的翼型的气动特性有着显著提高。     

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

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

基于深度学习驱动的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型高导热碳纤维丝可以将高温区域的热量沿纤维方向定向疏导至低温区域。     

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

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

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

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

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

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

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

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

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

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

后缘连续变弯度对跨声速翼型气动特性的影响

针对后缘连续变弯度对跨声速翼型气动特性的影响进行了研究。首先不考虑翼型后缘连续变弯度,基于搭建的优化设计系统对跨声速翼型进行气动减阻优化设计,通过添加不同的约束优化得到两种跨声速翼型:无激波翼型和超临界翼型。然后在这两种翼型的基础上,以后缘偏转角度为设计变量、以阻力系数最小为目标,针对不同的升力系数分别进行优化设计,并根据优化结果深入分析后缘连续变弯度对这两种翼型极曲线特性的影响机理。优化结果表明:无激波翼型与超临界翼型相比,其设计点处的气动特性较好,但鲁棒性较差;升力系数小于设计升力系数时,应用后缘连续变弯度后,无激波翼型的极曲线特性明显提高,减阻最高达到3.9%,而超临界翼型的极曲线特性提高不明显;升力系数大于设计升力系数时,应用后缘连续变弯度后,无激波翼型和超临界翼型的极曲线特性都明显提高,减阻分别达到2.4%~18.1%和1.7%~13.2%。