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

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

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

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

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

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

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

复杂外形粘性流动的计算网格生成一直是阻碍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核为基准）。     

变形飞机动态气动特性数值模拟研究

根据近年来国内外提出的折叠翼变形飞机方案,设计了一种折叠翼变形飞机的3D数模,利用作者发展的动态混合网格技术和非定常计算方法,对该折叠翼变形飞机机翼折叠和展开过程中的动态气动特性进行了数值模拟,并与准定常计算结果进行了对比分析.计算结果表明:机翼在展开/折叠的过程中,飞机的气动特性会发生剧烈变化,在变形飞行控制率设计时需要引起高度的重视.     

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

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

HyperFLOW 软件数值模拟TrapWing 高升力外形北大核心CSCD

采用High-Lift研讨会提供的梯形翼外形(TrapWing),利用自主研发的基于结构、非结构网格的通用CFD软件HyperFLOW进行了数值计算,以评估其对复杂外形低速流的模拟能力。分别采用了三套不同拓扑结构的计算网格,包括两套非结构/混合网格和一套多块结构网格,每套网格又分为粗、中、细三种密度不同的网格数量以考察其网格收敛性。利用Richardson插值法,对计算结果开展了可信度分析。结果表明,不管是结构网格还是非结构网格,HyperFLOW均建立了可接受的可信度;对于高升力外形数值模拟,SA湍流模型要比SST湍流模型模拟的更准确;在失速迎角附近,现有的二阶精度解算器仍需持续改进。