CFD软件自动化验证确认云平台设计与实现

针对大型计算流体力学(CFD)软件的验证与确认,为了减少人工成本,提高软件质量和开发效率,并适应于未来高性能计算发展的要求,基于云计算思想,提出了自动化测试云平台的解决方案。该方案采用模块化的浏览器/服务器(B/S)网络架构;利用LAMP(Linux+Apache+MySQL+PHP/Python)开发工具;建立了持续集成的专业数据库;构建了涵盖可靠集群监控、复杂作业调度及大规模并行计算功能的云端环境;实现了通过便捷的网络访问,自动加载测试算例、提交集群计算、监测实时进度、自动化后处理、输出分析结果,并提供丰富的实验对比图、误差分析报表和汇总报告等;完成了对大型CFD软件的自动化验证与确认。将该方案应用于某大规模并行计算的CFD软件,验证了该解决方案的可行性与实用性。     

DES与DDES在湍流分离中的原理与性能研究

随着工程上流动结构的日益复杂，兼具雷诺平均Navier-Stokes(RANS)方法高效率与大涡模拟（LES）高精度的分离涡模拟（DES）类混合方法成为现阶段工程中最有效的湍流模拟方法之一。围绕DES类混合方法中的DES与延迟分离涡模拟（DDES）方法开展工作，分析二者开关函数构造上的不同，研究延迟因子作用机理，并考察DES与DDES方法的求解能力。研究表明：DES与DDES方法在模拟表现上存在一定差异，DDES方法通过引入延迟因子，保护RANS求解区域，改善模化应力不足，降低了DDES方法对交界面系数CDES敏感程度；DDES方法在计算过程中容易出现过度保护，导致求解瞬时涡结构能力不如DES方法，分析与延迟因子引入比重及开关函数构造形式有关。     

基于图剖分的多块结构网格负载平衡方法

负载平衡是影响并行计算性能的重要因素。针对多块结构网格,给出了一种改进的多层次图剖分负载平衡方法。该方法设计了新的网格剖分算法,采用改进的子块分裂方法与图剖分算法的循环调用实现结构对接网格剖分,并通过建立不同物体重叠网格间的连接关系,实现了结构重叠网格的负载平衡。采用2个典型算例对方法进行了对比验证,数值结果表明,子块分裂方法对剖分结果具有重要影响,采用循环调用算法及改进的子块分裂方法能有效地实现计算负载均衡及通信量优化,同时显著减少了网格块数及因虚网格导致的内存需求,有利于提高并行效率。该负载平衡方法与网格拓扑无关,适用于多块结构对接网格及重叠网格,且整体型剖分方式对于多块结构重叠网格具有更好的剖分效果。     

带反馈控制的基本带翼外形的耦合气动/运动数值模拟研究

介绍了一种先进的、时间精度的多学科数值模拟平台的开发和应用,该平台可用于预测飞行器机动过程中的非定常气动特性和飞行力学特性.通过耦合求解非定常Navier-Stokes方程、刚体动力学方程及控制律,对基本带翼外形进行了虚拟飞行数值模拟研究.求解器采用基于结构网格的有限体积法,时间方法采用双时间步方法.使用反馈控制器获得舵偏控制指令,重叠网格方法用来模拟物体间的相对运动.模拟结果显示控制器的增益必须正确选择才能获得正确的响应.仿真结果表明本文发展的虚拟飞行数值模拟方法能够预测包含动边界的复杂非定常流动,具有较强的工程适用性.     

Roe格式中不同类型熵修正性能分析

针对Roe格式中较为流行的3类熵修正:Muller型、Harten-Yee型和Harten-Hyman型熵修正,从理论分析入手, 辅以Euler方程为控制方程的激波管问题,前台阶流动和运动激波的双马赫反射3个数值实验,对不同熵修正的性能做了深入研究,得出如下结论: Muller型和Harten-Yee型熵修正方法作用于激波和非物理的膨胀激波2种情况:激波情况下引入的数值耗散,能根本改善"Carbuncle"现象,膨胀激波时其数值粘性也足以使膨胀激波得以耗散;而类Harten-Hyman型熵修正只对膨胀过程起修正作用,对激波情况无效,不能改善"Carbuncle"现象,该类熵修正的数值粘性较小,不足以使膨胀过程正确求解;直接采用 δ 值代替特征值的特征值修正方法效果好于传统的特征值修正方法.     

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

吸气式高超声速飞行器多参数灵敏度分析

为实现吸气式高超声速飞行器多参数情况下的灵敏度分析及参数分类,降低设计的复杂度,在吸气式高超声速飞行器参数化建模的基础上,首先采用正交试验设计生成样本,再通过计算流体力学(CFD)进行高精度气动力性能计算,最后运用方差分析法进行气动性能的参数灵敏度分析,在运用较少的试验样本点的情况下,即可完成多参数、复杂构型条件下气动性能的参数灵敏度分析。结果表明,该方法可以正确地分析出参数对飞行器气动性能的敏感程度,得到参数对气动性能的影响规律。同时,通过灵敏度分析的计算样本,还可以初步选出气动性能较优的飞行器构型,为地面试验和优化设计提供参考。     

射流对高超声速进气道起动性能的影响

为探索利用射流技术降低进气道起动马赫数的可行性,对二元高超声速进气道二维流场进行了数值模拟,通过对比不同工况的流场结构、流量系数及总压恢复系数,分析了射流对高超声速进气道的作用效果,并研究了射流速度、压强及倾角对进气道起动性能的影响。分析结果表明:施加射流,激波与进气道边界层原有干扰形式发生改变,是降低进气道起动马赫数的主要原因。研究还表明,增大射流速度利于提高控制效果,但持续增大射流速度,会造成隔离段反压升高,并且这一现象与射流压强相关,降低射流压强能使进气道起动的射流速度区间扩大,同时在不同射流倾角下,上述规律表现一致。该研究揭示了进气道起动能力随射流参数变化的系统性规律,可用于指导工程设计及优化。     

适用于混合网格的并行GMRES+LU-SGS方法

给出了一种适用于混合网格的并行无矩阵GMRES+LU-SGS隐式时间格式。首先采用LU-SGS方法迭代若干步以获得一个合适的初场,然后切换到GMRES方法在每一时间步内近似求解,并将LU-SGS方法作为其预处理器。为加速收敛,将CFL数随着残差的降低逐步放大;为减少存储量和计算量,通量Jacobian采用无矩阵处理。在保证与串行执行一致的前提下,采用基于共享内存的OpenMP方法实现了并行计算,并通过对网格的分组避免了内存争夺。算例验证表明,方法极大地提高了计算收敛效率,并行结果与串行结果完全一致,计算结果与实验结果吻合较好。     

翼型的气动最优化设计方法和反设计方法

本文将数值优化方法同计算流体动力学(CFD)相结合,形成两种跨声速翼型的气动设计方法,即:最优化方法和反设计方法.采用求解Euler方程的有限体积法计算流场,通过结合优化算法,从正反两个方向对跨声速翼型进行气动优化设计.实例证明它们是翼型气动设计的有效方法,有较高的工程应用价值.