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

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

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

二维非结构网格的可压缩Euler方程流场解

讨论了非结构网格的生成方法,用Delaunay三角化方法和推进阵面法相结合的方法生成了非结构网格,提高了生成非结构网格的自动化程度和网格的质量.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼形和两段翼型的算例,并与实验结果进行了比较和分析.     

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

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