航空学报 > 2008, Vol. 29 Issue (5): 1118-1124

三维动态非结构重叠网格Navier-Stokes方程并行算法

夏健,田书玲,王江峰,伍贻兆   

  1. 南京航空航天大学 航空宇航学院
  • 收稿日期:2008-01-21 修回日期:2008-05-26 出版日期:2008-09-25 发布日期:2008-09-25
  • 通讯作者: 夏健

Parallel Computing Strategy for 3D Dynamic Overset Unstructured Navier-Stokes Solver

Xia Jian,Tian Shuling,Wang Jiangfeng,Wu Yizhao   

  1. College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics
  • Received:2008-01-21 Revised:2008-05-26 Online:2008-09-25 Published:2008-09-25
  • Contact: Xia Jian

摘要: 提出了一种三维动态非结构重叠网格Navier-Stokes(N-S)方程的并行计算方法。N-S方程的空间离散采用格点有限体积方法,时间离散采用隐式的双时间步长方法。应用一方程Spalart-Allmaras(S-A)模型来计算湍流黏性。并行计算采用动态的区域分裂方法,在每一物理时间步利用METIS网格分区系统对网格进行分区。为了实现各CPU之间的负载均衡,每块网格都按CPU个数进行分区并对活动节点和非活动节点进行了加权处理。最后,通过对外挂物投放无黏流动的数值模拟和内埋武器弹舱开启黏性流动的数值模拟,验证了该并行程序的准确性、高性能并行计算以及处理复杂几何外形的能力。

关键词: 非结构重叠网格, 并行计算, 区域分裂, Navier-Stokes 方程, 有限体积方法

Abstract: A parallel computing strategy for a 3D dynamic unstructured overset grids Navier-Stokes solver is developed in this article. The Navier-Stokes equations are discretized by a cell-vertex finite volume method in space and matched implicitly in time using a dual stepping method. The Spalart-Allmaras one-equation turbulence model is implemented to evaluate the turbulent viscosity. Parallel computation is based on the dynamic domain decomposition method which is performed using METIS system at each physical time step. Load balance is achievedby partitioning each component mesh according to the number of CPUs and weighing between active nodes and inactive nodes. Finally the parallel code is applied to inviscid flow simulation of store separation and viscous flow simulation of internal weapon bay opening to validate its accuracy, high parallel performance and the ability to deal with complex geometries.

Key words: unstructured overset grids, parallel computation, domain decomposition, Navier-Stokes equations, finite volume method

中图分类号: