Refine
Language
- English (23)
Keywords
- optimal control (6)
- state constraints (4)
- Dynamical contact problems (2)
- Newmark method (2)
- Signorini condition (2)
- adaptivity (2)
- coupling conditions (2)
- energy minimization (2)
- path-following (2)
- (visco-)elasticity (1)
The paper presents a new affine invariant theory on asymptotic mesh
independence of Newton’s method for discretized nonlinear operator equations.
Compared to earlier attempts, the new approach is both much simpler
and more intuitive from the algorithmic point of view. The theory
is exemplified at collocation methods for ODE boundary value problems
and at finite element methods for elliptic PDE problems.
We consider first order optimality conditions for state constrained optimal control problems. In particular we study the case where the state equation has not enough regularity to admit existence of a Slater point in function space. We overcome this difficulty by a special transformation. Under a density condition we show existence of Lagrange multipliers, which have a representation via measures and additional regularity properties.
In this paper we are concerned with the application of interior point methods in function space to gradient constrained optimal control problems, governed by partial differential equations. We will derive existence of solutions together with first order optimality conditions. Afterwards we show continuity of the central path, together with convergence rates depending on the interior point parameter.
Interior Point Methods in Function Space for State Constraints - Inexact Newton and Adaptivity
(2009)
We consider an interior point method in function space for PDE constrained optimal control problems with state constraints. Our emphasis is on the construction and analysis of an algorithm that integrates a Newton path-following method with adaptive grid refinement. This is done in the framework of inexact Newton methods in function space, where the discretization error of each Newton step is controlled by adaptive grid refinement in the innermost loop. This allows to perform most of the required Newton steps on coarse grids, such that the overall computational time is dominated by the last few steps. For this purpose we propose an a-posteriori error estimator for a problem suited norm.
A continuity result for Nemyckii Operators and some applications in PDE constrained optimal control
(2008)
This work explores two applications of a classical result on the continuity of Nemyckii operators to optimal control with PDEs. First, we present an alternative approach to the analysis of Newton's method for function space problems involving semi-smooth Nemyckii operators. A concise proof for superlinear convergence is presented, and sharpened bounds on the rate of convergence are derived. Second, we derive second order sufficient conditions for problems, where the underlying PDE has poor regularity properties. We point out that the analytical structure in both topics is essentially the same.
We propose and analyse an interior point path-following method in function space
for state constrained optimal control. Our emphasis is on proving convergence in
function space and on constructing a practical path-following algorithm. In particular, the introduction of a pointwise damping step leads to a very efficient method, as verified by numerical experiments.
This paper is intended to be a first step towards the continuous dependence of dynamical contact problems on the initial data as well as the uniqueness of a solution. Moreover, it provides the basis for a proof of the convergence of popular time integration schemes as the Newmark method.
We study a frictionless dynamical contact problem between both linearly elastic and viscoelastic bodies which is formulated via the Signorini contact conditions. For viscoelastic materials fulfilling the Kelvin-Voigt constitutive law, we find a characterization of the class of problems which satisfy a perturbation result in a non-trivial mix of norms in function space. This characterization is given in the form of a stability condition on the contact stresses at the contact boundaries.
Furthermore, we present perturbation results for two well-established approximations of the classical Signorini condition: The Signorini condition formulated in velocities and the model of normal compliance, both satisfying even a sharper version of our stability condition.
The minimization of an L^{\infty}-functional subject to an elliptic PDE and state constraints
(2008)
We study the optimal control of a maximum-norm objective functional subjected to an elliptic-type PDE and pointwise state constraints. The problem is transformed into a problem where the non-differentiable L^{\infty}-norm in the functional will be replaced by a scalar variable and additional state constraints. This problem is solved by barrier methods. We will show the existence and convergence of the central path for a class of barrier functions. Numerical experiments complete the presentation.
We study barrier methods for state constrained optimal control problems
with PDEs. In the focus of our analysis is the path of minimizers of the
barrier subproblems with the aim to provide a solid theoretical basis for function
space oriented path-following algorithms. We establish results on existence,
continuity, and convergence of this path. Moreover, we consider the structure of
barrier subdifferentials, which play the role of dual variables.
An extended mathematical framework for barrier methods for state constrained optimal control compared to [Schiela, ZIB-Report 07-07] is considered. This allows to apply the results derived there to more general classes of optimal control problems, in particular to boundary control and finite dimensional control.
We discuss first order optimality conditions for state constrained
optimal control problems. Our concern is the treatment of problems,
where the solution of the state equation is not known to be continuous,
as in the case of boundary control in three space dimensions or optimal
control with parabolic partial differential equations. We show existence of
measure valued Lagrangian multipliers, which have just enough additional
regularity to be applicable to all possibly discontinuous solutions of the
state equation.
The aim of this paper is to devise an adaptive timestep control in the contact--stabilized Newmark method (CONTACX) for dynamical contact problems between two viscoelastic bodies in the framework of Signorini's condition. In order to construct a comparative scheme of higher order accuracy, we extend extrapolation techniques. This approach demands a subtle theoretical investigation of an asymptotic error expansion of the contact--stabilized Newmark scheme. On the basis of theoretical insight and numerical observations, we suggest an error estimator and a timestep selection which also cover the presence of contact. Finally, we give a numerical example.
This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study.
We formulate the static mechanical coupling of a geometrically exact Cosserat rod
to a nonlinearly elastic continuum. In this setting, appropriate coupling conditions have
to connect a one-dimensional model with director variables to a three-dimensional
model without directors.
Two alternative coupling conditions are proposed,
which correspond to two different configuration trace spaces.
For both we show existence of solutions of the coupled problems, using the direct
method of the calculus of variations. From the first-order optimality conditions
we also derive the corresponding conditions for the dual variables. These are
then interpreted in mechanical terms.
We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented.
A primal interior point method for control constrained optimal control problems
with PDE constraints is considered. Pointwise elimination of the control
leads to a homotopy in the remaining state and dual variables, which is addressed
by a short step pathfollowing method. The algorithm is applied to the
continuous, infinite dimensional problem, where discretization is performed
only in the innermost loop when solving linear equations. The a priori elimination
of the least regular control permits to obtain the required accuracy with
comparatively coarse meshes. Convergence of the method and discretization
errors are studied, and the method is illustrated at two numerical examples.
A thorough convergence analysis of the Control Reduced Interior Point
Method in function space is performed. This recently proposed method is a
primal interior point pathfollowing scheme with the special feature, that the
control variable is eliminated from the optimality system. Apart from global
linear convergence we show, that this method converges locally almost quadratically,
if the optimal solution satisfies a certain non-degeneracy condition. In
numerical experiments we observe, that a prototype implementation of our
method behaves as predicted by our theoretical results.
The paper considers the time integration of frictionless dynamical contact problems between viscoelastic bodies in the frame of the Signorini condition. Among the numerical integrators, interest focuses on the contact-stabilized Newmark method recently suggested by Deuflhard et al., which is compared to the classical Newmark method and an improved energy dissipative version due to Kane et al. In the absence of contact, any such variant is equivalent to the Störmer-Verlet scheme, which is well-known to have consistency order 2. In the presence of contact, however, the classical approach to discretization errors would not show consistency at all because of the discontinuity at the contact. Surprisingly, the question of consistency in the constrained situation has not been solved yet. The present paper fills this gap by means of a novel proof technique using specific norms based on earlier perturbation results due to the authors. The corresponding estimation of the local discretization error requires the bounded total variation of the solution. The results have consequences for the construction of an adaptive timestep control, which will be worked out subsequently in a forthcoming paper.
The well known De Giorgi result on Hölder continuity for solutions of the Dirichlet problem is re-established for mixed boundary value problems, provided that the underlying domain is a Lipschitz domain and the border between the Dirichlet and the Neumann boundary part satisfies a very general geometric condition. Implications of this result for optimal control theory are presented.
Convergence Analysis of Smoothing Methods for Optimal Control of Stationary Variational Inequalities
(2012)
In the article an optimal control problem subject to a stationary variational inequality is investigated. The optimal control problem is complemented with pointwise
control constraints. The convergence of a smoothing scheme is analyzed. There, the variational inequality is replaced by a semilinear elliptic equation. It is shown that solutions of the regularized optimal control problem converge to solutions of the original one. Passing to
the limit in the optimality system of the regularized problem allows to prove C-stationarity of local solutions of the original problem. Moreover, convergence rates with respect to the regularization parameter for the error in the control are obtained. These rates coincide with
rates obtained by numerical experiments, which are included in the paper.