有相对运动多体动力学系统流动计算方法若干问题讨论

对应用非结构动网格技术数值模拟有相对运动多体动力学系统的流动现象时发现的若干问题进行了讨论.应用动网格技术求解均匀流场,验证了几何守恒律方程离散形式对流场计算精度的重要性,还发现格式时间精度越高、网格运动速度越接近声速,几何守恒律的作用越明显;采用类似于可控轨迹实验的准定常计算方法和非结构动网格技术,分别对超声速气流中物体运动的动力学过程进行了数值模拟,通过流场和运动特性比较和理论分析,指出准定常计算方法本质上无法反映非定常流动特征,即使时间步长趋向零也不会趋向于真解;还讨论了数值模拟非定常流动的双时间法,分析时间推进步长对解的影响,随着时间步长减小解的误差降低,但是这一方法与动网格技术结合解决多体分离问题有难度.     

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

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

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

径向基函数在三维Euler方程数值计算中的应用

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

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格,该方法的特点是在生成单元的同时引入结点,与其他预布结点的方法相比,其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构,有效地提高了程序运行效果,最后给出了Ｍ６机翼的数值算例。     

飞机外挂物离机过渡过程的有限元算法.

研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中，流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量，而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况，只要将固定网格的守恒通量项Fk用运动网格的Ek来代替，并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。