Cette thèse analyse le spectre d'opérateurs de Schrödinger avec champ magnétique constant dans des ouverts de type diédraux. Pour comprendre l'influence d'une arête courbe sur la première valeur propre de l'opérateur dans la limite semi-classique, il faut connaître le bas du spectre de l'opérateur de Schrödinger magnétique avec champ constant sur un dièdre infini. Par transformation de Fourier ce problème se ramène à l'étude d'une famille d'opérateurs à paramètre sur un secteur infini. On calcule le spectre essentiel de ces opérateurs sur le secteur et on montre que dans certains cas il y a des valeurs propres discrètes sous le spectre essentiel. Par comparaison avec des opérateurs de Sturm-Liouville singuliers sur le demi-axe on obtient des majorations du bas du spectre de l'opérateur sur le dièdre : pour un angle d'ouverture assez petit et certaines orientations du champ magnétique, celui-ci est strictement inférieur aux quantités spectrales issues du cas régulier. Finalement on applique ces résultats à l'opérateur de Schrödinger avec champ magnétique constant et petit paramètre dans des domaines bornés de l'espace possédant des arêtes courbes. Pour déterminer une asymptotique de la première valeur propre dans la limite semi-classique, on construit des quasi-modes près de l'arête à l'aide des fonctions propres du problème à paramètre sur le secteur. En utilisant une partition du domaine selon que l'on soit près de l'arête ou du bord régulier, on obtient le premier terme de l'asymptotique pour diverses orientations du champ magnétique et on montre dans certains cas que la première valeur propre est inférieure aux valeurs propres associées à des ouverts réguliers.

Ce travail comporte deux volets. Le premier concerne l'étude de défauts locaux dans des matériaux cristallins. Le chapitre 1 donne un bref panorama des principaux modèles utilisés en chimie quantique pour le calcul de structures électroniques. Dans le chapitre 2, nous présentons un modèle variationnel exact qui permet de décrire les défauts locaux d'un cristal périodique dans le cadre de la théorie de Thomas-Fermi-von Weiszäcker. Celui-ci est justifié à l'aide d'arguments de limite thermodynamique. On montre en particulier que les défauts modélisés par cette théorie ne peuvent pas être chargés électriquement. Les chapitres 3 et 4 de cette thèse traitent du phénomène de pollution spectrale. En effet, lorsqu'un opérateur est discrétisé, il peut apparaître des valeurs propres parasites, qui n'appartiennent pas au spectre de l'opérateur initial. Dans le chapitre 3, nous montrons que des méthodes d'approximation de Galerkin via une discrétisation en éléments finis pour approcher le spectre d'opérateurs de Schrödinger périodiques perturbés sont sujettes au phénomène de pollution spectrale. Par ailleurs, les vecteurs propres associés aux valeurs propres parasites peuvent être interprétés comme des états de surface. Nous prouvons qu'il est possible d'éviter ce problème en utilisant des espaces d'éléments finis augmentés, construits à partir des fonctions de Wannier associées à l'opérateur de Schrödinger périodique non perturbé. On montre également que la méthode dite de supercellule, qui consiste à imposer des conditions limites périodiques sur un domaine de simulation contenant le défaut, ne produit pas de pollution spectrale. Dans le chapitre 4, nous établissons des estimations d'erreur a priori pour la méthode de supercellule. En particulier, nous montrons que l'erreur effectuée décroît exponentiellement vite en fonction de la taille de la supercellule considérée. Un deuxième volet concerne l'étude d'algorithmes gloutons pour résoudre des problèmes de propagation d'incertitudes en grande dimension. Le chapitre 5 de cette thèse présente une introduction aux méthodes numériques classiques utilisées dans le domaine de la propagation d'incertitudes, ainsi qu'aux algorithmes gloutons. Dans le chapitre 6, nous prouvons que ces algorithmes peuvent être appliqués à la minimisation de fonctionnelles d'énergie fortement convexes non linéaires et que leur vitesse de convergence est exponentielle en dimension finie. Nous illustrons ces résultats par la résolution de problèmes de l'obstacle avec incertitudes via une formulation pénalisée

Cette thèse porte sur le calcul des corrections radiatives supersymétriques pour des processus entrants dans le calcul de la densité relique de matière noire, dans le MSSM et le scénario cosmologique standard, ainsi que sur l'influence du choix du schéma de renormalisation du secteur des neutralinos/charginos à partir de la mesure de trois masses physiques. Cette étude a été faite à l'aide d'un programme automatique de calcul à une boucle d'observables physiques dans le MSSM, appelé SloopS. Pour le calcul de la densité relique nous nous sommes penchés sur des scénarios où le candidat supersymétrique le plus étudié, le neutralino, se désintégrait en majoritairement en bosons de jauge. Nous avons couvert les cas où sa masse était de l'ordre de quelques centaines de GeV jusqu'à 2 TeV. Cela a nécessité la prise en compte complète des corrections électrofaibles et fortes, impliquées dans des processus sous-dominants impliquant des quarks. Dans le cas des neutralinos très lourds deux effets importants ont été mis à jour : les amplifications de type Sommerfeld dues aux bosons de jauge massifs et peut-être plus important encore des corrections de type Sudakov.

Ce travail propose des méthodes numériques pour la résolution du problème de Stokes en géométrie complexe pour des fluides non homogènes. Ce modèle décrit l'écoulement d'un fluide très visqueux, incompressible, dont la viscosité n'est pas uniforme mais dépend de la concentration d'un certain agent. D'un point de vue mathématique, il s'agit de résoudre un problème elliptique couplé à une équation de convection-diffusion, qui génèrent une dynamique non linéaire. L'algorithme de résolution est basé sur une discrétisation hybride utilisant une grille et des particules. Des algorithmes à pas fractionnaires permettent de séparer la résolution des différents phénomènes pour profiter des avantages spécifiques à ces discrétisations: méthodes lagrangiennes adaptées à la convection et méthodes eulériennes pour la diffusion. Une méthode de pénalisation permet de gérer efficacement l'interaction entre le fluide et la géométrie mobile du domaine. Une méthode de projection itérative est développée pour ce problème quasi-statique, cela permet d'utiliser des solveurs rapides propices aux calculs en grande dimension. Plusieurs cas tests viennent valider la convergence, la conservation et les performances de l'algorithme en 3D. Ce travail s'inscrit dans le contexte de l'étude de l'écoulement du mucus pulmonaire autour des cellules épithéliales ciliées tapissant les bronches. L'efficacité du transport du mucus, assurant la capture et l'expectoration des agents pathogènes, est étudiée en fonction des paramètres biologiques. D'autres simulations d'un micro-nageur et d'écoulements en milieux poreux complètent cette étude.

Les problèmes de complémentarité interviennent dans de nombreux domaines scientifiques : économie, mécanique des solides, mécanique des fluides. Ce n'est que récemment qu'ils ont commencé d'intéresser les chercheurs étudiant les écoulements et le transport en milieu poreux. Les problèmes de complémentarité sont un cas particulier des inéquations variationnelles. Dans cette thèse, on offre plusieurs contributions aux méthodes numériques pour résoudre les problèmes de complémentarité. Dans la première partie de cette thèse, on étudie les problèmes de complémentarité linéaires 0 6 x ⊥ (Mx+q) > 0 où, x l'inconnue est dans Rn et où les données sont q, un vecteur de Rn, et M, une matrice d'ordre n. L'existence et l'unicité de ce problème est obtenue quand la matrice M est une P-matrice. Une méthode très efficace pour résoudre les problèmes de complémentarité est la méthode de Newton-min, une extension de la méthode de Newton aux problèmes non lisses.Dans cette thèse on montre d'abord, en construisant deux familles de contre-exemples, que la méthode de Newton-min ne converge pas pour la classe des P-matrices, sauf si n= 1 ou 2. Ensuite on caractérise algorithmiquement la classe des P-matrices : c'est la classe des matrices qui sont telles que quel que, soit le vecteur q, l'algorithme de Newton-min ne fait pas de cycle de deux points. Enfin ces résultats de non-convergence nous ont conduit à construire une méthode de globalisation de l'algorithme de Newton-min dont nous avons démontré la convergence globale pour les P-matrices. Des résultats numériques montrent l'efficacité de cet algorithme et sa convergence polynomiale pour les cas considérés. Dans la deuxième partie de cette thèse, nous nous sommes intéressés à un exemple de problème de complémentarité non linéaire concernant les écoulements en milieu poreux. Il s'agit d'un écoulement liquide-gaz à deux composants eau-hydrogène que l'on rencontre dans le cadre de l'étude du stockage des déchets radioactifs en milieu géologique. Nous présentons un modèle mathématique utilisant des conditions de complémentarité non linéaires décrivant ces écoulements. D'une part, nous proposons une méthode de résolution et un solveur pour ce problème. D'autre part, nous présentons les résultats numériques que nous avons obtenus suite à la simulation des cas-tests proposés par l'ANDRA (Agence Nationale pour la gestion des Déchets Radioactifs) et le GNR MoMaS. En particulier, ces résultats montrent l'efficacité de l'algorithme proposé et sa convergence quadratique pour ces cas-tests

Nous étudions le problème de rupture des systèmes de couches minces soumis à contraintes de tension dues aux chargements mécaniques ou à d'autres phénomènes élastiques, associés e.g. à couplages thérmiques où humidité. Dans ces systèmes, chargements homogènes conduisent à la nucléation de fissures interagissantes transverses et de décollement, produisant l'auto-structuration de réseaux de fissures quasi-périodiques et la propagation de patterns complexes qui montrent caractéristiques morphologiques robustes. On s'intéresse à décrire l'évolution de ces fissures, en prenant en compte les phases de nucléation, sélection du trajet de fissure et évolution irreversible en espace et en temps. Les résultats disponibles en littérature se basent sur des modèles phénoménologiques, dépourvus d'une dérivation rigoureuse, et sont limités à des cas géométriquement simples. Dans ces derniers, le problème de nucléation, les mécanismes de sélection du chemin de fissuration et l'évolution non régulière en espace et en temps ne sont pas explorés, à cause des limitations de la théorie classique de la mécanique de la rupture. Nous proposons la dérivation d'une théorie variationnelle asymptotique, bidimensionnelle et globale, à partir d'un problème tridimensionnel d'élasticité fragile dans le cadre de l'approche variationnelle à la mécanique de la rupture, en faisant intervenir une notion de convergence variationnelle. Ensuite, nous introduisons une régularisation du problème faible de rupture par le moyen d'un modèle en gradient d'endommagement, adapté à la solution numérique via la méthode des éléments finis. Le travail proposé permet d'obtenir une compréhension des mécanismes couplés élastiques, de fissuration et décollement; d'établir un modèle asymptotique, réduit et variationnel, valable pour des systèmes de couches minces suffisamment riche pour capturer les mécanismes physiques essentiels; et d'aborder une étude détaillé des expériences numériques qui révèlent les patterns complexes de fissures observés dans les systèmes de couches minces.

Mes recherches portent essentiellement sur l''elaboration de m'ethodes de calcul formel pour l''etude constructive des 'equations diff'erentielles lin'eaires, plus particuli'erement autour de la th'eorie de Galois diff'erentielle. Celles-ci vont du d'eveloppement de la th'eorie sous-jacente aux algorithmes, en incluant leur implantation en Maple. Ces travaux ont en commun une approche exp'erimentale des math'ematiques o'u l'on met l'accent sur l'examen d'exemples les plus pertinents possibles. L''etude d'etaill'ee de cas provenant de la m'ecanique rationnelle ou de la physique th'eorique nourrit en retour le d'eveloppement de th'eories math'ematiques idoines. Mes travaux s'articulent suivant trois grands th'emes interd'ependants : la th'eorie de Galois diff'erentielle effective, ses applications 'a l'int'egrabilit'e de syst'emes hamiltoniens et des applications en physique th'eorique.

L'objet de la thèse est de contribuer à l'étude numérique des lois de conservation hyperboliques avec termes sources, ce qui est motivé par les applications aux équations de Saint-Venant pour les eaux peu profondes. La première partie traite des questions habituelles de l'analyse des approximations numériques des lois de conservation scalaires. On se concentre sur des schémas aux volumes finis semi-discrets, dans le cas général d'un maillage non-uniforme. Pour définir des discrétisations appropriées du terme source, on introduit le formalisme spécifique de la méthode "Upwind Interface Source" et on établit des conditions sur les fonctions numériques telles que le solveur discret préserve les solutions stationnaires. Une définition rigoureuse de consistance est ensuite formulée, adaptée aux "schémas équilibres", pour laquelle on est capable de prouver un théorème de convergence faible de type Lax-Wendroff. La méthode considérée dans un premier temps est essentiellement d'ordre un en espace. Pour améliorer la précision, on développe des approches à haute résolution pour la méthode "Upwind Interface Source" et on montre que celles-ci sont un moyen efficace de dériver des schémas d'ordre plus élevé avec des propriétés convenables. On prouve une estimation d'erreur dans $L^p$, $1\le p < +\infty$, qui est un résultat optimal dans le cas d'un maillage uniforme. On conclut alors que les mêmes taux de convergence $O(h)$ et $O(h^2)$ que pour les systèmes homogènes correspondants sont valables. La deuxième partie présente un schéma numérique pour approcher les équations de Saint-Venant, avec un terme source géométrique, qui vérifie les propriétés théoriques suivantes: il préserve les états stationnaires de l'eau au repos, vérifie une inégalité d'entropie discrète, préserve la positivité de la hauteur de l'eau et reste stable avec des profiles du fond discontinus. Cela est obtenu grâce à une approche cinétique au système; dans ce contexte, on utilise une description formelle du comportement microscopique du système pour définir les flux numériques aux interfaces d'un maillage non-structuré. On utilise aussi le concept de variables conservatives centrées (typique de la méthode des volumes finis) et des termes sources décentrés aux interfaces. Finalement, on présente des simulations numériques du système des équations de Saint-Venant modifiées pour prendre en compte le frottement et la viscosité, afin de retrouver les résultats de certaines études expérimentales. Une application à la modélisation des termes de frottement pour les avalanches de neige est discutée dans l'Appendice.