Funding agency | Austrian Science Fund (FWF) |
---|---|

Project number | J 3362 |

Principal investigator | Stefan Takacs |

First research stay | Andy Wathen Mathematical Institute University of Oxford Start: 2012-11-15End: 2013-11-16 |

Second research stay | Roland Herzog Fakultät für Mathematik TU Chemnitz Start: 2013-11-18End: 2014-11-26 |

Return phase | Ulrich Langer Institute of Computational Mathematics (NuMa) Universität Linz Start: 2014-12-01End: 2015-11-30 |

The main scope of the proposed project is the construction and the analysis of fast numerical methods for solving optimization problems constrained to partial differential equations (PDE-constrained optimization problems). To this class of problems belong optimal control problems, optimal design problems, shape and topology optimization problems, and many others. We will concentrate on optimal control problems. The proposer has considered in his previous work elliptic partial differential equations (PDEs) as constraints. The solution of such a problem is characterized by the optimality system, which is a system of PDEs. Such systems can be discretized using standard techniques. The resulting system is typically a large-scale linear system, which is sparse and indefinite. The construction of fast iterative solvers for such a system is of particular interest.

One issue in solving such linear systems is robustness. The optimal control problem itself depends on a parameter, which can be interpreted as regularization parameter or as cost parameter. If this parameter approaches zero, the condition number of the problem grows. Typically this leads to deteriorating convergence rates. A second parameter is introduced by discretizing the problem: the grid size or the grid level. Also here, if the grid is refined, the condition number of the problems grows. We are interested in an iterative method where the number of iterations is independent of these parameters. For standard elliptic problems, robustness in the grid parameter (optimal convergence order) can be achieved using hierarchical methods, like multigrid methods.

To obtain such a behavior also for an optimal control problem, several approaches are available. One possibility is the use of standard multigrid methods as a part of a block-preconditioner that is applied in the framework of the outer iteration, like a Krylov space method. An alternative is to apply the multigrid idea to the whole block-system (all-at-once approach or one-shot approach).

Both kinds of approaches have already been successfully applied to the elliptic optimal control problem mentioned above. The goal of the project is to develop such methods also for other PDE-constrained optimization problems, in particular for other optimal control problems. Such problems involve other state equations, like the Stokes control problem (velocity tracking), optimization of elastic deformation or optimization of Maxwell's equations. Also the consideration of box constraints for state or control, is of interest. The mentioned problems are linear but the constrained problems are non-linear due to the constraints. The methods introduced above have the potential to solve, both, the problems problems which are already linear and the linearizations of the non-linear problems. The linearizations may arise as subproblems to be solved in Newton-like methods. For these subproblems, robustness is an particular issue as additional parameters, like penalty parameters or active sets, appear.

- A robust multigrid method for Isogeometric Analysis in two dimensions using boundary correction.
- Computer Methods in Applied Mechanics and Engineering, vol. 316, p. 22 – 42, 2017.
- arXiv:1512.07091
- Approximation error estimates and inverse inequalities for B-splines of maximum smoothness.
- Mathematical Models and Methods in Applied Sciences, vol. 26 (7), p. 1411 – 1445, 2016.
- arXiv:1502.03733
- A robust all-at-once multigrid method for the Stokes control problem.
- Numerische Mathematik, vol. 130 (3), p. 517 – 540, 2015.
- arXiv:1502.04121
- A robust multigrid method for the time-dependent Stokes problem.
- SIAM Journal on Numerical Analysis, vol. 53 (6), p. 2634 – 2654, 2015.
- arXiv:1502.04070
- Using cylindrical algebraic decomposition and local Fourier analysis to study numerical methods: Two examples.
- In F. Winkler, V. Negru, T. Ida, T. Jebelean, D. Petcu, S. Watt, and D. Zaharie (eds.): 16th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing (SYNASC), p. 42 – 49, 2014.
- arXiv:1502.03926 Supplemental material
- Efficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems.
- In C. Pötzsche, C. Heuberger, B. Kaltenbacher, and F. Rendl (eds.): System Modeling and Optimization, p. 337 – 347, 2014.
- arXiv:1502.03917