基于数据挖掘的飞行器气动布局设计知识提取

为了更深入地理解飞行器气动布局设计优化中多目标/多设计变量间的影响关系,提高优化模型的科学性及优化效率,对基于数据挖掘技术的飞行器气动布局隐含设计知识提取问题开展了探索研究。以高升阻比滑翔飞行器布局设计优化问题为例,基于当前比较有代表性的方差分析、等度量映射、决策树、自组织映射4类机器学习算法对气动布局优化设计中产生的中间数据进行了挖掘分析。对不同方法得到的升阻比、横/侧向稳定性及容积率4种目标性能间的权衡关系,目标性能与设计变量间的敏感性关系及产生较优布局外形的设计变量取值规则进行了综合对比分析,凝练形成了适用于该类飞行器的设计知识,同时对4种方法的特点及适用性进行了总结分析,给出了相关结论。     

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

考虑数值离散误差的湍流模型选择引入的不确定度量化

在流体力学数值模拟过程中存在着多种来源的不确定因素,科学、定量地描述这些因素对模拟结果的影响对模型确认、工业产品设计优化和性能评估等过程十分重要。数值离散、模型选择和模型预测偏差是模拟过程中3种重要的不确定性来源,为将这3种不确定性因素对模拟目标量的影响统一考虑,发展了考虑数值离散误差的贝叶斯模型平均方法。首先,通过对数值离散解和网格尺度进行拟合完成数值离散误差估计,得到每个备选模型真实解的置信区间。然后,通过嵌套方法和条件优化算法,结合贝叶斯模型平均方法构建目标量的概率盒,定义目标变量累积分布函数的上、下限,以此分析其置信区间。最后,针对NACA0012低速绕流和CHN-T1跨声速绕流问题,给出了同时考虑上述3种不确定性因素之后升、阻力系数的置信区间分析示例。     

随机和认知不确定性框架下的CFD模型确认度量综述

随着计算机技术迅猛发展,计算流体力学(CFD)数值模拟技术已广泛应用于航空航天、船舶、风力、水利等领域,为装备的复杂流体分析、性能参数评估以及产品气动设计与优化提供重要验证与决策手段。现有的CFD技术都是在确定性的数学模型、物理参数、边界条件下进行数值模拟的。然而,由于物理过程的复杂性及人们的认知偏差,CFD技术中存在着模型参数、数值离散、模型形式等诸多不确定性因素且表征形式各异,给CFD数值模拟的可信度评估带来了极大挑战。文章阐述了CFD技术中面临的不确定性因素,探讨了主流模型确认度量方法的应用条件和适用范围。着重介绍在认知不确定性下的模型确认度量方法,包括几类重要的区间、概率盒下模型确认度量方法。通过NACA0012翼型绕流问题说明了CFD技术中考虑各种不确定性下的模型确认度量方法的有效性。     

烧蚀热响应计算中的不确定性传播分析方法研究

发展了基于代理模型的不确定度量化方法，以少量确定性模拟结果为样本，通过代理模型来构建目标变量和不确定性输入参数的响应关系，以此来分析目标变量的不确定性信息，包括均值、标准差、概率密度分布等统计信息以及各个输入参数对目标变量的敏感性大小。本文以烧蚀热响应计算为例，通过分析材料物性参数的不确定性对烧蚀热响应预测的影响，表明该方法具有准确度高、效率高的特点，为工程中的高维不确定度量化问题提供了很好的解决思路。     

基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。     

湍流模型系数不确定度对翼型绕流模拟的影响

使用非嵌入式多项式混沌方法研究了湍流模型系数的不确定度对RAE2822跨声速翼型绕流模拟的影响。计算中关注了数值模拟的积分量（升力系数、阻力系数）和局部量（壁面压力、摩擦系数和空间马赫数分布）的不确定度量化结果。首先,从单输入变量入手,研究卡门常数的不确定度对数值模拟的影响。然后,同时考虑Spalart-Allmaras模型中9个参数的不确定度带来的影响。通过多项式混沌展开,得到系统输出对不确定输入变量的响应,由此可以得到输出的统计特性,包括平均值、方差和极值等信息。最后,在多变量不确定度量化过程中,通过Sobol指标来量化每个输入变量的不确定度对输出不确定度的贡献程度。本文计算只考虑了RAE2822跨声速翼型模拟的单一计算状态,影响规律是否可以推及其他工况和算例需要进一步检验。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。