航空学报 > 2004, Vol. 25 Issue (3): 229-232

并行计算中一种非结构网格分割方法

周春华   

  1. 南京航空航天大学空气动力学系 江苏南京 210016
  • 收稿日期:2003-05-23 修回日期:2004-02-12 出版日期:2004-06-25 发布日期:2004-06-25

A Method of Non-structured Mesh Partition for Parallel Computation

ZHOU Chun-hua   

  1. Department of Aerodynamics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
  • Received:2003-05-23 Revised:2004-02-12 Online:2004-06-25 Published:2004-06-25

摘要: 将递归谱对剖分方法应用于流体力学并行计算中的非结构网格分割,以解决负载平衡和最小切割问题。为使用这种方法,计算了网格伴随图的离散Laplacian矩阵的第二特征矢,然后从该特征矢的分量引入网格的对剖分。特征矢计算中应用了Rayleigh商迭代,并进行了一些修正以使收敛强烈地偏向于第二特征矢及考虑逆迭代步中线性方程组的迭代求解。最后,通过非结构自适应网格上Euler方程分区计算的数值结果验证了所发展的网格分割方法。

关键词: 计算流体力学, 并行计算, 递归谱对剖分, 负载平衡, 网格分割, Rayleigh商迭代

Abstract: A recursive spectral bisection method is applied to partition non-structured triangular meshes in parallel CFD for the treatment of load balance and minimum cut graph bisection. To use this method, the second eigenvector of the discrete Laplacian matrix of the dual graph of a mesh is computed and the bisection is obtained from components of this eigenvector. A classical Rayleigh quotient iteration is applied in the computation of the second eigenvector, and some modifications are made both to bias convergence to the second eigenvector and to account for the solution of the linear system in the inverse iteration step by an iterative process. Finally, the mesh partition approach is validated in the numerical experiments for solving Euler equations in parallel on the adaptive non-structured meshes.

Key words: computational fluid dynamics(CFD), parallel computation, recursive spectral bisection, load balance, mesh partition, Rayleigh quotient iteration