Les méthodes barrières proposent de résoudre le problème non linéaire en résolvant une suite de problèmes pénalisés. Le lien entre la suite, dite externe, des solutions des fonctions pénalisées et la solution du problème initial a été établie dans les années soixante.<br /><br /> Dans cette thèse, nous avons utilisé une fonction barrière logarithmique. A chaque itération externe, la technique SQP se charge de produire une série de sous-problèmes quadratiques dont les solutions forment une suite, dite interne, de directions de descente pour résoudre le problème non linéaire pénalisé.<br /><br /> Nous avons introduit un changement de variable sur le pas de déplacement ce qui a permis d'obtenir des conditions d'optimalité plus stable numériquement.<br /><br /> Nous avons réalisé des simulations numériques pour comparer les performances de la méthode des gradients conjugués à celle de la méthode D.C., appliquées pour résoudre des problèmes quadratiques de région de confiance.<br /><br /> Nous avons adapté la méthode D.C. pour résoudre les sous-problèmes verticaux, ce qui nous a permis de ramener leurs dimensions de $n+m$ à $m+p$ ($ p < n $).<br /><br /> L'évolution de l'algorithme est contrôlée par la fonction de mérite. Des tests numériques permettent de comparer les avantages de différentes formes de la fonction de mérite. Nous avons introduit de nouvelles règles pour améliorer cette évolution.<br /><br /> Les expériences numériques montrent un gain concernant le nombre de problèmes résolus. L'étude de la convergence de notre méthode SDC, clôt ce travail.

Cette thèse propose des techniques numériques visant à optimiser les méthodes itératives d'évaluation de performances de modèles Markoviens. Ces techniques s'appliquent à des modèles où la matrice de transition de la chaîne de Markov associée est stockée sous un format tensoriel. Particulièrement, le formalisme des réseaux d'automates stochastiques est employé pour la description des modèles. L'évaluation de performances cherchée est la détermination de l'état stationnaire de la chaîne de Markov (\emph(résolution)). De ce fait, les propriétés de l'algèbre tensorielle généralisée sont proposées et démontrées de façon à établir la base nécessaire aux algorithmes de résolution introduits. Le principal apport de cette thèse réside dans l'efficacité des ces algorithmes, qui est obtenue avec l'accélération des méthodes itératives. Ceci est fait à deux niveaux: la réduction du coût de chaque itération; et la réduction du nombre d'itérations nécessaire à la convergence. La multiplication d'un vecteur par une matrice sous format tensoriel (produit vecteur-descripteur) est l'opération de base des itérations. L'efficacité de cette opération est le premier objectif à atteindre. Le deuxième objectif est l'implémentation des méthodes de la puissance, d'Arnoldi et GMRES dans ses versions standards et pré-conditionnées de façon a minimiser le nombre d'itérations sans trop augmenter le coût de chaque itération. La totalité des concepts introduits est alors utilisée dans le logiciel PEPS 2.0. Plusieurs exemples pratiques de modèles en réseaux d'automates stochastiques ont été mesurés sur PEPS 2.0 pour illustrer les résultats de cette thèse.

Ce travail est divisé en deux axes de recherches. Le premier axe concerne l'étude de quelques systèmes d'équations aux dérivées partielles non linéaires issus de la modélisation macroscopique des composants semi-conducteurs. Le deuxième axe de recherche est consacré à l'étude de quelques problèmes d'identification. Nous nous intéressons en particulier à deux types de problèmes d'identification. Le premier concerne la reconstruction des données sur le bord pour des problèmes elliptiques. Le deuxième type de problèmes auquel nous nous sommes intéressés est celui de l'identification des frontières dans des problèmes gouvernés par des équations elliptiques.

Dans cette thèse nous nous intéressons aux structures des algèbres quotients et plus particulièrement à l'apport de la dualité pour la représentation des algèbres de coordonnées. Une première partie de cette thèse est consacrée à la représentation des algèbres de dimension zéro et à des applications de la dualité à des problèmes d'interpolation. Nous généralisons les bases d'interpolation de Lagrange et d'Hermite pour lesquelles nous donnons des formules explicites. Cela nous permet de donner les relations entre les racines d'un système algébrique et ses coefficients avec des formules généralisant celles du cas univarié. Dans une deuxième partie, nous appliquons les résultats développés dans la première partie à la conception de méthodes itératives pour l'approximation simultanée de l'ensemble des solutions d'un système algébrique. La troisième partie est consacrée aux résidus algébriques. Nous rappelons les notions relatives aux algèbres de Gorenstein et à leurs représentations. Nous introduisons les bézoutiens et les résidus algébriques dont nous donnons des applications en géométrie. Dans la quatrième partie, nous nous intéressons à l'algorithmique associé aux matrices quasi-Toeplitz, quasi-Hankel, ..., telles que définies par B. Mourrain et V.Y. Pan. Nous en montrons des applications dans le cadre de l'algorithmique permettant des accélérations asymptotiques de méthodes de résolution de systèmes algébriques.

Le but de la thèse est avant tout d'ordre méthodologique : il s'agit de présenter et de valider une méthode de modélisation numérique pour l'équation des ondes élastiques (2D et 3D) dans des milieux réels. De tels milieux peuvent notamment présenter une forte topographie, ainsi que des couches de forme quelconque avec de forts contrastes de propriétés élastiques et de densité. Nous cherchons à modéliser l'ensemble du champ d'onde, c'est-à-dire tant les ondes de volume (compression et cisaillement) que les ondes de surface et d'interface. <br /><br />Dans l'introduction générale, nous définissons les différents types de modèles géophysiques que nous souhaitons étudier, et rappelons l'intérêt qu'il peut y avoir à être capable de calculer le comportement des ondes sismiques dans de telles structures. Nous passons brièvement en revue les différentes techniques de modélisation numérique ordinairement utilisées pour résoudre les problèmes de propagation d'ondes dans de tels milieux, et rappelons sommairement leurs principales caractéristiques ainsi que leurs principales limitations. Les avantages potentiels de l'approche que nous introduisons dans cette thèse sont déduits de ce rapide tour d'horizon. <br /><br />Dans une première partie, nous rappelons brièvement les lois essentielles de la mécanique des milieux continus et de l'élastodynamique. Nous mentionnons les principaux types d'ondes pouvant exister dans un milieu élastique. Les propriétés particulières des ondes de surface et d'interface dans de tels modèles sont également mentionnées. <br /><br />Dans une deuxième partie, nous introduisons et validons une méthode pseudo-spectrale globale fondée sur la formulation différentielle tensorielle des équations de l'élastodynamique, adaptée au traitement de géométries déformées et de surfaces non planes. Nous montrons son application à quelques cas simples, nous mettons en évidence des effets intéressants liés à la topographie, mais soulignons aussi ses limitations intrinsèques. Nous en déduisons la nécessité d'utiliser une méthode plus souple que la formulation différentielle classique (différences finies, méthodes spectrales ou pseudo-spectrales) pour le traitement de cas réalistes (typiquement, une méthode variationnelle avec décomposition de domaine).<br /><br />Dans une troisième partie, constituant le cœur de ce travail de thèse, nous introduisons et développons tant à 2D qu'à 3D une formulation variationnelle d'ordre élevé des équations de l'élastodynamique, dite " méthode des éléments spectraux ", et nous mettons en évidence ses propriétés sur différents problèmes classiques de complexité croissante (problème de Lamb, problème de Garvin, onde de Rayleigh sur une surface courbe...). Puis nous présentons l'application de cette méthode à des modèles 2D et 3D plus complexes, pour lesquels de forts effets liés notamment à la topographie sont mis en évidence. Nous obtenons en particulier de fortes amplifications locales du champ d'accélération et de déplacement dans de tels modèles (effets de site). Dans le cas d'un modèle tridimensionnel, nous montrons que ces effets dépendent fortement de la polarisation du champ incident. L'efficacité de la méthode sur un calculateur parallèle est également discutée et illustrée. <br /><br />Dans une quatrième partie, nous avons regroupé deux petites études que nous avons menées au cours de nos recherches pour résoudre l'équation des ondes à 1D (ondelettes et volumes finis).

Cette thèse s'articule autour de deux thématiques principales. La première traite de l'analyse de performance dans les contextes des systèmes utilisant les séquences d'apprentissage et des systèmes multi-utilisateurs. En nous basant sur une analyse asymptotique, nous donnons des expressions explicites des taux d'erreur binaire et de la probabilité de coupure et nous montrons leur précision même pour des systèmes de dimensions réduites. La deuxième partie de cette thèse traite des méthodes d'estimation aveugle et semi-aveugle. Nos contributions se situent sur deux plans: algorithmique et théorique. Sur le plan algorithmique, nous proposons de nouvelles approches qui permettent de résoudre certains problèmes des méthodes conventionnelles telque la sensibilité à la surestimation de l'ordre pour les méthodes aveugles et la recherche du paramètre de régularisation pour les méthodes semi-aveugles. Sur le plan théorique, nous apportons la démonstration de la quasi-convexité de l'erreur asymptotique quadratique moyenne des systèmes semi-aveugle basé sur la régularisation. Nous menons aussi une étude théorique qui témoigne de l'efficacité de l'utilisation de la norme lp pour robustifier les méthodes aveugles.

Une des difficultés fondamentales de la simulation numérique de plasmas magnétisés est l'existence d'échelles de temps et d'espace mettant en jeu plusieurs ordres de grandeurs très différents. Afin de réaliser des simulations numériques efficaces de ces phénomènes physiques, il est essentiel de développer des modèles et des méthodes numériques adaptés à ces problèmes. <br />A ce jour, la notion de convergence 2-échelles introduite par G. Allaire et G. Nguetseng est un des outils permettant de dériver rigoureusement des limites multi-échelles, ce qui nous permet d'obtenir des modèles limites qu'il est possible de discrétiser avec une méthode numérique usuelle : nous parlons alors d'une méthode numérique 2-échelles. <br />L'objectif de cette thèse est de développer une méthode semi-lagrangienne 2 échelles sur un modèle de type Vlasov gyrocinétique afin de simuler un plasma soumis à un champ magnétique fort du même type que ceux utilisés pour le projet ITER. Cependant, comme les phénomènes physiques à simuler sont assez complexes et comme nous ne savons que peu de choses sur le comportement d'une méthode numérique 2-échelles sur un modèle non-linéaire, il convient de procéder par étapes avant de développer une telle méthode sur un modèle gyrocinétique. <br />Dans une première partie, nous construisons une méthode de volumes finis 2-échelles sur les équations d'Euler 1D isentropiques faiblement compressibles. Bien que ce modèle soit assez différent d'un modèle de type Vlasov, il n'en est pas moins un cadre de travail relativement simple pour étudier le comportement d'une méthode numérique 2-échelles face à un modèle non-linéaire. <br />Dans une seconde partie, nous nous basons sur le modèle limite développé par E. Frénod, F. Salvarani et E. Sonnendrücker afin de construire une méthode semi-lagrangienne 2-échelles pour simuler des faisceaux de particules en géométrie axisymétrique. Même si le modèle de Vlasov axisymétrique utilisé est différent d'un modèle gyrocinétique, il constitue un contexte idéal pour établir les bases d'une méthode semi-lagrangienne 2 échelles.<br />Enfin, dans une troisième partie, nous utilisons la convergence 2-échelles afin d'améliorer les résultats de convergence faible-* établis par M. Bostan en 2007, et nous proposons une méthode semi-lagrangienne en avant permettant de valider numériquement ces résultats.

Dans cette thèse nous proposons de développer un système d'imagerie ultrasonore innovante de micro- défauts basé sur l'utilisation conjointe de techniques d'acoustique non linéaire et du concept de "transducteur à cavité chaotique". Ce transducteur correspond à la combinaison d'une céramique piézoélectrique collée sur une cavité de forme chaotique et du principe de retournement temporel. La faisabilité et les performances de ce nouveau système sont explorées par des simulations numériques. Des paramètres optimaux d'utilisation pour une implémentation expérimentale sont proposés. Une grande partie des travaux menés dans le cadre de cette thèse se concentre sur le développement d'outils numériques permettant l'amélioration de telles techniques d'imagerie. Un schéma d'éléments finis de type Galerkin Discontinu (GD) est étendu à l'élastodynamique non linéaire. Un type de zone absorbante parfaitement adaptée, appelée "Nearly Perfectly Matched Layer" (NPML) a aussi été développé. Dans le cas de matériaux orthotropes, comme des problèmes de stabilité apparaissent, un mélange de NPML et de zone atténuante, dont on contrôle la proportion respective, est introduit afin de stabiliser les NPML. Une validation expérimentale du concept de "transducteur à cavité chaotique" pour la focalisation dans un milieu solide, réverbérant ou non, en utilisant une seule source est réalisée. Les méthodes de retournement temporel et de filtre inverse sont présentées et comparées. La démonstration expérimentale qu'un "transducteur à cavité chaotique" peut être utilisé conjointement avec les méthodes d'inversion d'impulsion afin de réaliser une image de non linéarités localisées est présentée

Une antenne à réflecteur est un dispositif largement utilisé pour la communication satellite. La durée de vie d'un tel dispositif est étroitement liée à la fatigue due à la consommation d'énergie pour émettre le signal. Un des enjeux de la conception optimale d'une antenne revient donc à produire des systèmes dont le rendement est le meilleur possible par rapport à une tâche donnée. La particularité d'une antenne à réflecteur se traduit par la présence de surfaces rayonnantes dont la géométrie constitue le paramètre principal pour assumer cette tâche. Sur la base de la simulation de la propagation d'une onde électromagnétique en espace libre et en régime harmonique, on est capable de développer des méthodes d'optimisation numérique de la forme de surfaces rayonnantes. On cherche à minimiser un critère qui traduit en terme mathématique la tâche à effectuer d'un point de vue énergétique. Cependant, les méthodes utilisées sont souvent soumis à des difficultés liées au fait que ces problèmes sont mal posés et numériquement raides. Le contrôle étant géométrique, on a examiné dans cette thèse les contributions potentielles de représentations hiérarchiques afin d'étendre les performances d'algorithmes classiques d'optimisation. Ces extensions empruntent ses fondements aux méthodes multigrilles pour la résolution d'EDP. Un exemple théorique d'optimisation de forme permet d'assoir les stratégies appliquées à l'optimisation d'antennes. Puis des expériences numériques d'optimisation montrent que les algorithmes de bases sont améliorés en terme de robustesse comme en terme de vitesse de convergence.

Dans ce travail, nous nous intéressons à l'analyse numérique de problèmes aux valeurs propres non linéaires, comme on peut en trouver en chimie quantique ou en mécanique. La résolution de ces problèmes étant très coûteuse, l'idée est de proposer de nouvelles méthodes permettant de simplifier la résolution de ce type de problèmes et ainsi diminuer le coût de calcul. L'analyse numérique est nécessaire pour comprendre si l'impact positif sur le coût de calcul total n'a pas de mauvaise conséquence sur la précision des résultats. On propose un complément aux travaux existants sur les estimations d'erreur a priori, afin d'obtenir des résultats équivalents à ceux connus dans le cas de problèmes aux valeurs propres linéaires. Ces résultats ont été utilisés pour la mise en oeuvre et l'analyse numérique de nouveaux schémas à deux grilles pour l'approximation de problèmes aux valeurs propres non linéaires. Ensuite, on propose d'adapter ce type de méthode de sous-grilles, pour une utilisation associée à la méthode des bases réduites.