Ce travail est consacré à l'étude du couplage des équations de Navier-Stokes et de la chaleur posées dans un domaine polygonal non convexe. Après avoir analysé le comportement singulier des solutions de ces équations près des coins du domaine considéré, nous présentons une formulation mixte duale de ces équations basée sur l'introduction de deux nouvelles inconnues : a qui est le tenseur gradient de la vitesse et le champ vectoriel ^ qui désigne le gradient de la température. Ensuite, en considérant une famille de maillage Th du domaine fi, nous analysons une méthode d'éléments finis mixte duale basée sur cette dernière formulation en utilisant l'élément fini de Raviart-Thomas de plus bas degré pour approximer les nouvelles inconnues cr et <^ sur chaque triangle K de la triangulation Th. Tandis que les variables n, p et T seront approximées par des polynômes de degré zéro sur chaque triangle K. En particulier, nous montrons que l'on peut retrouver l'ordre de convergence quasi-optimal si le maillage est raffiné suivant certaines règles qui sont essentiellement celles introduites par G. Raugel et basées sur le fait que les solutions sont régulières dans des espaces de Sobolev à poids. Nous discutons les aspects d'implémentation de la méthode d'éléments finis mixte duale raffinée pour ces équations en utilisant un algorithme de point fixe combiné à une formulation hybride de deux systèmes issus du découplage du problème discret : l'un correspondant aux équations de Navier-Stokes et l'autre à l'équation de la chaleur. Nos résultats numériques obtenus, en plus de confirmer l'ordre de convergence optimal pour un problème test posé dans un domaine polygonal non convexe, sont tout-à-fait comparables avec ceux existants dans la littérature pour la convection naturelle dans une cavité carrée.

L'idée d'utiliser des bases d'ondelettes dans l'analyse numérique (résolution des équations elliptiques, aux dérivées partielles, intégrales) s'est imposée depuis que ces bases ont fait preuve de leur efficacité dans le traitement du signal. Deux problèmes se posent quant au calcul de la solution dans une base d'ondelettes : - problème 1 : l'étude de la structure de la matrice associée à un noyau K d'un opérateur intégral T dans une base d'ondelettes ; - problème 2 : l'adaptation des techniques de discrétisation de Galerkin aux bases d'ondelettes. Cette thèse contribue à l'étude de ces problèmes par l'introduction d'une nouvelle classe d'opérateurs définis par leur matrice représentative dans une base d'ondelettes et caractérisés par les dérivées fractionnaires de leurs noyaux.

On étudie dans cette thèse l'équation de Szegö sur la droite réelle ainsi que ses perturbations. Cette équation a été introduite il y a quelques années par Gérard et Grellier comme modèle mathématique d'une équation non linéaire totalement non dispersive.L'équation de Szegöapparait naturellement dans l'étude de l'équation de Schrödinger non linéaire (NLS) danscertaines situations sur-critiques où l'on constate un manque de dispersion, par exemplelorsque l'on considère NLS sur le groupe de Heisenberg. Par conséquent, une des motivationsde cette thèse est d'établir des résultats concernant l'équation de Szegö qui pourrontéventuellement être utilisés dans le contexte de l'équation de Schrödinger non linéaire.Le premier résultat de cette thèse est la classification des solitons de l'équation de Szegö.On montre que ce sont tous des fonctions rationnelles ayant un unique pôle qui est simple.De plus, on prouve que les solitons sont orbitalement stables.La propriété la plus remarquable de l'équation de Szegö est le fait qu'elle est complètement intégrable, ce qui permet notamment d'établir une formule explicite de sa solution.Comme applications de cette formule, on obtient les trois résultats suivants. (A) On montreque les solutions fonctions rationnelles génériques se décomposent en une somme de solitonset d'un reste qui est petit lorsque le temps tend vers l'infini. (B) On met en évidence unexemple de solution non générique dont les grandes normes de Sobolev tendent vers l'infiniavec le temps. (C) On détermine des coordonnées action-angle généralisées lorsque l'on restreintl'équation de Szegö à une sous-variété de dimension finie. En particulier, on en déduitqu'une grande partie des trajectoires de cette équation sont des spirales autour de cylindrestoroïdaux.Comme l'équation de Szegö est complètement intégrable, il est ensuite naturel d'étudierses perturbations et d'établir de nouvelles propriétés pour celles-ci à partir des résultatsconnus pour l'équation de Szegö. Une des perturbations de l'équation de Szegö est une équation desondes non linéaire (NLW) de donnée bien préparée.On prouve que si la donnée initiale de NLW est petite et à support dans l'ensemble desfréquences positives, la solution de NLW est alors approximée pour un temps long par lasolution de l'équation de Szegö. Autrement dit, on démontre ainsi que l'équation de Szegöest la première approximation de NLW. On construit ensuite une solution de NLW dont lesgrandes normes de Sobolev augmentent (relativement à la norme de la donnée initiale).Sur le tore T, Gérard et Grellier ont démontré un résultat analogue d'approximation deNLW. On améliore ce résultat en trouvant une approximation plus fine, de deuxième ordre.Dans une dernière partie, on s'intéresse à l'équation de Szegö perturbée par un potentielmultiplicatif petit. On étudie l'interaction de ce potentiel avec les solitons. Plus précisément,on montre que, si la donnée initiale est celle d'un soliton pour l'équation non perturbée, lasolution de l'équation perturbée garde la forme d'un soliton sur un long temps. De plus, ondéduit la dynamique effective, i.e. les équations différentielles satisfaites par les paramètresdu soliton. / In this Ph.D. thesis, we study the Szegö equation on the real lineas well as its perturbations.It was recently introduced by Gérard and Grellier as a toy model of a non-lineartotally non dispersive equation. The Szegö equation appears naturally in the study of thenon-linear Schrödinger equation (NLS) in super-critical situations where dispersion lacks,for example, when one considers NLS on the Heisenberg group. Consequently, one of themotivations of this Ph.D. thesis is fi nding new results for the Szegö equation in hope thatthey could be eventually used in the context of the non-linear Schrödinger equation.Our first result is a classification of the solitons of the Szegö equation. We show thatthey are all rational functions with one simple pole. In addition, we prove the orbitalstability of solitons.The Szegö equation has the remarkable property of being completely integrable. Thisallows us to find an explicit formula for solutions. We obtain three applications of thisformula. (A) We prove soliton resolution for solutions which are generic rational functions.(B) We construct an example of non-generic solution whose high Sobolev norms grow toinfinity over time. (C) We find generalized action-angle variables when restricting the Szegöequation to a finite dimensional sub-manifold. In particular, this yields that most of thetrajectories of the Szegö equation are spirals around toroidal cylinders.Since the Szegö equation is completely integrable, it is natural to study its perturbationsand deduce new properties of such perturbations from the known results for the Szegöequation. One perturbation of the Szegö equation is a non-linear wave equation(NLW) with small initial data.We prove that the Szegö equation is the first order approximation of NLW. More precisely,if an initial condition of NLW is small and supported only on non-negative frequencies, thenthe corresponding solution can be approximated by the solution of the Szegö equation, fora long time. We then construct a solution of NLW whose high Sobolev norms grow.On the torus T, Gérard and Grellier proved an analogous first order approximationresult for NLW. By considerning the second order approximation, we obtain an improvedresult with a smaller error.Lastly, we consider the Szegö equation perturbed by a small multiplicative potential.We study the interaction of this potential with solitons. More precisely, we show that, if theinitial condition is that of a soliton for the unperturbed Szegö equation, then the solutionpreserves the shape of a soliton for a long time. In addition, we prescribe the effectivedynamics, i.e. we derive the differential equations satisfied by the parameters of the soliton.

No description available.

Cette thèse regroupe une série de travaux ayant tous trait à des systèmes hamiltoniens non linéaires spatialement étendus présentant une bifurcation nœud-col. Elle est constituée de deux parties. Nous étudions dans une première partie la transition à la dissipation de systèmes unidimensionnels soumis à un forçage local et régis par des équations de type sine-Gordon ou Schrödinger non linéaire (ESNL). Nous en calculons analytiquement les solutions stationnaires et caractérisons le comportement dynamique au voisinage de celles-ci près de la bifurcation. Lorsque la relation de dispersion des systèmes possède une fréquence de coupure, le comportement dynamique est caractéristique de systèmes hamiltoniens. A contrario, lorsque la relation de dispersion ne possède pas de fréquence de coupure, la dynamique du système se couple avec l'émission d'ondes sonores qui joue le rôle d'un amortissement effectif. Elle devient alors typique de systèmes dissipatifs. En outre, les modes propres temporels du système subissent une délocalisation spatiale. La seconde partie de la thèse concerne l'étude de deux types d'écoulements bidimensionnels de fluides parfaits barotropes autour d'un obstacle : un écoulement décrit par l'ESNL et un écoulement à surface libre dans l'approximation eau peu profonde, où sont pris en compte les effets dispersifs dus aux effets de tension de surface. Lorsque la longueur caractérisant la dispersion des ondes sonores tend vers zéro, ces deux écoulements se réduisent à l'écoulement autour d'un disque d'un fluide eulérien compressible, auquel se superpose une couche limite que nous calculons analytiquement. Par des méthodes de suivi de branches fondés sur des développements pseudo-spectraux, nous calculons le diagramme de bifurcation complet des deux écoulements. En étudiant la dynamique des deux systèmes au-delà de la bifurcation, nous mettons en évidence une émission d'excitations (dans le cas de l'ESNL) dont la nature dépend du rapport de la longueur de cohérence sur la taille de l'obstacle. Dans le cadre de l'écoulement en eau peu profonde, cette émission est remplacée par une singularité à temps fini de démouillage.

Dans cette thèse, nous nous intéressons à la modélisation et la simulation de dispositifs nanoélectroniques innovants. Premièrement, nous dérivons formellement un modèle avec masse effective pour décrire le transport quantique des électrons dans des nanostructures très fortement confinées. Des simulations numériques illustrent l'intérêt du modèle obtenu pour un dispositif simplifié mais déjà significatif. La deuxième partie est consacrée à l'étude du transport non ballistique dans ces mêmes structures confinées. Nous analysons rigoureusement un modèle de drift-diffusion et puis nous décrivons et implémentons une approche de couplage spatial classique-quantique. Enfin, nous modélisons et simulons un nanodispositif de spintronique. Plus précisement, nous étudions le renversement d'aimantation dans un matériau ferromagnétique multi-couches sous l'effet d'un courant de spin.

L'objectif de cette thèse est, d'une part, de fournir des méthodes de planification de mouvements pour les systèmes non-holonomes, et d'autre part, d'étudier la contrôlabilité spectrale pour les équations de Schrödinger linéarisées. Nous avons apporté une double contribution au problème de la planification de mouvements pour les systèmes non-holonomes. Fondé sur la géométrie sous-riemannienne, nous avons conçu un nouvel algorithme qui résout complètement le problème dans un cadre général. Nous avons également proposé une implémentation numérique de la méthode de continuation qui fournit des solutions satisfaisantes au problème de la planification du roulement sur le plan, un exemple classique de systèmes non-holonomes à deux entrées. Nous avons donné des conditions nécessaires et suffisantes de contrôlabilité spectrale en temps fini des équations de Schrödinger linéarisées en dimension 2 et 3. Leur généricité par rapport au domaine a été étudiée par une technique originale basée sur les équations intégrales.

Le sujet de la thèse porte sur l'étude d'une nouvelle classe de modèles de transport quantique: les modèles fluides quantiques issus du principe de minimisation d'entropie. Ces modèles ont été dérivés dans deux articles publiés en 2003 et 2005 par Degond, Méhats et Ringhofer dans Journal of Statistical Physics en adaptant au cadre de la théorie quantique la méthode des moments développée par Levermore dans le cadre classique. Cette méthode consiste à prendre les moments de l'équation de Liouville quantique et à fermer ce système par un équilibre local (ou Maxwellienne quantique) défini comme minimiseur d'une certaine entropie quantique sous contrainte de conservation de certaines quantités physiques comme la masse, le courant, et l'énergie. Le principal intérêt des modèles quantiques ainsi obtenus provient du fait qu'étant macroscopiques, ils sont biens moins coûteux numériquement que des modèles microscopiques comme l'équation de Schrödinger ou l'équation de Wigner, et de plus, ils prennent en compte implicitement des effets de collision bien plus difficiles à modéliser à un niveau microscopique. Le but de cette thèse est donc de proposer des méthodes numériques pour implémenter ces modèles et de les tester sur des dispositifs physiques adéquats.<br />Nous avons donc commencé dans le chapitre I par proposer une discrétisation du plus simple de ces modèles qu'est le modèle de Dérive-Diffusion Quantique sur un domaine fermé. Puis nous avons décidé dans le chapitre II et III d'appliquer ce modèle au transport d'électrons dans les semiconducteurs en choisissant comme dispositif ouvert la diode à effet tunnel résonnant. Ensuite nous nous sommes intéressés au chapitre IV à l'étude et l'implémentation du modèle d'Euler Quantique Isotherme, avant de s'attaquer aux modèles non isothermes dans le chapitre V avec l'étude des modèles d'Hydrodynamique Quantique et de Transport d'Énergie Quantique. Enfin, le chapitre VI s'intéresse à un problème un petit peu différent en proposant un schéma asymptotiquement stable dans la limite semi-classique pour l'équation de Schrödinger écrite dans sa formulation fluide: le système de Madelung.

Cette thèse porte sur la méthode dite « méthode multipôle rapide » qui résout hiérarchiquement le problème à N-corps avec une complexité linéaire pour n'importe quelle précision. Dans le cadre de l'équation de Laplace, nous souhaitons pouvoir traiter efficacement toutes les distributions de particules rencontrées en astrophysique et en dynamique moléculaire.<br /> Nous étudions tout d'abord deux expressions distinctes du principal opérateur (« multipôle-to-local ») ainsi que les bornes d'erreur associées. Pour ces deux expressions, nous présentons une formulation matricielle dont l'implémentation avec des routines BLAS (Basic Linear Algebra Subprograms) permet d'améliorer fortement l'efficacité de calcul. Dans la gamme de précisions qui nous intéresse, cette approche se révèle plus performante que les améliorations existantes (FFT, rotations et ondes planes), pour des distributions uniformes ou non.<br /> Outre une nouvelle structure de données pour l'octree sous-jacent et des contributions algorithmiques à la version adaptative, nous avons aussi efficacement parallélisé notre méthode en mémoire partagée et en mémoire distribuée. Enfin, des comparaisons avec des codes dédiés justifient l'intérêt de notre code pour des simulations en astrophysique.

On considère le problème d'advection-diffusion stationnaire v(∇u, ∇v)+( β*∇u, v) = (f, v) et non stationnaire d/dt (u(t), v) + v(∇u, ∇v)+( β*∇u, v) = (g(t), v), ainsi que le problème d'advection (β*∇u, v) = (f, v) sur un domaine polygonal borné du plan. Le terme de diffusion est approché par des éléments de Crouzeix Raviart et le terme de convection par une méthode upwind sur des volumes barycentriques finis avec un maillage triangulaire. Pour le problème stationnaire d'advection-diffusion, la L²-stabilité (c'est-à-dire indépendante du coefficient de diffusion v) est démontrée pour la solution du problème approché obtenue par cette méthode d'éléments finis et de volumes finis. Pour cela une condition sur la géométrie doit être satisfaite. Des exemples de maillages sont donnés. Toujours avec cette condition géométrique sur le maillage, une inégalité de stabilité (où la discrétisation en temps n'est pas couplée à une condition sur la finesse du maillage) est obtenue pour le cas non-stationnaire. La discrétisation en temps y est faite par un schéma d'Euler implicite. Une majoration de l'erreur, proportionnelle au pas en temps et à la finesse du maillage, est ensuite proposée et exprimée explicitement en fonction des données du problème. Pour le problème d'advection, une approche utilisant la théorie des graphes est utilisée pour obtenir l'existence et l'unicité de la solution, ainsi que le résultat de stabilité. Comme pour la stabilité du problème d'advection-diffusion, une condition géométrique - qui est équivalente pour les points intérieurs du maillage à celle du problème d'advection-diffusion - est nécessaire.