The dual variable method for solving fluid flow difference equations on Delaunay triangulations

Abstract

A complementary volume discretization of the transient Navier-Stokes equations on a triangular mesh is viewed as a system defining flows on an associated network. The dual variable method transforms this primitive system into an equivalent system which is one-fifth the size of the primitive system, eliminates the pressures and yields velocities that are exactly discretely divergence free.