Dans cette thèse, nous nous intéressons à un problème mathématique issu de la modélisation de taux d'érosion maximale dans la stratigraphie géologique. Une contrainte globale sur $\partial_t u$, la dérivée par rapport au temps de la solution, est la principale caractéristique de ce modèle. Ce qui nous amène à considérer une équation non linéaire pseudo-parabolique avec un coefficient de diffusion qui est une fonction non-linéaire de $\partial_t u$. En outre, le problème dégénère de telle sorte de tenir compte implicitement de la contrainte. Nous présentons un résultat de l'existence d'une solution au problème continu. Ensuite, une méthode DgFem (discontinuous Galerkin finite element method) pour son approximation numérique est développée. Notre objectif est d'utiliser les propriétéess d'approximation constante par morceaux pour tenir compte implicitement de la contrainte.

L'objet de cette thèse est d'étudier un problème de la magnétostatique tridimensionnel. On propose trois formulations mixtes couplant une méthode d'éléments finis pour tenir compte du milieu hétérogène et une méthode éléments de frontière pour le milieu extérieur homogène. Pour la méthode intégrale on a utilisé les équations de Calderon, l'opérateur de Neumann-Dirichlet ou d'autres opérateurs intégraux. L'utilisation des éléments d'arête de Nédélec pour le champ magnétique, et les éléments de face de Raviart pour l'induction magnétique permet d'utiliser des méthodes éléments finis conformes. Des résultats numériques ont permis de valider ces méthodes. La deuxième partie a porté sur la comparaison de diverses discrétisations pour l'opérateur de Poincaré-Steklov. Ces méthodes ont été comparées sur une formulation de la magnétostatique. Enfin, on propose des formulations discontinues du problème de la magnétostatique avec des conditions aux limites. On montre que ces formulations sont consistantes et des estimations d'erreur sont obtenues.

Ce travail de thèse a pour objet la mise en oeuvre d'une méthode de calcul des états liés dans les structures à multipuits quantiques. Il participe ainsi à l'amélioration des outils de simulation permettant d'optimiser les composants avant leur réalisation. Nous présentons le modèle physique utilisé ainsi que les différentes méthodes couramment employées pour le calcul de ces états. Une comparaison avec le calcul par éléments finis du premier ordre montre un avantage majeur de notre approche dans des cas limites usuels comme le couplage de deux puits identiques ou le calcul des bandes de valence d'un puits quantique large, ainsi qu'en terme de rapidité. La finalité de ce calcul est l'évaluation du gain matériau, élément de base de la simulation des composants. Ce nouveau module vient compléter le simulateur BCBV dont nous rappelons les principaux modèles tels que celui de dérive-diffusion et du couplage électrique-optique en semi-classique. Cependant, la présence de zones quantiques peut nécessiter une approche par la matrice de densité pour rendre compte, de manière plus précise, des phénomènes de transport. Pour finir, nous tentons de comparer les résultats de la simulation du gain avec des mesures effectuées à partir de lasers de type Fabry-Pérot.

Les travaux présentés concernent l'étude de problèmes d'existence, d'unicité et de régularité de problèmes d'écoulement de fluides géophysiques. On s'intéresse en particulier aux écoulements des masses d'eau à travers des modèles bi-dimensionnels et tri-dimensionnels. Dans une première partie est présentée différentes approches de la simulation de ces écoulements par les méthodes de Galerkin linéaire et non linéaire. Les résultats de convergence de ces méthodes numériques sont présentés dans le cas bidimensionnel. Dans une seconde partie l'impact de l'hypothèse du toit rigide sur le comportement des solutions est étudié. Différents schémas numériques de résolution sont proposés, permettant de traiter efficacement les différentes contraintes globales liées à cette hypothèse. Des résultats de convergence de la méthode sont présentés pour des modèles de shallow water mono-couche et bi-couche. La dernière partie de ces travaux est consacrée à l'étude d'un couplage entre les vibrations de la croûte terrestre et la masse d'eau océanique. Cette approche utilisant un opérateur de plaque pour la modélisation de la croûte terrestre permet d'obtenir des résultats d'existence et d'unicité des solutions. Ce type de modèle est adaptée, entre autre, la modélisation de phénomènes de tsunami.

Nous avons étudié l'influence de fluctuations de grandes échelles sur le seuil de l'instabilité dynamo. Pour cela, nous avons résolu le problème cinématique pour un champ de vitesse hélicoïdal, auquel nous avons rajouté à sa partie stationnaire, une modulation périodique dans le temps, également hélicoïdale. Pour un champ de vitesse hélicoïdal stationnaire des études précédentes ont montré que pour de grands nombres de Reynolds magnétique le champ magnétique était généré au niveau d'une surface caractérisée par une condition de résonance sur le champ de vitesse. Pour un champ de vitesse modulé, nous avons montré que pour une faible amplitude de modulation, c'est la condition de résonance portant sur la partie stationnaire qui gouverne la génération du champ magnétique. Pour une grande amplitude de modulation, c'est la condition de résonance portant sur la modulation qui contrôle la génération du champ magnétique. Dans la plupart des cas et si la condition de résonance est vérifiée pour les deux parties du champ de vitesse, on trouve que le seuil augmente en fonction de l'intensité de la modulation, puis diminue tout en restant supérieur au seuil de l'écoulement stationnaire de même géométrie. Si la condition de résonance n'est pas vérifiée pour la modulation, alors le seuil augmente drastiquement avec son intensité.<br />Cette étude suggère que l'optimisation des expériences dynamo dépend non seulement de la partie stationnaire du champ de vitesse, mais aussi de ses fluctuations de grande échelle. Si ces dernières ne sont pas optimisées alors le seuil dynamo peut augmenter drastiquement, même si celles-ci sont de faible intensité.

Le givrage a été identifié comme un danger important dès le début de l'aéronautique.L'accrétion de givre sur les ailes d'avion, due à la présence de gouttelettes surfonduesdans les nuages, cause parmi d'autres conséquences néfastes une dégradation des performancesaérodynamiques pouvant conduire au décrochage. C'est pourquoi les avionneursdéveloppent depuis longtemps des systèmes de protection. Comme les essais en vols ou ensoufflerie sont souvent complexes à mettre en oeuvre et onéreux, la simulation numériqueest devenue un outil efficace et complémentaire pour dimensionner ces systèmes.Cette thèse s'inscrit dans le contexte de la modélisation de l'accrétion de givre, duruissellement et des systèmes de protection thermique. Elle s'articule en sept chapitres.Après avoir présenté les enjeux et contexte, on introduit une approche tricouche permettantde modéliser l'accrétion de givre et le ruissellement de manière instationnaire. Les troischapitres suivants traitent des méthodes de discrétisation de ce modèle ainsi que de soncouplage avec un modèle du système de protection thermique. Les deux derniers sontconsacrés à la présentation des résultats de simulations numériques montrant l'intérêt del'approche développée et la faisabilité de simulations complètes de phénomènes d'accrétionde givre sur une paroi chauffée ou non. / Icing has been identified as a serious issue since the start of aeronautics. Ice accretion onwings, due to supercooled droplets inside clouds, leads to severe degradation of aerodynamicperformances, among other undesirable effects. Therefore, aircraft manufacturers have sincea long time developed ice protection systems. As flight tests or wind tunnel experimentsare often complicated to implement and expensive, numerical modeling is an effective andcomplementary tool to design those systems. This thesis concerns the modeling of ice accretion, runback and thermal ice protectionsystems. It consists of seven chapters. The first one is dedicated to the presentation of theconcerns and the context. Then a three layer approach allowing to model in an unsteadyway ice accretion and runback is presented. Following three chapters deal with this modeldiscretization as well as a method to couple it with a thermal ice protection system model.Two last chapters are dedicated to numerical simulations showing the sake of the approachand the feasibility of a whole simulation of ice accretion on a heated or unheated surface.

La tomographie optique diffuse de fluorescence permet de développer des systèmes d'imagerie fonctionnelle. Elle permet le suivi in vivo de la biodistribution de marqueurs spécifiques fluorescents, dans le but de développer de nouvelles thérapies, en particulier dans le domaine de l'étude de cancers. Le travail de cette thèse est appliqué au problème mathématique de la tomographie, qui peut être formulé en termes d'évaluation de paramètres physiques dans un ensemble d'équations aux dérivées partielles (EDPs).<br />La Méthode des Eléments Finis (MEF) a été choisie ici pour résoudre ces EDPs, car elle permet de ne pas avoir besoin d'hypothèses sur la géométrie et les inhomogénéités du système. Mais elle consomme beaucoup de temps et de mémoire de calcul, principalement en raison des grandes dimensions des matrices impliquées.<br />Notre objectif est d'accélérer la résolution du système.<br />Pour cela, dans un premier temps, une analyse multirésolution a été choisie : les matrices apparaissant dans la version discrétisée des EDPs sont projetées sur une base orthonormale d'ondelettes de Haar ; on obtient alors des approximations des matrices, de taille réduite, contenant l'essentiel de l'information. En utilisant ces approximations à la place des matrices originales, le temps de calcul pour la résolution des EDPs est alors sensiblement réduit.<br />Dans un deuxième temps, la méthode de Galerkin-ondelettes est choisie à la place de la MEF standard : les polynômes d'interpolation classiquement utilisés par la MEF sont directement remplacés par des fonctions d'ondelettes de Daubechies. On s'affranchit alors de l'étape du maillage, coûteuse en temps de calcul, ce qui permet de diminuer également du temps de calcul tout en gardant une qualité comparable de la solution.

La prise en compte des nuages est essentielle dans l'étude du climat, en particulier leur rôle radiatif. A cause de ses aspects non-linéaire et multiéchelle, une des difficultées du transfert radiatif tient à la prise en compte des hétérogénéitées spatiales des propriétées microphysiques (variabilité horizontale, hétérogénéitée géométrique, ...), tant d'un point de vue théorique que numérique. Afin de répondre à cette problématique, nous utilisons ici l'analyse multirésolution de type ondelette dans une méthode de Galerkin appliquée à l'Equation du Transfert Radiatif. Cette méthode permet de formaliser le problème du transfert radiatif en milieu hétérogène en distinguant les processus élémentaires d'interaction entre le champ de rayonnement et le milieu nuageux, décrits par des coefficients de connexion. Cette réécriture a conduit à la mise au point d'un code de calcul, conçu comme un "modèle instrument", un outil d'analyse des problèmes d'échelle en transfert radiatif. Les résultats ont été validés par rapport aux codes Monte Carlo et SHDOM. Des premières simulations ont montré l'intérêt de cette approche pour quantifier les sources radiatives additionnelles à l'échelle de pixels dits d'approximation dues aux couplages locaux entre les echelles, et étudier les effets radiatifs 3D. Cette nouvelle approche ouvre des possibilités d'application aux études en cours sur les limitations des codes de calcul existants et des techniques d'inversion de mesures satellites, et peut amener à la définition de paramétrisation multichelle des hétérogénéités nuageuses.

La compréhension et l'amélioration de l'environnement vibratoire des charges utiles demande la mise au point de démarches prédictives maîtrisées qui permettent de comprendre les phénomènes de transmission des ondes de chocs d'origine pyrotechnique dans le lanceur Ariane5. Plus particulièrement, la maîtrise du comportement transitoire des coques sandwichs en nid d'abeilles, principaux constituants de l'Adaptateur de Charges Utiles – structure porteuse des satellites, est nécessaire pour prédire les vibrations au pied des équipements électroniques des satellites et des lanceurs. Cette problématique présente un caractère multi-échelle tant d'un point de vue temporel (charge mobile supersonique, temps d'analyse) que spatial (dimensions des structures du lanceur, taille des cellules en nid d'abeilles, longueurs d'ondes liées aux hautes fréquences). Celui-ci a été traité dans cette thèse en s'appuyant d'une part, sur une qualification à la fois analytique et numérique des modèles classiques homogénéisés des plaques sandwichs en nid d'abeilles pour la gamme de fréquence mise en jeu et d'autre part, sur une application des stratégies de remaillage adaptatif pour la propagation des ondes développées dans le cadre de la méthode de Galerkin espace-temps discontinue en temps. Deux catégories de modèles de plaques épaisses ont été ainsi construites dans le but d'enrichir la cinématique classique de plaques épaisses de Mindlin-Reissner qui s'est avérée être insuffisante pour correctement représenter le comportement dynamique hors-plan des plaques sandwich en nid d'abeilles. Ainsi ont été analysés les modèles dits monocouches basés sur un enrichissement de la cinématique par ajout de degrés de liberté dans l'épaisseur, et les modèles multicouches composés d'une superposition de trois plaques avec une homogénéisation séparée des matériaux. Il a été montré que ces deux sortes de modèles améliorent la description des phénomènes de hautes fréquences, notamment ceux de flexion et de cisaillement transverse qui sont plus délicats à retranscrire. Toutes les études numériques ont été effectuées avec un code éléments finis qui emploie des solveurs adaptatifs dynamiques basés sur la méthode de Galerkin espace-temps discontinue en temps. Cette méthode d'intégration en temps introduit un amortissement numérique dépendant du pas de temps et qui peut interférer avec un amortissement physique susceptible d'être introduit dans un modèle numérique et conduire au final à un amortissement total différent de celui qui est attendu. Cette interaction a été analysée et mise en évidence dans ce travail à travers l'introduction de l'amortissement de Rayleigh dans les modèles de propagation de chocs. Les outils et les modèles de propagation ainsi développés ont été validés sur plusieurs structures académiques et industrielles. Des comparaisons avec des données expérimentales sur des structures industrielles de grande taille, plus particulièrement sur un Adaptateur de Charges Utiles d'Ariane5, sont effectuées et soulignent la cohérence de notre approche ainsi que la fiabilité et l'efficacité des modèles de propagation proposés.

Le principe du contrôle par flux de fuite magnétique (Magnetic Flux Leakage MFL) consiste à aimanter la pièce à contrôler par un champ magnétique et à détecter à l'aide d'un capteur magnétique les fuites des lignes du champ qui résultent de la présence d'un défaut dans la pièce. Dans le but d'améliorer les performances d'un dispositif de détection, le CEA et la société Vallourec collaborent pour développer un modèle numérique dédié au contrôle virtuel des défauts longitudinaux dans les tubes ferromagnétiques. Le dispositif expérimental comprend un circuit magnétique tournant à une vitesse constante autour du tube qui défile. Dans le cadre de cette thèse, on débute le problème de la modélisation sans tenir compte des effets de la vitesse de rotation, il s'agit donc de résoudre un problème d'électromagnétisme en régime magnétostatique.Pour résoudre ce problème, on propose de comparer une approche semi-analytique basée sur le formalisme des équations intégrales (EI) et une approche purement numérique utilisant les éléments finis (EF).Dans la première partie de cette thèse, après avoir établi le formalisme théorique par EI, un premier modèle considérant des matériaux ferromagnétiques à perméabilité magnétique constante (régime linéaire) a été mis en œuvre en 2D. Ce modèle a été appliqué pour un exemple de système extrait de la littérature et validé numériquement par une comparaison des résultats EI/EF. Pour une meilleure détection, il est opportun de saturer magnétiquement la pièce. Le matériau ferromagnétique est alors caractérisé par une courbe B(H) non-linéaire. Par conséquent, la deuxième partie de la thèse a été consacrée à la mise en œuvre du modèle en régime non linéaire qui tient compte de cette caractéristique.Différentes méthodes de discrétisation ont été étudiées afin de réduire le nombre d'inconnues et le temps de calcul. L'originalité de la thèse réside dans l'utilisation des fonctions d'interpolation d'ordre élevé (polynôme de Legendre) pour une discrétisation des équations intégrales par une approche de type Galerkin. Les premiers essais de validation numérique de ce modèle ont été effectués sur un système MFL simplifié. Des premiers essais de validation expérimentale pour des données obtenues par EF ont été effectués en deux phases : La première a consisté à vérifier le distribution du champ magnétique pour un tube sain et en régime magnétostatique. La deuxième phase a consisté à calculer la réponse d'un défaut dans le tube ferromagnétique en tenant en compte les effets éventuels de la rotation du circuit magnétique par rapport au tube.