Software for simplified Lanczos and QMR algorithms
Applied Numerical Mathematics - Special issue on iterative methods for linear equations
Domain decomposition: parallel multilevel methods for elliptic partial differential equations
Domain decomposition: parallel multilevel methods for elliptic partial differential equations
Computational Modeling in Semiconductor Processing
Computational Modeling in Semiconductor Processing
Perspectives in Flow Control and Optimization
Perspectives in Flow Control and Optimization
A variational finite element method for source inversion for convective-diffusive transport
Finite Elements in Analysis and Design - Special issue: 14th Robert J. Melosh competition
Journal of Computational and Applied Mathematics
SIAM Journal on Scientific Computing
Numerical Approximation of Partial Differential Equations
Numerical Approximation of Partial Differential Equations
Hi-index | 7.29 |
We present a non-overlapping spatial domain decomposition method for the solution of linear-quadratic parabolic optimal control problems. The spatial domain is decomposed into non-overlapping subdomains. The original parabolic optimal control problem is decomposed into smaller problems posed on space-time cylinder subdomains with auxiliary state and adjoint variables imposed as Dirichlet boundary conditions on the space-time interface boundary. The subdomain problems are coupled through Robin transmission conditions. This leads to a Schur complement equation in which the unknowns are the auxiliary state adjoint variables on the space-time interface boundary. The Schur complement operator is the sum of space-time subdomain Schur complement operators. The application of these subdomain Schur complement operators is equivalent to the solution of an subdomain parabolic optimal control problem. The subdomain Schur complement operators are shown to be invertible and the application of their inverses is equivalent to the solution of a related subdomain parabolic optimal control problem. We introduce a new family of Neumann-Neumann type preconditioners for the Schur complement system including several different coarse grid corrections. We compare the numerical performance of our preconditioners with an alternative approach recently introduced by Benamou.