An Adaptive Method for Optimal Control Problems Constrained by Parabolic Differential Equations

TL;DR

This paper introduces an adaptive collocation method combining hp-Galerkin finite elements and orthogonal collocation for efficiently solving optimal control problems constrained by parabolic PDEs, with integrated error estimation and mesh adaptivity.

Contribution

It develops a novel adaptive method that integrates variational reformulation, integral linearization, residual-based error estimation, and solution regularity assessment for parabolic optimal control.

Findings

Error reduction up to five orders of magnitude.

Effective adaptive mesh refinement based on solution regularity.

Demonstrated efficiency on numerical examples.

Abstract

An adaptive direct collocation method is developed for solving optimal control problems constrained by parabolic partial differential equations. The partial differential equation is first reformulated in a variational setting, where the spatial domain is discretized using the hp-Galerkin finite element method. To address nonlinearities in the variational form, a Kirchhoff-like integral transformation is applied to linearize the dynamics. In the temporal dimension, an orthogonal collocation scheme, the hp-flipped Legendre-Gauss-Radau method, is employed to fully discretize the problem, yielding a large, sparse nonlinear programming problem. Upon solving the nonlinear programming problem, solution accuracy is assessed through an implicit residual estimation procedure. This approach evaluates the local error by solving auxiliary residual problems over selected subdomains, providing a novel means of error estimation within an orthogonal collocation framework for optimal control. Based on the computed error estimate, the mesh is adaptively refined or coarsened to meet a prescribed error tolerance. Mesh refinement is guided by the estimated regularity of the solution which is determined via the decay rate of the coefficients of a Legendre polynomial expansion. In overcollocated regions, a mesh reduction strategy is adapted from orthogonal collocation methods for application within the finite element framework. Numerical examples demonstrate that the proposed method can reduce the error by up to five orders of magnitude in both spatial and temporal dimensions.