航空学报 > 2016, Vol. 37 Issue (8): 2395-2407   doi: 10.7527/S1000-6893.2016.0141

非结构变形网格和离散几何守恒律

刘君1, 刘瑜2, 陈泽栋3   

  1. 1. 大连理工大学 航空航天学院, 大连 116024;
    2. 装备学院 航天装备系, 北京 101416;
    3. 大连理工大学 工程力学系, 大连 116024
  • 收稿日期:2016-01-16 修回日期:2016-05-03 出版日期:2016-08-15 发布日期:2016-05-10
  • 通讯作者: 刘君,Tel.:0411-84707176。E-mail:liujun65@dlut.edu.cn E-mail:liujun65@dlut.edu.cn
  • 作者简介:刘君,男,博士,教授,博士生导师。主要研究方向:空气动力学。Tel:0411-84707176。E-mail:liujun65@dlut.edu.cn;刘瑜,男,博士,讲师。主要研究方向:燃烧数值模拟。E-mail:liuyu@nudt.edu.cn;陈泽栋,男,博士研究生。主要研究方向:计算流体力学,空气动力学。E-mail:chenzd_dut@163.com
  • 基金资助:

    国家自然科学基金(91541117)

Unstructured deforming mesh and discrete geometric conservation law

LIU Jun1, LIU Yu2, CHEN Zedong3   

  1. 1. School of Aeronautics and Astronautics, Dalian University of Technology, Dalian 116024, China;
    2. Department of Space Equipment, Academy of Equipment, Beijing 101416, China;
    3. Department of Engineering Mechanics, Dalian University of Technology, Dalian 116024, China
  • Received:2016-01-16 Revised:2016-05-03 Online:2016-08-15 Published:2016-05-10
  • Supported by:

    National Natural Science Foundation of China (91541117)

摘要:

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉(ALE)方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格(MSA)及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数(RBFs)和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律(DGCL)概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。

关键词: 非结构网格, 网格变形, 非定常流动, 几何守恒律, 界面耦合

Abstract:

A popular method for simulating unsteady flow of fluid-structure interaction or multi-body separation is finite volume method based on arbitrary Lagrangian-Eulerian (ALE) equations, which involves mesh deformation and discrete geometric conservation law. The fundamental principle, state of the art and boundaries of validity of mesh deformation are reviewed following 5 categories of constructing ideas, e.g., physics analogy, ellipse smoothing, interpolation, moving submesh approach (MSA) and hybrid method. A hybrid method which combines the benefits of MSA and radial basis functions (RBFs) interpolation is proved to be robust and efficient via several numerical examples. After the concept of discrete geometric conservation law (DGCL) is introduced, the intrinsic mechanism of DGCL is analyzed through 2D model, which is the inequality between volume increment and the volume sweeping by the surfaces enclosing of mesh cell. The different implementations of DGCL which could be mainly categorized as area correction, area correction via assuming velocity of surface, velocity correction and volume correction are surveyed and their range of application and the existing problems are discussed. We found that the proposed volume correction method can satisfy the fluid-structure interface condition, and is also appropriate for multi-step temporal discrete schemes.

Key words: unstructured mesh, mesh deformation, unsteady flow, geometric conservation law, interface coupling

中图分类号: