In this paper, we study the finite volume element method for solving the Neumann boundary control problems governed by elliptic partial differential equations. The following control problem will be considered
where is a plane polygonal domain with piecewise smooth boundary , is a 2 × 2 symmetric and uniformly positive matrix, is a sufficient smooth function defined on Ω, B denotes the linear and continuous control operator which makes , and have enough regularity so that this problem has a unique solution.
The control problem (1)-(2) has many applications, see for example , for a brief review of applications. Finite element method is an important numerical method for the problems of partial differential equations and is widely used in the numerical solution of distributed optimal control problems    . Compared with distributed control problems, the Neumann boundary control problem is more difficult. Error estimates for finite element discretization of Neumann boundary control problems governed by elliptic equations are discussed in some publications. In , piecewise constant functions are used to discretize the control and obtain the rates of convergence in . The continuous piecewise linear functions are used to discretize the control and the corresponding error estimates are obtained in . The variational discretization of Neumann type elliptic optimal control problems are considered by Hinze . Error estimates for the postprocessing approach applied to the Neumann boundary control problems in polyhedral domains are considered  .
The finite volume element (FVE) method has been one of the most commonly used numerical methods for solving partial differential equations. The main advantage of the method is that the local physical conservation law can be maintained. So it has been extensively used in computational fluid dynamics. We can refer to  for groundwater flow,  for weather prediction,  for shallow water wave, and to  for sedimentation problem. However, there are only a few published results on the finite volume element method for the distributed optimal control problems. In  , the authors discussed distributed optimal control problems governed by elliptic equations by using the finite volume element methods. The variational discretization approach is used to deal with the control and the error estimates are obtained in some norms. In , the authors considered the convergence analysis of discontinuous finite volume methods applied to distributed optimal control problems governed by a class of second-order elliptic equations.
In this paper, we consider the finite volume element method for solving the elliptic Neumann boundary control problems. Firstly, we introduce the Neumann boundary optimal control problems and their optimal conditions. To solve the optimal control problems, the associated FVE schemes are constructed. The variational discretization approach is used to deal with the control, which can avoid explicit discretization of the control and improve the approximation of the control.
The rest of the paper is organized as follows. In Section 2, we introduce some notations and the associated optimality conditions for the Neumann boundary optimal control problems. Section 3 presents the finite volume element schemes for the Neumann boundary optimal control problems. In Section 4, numerical results are presented to illustrate the effectiveness. Brief conclusions are given in Section 5.
2. Some Notations and Optimality Conditions
In the following, we use the following notations for the inner products and norms on and :
The corresponding weak formulation for (2) is: Find such that
Now, we consider the following optimal control problem for the state and the control :
over all subject to elliptic state problem (3) and the control constraints
where is a given desired state, is a regularization parameter and . We define the set of admissible controls by
where is a nonempty, closed and convex subset of .
From standard arguments for elliptic equations, we can obtain the following propositions.
Proposition 1 For fixed control , and , the state Equation (3) admits a unique solution and the following a priori estimate holds:
Proposition 2 Let be a nonempty, closed, bounded and convex set, in and , then the optimal control problem (6) admits a unique solution .
This proof follows standard techniques .
The adjoint state equation for is given by
Proposition 3 The necessary and sufficient optimality conditions for (6) and (7) can be expressed as the variational inequality
Further, the variational inequality is equivalent to
where denotes the orthogonal projection in onto the admissible set of the control, and is the adjoint operator of B.
3. FVE Method
Now we describe the finite volume element discretization of the optimal control problem (6). We consider a quasi-uniform triangulation . Divide into a sum of finite number of small triangles K such that they have no overlapping internal region and a vertex of any triangle does not belong to a side of any other triangle. At last, we can obtain a triangulation such that . The set of all nodes are denoted by
In addition, we denote by the index set of those vertices that, along with , are in some element of .
The dual partition corresponding to the primal partition is constructed as follows: In each element consisting of vertices , and , take the barycenter q in K and select the midpoint on each of the three edges . Then connect q to the points by straight lines . In this way, we obtain the dual element whose edges are . Figure 1 gives a sketch of a dual element centered at a vertex .
Figure 1. Dual element/control volume with barycenter as internal point.
We now formulate the discrete form of the optimal control problem. Let be the trial function space defined on the triangulation ,
and let be the test function space defined on the dual mesh ,
In this way, we have
where are the standard node basis functions with the nodes , and are the characteristic functions of the control volume . Let and be the usual interpolation operators, i.e.,
By the interpolation theorem of Sobolev spaces, we have for and
Then the finite volume element schemes for (3), (10) and (11) are defined as follows:
We use the following fixed-point iteration algorithm to compute the above discrete optimality problem.
4. Numerical Examples
In this section, we report some numerical results of finite volume element schemes for the Neumann boundary optimal control problems. To illustrate the theoretical analysis, the following rate of convergence r is defined,
where is the numerical solution with space step size h and u is the analytical solution.
4.1. Example 1
To validate the finite volume element schemes for the solution of elliptic problem with the Neumann boundary condition, test example is needed for which the exact solutions are known in advance. We consider the following problem with inhomogeneous Neumann boundary condition,
where Ω denotes unit square , is the outer unit normal vector. We consider the following exact solution to the above boundary value problem
In this example, we choose . The source term and the boundary condition can be derived from the exact solution.
We compute the , and errors for . They are displayed in Table 1. Examination of the table shows that the error measures of the FVE schemes diminish approximately for the error in the norm and h in the norm, which is optimal from the the viewpoint of polynomial degrees. And the error in the norm is quasi-optimal.
4.2 Example 2
This example is taken from . The optimization problem reads
Table 1. Errors results for Example 1.
where Ω denotes unit square , , , , , ,
where to are the four sides of the square, starting at the bottom side and going counterclockwise,
This problem has the following exact solution: , and
The convergence results of the control in Table 2 demonstrate second-order accuracy in the and norm. The , and errors for state and adjoint state have been computed by the finite volume element schemes. They are displayed in Table 3 and Table 4, respectively. Examination of the tables shows that the error measures of the schemes diminish approximately quadratically for the error in and norm, and linearly for the error in norm except the error of adjoint state . This may be due to its better smoothness. Figure 2 shows the computed optimal control with .
Table 2. Errors of the control for Example 2.
Table 3. Errors of the state for Example 2.
Table 4. Errors of the adjoint state for Example 2.
Figure 2. The computed optimal control .
In this article, we have considered the finite volume element approximation of the Neumann boundary optimal control problem which is governed by the elliptic partial differential equations. To obtain optimal convergence order, the variational discretization approach is used to deal with the control. Numerical results demonstrate that the proposed FVE method can effectively solve the the Neumann boundary optimal control problem. Second-order accuracy is obtained for the control in the and norm. In the future, we will consider the error estimates for the present method.
This project is partially supported by the Fundamental Research Funds for the Central Universities No. KJQN201839.
 Zhang, Z. (2009) Error Estimates of Finite Volume Element Method for the Pollution in Groundwater Flow. Numerical Methods for Partial Differential Equations, 25, 259-274.
 Arada, N., Casas, E. and Tröltzsch, F. (2002) Error Estimates for the Numerical Approximation of a Semilinear Elliptic Control Problem. Computational Optimization and Applications, 23, 201-229.
 Casas, E. and Mateos, M. (2008) Error Estimates for the Numerical Approximation of Neumann Control Problems. Computational Optimization and Applications, 39, 265-295.
 Apel, T., Pfefferer, J. and Rösch, A. (2015) Finite Element Error Estimates on the Boundary with Application to Optimal Control. Mathematics of Computation, 84, 33-70.
 Casas, E., Mateos, M. and Tröltzsch, F. (2005) Error Estimates for the Numerical Approximation of Boundary Semilinear Elliptic Control Problems. Computational Optimization and Applications, 31, 193-219.
 Hou, T. and Li, L. (2016) Error Estimates of Mixed Methods for Optimal Control Problems Governed by General Elliptic Equations. Advances in Applied Mathematics and Mechanics, 8, 1050-10711.
 Becker, R., Kapp, H. and Rannacher, R. (2000) Adaptive Finite Element Methods for Optimal Control of Partial Differential Equations: Basic Concept. SIAM Journal on Control and Optimization, 39, 113-132.
 Wang, Q., Zhao, T. and Zhang, Z. (2018) Finite Volume Element Approximation for the Elliptic Equation with Distributed Control. International Journal of Differential Equations, 2018, Article ID: 4753792.
 Sandilya, R. and Kumar, S. (2016) Convergence Analysis of Discontinuous Finite Volume Methods for Elliptic Optimal Control Problems. International Journal of Computational Methods, 13, Article ID: 1640012.
 Wang, Q., Zhang, Z. and Li, Z. (2013) A Fourier Finite Volume Element Method for Solving Two-Dimensional Quasi-Geostrophic Equations on a Sphere. Applied Numerical Mathematics, 71, 1-13.
 Apel, T., Winkler, M. and Pfefferer, J. (2018) Error Estimates for the Postprocessing Approach Applied to Neumann Boundary Control Problems in Polyhedral Domains. IMA Journal of Numerical Analysis, 38, 1984-2025.
 Luo, X., Chen, Y. and Huang, Y. (2013) Some Error Estimates of Finite Volume Element Approximation for Elliptic Optimal Control Problems. International Journal of Numerical Analysis and Modeling, 10, 697-711.
 Wang, Q., Zhang, Z., Zhang, X. and Zhu, Q. (2014) Energy-Preserving Finite Volume Element Method for the Improved Boussinesq Equation. Journal of Computational Physics, 270, 58-69.
 Ruiz-Baier, R. and Torres, H. (2014) Numerical Solution of a Multidimensional Sedimentation Problem Using Finite Volume-Element Methods. Applied Numerical Mathematics, 95, 280-291.