Non-finite element implementation of the finite element projection in two dimensions

摘要

Embodiments of the present invention provide a finite difference method, a finite element projection, and a non-finite-element implementation of the finite element projection for solving the incompressible continuity equation on quadrilateral grids to obtain pressure and fluid velocity. An exemplary apparatus for simulating and analyzing the flow of an incompressible fluid is also provided. More particularly, embodiments of the present invention provide a non-finite-element-method (non-FEM) to implement the finite element projection to obtain a linear system whose coefficient matrix has an average bandwidth of five, in comparison to the finite difference scheme using conventional finite element projection and implementation with a bandwidth of nine. A two-dimensional-five-point-stencil-based non-FEM discretization is formulated to enforce the incompressible continuity equation. A shape function that is different from a weighting function used in the incompressible continuity equation is utilized to guarantee the five-point-stencil discretization is possible. The smaller bandwidth of the non-FEM scheme requires less computing resources and can compute the results in a reduced time duration, due to the simpler matrix to solve.