Source author record

Emmanuel Trélat

Emmanuel Trélat appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2024arXiv

Approximate control of parabolic equations with on-off shape controls by Fenchel duality

We consider the internal control of linear parabolic equations through on-off shape controls, i.e., controls of the form $M(t)χ_{ω(t)}$ with $M(t) \geq 0$ and $ω(t)$ with a prescribed maximal measure. We establish small-time approximate controllability towards all possible final states allowed by the comparison principle with nonnegative controls. We manage to build controls with constant amplitude $M(t) \equiv M$. In contrast, if the moving control set $ω(t)$ is confined to evolve in some region of the whole domain, we prove that approximate controllability fails to hold for small times. The method of proof is constructive. Using Fenchel-Rockafellar duality and the bathtub principle, the on-off shape control is obtained as the bang-bang solution of an optimal control problem, which we design by relaxing the constraints. Our optimal control approach is outlined in a rather general form for linear constrained control problems, paving the way for generalisations and applications to other PDEs and constraints.

preprint2023arXiv

Linear turnpike theorem

The turnpike phenomenon stipulates that the solution of an optimal control problem in large time, remains essentially close to a steady-state of the dynamics, itself being the optimal solution of an associated static optimal control problem. Under general assumptions, it is known that not only the optimal state and the optimal control, but also the adjoint state coming from the application of the Pontryagin maximum principle, are exponentially close to a steady-state, except at the beginning and at the end of the time frame. In such results, the turnpike set is a singleton, which is a steady-state. In this paper, we establish a turnpike result for finite-dimensional optimal control problems in which some of the coordinates evolve in a monotone way, and some others are partial steady-states of the dynamics. We prove that the discrepancy between the optimal trajectory and the turnpike set is then linear, but not exponential: we thus speak of a linear turnpike theorem.

preprint2022arXiv

Optimization of vaccination for COVID-19 in the midst of a pandemic

During the Covid-19 pandemic a key role is played by vaccination to combat the virus. There are many possible policies for prioritizing vaccines, and different criteria for optimization: minimize death, time to herd immunity, functioning of the health system. Using an age-structured population compartmental finite-dimensional optimal control model, our results suggest that the eldest to youngest vaccination policy is optimal to minimize deaths. Our model includes the possible infection of vaccinated populations. We apply our model to real-life data from the US Census for New Jersey and Florida, which have a significantly different population structure. We also provide various estimates of the number of lives saved by optimizing the vaccine schedule and compared to no vaccination.

preprint2022arXiv

Quantum Limits on product manifolds

We establish some properties of quantum limits on a product manifold, proving for instance that, under appropriate assumptions, the quantum limits on the product of manifolds are absolutely continuous if the quantum limits on each manifolds are absolutely continuous. On a product of Riemannian manifolds satisfying the minimal multiplicity property, we prove that a periodic geodesic can never be charged by a quantum limit.

preprint2021arXiv

Spiraling of sub-Riemannian geodesics around the Reeb flow in the 3D contact case

We consider a closed three-dimensional contact sub-Riemannian manifold. The objective of this note is to provide a precise description of the sub-Riemannian geodesics with large initial momenta: we prove that they "spiral around the Reeb orbits", not only in the phase space but also in the configuration space. Our analysis is based on a normal form along any Reeb orbit due to Melrose.

preprint2020arXiv

Observability for generalized Schrödinger equations and quantum limits on product manifolds

Given a closed product Riemannian manifold N = M x M equipped with the product Riemannian metric g = h + h , we explore the observability properties for the generalized Schr{ö}dinger equation i$\partial$ t u = F (g)u, where g is the Laplace-Beltrami operator on N and F : [0, +$\infty$) $\rightarrow$ [0, +$\infty$) is an increasing function. In this note, we prove observability in finite time on any open subset $ω$ satisfying the so-called Vertical Geometric Control Condition, stipulating that any vertical geodesic meets $ω$, under the additional assumption that the spectrum of F (g) satisfies a gap condition. A first consequence is that observability on $ω$ for the Schr{ö}dinger equation is a strictly weaker property than the usual Geometric Control Condition on any product of spheres. A second consequence is that the Dirac measure along any geodesic of N is never a quantum limit.

preprint2020arXiv

PI regulation control of a 1-D semilinear wave equation

This paper is concerned with the Proportional Integral (PI) regulation control of the left Neu-mann trace of a one-dimensional semilinear wave equation. The control input is selected as the right Neumann trace. The control design goes as follows. First, a preliminary (classical) velocity feedback is applied in order to shift all but a finite number of the eivenvalues of the underlying unbounded operator into the open left half-plane. We then leverage on the projection of the system trajectories into an adequate Riesz basis to obtain a truncated model of the system capturing the remaining unstable modes. Local stability of the resulting closed-loop infinite-dimensional system composed of the semilinear wave equation, the preliminary velocity feedback, and the PI controller, is obtained through the study of an adequate Lyapunov function. Finally, an estimate assessing the set point tracking performance of the left Neumann trace is derived.

preprint2020arXiv

Shape turnpike for linear parabolic PDE models

We introduce and study the turnpike property for time-varying shapes, within the viewpoint of optimal control. We focus here on second-order linear parabolic equations where the shape acts as a source term and we seek the optimal time-varying shape that minimizes a quadratic criterion. We first establish existence of optimal solutions under some appropriate sufficient conditions. We then provide necessary conditions for optimality in terms of adjoint equations and, using the concept of strict dissipativity, we prove that state and adjoint satisfy the measure-turnpike property, meaning that the extremal time-varying solution remains essentially close to the optimal solution of an associated static problem. We show that the optimal shape enjoys the exponential turnpike property in term of Hausdorff distance for a Mayer quadratic cost. We illustrate the turnpike phenomenon in optimal shape design with several numerical simulations.

preprint2020arXiv

Small-time asymptotics of hypoelliptic heat kernels near the diagonal, nilpotentization and related results

We establish small-time asymptotic expansions for heat kernels of hypoelliptic Hörmander operators in a neighborhood of the diagonal, generalizing former results obtained in particular by Métivier and by Ben Arous. The coefficients of our expansions are identified in terms of the nilpotentization of the underlying sub-Riemannian structure. Our approach is purely analytic and relies in particular on local and global subelliptic estimates as well as on the local nature of small-time asymptotics of heat kernels. The fact that our expansions are valid not only along the diagonal but in an asymptotic neighborhood of the diagonal is the main novelty, useful in view of deriving Weyl laws for subelliptic Laplacians. In turn, we establish a number of other results on hypoelliptic heat kernels that are interesting in themselves, such as Kac's principle of not feeling the boundary, asymptotic results for singular perturbations of hypoelliptic operators, global smoothing properties for selfadjoint heat semigroups.

preprint2020arXiv

Unified Riccati theory for optimal permanent and sampled-data control problems in finite and infinite time horizons

We revisit and extend the Riccati theory, unifying continuous-time linear-quadratic optimal permanent and sampled-data control problems, in finite and infinite time horizons. In a nutshell, we prove that:-- when the time horizon T tends to $+\infty$, one passes from the Sampled-Data Difference Riccati Equation (SD-DRE) to the Sampled-Data Algebraic Riccati Equation (SD-ARE), and from the Permanent Differential Riccati Equation (P-DRE) to the Permanent Algebraic Riccati Equation (P-ARE);-- when the maximal step of the time partition $Δ$ tends to $0$, one passes from (SD-DRE) to (P-DRE), and from (SD-ARE) to (P-ARE).Our notations and analysis provide a unified framework in order to settle all corresponding results.

preprint2016arXiv

Asymptotic analysis and optimal control of an integro-differential system modelling healthy and cancer cells exposed to chemotherapy

We consider a system of two coupled integro-differential equations modelling populations of healthy and cancer cells under therapy. Both populations are structured by a phenotypic variable, representing their level of resistance to the treatment. We analyse the asymptotic behaviour of the model under constant infusion of drugs. By designing an appropriate Lyapunov function, we prove that both densities converge to Dirac masses. We then define an optimal control problem, by considering all possible infusion protocols and minimising the number of cancer cells over a prescribed time frame. We provide a quasi-optimal strategy and prove that it solves this problem for large final times. For this modelling framework, we illustrate our results with numerical simulations, and compare our optimal strategy with periodic treatment schedules.

preprint2016arXiv

Interaction Network, State Space and Control in Social Dynamics

In the present chapter we study the emergence of global patterns in large groups in first and second-order multi-agent systems, focusing on two ingredients that influence the dynamics: the interaction network and the state space. The state space determines the types of equilibrium that can be reached by the system. Meanwhile, convergence to specific equilibria depends on the connectivity of the interaction network and on the interaction potential. When the system does not satisfy the necessary conditions for convergence to the desired equilibrium, control can be exerted, both on finite-dimensional systems and on their mean-field limit.

preprint2016arXiv

Linear-quadratic optimal sampled-data control problems: convergence result and Riccati theory

We consider a general linear control system and a general quadratic cost, where the state evolves continuously in time and the control is sampled, i.e., is piecewise constant over a subdivision of the time interval. This is the framework of a linear-quadratic optimal sampled-data control problem. As a first result, we prove that, as the sampling periods tend to zero, the optimal sampled-data controls converge pointwise to the optimal permanent control. Then, we extend the classical Riccati theory to the sampled-data control framework, by developing two different approaches: the first one uses a recently established version of the Pontryagin maximum principle for optimal sampled-data control problems, and the second one uses an adequate version of the dynamic programming principle. In turn, we obtain a closed-loop expression for optimal sampled-data controls of linear-quadratic problems.

preprint2016arXiv

Low-Thrust Lyapunov to Lyapunov and Halo to Halo with $L^2$-Minimization

In this work, we develop a new method to design energy minimum low-thrust missions (L2-minimization). In the Circular Restricted Three Body Problem, the knowledge of invariant manifolds helps us initialize an indirect method solving a transfer mission between periodic Lyapunov orbits. Indeed, using the PMP, the optimal control problem is solved using Newton-like algorithms finding the zero of a shooting function. To compute a Lyapunov to Lyapunov mission, we first compute an admissible trajectory using a heteroclinic orbit between the two periodic orbits. It is then used to initialize a multiple shooting method in order to release the constraint. We finally optimize the terminal points on the periodic orbits. Moreover, we use continuation methods on position and on thrust, in order to gain robustness. A more general Halo to Halo mission, with different energies, is computed in the last section without heteroclinic orbits but using invariant manifolds to initialize shooting methods with a similar approach.

preprint2016arXiv

Minimal time spiking in various ChR2-controlled neuron models

We use conductance based neuron models and the mathematical modeling of Optogenetics to define controlled neuron models and we address the minimal time control of these affine systems for the first spike from equilibrium. We apply tools of geometric optimal control theory to study singular extremals and we implement a direct method to compute optimal controls. When the system is too large to theoretically investigate the existence of singular optimal controls, we observe numerically the optimal bang-bang controls.

preprint2016arXiv

Optimal control of infinite-dimensional piecewise deterministic Markov processes and application to the control of neuronal dynamics via Optogenetics

In this paper we define an infinite-dimensional controlled piecewise deterministic Markov process (PDMP) and we study an optimal control problem with finite time horizon and unbounded cost. This process is a coupling between a continuous time Markov Chain and a set of semilinear parabolic partial differential equations, both processes depending on the control. We apply dynamic programming to the embedded Markov decision process to obtain existence of optimal relaxed controls and we give some sufficient conditions ensuring the existence of an optimal ordinary control. This study, which constitutes an extension of controlled PDMPs to infinite dimension, is motivated by the control that provides Optogenetics on neuron models such as the Hodgkin-Huxley model. We define an infinite-dimensional controlled Hodgkin-Huxley model as an infinite-dimensional controlled piecewise deterministic Markov process and apply the previous results to prove the existence of optimal ordinary controls for a tracking problem.

preprint2016arXiv

Regularization of chattering phenomena via bounded variation control

In control theory, the term chattering is used to refer to strong oscillations of controls, such as an infinite number of switchings over a compact interval of times. In this paper we focus on three typical occurences of chattering: the Fuller phenomenon, referring to situations where an optimal control switches an infinite number of times over a compact set; the Robbins phenomenon, concerning optimal control problems with state constraints, meaning that the optimal trajectory touches the boundary of the constraint set an infinite number of times over a compact time interval; the Zeno phenomenon, referring as well to an infinite number of switchings over a compact set, for hybrid optimal control problems. From the practical point of view, when trying to compute an optimal trajectory, for instance by means of a shooting method, chattering may be a serious obstacle to convergence. In this paper we propose a general regularization procedure, by adding an appropriate penalization of the total variation. This produces a quasi-optimal control, and we prove that the family of quasi-optimal solutions converges to the optimal solution of the initial problem as the penalization tends to zero. Under additional assumptions, we also quantify the quasi-optimality property by determining a speed of convergence of the costs.

preprint2015arXiv

Coupled attitude and trajectory optimization for a launcher tilting maneuver

We study the minimum time control problem of the launchers. The optimal trajectories of the problem may contain singular arcs, and thus chattering arcs. The motion of the launcher is described by its attitude kinematics and dynamics and also by its trajectory dynamics. Based on the nature of the system, we implement an efficient indirect numerical method, combined with numerical continuation, to compute numerically the optimal solutions of the problem. Numerical experiments show the efficiency and the robustness of the proposed method.

preprint2015arXiv

Finite-dimensional predictor-based feedback stabilization of a 1D linear reaction-diffusion equation with boundary input delay

We consider a one-dimensional controlled reaction-diffusion equation, where the control acts on the boundary and is subject to a constant delay. Such a model is a paradigm for more general parabolic systems coupled with a transport equation. We prove that this is possible to stabilize (in H 1 norm) this process by means of an explicit predictor-based feedback control that is designed from a finite-dimensional subsystem. The implementation is very simple and efficient and is based on standard tools of pole-shifting. Our feedback acts on the system as a finite-dimensional predictor. We compare our approach with the backstepping method.

preprint2015arXiv

High order variational integrators in the optimal control of mechanical systems

In recent years, much effort in designing numerical methods for the simulation and optimization of mechanical systems has been put into schemes which are structure preserving. One particular class are variational integrators which are momentum preserving and symplectic. In this article, we develop two high order variational integrators which distinguish themselves in the dimension of the underling space of approximation and we investigate their application to finite-dimensional optimal control problems posed with mechanical systems. The convergence of state and control variables of the approximated problem is shown. Furthermore, by analyzing the adjoint systems of the optimal control problem and its discretized counterpart, we prove that, for these particular integrators, dualization and discretization commute.

preprint2015arXiv

Impulse and sampled-data optimal control of heat equations, and error estimates

We consider the optimal control problem of minimizing some quadratic functional over all possible solutions of an internally controlled multi-dimensional heat equation with a periodic terminal state constraint. This problem has a unique optimal solution, which can be characterized by an optimality system derived from the Pontryagin maximum principle. We define two approximations of this optimal control problem. The first one is an impulse approximation, and consists of considering a system of linear heat equations with impulse control. The second one is obtained by the sample-and-hold procedure applied to the control, resulting into a sampled-data approximation of the controlled heat equation. We prove that both problems have a unique optimal solution, and we establish precise error estimates for the optimal controls and optimal states of the initial problem with respect to its impulse and sampled-data approximations.

preprint2015arXiv

Minimum time control of the rocket attitude reorientation associated with orbit dynamics

In this paper, we investigate the minimal time problem for the guidance of a rocket, whose motion is described by its attitude kinematics and dynamics but also by its orbit dynamics. Our approach is based on a refined geometric study of the extremals coming from the application of the Pontryagin maximum principle. Our analysis reveals the existence of singular arcs of higher-order in the optimal synthesis, causing the occurrence of a chattering phenomenon, i.e., of an infinite number of switchings when trying to connect bang arcs with a singular arc. We establish a general result for bi-input control-affine systems, providing sufficient conditions under which the chattering phenomenon occurs. We show how this result can be applied to the problem of the guidance of the rocket. Based on this preliminary theoretical analysis, we implement efficient direct and indirect numerical methods, combined with numerical continuation, in order to compute numerically the optimal solutions of the problem.

preprint2015arXiv

Multiple Shape Registration using Constrained Optimal Control

Lagrangian particle formulations of the large deformation diffeomorphic metric mapping algorithm (LDDMM) only allow for the study of a single shape. In this paper, we introduce and discuss both a theoretical and practical setting for the simultaneous study of multiple shapes that are either stitched to one another or slide along a submanifold. The method is described within the optimal control formalism, and optimality conditions are given, together with the equations that are needed to implement augmented Lagrangian methods. Experimental results are provided for stitched and sliding surfaces.

preprint2015arXiv

Nonlinear damped partial differential equations and their uniform discretizations

We establish sharp energy decay rates for a large class of nonlinearly first-order damped systems, and we design discretization schemes that inherit of the same energy decay rates, uniformly with respect to the space and/or time discretization parameters, by adding appropriate numerical viscosity terms. Our main arguments use the optimal-weight convexity method and uniform observability inequalities with respect to the discretization parameters. We establish our results, first in the continuous setting, then for space semi-discrete models, and then for time semi-discrete models. The full discretization is inferred from the previous results. Our results cover, for instance, the Schrödinger equation with nonlinear damping, the nonlinear wave equation, the nonlinear plate equation, as well as certain classes of equations with nonlocal terms.

preprint2015arXiv

Optimal Neumann control for the 1D wave equation: finite horizon, infinite horizon, boundary tracking terms and the turnpike property

We consider a vibrating string that is fixed at one end with Neumann control action at the other end. We investigate the optimal control problem of steering this system from given initial data to rest, in time T , by minimizing an objective functional that is the convex sum of the L 2-norm of the control and of a boundary Neumann tracking term. We provide an explicit solution of this optimal control problem, showing that if the weight of the tracking term is positive, then the optimal control action is concentrated at the beginning and at the end of the time interval, and in-between it decays exponentially. We show that the optimal control can actually be written in that case as the sum of an exponentially decaying term and of an exponentially increasing term. This implies that, if the time T is large the optimal trajectory approximately consists of three arcs, where the first and the third short-time arcs are transient arcs, and in the middle arc the optimal control and the corresponding state are exponentially close to 0. This is an example for a turnpike phenomenon for a problem of optimal boundary control. If T = +$\infty$ (infinite horizon time problem), then only the exponentially decaying component of the control remains, and the norms of the optimal control action and of the optimal state decay exponentially in time. In contrast to this situation if the weight of the tracking term is zero and only the control cost is minimized, then the optimal control is distributed uniformly along the whole interval [0, T ] and coincides with the control given by the Hilbert Uniqueness Method.

preprint2015arXiv

Optimal sampled-data control, and generalizations on time scales

In this paper, we derive a version of the Pontryagin maximum principle for general finite-dimensional nonlinear optimal sampled-data control problems. Our framework is actually much more general, and we treat optimal control problems for which the state variable evolves on a given time scale (arbitrary non-empty closed subset of R), and the control variable evolves on a smaller time scale. Sampled-data systems are then a particular case. Our proof is based on the construction of appropriate needle-like variations and on the Ekeland variational principle.

preprint2015arXiv

Planar tilting maneuver of a spacecraft: singular arcs in the minimum time problem and chattering

In this paper, we study the minimum time planar tilting maneuver of a spacecraft, from the theoretical as well as from the numerical point of view, with a particular focus on the chattering phenomenon. We prove that there exist optimal chattering arcs when a singular junction occurs. Our study is based on the Pontryagin Maximum Principle and on results by M.I. Zelikin and V.F. Borisov. We give sufficient conditions on the initial values under which the optimal solutions do not contain any singular arc, and are bang-bang with a finite number of switchings. Moreover, we implement sub-optimal strategies by replacing the chattering control with a fixed number of piecewise constant controls. Numerical simulations illustrate our results.

preprint2015arXiv

Pontryagin maximum principle for optimal sampled-data control problems

In this short communication, we first recall a version of the Pontryagin maximum principle for general finite-dimensional nonlinear optimal sampled-data control problems. This result was recently obtained in [L. Bourdin and E. Tr{é}lat , Optimal sampled-data control, and generalizations on time scales, arXiv:1501.07361, 2015]. Then we discuss the maximization condition for optimal sampled-data controls that can be seen as an average of the weak maximization condition stated in the classical Pontryagin maximum principle for optimal (permanent) controls. Finally, applying this theorem, we solve a linear-quadratic example based on the classical parking problem.

preprint2015arXiv

Quantum ergodicity and quantum limits for sub-Riemannian Laplacians

This paper is a proceedings version of \cite{CHT-I}, in which we state a Quantum Ergodicity (QE) theorem on a 3D contact manifold, and in which we establish some properties of the Quantum Limits (QL). We consider a sub-Riemannian (sR) metric on a compact 3D manifold with an oriented contact distribution. There exists a privileged choice of the contact form, with an associated Reeb vector field and a canonical volume form that coincides with the Popp measure. We state a QE theorem for the eigenfunctions of any associated sR Laplacian, under the assumption that the Reeb flow is ergodic. The limit measure is given by the normalized canonical contact measure. To our knowledge, this is the first extension of the usual Schnirelman theorem to a hypoelliptic operator. We provide as well a decomposition result of QL's, which is valid without any ergodicity assumption. We explain the main steps of the proof, and we discuss possible extensions to other sR geometries.

preprint2015arXiv

Stabilization of semilinear PDE's, and uniform decay under discretization

These notes are issued from a short course given by the author in a summer school in Chamb{é}ry in June 2015. We consider general semilinear PDE's and we address the following two questions: 1) How to design an efficient feedback control locally stabilizing the equation asymptotically to 0? 2) How to construct such a stabilizing feedback from approximation schemes? To address these issues, we distinguish between parabolic and hyperbolic semilinear PDE's. By parabolic, we mean that the linear operator underlying the system generates an analytic semi-group. By hyperbolic, we mean that this operator is skew-adjoint. We first recall general results allowing one to consider the nonlinear term as a perturbation that can be absorbed when one is able to construct a Lyapunov function for the linear part. We recall in particular some known results borrowed from the Riccati theory. However, since the numerical implementation of Riccati operators is computationally demanding, we focus then on the question of being able to design "simple" feedbacks. For parabolic equations, we describe a method consisting of designing a stabilizing feedback, based on a small finite-dimensional (spectral) approximation of the whole system. For hyperbolic equations, we focus on simple linear or nonlinear feedbacks and we investigate the question of obtaining sharp decay results. When considering discretization schemes, the decay obtained in the continuous model cannot in general be preserved for the discrete model, and we address the question of adding appropriate viscosity terms in the numerical scheme, in order to recover a uniform decay. We consider space, time and then full discretizations and we report in particular on the most recent results obtained in the literature. Finally, we describe several open problems and issues.

preprint2015arXiv

Sub-Riemannian structures on groups of diffeomorphisms

In this paper, we define and study strong right-invariant sub-Riemannian structures on the group of diffeomorphisms of a manifold with bounded geometry. We derive the Hamiltonian geodesic equations for such structures, and we provide examples of normal and of abnormal geodesics in that infinite-dimensional context. The momentum formulation gives a sub-Riemannian version of the Euler-Arnol'd equation. Finally, we establish some approximate and exact reachability properties for diffeomorphisms, and we give some consequences for Moser theorems.

preprint2014arXiv

Control to flocking of the kinetic Cucker-Smale model

The well-known Cucker-Smale model is a macroscopic system reflecting flocking, i.e. the alignment of velocities in a group of autonomous agents having mutual interactions. In the present paper, we consider the mean-field limit of that model, called the kinetic Cucker-Smale model, which is a transport partial differential equation involving nonlocal terms. It is known that flocking is reached asymptotically whenever the initial conditions of the group of agents are in a favorable configuration. For other initial configurations, it is natural to investigate whether flocking can be enforced by means of an appropriate external force, applied to an adequate time-varying subdomain. In this paper we prove that we can drive to flocking any group of agents governed by the kinetic Cucker-Smale model, by means of a sparse centralized control strategy, and this, for any initial configuration of the crowd. Here, "sparse control" means that the action at each time is limited over an arbitrary proportion of the crowd, or, as a variant, of the space of configurations; "centralized" means that the strategy is computed by an external agent knowing the configuration of all agents. We stress that we do not only design a control function (in a sampled feedback form), but also a time-varying control domain on which the action is applied. The sparsity constraint reflects the fact that one cannot act on the whole crowd at every instant of time. Our approach is based on geometric considerations on the velocity field of the kinetic Cucker-Smale PDE, and in particular on the analysis of the particle flow generated by this vector field. The control domain and the control functions are designed to satisfy appropriate constraints, and such that, for any initial configuration, the velocity part of the support of the measure solution asymptotically shrinks to a singleton, which means flocking.

preprint2014arXiv

Optimal shape and location of sensors for parabolic equations with random initial data

In this article, we consider parabolic equations on a bounded open connected subset $Ω$ of $\R^n$. We model and investigate the problem of optimal shape and location of the observation domain having a prescribed measure. This problem is motivated by the question of knowing how to shape and place sensors in some domain in order to maximize the quality of the observation: for instance, what is the optimal location and shape of a thermometer? We show that it is relevant to consider a spectral optimal design problem corresponding to an average of the classical observability inequality over random initial data, where the unknown ranges over the set of all possible measurable subsets of $Ω$ of fixed measure. We prove that, under appropriate sufficient spectral assumptions, this optimal design problem has a unique solution, depending only on a finite number of modes, and that the optimal domain is semi-analytic and thus has a finite number of connected components. This result is in strong contrast with hyperbolic conservative equations (wave and Schrödinger) studied in [56] for which relaxation does occur. We also provide examples of applications to anomalous diffusion or to the Stokes equations. In the case where the underlying operator is any positive (possible fractional) power of the negative of the Dirichlet-Laplacian, we show that, surprisingly enough, the complexity of the optimal domain may strongly depend on both the geometry of the domain and on the positive power. The results are illustrated with several numerical simulations.

preprint2014arXiv

Shape deformation analysis from the optimal control viewpoint

A crucial problem in shape deformation analysis is to determine a deformation of a given shape into another one, which is optimal for a certain cost. It has a number of applications in particular in medical imaging. In this article we provide a new general approach to shape deformation analysis, within the framework of optimal control theory, in which a deformation is represented as the flow of diffeomorphisms generated by time-dependent vector fields. Using reproducing kernel Hilbert spaces of vector fields, the general shape deformation analysis problem is specified as an infinite-dimensional optimal control problem with state and control constraints. In this problem, the states are diffeomorphisms and the controls are vector fields, both of them being subject to some constraints. The functional to be minimized is the sum of a first term defined as geometric norm of the control (kinetic energy of the deformation) and of a data attachment term providing a geometric distance to the target shape. This point of view has several advantages. First, it allows one to model general constrained shape analysis problems, which opens new issues in this field. Second, using an extension of the Pontryagin maximum principle, one can characterize the optimal solutions of the shape deformation problem in a very general way as the solutions of constrained geodesic equations. Finally, recasting general algorithms of optimal control into shape analysis yields new efficient numerical methods in shape deformation analysis. Overall, the optimal control point of view unifies and generalizes different theoretical and numerical approaches to shape deformation problems, and also allows us to design new approaches. The optimal control problems that result from this construction are infinite dimensional and involve some constraints, and thus are nonstandard. In this article we also provide a rigorous and complete analysis of the infinite-dimensional shape space problem with constraints and of its finite-dimensional approximations.

preprint2014arXiv

Sparse Stabilization and Control of Alignment Models

From a mathematical point of view self-organization can be described as patterns to which certain dynamical systems modeling social dynamics tend spontaneously to be attracted. In this paper we explore situations beyond self-organization, in particular how to externally control such dynamical systems in order to eventually enforce pattern formation also in those situations where this wished phenomenon does not result from spontaneous convergence. Our focus is on dynamical systems of Cucker-Smale type, modeling consensus emergence, and we question the existence of stabilization and optimal control strategies which require the minimal amount of external intervention for nevertheless inducing consensus in a group of interacting agents. We provide a variational criterion to explicitly design feedback controls that are componentwise sparse, i.e. with at most one nonzero component at every instant of time. Controls sharing this sparsity feature are very realistic and convenient for practical issues. Moreover, the maximally sparse ones are instantaneously optimal in terms of the decay rate of a suitably designed Lyapunov functional, measuring the distance from consensus. As a consequence we provide a mathematical justification to the general principle according to which "sparse is better" in the sense that a policy maker, who is not allowed to predict future developments, should always consider more favorable to intervene with stronger action on the fewest possible instantaneous optimal leaders rather than trying to control more agents with minor strength in order to achieve group consensus. We then establish local and global sparse controllability properties to consensus and, finally, we analyze the sparsity of solutions of the finite time optimal control problem where the minimization criterion is a combination of the distance from consensus and of the l1-norm of the control.

preprint2014arXiv

The turnpike property in finite-dimensional nonlinear optimal control

Turnpike properties have been established long time ago in finite-dimensional optimal control problems arising in econometry. They refer to the fact that, under quite general assumptions, the optimal solutions of a given optimal control problem settled in large time consist approximately of three pieces, the first and the last of which being transient short-time arcs, and the middle piece being a long-time arc staying exponentially close to the optimal steady-state solution of an associated static optimal control problem. We provide in this paper a general version of a turnpike theorem, valuable for nonlinear dynamics without any specific assumption, and for very general terminal conditions. Not only the optimal trajectory is shown to remain exponentially close to a steady-state, but also the corresponding adjoint vector of the Pontryagin maximum principle. The exponential closedness is quantified with the use of appropriate normal forms of Riccati equations. We show then how the property on the adjoint vector can be adequately used in order to initialize successfully a numerical direct method, or a shooting method. In particular, we provide an appropriate variant of the usual shooting method in which we initialize the adjoint vector, not at the initial time, but at the middle of the trajectory.

preprint2013arXiv

Pontryagin Maximum Principle for finite dimensional nonlinear optimal control problems on time scales

In this article we derive a strong version of the Pontryagin Maximum Principle for general nonlinear optimal control problems on time scales in finite dimension. The final time can be fixed or not, and in the case of general boundary conditions we derive the corresponding transversality conditions. Our proof is based on Ekeland's variational principle. Our statement and comments clearly show the distinction between right-dense points and right-scattered points. At right-dense points a maximization condition of the Hamiltonian is derived, similarly to the continuous-time case. At right-scattered points a weaker condition is derived, in terms of so-called stable $Ω$-dense directions. We do not make any specific restrictive assumption on the dynamics or on the set $Ω$ of control constraints. Our statement encompasses the classical continuous-time and discrete-time versions of the Pontryagin Maximum Principle, and holds on any general time scale, that is any closed subset of $\R$.

preprint2012arXiv

General Cauchy-Lipschitz theory for shifted and non shifted Delta-Cauchy problems on time scales

This article is devoted to completing some aspects of the classical Cauchy-Lipschitz (or Picard-Lindelöf) theory for general nonlinear systems posed on time scales, that are closed subsets of the set of real numbers. Partial results do exist but do not cover the framework of general dynamics on time-scales encountered e.g. in applications to control theory. In the present work, we first introduce the notion of absolutely continuous solution for shifted and non shifted Delta-Cauchy problems, and then the notion of a maximal solution. We state and prove a Cauchy-Lipschitz theorem, providing existence and uniqueness of the maximal solution of a given Delta-Cauchy problem under suitable assumptions like regressivity and local Lipschitz continuity, and discuss some related issues like the behavior of maximal solutions at terminal points.

preprint2012arXiv

Optimal observability of the multi-dimensional wave and Schrödinger equations in quantum ergodic domains

We consider the wave and Schrödinger equations on a bounded open connected subset $Ω$ of a Riemannian manifold, with Dirichlet, Neumann or Robin boundary conditions whenever its boundary is nonempty. We observe the restriction of the solutions to a measurable subset $ω$ of $Ω$ during a time interval [0,T] with T>0. We address the problem of the optimal location of the observation subset ωamong all possible subsets of a given measure or volume fraction. We solve it in two different situations. First, when a specific choice of the initial data is given and therefore we deal with a particular solution, we show that the problem always admits at least one solution that can be regular or of fractal type depending on the regularity of the initial data. We then introduce a relevant criterion, not depending on the initial conditions and to choose the observation set in an uniform way, independent of the data and solutions under consideration. Through spectral decompositions, this leads to a second problem which consists of maximizing a spectral functional that can be viewed as a measure of eigenfunction concentration. This spectral criterion can be interpreted in two ways: it corresponds either to a time asymptotic observability constant as the observation time interval tends to infinity, or to a randomized version of the deterministic observability inequality. Under appropriate quantum ergodicity assumptions on $Ω$, we compute the optimal value. We then provide spectral approximations and present some numerical simulations that fully confirm the theoretical results in the paper and support our conjectures. Our results highlight precise connections between optimal observability issues and quantum ergodic properties of the domain under consideration.

preprint2007arXiv

Nonlinear optimal control via occupation measures and LMI-relaxations

We consider the class of nonlinear optimal control problems (OCP) with polynomial data, i.e., the differential equation, state and control con- straints and cost are all described by polynomials, and more generally for OCPs with smooth data. In addition, state constraints as well as state and/or action constraints are allowed. We provide a simple hierarchy of LMI (lin- ear matrix inequality)-relaxations whose optimal values form a nondecreasing sequence of lower bounds on the optimal value. Under some convexity assump- tions, the sequence converges to the optimal value of the OCP. Preliminary results show that good approximations are obtained with few moments.

preprint2006arXiv

Asymptotics of accessibility sets along an abnormal trajectory

We describe precisely, under generic conditions, the contact of the accessibility set at time $T$ with an abnormal direction, first for a single-input affine control system with constraint on the control, and then as an application for a sub-Riemannian system of rank 2. As a consequence we obtain in sub-Riemannian geometry a new splitting-up of the sphere near an abnormal minimizer $γ$ into two sectors, bordered by the first Pontryagin's cone along $γ$, called the $\xLinfty$-sector and the $\xLtwo$-sector. Moreover we find again necessary and sufficient conditions of optimality of an abnormal trajectory for such systems, for any optimization problem.

preprint2006arXiv

Control of travelling walls in a ferromagnetic nanowire

We investigate the problem of controlling the magnetic moment in a ferromagnetic nanowire submitted to an external magnetic field in the direction of the nanowire. The system is modeled with the one dimensional Landau-Lifschitz equation. In the absence of control, there exist particular solutions, which happen to be relevant for practical issues, called travelling walls. In this paper, we prove that it is possible to move from a given travelling wall profile to any other one, by acting on the external magnetic field. The control laws are simple and explicit, and the resulting trajectories are shown to be stable.

preprint2006arXiv

Quasi-optimal robust stabilization of control systems

In this paper, we investigate the problem of semi-global minimal time robust stabilization of analytic control systems with controls entering linearly, by means of a hybrid state feedback law. It is shown that, in the absence of minimal time singular trajectories, the solutions of the closed-loop system converge to the origin in quasi minimal time (for a given bound on the controller) with a robustness property with respect to small measurement noise, external disturbances and actuator noise.

preprint2006arXiv

Singular trajectories of control-affine systems

When applying methods of optimal control to motion planning or stabilization problems, some theoretical or numerical difficulties may arise, due to the presence of specific trajectories, namely, singular minimizing trajectories of the underlying optimal control problem. In this article, we provide characterizations for singular trajectories of control-affine systems. We prove that, under generic assumptions, such trajectories share nice properties, related to computational aspects; more precisely, we show that, for a generic system -- with respect to the Whitney topology --, all nontrivial singular trajectories are of minimal order and of corank one. These results, established both for driftless and for control-affine systems, extend previous results. As a consequence, for generic systems having more than two vector fields, and for a fixed cost, there do not exist minimizing singular trajectories. We also prove that, given a control system satisfying the LARC, singular trajectories are strictly abnormal, generically with respect to the cost. We then show how these results can be used to derive regularity results for the value function and in the theory of Hamilton-Jacobi equations, which in turn have applications for stabilization and motion planning, both from the theoretical and implementation issues.

preprint2006arXiv

Some properties of the value function and its level sets for affine control systems with quadratic cost

Let $T>0$ fixed. We consider the optimal control problem for analytic affine systems: $\ds{\dot{x}=f\_0(x)+\sum\_{i=1}^m u\_if\_i(x)}$, with a cost of the form: $\ds{C(u)=\int\_0^T \sum\_{i=1}^m u\_i^2(t)dt}$. For this kind of systems we prove that if there are no minimizing abnormal extremals then the value function $S$ is subanalytic. Secondly we prove that if there exists an abnormal minimizer of corank 1 then the set of end-points of minimizers at cost fixed is tangent to a given hyperplane. We illustrate this situation in sub-Riemannian geometry.