We derive local error estimates for the discretization of optimal control problems governed by linear advection-diffusion partial differential equations (PDEs) using the streamline upwind/Petrov Galerkin (SUPG) stabilized finite element method. We show that if the SUPG method is used to solve optimization problems governed by an advection-dominated PDE the convergence properties of the SUPG method is substantially different from the convergence properties of the SUPG method applied for the solution of an advection-dominated PDE. The reason is that the solution of the optimal control problem involves another advection dominated PDE, the so-called adjoint equation, whose advection field is just the negative of the advection of the governing PDEs. For the solution of the optimal control problem, a coupled system involving both the original governing PDE as well as the adjoint PDE must be solved. We show that in the presence of a boundary layer, the local error between the solution of the SUPG discretized optimal control problem and the solution of the infinite dimensional problem is only of first order even if the error is computed locally in a region away from the boundary layer. In the presence of interior layers, we prove optimal convergence rates for the local error in a region away from the layer between the solution of the SUPG discretized optimal control problems and the solution of the infinite dimensional problem. Numerical examples are presented to illustrate some of the theoretical results.

In this paper we address several issues arising from a singularly perturbed fourth order problem with small parameter ε. First, we introduce a new family of nonconforming elements. We then prove that the corresponding finite element method is robust with respect to the parameter ε and uniformly convergent to order h 1/2 . In addition, we analyze the effect of treating the Neumann boundary condition weakly by Nitsche's method. We show that such treatment is superior when the parameter ε is smaller than the mesh size h and obtain sharper error estimates. Such error analysis is not restricted to the proposed elements and can easily be carried out to other elements as long as the Neumann boundary condition is imposed weakly. Finally, we discuss the local error estimates and the pollution effect of the boundary layers in the interior of the domain.

Abstract. A model second-order elliptic equation on a general convex polyhedral domain in three dimensions is considered. The aim of this paper is twofold: First sharp Hölder estimates for the corresponding Green's function are obtained. As an applications of these estimates to finite element methods, we show the best approximation property of the error in W 1 ∞ . In contrast to previously known results, W 2 p regularity for p > 3, which does not hold for general convex polyhedral domains, is not required. Furthermore, the new Green's function estimates allow us to obtain localized error estimates at a point.

The aim of the paper is to show the stability of the finite element solution for the Stokes system in W 1 ∞ norm on general convex polyhedral domain. In contrast to previously known results, W 2 r regularity for r > 3, which does not hold for general convex polyhedral domains, is not required. The argument uses recently available sharp Hölder pointwise estimates of the corresponding Green's matrix together with novel local energy error estimates, which do not involve an error of the pressure in a weaker norm.

In this paper we consider a parabolic optimal control problem with a pointwise (Dirac type) control in space, but variable in time, in two space dimensions. To approximate the problem we use the standard continuous piecewise linear approximation in space and the piecewise constant discontinuous Galerkin method in time. Despite low regularity of the state equation, we show almost optimal h 2 + k convergence rate for the control in L 2 norm. This result improves almost twice the previously known estimate in [W.

The main goal of the paper is to establish time semidiscrete and space-time fully discrete maximal parabolic regularity for the time discontinuous Galerkin solution of linear parabolic equations. Such estimates have many applications. They are essential, for example, in establishing optimal a priori error estimates in nonHilbertian norms without unnatural coupling of spatial mesh sizes with time steps.Keywords maximal parabolic regularity · finite elements · maximum norm · fully discrete · resolvent estimates · resolvent estimates · optimal error estimates · parabolic smoothing

Abstract. We consider finite element methods for a model second-order elliptic equation on a general bounded convex polygonal or polyhedral domain. Our first main goal is to extend the best approximation property of the error in the W 1 ∞ norm, which is known to hold on quasi-uniform meshes, to more general graded meshes. We accomplish it by a novel proof technique. This result holds under a condition on the grid which is mildly more restrictive than the shape regularity condition typically enforced in adaptive codes. The second main contribution of this work is a discussion of the properties of and relationships between similar mesh restrictions that have appeared in the literature.

In this paper we consider a problem of initial data identification from the final time observation for homogeneous parabolic problems. It is well-known that such problems are exponentially ill-posed due to the strong smoothing property of parabolic equations. We are interested in a situation when the initial data we intend to recover is known to be sparse, i.e. its support has Lebesgue measure zero. We formulate the problem as an optimal control problem and incorporate the information on the sparsity of the unknown initial data into the structure of the objective functional. In particular, we are looking for the control variable in the space of regular Borel measures and use the corresponding norm as a regularization term in the objective functional. This leads to a convex but non-smooth optimization problem. For the discretization we use continuous piecewise linear finite elements in space and discontinuous Galerkin finite elements of arbitrary degree in time. For the general case we establish error estimates for the state variable. Under a certain structural assumption, we show that the control variable consists of a finite linear combination of Dirac measures. For this case we obtain error estimates for the locations of Dirac measures as well as for the corresponding coefficients. The key to the numerical analysis are the sharp smoothing type pointwise finite element error estimates for homogeneous parabolic problems, which are of independent interest. Moreover, we discuss an efficient algorithmic approach to the problem and show several numerical experiments illustrating our theoretical results.