Necessary optimality conditions for discrete systems with state-dependent control region
New regularity results and improved error estimates for optimal control problems with state constraints
In this paper we are concerned with a distributed optimal control problem governed by an elliptic partial differential equation. State constraints of box type are considered. We show that the Lagrange multiplier associated with the state constraints, which is known to be a measure, is indeed more regular under quite general assumptions. We discretize the problem by continuous piecewise linear finite elements and we are able to prove that, for the case of a linear equation, the order of convergence...
New sufficient convergence conditions for the secant method
We provide new sufficient conditions for the convergence of the secant method to a locally unique solution of a nonlinear equation in a Banach space. Our new idea uses “Lipschitz-type” and center-“Lipschitz-type” instead of just “Lipschitz-type” conditions on the divided difference of the operator involved. It turns out that this way our error bounds are more precise than the earlier ones and under our convergence hypotheses we can cover cases where the earlier conditions are violated.
New system of general nonconvex variational inequalities.
Newton and conjugate gradient for harmonic maps from the disc into the sphere
We compute numerically the minimizers of the Dirichlet energyamong maps from the unit disc into the unit sphere that satisfy a boundary condition and a degree condition. We use a Sobolev gradient algorithm for the minimization and we prove that its continuous version preserves the degree. For the discretization of the problem we use continuous finite elements. We propose an original mesh-refining strategy needed to preserve the degree with the discrete version of the algorithm (which is a preconditioned...
Newton's method in the context of gradients.
Newton's methods from a geometric point of view.
Non-local approximation of free-discontinuity problems with linear growth
We approximate, in the sense of Γ-convergence, free-discontinuity functionals with linear growth in the gradient by a sequence of non-local integral functionals depending on the average of the gradients on small balls. The result extends to higher dimension what we already proved in the one-dimensional case.
Non-parametric approximation of non-anticipativity constraints in scenario-based multistage stochastic programming
We propose two methods to solve multistage stochastic programs when only a (large) finite set of scenarios is available. The usual scenario tree construction to represent non-anticipativity constraints is replaced by alternative discretization schemes coming from non-parametric estimation ideas. In the first method, a penalty term is added to the objective so as to enforce the closeness between decision variables and the Nadaraya–Watson estimation of their conditional expectation. A numerical application...
Numerical analysis for optimal shape design in elliptic boundary value problems
Shape optimization problems are optimal design problems in which the shape of the boundary plays the role of a design, i.e. the unknown part of the problem. Such problems arise in structural mechanics, acoustics, electrostatics, fluid flow and other areas of engineering and applied science. The mathematical theory of such kind of problems has been developed during the last twelve years. Recently the theory has been extended to cover also situations in which the behaviour of the system is governed...
Numerical analysis of oscillations in nonconvex problems.
Numerical analysis of some optimal control problems governed by a class of quasilinear elliptic equations
In this paper, we carry out the numerical analysis of a distributed optimal control problem governed by a quasilinear elliptic equation of non-monotone type. The goal is to prove the strong convergence of the discretization of the problem by finite elements. The main issue is to get error estimates for the discretization of the state equation. One of the difficulties in this analysis is that, in spite of the partial differential equation has a unique solution for any given control, the uniqueness...
Numerical analysis of some optimal control problems governed by a class of quasilinear elliptic equations*
In this paper, we carry out the numerical analysis of a distributed optimal control problem governed by a quasilinear elliptic equation of non-monotone type. The goal is to prove the strong convergence of the discretization of the problem by finite elements. The main issue is to get error estimates for the discretization of the state equation. One of the difficulties in this analysis is that, in spite of the partial differential equation has a unique solution for any given control, the uniqueness...
Numerical approximation of relaxed variational problems.
Numerical procedure to approximate a singular optimal control problem
In this work we deal with the numerical solution of a Hamilton-Jacobi-Bellman (HJB) equation with infinitely many solutions. To compute the maximal solution – the optimal cost of the original optimal control problem – we present a complete discrete method based on the use of some finite elements and penalization techniques.
Numerical realization of a fictitious domain approach used in shape optimization. Part I: Distributed controls
We deal with practical aspects of an approach to the numerical realization of optimal shape design problems, which is based on a combination of the fictitious domain method with the optimal control approach. Introducing a new control variable in the right-hand side of the state problem, the original problem is transformed into a new one, where all the calculations are performed on a fixed domain. Some model examples are presented.
Numerical resolution of an “unbalanced” mass transport problem
We introduce a modification of the Monge–Kantorovitch problem of exponent 2 which accommodates non balanced initial and final densities. The augmented lagrangian numerical method introduced in [6] is adapted to this “unbalanced” problem. We illustrate the usability of this method on an idealized error estimation problem in meteorology.
Numerical resolution of an “unbalanced” mass transport problem
We introduce a modification of the Monge–Kantorovitch problem of exponent 2 which accommodates non balanced initial and final densities. The augmented Lagrangian numerical method introduced in [6] is adapted to this “unbalanced” problem. We illustrate the usability of this method on an idealized error estimation problem in meteorology.
Numerical study of discretizations of multistage stochastic programs
This paper presents a numerical study of a deterministic discretization procedure for multistage stochastic programs where the underlying stochastic process has a continuous probability distribution. The discretization procedure is based on quasi-Monte Carlo techniques originally developed for numerical multivariate integration. The solutions of the discretized problems are evaluated by statistical bounds obtained from random sample average approximations and out-of-sample simulations. In the numerical...