General Order Virtual Element Methods for Neumann Boundary Optimal Control Problems in Saddle Point Formulation

Andrea Borio,Francesca Marcon,Maria Strazzullo

2024-11-13

Abstract:In this work, we explore the application of the Virtual Element Methods for Neumann boundary Optimal Control Problems in saddle point formulation. The method is proposed for arbitrarily polynomial order of accuracy and general polygonal meshes. Our contribution includes a rigorous a priori error estimate that holds for general polynomial degree. On the numerical side, we present (i) an initial convergence test that reflects our theoretical findings, and (ii) a second test case based on a more application-oriented experiment. For the latter test, we focus on the role of VEM stabilization, conducting a detailed experimental analysis, and proposing an alternative structure-preserving strategy to circumvent issues related to the choice of the stabilization parameter.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to apply the Virtual Element Methods (VEM) to solve Neumann boundary Optimal Control Problems (OCPs), especially the numerical approximation with arbitrary polynomial accuracy on general polygonal meshes under the saddle - point structure**.

Specifically, this paper focuses on the VEM discretization of the unconstrained Neumann boundary control problem under the saddle - point structure and analyzes its stability and error estimates. The main objectives include:

1. **Propose a VEM method applicable to any polynomial order** to deal with Neumann boundary optimal control problems.
2. **Prove the well - posedness and a priori error estimates of the proposed VEM method** to ensure the theoretical reliability of the method.
3. **Verify the theoretical results through numerical experiments**, including the initial convergence test and a more application - oriented test case.
4. **Study the influence of the choice of VEM stabilization parameters on the accuracy of the method** and propose a structure - preserving strategy to avoid problems related to the choice of stabilization parameters.

The main contributions of the paper are:
- A VEM method with arbitrary polynomial accuracy is proposed.
- Neumann boundary optimal control problems are analyzed under the saddle - point structure.
- Well - posedness and a priori error estimates for general polynomial orders and variable - coefficient cases are proved.
- Numerical experiments cover high - order schemes and the influence of stabilization parameters is studied in detail.
- A structure - preserving strategy based on the latest VEM theory development is proposed to avoid problems brought by the choice of stabilization parameters.

Through these works, the paper provides a new, flexible and efficient numerical solution method for Neumann boundary optimal control problems.