TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟

 采用“亚跨超声速计算流体力学软件平台”（TRIP2.0）数值模拟了阻力预测小组（AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ）翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果。重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致。采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响。     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

DPWⅢ机翼和翼身组合体构型数值模拟

采用亚跨超CFD软件平台(TRIP)数值模拟了DPW Ⅲ提供的DPW_W1/W2两种机翼构型和DPW-F6/F6_FX2B两种翼身组合体构型,主要目的是通过两种机翼构型和两种翼身组合体构型的数值模拟,研究网格密度对运输机构型气动特性计算结果的影响。数值模拟采用的多块对接网格来自AIAA CFD Drag Prediction Workshop Ⅲ(DPW Ⅲ),采用National Transonic Facility(NTF)的试验结果和CFL3D的计算结果作对比。详细研究了网格密度对两种机翼构型和翼身组合体的总体气动特性和压力分布的影响。采用SST两方程模型计算两种构型均得到了网格收敛结果,网格密度主要影响压差阻力而对摩擦阻力影响较小,计算结果较好地预测了机翼和翼身组合体外形优化前后总体气动特性的变化量。     

NASA桘CRM阻力预测模型的计算验证

采用自研的非结构网格解算器UNSMB 进行了AIAA第四届阻力会议提供的NASA-CRM 翼身组合体（ WB）以及翼身组合加平尾（ WBT）两种构型的计算验证。重点分析了WBT模型的网格收敛特性、升阻力曲线、雷诺数效应、下洗效应以及压力分布等,并把计算结果与阻力预测会议上ONERA的计算结果进行了对比。分析结果显示,非结构混合网格解算器的计算结果与ONERA的计算结果吻合度较好,同时在一定程度上验证与确认了解算器的阻力预测精度。     

NASA共同研究模型翼身组合体数值模拟研究

采用结构网格求解器WiseManPlus软件,对DPW5提供的NASA共同研究模型CRM翼身组合体构型进行数值模拟研究.采用雷诺平均N-S方程,选择SA一方程湍流模型和SST两方程湍流模型,开展了网格收敛性分析及抖振分析计算研究,评估了该软件对此类民机构型的阻力预测能力.按照DPW5的要求,采用会议统一提供的基准结构网格进行计算分析,研究结果表明,计算的巡航设计点的阻力及大攻角下的气动力系数与试验值吻合较好,达到会议统计分析平均水平.计算结果表明,所采用的结构网格求解器WiseManPlus软件计算精度较高,适用于大型民机高速巡航及抖振现象的数值模拟研究.     

多段翼型粘性绕流计算研究

针对多段翼型外形结构复杂,结构网格生成难度大的缺点,运用割补法对重叠网格技术进行了研究,生成贴体及与边界正交的高质量的结构化网格,在识别贡献单元方面,介绍了一种基于kd树的找重方法.流场求解应用雷诺平均Navier-Stokes方程结合Spalart-Allmaras一方程湍流模型在所生成的重叠网格上对多段翼型的粘性绕流进行了模拟,为了兼顾流场计算效率问题,在重叠网格的各个子块上采用多重网格方法.计算结果表明,用该方法可以很好地预计压力分布,并且计算效率得到很大的提高.     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

机体/动力装置一体化分析中的动力影响效应数值模拟

采用合适的进、排气边界条件, 对带动力的单独涡扇发动机模型和带动力的DLR-F6翼身架舱复杂组合体(WBPN)模型的流场进行了数值模拟和分析.运用Roe三阶迎风偏置通量差分裂方法和隐式近似因子分解方法求解雷诺平均Navier-Stokes方程, 基于点搭接多块网格技术生成高质量计算网格, 通过多重网格技术来加速收敛.着重对动力发动机进、排气边界条件的给定进行了推导, 以准确地模拟发动机动力效应.计算结果与实验值吻合良好, 验证了所采用计算方法、网格分块策略和边界条件处理方法的正确性.      

旋翼流场计算嵌套网格并行装配方法改进研究

针对嵌套网格系统下进行旋翼流场计算时的网格装配问题,给出了一种改进的并行化网格装配模型。在该模型中,使用压缩存储的辅助结构网格进行并行"洞切割",生成插值边界;应用基于插值点坐标范围的局部可自适应辅助结构网格定位,缩小了处理器之间进行并行贡献单元搜索时的交互数据范围及生成辅助结构网格的内存消耗,结合局部范围的Neighbor to Neighbor(N2N)搜索重启和穷举方法,消除了网格分区、物面、外场等边界面对贡献单元搜索过程的影响,且方法对网格的拓扑没有限制,适合于结构、非结构、直角网格等网格块的装配过程。在流场计算方面,采用并行化的隐式双时间格式数值求解Navier-Stokes方程,用于旋翼数值模拟问题。理论分析和计算测试表明,在使用粗化的辅助网格时,改进的计算方法能够达到传统方法贡献单元搜索速度近似10倍,在低内存开销情况下降低网格装配过程对流场计算时间占比。在内迭代大于7次时,网格装配时间小于流场计算的1%。该方法是对旋翼多周期计算时反复网格装配过程的很好的效率改进。     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

Enhancement on parallel unstructured overset grid method for complex aerospace engineering applications

In the present study, an efficient overset grid method by means of parallel implicit hole-cutting is proposed for the sake of simulating unsteady flows in aerospace engineering involving multiple bodies in relative movement. In view of the degraded computational efficiency and robustness for conventional overset grid assembly, several innovative techniques are developed within the overset grid assembly process, viz., a bookkeeping alternative digital tree method to speed up the donor-cell search...     

运输类飞机巡航阻力CFD计算分析

针对运输类飞机设计过程中的巡航阻力准确预测问题,采用基于多块结构网格求解三维雷诺平均Navier-Stokes方程的并行流动求解器“CCFD-MB”,对AIAA阻力预测小组提供的典型运输类飞机DLR—F6翼身构型进行了CFD计算分析,比较了小同规模网格及不同湍流模型对计算结果的影响;并通过与DLR-F4翼身构型计算结果...     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切.使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格.应用特征恢复技术恢复凹面特征并改善凸面处的网格质量.流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量.加入Spalart-Allmaras一方程湍流模型以获得湍流解.上下对称高斯-赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算.     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

非结构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核为基准）。     

CRM翼/身/平尾组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼/身/平尾组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测方面的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细3套网格,对应的y⁺分别为1、2/3和4/9。通过对CRM模型的计算和分析,研究了网格规模对气动特性、压力分布和翼根后缘局部分离区的影响。通过与DPW4统计结果和部分实验结果的对比,高阶精度数值模拟结果表明,本文的气动特性计算结果与统计平均结果吻合较好;网格密度对机翼内侧的展向站位压力分布影响较小,对机翼外侧展向压力分布影响较大;网格密度对翼根后缘局部分离区略有影响。