Rotated block diagonal preconditioners for Navier-Stokes control problems

Abstract

Fluid flow control problems play a crucial role in various industrial applications. The optimal control of the Navier-Stokes equations poses a significant challenge. By employing Oseen's approximation to linearize the problem, we encounter a series of large sparse structured linear systems. These coefficient matrices exhibit a two-by-two block structure with square blocks. Leveraging this block structure, we exploit matrix splitting preconditioners. Theoretical analysis indicates that the real parts of all eigenvalues of preconditioned matrix are equal to 1/2. A large number of eigenvalues are clustered in (1±i)/2. The imaginary part of the eigenvalues is bounded explicitly. The eigenvalue distribution predicts the fast convergence of Krylov subspace methods. To avoid solving the saddle point subsystems in the preconditioning procedure, we propose a practical variant of the preconditioner. The theoretical analysis demonstrates that if the approximation is sufficiently close to Schur complement, the eigenvalues of the modified preconditioned matrix will lie within a circle centered at the original eigenvalues with a radius less than 1. This implies that the modified preconditioner exhibits excellent performance in terms of preconditioning. The numerical experiments demonstrate that the generalized minimal residual method, when combined with the proposed preconditioners, proves to be an efficient and effective solution for Oseen's control optimization with a variable viscosity coefficient. The number of iterations required by the preconditioned methods remains unaffected by the mesh size used in finite element discretization and only marginally depends on the regularization parameter.