Cette thèse porte sur l'étude d'une méthode de type Galerkin discontinu en domaine temporel (GDDT), afin de résoudre numériquement les équations de Maxwell instationnaires sur des maillages hybrides tétraédriques/hexaédriques en 3D (triangulaires/quadrangulaires en 2D) et non-conformes, que l'on note méthode GDDT-PpQk. Comme dans différents travaux déjà réalisés sur plusieurs méthodes hybrides (par exemple des combinaisons entre des méthodes Volumes Finis et Différences Finies, Éléments Finis et Différences Finies, etc.), notre objectif principal est de mailler des objets ayant une géométrie complexe à l'aide de tétraèdres, pour obtenir une précision optimale, et de mailler le reste du domaine (le vide environnant) à l'aide d'hexaèdres impliquant un gain en terme de mémoire et de temps de calcul. Dans la méthode GDDT considérée, nous utilisons des schémas de discrétisation spatiale basés sur une interpolation polynomiale nodale, d'ordre arbitraire, pour approximer le champ électromagnétique. Nous utilisons un flux centré pour approcher les intégrales de surface et un schéma d'intégration en temps de type saute-mouton d'ordre deux ou d'ordre quatre. Après avoir introduit le contexte historique et physique des équations de Maxwell, nous présentons les étapes détaillées de la méthode GDDT-PpQk. Nous réalisons ensuite une analyse de stabilité L2 théorique, en montrant que cette méthode conserve une énergie discrète et en exhibant une condition suffisante de stabilité de type CFL sur le pas de temps, ainsi que l'analyse de convergence en h (théorique également), conduisant à un estimateur d'erreur a-priori. Ensuite, nous menons une étude numérique complète en 2D (ondes TMz), pour différents cas tests, des maillages hybrides et non-conformes, et pour des milieux de propagation homogènes ou hétérogènes. Nous faisons enfin de même pour la mise en oeuvre en 3D, avec des simulations réalistes, comme par exemple la propagation d'une onde électromagnétique dans un modèle hétérogène de tête humaine. Nous montrons alors la cohérence entre les résultats mathématiques et numériques de cette méthode GDDT-PpQk, ainsi que ses apports en termes de précision et de temps de calcul.

Mes travaux de recherche portent sur l'analyse numérique des intégrateurs géométriques et multi-échelles pour les équations différentielles déterministes ou stochastiques. Les modèles d'équations différentielles issus de la physique ou la chimie possèdent souvent une structure géométrique ou multi-échelles particulière (par exemple, les structures hamiltoniennes, les intégrales premières, les structures multi-échelles en temps ou en espace, les systèmes hautement oscillatoires), mais leur complexité est souvent telle qu'une solution satisfaisante est hors de portée en utilisant seulement des méthodes numériques standards à usage général. L'objectif est donc d'identifier les propriétés géométriques ou multi-échelles pertinentes de ces problèmes, et d'en tirer avantage pour concevoir et analyser de nouveaux intégrateurs efficaces, fiables et précis, reproduisant fidèlement le comportement qualitatif de la solution exacte des modèles considérés.

Le but général de cette thèse est d'exploiter des modélisations asymptotiques pour la résolution de problèmes de diffraction inverse en électromagnétisme. Nous nous intéressons plus particulièrement au cas des conditions d'impédance généralisée qui modélisent notamment des matériaux fortement absorbants ou des revêtements de faible épaisseur. L'expression "impédance généralisée" signifie que la condition au bord fait intervenir un opérateur surfacique. Les conditions dites d'impédance classique entrent dans cette famille de conditions aux bord, dans ce cas, l'opérateur surfacique se réduit à la multiplication par une fonction. Dans le cadre des problèmes inverses, l'utilisation de modèles approchés permet de simplifier aussi bien la résolution numérique que l'analyse mathématique. De nombreux travaux ont été menés en diffraction inverse sur l'utilisation d'une condition d'impédance classique, nous les avons étendus pour des opérateurs surfaciques plus complexes faisant intervenir des dérivées tangentielles. Une partie importante de la thèse est consacrée à la mise en oeuvre des méthodes d'optimisation pour retrouver un obstacle ainsi que les paramètres définissant l'opérateur d'impédance. Nous présentons en particulier un calcul de dérivée de forme dans le cas où les équations de l'électromagnétisme se simplifient en une équation scalaire, nous étendons ensuite ce calcul aux équations de Maxwell vectorielles. Des exemples numériques de reconstruction de forme et de paramètres d'impédance viennent illustrer l'applicabilité des méthodes d'optimisation à notre problème inverse. Afin de compléter cette étude, nous avons utilisé une méthode qualitative - la méthode de factorisation - pour identifier un objet diffractant caractérisé par une condition d'impédance généralisée. Enfin, en relation avec les méthodes qualitatives, nous nous sommes penché sur l'utilisation des valeurs propres de transmission associées au problème de diffraction par des couches minces pour obtenir des informations sur la couche. Dans ce but, nous avons calculé et justifié le développement asymptotique de la première valeur propre de transmission intérieure par rapport à la faible épaisseur du revêtement. Ce développement donne une manière simple de calculer l'épaisseur du revêtement à partir du champ diffracté pour plusieurs fréquences.

Les écoulements de mousse liquide sont directement impliqués dans de nombreuses applications dans des domaines aussi variés que les industries agro-alimentaire et cosmétique, l'extraction pétrolière, ou encore la décontamination nucléaire. Par ailleurs, l'étude des mousses apporte des connaissances fondamentales : plus facile à manipuler et analyser, la mousse est un fluide modèle pour comprendre des matériaux tels que les émulsions, les polymères, les pâtes, ou les agrégats de cellules, qui possèdent à la fois des propriétés liquides et solides. Les expériences systématiques réalisées par l'équipe de F. Graner ont fourni une série de données précises qui mettent l'accent sur les propriétés non newtoniennes de la mousse. Dans le même temps, P. Saramito a proposé un modèle visco-élasto-plastique (VEP) continu et tensoriel, à même de prédire le comportement de la mousse. L'objectif de cette thèse est de comprendre ce comportement complexe en s'appuyant sur ces deux éléments. Nous avons élaboré et validé un algorithme de résolution basé sur une méthode d'éléments finis bidimensionnelle. Les solutions numériques sont en excellent accord avec tous les champs mesurés (vitesse, déformation élastique, taux de déformation plastique), et nous avons confirmé le caractère prédictif du modèle. Nous avons identifié les paramètres dominants et établi la nécessité de traiter simultanément les contributions visqueuse, élastique, et plastique dans un formalisme tensoriel. Notre travail apporte une contribution substantielle à la compréhension des mousses et ouvre la voie à la simulation réaliste d'écoulements complexes de fluides VEP en vue d'applications industrielles.

Cette thèse s'intéresse aux méthodes de champs de phase pour l'approximation de mouvements par courbure moyenne. En particulier, nous proposons de nouveaux modèles pour prendre en compte des contraintes de volumes, des tensions de surfaces anisotropes et des angles de contact au bord du domaine d'évolution.

Cette thèse est dédiée à l'étude mathématique et numérique d'une classe d'équations cinétiques collisionnelles, de type équation de Boltzmann. Nous avons porté un intérêt tout particulier à l'équation des milieux (ou gaz) granulaires, initialement introduite dans la littérature physique pour décrire le comportement hors équilibre de matériaux composés d'un grand nombre de grains, ou particules, non nécessairement microscopiques, et interagissant par des collisions dissipant l'énergie cinétique. Ces modèles se sont révélés avoir une structure mathématique très riche. Cette thèse se structure en trois partie pouvant être lues de manière indépendante, mais néanmoins en rapport avec des équations cinétiques collisionnelles en général, et l'équation des milieux granulaires en particulier. La première partie est dédiée à l'étude mathématique du comportement asymptotique de certaines équations cinétiques collisionnelles dans un cadre homogène en espace. Nous y montrons des résultats de type explosion et convergence vers la solution autosimilaire avec calcul explicite des taux, pour des opérateurs de type Boltzmann, grâce à l'utilisation (entre autre) d'une nouvelle méthode de changement de variables dépendant directement de la solution de l'équation considérée. En particulier, nous démontrons que pour un modèle de gaz granulaire - dit anormal - il est possible d'observer une explosion en temps fini. Dans la deuxième partie, orientée analyse numérique et calcul scientifique, nous nous intéressons développement et à l'étude de méthodes spectrales pour la résolution de problèmes multi-échelles, issus de la théorie des équations cinétiques collisionnelles. Les méthodes de changement de variables tiennent aussi une place importante dans cette partie, et permettent d'observer numériquement des phénomènes non triviaux qui apparaissent lors de l'étude de gaz granulaires, comme la création d'amas de matière ou la caractérisation précise du retour vers l'équilibre. La troisième et dernière partie est dédiée à l'étude spectrale de l'opérateur des milieux granulaires avec bain thermique, linéarisé au voisinage d'un équilibre homogène en espace, afin d'établir des résultats de type stabilité et convergence vers une limite hydrodynamique. Ce travail est en fait la généralisation d'un résultat célèbre dans la théorie de l'équation de Boltzmann, dû à R. Ellis et M. Pinsky, et établissant rigoureusement la première limite hydrodynamique vers les équations d'Euler compressibles linéaires puis Navier-Stokes de cette équation.

La thèse est portée essentiellement sur la stabilisation et la contrôlabilité de deux équations des ondes moyennant un seul contrôle agissant sur le bord du domaine. Dans le cas du contrôle dynamique, le contrôle est introduit dans le système par une équation différentielle agissant sur le bord. C'est en effet un système hybride. Le contrôle peut être aussi applique directement sur le bord d'une équation, c'est le cas du contrôle indirecte mais non borne. La nature du système ainsi coupledépend du couplage des équations, et ceci donne divers résultats par la stabilisation (exponentielle et polynomiale) et la contrôlabilité exacte (espace contrôlable). Des nouvelles inégalités d'énergie permettent de mettre en oeuvre la Méthode fréquentielle et la Méthode d'Unicité de Hilbert.

Dans une première partie, on rassemble plusieurs résultats en théorie du contrôle autour des inégalités d'Ingham, généralisations de l'égalité de Parseval, qui inter- viennent pour montrer l'observabilité, la contrôlabilité ou la stabilisation frontière ou interne de l'équation des ondes ou d'équations similaires dans certains cas parti- culiers. On s'intéresse dans un premier temps à l'optimalité de ce type d'inégalités en généralisant un résultat précédent au cas vectoriel. On développe ensuite un théo- rème de type Ingham adapté pour traiter le cas d'une géométrie cartésienne. Enfin, on donne des résultats d'observabilité dans le cas d'approximations numériques. Dans une seconde partie, on présente les méthodes semi-Lagrangiennes qui sont composées essentiellement de deux ingrédients : calcul des caractéristiques le long desquelles la fonction de distribution est constante et étape d'interpolation. On ana- lyse des schémas d'ordre élevé en temps pour le système de Vlasov-Poisson 1D×1D, basés sur le splitting directionnel, qui est une succession d'étapes de transport li- néaire. On étudie alors les méthodes semi-Lagrangiennes dans ce cas particulier et on fait le lien entre différentes formulations. On obtient également un théorème de convergence pour le système de Vlasov-Poisson dans ce cadre, qui reste valable pour des petits déplacements. On développe ensuite ce type de méthodes dans un cadre plus général, en se basant sur le splitting uni-dimensionnel conservatif, avec une variante de type Galerkin discontinu. Dans une dernière partie, on étudie l'opérateur de gyromoyenne qui intervient en physique des plasmas pour prendre en compte des corrections de rayon de Larmor fini. Enfin, on discute de la problématique de la divergence discrète nulle qui donne une compatibilité entre le calcul du champ et la méthode numérique de transport.

Ce travail présente quelques résultats sur les systèmes de particules en interaction pour l'interprétation probabiliste des équations aux dérivées partielles, avec des applications à des questions de dynamique moléculaire et de chimie quantique. On présente notamment une méthode particulaire permettant d'analyser le processus de la force biaisante adaptative, utilisé en dynamique moléculaire pour le calcul de différences d'énergies libres. On étudie également la sensibilité de dynamiques stochastiques par rapport à un paramètre, en vue du calcul des forces dans l'approximation de Born-Oppenheimer pour rechercher l'état quantique fondamental de molécules. Enfin, on présente un schéma numérique basé sur un système de particules pour résoudre des lois de conservation scalaires, avec un terme de diffusion anormale se traduisant par une dynamique de sauts sur les particules

Les modèles mathématiques sont importants pour la compréhension de la dynamique de l'atmosphère et de l'océan. Mais si ils étaient notre seule source d'information aucune prévision ne serait possible faute, notamment au manque de la connaissance d'une condition initiale cohérente. On dispose également d'observations de ces systèmes en nombre de plus en plus important, notamment grâce aux nombreux satellites d'observation qui croisent maintenant au large de notre planète. Ces observations sont souvent indirectes et incomplètes, et de ce fait ne fournissent pas non plus, à elles seules, une connaissance approfondie de l'état du milieu considéré. Et pour finir, on dispose de statistiques sur les champs des variables atmosphériques, leur variabilité, leur cohérence en temps et en espace. Je présente donc, dans ce document, des méthodes permettant de combiner tout ou partie de ces informations afin d'améliorer la prévision et la connaissance du fonctionnement de ces systèmes. Ces méthodes se basent le plus souvent sur une théorie mathématique solide, mais les appliquer dans un contexte réaliste n'est pas toujours chose aisée. C'est pourquoi on gardera le souci d'accompagner les développements que nous effectuons jusqu'à des applications opérationnelles ou quasi opérationnelles afin de démontrer la faisabilité de ceux ci.