Abstract:We present a scalable approach to solve a class of elliptic partial differential equation (PDE)-constrained optimization problems with bound constraints. This approach utilizes a robust full-space interior-point (IP)-Gauss-Newton optimization method. To cope with the poorly-conditioned IP-Gauss-Newton saddle-point linear systems that need to be solved, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE-constrained optimization problems. A block Gauss-Seidel preconditioner is proposed for the GMRES-based solution of the IP-Gauss-Newton linear systems. It is shown, for a large-class of PDE- and bound-constrained optimization problems, that the spectrum of the block Gauss-Seidel preconditioned IP-Gauss-Newton matrix is asymptotically independent of discretization and is not impacted by the ill-conditioning that notoriously plagues interior-point methods. We propose a regularization and log-barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)-based solution of the related IP-Gauss-Newton-Schur complement linear systems. The scalability of the approach is demonstrated on an example problem with bound and nonlinear elliptic PDE constraints. The numerical solution of the optimization problem is shown to require a discretization independent number of IP-Gauss-Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill-conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of preconditioner and linear system matrix applies, achieved with algebraic multigrid based solvers, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE- and bound-constrained problems.

Approximate inner solvers for block preconditioning of the incompressible Navier–Stokes problems discretized by isogeometric analysis

Rotated block diagonal preconditioners for Navier-Stokes control problems

Comparative study of inner-outer Krylov solvers for linear systems in structured and high-order unstructured CFD problems

An augmented Lagrangian preconditioner for the control of the Navier--Stokes equations

Performance of preconditioned iterative linear solvers for cardiovascular simulations in rigid and deformable vessels

Effective Block Preconditioners for Fluid Dynamics Coupled to Reduced Models of a Non-Local Nature

Preconditioning of the generalized Stokes problem arising from the approximation of the time-dependent Navier-Stokes equations

Preconditioned iterative solvers for constrained high-order implicit shock tracking methods

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

A monolithic ALE Newton-Krylov solver with Multigrid-Richardson-Schwarz preconditioning for incompressible Fluid Structure Interaction

Novel Approach for solving the discrete Stokes problems based on Augmented Lagrangian and Global Techniques: Application to Saddle-Point Linear Systems from Incompressible flow

A robust and time-parallel preconditioner for parabolic reconstruction problems using Isogeometric Analysis

A preconditioning for the spectral solution of incompressible variable-density flows

A novel block non-symmetric preconditioner for mixed-hybrid finite-element-based flow simulations

Convergence analysis of iterative solution with inexact block preconditioning for weak Galerkin finite element approximation of Stokes flow

Preconditioning a mass‐conserving discontinuous Galerkin discretization of the Stokes equations

A Scalable Interior-Point Gauss-Newton Method for PDE-Constrained Optimization with Bound Constraints

A Symmetric Multigrid-Preconditioned Krylov Subspace Solver for Stokes Equations

Hybrid discontinuous Galerkin discretisation and domain decomposition preconditioners for the Stokes problem

A robust and efficient iterative method for hyper-elastodynamics with nested block preconditioning

Inexact block LU preconditioners for incompressible fluids with flow rate conditions