Application of the GSMAC-CIP Method to Incompressible Navier-Stokes Equations at High Reynolds Numbers

The conventional shape function for the finite-element method (FEM) is linear, and it is thus inadequate for analyzing numerically complex flows at high Reynolds numbers. In this study, we propose a new scheme, GSMAC-CIP, using the third-order shape function, which requires continuity of the value of the function and its first space derivative in the whole space and is formulated by a finite element method for the cubic interpolated pseudo-particle (CIP) method. We verified the effectiveness of this new scheme by analyzing the forced-driven convection in a square cavity at Re = 1000, 5000 and 10000. The numerical results obtained by the present scheme are compared with those of GSMAC-FEM using coarser meshes, and it is shown that the present scheme is superior to GSMAC-FEM in terms of space accuracy. Moreover, it is shown that the numerical results obtained by the present scheme using fine meshes were in precise agreement with those obtained by Ghia et al.