On a network method for unsteady incompressible fluid flow on triangular grids

The dual variable method for Delaunay triangulations is a network‐theoretic method that transforms a set of primitive variable finite difference or finite element equations for incompressible flow into an equivalent system which is one‐fifth the size of the original. Additionally, it eliminates the pressures from the system and produces velocities that are exactly discretely divergence‐free. In this paper new discretizations of the convection term are presented for Delaunay triangulations, the dual variable method is extended to tessellations that contain obstacles, and an efficient algorithm for the solution of the dual variable system is described.