航空学报 > 2015, Vol. 36 Issue (5): 1480-1490   doi: 10.7527/S1000-6893.2014.0331

飞翼气动优化中参数化和网格变形技术

唐静1, 邓有奇2, 马明生1,2, 李彬2   

  1. 1. 西北工业大学 航空学院, 西安 710072;
    2. 中国空气动力研究与发展中心 计算空气动力研究所, 绵阳 621000
  • 收稿日期:2014-07-16 修回日期:2014-12-01 出版日期:2015-05-15 发布日期:2014-12-11
  • 通讯作者: 唐静Tel.: 0816-2463270 E-mail: tangjingn@hotmail.com E-mail:tangjingn@hotmail.com
  • 作者简介:唐静 男, 博士研究生, 助理研究员。主要研究方向: 飞行器气动外形优化设计, 非结构网格技术。 Tel: 0816-2463270 E-mail: tangjingn@hotmail.com;邓有奇 男, 博士, 研究员, 博士生导师。主要研究方向: 计算空气动力学, 飞行器优化设计。 Tel: 0816-2463270 E-mail: cai@cardc.cn;马明生 男, 博士, 研究员, 博士生导师。主要研究方向: 飞行器气动总体设计。 Tel: 0816-2463090 E-mail: cai@cardc.cn

Parameterization and grid deformation techniques for flying-wing aerodynamic optimization

TANG Jing1, DENG Youqi2, MA Mingsheng1,2, LI Bin2   

  1. 1. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China;
    2. Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China
  • Received:2014-07-16 Revised:2014-12-01 Online:2015-05-15 Published:2014-12-11

摘要:

几何参数化和网格变形是飞行器气动外形数值优化迭代过程中的两个关键技术。基于非均匀有理B样条(NURBS)的自由型面变形(NFFD)技术对几何表示形式具有普适性,距离权函数(DWF)网格变形技术具有计算快速和网格拓扑无关性,两者广泛应用于曲面优化。基于NFFD技术改进了反求参数的Newton迭代算法,并通过提高物面附近网格刚度改进的距离权函数(IDWF)技术使其适用于更大程度的网格变形。还提出了改进后的参数化和网格变形两种技术并行计算的具体实现算法。结合离散伴随方法,使用参数化和网格变形技术,实现了由NACA0012初始翼型到飞翼标准翼型EH1590的反设计;针对某飞翼标模完成了单点全机升阻比优化,升阻比提高约18%。数值结果表明,建立的NFFD和IDWF动网格技术可满足飞翼气动外形优化参数化和快速网格变形的需求。

关键词: 外形参数化, 自由型面变形, 距离权函数, 动网格, 飞翼, 外形优化

Abstract:

Geometry shape parameterization and grid deformation techniques are two key issues for aerodynamic shape optimization of aircraft by CFD. The non-uniform rational B-splines (NURBS)-based free-form deformation (NFFD), which is universal for arbitrary representation of geometry, and the dynamic grid technique based on distance weighted function (DWF), which is independent from grid topology and fast on calculation, are widely used for shape optimization. In this paper, the NFFD technique is introduced in detail and the Newton iteration algorithms for inverse local coordinates calculation are developed. The DWF is improved (IDWF) to satisfy larger range grid deformation by increasing the stiffness of grid cells near surface. The parallel implementation algorithms for both techniques are designed and introduced in detail. With an optimization solver based on gradient calculated by discrete adjoint method, the standard flying-wing airfoil EH1590 is inversely designed form the initial airfoil NACA0012, and the ratio of lift and drag is optimized and increased by 18% for a whole flying-wing aircraft in a single design state. The result indicates that for flying-wing aerodynamic shape optimization, the NFFD technique is sufficient for shape parameterization and IDWF technique is efficient for grid deformation.

Key words: shape parameterization, free-form deformation, distance weighted function, dynamic grid, flying-wing, shape optimization

中图分类号: