Unsteady solution of incompressible Navier-Stokes equations

The numerical scheme used by the present time-accurate FEM numerical method for incompressible Navier-Stokes equations, using primitive variables as the unknowns, is a Crank-Nicholson implicit treatment of all equation terms with central differencing for space derivatives. The introduction of a continuous auxilliary system in pseudo-time, with artificial compressibility, yields the incompressible solution at the advanced time level; time-accurate solutions are thereby obtained for two-dimensional fluid flows in a square cavity, in the cases of both an impulsively starting lid and an oscillating lid.