On the convergence order of the finite element error in the kinetic energy for high Reynolds number incompressible flows

Abstract

The kinetic energy of a flow is proportional to the square of the L-2 (Omega) norm of the velocity. Given a sufficient regular velocity field and a velocity finite element space with polynomials of degree r, then the best approximation error in L-2 (Omega) is of order r + 1. In this survey, the available finite element error analysis for the velocity error in L-infinity(0, T; L-2 (Omega)) is reviewed, where T is a final time. Since in practice the case of small viscosity coefficients or dominant convection is of particular interest, which may result in turbulent flows, robust error estimates are considered, i.e., estimates where the constant in the error bound does not depend on inverse powers of the viscosity coefficient. Methods for which robust estimates can be derived enable stable flow simulations for small viscosity coefficients on comparatively coarse grids, which is often the situation encountered in practice. To introduce stabilization techniques for the convection-dominated regime and tools used in the error analysis, evolutionary linear convection-diffusion equations are studied at the beginning. The main part of this survey considers robust finite element methods for the incompressible Navier-Stokes equations of order r -1, r, and r 1/2 for the velocity error in L-infinity(0, T; L-2 (Omega)). All these methods are discussed in detail. In particular, a sketch of the proof for the error bound is given that explains the estimate of important terms which determine finally the order of convergence. Among them, there are methods for inf-sup stable pairs of finite element spaces as well as for pressure-stabilized discretizations. Numerical studies support the analytic results for several of these methods. In addition, methods are surveyed that behave in a robust way but for which only a non-robust error analysis is available. The conclusion of this survey is that the problem of whether or not there is a robust method with optimal convergence order for the kinetic energy is still open.