Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems

The governing dynamics of fluid flow is stated as a system of partial differential equations referred to as the Navier-Stokes system. In industrial and scientific applications, fluid flow control becomes an optimization problem where the governing partial differential equations of the fluid flow are stated as constraints. When discretized, the optimal control of the Navier-Stokes equations leads to large sparse saddle point systems in two levels. In this paper, we consider distributed optimal control for the Stokes system and test the particular case when the arising linear system can be compressed after eliminating the control function. In that case, a system arises in a form which enables the application of an efficient block matrix preconditioner that previously has been applied to solve complex-valued systems in real arithmetic. Under certain conditions, the condition number of the so preconditioned matrix is bounded by 2. The numerical and computational efficiency of the method in terms of number of iterations and execution time is favorably compared with other published methods.

[1]  Matthias Heinkenschloss,et al.  Local Error Estimates for SUPG Solutions of Advection-Dominated Elliptic Linear-Quadratic Optimal Control Problems , 2010, SIAM J. Numer. Anal..

[2]  A. Wathen,et al.  FAST ITERATIVE SOLVERS FOR CONVECTION-DIFFUSION CONTROL PROBLEMS ∗ , 2013 .

[3]  Howard C. Elman,et al.  Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics , 2014 .

[4]  J. Lions.  Optimal Control of Systems Governed by Partial Differential Equations , 1971 .

[5]  G. Golub,et al.  Inexact and preconditioned Uzawa algorithms for saddle point problems , 1994 .

[6]  J. Cahouet,et al.  Some fast 3D finite element solvers for the generalized Stokes problem , 1988 .

[7]  Andrei Draganescu,et al.  Multigrid solution of a distributed optimal control problem constrained by the Stokes equations , 2011, Appl. Math. Comput..

[8]  Andrew J. Wathen,et al.  Preconditioning Iterative Methods for the Optimal Control of the Stokes Equations , 2011, SIAM J. Sci. Comput..

[9]  Roland Becker,et al.  Optimal control of the convection-diffusion equation using stabilized finite element methods , 2007, Numerische Mathematik.

[10]  Anne Greenbaum,et al.  Any Nonincreasing Convergence Curve is Possible for GMRES , 1996, SIAM J. Matrix Anal. Appl..

[11]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[12]  W. Bangerth,et al.  deal.II—A general-purpose object-oriented finite element library , 2007, TOMS.

[13]  Walter Zulehner,et al.  Efficient Solvers for Saddle Point Problems with Applications to PDE–Constrained Optimization , 2013 .

[14]  Owe Axelsson,et al.  A comparison of iterative methods to solve complex valued linear algebraic systems , 2014, Numerical Algorithms.

[15]  Zhong-Zhi Bai,et al.  Block preconditioners for elliptic PDE-constrained optimization problems , 2011, Computing.

[16]  Apostol T. Vassilev,et al.  Analysis of the Inexact Uzawa Algorithm for Saddle Point Problems , 1997 .

[17]  A. Wathen,et al.  Chebyshev semi-iteration in preconditioning for problems including the mass matrix. , 2008 .

[18]  Stefan Ulbrich,et al.  Optimization with PDE Constraints , 2008, Mathematical modelling.

[19]  Maxim A. Olshanskii,et al.  Uniform preconditioners for a parameter dependent saddle point problem with application to generalized Stokes interface equations , 2006, Numerische Mathematik.

[20]  J. Pearson.  On the development of parameter-robust preconditioners and commutator arguments for solving Stokes control problems , 2015 .

[21]  Owe Axelsson,et al.  Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems , 2016, Numerical Algorithms.

[22]  O. Axelsson,et al.  Real valued iterative methods for solving complex symmetric linear systems , 2000 .

[23]  Maya Neytcheva,et al.  Efficient Preconditioners for Large Scale Binary Cahn-Hilliard Models , 2012, Comput. Methods Appl. Math..

[24]  Walter Zulehner,et al.  Nonstandard Norms and Robust Estimates for Saddle Point Problems , 2011, SIAM J. Matrix Anal. Appl..

[25]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

[26]  X. Wu,et al.  Numerical and computational efficiency of solvers for two-phase problems , 2013, Comput. Math. Appl..

[27]  Ralf Hiptmair,et al.  Operator Preconditioning , 2006, Comput. Math. Appl..

[28]  F. Tröltzsch.  Optimal Control of Partial Differential Equations: Theory, Methods and Applications , 2010 .

[29]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

[30]  Kent-André Mardal,et al.  Uniform preconditioners for the time dependent Stokes problem , 2006, Numerische Mathematik.

[31]  Andrew J. Wathen,et al.  A new approximation of the Schur complement in preconditioners for PDE‐constrained optimization , 2012, Numer. Linear Algebra Appl..

[32]  Alfio Borzì,et al.  Computational Optimization of Systems Governed by Partial Differential Equations , 2012, Computational science and engineering.

[33]  Michael A. Heroux,et al.  Trilinos users guide. , 2003 .