计算流体力学非结构混合网格高精度格式研究现状

介绍了航空CFD技术中的非结构混合网格高精度格式的基本理论和方法,国内外研究现状,几种非结构混合网格高精度格式。     

基于非结构嵌套网格和逆风格式的旋翼悬停流场数值模拟

为更好地计及桨叶间相互影响,同时捕捉桨尖附近的流场细节,发展了一套非结构运动嵌套网格方法,并针对关键的贡献单元搜寻问题提出了一种高效的搜索方法.为了可方便地扩展到前飞状态的计算,同时为提高旋翼流场涡量的捕捉精度和效率,采用双时间法求解惯性坐标系下的非定常N-S主控方程,空间方向上使用二阶迎风格式,在此基础上,将不同精度格式的计算结果与试验值进行对比,并讨论了桨尖形状变化对旋翼气动特性的影响,结果表明本文建立的CFD方法对旋翼悬停流场的模拟是很有效的.     

非结构网格、混合网格下计算方法研究

为了提高计算精度，本文研究了计算格式和网格技术。通过比较非结构网格下的NND格式、MUSCL格式、ENO格式及最小二乘格式的计算精度和波动，选择利用NND方法中的守恒变量法来建立非结构网格下的计算格式，通过改进NND方法中的限制器来提高计算精度和分辨率，计算结果比较显示出该格式具有很高的分辨率和计算精度。本文也研究了非定常网格自适应技术，计算给出了双马赫反射问题的新旧限制器及网格自适应计算结果比较，新限制器网格自适应后计算得到的滑移线出现了剪切涡现象，这表明，数值计算中取足够密的网格和构造高精度格式是非常必要的。     

基于迎风格式的双曲型网格生成方法

双曲型网格生成方法是一种十分适合于重叠网格方法应用的网格生成方法。为了避免双曲型方程将初始外形的不光滑性继续向外传播，本文采用了自动隐含数值耗散项的迎风格式进行数值离散。在法向推进方向，采用了隐式格式，使得网格推进步长的大小不受稳定性条件的限制。对实际应用中可能出现的边界条件进行了分类和处理。几个典型外形网格生成的质量是令人满意的，说明本文方法是成功的。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

二维非结构网格上的高精度有限体积WENO格式

构造了非结构网格上二维双曲型守恒律的一类新的高精度有限体积WENO格式。其主要思想是:根据格式精度的要求,按照谱体积方法对三角形单元网格进行剖分,通过选取适当的子单元组成模板,利用WENO重构方法重构二阶和三阶多项式,利用有限体积公式和高阶Runge-Kutta TVD时间离散方法,构造了非结构网格上二维双曲型守恒律的一致二阶和三阶精度的有限体积WENO格式。然后,推广到二维Euler方程组。最后,给出几个数值算例,验证了格式的稳定性、高阶精度和高分辨捕捉激波等间断的能力。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

利用紧致格式捕捉间断的数值方法研究

将通量限制或加权的思想引入到紧致格式中,构造了两类传统方法与紧致格式混合组成的差分格式:通量限制与加权型差分格式。通过方波、组合波、定常激波、非定常Sod问题、Shu问题和Lax问题上的计算,以及与精确解的比较,结果表明这两种方法在间断的捕捉上具有高精度、高分辨率,而在计算方波、组合波或Shu问题和Lax问题上,加权格式具有更大的优势。     

非结构网格紧致高精度气体动理学格式

简要介绍了非结构网格上基于内自由度高效紧致重构方法的几种高精度气体动理学格式(GKS),包括结合三阶GKS通量求解器与子单元有限体积法(SCFV)及通量重构方法(CPR)构造的两种单步时空三阶格式SCFV-GKS和CPR-GKS,以及结合二阶GKS通量和两步四阶方法构造的时空四阶CPR-GKS.进而在CPR-GKS中混...     

非结构网格及混合网格复杂无粘流场并行计算方法研究

在分布式存储网络异构编程PVM环境下,本文将原串行计算程序改造为并行计算程序,提高了计算效率.利用非结构网格,对球体及航天飞机外形绕流流场进行了计算,取得了较好的性能结果.最后将该方法应用于混合网格的计算,求解了多种复杂外形的无粘绕流,如翼身组合体和DDF1C导弹,得到了满意的结果.     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

基于非结构自适应网格和嵌入边界法的残差格式研究

嵌入边界法由于在求解NS方程时能够简化网格生成问题而在计算流体领域受到越来越广泛的关注.简言之,嵌入边界法能够简化大变形和运动条件下多物理流动模拟、流固相互作用耦合问题,然而壁面边界条件的精确处理仍旧是该方法需要解决的问题.在本文工作中,为考虑壁面边界条件而在NS方程中增加了补偿项,同时采用非结构网格自适应技术保持了壁面边界条件的精度.     

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

一种非结构网格上基于径向基函数重构的ENO格式

基于二维Euler方程,对非结构三角形网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个三角形单元构造插值径向基函数,而在计算交界面的流通量采用两点高斯积分公式以保证格式的整体精度。时间离散采用三阶TVD Runge—Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

运输机组合体绕流非结构混合网格N-S方程数值模拟

本文利用通用网格生成软件生成运输类飞机翼身组合体（DLR-F6）非结构混合网格，在边界附面层区域内，生成三棱柱形网格，在其他流动区域采用四面体网格进行填充。粘性流场求解器采用基于N-S方程的CFX5计算流体软件求解，通过计算结果与风洞试验所得到的数据对比分析，验证了网格生成方法及网格数据的正确性和实用性。本文研究结果也为运输类飞机高速风洞试验提供了对比数据。     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。