二维多段翼型网格生成与数值模拟研究

为提高多段翼型的网格生成效率和数值模拟精度,发展了一套自适应混合笛卡尔网格(AHCG)生成方法和基于有限体积方法的雷诺数平均Navier-Stokes (RANS)的数值求解技术.混合笛卡尔网格由围绕物体几何外形的贴体结构网格和填充流场其他区域的笛卡尔网格构成,两套网格之间的信息传递由“贡献单元”提供,且“贡献单元”由基于ADT(Alternating digital tree)技术的搜寻方法获得.为更准确地捕捉流场信息,采用了基于流场特征的网格自适应技术.数值模拟结果显示,AHCG方法能够准确且高效地模拟高升力多段翼型绕流问题.     

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

国家数值风洞(NNW)工程中的黏性自适应笛卡尔网格方法研究进展

相比传统的结构和非结构网格,笛卡尔网格具有自动化和高质量兼顾的优势,是未来网格技术发展的重要方向。面向笛卡尔网格生成软件开发需求,依托国家数值风洞（NNW）工程,介绍了黏性自适应笛卡尔网格技术相关课题的研究工作进展,重点从网格生成技术、自适应方法、黏性物面边界处理方法等方面展开。网格生成技术方面,基于全线程树数据结构的思想,开展了叉树数据结构的优化;针对物面离散单元信息的快速检索,构建了稳定、平衡、高效的叉树式物面网格数据结构;对于网格单元类型的判断,发展了染色方法和包围盒方法以提高效率;针对大规模网格量,发展了高性能并行计算技术。自适应技术方面,发展了基于几何特征和流场解特征的自适应方法,建立了流场结构捕捉判据,并开展了三维复杂构型的应用;为了降低全局网格规模,发展了各向异性自适应方法,构建了各向异性特征识别的综合判据;发展了边界层法向射线加密技术,以减少边界层内的网格数目。物面边界处理方面,发展了重叠笛卡尔网格方法,构建了重叠网格尺度匹配技术和数据交互方法;开展了壁面函数方法在笛卡尔网格框架下的适应性应用,以放宽壁面网格尺度限制、降低边界层网格规模,并通过典型算例验证了方法的可靠性。     

基于全线程树数据结构的笛卡尔网格高效生成技术

笛卡尔网格方法具有易于自适应、自动化程度高、网格质量好等优势。由于其非贴体特性,在处理复杂构型或者复杂流动问题时候往往存在网格量过大的问题,不容忽视。笛卡尔网格生成的效率直接影响了整个计算周期的长短,有必要发展高效的生成技术。对于笛卡尔网格而言,决定网格生成效率的关键在于网格数据结构,其直接影响计算量和存储量。针对三维构型进行笛卡尔网格生成,发展了邻居查询更快捷、内存利用率更高的全线程树数据结构,并在本文的方法框架下进行了适应性应用和改进。同时,为了高效地判断网格单元类型,构建了物面单元的快速检索方式,并引入了染色方法,进一步提高网格生成效率。还提出了一种鲁棒的奇异性检测算法,保证网格单元类型判断的鲁棒性。在流场解自适应方面,采用的是速度散度和旋度相结合的三维判据,以保证对于多种流动特征的捕捉能力。通过圆球、导弹、翼身组合体、机翼-副翼等算例进行了考核验证,经对比,网格自适应位置与理论解吻合较好,且网格单元生成耗时短、平均耗时受物面网格分布影响小,证明了方法的可靠和高效性。     

一种适合于旋翼涡流场计算的非结构自适应嵌套网格方法

针对旋翼涡流场的特点,提出了一种新的非结构自适应运动嵌套网格方法。在该运动嵌套网格方法中,贴体网格采用非结构网格,背景网格使用可自适应的笛卡尔网格,以便发挥非结构网格的形状描述能力和笛卡尔网格的空间自适应优势,方便网格生成过程并提高计算精度,同时克服在物体表面附近自适应笛卡尔网格生成工作量大和过程烦琐的不足。在此基础上,使用求解N-S主控方程的方法,分别对旋翼的悬停和前飞流场进行了数值模拟,计算结果表明本文的非结构自适应运动嵌套网格方法对提高旋翼涡流场的细节捕捉能力和计算精度是很有效的。     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

基于粘性笛卡尔网格的N-S方程计算

本文采用基于叉树数据结构的粘性笛卡尔网格求解N-S方程。为了使网格的分布更为合理,采用基于物体几何外形的方法对网格进行自适应加密。在物面处采用光顺、投影、分割加密的方法保证了物面处网格的正交性与贴体性,从而生成高质量粘性笛卡尔网格。将有限体积法应用于粘性笛卡尔网格求解N-S方程,采用了SA紊流模型。算例计算结果与实验数据符合良好。     

基于参数坐标的曲面网格生成方法研究

体网格的品质很大程度上依赖于其边界由西阿格的品质．基于参数坐标的曲面网格生成方法是生成曲面网格的最常用的方法。本文分析比较了两种有代表性的基于参数坐标的曲面网格生成方法；讨论了两种方法的优缺点而适用范围；同时，介绍了在CAEGrid网格生成系统中采用的一种通用的曲面网格生成策略．     

一种新的边界处理方法在笛卡尔网格中的应用

采用自适应笛卡尔网格方法求解Euler方程,固壁边界条件通过一种离散强迫浸入边界方法引入,边界切割网格与流场网格保持一致,从而有效解决了小切割网格单元的时间步长限制问题.针对自适应笛卡尔网格的特点提出了新的虚拟点反射方法和插值方法,解决了数值吸力峰值问题.对RAE2822翼型、双NACA0012翼型和Suddhoo三元翼型的流动状况进行了无黏数值模拟,并与现有的理论解、单域以及分区结构网格解进行了对比.结果表明:该离散强迫浸入边界方法有效地引进了固壁边界条件,结合自适应笛卡尔网格技术能够准确模拟复杂几何外形下的流动状况.      

基于RANS方法的民用飞机网格分析

基于RANS方程分别对代表民机巡航的DLR—F6模型和增升的NASATRAPWING模型进行多块结构化网格生成技术研究,通过网格拓扑、分布、质量参数等多项研究,得到影响民机外流场粘性计算需要的网格生成因子,形成一套适用高精度阻力预测的民机CFD网格生成指南。同时,为验证网格设计方法的可靠性及准确性,分析了不同网格生成准则对于CFD计算模拟精准度的影响,为民机计算策略研究及气动力设计提供了参考依据,具有重要的工程实用价值。     

二维Euler方程的自适应非结构网格生成方法

发展了一种用于数值模拟二维Euler方程的自适应非结构网格的生成方法。这一方法根据背景网格提供的信息直接确定内部结点的位置，然后依照Delaunay准则，将新结点加入到已有网格划分中。用于平观流及管道流的实际计算中，结果表明，该方法不受流场形状及流动特征的限制，具有较好的通用性。     

鲁棒的结构网格自动化重叠方法

网格技术是目前数值模拟中的关键技术之一。重叠网格是一种放宽拓扑要求、减小网格生成难度的网格技术。本文以结构重叠网格为基础,分别针对挖洞、寻点以及洞面优化方法进行了研究和改进,同时完成物面网格重叠,形成了一套鲁棒的、自动化的网格重叠系统。在挖洞方面,结合“最小洞映射”方法,提出“复合式挖洞”方法,节省内存开销;在寻点方面,通过构建格心虚网格,保证搜索空间的连续性,同时结合“有效搜索”思想,排除部分对寻点无贡献的网格点,进而减少ADT叉树节点;在洞面优化上,改变填补判别法则并引入两类受保护洞内点,确保两层插值边界建立,提高鲁棒性;在物面网格重叠上,利用物面投影法完成坐标修正,实现物面附近网格流动变量的准确传递。为验证本文方法,分别对定常翼身组合体DLR-F6绕流和非定常机翼挂载分离过程进行了数值模拟,计算结果与实验结果吻合良好,表明该结构重叠网格系统对多物体间定常、非定常扰流具有较好的数值模拟能力和较高的模拟精度,具有较高的工程应用价值。     

基于STL文件的笛卡尔网格生成方法研究

采用从STL文件中直接获取模型几何数据的方法,利用ADT数据结构快速剔除相同顶点,提出一种新的“推进查找”方法替代以往网格单元是否在物面内部的判定方法,避免了对物面的自封闭要求,无需进行几何重构,直接完成网格生成,更进一步提高了自适应笛卡尔网格生成速度和自动化程度,为实现CAD与CFD的之间的数据对接提供了一种新的方式。     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

多部件组合体分块网格生成技术及应用

给出多部件组合体分块网格生在的一种方案及实践,包括交界面网格生成技术、空间网格生成方法及质量控制。以一个平尾和机翼具有一定高度差的机身－机翼－平尾－立尾组合体外形为例,生成与公共交界面或边界面正交或接近正交的分块网格。用于Ｅｕｌｅｒ方程跨声速来流几个迎角的计算,结果与实验符合较好。     

网格生成和数值模拟的讨论

讨论了多种网格生成方法及其数值模拟。给出一种生成多部件复杂外形飞行器分区网格和Ｅｕｌｅｒ方程分区解。讨论了前缘分离涡为主体的旋涡流场。作了跨声速流分区的并行计算的初步尝试。给出一种二维非结构网格的生成方法。     

非结构网格生成技术的研究

本文描述了非结构网格生成的过程。在生成方法上综合Ｄｅｌａｕｎｅｙ和Ａｄｖａｎｃｉｎｇ Ｆｒｏｎｔ两种方法之长，构造了新的网格生成方法，使非结构网格生成更方便、快捷。网格生成过程分三个部分：１．表面网格的构造和最佳逼近；２．辅助点的生成与确定；３．空间四面体网格的生成。使用上述方法，本文构造限复杂构形物体的非结构网格。     

机翼机身跨声速绕流的计算

本文采用变换与数值网格生成方法，对由椭圆抛物而描述的每一个展向曲面进行网格生成，引入一种用于控制网格疏密分布的函数变换方法，其特点是易于对不同族网格线分别控制，根据物面梯度自适应地调整网格的疏密。用有限体积法计算守恒型全速势方程跨声速绕流，对M-100翼身组合体进行计算。理论计算结果与实验结果表明，本文方法有效。     

直接统计模拟位置元方法中的分子表面反射确定论判据

在直接模拟Ｍｏｎｔｅ－Ｃａｒｌｏ（ＤＳＭＣ）方法的位置元方法中准确判断分子是否与某一表面元相撞,直接关系着气动力热的计算精度,是方法需要解决的关键问题。本文发展的分子表面反射判据,使这个问题得到了有效的解决。此外,还讨论了分子空间位置坐标的记录方法、分子碰撞对的极限距离、网格自适应等实际模拟技术。与圆球自由分子流准确解和过渡领域贴体网格ＤＳＭＣ模拟结果的比较,证明了上述判据和技术的有效性。     

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。