Academic literature on the topic 'EDPs elliptiques non linéaires'

Les travaux présentés dans ce mémoire portent sur des problèmes d'équations ou de systèmes d'équations aux dérivées partielles (EDPs) elliptiques non linéaires. Ils apparaissent comme des équations d'Euler-Lagrange de problèmes de minimisation sous contrainte avec perte de compacité à l'infini. Ces problèmes sont de plus tous liés à des modèles de physique : strucure électronique des solides et (hyper)-élasticité non linéaire (chapitres 1,2 et 3 d'une part, et condensats de Bose-Einstein (chapitre 4) d'autre part.<br /><br />La base de travail des chapitres 1, 2 et 3 est le modèle de Thomas-Fermi-von Weizsäcker (TFW), ou certaines de ses extensions. Dans ce modèle, un système moléculaire est décrit par N noyaux, qui sont des particules classiques ponctuelles, et N électrons, qui sont des particules quantiques définies par leur densité collective. L'énergie TFW, qui dépend des positions des noyaux et de la densité électronique, est minimisée par rapport à cette dernière. Ce modèle est défini au départ pour un nombre fini de noyaux et d'électrons, et sa définition pour une infinité de particules est un problème non trivial. Ce problème, dit de limite thermodynamique, consiste à faire tendre conjointement le nombre de noyaux et d'électrons vers l'infini, en imposant une certaine géométrie (typiquement la périodicité) aux noyaux, et à obtenir la convergence de la densité d'électrons, ainsi que de l'énergie moyenne du système. Ce problème a été résolu dans le cas périodique par I. Catto, C. Le Bris et P.-L. Lions.<br /><br />Le chapitre 1 aborde le problème de la justification de la périodicité supposée dans l'ouvrage de Catto, Le Bris et Lions. Dans la section 1.3, on considère l'énergie TFW d'un cristal comme une fonction du réseau périodique définissant la position des noyaux, et on étudie l'existence d'un minimiseur. Un préliminaire à ce travail, présenté dans la section 1.2, est l'étude des cas dégénérés de réseaux périodiques, à savoir le cas où les noyaux sont répartis périodiquement sur un plan d'une part, et celui où les noyaux sont répartis périodiquement sur une droite d'autre part.<br /><br />Les sections 1.4 et 1.5 abordent le problème sans supposer la périodicité : on minimise l'énergie TFW par rapport à la densité électronque et par rapport à la position des noyaux, à N fixé, et on démontre alors que quand N tend vers l'infini, la configuration minimisante devient périodique. Ce problème est traité théoriquement pour le cas 1D (section 1.4), puis une étude numérique est faite sur le cas 2D (section 1.5), indiquant que le résultat est aussi vrai dans ce cas.<br /><br />Bien que la périodicité soit une bonne approximation pour les cristaux simples, il arrive souvent (dans le cas des polycristaux, des solides amorphes ou de solides cristallins présentant des dislocations par exemple) que cette hypothèse ne soit pas valable. C'est pourquoi on étudie dans le chapitre 2 les problèmes de définition du modèle TFW, pour des solides dont les positions de noyaux ne sont pas périodiques. Un cas déterministe est présenté dans la section 2.1.1, où l'on construit le cadre fonctionnel nécessaire à la définition du modèle, puis on résout le problème de limite thermodynamique associé. La section 2.1.2 présente un cas où les positions des noyaux sont stochastiques. Là aussi, on commence par construire un cadre stochastique (stationnaire ergodique) nécessaire, puis on résout le problème de limite thermodynamique correspondant.<br /><br />Outre ces problèmes de limite thermodynamique, qui font le lien entre un modèle moléculaire et le modèle de théorie des solides correspondant, on étudie dans la section 2.2 des modèles (dits "orbital-free'') plus élaborés utilisés dans certains codes de chimie, sans chercher à les justifier par limite thermodynamique. Cette étude montre que le problème variationnel est mal posé, et que le "minimum'' calculé est un minimum local vraisemblablement dépendant de la discrétisation utilisée et du point de départ de l'algorithme de minimisation.<br /><br />Le modèle TFW est un modèle microscopique. Il est cependant naturel, après l'avoir défini pour des solides (cristallins ou non), d'étudier le lien de ce modèle avec des modèles d'élasticité non linéaire. Ce problème est évoqué dans le chapitre 3, où on considère l'énergie d'un système atomique déformé par un diffémorphisme u, et on passe à la limite quand la distance inter-atomique tend vers 0. On obtient ainsi une énergie hyperélastique qui a la forme de celles utilisées en mécanique. La section 3.1 présente ce travail dans un cadre déterministe, la section 3.2 le même type de résultat dans le cas où les positions des noyaux sont stochastiques.<br /><br />La section 3.3 présente une étude similaire, mais dans le cas d'un joint collé, c'est-à-dire d'une interface d'épaisseur nulle au niveau macroscopique (mais infinie au niveau microscopique). Ce cas est particulier car il doit autoriser un saut de la déformation à travers l'interface, ce qui lui impose une régularité moindre que précédemment.<br /><br />Dans le même esprit, la section 3.4 présente l'analyse du couplage entre un modèle de mécanique des milieux continus et le modèle discret correspondant. L'idée est ici d'étudier la déformation d'un solide qui est régulière dans une partie du solide, mais présente des singularités. Là où la déformation est régulière, on utilise un modèle d'élasticité standard, et là où la déformation est singulière, on revient au modèle discret mettant en jeu les atomes et leurs interactions. Comme à notre connaissance aucune étude théorique n'existait sur ce type de théorie, nous avons étudié un cas très simple de dimension 1, et obtenu des résultats qui laissent penser que le modèle est "bon'' dans le cas convexe (i.e si le potentiel d'interaction des atomes est convexe), mais beaucoup plus douteux dans le cas contraire.<br /><br />Le chapitre 4 présente des travaux sur les condensats de Bose-Einstein. La première section porte sur l'écoulement d'un condensat autour d'un obstacle (physiquement, un laser). Nous établissons l'existence d'une solution sans vortex si la vitesse de translation de l'obstacle est suffisamment faible. Ce résultat avait déjà été établi pour un modèle de dimension 2, et nous l'avons étendu au cas plus réaliste de dimension 3, en étudiant en particulier la zone du bord du condensat où le modèle 2D n'est pas valable (contrairement au coeur du condensat).<br /><br />La section 4.3 concerne l'étude de condensats en rotation, et en particulier des vortex nucléés par cette rotation. Les résultats présentés portent sur la rotation rapide : si Omega est la vitesse de rotation, le système n'a de minimum d'énergie que si Omega < 1. La rotation rapide correspond à la limite Omega tend vers 1. Dans ce régime, la fonction d'onde peut être approximée avec une bonne précision par une fonction analytique multipliée par une gaussienne. Les vortex sont alors les zéros de cette fonction. Nous établissons une borne supérieure de l'énergie en utilisant une fonction test dont les zéros forment un réseau distordu sur les bords du condensat. Ceci est en accord avec les observations expérimentales et numériques.