×

Distributed solution of optimal control problems governed by parabolic equations. (English) Zbl 1120.49024

Kurdila, Andrew J. (ed.) et al., Robust optimization-directed design. Papers based on the 1st Conference on Robust Optimization-Directed Design (RODD), Shalimar, FL, USA, April 19–21, 2004. New York, NY: Springer (ISBN 0-387-28263-7/hbk). Nonconvex Optimization and Its Applications 81, 73-94 (2006).
Summary: We present a spatial domain decomposition (DD) method for the solution of discretized parabolic linear-quadratic optimal control problems. Our DD preconditioners are extensions of Neumann-Neumann DD methods, which have been successfully applied to the solution of single elliptic partial differential equations and of linear-quadratic optimal control problems governed by elliptic equations.
We use a decomposition of the spatial domain into non-overlapping subdomains. The optimality conditions for the parabolic linear-quadratic optimal control problem are split into smaller problems restricted to spatial subdomain-time cylinders. These subproblems correspond to parabolic linear-quadratic optimal control problems on subdomains with Dirichlet data on interfaces. The coupling of these subdomain problems leads to a Schur complement system in which the unknowns are the state and adjoint variables on the subdomain interfaces in space and time.
The Schur complement system is solved using a preconditioned Krylov subspace method. The preconditioner is obtained from the solution of appropriate subdomain parabolic linear-quadratic optimal control problems. The dependence of the performance of these preconditioners on mesh size and subdomain size is studied numerically. Our tests indicate that their dependence on mesh size and subdomain size is similar to that of its counterpart applied to elliptic equations only. Our tests also suggest that the preconditioners are insensitive to the size of the control regularization parameter.
For the entire collection see [Zbl 1087.93006].

MSC:

49M25 Discrete approximations in optimal control
49M20 Numerical methods of relaxation type
65M55 Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs
49M27 Decomposition methods

Software:

QMRPACK
PDFBibTeX XMLCite
Full Text: DOI