CHARACTERISTICS OF MATRICES IN THE P2P1 FINITE ELEMENT METHODS FOR SOLVING THE INCOMPRESSIBLE NAVIER-STOKES EQUATION

P2P1 유한요소를 이용한 비압축성 Navier-Stokes 방정식 해법들의 행렬 특성

Numerical algorithms for solving the incompressible Navier-Stokes equations using P2P1 finite element are compared regarding the eigenvalues of matrices. P2P1 element allocates pressure at vertex nodes and velocity at both vertex and mid nodes. Therefore, compared to the P1P1 element, the number of pressure variables in the P2P1 element decreases to 1/4 in the case of two-dimensional problems and to 1/8 in the three-dimensional problems. Fully-implicit-integrated, semi-implicit- integrated and semi-segregated finite element formulations using P2P1 element are compared in terms of elapsed time, accuracy and eigenvlue distribution (condition number). For the comparison,they have been applied to the well-known benchmark problems. That is, the two-dimensional unsteady flows around a fixed circular cylinder and decaying vortex flow are adopted to check spatial accuracy.