航空学报 > 2014, Vol. 35 Issue (3): 727-735   doi: 10.7527/S1000-6893.2013.0300

基于径向基函数改进的Delaunay图映射动网格方法

孙岩1, 邓小刚1,2, 王光学1, 王运涛1, 毛枚良1   

  1. 1. 中国空气动力研究与发展中心 空气动力学国家重点实验室, 四川 绵阳 621000;
    2. 国防科学技术大学 计算机学院, 湖南 长沙 410073
  • 收稿日期:2013-04-19 修回日期:2013-06-06 出版日期:2014-03-25 发布日期:2013-07-20
  • 通讯作者: 孙岩,Tel.:0816-7067914 E-mail:supersunyan@163.com E-mail:supersunyan@163.com
  • 作者简介:孙岩 男,博士研究生,助理研究员。主要研究方向:计算流体力学。Tel:0816-7067914 E-mail:supersunyan@163.com;邓小刚 男,博士,研究员,博士生导师。主要研究方向:计算流体力学。E-mail:xgdeng@skla.cardc.cn;王光学 男,硕士,副研究员。主要研究方向:计算空气动力学和大迎角分离流。E-mail:wgx111@sina.com
  • 基金资助:

    国家“863”计划(2012AA051304);国家自然科学基金(11372337)

Improvement on Delaunay Graph Mapping Dynamic Grid Method Based on Radial Basis Functions

SUN Yan1, DENG Xiaogang1,2, WANG Guangxue1, WANG Yuntao1, MAO Meiliang1   

  1. 1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    2. School of Computer Science, National University of Defense Technology, Changsha 410073, China
  • Received:2013-04-19 Revised:2013-06-06 Online:2014-03-25 Published:2013-07-20
  • Supported by:

    National High-tech Research and Development Program of China (2012AA051304); National Natural Science Foundation of China (11372337)

摘要:

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。

关键词: 径向基函数, 多项式, Delaunay, 动网格, 衰减函数, 翼型

Abstract:

An improvement on the Delaunay graph mapping dynamic gird method based on radial basis function interpolation is presented in this paper to enhance the efficiency and quality of dynamic grid generation. First, the displacements of far-field background grid nodes are calculated by means of radial basis function interpolation with polynomial terms, and then the displacements of computational grid nodes are acquired with the interpolation of background grid nodes displacements. Finally, the displacements of the far-field computational grid nodes are decreased to zero with the use of damp functions. Differences between the improved method and the original Delaunay graph mapping method are investigated based on a grid deformation example of rectangle rotation. The effect of the number of radial basis functions'center points and damp functions on deformed grid quality are also investigated. The grid deformation for rectangle rotation shows that the improved method can accurately recover the rotation of the Delaunay background grid. Dynamic grid generation for NACA 0012 airfoil, NLR 7301 bi-element airfoil and M6 wing further demonstrates that with the insertion of control points the improved method can generate high quality dynamic grids for large deformation without the regeneration of the background grid.

Key words: radial basis functions, polynomial, Delaunay, dynamic grid, damp function, airfoil

中图分类号: