On saturation effects in the Neumann boundary control of elliptic optimal control problems

A Neumann boundary control problem for a linear-quadratic elliptic optimal control problem in a polygonal domain is investigated. The main goal is to show an optimal approximation order for discretized problems after a postprocessing process. It turns out that two saturation processes occur: The regularity of the boundary data of the adjoint is limited if the largest angle of the polygon is at least 2π/3. Moreover, piecewise linear finite elements cannot guarantee the optimal order, if the largest angle of the polygon is greater than π/2. We will derive error estimates of order h ^α with α∈[1,2] depending on the largest angle and properties of the finite elements. Finally, numerical test illustrates the theoretical results.     

A Priori Error Estimates for Finite Element Methods for H (2,1)-Elliptic Equations

The convergence of finite element methods for linear elliptic boundary value problems of second and forth order is well understood. In this article, we introduce finite element approximations of some linear semi-elliptic boundary value problem of mixed order on a two-dimensional rectangular domain Q. The equation is of second order in one direction and forth order in the other and appears in the optimal control of parabolic partial differential equations if one eliminates the control and the state (or the adjoint state) in the first order optimality conditions. We establish a regularity result and estimate for the finite element error of conforming approximations of this equation. The finite elements in use have a tensor product structure, in one dimension we use linear, quadratic or cubic Lagrange elements in the other dimension cubic Hermite elements. For these elements, we prove the error bound O(h ² + τ ^k ) in the energy norm and O((h ² + τ ^k )(h ² + τ)) in the L ²(Q)-norm.     

The Mortar Element Method with Locally Nonconforming Elements

We consider a discretization of linear elliptic boundary value problems in 2-D by the new version of the mortar finite element method which uses locally nonconforming Crouzeix-Raviart elements. We show that if a solution of the original differential problem belongs to the space H ²(), then an error is of the same order as in the standard nonconforming finite element method. We also propose an additive Schwarz method of solving the discrete problem and show that its rate of convergence is almost optimal.     

Local mesh refinement for the discretization of Neumann boundary control problems on polyhedra

This paper deals with optimal control problems constrained by linear elliptic partial differential equations. The case where the right‐hand side of the Neumann boundary is controlled, is studied. The variational discretization concept for these problems is applied, and discretization error estimates are derived. On polyhedral domains, one has to deal with edge and corner singularities, which reduce the convergence rate of the discrete solutions, that is, one cannot expect convergence order two for linear finite elements on quasi‐uniform meshes in general. As a remedy, a local mesh refinement strategy is presented, and a priori bounds for the refinement parameters are derived such that convergence with optimal rate is guaranteed. As a by‐product, finite element error estimates in the H¹(Ω)‐norm, L²(Ω)‐norm and L²(Γ)‐norm for the boundary value problem are obtained, where the latter one turned out to be the main challenge. Copyright © 2015 John Wiley & Sons, Ltd.     

A Priori Error Estimates of Crank–Nicolson Finite Volume Element Method for a Hyperbolic Optimal Control Problem

In this article, a Crank–Nicolson linear finite volume element scheme is developed to solve a hyperbolic optimal control problem. We use the variational discretization technique for the approximation of the control variable. The optimal convergent order O(h² + k²) is proved for the numerical solution of the control, state and adjoint‐state in a discrete L²‐norm. To derive this result, we also get the error estimate (convergent order O(h² + k²)) of Crank–Nicolson finite volume element approximation for the second‐order hyperbolic initial boundary value problem. Numerical experiments are presented to verify the theoretical results.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1331–1356, 2016     

Two-grid methods for –P1 mixed finite element approximation of general elliptic optimal control problems with low regularity

In this paper, we present a two-grid mixed finite element scheme for distributed optimal control governed by general elliptic equations. –P₁ mixed finite elements are used for the discretization of the state and co-state variables, whereas piecewise constant function is used to approximate the control variable. We first use a new approach to obtain the superclose property between the centroid interpolation and the numerical solution of the optimal control u with order h² under the low regularity. Based on the superclose property, we derive the optimal a priori error estimates. Then, using a postprocessing projection operator, we get a second-order superconvergent result for the control u. Next, we construct a two-grid mixed finite element scheme and analyze a priori error estimates. In the two-grid scheme, the solution of the elliptic optimal control problem on a fine grid is reduced to the solution of the elliptic optimal control problem on a much coarser grid and the solution of a linear algebraic system on the fine grid and the resulting solution still maintains an asymptotically optimal accuracy. Finally, a numerical example is presented to verify the theoretical results.     

Error estimates for the discretization of elliptic control problems with pointwise control and state constraints

A family of elliptic optimal control problems with pointwise constraints on control and state is considered. We are interested in approximation of the optimal solution by a finite element discretization of the involved partial differential equations. The discretization error for a problem with mixed state constraints is estimated in the semidiscrete case and in the fully discrete scheme with the convergence of order h|ln h| and h ^1/2, respectively. However, considering the unregularized continuous problem and the discrete regularized version, and choosing suitable relation between the regularization parameter and the mesh size, i.e., ε∼h ², a convergence order arbitrary close to 1, i.e., h ^1−β is obtained. Therefore, we benefit from tuning the involved parameters.     

Boundary concentrated finite elements for optimal boundary control problems of elliptic PDEs

We investigate the discretization of optimal boundary control problems for elliptic equations on two-dimensional polygonal domains by the boundary concentrated finite element method. We prove that the discretization error ||u^*-u_h^*||_L²(G)\|u^{*}-u_{h}^{*}\|_{L^{2}(\Gamma)} decreases like N ⁻¹, where N is the total number of unknowns. This makes the proposed method favorable in comparison to the h-version of the finite element method, where the discretization error behaves like N ^−3/4 for uniform meshes. Moreover, we present an algorithm that solves the discretized problem in almost optimal complexity. The paper is complemented with numerical results.     

Finite groups in which any two elements of the same order are either fused or inverse-fused

Elements a, b of a group G are said to be fused or inverse-fused if there exists σεAut(G) such that a = b^σ or a = (b^-1)^σ respectively. This paper gives a classification of all finite groups in which any two elements of the same order are fused orinverse-fused.     

A finite element method for the solution of a potential theory integral equation

This paper discusses a finite element approximation for an integral equation of the second kind deduced from a potential theory boundary value problem in two variables. The equation is shown to admit a unique solution, to be variational and coercive in the Hilbert space of functions σ ε H^1/2(Γ), f_rd γ = 0. The Galerkin method with finite elements as trial functions is shown to lead to an optimal rate of convergence.     

Composite finite element approximation for nonlinear parabolic problems in nonconvex polygonal domains

We study a new class of finite elements so‐called composite finite elements (CFEs), introduced earlier by Hackbusch and Sauter, Numer. Math., 1997; 75:447‐472, for the approximation of nonlinear parabolic equation in a nonconvex polygonal domain. A two‐scale CFE discretization is used for the space discretizations, where the coarse‐scale grid discretized the domain at an appropriate distance from the boundary and the fine‐scale grid is used to resolve the boundary. A continuous, piecewise linear CFE space is employed for the spatially semidiscrete finite element approximation and the temporal discretizations is based on modified linearized backward Euler scheme. We derive almost optimal‐order convergence in space and optimal order in time for the CFE method in the L^∞(L²) norm. Numerical experiment is carried out for an L‐shaped domain to illustrate our theoretical findings.     

An adaptive algorithm for the Crank–Nicolson scheme applied to a time-dependent convection–diffusion problem

An a posteriori upper bound is derived for the nonstationary convection–diffusion problem using the Crank–Nicolson scheme and continuous, piecewise linear stabilized finite elements with large aspect ratio. Following Lozinski et al. (2009) [13], a quadratic time reconstruction is used.A space and time adaptive algorithm is developed to ensure the control of the relative error in the L²(H¹) norm. Numerical experiments illustrating the efficiency of this approach are reported; it is shown that the error indicator is of optimal order with respect to both the mesh size and the time step, even in the convection dominated regime and in the presence of boundary layers.     

The finite simple groups with at most two fusion classes of every order

Elements a,b of a group G are said to be fused if a = b^σ and to be inverse-fused if a =(b^-1)^σ for some σ ? Aut(G). The fusion class of a ? G is the set {a^σ | σ ? Aut(G)}, and it is called a fusion class of order i if a has order iThis paper gives a complete classification of the finite nonabelian simple groups G for which either (i) or (ii) holds, where:

(i) G has at most two fusion classes of order i for every i (23 examples); and

(ii) any two elements of G of the same order are fused or inversenfused.

The examples in case (ii) are: A₅, A₆,L₂(7),L₂(8), L₃(4), Sz(8), M₁₁ and M₂₃An application is given concerning isomorphisms of Cay ley graphs.     

Optimization of a time-variable,constrained control matrix

The problem posed is to choose, in a optimal manner, a time-variable, bounded, linear transformation defining the velocity of a state point inn-dimensional space in terms of the state. The two-point boundary-value problem which arises from an application of the Pontryagin maximum principle is explicitly solvable; hence, a formula is derived showing that the optimal trajectories in state space are equiangular spirals in two-dimensional subspaces ofR ⁿ and also describing the boundary of the set of attainability. This formula is used to solve the problem of minimal-time transfer between any two given points, and the optimal control is specified both as an open-loop and a closed-loop controller. The solutions to the problem of maximizing a linear payoff function of the final state and of maximizing the angle of rotation of the state vector about the origin are also given.     

An implicit–explicit residual error estimator for the coupling of dual‐mixed finite elements and boundary elements in elastostatics

We consider the coupling of dual‐mixed finite elements and boundary elements to solve a mixed Dirichlet–Neumann problem of plane elasticity. We derive an a‐posteriori error estimate that is based on the solution of local Dirichlet problems and on a residual term defined on the coupling interface. The general error estimate does not make use of any special finite element or boundary element spaces. Here the residual term is given in a negative order Sobolev norm. In practical applications, where a certain boundary element subspace is used, this norm can be estimated by weighted local L²‐norms. Copyright © 2001 John Wiley & Sons, Ltd.     

A priori error estimates for optimal control problems with pointwise constraints on the gradient of the state

We analyze a finite element approximation of an elliptic optimal control problem with pointwise bounds on the gradient of the state variable. We derive convergence rates if the control space is discretized implicitly by the state equation. In contrast to prior work we obtain these results directly from classical results for the W ^1,∞-error of the finite element projection, without using adjoint information. If the control space is discretized directly, we first prove a regularity result for the optimal control to control the approximation error, based on which we then obtain analogous convergence rates.     

A family of higher order mixed finite element methods for plane elasticity

Summary The Dirichler problem for the equations of plane elasticity is approximated by a mixed finite element method using a new family of composite finite elements having properties analogous to those possessed by the Raviart-Thomas mixed finite elements for a scalar, second-order elliptic equation. Estimates of optimal order and minimal regularity are derived for the errors in the displacement vector and the stress tensor inL ²(), and optimal order negative norm estimates are obtained inH ^s () for a range ofs depending on the index of the finite element space. An optimal order estimate inL () for the displacement error is given. Also, a quasioptimal estimate is derived in an appropriate space. All estimates are valid uniformly with respect to the compressibility and apply in the incompressible case. The formulation of the elements is presented in detail.This work was performed while Professor Arnold was a NATO Postdoctoral Fellow     

A finite difference method for the wide‐angle “parabolic” equation in a waveguide with downsloping bottom

We consider the third‐order wide‐angle “parabolic” equation of underwater acoustics in a cylindrically symmetric fluid medium over a bottom of range‐dependent bathymetry. It is known that the initial‐boundary‐value problem for this equation may not be well posed in the case of (smooth) bottom profiles of arbitrary shape, if it is just posed e.g. with a homogeneous Dirichlet bottom boundary condition. In this article, we concentrate on downsloping bottom profiles and propose an additional boundary condition that yields a well‐posed problem, in fact making it L² ‐conservative in the case of appropriate real parameters. We solve the problem numerically by a Crank–Nicolson‐type finite difference scheme, which is proved to be unconditionally stable and second‐order accurate and simulates accurately realistic underwater acoustic problems. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

The first Steklov eigenvalue, conformal geometry, and minimal surfaces

We consider the relationship of the geometry of compact Riemannian manifolds with boundary to the first nonzero eigenvalue σ₁ of the Dirichlet-to-Neumann map (Steklov eigenvalue). For surfaces Σ with genus γ and k boundary components we obtain the upper bound σ₁L(∂Σ)?2(γ+k)π. For γ=0 and k=1 this result was obtained by Weinstock in 1954, and is sharp. We attempt to find the best constant in this inequality for annular surfaces (γ=0 and k=2). For rotationally symmetric metrics we show that the best constant is achieved by the induced metric on the portion of the catenoid centered at the origin which meets a sphere orthogonally and hence is a solution of the free boundary problem for the area functional in the ball. For a general class of (not necessarily rotationally symmetric) metrics on the annulus, which we call supercritical, we prove that σ₁(Σ)L(∂Σ) is dominated by that of the critical catenoid with equality if and only if the annulus is conformally equivalent to the critical catenoid by a conformal transformation which is an isometry on the boundary. Motivated by the annulus case, we show that a proper submanifold of the ball is immersed by Steklov eigenfunctions if and only if it is a free boundary solution. We then prove general upper bounds for conformal metrics on manifolds of any dimension which can be properly conformally immersed into the unit ball in terms of certain conformal volume quantities. We show that these bounds are only achieved when the manifold is minimally immersed by first Steklov eigenfunctions. We also use these ideas to show that any free boundary solution in two dimensions has area at least π, and we observe that this implies the sharp isoperimetric inequality for free boundary solutions in the two-dimensional case.     

Linear and Discontinuous Approximations for Optimal Control Problems

Abstract

An optimal control problem for 2D and 3D elliptic equations is investigated with pointwise control constraints. This paper is concerned with the discretization of the control by piecewise linear but discontinuous functions. The state and the adjoint state are discretized by linear finite elements. The paper is focused on similarities and differences to piecewise constant and piecewise linear (continuous) approximation of the controls. Approximation of order h in the L ^∞-norm is proved in the main result.