Abstract
A mesh-independent and second-order accurate multigrid strategy to solve control-constrained parabolic optimal control problems is presented. The resulting algorithms appear to be robust with respect to change of values of the control parameters and have the ability to accommodate constraints on the control also in the limit case of bang-bang control. Central to the development of these multigrid schemes is the design of iterative smoothers which can be formulated as local semismooth Newton methods. The design of distributed controls is considered to drive nonlinear parabolic models to follow optimally a given trajectory or attain a final configuration. In both cases, results of numerical experiments and theoretical twogrid local Fourier analysis estimates demonstrate that the proposed schemes are able to solve parabolic optimality systems with textbook multigrid efficiency. Further results are presented to validate second-order accuracy and the possibility to track a trajectory over long time intervals by means of a receding-horizon approach.
Similar content being viewed by others
References
Allgöwer, F., Badgwell, T.A., Qin, J.S., Rawlings, J.B., Wright, S.J.: Nonlinear predictive control and moving horizon estimation—An introduction overview. In: Frank, P.M. (ed.) Advances in Control, Highlights of ECC’99, pp. 391–449 (1999), Chap. 5
Ascher, U.M., Petzold, L.R.: Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations. SIAM, Philadelphia (1998)
Batista, M.: A method for solving cyclic block penta-diagonal systems of linear equations (2008). arXiv:0803.0874v3
Batista, M., Karawia, A.A.: The use of the Sherman-Morrison-Woodbury formula to solve cyclic block tri-diagonal and cyclic block penta-diagonal linear systems of equations. Appl. Math. Comput. 210, 558–563 (2009)
Benkert, K., Fischer, R.: An efficient implementation of the Thomas-Algorithm for block penta-diagonal systems on vector computers. In: Shi, Y., van Albada, G.D., Dongarra, J., Sloot, P.M.A. (eds.) Proceedings of the 7th International Conference on Computer Science, ICCS (2007)
Biegler, L.T., Ghattas, O., Heinkenschloss, M., Keyes, D., van Bloemen Waanders, B. (eds.) Real-Time PDE-Constrained Optimization, Computational Science and Engineering, vol. 3. SIAM, Philadelphia (2007)
Borzì, A.: High-order discretization and multigrid solution of elliptic nonlinear constrained optimal control problems. J. Comput. Appl. Math. 200, 67–85 (2007)
Borzì, A.: Multigrid methods for parabolic distributed optimal control problems. J. Comput. Appl. Math. 157, 365–382 (2003)
Borzì, A.: Space-time multigrid methods for solving unsteady optimal control problems. In: Biegler, L.T., Ghattas, O., Heinkenschloss, M., Keyes, D., van Bloemen Waanders, B. (eds.) Real-Time PDE-Constrained Optimization, Computational Science and Engineering, vol. 3. SIAM, Philadelphia (2007), Chap. 5
Borzì, A., Griesse, R.: Distributed optimal control of lambda-omega systems. J. Numer. Math. 14, 17–40 (2006)
Borzì, A., Kunisch, K.: A multigrid scheme for elliptic constrained optimal control problems. Comput. Optim. Appl. 31, 309–333 (2005)
Borzì, A., Schulz, V.: Multigrid methods for PDE optimization. SIAM Rev. 51, 361–395 (2009)
Borzì, A., von Winckel, G.: Multigrid methods and sparse-grid collocation techniques for parabolic optimal control problems with random coefficients. SIAM J. Sci. Comput. 31, 2172–2192 (2009)
Brandt, A.: Multi-level adaptive solutions to boundary-value problems. Math. Comput. 31, 333–390 (1977)
Brézis, H., Crandall, M.-G., Pazy, A.: Perturbations of nonlinear maximal monotone sets in Banach spaces. Commun. Pure Appl. Math. 23, 123–144 (1970)
Dreyer, Th., Maar, B., Schulz, V.: Multigrid optimization in applications. J. Comput. Appl. Math. 120, 67–84 (2000)
Emmrich, E.: Two-step BDF time discretisation of nonlinear evolution problems governed by monotone operators with strongly continuous perturbations. Comput. Methods Appl. Math. 9, 37–62 (2009)
Glashoff, K., Sachs, E.: On theoretical and numerical aspects of the bang-bang principle. Numer. Math. 29, 93–113 (1977)
Goldberg, H., Tröltzsch, F.: Second order sufficient optimality conditions for a class of non-linear parabolic boundary control problems. SIAM J. Control Optim. 31, 1007–1027 (1993)
Goldberg, H., Tröltzsch, F.: On a SQP–multigrid technique for nonlinear parabolic boundary control problems. In: Hager, W.W., Pardalos, P.M. (eds.) Optimal Control: Theory, Algorithms, and Applications, pp. 154–174. Kluwer Academic, Dordrecht (1998)
Hackbusch, W.: Parabolic multigrid methods. In: Glowinski, R., Lions, J.-L. (eds.) Computing Methods in Applied Sciences and Engineering VI. North-Holland, Amsterdam (1984)
Hackbusch, W.: Elliptic Differential Equations. Springer, New York (1992)
Hackbusch, W.: On the fast solving of parabolic boundary control problems. SIAM J. Control Optim. 17, 231–244 (1979)
Hackbusch, W.: Numerical Solution of Linear and Nonlinear Parabolic Optimal Control Problems. Lecture Notes in Control and Information Science, vol. 30. Springer, Berlin (1981)
Hintermüller, M., Ito, K., Kunisch, K.: The primal-dual active set strategy as a semi-smooth Newton method. SIAM J. Optim. 13, 865–888 (2003)
Hintermüller, M., Hoppe, R.H.W.: Goal-oriented adaptivity in control constrained optimal control of partial differential equations. SIAM J. Control Optim. 47, 1721–1743 (2008)
Hinze, M.: A variational discretization concept in control constrained optimization: the linear-quadratic case. Comput. Optim. Appl. 30, 45–63 (2005)
Horton, G., Vandewalle, S.: A space-time multigrid method for parabolic partial differential equations. SIAM J. Sci. Comput. 16(4), 848–864 (1995)
Ito, K., Kunisch, K.: Asymptotic properties of receding horizon optimal control problems. SIAM J. Control Optim. 40(5), 1585–1610 (2002)
Kunisch, K., Volkwein, S.: Proper orthogonal decomposition for optimality systems. ESAIM: Math. Model. Numer. Anal. 42, 1–23 (2008)
Lions, J.L.: Optimal Control of Systems Governed by Partial Differential Equations. Springer, Berlin (1971)
Malanowski, K.: Convergence of approximations vs. regularity of solutions for convex, control-constrained optimal-control problems. Appl. Math. Optim. 8, 69–95 (1981)
Meidner, D., Vexler, B.: Adaptive space-time finite element methods for parabolic optimization problems. SIAM J. Control Optim. 46, 116–142 (2007)
Nagaiah, Ch., Kunisch, K., Plank, G.: Numerical solution for optimal control of the reaction-diffusion equations in cardiac electrophysiology. Comput. Optim. Appl. doi:10.1007/s10589-009-9280-3
Neittaanmäki, P., Tiba, D.: Optimal Control of Nonlinear Parabolic Systems. Dekker, New York (1994)
Qi, L., Sun, J.: A nonsmooth version of Newton’s method. Math. Program. 58, 353–368 (1993)
Rösch, A.: Error estimates for linear-quadratic control problems with control constraints. Optim. Methods Softw. 21, 121–134 (2006)
Showalter, R.E.: Monotone Operators in Banach Space and Nonlinear Partial Differential Equations. Mathematical Surveys and Monographs AMS. AMS, Providence (1997)
Stadler, G.: Semismooth Newton and augmented Lagrangian methods for a simplified friction problem. SIAM J. Optim. 15, 39–62 (2004)
Thomas, J.W.: Numerical Partial Differential Equations: Finite Difference Methods. Springer, Berlin (1995)
Trottenberg, U., Oosterlee, C., Schüller, A.: Multigrid. Academic Press, London (2001)
Ulbrich, M.: Semismooth Newton methods for operator equations in function spaces. SIAM J. Optim. 13, 805–842 (2002)
Varga, S.R.: Matrix Iterative Analysis. Prentice Hall, New York (1962)
Zuazua, E.: Switching control. J. Eur. Math. Soc. (to appear)
Author information
Authors and Affiliations
Corresponding author
Additional information
Supported in part by the Austrian Science Fund SFB Project F3205-N18 “Fast Multigrid Methods for Inverse Problems”.
Rights and permissions
About this article
Cite this article
González Andrade, S., Borzì, A. Multigrid second-order accurate solution of parabolic control-constrained problems. Comput Optim Appl 51, 835–866 (2012). https://doi.org/10.1007/s10589-010-9358-y
Received:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10589-010-9358-y