Optimal solvers for PDE-Constrained Optimization

Optimization problems with constraints which require the solution of a partial differential equation arise widely in many areas of the sciences and engineering, in particular in problems of design. The solution of such PDE-constrained optimization problems is usually a major computational task. Here...

Description complète

Détails bibliographiques
Auteurs principaux: Rees, T, Dollar, H, Wathen, A
Format: Report
Publié: Unspecified 2008
_version_ 1826267794775212032
author Rees, T
Dollar, H
Wathen, A
author_facet Rees, T
Dollar, H
Wathen, A
author_sort Rees, T
collection OXFORD
description Optimization problems with constraints which require the solution of a partial differential equation arise widely in many areas of the sciences and engineering, in particular in problems of design. The solution of such PDE-constrained optimization problems is usually a major computational task. Here we consider simple problems of this type: distributed control problems in which the 2- and 3-dimensional Poisson problem is the PDE. The large dimensional linear systems which result from discretization and which need to be solved are of saddle-point type. We introduce two optimal preconditioners for these systems which lead to convergence of symmetric Krylov subspace iterative methods in a number of iterations which does not increase with the dimension of the discrete problem. These preconditioners are block structured and involve standard multigrid cycles. The optimality of the preconditioned iterative solver is proved theoretically and verified computationally in several test cases. The theoretical proof indicates that these approaches may have much broader applicability for other partial differential equations.
first_indexed 2024-03-06T20:59:41Z
format Report
id oxford-uuid:3a772dfd-c1f4-4052-acd0-77b25f94b0e3
institution University of Oxford
last_indexed 2024-03-06T20:59:41Z
publishDate 2008
publisher Unspecified
record_format dspace
spelling oxford-uuid:3a772dfd-c1f4-4052-acd0-77b25f94b0e32022-03-26T14:01:50ZOptimal solvers for PDE-Constrained OptimizationReporthttp://purl.org/coar/resource_type/c_93fcuuid:3a772dfd-c1f4-4052-acd0-77b25f94b0e3Mathematical Institute - ePrintsUnspecified2008Rees, TDollar, HWathen, AOptimization problems with constraints which require the solution of a partial differential equation arise widely in many areas of the sciences and engineering, in particular in problems of design. The solution of such PDE-constrained optimization problems is usually a major computational task. Here we consider simple problems of this type: distributed control problems in which the 2- and 3-dimensional Poisson problem is the PDE. The large dimensional linear systems which result from discretization and which need to be solved are of saddle-point type. We introduce two optimal preconditioners for these systems which lead to convergence of symmetric Krylov subspace iterative methods in a number of iterations which does not increase with the dimension of the discrete problem. These preconditioners are block structured and involve standard multigrid cycles. The optimality of the preconditioned iterative solver is proved theoretically and verified computationally in several test cases. The theoretical proof indicates that these approaches may have much broader applicability for other partial differential equations.
spellingShingle Rees, T
Dollar, H
Wathen, A
Optimal solvers for PDE-Constrained Optimization
title Optimal solvers for PDE-Constrained Optimization
title_full Optimal solvers for PDE-Constrained Optimization
title_fullStr Optimal solvers for PDE-Constrained Optimization
title_full_unstemmed Optimal solvers for PDE-Constrained Optimization
title_short Optimal solvers for PDE-Constrained Optimization
title_sort optimal solvers for pde constrained optimization
work_keys_str_mv AT reest optimalsolversforpdeconstrainedoptimization
AT dollarh optimalsolversforpdeconstrainedoptimization
AT wathena optimalsolversforpdeconstrainedoptimization