A New Control Volume Finite Element Method for Three Dimensional Analysis of Polymer Flow

고분자 유동의 3차원 해석을 위한 새로운 검사 체적 유한 요소법

A new control volume finite element method is proposed for three dimensional analysis of polymer flow. Tetrahedral finite element is employed and co-located interpolation procedure for pressure and velocity is implemented. Inclusion of pressure gradient term in the velocity shape functions prevents the checkerboard pressure field from being developed. Vectorial nature of pressure gradient is considered in the velocity shape function so that velocity profile in the limit of very small Reynolds number becomes physically meaningful. The proposed method was verified through three dimensional simulation of pipe flow problem for Newtonian and power-law fluid. Calculated pressure and velocity field showed an excellent agreement with analytic solutions for pressure and velocity. Driven-cavity problem, which is reported to yield checkerboard pressure filed when conventional finite element method is applied, could be solved without yielding checkerboard pressure field when the proposed control volume finite element method was applied. The proposed method could be successfully applied to the three dimensional mold filling problem.