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

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

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较.     

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。      

基于非结构网格的离子推力器羽流数值模拟方法研究

为适应复杂的航天器外形,文章采用基于非结构网格的数值模拟方法研究离子推力器真空羽流的基本特征。首先,介绍了非结构网格划分的方法,并分析了建模区域生成各种拓扑关系的特点;然后,提出了粒子穿越网格的具体算法,包括有利于减少计算量的结合矩阵重标号的带宽压缩算法;最后,开展了模拟结果和试验结果的对比研究,验证了以上数值模拟方法的可行性。     

三维分块非结构化网格上的可压缩流计算方法

提出了分块非结构化网格生成方法,通过子块划分、子块网格划分、整体网格拓扑连接生成可局部加密的分块非结构化网格,既具备了非结构化网格对复杂区域的处理能力,又避免了分块结构化网格块间信息传递的复杂计算过程.通过将密度修正值引入压力修正方程,发展了非结构化网格上可压缩流动的SIMPLE算法.对孤立翼型及透平叶栅的流动进行了数值计算,压力系数、速度、折转角的数值结果与已有的实验值吻合较好,验证了数值方法的精度.      

基于非结构网格有限差分法的扎染算法

应用有限差分法时遵循贴体坐标系下离散等价方程的离散准则,格式中仅用到当地网格点的坐标变换系数,据此可以构建新的非结构网格有限差分法,数值算例表明在空间二维情况下用连接离散点的3条网格线构造出来的一阶迎风格式可以稳定收敛。把这种非结构网格有限差分法推广到三维,在4条网格线基础上进行离散计算,然后提出一种局部区域多重网格、重复计算的新算法。首先在包含物体的整个区域采用直角坐标系下的均匀网格进行计算,然后消除物体内部点对物体外部流场区域的影响,最后在物面和直角网格之间很小的局部区域填补非结构网格进行计算,计算过程类似于中国传统的扎染工艺。处理流场内存在曲面物体边界问题时,相较于常规结构网格差分算法,扎染算法尽管在计算中包含需要特殊处理的无用网格点,也存在重复计算的局部区域,但是以这些计算能力的浪费为代价,换取了编程简便、内存量小、网格生成快捷、易于扩展网格规模等优势。定性分析表明,这种扎染算法非常适合研制大规模并行计算,建立的计算流体力学应用软件可以充分发挥数十万至百万处理器核心的超级计算机的效能。     

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

基于非结构动网格的跨声速非定常流场求解

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.采用改进后的方法并通过耦合求解了基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的刚性无后掠机翼及某翼身组合体的跨声速绕流.结果表明,改进后的方法可以大大提高网格变形能力和网格质量,且计算结果与实验结果吻合良好.     

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

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