Stable finite-element solution of the incompressible Navier-Stokes equations using linear interpolations for velocity and pressure

A velocity decomposition least squares method is reexamined and the results of several test cases are presented. All the present calculations are based on linear triangular elements, with a new treatment of boundary conditions. For high Reynolds number flows, second-order upwind schemes are adopted and the Lanczos method is applied to accelerate the convergence of the iterative scheme which is based on the symmetric Gauss-Seidel algorithm. Extension to unsteady flows is also presented. (C) 2001 Elsevier Science B.V. All rights reserved.