UPThèses

Ce rapport de thèse est consacré à l'étude de modèles de champ de phase de type Caginalp. Nous considérons ici, deux modèles : le premier étant une généralisation du modèle de champ de phase de Caginalp basée sur une généralisation de la loi de Maxwell-Cattaneo et le second une généralisation provenant de la théorie de la conduction de chaleur introduite par Chen et Gurtin. L'étude du premier modèle est faite aussi bien dans un domaine borné (avec un potentiel régulier puis dans le cas d'un potentiel non régulier), que dans un domaine non borné, en l'occurrence R3. Le second modèle est un problème de champ de phase avec un couplage (linéaire et non linéaire). Tout d'abord, l'existence, l'unicité et la régularité des solutions sont analysées aux moyens d'arguments classiques. Ensuite, l'existence d'ensembles bornés absorbants et compacts attractifs est établie, assurant ainsi l'existence de l'attracteur global. Enfin, dans certains cas, l'existence d'attracteurs exponentiels, ainsi que le comportement spatial des solutions lorsque le domaine spatial est un cylindre semi-infini tri-dimensionnel, sont analysés.

Cette thèse est divisée en deux parties principales : La première partie concerne la déformation plastique d'un matériau contraint. Nous commençons cette partie par une introduction physique sur la dislocation et son rôle dans l'étude de la déformation plastique. Nous exposons ensuite deux types de modélisation de la déformation plastique ce qui nous conduit à deux équations différentielles à retard de Mecking-Lüke-Grilhé. Nous présentons une analyse mathématique complète des deux modèles linéaire et non linéaire. Nous terminons cette partie par des tests numériques et une comparaison des deux modèles. La deuxième partie de la thèse traite l'instabilité de Rayleigh-Plateau. Cette étude porte sur les instabilités de surface d'un pore cylindrique sans contraintes. Nous nous intéressons à une EDP parabolique non linéaire d'ordre quatre, obtenue à partir d'une équation d'évolution des films minces. Le résultat principal est l'existence globale de la solution et la convergence vers la valeur moyenne de la donnée initiale en temps long. L'étude théorique est aussi appuyée comme dans la première partie par une validation numérique.

L'objectif de notre recherche est d'étudier l'existence et la régularité des solutions pour des problèmes de valeurs propres faisant intervenir un opérateur →p-multivoque A : V → P(V*) sur un domaine régulier Ω C Rᶰ. Par l'intermédiaire des N-fonctions, nous construisons un opérateur →p-multivoque de Leray-Lions "fortement monotone" sur un espace d'Orlicz-Sobolev anisotrope. Nous signalons que la formulation théorique des problèmes associés à cet opérateur repose essentiellement sur la notion de sous-différentielle de Clarke, pour cela, nous donnons des nouvelles méthodes variationelles qui correspondent à la résolution de ces problèmes dans le cas "sous-critique" dans lequel la compacité joue un rôle important puis dans le cas critique lorsque nous perdons la compacité. Différentes applications sont données pour illustrer nos résultats abstraits, par exemple, un opérateur anisotrope aux exposants variables et un opérateur avec un poids de type Hardy.

L'objectif de notre recherche est d'étudier l'existence d'attracteurs de dimension finie pour trois modèles d'équations aux dérivées partielles non linéaires : Navier-Stokes, Cahn-Hilliard et Cahn-Hilliard visqueux. En outre, nous étudions ces modèles dans les cas autonome et non autonome. Dans le but de définir les attracteurs pour chaque modèle nous avons démontré l'existence et l'unicité de la solution, puis l'existence d'un ensemble borné absorbant. Nous avons démontré également l'existence d'un attracteur global et d'attracteurs exponentiels dans le cas autonome, ainsi que les attracteurs exponentiels, uniforme, rétrograde et exponentiels rétrogrades dans le cas non autonome. Par ailleurs, nous avons démontré que pour les deux derniers modèles, l'attracteur pullback est unique et de dimension finie. L'existence d'un paramètre ε, dans le modèle de Cahn-Hilliard visqueux, nous a conduit à construire une famille robuste d'attracteurs exponentiels en étudiant la limite ε tend vers zéro.

Dans cette thèse, nous résolvons un problème inverse de type Cauchy associé à l'opérateur biharmonique. Pour des données compatibles, comme ce problème est mal posé au sens d'Hadamard, nous utilisons la méthode de régularisation évanescente. Elle est itérative. Son originalité est de faire intervenir, à chaque itération, un problème d'optimisation bien posé qui dépend d'un terme de régularisation dont les effets se dissipent à la limite du processus itératif. Cette limite est la solution du problème de Cauchy. Pour adapter des algorithmes élaborés pour les problèmes de Cauchy associés au laplacien, nous factorisons le problème initial en deux problèmes inverses de Cauchy pour l'opérateur harmonique. Les résultats principaux sont la convergence de la solution discrète vers la solution continue et l'efficacité de la méthode à gérer numériquement, via les éléments finis, le problème factorisé sur différents domaines, même lorsque les données sont bruitées.

Cette thèse traite des fonctionnelles exponentielles du mouvement brownien et porte en particulier sur des calculs explicites de prix de bonds zérocoupon associés au modèle de taux d'intérêt de Dothan. En utilisant des méthodes de noyaux de la chaleur et de résolution d'équations de Fokker-Planck, nous donnons des formules explicites de densités de probabilités ou de leurs transformées de Laplace. Les différentes formules intégrales obtenues complètent celles de l'article original "On the Term Structure of Interest Rates" (L. U.Dothan). La méthode utilisée est directe et implique notamment une nouvelle représentation intégrale pour le module au carré de la fonction Gamma. Nous étudions ensuite les applications à la physique et aux mathématiques financières des résultats obtenus pour les fonctionnelles périodiques et hyperboliques du mouvement brownien. Nous traitons aussi de calculs de sensibilités d'options par le calcul de Malliavin. Nous donnons des expressions explicites de l'indicateur delta pour des prix d'options asiatiques et des obligations reposant sur des taux courts traités dans la première partie de la thèse.

La notion de l'approximation rationnelle est normalement conçue pour la discrétisation en temps. Dans cette thèse nous mélangeons cette notion avec la notion de la convergence au sens de Kato qui découle d'une discrétisation en espace pour l'équation de transport neutronique. Nous appliquons cette procedure aux schémas d'Euler explicite et implicite, Crank-Nicolson et Prédicateur-Correcteur qui ont le degré de convergence 1,2 et 3 au sens de l'approximation rationnelle. Pour démontrer la convergence nous utiliserons le théorème de Cherno et nous donnons aussi des illustrations numérique pour justifier ces degrés de convergence. Dans le dernier chapitre nous donnons quelques nouvelles généralisations des théorèmes de point fixe de type Schauder et de type Krasnoselskii qui se basent sur la notion de la compacité faible sur des espaces Fréchet ayant la propriété de Dunford- Pettis et sur la notion de la U-équicontraction.

Cette thèse a pour objectif de reconstruire la conductivité isotrope d'un domaine borné à partir des données sur le bord. Pour cela, nous considérons le problème de conductivité inverse et on utilise les méthodes du Dbar de Nachman et de Brown-Uhlmann. A partir des données aux bords, ces méthodes consistent à calculer une fonction complexe dite la transformée scattering, ensuite calculer la conductivité en résolvant une équation (qui contient la transformée scattering) dite l'équation Dbar.Pour des raisons de stabilité, nous approchons les transformées scattering par plusieurs façons et nous étudions l'erreur de ces approximations. Nous montrons la stabilité des méthodes Dbar via les approximations des transformées scattering. Nous étudions en détails le cas des conductivités radiales et nous obtenons des expressions explicites des approximations des transformées scattering. Nous utilisons la méthode de Vainikko pour la résolution numérique de l'équation Dbar et nous introduisons un schéma de point fixe et nous étudions sa convergence. Les résultats numériques obtenus montrent que la méthode est efficace et justifient nos résultats théoriques.

Dans cette thèse, nous étudions d'un point de vue numérique des équations de Cahn-Hilliard non isotrope et non isotherme modélisées d'après une approche de Gurtin. Concernant l'équation de Cahn-Hilliard-Gurtin non isotrope, dont la structure est proche d'un flot de gradient, nous proposons une discrétisation en espace par éléments finis mixtes, et en même temps par le schéma d'Euler implicite. Nous établissons la stabilité du schéma semi-discrétisé en espace et du schéma complètement discrétisé pour une non linéarité polynômiale. Nous établissons également, pour ces mêmes schémas, des estimations d'erreurs optimales en norme H1 et en norme L2. Ces résultats sont illustrés par des simulations numériques en dimension un et deux d'espace, qui permettent d'étudier l'influence des différents paramètres. Concernant le modèle de Cahn-Hilliard-Gurtin non isotherme, pour lequel il n'existe pas de résultat d'existence locale, nous proposons un schéma totalement discrétisé qui se révèle stable en pratique. Des simulations numériques en dimension un d'espace permettent d'observer des comportements asymptotiques proches du cas isotherme.

Nous nous intéressons dans cette thèse à l'étude d'une équation parabolique quasilinéaire dans laquelle la diffusion est paramétrée par la longueur des différentes interactions non locales. Pour ce qui est du problème stationnaire associé, après avoir montré des résultats d'existence, d' unicité et de continuité. Nous présentons ensuite un critère général d' inversibilité dépendant du paramètre, ce critère très important va par la suite nous permettre en exemple d'application, de retrouver des résultats d'inversibilité déjà connus lorsque le paramètre est égale au diamètre du domaine. Nous donnons ensuite un résultat de principe de comparaisons de solutions symétriques radiales et une généralisation du compte du nombre de solutions. Enfin nous donnons quelques applications numériques utilisant une méthode de point fixe et de Newton pour illustrer ces résultats. Pour le problème d'évolution, après avoir montré l'existence d'un attracteur global associé à notre problème, nous démontrons une estimation L∞ de la solution en fonction d'estimations Lq, q > 1 utilisant des itérations de type Moser