Resolution of the implicit Euler scheme for the Navier–Stokes equation through a least-squares method

Abstract

This work analyzes a least-squares method in order to solve implicit time schemes associated to the 2D and 3D Navier–Stokes system, introduced in 1979 by Bristeau, Glowinksi, Periaux, Perrier and Pironneau. Implicit time schemes reduce the numerical resolution of the Navier–Stokes system to multiple resolutions of steady Navier–Stokes equations. We first construct a minimizing sequence (by a gradient type method) for the least-squares functional which converges strongly and quadratically toward a solution of a steady Navier–Stokes equation from any initial guess. The method turns out to be related to the globally convergent damped Newton approach applied to the Navier–Stokes operator. Then, we apply iteratively the analysis on the fully implicit Euler scheme and show the convergence of the method uniformly with respect to the time discretization. Numerical experiments for 2D examples support our analysis.