航空学报 > 2016, Vol. 37 Issue (12): 3605-3614   doi: 10.7527/S1000-6893.2016.0175

基于自适应混合网格的脱体涡模拟

张扬1,3, 张来平1,2, 赫新1,2, 邓小刚4   

  1. 1. 中国空气动力研究与发展中心 空气动力学国家重点实验室, 绵阳 621000;
    2. 中国空气动力研究与发展中心 计算空气动力研究所, 绵阳 621000;
    3. 中国空气动力研究与发展中心 低速空气动力研究所, 绵阳 621000;
    4. 国防科技大学, 长沙 410073
  • 收稿日期:2016-01-18 修回日期:2016-06-02 出版日期:2016-12-15 发布日期:2016-06-14
  • 通讯作者: 张来平,Tel.:0816-2463292,E-mail:zhanglp_cardc@126.com E-mail:zhanglp_cardc@126.com
  • 作者简介:张扬,男,博士研究生,工程师。主要研究方向:低速空气动力学计算与试验。Tel.:0816-2463205,E-mail:zhangy29v@sina.com;张来平,男,博士,研究员,博士生导师。主要研究方向:计算流体力学、非定常流动机理。Tel.:0816-2463292,E-mail:zhanglp_cardc@126.com
  • 基金资助:

    国家自然科学基金(11532016);国家科技支撑计划(2016YFB0200700)

Detached eddy simulation based on adaptive hybrid grids

ZHANG Yang1,3, ZHANG Laiping1,2, HE Xin1,2, DENG Xiaogang4   

  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;
    3. Low Speed Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    4. National University of Defense Technology, Changsha 410073, China
  • Received:2016-01-18 Revised:2016-06-02 Online:2016-12-15 Published:2016-06-14
  • Supported by:

    National Natural Science Foundation of China (11532016); National Key Technology Research and Development Program (2016YFB0200700)

摘要:

基于混合网格和CGNS(CFD General Notation System)数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟(DES)算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明:采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。

关键词: 网格自适应, 脱体涡模拟, 非结构网格, 混合网格, 三角翼, 有限体积法

Abstract:

An adaptive mesh technique with isotropic refining/coarsening approach based on CFD general notation system (CGNS) data structure is presented and implemented for hybrid grids. In order to simplify the possible refinement cases, elements with hanging nodes are changed into polyhedron and meanwhile the compatibility of the grid is maintained, so flow solvers that allow polyhedrons can operate on the adapted meshes without any modifications. Then, detached eddy simulation (DES) method combined with the grid adaptation technique is applied to simulate the flow over a 65° sweep delta wing at high angle of attack. The comparison of the DES results on initial grid and adaptive grid, as well as the experimental data, is carried out. The numerical results demonstrate that with the use of the adaptive technique the spatial resolution in the leeward side of the delta wing can be improved effectively due to the increment of cell number in some local region and the capability of present DES solver to resolve the small scale turbulent flow structure is effectively enhanced, consequently the issue of resolution is alleviated in DES simulation based on hybrid grid with a commonly used second-order scheme.

Key words: grid self-adaptive, detached eddy simulation, unstructured grid, hybrid grid, delta wing, finite-volume method

中图分类号: