航空学报 > 2022, Vol. 43 Issue (5): 125170-125170   doi: 10.7527/S1000-6893.2021.25170

基于全线程树数据结构的笛卡尔网格高效生成技术

陈浩1,2, 毕林1,2, 华如豪1,2, 周清清1,2, 唐志共2, 袁先旭1,2   

  1. 1. 中国空气动力研究与发展中心 空气动力学国家重点实验室, 绵阳 621000;
    2. 中国空气动力研究与发展中心 计算空气动力研究所, 绵阳 621000
  • 收稿日期:2020-12-28 修回日期:2021-01-11 发布日期:2021-02-02
  • 通讯作者: 袁先旭 E-mail:yuanxianxu@cardc.cn
  • 基金资助:
    国家自然科学基金(12002358)

Efficient Cartesian mesh generation method based on fully threaded tree data structure

CHEN Hao1,2, BI Lin1,2, HUA Ruhao1,2, ZHOU Qingqing1,2, TANG Zhigong2, YUAN Xianxu1,2   

  1. 1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    2. Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China
  • Received:2020-12-28 Revised:2021-01-11 Published:2021-02-02
  • Supported by:
    National Natural Science Foundation of China (12002358)

摘要: 笛卡尔网格方法具有易于自适应、自动化程度高、网格质量好等优势。由于其非贴体特性,在处理复杂构型或者复杂流动问题时候往往存在网格量过大的问题,不容忽视。笛卡尔网格生成的效率直接影响了整个计算周期的长短,有必要发展高效的生成技术。对于笛卡尔网格而言,决定网格生成效率的关键在于网格数据结构,其直接影响计算量和存储量。针对三维构型进行笛卡尔网格生成,发展了邻居查询更快捷、内存利用率更高的全线程树数据结构,并在本文的方法框架下进行了适应性应用和改进。同时,为了高效地判断网格单元类型,构建了物面单元的快速检索方式,并引入了染色方法,进一步提高网格生成效率。还提出了一种鲁棒的奇异性检测算法,保证网格单元类型判断的鲁棒性。在流场解自适应方面,采用的是速度散度和旋度相结合的三维判据,以保证对于多种流动特征的捕捉能力。通过圆球、导弹、翼身组合体、机翼-副翼等算例进行了考核验证,经对比,网格自适应位置与理论解吻合较好,且网格单元生成耗时短、平均耗时受物面网格分布影响小,证明了方法的可靠和高效性。

关键词: 笛卡尔网格, 全线程树, 奇异性检测算法, 流场解自适应, 染色算法

Abstract: The Cartesian mesh method has the advantages of easy adaptation, high degree of automation and good quality. However, its non-body fitted characteristics leads to unnegligible massive amount of grids when addressing complex configuration or complex flow problems. The efficiency of Cartesian mesh generation directly affects the whole computation cycle, entailing the development of efficient generation technology. The key to Cartesian grid generation efficiency is the grid data structure, which directly influences the amount of computation and storage. In this paper, Cartesian grid generation is conducted for three-dimensional configurations, and the Fully Threaded Tree data structure with faster neighbor query and higher memory utilization is developed, which is applied and further improved in this framework. Meanwhile, to effectively determine the type of grid cells, a fast retrieval method of facets is constructed, and the painting algorithm introduced to further improve the grid generation efficiency. In addition, a robust singularity detection algorithm is proposed to ensure the robustness of the grid cell type determination. With respect to mesh adaption based on the flow field solution, a three-dimensional criterion combining velocity divergence and curl is adopted to ensure the ability to capture various flow characteristics. Some test examples such as spheres, missiles, wing-body and wing-aileron configurations are used to verify the method. Comparison shows that the adaptive position of the mesh is in good agreement with the theoretical solution, the generation time cost of the mesh is low, and the average time is less affected by the distribution of the facets, thereby proving the reliability and efficiency of the proposed method.

Key words: Cartesian mesh, fully threaded tree, singularity check algorithm, flow feature based adaptation, painting algorithm

中图分类号: