计算机应用 ›› 2017, Vol. 37 ›› Issue (1): 42-47.DOI: 10.11772/j.issn.1001-9081.2017.01.0042

• 2016年全国开放式分布与并行计算学术年会(DPCS2016)论文 • 上一篇    下一篇

二维不可压缩Navier-Stokes方程的并行谱有限元法求解

胡园园1, 谢江1, 张武2   

  1. 1. 上海大学 计算机工程与科学学院, 上海 200444;
    2. 上海大学 高性能计算中心, 上海 200444
  • 收稿日期:2016-08-16 修回日期:2016-08-26 出版日期:2017-01-10 发布日期:2017-01-09
  • 通讯作者: 谢江
  • 作者简介:胡园园(1991-),女,江苏淮安人,CCF会员,主要研究方向:生物信息学、高性能计算;谢江(1971-),女,湖北恩施人,副教授,博士,CCF会员,主要研究方向:生物信息学、高性能计算;张武(1957-)男,江西武宁人,教授,博士,CCF会员,主要研究方向:高性能计算、生物信息学、计算流体学。
  • 基金资助:
    国家自然科学基金重大研究计划项目(91330116)。

Solution of two dimensional incompressible Navier-Stokes equation by parallel spectral finite element method

HU Yuanyuan1, XIE Jiang1, ZHANG Wu2   

  1. 1. School of Computer Engineering and Science, Shanghai University, Shanghai 200444, China;
    2. High Performance Computing Center, Shanghai University, Shanghai 200444, China
  • Received:2016-08-16 Revised:2016-08-26 Online:2017-01-10 Published:2017-01-09
  • Supported by:
    This work is supported by the Major Research Plan of the National Natural Science Foundation of China (91330116).

摘要: 针对不可压缩Navier-Stokes (N-S)方程求解过程中的有限元法存在计算网格量大、收敛速度慢的缺点,提出了基于面积坐标的三角网格剖分谱有限元法(TSFEM)并进一步给出了利用OpenMP对其并行化的方法。该算法结合谱方法和有限元法思想,选取具有无限光滑特性的指数函数取代传统有限元法中的多项式函数作为基函数,能够有效减少计算网格数量,提高算法的精度和收敛速度;利用面积坐标便于三角形单元计算的特点,选取三角单元作为计算单元,增强了适用性;在顶盖方腔驱动流问题上对该算法进行验证。实验结果表明,TSFEM较传统有限元法(FEM)无论是收敛速度还是计算效率都有了显著提高。

关键词: 不可压缩N-S方程, OpenMP, 方腔驱动流, 高精度, 无穷收敛性

Abstract: Due to a large number of computational grids and slow convergence existed in the numerical simulation of Navier-Stokes (N-S) equation, Triangular mesh Spectral Finite Element Method based on area coordinate (TSFEM) was proposed. And further, TSFEM was paralleled with OpenMP. Spectral method was combined with finite element method, and the exponential function with infinite smoothness was selected as the basis function to replace the polynomial function in the traditional finite element method, which can efficiently reduce the amount of computational grids as well as improve the convergence and accuracy of the proposed algorithm. Because area coordinates can facilitate the calculation of triangular units, which were selected as the computing units to enhance the applicability of the algorithm. The lid-driven cavity flow was used to verify the TSFEM. The experimental results show that, compared with the traditional Finite Element Method (FEM), the TSFEM greatly improves the convergence rate and the calculation efficiency.

Key words: incompressible Navier-Stokes (N-S) equation, OpenMP, lid-driven cavity flow, high-precision, infinite convergence

中图分类号: