Convergence analysis of a finite element projection/Lagrange-Galerkin method for the incompressible Navier-Stokes equations

This paper provides a convergence analysis of a fractional-step method to compute incompressible viscous flows by means of finite element approximations. In the proposed algorithm, the convection, the diffusion, and the incompressibility are treated in three different substeps. The convection is treated first by means of a Lagrange-Galerkin technique, whereas the diffusion and the incompressibility are treated separately in two subsequent substeps by means of a projection method. It is shown that provided the time step, delta t, is of O(h(d/4)), where h is the meshsize and d is the space dimension (2 less than or equal to d less than or equal to 3), the proposed method yields for finite time T an error of O(h(l+1) + delta t) in the L-2 norm for the velocity and an error of O(h(l) + delta t) in the H-1 norm (or the L-2 norm for the pressure), where l is the polynomial degree of the approximate velocity.