​

We describe a family of low-order finite element solvers for incompressible flows. Both advection and divergence terms are treated using consistent numerical fluxes along edges. Several techniques to accelerate convergence to steady state are explored and compared. The techniques are then used in a fully implicit time-marching scheme that solves a steady problem at every time step. Several examples demonstrate the usefulness of the developed schemes.

​