Approximation of boundary control problems on curved domains

In this paper we consider boundary control problems associated to a semilinear elliptic equation defined in a curved domain Ω. The Dirichlet and Neumann cases are analyzed. To deal with the numerical analysis of these problems, the approximation of Ω by an appropriate domain Ωh (typically polygonal) is required. Here we do not consider the numerical approximation of the control problems. Instead, we formulate the corresponding infinite dimensional control problems in Ωh, and we study the influence of the replacement of Ω by Ωh on the solutions of the control problems. Our goal is to compare the optimal controls defined on Γ = ∂Ω with those defined on Γh = ∂Ωh and to derive some error estimates. The use of a convenient parametrization of the boundary is needed for such estimates

Error estimates for the numerical approximation of Dirichlet boundary control for semilinear elliptic equations

We study the numerical approximation of boundary optimal control problems governed by semilinear elliptic partial differential equations with pointwise constraints on the control. The control is the trace of the state on the boundary of the domain, which is assumed to be a convex, polygonal, open set in R2. Piecewise linear finite elements are used to approximate the control as well as the state. We prove that the error estimates are of order O(h1−1/p) for some p > 2, which is consistent with the W1−1/p,p(Γ)-regularity of the optimal control

Error estimates for the numerical approximation of boundary semilinear elliptic control problems

We study the numerical approximation of boundary optimal control problems governed by semilinear elliptic partial differential equations with pointwise constraints on the control. The analysis of the approximate control problems is carried out. The uniform convergence of discretized controls to optimal controls is proven under natural assumptions by taking piecewise constant controls. Finally, error estimates are established

Superconvergence for Neumann boundary control problems governed by semilinear elliptic equations

This paper is concerned with the discretization error analysis of semilinear Neumann boundary control problems in polygonal domains with pointwise inequality constraints on the control. The approximations of the control are piecewise constant functions. The state and adjoint state are discretized by piecewise linear finite elements. In a postprocessing step approximations of locally optimal controls of the continuous optimal control problem are constructed by the projection of the respective discrete adjoint state. Although the quality of the approximations is in general affected by corner singularities a convergence order of $h^2|\ln h|^{3/2}$ is proven for domains with interior angles smaller than $2\pi/3$ using quasi-uniform meshes. For larger interior angles mesh grading techniques are used to get the same order of convergence

Discontinuous Galerkin Methods for Mass Transfer through Semi-Permeable Membranes

A discontinuous Galerkin (dG) method for the numerical solution of initial/boundary value multi-compartment partial differential equation (PDE) models, interconnected with interface conditions, is presented and analysed. The study of interface problems is motivated by models of mass transfer of solutes through semi-permeable membranes. More specifically, a model problem consisting of a system of semilinear parabolic advection-diffusion-reaction partial differential equations in each compartment, equipped with respective initial and boundary conditions, is considered. Nonlinear interface conditions modelling selective permeability, congestion and partial reflection are applied to the compartment interfaces. An interior penalty dG method is presented for this problem and it is analysed in the space-discrete setting. The a priori analysis shows that the method yields optimal a priori bounds, provided the exact solution is sufficiently smooth. Numerical experiments indicate agreement with the theoretical bounds and highlight the stability of the numerical method in the advection-dominated regime