Abstract:In this paper we present an all-at-once multigrid method for a distributed Stokes control problem (velocity tracking problem). For solving such a problem, we use the fact that the solution is characterized by the optimality system (Karush-Kuhn-Tucker-system). The discretized optimality system is a large-scale linear system whose condition number depends on the grid size and on the choice of the regularization parameter forming a part of the problem. Recently, block-diagonal preconditioners have been proposed, which allow to solve the problem using a Krylov space method with convergence rates that are robust in both, the grid size and the regularization parameter or cost parameter. In the present paper, we develop an all-at-once multigrid method for a Stokes control problem and show robust convergence, more precisely, we show that the method converges with rates which are bounded away from one by a constant which is independent of the grid size and the choice of the regularization or cost parameter.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the effective numerical solution to the Stokes control problem (especially the velocity tracking problem). Specifically, the author aims to develop a robust all - at - once multigrid method to solve the optimization problem constrained by the Stokes equations. The characteristic of such problems lies in that the optimality system (KKT system) after discretization is a large - scale linear system, and its condition number depends on the choice of grid size and regularization parameter. ### Core challenges of the problem 1. **Dependence of condition number**: The condition number of the discretized linear system depends not only on the grid size \(h\), but also on the regularization parameter \(\alpha\). This makes traditional iterative methods inefficient when dealing with problems of different scales. 2. **Computational complexity**: In order to ensure that the computational complexity is linear with the number of unknowns and does not depend on the grid level and parameter \(\alpha\), an efficient solver needs to be designed. ### Solutions in the paper The author proposes an all - at - once multigrid method, which can: - **Robust convergence**: It is proved that the convergence rate of this method is independent of the grid size and regularization parameter, that is, the convergence rate is bounded by a constant independent of these parameters. - **Efficient solution**: By combining the smoothing property and the approximation property, the efficiency and robustness of the algorithm in dealing with large - scale problems are ensured. ### Mathematical model Consider the following constrained optimization problem: \[ \min_{u, f} J(u, f)=\frac{1}{2}\|u - u_D\|_{L^2(\Omega)}^2+\frac{\alpha}{2}\|f\|_{L^2(\Omega)}^2 \] where \(u\in[H^1(\Omega)]^d\) is the velocity field, \(p\in L^2(\Omega)\) is the pressure distribution, \(f\in[L^2(\Omega)]^d\) is the control variable (force field), and satisfies the Stokes equations: \[ \begin{cases} -\Delta u+\nabla p = f & \text{in }\Omega,\\ \nabla\cdot u = 0 & \text{in }\Omega,\\ u = 0 & \text{on }\partial\Omega. \end{cases} \] ### Main contributions 1. **All - at - once multigrid method**: It is directly applied to the unreduced or reduced KKT system, avoiding the complexity in the block - iterative method. 2. **Convergence proof**: Based on the analysis of the classical smoothing property and approximation property, the robust convergence of this method under different grid levels and parameter selections is proved. 3. **Numerical verification**: The theoretical results are verified by numerical experiments, showing the effectiveness and robustness of the method. In conclusion, by developing a new all - at - once multigrid method, this paper solves the challenges brought by the dependence of condition number and computational complexity in the Stokes control problem, providing a powerful tool for research in related fields.

A robust all-at-once multigrid method for the Stokes control problem

Efficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems

Achieving $h$- and $p$-robust monolithic multigrid solvers for the Stokes equations

A staggered-grid multilevel incomplete LU for steady incompressible flows

A Symmetric Multigrid-Preconditioned Krylov Subspace Solver for Stokes Equations

Multigrid methods for the Stokes problem on GPU systems

Symbol based Convergence Analysis in Block Multigrid Methods with Applications for Stokes Problems

Monolithic Algebraic Multigrid Preconditioners for the Stokes Equations

Solution Techniques for the Stokes System: A priori and a posteriori modifications, resilient algorithms

Coarse-grid operator optimization in multigrid reduction in time for time-dependent Stokes and Oseen problems

RBF collocation and hybrid-LHI methods for Stokes systems and its application to controllability problems

A Multigrid Discretization of Discontinuous Galerkin Method for the Stokes Eigenvalue Problem

Multigrid preconditioning for discontinuous Galerkin discretizations of an elliptic optimal control problem with a convection-dominated state equation

A two-grid Adaptive Finite Element Method for the Dirichlet Boundary Control Problem Governed by Stokes Equation

An augmented Lagrangian preconditioner for the control of the Navier--Stokes equations

Scalable robust solvers for unstructured FE geodynamic modeling applications: Solving the Stokes equation for models with large localized viscosity contrasts

Optimal Control of the Multiphase Stefan Problem

Robust finite element solvers for distributed hyperbolic optimal control problems

Novel Approach for solving the discrete Stokes problems based on Augmented Lagrangian and Global Techniques: Application to Saddle-Point Linear Systems from Incompressible flow

Exploiting mesh structure to improve multigrid performance for saddle-point problems

Rotated block diagonal preconditioners for Navier-Stokes control problems