An Error Analysis of Second Order Elliptic Optimal Control Problem via Hybrid Higher Order Methods
Gouranga Mallik, Ramesh Chandra Sau
TL;DR
The paper develops Hybrid High-Order (HHO) methods for distributed optimal control of the Poisson equation, proposing three unconstrained schemes (UC1, UC2, UC3.1/UC3.2) and two constrained schemes (WC1, WC2). By leveraging a reconstruction-based approach, UC3.1/UC3.2 achieve improved control accuracy (up to \\(k+2\\) order) compared to standard FEM, while WC1 attains linear convergence with minimal element order and WC2 attains cubic convergence for the control under suitable regularity. The analysis combines the HHO framework with variational discretization concepts, providing rigorous energy and L2-error estimates for states, adjoints, and controls. Numerical experiments on rectangular and polygonal meshes validate the theoretical rates and demonstrate the practical robustness and efficiency of the proposed schemes. The results position HHO methods as powerful tools for high-accuracy optimal control in polytopal meshes, with clear advantages over traditional discretizations.
Abstract
This paper presents the design and analysis of a Hybrid High-Order (HHO) approximation for a distributed optimal control problem governed by the Poisson equation. We propose three distinct schemes to address unconstrained control problems and two schemes for constrained control problems. For the unconstrained control problem, while standard finite elements achieve a convergence rate of \( k+1 \) (with \( k \) representing the polynomial degree), our approach enhances this rate to \( k+2 \) by selecting the control from a carefully constructed reconstruction space. For the box-constrained problem, we demonstrate that using lowest-order elements (\( \mathbb{P}_0 \)) yields linear convergence, in contrast to finite element methods (FEM) that require linear elements to achieve comparable results. Furthermore, we derive a cubic convergence rate for control in the variational discretization scheme. Numerical experiments are provided to validate the theoretical findings.
