Abstract
We are concerned with efficient solutions of the time-periodic parabolic optimal control problems. By using the multiharmonic FEM, the linear algebraic equations characterizing the first-order optimality conditions can be decoupled into a series of parallel solvable block 4 × 4 linear systems with respect to the cosine and sine Fourier coefficients of the state and scaled control variables for different frequencies. Parameter robust preconditioners are proposed for solving these linear systems along with information on practical algorithm implementation and detailed spectral analysis. Problem independent eigenvalue bounds and upper bound approximations of the condition numbers of the eigenvector matrices are obtained for the preconditioned matrices. Such results ensure efficient Krylov subspace acceleration methods and a parameter-free Chebyshev acceleration method, which are both robust in view of all discretization and model parameters. Numerical experiments are presented to demonstrate the robustness and effectiveness of the proposed preconditioners within both Krylov subspace and Chebyshev accelerations compared with some already available preconditioned Krylov subspace methods.
Similar content being viewed by others
References
Abbeloos, D, Diehl, M., Hinze, M., Vandewalle, S.: Nested multigrid methods for time-periodic, parabolic optimal control problems. Comput. Vis. Sci. 14, 27–38 (2011)
Anjam, I., Valdman, J.: Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements. Appl. Math. Comput. 267, 252–263 (2015)
Axelsson, O.: Iterative Solution Methods. Cambridge University Press, Cambridge (1994)
Axelsson, O., Farouq, S., Neytcheva, M.: Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems Poisson and convection–diffusion control. Numer. Algorithms 73(3), 631–663 (2016)
Axelsson, O., Liang, Z.-Z.: A note on preconditioning methods for time-periodic eddy current optimal control problems. J. Comput. Appl. Math. 352, 262–277 (2019)
Axelsson, O., Lukáš, D.: Preconditioning methods for eddy-current optimally controlled time-harmonic electromagnetic problems. J. Numer. Math. 27(1), 1–21 (2019)
Axelsson, O., Neytcheva, M., Ahmad, B.: A comparison of iterative methods to solve complex valued linear algebraic systems. Numer. Algorithm. 66 (4), 811–841 (2014)
Bai, Z. -Z.: Block preconditioners for elliptic PDE-constrained optimization problems. Computing 91(4), 379–395 (2011)
Bai, Z. -Z., Benzi, M., Chen, F., Wang, Z.-Q.: Preconditioned MHSS Iteration methods for a class of block two-by-two linear systems with applications to distributed control problems. IMA. J. Numer. Anal. 33(1), 343–369 (2013)
Bai, Z. -Z., Chen, F., Wang, Z.-Q.: Additive block diagonal preconditioning for block two-by-two linear systems of skew-H,amiltonian coefficient matrices. Numer. Algorithm. 62, 655–675 (2013)
Barker, A.T., Stoll, M.: Domain decomposition in time for PDE,-constrained optimization. Comput. Phys Commun. 197, 136–143 (2015)
Benamou, J.D.: Domain decomposition, optimal control of systems governed by partial differential equations, and synthesis of feedback laws. J. Optim. Theory Appl. 102(1), 15–36 (1999)
Benamou, J.D., Desprès, B.: A domain decomposition method for the Helmholtz equation and related optimal control problems. J. Comput. Phy. 136(1), 68–82 (1997)
Benzi, M., Golub, G.H., Liesen, J.: Numerical solution of saddle point problems. Acta Numer. 14, 1–137 (2005)
Borzì, A., Schulz, V.: Multigrid methods for PDE optimization. SIAM Rev. 51(2), 361–395 (2009)
Bramble, J.H., Pasciak, J.E.: A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems. Math. Comp. 51, 1–17 (1988)
Elman, H.C., Silverster, D.J., Wathen, A.J.: Finite elements and fast iterative solvers: with applications in incompressible fluid dynamics. Oxford university press, New York (2005)
Elman, H.C., Ramage, A., Silvester, D.J.: Algorithm 866 IFISS, a Matlab toolbox for modelling incompressible flow. ACM Trans. Math. Softw. 33(2), 14 (2007)
Elman, H.C., Ramage, A., Silvester, D.J.: IFISS A computational laboratory for investigating incompressible flow problems. SIAM Rev. 56(2), 261–273 (2014)
Gould, N.I.M., Hribar, M.E., Nocedal, J.: On the solution of equality constrained quadratic programming problems arising in optimization. SIAM J. Sci. Comput. 23(4), 1376–1395 (2001)
Greenbaum, A.: Iterative Methods for Solving Linear Systems, vol. 17. SIAM, Philadelphia (1997)
Güttel, S., Pearson, J.W.: A rational deferred correction approach to parabolic optimal control problems. IMA J. Numer. Anal. 38(4), 1861–1892 (2017)
Herzog, R., Sachs, E.: Preconditioned conjugate gradient method for optimal control problems with control and state constraints. SIAM. J. Matrix Anal. Appl. 31(5), 2291–2317 (2010)
Hinze, M., Pinnau, R., Ulbrich, M., Ulbrich, S.: Optimization with PDE Constraints, Mathematical Modelling: Theory and Applications, vol. 23. Springer, New York (2008)
Kollmann, M.: Efficient iterative solvers for saddle point systems arising in PDE-constrained optimization problems with inequality constraints. PhD thesis Linz: Johannes Kepler University (2013)
Kollmann, M., Kolmbauer, M.: A preconditioned minres solver for time-periodic parabolic optimal control problems. Numer Linear Algebra Appl. 20(5), 761–784 (2013)
Kollmann, M., Kolmbauer, M., Langer, U., Wolfmayr, M., Zulehner, W.: A robust finite element solver for a multiharmonic parabolic optimal control problem. Comput. Math. Appl. 65(3), 469–486 (2013)
Kolmbauer, M.: The multiharmonic finite element and boundary element method for simulation and control of eddy current problems. PhD thesis Linz: Johannes Kepler University (2012)
Krendl, W., Simoncini, V., Zulehner, W.: Stability estimates and structural spectral properties of saddle point problems. Numer. Math. 124(1), 183–213 (2013)
Langer, U., Wolfmayr, M.: Multiharmonic finite element analysis of a time-periodic parabolic optimal control problem. J. Numer. Math. 21 (4), 265–300 (2013)
Liang, Z. -Z., Axelsson, O., Neytcheva, M.: A robust structured preconditioner for time-harmonic parabolic optimal control problems. Numer. Algorithm. 79, 575–596 (2018)
Logist, F., Wouwer, A.V., Smets, I., Impe, J.V.: Optimal temperature profiles for tubular reactors implemented through a flow reversal strategy. Chem. Eng Process. 62, 4675–4688 (2007)
Mardal, K.A., Nielsen, B.F., Nordaas, M.: Robust preconditioners for PDE,-constrained optimization with limited observations. BIT Numer. Math. 57, 405–431 (2017)
Napov, A., Notay, Y.: An algebraic multigrid method with guaranteed convergence rate. SIAM J. Sci. Comput. 34(2), A1079–A1109 (2012)
Nédélec, J. -C.: Mixed finite elements in \(\mathbb {R}^{3}\). Numer. Math. 35(3), 315–341 (1980)
Nédélec, J. C.: A new family of mixed finite elements in \(\mathbb {R}^{3}\). Numer. Math. 50(1), 57–81 (1986)
Nilchan, S., Pantelides, C.: On the optimisation of periodic adsorption processes. Adsorption 4, 113–147 (1998)
Notay, Y.: AGMG software and documentation; see http://agmg.eu/
Notay, Y.: An aggregation-based algebraic multigrid method. Electron. Trans. Numer. Anal. 37, 123–146 (2010)
Notay, Y.: Aggregation-based algebraic multigrid for convection-diffusion equations. SIAM J. Sci. Comput. 34(4), A2288–A2316 (2012)
Paige, C.C., Saunders, M.A.: Solution of sparse indefinite systems of linear equations. SIAM J. Numer. Anal. 12(4), 617–629 (1975)
Pearson, J.W., Gondzio, J.: Fast interior point solution of quadratic programming problems arising from PDE-constrained optimization. Numer. Math. 137(4), 959–999 (2017)
Pearson, J.W., Stoll, M., Wathen, A.J.: Robust iterative solution of a class of time-dependent optimal control problems. PAMM 12(1), 3–6 (2012)
Pearson, J.W., Wathen, A.J.: Matching Schur complement approximations for certain saddle-point systems. In Contemporary Computational Mathematics-A Celebration of the 80th Birthday of Ian Sloan, pp. 1001–1016. Springer (2018)
Pearson, J.W., Stoll, M., Wathen, A.J.: Regularization-robust preconditioners for time-dependent PDE,-constrained optimization problems. SIAM J Matrix Anal. Appl. 33(4), 1126–1152 (2012)
Pearson, J.W., Wathen, A.J.: A new approximation of the S,chur complement in preconditioners for PDE-constrained optimization. Numer. Linear Algebra Appl. 19(5), 816–829 (2012)
Rees, T., Dollar, H.S., Wathen, A.J.: Optimal solvers for PDE,-constrained optimization. SIAM J. Sci. Comput. 32(1), 271–298 (2010)
Ress, T., Stoll, M.: Block-triangular preconditioners for PDE,-constrained optimization. Numer. Linear Algebra Appl. 17(6), 977–996 (2010)
Saad, Y.: A flexible inner-outer preconditioned GMRES Algorithm. SIAM J. Sci. Comput. 14(2), 461–419 (1993)
Saad, Y.: Iterative methods for sparse linear systems. SIAM (2003)
Saad, Y., Schultz, M.h.: GMRES A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 7(3), 856–869 (1986)
Schöberl, J., Zulehner, W.: Symmetric indefinite preconditioners for saddle point problems with applications to PDE-constrained optimization problems. SIAM J. Matrix Anal. Appl. 29(3), 752–773 (2007)
Simoncini, V.: Reduced order solution of structured linear systems arising in certain PDE-constrained optimization problems. Comput. Optim. Appl. 53(2), 591–617 (2012)
Stoll, M.: One-shot solution of a time-dependent time-periodic PDE-constrained optimization problem. IMA J. Numer. Anal. 34(4), 1554–1577 (2013)
Stoll, M., Breiten, T.: A low-rank in time approach to PDE-constrained optimization. SIAM J. Sci. Comput. 37(1), B1–B29 (2015)
Stoll, M., Wathen, A.J.: All-at-once solution of time-dependent PDE-constrained optimization problems. Technical Report 1017. The Mathematical Institute, University of Oxford (2010)
Toumi, A., Engell, S., Diehl, M., Bock, H.G., Schloder, J.P.: Efficient optimization of simulated moving bed processes. Chem. Eng Process. 46, 1067–1084 (2007)
Tröltzsch, F.: Optimal Control of Partial Differential Equations: Theory, Methods and Applications, vol. 112. American Mathematical Society, Providence (2010)
Wolfmayr, M.: Multiharmonic finite element analysis of parabolic time-periodic simulation and optimal control problems. Phd thesis Linz: Johannes Kepler University (2014)
Zulehner, W.: Nonstandard norms and robust estimates for saddle point problems. SIAM J Matrix Anal. Appl. 32(2), 536–560 (2011)
Acknowledgements
We would like to express our sincere thanks to the two anonymous reviewers for their insightful comments and valuable suggestions, which greatly improve the quality of the paper. We would also thank one of them for pointing us to [12, 13, 15] and references therein, which are unnoticed by us previously, about the parallelizable domain decomposition and multigrid methods for solving the time-dependent PDE-constrained optimal control problems.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 11801242 and 11771193).
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by: Stefan Volkwein
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Liang, ZZ., Zhang, GF. Robust preconditioning techniques for multiharmonic finite element method with application to time-periodic parabolic optimal control problems. Adv Comput Math 47, 67 (2021). https://doi.org/10.1007/s10444-021-09887-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10444-021-09887-2
Keywords
- Parabolic optimal control problem
- Multiharmonic FEM
- Preconditioning
- Spectral analysis
- Krylov subspace acceleration
- Chebyshev acceleration