三维复杂外形的双曲型非结构网格生成方法

为了实现棱柱网格的自动化生成,发展了使用双曲型偏微分方程生成棱柱网格的方法。利用网格正交性条件和单元体积约束相结合的双曲偏微分方程,从物面网格向外逐步推进生成棱柱形网格。通过类比结构网格双曲型网格生成方法,在非结构网格局部格点坐标系上建立了双曲型网格生成方程,通过格点有限体积法、迎风格式和隐式的方法离散控制方程,使用GMRES求解器求解线性方程组,进而自动生成棱柱网格。该方法应用于球外形、具有尖锐边缘、深凹高凸的旋成体外形、双椭球外形以及X38复杂飞行器,生成的网格表现出良好的正交性和光滑性,表明本文的方法是成功的。数值仿真结果表明本文的网格适用于气动力热的数值模拟。     

基于迎风格式数值模拟三角翼涡流场

本文应用基于有限体积方法的高精度迎风格式模拟了三角翼在大迎角条件下的绕流流场。选用不同几何外形的三角翼，由计算结果可以清楚的看到前缘自由剪切层的发展，涡的形成及其复杂的相互作用。     

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

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

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

复杂外形网格生成技术

复杂外形网格生成是进行流场数值模拟的一个关键技术及前提。本文论述了采用保角变换方法,求解偏微分方程方法以及它们在代数方法结合生成二维及三维具有多块网格结构的贴体且与边界正交网格的技术。     

复杂空气动力外形网格生成技术

通过对用于网格生成的Ｌａｐｌａｃｅ方程进行抛物化差分离散,发展了一种快速抛物型方程贴体及与边界正交网格生成方法,另外,采用源项修正技术求解用于网格生成的Ｐｏｉｓｓｏｎ方程,发展了一种可实现对网格线与边界间距及正交性进行双重控制的椭圆型方程方法,运用这种两种基本方法以及这两种方法的混合,分别生成了真实歼击机全机外形,带十字叉翼导弹以及双发,下单翼运输机外形的计算网格。并已成功地秀于跨声速得杂流场的娄     

飞机外形网格拓扑和多块网络拼接技术

多块网格方法成功的关键在于网格拓扑的选择和块与块之间的拼接方法,针对真实飞机外形多块网格生成,本文分析了研究了翼类和体类部件的各种网格拓扑,论述了翼身组合体的多块网格非结构连续拼接技术,指出连续拼接方法的缺点,以及相应的解决方法;同时还讨论了分块策略和确保块与块之间网格线安全连续拼接的方法,通过实例进一步阐述了本文工作的另外一些独到之处。     

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分(PFD)格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过类似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分(FPFD)格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。     

求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

上风格式的若干性能分析

采用理论分析和数值试验相结合的方法，对上风格式的三个典型代表：FDS、FVS及AUSM格式的性能进行了分析比较：FDS格式具有高的间断分辨率但可能出现非物理懈；FVS格式捕捉强非线性波的能力很强，计算效率高，但耗散比较大；AUSM格式兼有高间断分辨率和高计算效率且克服了两者的弱点。以一维无粘激波管流场为算例考察了这些格式性能与间断分辨率。     

基于高阶逆风通量差分裂格式的直升机旋翼前飞流场模拟

描述了三维雷诺平均NS方程数值模拟前飞状态旋翼非定常流场的方法和模型.在该模型中,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,无需添加人工粘性,因而可有效地减少Jameson二阶中心差分有限体积法导致的旋翼尾迹数值耗散.为充分考虑旋翼尾迹对流场的影响和计入桨叶的旋转、挥舞、变距运动,发展了一种运动嵌套网格.针对嵌套网格中关键的贡献单元搜寻问题,提出了一种新的"伪贡献单元搜寻法".应用上述方法,对Caradonna模型旋翼在周期变距和挥舞运动条件下的前飞流场进行了计算,与采用Jameson二阶中心差分方法相比,计算结果更贴近实验值.并针对几种不同桨尖形状的旋翼前飞流场进行计算,表明后掠桨尖可有效减弱超临界流动.     

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

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

一种基于约束框架的棱柱网格生成方法

提出了一种基于约束框架的棱柱网格生成方法,即利用物面网格节点与背景结构框架网格之间的单一映射关系线性插值生成棱柱网格。背景结构框架网格基于TFI方法生成。物面网格节点与背景框架网格之间的映射关系通过一种线性投影方法建立。然后通过V型槽、圆球、M6机翼和F6翼身组合体四种不同外形对该棱柱网格生成方法进行了测试,测试结果表明该方法能够避免棱柱网格生成中的网格交叉现象,保证棱柱网格的成功生成。该方法在棱柱网格质量控制、棱柱网格局部修改等方面具有较好的效率和优势,但背景框架网格的生成依然是该方法的一个瓶颈,需要进一步的研究。     

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

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

双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。     

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

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

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