Efficient Time Domain Decomposition Algorithms for Parabolic PDE-Constrained Optimization Problems

Research paper by Jun Liu, Zhu Wang

Indexed on: 29 Dec '16Published on: 29 Dec '16Published in: arXiv - Mathematics - Numerical Analysis


Optimization with time-dependent partial differential equations (PDEs) as constraints appear in many engineering applications. The associated first-order necessary optimality system consists of one forward and one backward time-dependent PDE coupled with optimality conditions. An optimization process by using the one-shot method determines the optimal control, state and adjoint state at once, while, with the cost of solving a large scale, fully discrete optimality system. Hence, such one-shot method could easily become prohibitive when the time span is long or a small time step is taken. To overcome this difficulty, we propose in this paper several time domain decomposition algorithms for improving its computational efficiency. In these algorithms, the optimality system is split into many small subsystems over a much smaller time interval, which are coupled by appropriate continuity matching conditions. Both one-level and two-level multiplicative and additive Schwarz algorithms are developed for iteratively solving the decomposed subsystems in parallel. In particular, the convergence of the one-level multiplicative and additive Schwarz algorithms without overlap are proved. The effectiveness of our proposed algorithms is demonstrated by both 1D and 2D numerical experiments, where the developed two-level algorithms show very scalable convergence rates with respect to the number of subdomains.