SIAM Journal on Control and Optimization, Vol.49, No.5, 1998-2007, 2011
A PARADOX IN THE APPROXIMATION OF DIRICHLET CONTROL PROBLEMS IN CURVED DOMAINS
In this paper, we study the approximation of a Dirichlet control problem governed by an elliptic equation defined on a curved domain Omega. To solve this problem numerically, it is usually necessary to approximate Omega by a (typically polygonal) new domain Omega(h). The difference between the solutions of both infinite-dimensional control problems, one formulated in Omega and the second in Omega(h), was studied in [E. Casas and J. Sokolowski, SIAM J. Control Optim., 48 (2010), pp. 3746-3780], where an error of order O(h) was proved. In [K. Deckelnick, A. Gunther, and M. Hinze, SIAM J. Control Optim., 48 (2009), pp. 2798-2819], the numerical approximation of the problem defined in Omega was considered. The authors used a finite element method such that Omega(h) was the polygon formed by the union of all triangles of the mesh of parameter h. They proved an error of order O(h(3/2)) for the difference between continuous and discrete optimal controls. Here we show that the estimate obtained in [E. Casas and J. Sokolowski, SIAM J. Control Optim., 48 (2010), pp. 3746-3780] cannot be improved, which leads to the paradox that the numerical solution is a better approximation of the optimal control than the exact one obtained just by changing the domain from Omega to Omega(h).