Ce manuscrit se décompose en deux parties principales, associées à deux types d'applications assez différentes. Dans la première partie qui comprend les deux premiers chapitres, je m'intéresse à des systèmes issus de problèmes de contrôle et d'estimation en physique quantique; dans la deuxième partie (troisième chapitre du manuscrit), j'étudie la propagation d'ondes électriques le long des fils classiques dans un réseau de lignes de transmission et je considère certains problèmes d'estimation de paramètres. Dans le premier chapitre nous étudions le problème de la planification de trajectoires pour des systèmes quantiques fermés modélisés par des équations de Schrödinger bilinéaire. Nous démontrons alors des résultats de la stabilisation approchée pour le cas d'une boite quantique infinie ainsi que pour le cas d'un potentiel décroissant. Dans les deux cas, le manque de pré-compacité des trajectoires dans des espaces fonctionnels appropriés nous oblige à proposer des méthodes de Lyapunov qui évitent des phénomènes de perte de masse à l'infini. Dans le deuxième chapitre nous étudions le problème de stabilisation de systèmes quantiques en observation. Cette observation nécessite l'ouverture du système à son environnement. Les modèles pertinents pour l'évolution de ce type de systèmes sont des modèles stochastiques basés sur des trajectoires de Monte-Carlo quantiques. Nous étudions alors certains problèmes de stabilisation qui parviennent de vraies expériences physiques. Enfin, dans le chapitre 3 nous considérons le problème d'estimation de paramètres pour un réseau de fils de câblage électrique. Dans ce but, nous étudions deux approches : l'approche temporelle et l'approche fréquentielle. Dans l'approche temporelle, nous considérons le réseau le plus simple qui consiste d'une seule ligne de transmission et nous proposons un algorithme d'identification pour l'équation d'onde associé qui est basé sur l'application des observateurs asymptotiques. Dans l'approche fréquentielle, nous considérons un réseau plus compliqué de la forme étoile. Nous proposons alors des résultats d'identifiabilité basés sur des techniques de l'inverse scattering.

Mes travaux de thèse portent sur la construction de conditions aux limites absorbantes (CLAs) pour des problèmes de propagation d'ondes posés dans des milieux limités par des surfaces régulières. Ces conditions sont nouvelles car elles prennent en compte non seulement les ondes proagatives (comme la plupart des CLAs existantes) mais aussi les ondes évanescentes et rampantes. Elles sont donc plus performantes que les conditions existantes. De plus, elles sont facilement implémentables dans un schéma d'éléments finis de type Galerkine Discontinu (DG) et ne modifie pas la condition de stabilité de Courant-Friedrichs-Lewy (CFL). Ces CLAs ont été implémentées dans un code simulant la propagation des ondes acoustiques ainsi que dans un code simulant la propagation des ondes en régime harmonique. Les comparaisons réalisées entre les nouvelles conditions et celles qui sont les plus utilisées dans la littérature montrent que prendre en compte les ondes évanescentes et les ondes rampantes permet de diminuer les réflexions issues de la frontière artificielle et donc de rapprocher la frontière artificielle du bord de l'obstacle. On limite ainsi les coûts de calcul, ce qui est un des avantages de mes travaux. De plus, compte tenu du fait que les nouvelles CLAs sont écrites pour des frontières quelconques, elles permettent de mieux adapter le domaine de calcul à la forme de l'obstacle et permettent ainsi de diminuer encore plus les coûts de calcul numérique.

Nous étudions dans ce travail une méthode de type éléments finis dans le but de simuler le mouvement de particules rigides immergées. La méthode développée ici est une méthode de type domaine fictif. L'idée est de chercher un prolongement régulier de la solution exacte à tout le domaine fictif afin d'obtenir une solution régulière sur tout le domaine et retrouver l'ordre optimal de l'erreur avec des éléments d'ordre 1. Le prolongement régulier est cherché en minimisant une fonctionnelle dont le gradient est donné par la solution d'un nouveau problème fluide faisant intervenir une distribution simple couche dans le second membre. Nous faisons une analyse numérique, dans le cas scalaire, de l'approximation de cette distribution par une combinaison de masse de Dirac. Un des avantages de cette méthode est de pouvoir utiliser des solveurs rapides sur maillages cartésiens tout en conservant l'ordre optimal de l'erreur. Un autre avantage de la méthode vient du fait que les opérateurs ne sont pas modifiés, seul les seconds membres dépendent de la géométrie du domaine initial. Nous avons de plus écrit un code C++ parallèle en deux et trois dimensions, permettant de simuler des écoulements fluide/particules rigides avec cette méthode. Nous présentons ainsi une description des principales composantes de ce code.

Les méthodes classiques de résolution des équations de l'élasticité linéaire sont les méthodes éléments finis. Ces méthodes produisent de très bons résultats et sont très largement analysées mathématiquement pour l'étude des déformations solides. Pour des problèmes de couplage solide/fluide, pour des situations réalistes en présence de discontinuités (modélisation des fronts de gel dans les sols humides), ou bien encore pour des domaines de calcul mieux adaptés aux maillages non conformes, il parait intéressant de disposer de solveurs Volumes Finis. Les méthodes Volumes Finis sont très largement utilisées en mécanique des fluides. Appliquées aux problèmes de convection, elles sont bien adaptées à la capture de solutions présentant des discontinuités et ne nécessitent pas de maillages conformes. De plus, elles présentent l'avantage de conserver au niveau discret les flux à travers les interfaces du maillage. C'est pourquoi sont développées et testées, dans cette thèse, plusieurs méthodes de volumes finis, qui permettent de traiter le problème de l'élasticité. On a, dans un premier temps, mis en œuvre la méthode LSGR (Least Squares Gradient Reconstruction), qui reconstruit des gradients par volumes à partir d'une formule de moindres carrés pondérés sur les volumes voisins. Elle est testée pour des maillages tétraédriques non structurés, et montre un ordre 1 de convergence. La méthode des Volumes Finis mixtes est ensuite présentée, basée sur la conservation d'un flux "pénalisé" à travers les interfaces. Cette pénalisation impose une contrainte sur le type de maillage utilisé, et des tests sont réalisés en 2d avec des maillages structurés et non structurés de quadrangles. On étend ensuite la méthode des Volumes Finis diamants à l'élasticité. Cette méthode détermine un gradient discret sur des sous volumes associés aux interfaces à partir de l'interpolation de la solution aux sommets du maillage. La convergence théorique est prouvée sous réserve de vérifier une condition de coercivité. Les résultats numériques, en 2d pour des maillages non structurés, conduisent à un ordre de convergence meilleur que celui prouvé. Enfin, la méthode DDFV (Discrete Duality Finite Volume), qui est une extension de la méthode Diamant, est présentée. Elle est basée sur une correspondance entre plusieurs maillages afin d'y construire des opérateurs discrets en "dualité discrète". On montre que la méthode est convergente d'ordre 1. Les illustrations numériques, réalisées en 2d et en 3d pour des maillages non structurés, montrent une convergence d'ordre 2, ce qui est fréquemment observé pour cette méthode.

Dans cette thèse on étudie le couplage de modèles à dimensions spatiales hétérogènes, avec une application en hydraulique fluviale. On commence par donner un cadre mathématique général à cette problématique. Ensuite, on étudie un cas académique de couplage 1-D/2-D dans le cadre elliptique. On définit les opérateurs de restriction et d'extension nécessaires à l'analyse en se basant sur la dérivation du modèle 1-D à partir du modèle 2-D. Après cela, on met en oeuvre un algorithme de couplage de type Schwarz avec des conditions de type Robin. On montre alors la convergence de cet algorithme, plus particulièrement sa convergence optimale en utilisant l'opérateur absorbant exact 1-D. On termine cette partie en établissant une majoration de l'erreur entre la solution couplée et la solution globale de référence en fonction du rapport d'aspect du domaine d'étude et de la position de l'interface de couplage. Ces résultats sont illustrés numériquement. Dans la deuxième partie, on généralise cette analyse mathématique au cas du couplage des systèmes linéaires de Saint-Venant 2-D et de Navier-Stokes hydrostatiques 3-D. En faisant l'hypothèse d'une friction nulle au fond, on montre que la convergence de l'algorithme de couplage est équivalente à celle de l'algorithme usuel de décomposition de domaine du Système de Saint-Venant. On calcule une approximation des opérateurs absorbants exacts du système de Saint-Venant, ce qui nécessite de faire des hypothèses restrictives. Comme alternative, on propose un algorithme avec des conditions de type Robin, dont on montre la convergence. Enfin, on présente une première étude d'un cas test réel de couplage des systèmes de Saint-Venant 1-D et Navier-Stokes 3-D en utilisant les codes numériques Mascaret 1-D et Telemac 3-D développés par EDF R&D.

Cette thèse est consacrée à l'étude mathématique des propriétés de stabilité de systèmes quantiques infinis, décrits par des modèles non linéaires. Dans les chapitres 1 et 2, on étudie l'instabilité du vide relativiste menant au phénomène de création de paires électron-positron. Dans le chapitre 3, on considère la dynamique de ce même vide relativiste couplé à un champ scalaire. Les chapitres 4 et 5 sont consacrés au caractère dispersif de la dynamique non linéaire de Hartree pour des perturbations de la mer de Fermi, et en particulier à sa stabilité orbitale et asymptotique. Enfin, le chapitre 6 introduit une notion générale d'entropie relative entre deux états comportant une infinité de particules.

Nous présentons des modèles macroscopiques du comportement mécanique du parenchyme pulmonaire humain obtenus par homogénéisation double--échelle. Dans une première partie consacrée au couplage entre parenchyme et arbre bronchique, nous commençons par proposer un modèle de la déformation du parenchyme. Nous modélisons (i) le parenchyme par un matériau élastique linéaire, (ii) les alvéoles comme des cavités de taille epsilon réparties périodiquement dans le domaine macroscopique et (iii) l'arbre bronchique par un arbre dyadique résistif en supposant la loi de Poiseuille valide pour chaque voie aérienne. Nous obtenons en faisant tendre epsilon vers zéro, par homogénéisation, une description macroscopique du parenchyme comme un matériau viscoélastique, sous certaines conditions sur l'irrigation du domaine par l'arbre que nous étudions ensuite. L'arbre induit une dissipation non--locale en espace. Nous illustrons ces résultats par des résultats numériques. Dans une deuxième partie, nous étudions la propagation d'ondes sonores dans le parenchyme, sans prendre en compte l'arbre bronchique. Nous homogénéisons dans le domaine fréquentiel un premier modèle couplant l'élasticité linéarisée dans le parenchyme avec l'équation acoustique dans l'air. Nous déduisons une loi de comportement macroscopique élastique linéaire en dehors d'un ensemble de résonances. Ensuite, nous homogénéisons un deuxième modèle qui prend en compte le caractère viscoélastique et inhomogène du parenchyme au niveau microscopique. Le matériau macroscopique présente des effets de mémoire nouveaux par rapport aux composants microscopiques que nous étudions numériquement.

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.

En simulation de réservoir, les variations spatiales des propriétés des roches sont données par un modèle géologique défini sur une grille pouvant comporter plusieurs centaines de millions de mailles. Simuler des écoulements polyphasiques sur ce type de grille peut être très coûteux en temps de calculs. Les ingénieurs de gisement procèdent souvent à une mise à l'échelle de ces propriétés sur une grille plus grossière, utilisée ensuite pour la simulation. Cette façon de procéder ne permet pas toujours de bien prendre en compte l'influence sur l'écoulement des hétérogénéités présentes à l'échelle fine. Plusieurs méthodes multi-échelles ont donc été proposées ces dernières années afin d'obtenir un meilleur compromis entre les temps de calculs et la précision des solutions obtenues. Actuellement, la plupart de ces méthodes permettent de résoudre plus efficacement l'équation en pression. Les équations de transport sont, quant à elles, toujours résolues sur la grille fine. La première partie de cette thèse présente la mise en œuvre et les résultats de tests réalisés sur la méthode des éléments finis mixtes multi-échelles dans un environnement de calcul parallèle. Le second et principal objectif de ce travail est de proposer une nouvelle méthode multi-échelle pour la simulation du transport d'un traceur dans un milieu poreux. Cette méthode est construite à partir de résultats rigoureux d'homogénéisation en milieux périodiques. Cette thèse fournit une vitesse de convergence explicite pour ce procédé d'homogénéisation ainsi qu'une estimation d'erreur pour la méthode numérique multi-échelle associée. Des exemples numériques sont ensuite présentés.

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.