La distribution de marchandises est un secteur en constant développement et constitue un facteur économique important. Par contre, dans les villes, il contribue notamment aux problèmes de congestion, pollution, bruit et d'autres dérangements à la population des villes. Pour faire face à ces problèmes, une nouvelle discipline est née à la fin du XXe siècle, la " City Logistics ", qui a comme objectifs principaux la réduction de la congestion, la pollution et le bruit occasionné par le transport de marchandises en ville. Dans les dernières années, plusieurs études et expériences se sont développées en toute l'Europe, mais pour l'instant une politique commune en matière de logistique urbaine n'a pas encore été proposée par l'Union Européenne. En Italie, seulement certaines villes de petite taille ont expérimenté des politiques de " city logistics " avec succès, mais sans un lien entre elles. Nous observons que ces expériences utilisent des centres urbains de distribution de marchandises, ce qui peut se traduire en un système de transport à deux ou plus niveaux. Plusieurs études en recherche opérationnelle ont traité des problématiques liées à des systèmes à niveaux multiples pour la distribution de marchandise. Néanmoins, l'optimisation des coûts de transport est en générale réalisé en considérant chaque niveau indépendant des autres, ou en approximant les coûts du transport dans certains niveaux pour simplifier. Un autre problème est le manque d'une unification de la terminologie utilisée dans ces études, qui difficulte la recherche bibliographique. Le but de cette recherche est, d'un coté, proposer des lignes guide d'accion en matière de planification de la distribution urbaine de marchandises, en unifiant certains termes, et d'un autre coté présenter une famille de problèmes d'optimisation de routes des véhicules qui considère les systèmes à niveaux multiples dans son ensemble et pas comme une somme de systèmes indépendants. Dans un premier temps, nous présentons les principales expériences de " city logistics " en Italie, ainsi que des lignes d'action dans la planification des systèmes de distribution urbaine des marchandises qui puissent devenir opérationnels et efficients. Ensuite nous présentons les principales problématiques et limites de l'optimisation de systèmes de transports à niveaux multiples, en unifiant les concepts et la notation. Nous proposons une nouvelle famille de problèmes d'optimisation de routes de véhicules pour des systèmes à niveaux multiples, en détaillant le cas basique : le problème de routes de véhicules à deux niveaux. Nous proposons des modèles mathématiques pour ce problème et des résultats numériques pour illustrer les avantages et les limites de la modélisation de ces systèmes.

Les Algorithmes Évolutionnaires (AEs) ont été très étudiés en raison de leur capacité à résoudre des problèmes d'optimisation complexes en utilisant des opérateurs de variation adaptés à des problèmes spécifiques. Une recherche dirigée par une population de solutions offre une bonne robustesse par rapport à un bruit modéré et la multi-modalité de la fonction optimisée, contrairement à d'autres méthodes d'optimisation classiques telles que les méthodes de quasi-Newton. La principale limitation de AEs, le grand nombre d'évaluations de la fonction objectif, pénalise toutefois l'usage des AEs pour l'optimisation de fonctions chères en temps calcul. La présente thèse se concentre sur un algorithme évolutionnaire, Covariance Matrix Adaptation Evolution Strategy (CMA-ES), connu comme un algorithme puissant pour l'optimisation continue boîte noire. Nous présentons l'état de l'art des algorithmes, dérivés de CMA-ES, pour résoudre les problèmes d'optimisation mono- et multi-objectifs dans le scénario boîte noire. Une première contribution, visant l'optimisation de fonctions coûteuses, concerne l'approximation scalaire de la fonction objectif. Le meta-modèle appris respecte l'ordre des solutions (induit par la valeur de la fonction objectif pour ces solutions) ; il est ainsi invariant par transformation monotone de la fonction objectif. L'algorithme ainsi défini, saACM-ES, intègre étroitement l'optimisation réalisée par CMA-ES et l'apprentissage statistique de meta-modèles adaptatifs ; en particulier les meta-modèles reposent sur la matrice de covariance adaptée par CMA-ES. saACM-ES préserve ainsi les deux propriété clé d'invariance de CMA-ES~: invariance i) par rapport aux transformations monotones de la fonction objectif; et ii) par rapport aux transformations orthogonales de l'espace de recherche. L'approche est étendue au cadre de l'optimisation multi-objectifs, en proposant deux types de meta-modèles (scalaires). La première repose sur la caractérisation du front de Pareto courant (utilisant une variante mixte de One Class Support Vector Machone (SVM) pour les points dominés et de Regression SVM pour les points non-dominés). La seconde repose sur l'apprentissage d'ordre des solutions (rang de Pareto) des solutions. Ces deux approches sont intégrées à CMA-ES pour l'optimisation multi-objectif (MO-CMA-ES) et nous discutons quelques aspects de l'exploitation de meta-modèles dans le contexte de l'optimisation multi-objectif. Une seconde contribution concerne la conception d'algorithmes nouveaux pour l'optimi\-sation mono-objectif, multi-objectifs et multi-modale, développés pour comprendre, explorer et élargir les frontières du domaine des algorithmes évolutionnaires et CMA-ES en particulier. Spécifiquement, l'adaptation du système de coordonnées proposée par CMA-ES est couplée à une méthode adaptative de descente coordonnée par coordonnée. Une stratégie adaptative de redémarrage de CMA-ES est proposée pour l'optimisation multi-modale. Enfin, des stratégies de sélection adaptées aux cas de l'optimisation multi-objectifs et remédiant aux difficultés rencontrées par MO-CMA-ES sont proposées.

Cette thèse a pour cadre le développement de méthodes numériques pour la conception optimale de forme de voilure en aérodynamique compressible. Ce travail a donné lieu au développement et à la validation d'un adjoint discret aéro-structure pour l'analyse de sensibilité par rapport aux paramètres de forme et de structure interne de l 'aile dont dépend la fonction d'intérêt, qu' elle soit aérodynamique ou structurale. Les développements logiciels ont été réalisés dans le code de simulation numérique de mécanique des fluides elsA et font suite aux travaux de Marcelet portant sur l'adjoint aéroélastique et dont ils sont une extension. Alors que pour l'adjoint aéro-élastique, on considère une aile flexible, de caractéristiques structurales constantes, pour l'adjoint aérostructure, leur variations sont prises en compte. Pour cela, l'extension de la méthode adjointe s' est accompagnée du développement d' un module de modélisation de la structure interne de l'aile. Ce module est linéarisé et vient donc alimenter le système adjoint. Il a été validé par le dimensionnement de la structure primaire d'une configuration de recherche fournie par l' avionneur Airbus. Dans l'état actuel de développement de la méthode adjointe dans le code elsA, on peut donc désormais calculer les sensibilités d'une fonction aérodynamique ou d'une fonction structure par rapport aux paramètres de forme aérodynamique ou de structure interne de l'aile. Le calcul des gradients ainsi obtenus a été validé par des comparaisons systématiques aux différences finies.

Cette thèse présente des algorithmes permettant la parallélisation dans la direction temporelle de la simulation des systèmes régis par des équations aux dérivées partielles issues de trois domaines d'application proposant des modèles complexes: \textbf{1. Contrôle optimal d'équations paraboliques:}\\ Nous développons deux algorithmes parallèles (SITPOC et PITPOC). Ces deux algorithmes sont basés sur une méthode générale de décomposition en temps des problèmes de contrôle optimal. Un résultat de convergence est obtenu pour SITPOC ainsi que des interprétations matricielles des deux algorithmes. \textbf{2. Cinétique de la population de neutrons dans un réacteur nucléaire:}\\ Nous concevons un solveur parallélisé en temps regroupant toutes les variables issues de la modélisation neutronique d'un réacteur. Notre méthode est adaptée aux différents scénarios possibles réduisant la physique. Les résultats numériques sont comparables à ceux obtenus par le code MINOS du CEA. La parallélisation repose sur un schéma de type pararéel pour la résolution en temps. Nous considérons plusieurs modèles physiques de la cinétique des neutrons que nous traitons par notre algorithme dans lequel le solveur grossier utilisé correspond à une simulation par réduction du modèle physique. Cette réduction est bénéfique puisqu'elle permet une accélération importante du traitement en temps machine. \textbf{3. Contrôle par résonance magnétique nucléaire et information quantique:}\\ Ce chapitre présente un travail d'ouverture sur une méthode parallèle en temps pour la résolution d'un problème de contrôle optimal en résonance magnétique nucléaire. Notre méthode produit une accélération importante par rapport à l'algorithme de référence non parallèle. De plus, les champs de contrôle calculés par notre méthode sont lisses, ce qui permet une mise en œuvre expérimentale plus simple dans les instruments. Les tests numériques prouvent l'efficacité de notre approche. Sur des exemples académiques et sans faire usage de techniques de programmation avancées, nous obtenons des accélérations de résolution significatives.

Cette thèse s'intéresse au problème d'estimation et d'optimisation distribuée dans les réseaux asynchrones, c'est à dire en n'utilisant que des communications locales et asynchrones. A partir de multiples applications allant de l'apprentissage automatique aux réseaux de capteurs sans-fils, nous concevons et analysons théoriquement de nouveaux algorithmes résolvant trois problèmes de nature très différentes : la propagation de la plus grande des valeurs initiales, l'estimation de leur moyenne et enfin l'optimisation distribuée.

Ce mémoire présente les travaux réalisés au cours de ma thèse sur le contrôle d'équations aux dérivées partielles. La première partie est consacrée à l'étude d'équations de Schrödinger bilinéaires unidimensionnelles autour de deux axes : la non contrôlabilité en temps petit avec des contrôles petits et la contrôlabilité simultanée. On établit un cadre pour lequel, bien que la vitesse de propagation du système soit infinie, la contrôlabilité exacte locale avec des contrôles petits est vérifiée si et seulement si le temps est suffisamment grand. Ces résultats, basés sur la coercivité d'une forme quadratique associée au développement de l'état à l'ordre deux, sont étendus dans le contexte de la contrôlabilité simultanée. On montre alors, en utilisant la méthode du retour de J.-M.~Coron, des résultats de contrôle exact local simultané pour deux ou trois équations, à phase globale et/ou à retard global près. La trajectoire de référence utilisée est construite via des résultats de contrôle partiel. En utilisant un argument de perturbation, on étend cette idée pour montrer la contrôlabilité exacte globale d'un nombre quelconque d'équations sans hypothèse sur le potentiel. Dans la deuxième partie, on prend en compte dans le modèle un terme supplémentaire quadratique en le contrôle. Ce terme, dit de polarisabilité, généralement négligé, présente un intérêt physique dans la modélisation, mais aussi mathématique dans le cas où le terme bilinéaire est insuffisant pour conclure à la contrôlabilité. En dimension quelconque, on construit des contrôles explicites réalisant la contrôlabilité approchée de l'état fondamental. En adaptant conjointement l'argument de perturbation précédent et certains résultats du contrôle bilinéaire, on prouve la contrôlabilité globale exacte de l'équation de Schrödinger avec polarisabilité 1D. La dernière partie de ce mémoire est consacrée à l'étude de la continuation unique pour un opérateur de type Grushin sur un rectangle 2D. Cet opérateur présente une singularité et une dégénérescence sur un segment séparant le domaine en deux composantes. On donne une condition nécessaire et suffisante sur le coefficient du potentiel singulier pour obtenir la continuation unique.

Cette these presente les techniques d'optimisation combinatoire et de programmation entiere transversales a nos differents projets de recherche.

Cette thèse présente une étude des systèmes dynamiques polynomiaux motivée à la fois par le grand spectre d'applications de cette classe (modèles de réactions chimiques, modèles de circuits électriques ainsi que les modèles biologiques) et par la difficulté (voire incapacité) de la résolution théorique de tels systèmes. Dans une première partie préliminaire, nous présentons les polynômes multivariés et nous introduisons les notions de forme polaire d'un polynôme (floraison) et de polynômes de Bernstein qui seront d'un grand intérêt par la suite. Dans une deuxième partie, nous considérons le problème d'optimisation polynomial dit POP. Nous décrivons dans un premier temps les principales méthodes existantes permettant de résoudre ou d'approcher la solution d'un tel problème. Puis, nous présentons deux relaxations linéaires se basant respectivement sur le principe de floraison ainsi que les polynômes de Bernstein permettant d'approcher la valeur optimale du POP. La dernière partie de la thèse sera consacrée aux applications de nos deux méthodes de relaxation dans le cadre des systèmes dynamiques polynomiaux. Une première application s'inscrit dans le cadre de l'analyse d'atteignabilité : en effet, on utilisera notre relaxation de Bernstein pour pouvoir construire un algorithme permettant d'approximer les ensembles atteignables d'un système dynamique polynomial discrétisé. Une deuxième application sera la vérification et le calcul d'invariants pour un système dynamique polynomial. Une troisième application consiste à calculer un contrôleur et un invariant pour un système dynamique polynomial soumis à des perturbations. Dans le contexte de l'invariance, on utilisera la relaxation se basant sur le principe de floraison. Enfin, une dernière application sera d'exploiter les principales propriétés de la forme polaire pour pouvoir étudier des systèmes dynamiques polynomiaux dans des rectangles.

La gestion du trafic aérien est un système complexe. Actuellement en pleine mutation, une des problématiques essentielles à l'évolution du système est la recherche de méthodes automatiques de résolution de conflits. Nous présentons d'abord un nouveau modèle de trajectoire courbe basé sur les B-splines et permettant de définir une trajectoire à l'aide d'un nombre très limité de paramètres. À partir de cette modélisation, nous arrêtons une nouvelle formulation du problème de résolution de conflits pour obtenir un problème d'optimisation continue. Celle-ci repose sur une formulation dite semi-infinie de la contrainte de séparation entre deux avions. La manière dont nous avons défini la fonction-objectif et les fonctions contraintes nous permettent également d'en calculer les gradients. Nous utilisons trois différentes méthodes d'optimisation pour résoudre notre problème. Une méthode globale stochastique est d'abord testée : les algorithmes génétiques, couramment utilisés pour le problème de résolution de conflits. Deux méthodes d'optimisation locale sont aussi mises en œuvre, une méthode de points intérieurs et une méthode d'optimisation sans dérivées. Enfin, nous présentons des résultats numériques prometteurs montrant la fiabilité de l'optimisation locale pour le problème de résolution de conflits. Notre méthodologie, alliant une modèle de trajectoire courbe parcimonieux et une méthode d'optimisation locale appliquée à notre formulation mathématique du problème, est une option crédible pour le problème de résolution de conflits aériens.

Les moteurs de recherche jouent un rôle essentiel sur le Web. Ils rassemblent des informations sur les pages web et pour chaque requête d'un internaute, ils donnent une liste ordonnée de pages pertinentes. Ils utilisent divers algorithmes pour classer les pages en fonction de leur contenu textuel ou de la structure d'hyperlien du Web. Ici, nous nous concentrons sur les algorithmes qui utilisent cette structure d'hyperliens, comme le PageRank, HITS, SALSA et HOTS. La notion fondamentale pour tous ces algorithmes est le graphe du web. C'est le graphe orienté qui a un noeud pour chaque page web et un arc entre les noeuds i et j si il y a un hyperlien entre les pages i et j. Le problème original considéré dans cette thèse est l'optimisation du référencement des pages d'un site web donné. Il consiste à trouver une stratégie optimale de liens qui maximise une fonction scalaire d'un classement donné sous des contraintes de design. Le PageRank, HITS et SALSA classent les pages par un vecteur de Perron, c'est-à-dire qu'ils correspondent au vecteur propre principal d'une matrice à coefficients positifs. Quand on optimise le référencement, on optimise donc une fonction scalaire du vecteur propre de Perron sur un ensemble de matrices positives irréductibles. La matrice est construite à partir du graphe du web, donc commander les hyperliens revient à commander la matrice elle-même. Nous étudions d'abord un problème général d'optimisation du PageRank avec une fonction d'utilité correspondant au revenu total du site et des contraintes de design. Ce cas est d'un intérêt particulier car pour de nombreux sites la valeur du PageRank est corrélée au chiffre d'affaires. Nous avons réduit le problème d'optimisation du PageRank à des problèmes de décision markoviens dont les ensembles d'action sont définis implicitement comme étant les points extrêmes de polytopes qui ont un oracle de séparation polynomial. Nous montrons que de tels problèmes de décision markoviens sont solubles en temps polynomial et nous donnons un algorithme qui passe à l'échelle pour la résolution effective du problème d'optimisation du PageRank sur de grandes bases de données. Ensuite, nous étudions le problème général de l'optimisation d'une fonction scalaire du vecteur propre de Perron sur un ensemble de matrices positives irréductibles. Cela couvre tous les classements par vecteur de Perron, HITS et SALSA compris. Nous montrons que la matrice des dérivées partielles de la fonction objectif a un petit rang et peut être calculée par un algorithme qui a les mêmes propriétés de convergence que la méthode de la puissance utilisée pour calculer le classement. Nous donnons un algorithme d'optimisation qui couple les itérations puissance et gradient et nous prouvons sa convergence vers un point stationnaire du problème d'optimisation. En considérant HOTS comme un vecteur de Perron non linéaire, nous montrons que l'algorithme HOTS converge géométriquement et nous résolvons l'optimisation locale de HOTS. Finalement, nous étendons le domaine d'application des méthodes d'optimisation du vecteur propre et de la valeur propre de Perron à l'optimisation de la chimiothérapie, sous l'hypothèse que les cellules se comportent différemment suivant l'heure de la journée. Nous voulons profiter de cette caractéristique pour minimiser le taux de croissance des cellules cancéreuses tout en maintenant le taux de croissance des cellules saines au dessus d'un seuil de toxicité donné. L'objectif et la contrainte peuvent s'écrire comme les valeurs propres de Floquet de modèles d'EDP structurés en âge avec des coefficients périodiques, qui sont approchés par des valeurs propres de Perron dans le problème discrétisé. Nous cherchons des stratégies d'injection de médicament localement optimales par une méthode des multiplicateurs où les minimisations sans contrainte sont faites en utilisant l'algorithme couplant les itérations puissance et gradient développé dans le cadre de l'optimisation du référencement.