航空学报 > 2016, Vol. 37 Issue (10): 2980-2991   doi: 10.7527/S1000-6893.2016.0034

鲁棒的结构网格自动化重叠方法

王文1, 阎超1, 袁武2, 黄宇1, 席柯3   

  1. 1. 北京航空航天大学 航空科学与工程学院, 北京 100083;
    2. 中国科学院 计算机网络信息中心 超级计算中心, 北京 100190;
    3. 中国兵器工业导航与控制技术研究所, 北京 100089
  • 收稿日期:2016-01-03 修回日期:2016-01-18 出版日期:2016-10-15 发布日期:2016-02-22
  • 通讯作者: 阎超,Tel.:010-82317019 E-mail:yanchao@buaa.edu.cn E-mail:yanchao@buaa.edu.cn
  • 作者简介:王文 男,博士研究生。主要研究方向:计算流体力学,重叠网格方法。Tel:010-82338071 E-mail:wangwenbuaa@126.com;阎超 男,博士,教授,博士生导师。主要研究方向:空气动力学,计算流体力学。Tel:010-82317019 E-mail:yanchao@buaa.edu.cn

A robust and automatic structured overlapping grid approach

WANG Wen1, YAN Chao1, YUAN Wu2, HUANG Yu1, XI Ke3   

  1. 1. School of Aeronautic Science and Engineering, Beihang University, Beijing 100083, China;
    2. Supercomputing Center, Computer Network Information Center, Chinese Academy of Sciences, Beijing 100190, China;
    3. Institute of Ordnance Industry Navigation and Control Technology, Beijing 100089, China
  • Received:2016-01-03 Revised:2016-01-18 Online:2016-10-15 Published:2016-02-22

摘要:

网格技术是目前数值模拟中的关键技术之一。重叠网格是一种放宽拓扑要求、减小网格生成难度的网格技术。本文以结构重叠网格为基础,分别针对挖洞、寻点以及洞面优化方法进行了研究和改进,同时完成物面网格重叠,形成了一套鲁棒的、自动化的网格重叠系统。在挖洞方面,结合“最小洞映射”方法,提出“复合式挖洞”方法,节省内存开销;在寻点方面,通过构建格心虚网格,保证搜索空间的连续性,同时结合“有效搜索”思想,排除部分对寻点无贡献的网格点,进而减少ADT叉树节点;在洞面优化上,改变填补判别法则并引入两类受保护洞内点,确保两层插值边界建立,提高鲁棒性;在物面网格重叠上,利用物面投影法完成坐标修正,实现物面附近网格流动变量的准确传递。为验证本文方法,分别对定常翼身组合体DLR-F6绕流和非定常机翼挂载分离过程进行了数值模拟,计算结果与实验结果吻合良好,表明该结构重叠网格系统对多物体间定常、非定常扰流具有较好的数值模拟能力和较高的模拟精度,具有较高的工程应用价值。

关键词: 计算流体力学, 结构重叠网格, 洞映射, 寻点, 洞面优化方法, 翼身组合体, 机翼挂载分离

Abstract:

The grid technology is a crucial aspect of numerical simulation up till the present moment. Overlapping grid technology could relax the restriction of structure topology, thus reducing the difficulty of grid generation. Based on the structured overlapping grid, hole-cutting, donor search, and overlapping optimization are explored and improved, and surface grid assemblage is completed in this thesis. Simultaneously, a robust and automatic grid overlapping system is constructed. In order to save memory, based on the minimum hole mapping method, a composite hole-cutting method is proposed. Regarding donor search, the continuity of the search space is preserved with the virtual cell-centered grid. Moreover, the nodes of alternating digital tree (ADT) are decreased by an efficient search algorithm which excludes some non-contributing grid points. With respect to overlapping optimization, two kinds of protected hole-points are inserted and a new principle is developed to ensure the construction of two interpolated layers in paste phase, intensifying the robustness of overlapping optimization. The coordinates of near-body grids are corrected by projection method to make sure the accurate transmission of flow variables. With the enhanced algorithm, a steady DLR-F6 wing-body flow and an unsteady wing/pylon/store separation flow are performed, and excellent agreement of computational results compared with experimental data has been achieved. The algorithm shows remarkable capability and accuracy in the simulation of steady or unsteady multiple bodies flow and provides great application value for engineering.

Key words: computational fluid dynamics, structured overlapping grid, hole-mapping, donor search, hole optimization method, wing-body, wing/pylon/store separation

中图分类号: