MUSCOP swMATH ID: 6143 Software Authors: Potschka, Andreas Description: A direct method for the numerical solution of optimization problems with time-periodic PDE constraints. In this thesis we develop a numerical method based on direct multiple shooting for optimal control problems (OCPs) constrained by time-periodic partial differential equations (PDEs). The proposed method features asymptotically optimal scale-up of the numerical effort with the number of spatial discretization points. It consists of a linear iterative splitting approach (LISA) within a Newton-type iteration with globalization on the basis of natural level functions. We investigate the LISA-Newton method in the framework of Bock’s kappa-theory and develop reliable a-posteriori kappa-estimators. Moreover we extend the inexact Newton method to an inexact sequential quadratic programming (SQP) method for inequality constrained problems and provide local convergence theory. In addition we develop a classical and a two-grid Newton-Picard preconditioner for LISA and prove grid independent convergence of the classical variant for a model problem. Based on numerical results we can claim that the two-grid version is even more efficient than the classical version for typical application problems. Moreover we develop a two-grid approximation for the Lagrangian Hessian which fits well in the two-grid Newton-Picard framework and yields a reduction of 68 For the solution of the occurring large-scale quadratic programming problems (QPs) we develop a structure exploiting two-stage approach. In the first stage we exploit the multiple shooting and Newton-Picard structure to reduce the large-scale QP to an equivalent QP whose size is independent of the number of spatial discretization points. For the second stage we develop extensions for a parametric active set method (PASM) to achieve a reliable and efficient solver for the resulting, possibly nonconvex QP. Furthermore we construct three illustrative, counter-intuitive toy examples which show that convergence of a one-shot one-step optimization method is neither necessary nor sufficient for the convergence of the forward problem method. For three regularization approaches to recover convergence our analysis shows that de-facto loss of convergence cannot be avoided with these approaches. We have further implemented the proposed methods within a code called MUSCOP which features automatic derivative generation for the model functions and dynamic system solutions of first and second order, parallelization on the multiple shooting structure, and a hybrid language programming paradigm to minimize setup and solution time for new application problems. We demonstrate the applicability, reliability, and efficiency of MUSCOP and thus the proposed numerical methods and techniques on a sequence of PDE OCPs of growing difficulty ranging from linear academic problems, over highly nonlinear academic problems of mathematical biology to a highly nonlinear real-world chemical engineering problem in preparative chromatography: The simulated moving bed process. Homepage: http://www.iwr.uni-heidelberg.de/~Andreas.Potschka/software.html Keywords: direct multiple shooting; optimal control; linear iterative splitting approach; Newton-type iteration; sequential quadratic programming; convergence; Newton-Picard preconditioner; numerical results; grid refinement; parametric active set method; automatic derivative generation; parallelization Related Software: NewtonLib; ADOL-C; L-BFGS-B; LBFGS-B; L-BFGS; UFL; FEniCS; DOLFIN; MUMPS; PETSc; Chebfun; deal.ii; DAESOL-II; MINRES; CUTEst; MinRes; Trilinos; LiftOpt; SolvIND; NETLIB LP Test Set Cited in: 8 Documents Standard Articles 1 Publication describing the Software, including 1 Publication in zbMATH Year A direct method for the numerical solution of optimization problems with time-periodic PDE constraints. Zbl 1237.65062Potschka, Andreas 2001 all top 5 Cited by 7 Authors 7 Potschka, Andreas 1 Bock, Hans Georg 1 Fang, Liang 1 Hante, Falk M. 1 Meyers, Johan 1 Mommer, Mario S. 1 Vandewalle, Stefan G. Cited in 4 Serials 2 SIAM Journal on Numerical Analysis 1 Journal of Computational Physics 1 Numerical Algorithms 1 Mathematical Programming. Series A. Series B all top 5 Cited in 8 Fields 8 Numerical analysis (65-XX) 5 Calculus of variations and optimal control; optimization (49-XX) 5 Operations research, mathematical programming (90-XX) 3 Global analysis, analysis on manifolds (58-XX) 2 Partial differential equations (35-XX) 1 Operator theory (47-XX) 1 Mechanics of deformable solids (74-XX) 1 Biology and other natural sciences (92-XX) Citations by Year