Preconditioned SQP Solvers for Nonlinear Optimization Problems with PDEs


Project Description

Partial differential equations (PDEs) are the language to describe countless phenomena in the natural sciences and engineering. The optimization of such processes, or the identification of unknown model parameters, lead to optimization problems with PDEs. Sequential quadratic programming (SQP) algorithms are powerful and widely used solution methods for nonlinear problems of this type.

The overall effectiveness of an SQP algorithm depends on its global and local convergence properties, as well as on the fast solution of the quadratic programming (QP) subproblem in every iteration. In the proposed project, we will investigate preconditioners which are especially well suited for the efficient solution of the subproblems occuring in the popular composite-step trust-region SQP methods. Our research will lead to new preconditioned matrix-free SQP solvers for nonlinear large-scale optimization problems. Challenging applications governed by nonlinear, coupled and time-dependent PDEs will demonstrate the potential and limitations of these methods.

Associated Publications

  • Preconditioned solution of state gradient constrained elliptic optimal control problems
    SIAM Journal on Numerical Analysis 54(2), p.688-718, 2016
    bibtex
    @ARTICLE{HerzogMach:2016:1,
      AUTHOR = {Herzog, Roland and Mach, Susann},
      DATE = {2016},
      DOI = {10.1137/130948045},
      JOURNALTITLE = {SIAM Journal on Numerical Analysis},
      NUMBER = {2},
      PAGES = {688--718},
      TITLE = {Preconditioned solution of state gradient constrained elliptic optimal control problems},
      VOLUME = {54},
    }
  • Roland Herzog and Ekkehard Sachs
    Superlinear convergence of Krylov subspace methods for self-adjoint problems in Hilbert space
    SIAM Journal on Numerical Analysis 53(3), p.1304-1324, 2015
    bibtex
    @ARTICLE{HerzogSachs:2015:1,
      AUTHOR = {Herzog, Roland and Sachs, Ekkehard},
      DATE = {2015},
      DOI = {10.1137/140973050},
      JOURNALTITLE = {SIAM Journal on Numerical Analysis},
      NUMBER = {3},
      PAGES = {1304--1324},
      TITLE = {Superlinear convergence of Krylov subspace methods for self-adjoint problems in Hilbert space},
      VOLUME = {53},
    }

Agency logo