大型 CFD 软件自动化测试平台的初步设计与实现

针对大型通用 CFD 软件研制的需求，为了提升 CFD 软件的开发效率，初步设计与开发了一个 CFD 软件自动化测试平台。该平台基于 MySQL 数据库，通过前台界面和远程集群后台管理程序的网络信息交互，实现了自主研发的通用 CFD 软件平台(HyperFLOW)的自动测试及验证与确认。根据需要适当修改前置处理接口，该平台亦可推广应用于其他 CFD 软件的自动测试。本文重点介绍了平台的基本框架和总体设计思路、平台的数据库和前台管理 GUI 设计、网络交互和后台管理等，并利用典型算例对计算结果的验证与确认过程进行了简要介绍，最后对平台构建的后续工作进行了展望。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

基于“各向异性”四面体网格聚合的复杂外形混合网格生成方法

复杂外形粘性流动的计算网格生成一直是阻碍CFD实际工程应用的瓶颈问题之一,虽然三棱柱/四面体混合网格的出现极大地缓解了这一问题,但是在飞行器几何外形非常复杂的情况下,在粘性边界层内生成高质量的贴体三棱柱网格仍然是非常困难的事情。本文从"各向异性"四面体网格的生成算法中得到启示,发展了一种基于"各向异性"四面体网格聚合的混合网格生成方法。通过对DLR-F6翼身组合体和F16战斗机外形的混合网格生成及数值模拟表明,该方法生成的混合网格具有网格质量好、自动化程度高等特点,对提高复杂外形飞行器关键气动参数的模拟精度具有重要意义。     

前体涡非对称分离机理及前缘吹气控制研究

通过设计对称性算法,求解层流Navier-Stokes方程,数值模拟了细长体在低超声速情况下前体背风涡随攻角演化的规律,在此基础上,进一步研究了前缘吹气对前体涡演化和侧向力特性的控制机理.根据数值模拟结果及分析,倾向于支持在层流框架内,前体涡的非对称失稳是一种对流不稳定机制,要想根据需要产生对称或不对称的前体涡,就必须外加持续的扰动.在约16°～48°攻角区间内,前缘吹气可产生规律性较好的侧向力,有可能直接利用前体涡进行横侧向控制.为工程实用化,需提高前缘吹气的激励收益效费比.     

任意平面域的三角形网格和混合网格生成

本文以阵面推进法为基础,给出了一种可对任意平面域进行三角形网格和混合网格自动生成的方法。该方法包含了易于生成与调控的矩形背景网格技术,简便的外形输入和初始阵面划分技术,加速查寻的数据结构,以及节点松弛和基于Ｄｅｌａｕｎａｙ准则的对角线交换技术。网格生成的实例表明,方法能快速可靠地生成任意平面域的高质量三角形网格和适于粘性计算的三角形／四边形混合网格,具有较好的通用性与实用性。     

基于HyperFLOW平台的客机标模CHN-T1气动性能预测及可信度研究

近年来,针对实际飞行器外形的CFD气动性能预测及可信度研究逐步得到重视,国内也召开了第一届航空CFD可信度研讨会(AeCW-1)。本文首先基于自主研发的CFD软件平台HyperFLOW对NACA0012翼型低速绕流进行了网格收敛性研究,验证了软件对简单湍流问题的模拟能力且具备良好的网格收敛性。其次,针对AeCW-1提供的客机标模CHN-T1,选用其中的两个算例:(1)定升力系数的网格收敛性研究;(2)考虑模型支撑和模型静气动弹性变形的抖振特性研究,研究了计算结果的网格收敛性及模型支撑、静气动弹性变形和湍流模型等对气动特性预测精度的影响。结果表明:观测精度阶和网格收敛性指数显示数值结果具有良好的网格收敛性和可信度;是否考虑模型支撑对力矩的预测精度影响较大,引入尾撑和弹性变形后,数值结果与实验结果吻合较好;对于CHN-T1标模,采用QCR关系式对原始SA模型进行修正对标模力矩特性有一定影响。     

非结构CFD软件MPI+OpenMP混合并行及超大规模非定常并行计算的应用

常规工程应用中，非定常数值模拟（如多体分离）的计算量十分巨大，如果为了达到更高的计算精度，加密网格或者采用高精度方法将会使得计算量进一步增大，导致非定常数值模拟在CFD工程应用中成为十分耗时和昂贵的工作，因此，提高非定常数值模拟的可扩展性和计算效率十分必要。为充分发挥既有分布内存又有共享内存的多核处理器的性能和效率优势，对作者团队开发的非结构网格二阶精度有限体积CFD软件（HyperFLOW）进行了混合并行改造，在计算节点间采用MPI消息传递机制，在节点内采用OpenMP共享内存的MPI+OpenMP混合并行策略。首先分别实现了两种粒度（粗粒度和细粒度）的混合并行，并基于国产in-house集群采用CRM标模（约4 000万网格单元）定常湍流算例对两种混合并行模式进行了测试和比较。结果表明，粗粒度在进程数和分区数较少的小规模并行时具有效率优势，16线程时效率较高；而细粒度混合并行在大规模并行计算时具有优势，8线程时效率较高。其次，验证了混合并行在非定常计算情况下的可扩展性，采用机翼外挂物投放标模算例，分别生成3.6亿和28.8亿非结构重叠网格，采用对等的（P2P）网格读入模式和优化的重叠网格隐式装配策略，网格读入和重叠网格装配耗时仅需数十秒；采用3.6亿网格，完成了非定常状态效率测试及非定常分离过程的湍流流场计算，在in-house集群上12 288核并行效率达到90%（以768核为基准），在天河2号上12 288核并行效率达到70%（以384核为基准），数值模拟结果与试验结果符合良好。最后，在in-house集群上采用28.8亿非结构重叠网格进行了4.9万核的并行效率测试，结果显示，4.9万核并行效率达到55.3%（以4 096核为基准）。     

中国CFD史

本文简要回顾了中国的CFD史,早在20世纪70年代,在钱学森教授的建议下,一些中国学者对计算流体动力学(CFD)这一新型学科产生了浓厚兴趣并开始研究.从那时起,中国的CFD逐渐发展繁荣起来,特别是从世纪之交,中国CFD非常活跃,为航空航天和其他民用领域做出了巨大贡献.在中国,特别注重CFD的研究与物理分析的结合.本文首先介绍了M5A的概念,它概括了中国CFD的主要研究领域;随后,介绍了设计格式的原则,以此原则为指导发展了一系列格式,比如NND格式、ENN格式和紧致格式.最近,发展了一系列高阶精度数值格式并应用于实际外形的复杂流动计算.再后,讨论了求解气动方程的网格尺寸准则.最后,同顾了如流动分离、流动拓扑结构和结构稳定性理论、旋涡沿其轴向演化和Hopf分叉理论、动导数和飞行稳定性等代表性工作.     

矩形／非结构混合网格技术及在二维／三维复杂无粘流场数值模拟 …

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术，其中非结构网格仅限于物面附近，发挥非结构网的几何灵活性，以少量的网格模拟复杂外形，同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。     

基于动态混合网格的鱼体巡游数值模拟

对鱼类等柔性变形体的数值模拟需要解决如下几个问题(1) 柔性体动态网格生成;(2) 基于动态网格的不可压流非定常算法.本文基于Delaunay背景网格插值方法,发展了一套快速生成柔性体变形的动态混合网格生成方法;对于不可压流问题,采用虚拟压缩方法进行求解,结合双时间步方法来处理非定常问题;最后数值模拟了鱼体巡游的非定常流动,分析了鱼体的摆动频率对力、功耗等流体动力学特性的影响.