An inexact l1 penalty SQP algorithm for PDE-constrained optimization with an application to shape optimization in linear elasticity. Hess, W. & Ulbrich, S. Optimization Methods and Software, 2012.
doi  abstract   bibtex   
We develop an optimization algorithm which is able to deal with inexact evaluations of the objective function. The proposed algorithm employs sequential quadratic programming with a line search that uses the ℓ1 penalty function for an Armijo-like condition. Both the objective gradient computations for the quadratic subproblems and the objective function computations for the line search admit some inexactness which is controlled by the algorithm. We prove convergence results for the presented algorithm under suitable assumptions. The kind of optimization problems handled by the algorithm arises, for example, when we apply iterative solvers in the evaluation of the reduced objective function of a shape optimization problem subject to the linear elasticity equations. In the second part of this paper, we present an application of the algorithm to the shape optimization of steel profiles using nested iterations with adaptively refined meshes. It is shown how the inexactness in the algorithm translates to an acceptable residual of the state and adjoint state solutions in the iterative solver. We look at three different stopping criteria for the outer iterations. Numerical results are presented.
@Article{         Hess_2012aa,
  abstract      = { We develop an optimization algorithm which is able to deal with inexact evaluations of the objective function. The proposed algorithm employs sequential quadratic programming with a line search that uses the ℓ1 penalty function for an Armijo-like condition. Both the objective gradient computations for the quadratic subproblems and the objective function computations for the line search admit some inexactness which is controlled by the algorithm. We prove convergence results for the presented algorithm under suitable assumptions. The kind of optimization problems handled by the algorithm arises, for example, when we apply iterative solvers in the evaluation of the reduced objective function of a shape optimization problem subject to the linear elasticity equations. In the second part of this paper, we present an application of the algorithm to the shape optimization of steel profiles using nested iterations with adaptively refined meshes. It is shown how the inexactness in the algorithm translates to an acceptable residual of the state and adjoint state solutions in the iterative solver. We look at three different stopping criteria for the outer iterations. Numerical results are presented. },
  author        = {Hess, Wolfgang and Ulbrich, Stefan},
  doi           = {10.1080/10556788.2011.651082},
  issn          = {1055-6788},
  journal       = {Optimization Methods and Software},
  langid        = {english},
  pages         = {1--26},
  title         = {An inexact l1 penalty {SQP} algorithm for {PDE}-constrained optimization with an application to shape optimization in linear elasticity},
  year          = {2012},
  shortjournal  = {Optim Meth Software}
}

Downloads: 0