La grande extension méridienne et la faible profondeur du plateau de Kerguelen en font un obstacle majeur à l'écoulement zonal du Courant Circumpolaire Antarctique (CCA) dans le secteur indien de l'océan austral. Tandis que la majorité du transport du CCA est dévié au nord des îles Kerguelen, le reste ( 50 x 106 m3 s1) doit passer plus au sud, probablement par les passages profonds du Fawn Trough (56°S, 77°E, 2650 m) et du Princess Elizabeth Trough (64°S, 82°E, 3650 m). Pourtant, le détail de la circulation autour du plateau est longtemps resté méconnu, en raison des di cultés à récolter des données dans cette région éloignée de toute route commerciale, d'un climat particulièrement rude et de la présence de la banquise recouvrant la moitié du plateau en hiver. L'objectif de cette thèse est d'améliorer notre connaissance de la circulation moyenne autour du plateau de Kerguelen. Pour cela, un large éventail d'outils allant de l'analyse d'observations récentes à l'utilisation de modèles numériques de circulation a été utilisé. Les données obtenues à l'aide d'éléphants de mer instrumentés sont présentées en détail. Une procédure de calibration et de validation de ces pro ls hydrologiques est proposée. Ce jeu de données, constitué d'un grand nombre de sections à haute résolution spatiale, est ensuite combiné avec d'autres observations plus conventionnelles (in situ et satellites) pour construire un nouveau schéma de circulation dans la région. L'existence d'un courant intense traversant le passage du Fawn Trough est alors mise en évidence. Les observations directes du courant du Fawn Trough pendant la campagne récente TRACK (resp. : Y.-H. Park) con rment l'importance de ce courant, qui concentre plus de 40 x 106 m3 s1, soit 30 % du transport total du CCA. La connaissance de la circulation développée à partir des observations est ensuite utilisée pour valider une simulation au 1/4° de résolution (ORCA025-G70, Barnier et al., 2006). A l'aide de tests de sensibilité sur une con guration régionale de l'Océan Indien Sud, l'importance d'une bonne représentation des caractéristiques des masses d'eau profonde et de la bathymétrie pour obtenir une simulation réaliste est mise en évidence. Finalement, une analyse dynamique de la simulation est menée, montrant comment la circulation résulte d'une combinaison des e ets de la bathymétrie et du vent au travers de la balance de Sverdrup topographique. En forçant le passage d'une partie du CCA au niveau de la divergence antarctique (entre 60°S et 65°S), là où le pompage d'Ekman est maximal, le plateau de Kerguelen favorise une intense remontée vers la surface des isopycnes sud du CCA. Ainsi, paradoxalement, le plateau de Kerguelen pourrait être une des principales zones d'accélération du CCA.

Des aspects variés de la modélisation de l'écoulement glaciaire sont abordés. <br />Le modèle de Le Meur et Vincent (2003) est utilisé à plusieurs reprises. Plusieurs améliorations sont effectuées (conservation de la masse, traitement des épaisseurs négatives, traitement de la glace située en amont de la rimaye).<br />Sur le glacier de St. Sorlin (France) une comparaison entre le champ de bilan issu des mesures et celui issu d'un modèle est effectuée. Le comportement du glacier sous un scénario climatique futur est prédit.<br />Une intercomparaison entre différents types de modèles est effectuée avec pour objectif la détermination du type de modèle le plus approprié en fonction du type de glacier. Les modèles testés sont le modèle SIA de Le Meur et Vincent (2003), le modèle SIA et le modèle d'ordre supérieur de Pattyn (2003) et un modèle Full Stokes (Elmer). Des géométries synthétiques sont utilisées ainsi qu'un cas réel. Les tests synthétiques montrent les limites de l'applicabilité de la SIA. Par contre, une comparaison rapide montre le gain considérable en temps CPU. D'un autre côté, l'augmentation du coût en terme de temps CPU ne s'avère pas très importante lors du passage d'un modèle d'ordre supérieur à un modèle Full Stokes.<br />Les simulations effectuées sur le glacier de St. Sorlin donnent un aperçu des limites de la validité de la SIA sur ce glacier. Même si elle reproduit globalement l'évolution observée du glacier, elle ne reproduit pas correctement le champ de vitesse ni certaines structures de la géométrie.<br />Les simulations sont comparées avec les observations pour la position du front et les vitesses de surface. <br />Un dernier chapitre est consacré au glacier Cotopaxi (Andes).

L'objectif de ce travail est d'estimer les apports d'une topographie dynamique moyenne (MDT) déduite des données gravimétriques GRACE, pour la reconstruction de la circulation océanique du Pacifique Tropical par assimilation de données.<br />Les récentes avancées de notre connaissance du géoïde nous ont conduit à étudier l'impact de l'utilisation d'un signal altimétrique absolu. Un modèle (OPA) de l'Océan Pacifique Tropical, où des observations in-situ et quasi-synoptiques sont disponibles (réseau TAO), et un filtre de Kalman en rang réduit (SEEK) ont été mis au point et utilisés dans différentes configurations. La première suppose une situation pré-GRACE et utilise une MDT artificielle. La deuxième utilise une MDT observée déduite du géoïde GRACE. Conjointement à l'altimétrie, les profils de température TAO sont assimilés.<br />Ce travail montre l'importance d'une bonne référence pour les résidus altimétriques. Le résultat le plus important concerne la capacité du système d'assimilation utilisant GRACE à mieux composer avec des données mixtes: satellites et in-situ. Ici, l'assimilation conjointe d'altimétrie et de données TAO est plus performante grâce à la meilleure compatibilité des données.<br />En outre, une analyse physique, qui considère l'apport de l'assimilation à l'amélioration de la représentation de la dynamique du Pacifique Tropical, a été conduite. L'originalité de ce travail est de montrer comment l'assimilation contribue à améliorer notre compréhension des mécanismes physiques en action dans ce bassin.<br />De manière intéressante et rétrospective, cette analyse révèle aussi une zone (8°N) où les données GRACE semblent avoir des faiblesses qu'il serait judicieux de corriger.

La connaissance et la maîtrise de nombreux régimes transitoires qui apparaissent et se propagent sur les réseaux et dans les matériels qui y sont connectés (transformateurs, parafoudres, ...) est un problème d'actualité. Un transformateur de puissance relié à un réseau moyenne tension ou très haute tension est soumis à différents types de surtension dont la connaissance permet d'évaluer les contraintes auxquelles sont soumises les isolations internes. D'une part, nous avons élaboré un nouveau modèle de transformateur cuirassé valable sur une très large gamme de fréquences, issu du principe de la dualité magnétique électrique afin de connaître son comportement interne (fréquences de résonance, coefficients de surtension). Ce modèle est destiné à être introduit dans des codes de calcul de réseaux, par exemple EMTP. Le modèle prend en compte les pertes dans le cuivre et dans le noyau magnétique. Des méthodes ont été également développées pour un calcul fiable des paramètres du modèle. D'autre part, une méthode d'analyse modale permettant d'illustrer le comportement oscillatoire d'un transformateur a été développée. Cette méthode conduit à un schéma équivalent qui sert à étudier les interactions du transformateur avec le réseau. Dans ces précédents modèles, les résultats calculés de simulation sont en bon accord avec ceux obtenus par la mesure, tant sur modèle réduit que sur transformateur THT réel. Par ailleurs, un algorithme de calcul de réseaux électriques fondé sur I'inversion de la transformée de Laplace a été développé. II permet l'inversion de la transformée de Laplace de plusieurs fonctions à la fois, en y adaptant la méthode nodale. Les résultats de calcul sont en très bon accord avec ceux des mesures et permettent d'espérer une meilleure simulation du système ligne-transformateur.

Dans ce travail nous avons ouvert des voies nouvelles pour la méthode des équations intégrales de frontière appliquée à la résolution des problèmes tridimensionnels de l'électromagnétisme en basse fréquence. Elles permettent de réaliser un maillage adaptatif associé à des résolutions approchées, locales, du système intégral. Dans le premier chapitre nous discutons brièvement les différentes méthodes numériques appliquées à l'électromagnétisme. Le deuxième chapitre est consacré à la méthode des équations intégrales de frontière. Nous y justifions nos choix particuliers (type et ordre des éléments discrétisant les frontières) et décrivons les algorithmes réalisés, par exemple pour la prise en compte des symétries ou périodicités des systèmes étudiés. Dans le troisième chapitre, nos méthodes de calcul des intégrales élémentaires, combinant l'analytique et le numérique, sont décrites et validées. Dans le chapitre quatre, nous présentons un estimateur d'erreur basé sur l'écart entre la valeur interpolée du potentiel et sa valeur calculée par la méthode intégrale elle-même. Pour finir, nous utilisons dans le chapitre cinq cet estimateur d'erreur pour réaliser automatiquement un. maillage adapté à chaque problème traité, en autorisant aussi un maillage « non-conforme ». Une méthode de résolution locale permet un gain important en temps de calcul pendant cette phase de maillage adaptatif. Le maillage final est séparé en domaines, ce qui permet d'améliorer le résultat par itérations entre ces domaines, et sans recourir à une résolution globale, très coûteuse pour des problèmes de grande taille.

Quels sont les mécanismes qui contrôlent l'exhumation des roches métamorphiques de HP(UHP)-BT au niveau des différents types de zones de convergence? Quels sont les conditions indispensables à l'obtention de l'exhumation ? Ces mécanismes sont-ils différents selon la nature des roches impliquées ? Ces questions ont déjà fait l'objet de nombreuses études mais ne sont toujours pas complètement résolues.<br /><br /> Les modèles numériques constituent l'outil idéal pour une telle étude, permettant de réaliser des études paramétriques dont les résultats peuvent être évalués avec les données naturelles. L'étude des résultats P-T-t synthétiques obtenus à partir des pointeurs intégrés au code numérique thermomécanique PARA(O)VOZ, comparés aux données naturelles du métamorphisme pour les Alpes nous a ainsi permis d'étudier les processus d'exhumation lors de la convergence, à la fois en contexte de subduction océanique et continentale. Le cas particulier de l'obduction a été traité dans le cas de l'Oman, pour lequel des estimations P-T complémentaires de haute-résolution ont été spécifiquement réalisées.<br /><br /> Ce travail, par la réalisation de modélisations paramétriques, met en lumière l'influence de paramètres tels que la vitesse de convergence, le taux d'érosion ou encore la nature des lithosphères impliquées et la densité des sédiments sur les processus d'exhumation. Il montre que de bonnes corrélations existent entre les résultats des expériences numériques et les observations de terrain (topographie, morphologie, chemins P-T-t, vitesses d'exhumation).<br /> <br /> S'agissant des mécanismes d'exhumation, envisagés en fonction du contexte de convergence et du type de matériel mis en jeu, cette étude révèle l'importance des forces de volumes responsables des premiers stades rapides d'exhumation, tant au niveau du prisme d'accrétion qu'à l'échelle lithosphérique pour la croûte continentale. Elle confirme, en outre, que les processus d'érosion et de tectonique de surface contrôlent l'exhumation des stades finaux, plus lents. Enfin, elle permet de caractériser plus précisément les problèmes du mécanisme d'exhumation de la croûte océanique qui, une fois éclogitisée, ne peut être exhumée qu'avec l'aide de matériaux moins denses comme la serpentinite ou la croûte continentale.<br /><br /> Un bilan de l'exhumation du matériel sédimentaire, océanique et continental au niveau des zones de convergence est finalement proposé.

Dans la présente étude, la conductivité thermique et le module d'élasticité de revêtementsd'YPSZ élaborés par projection plasma ont été prédits par modélisations numériques 2D et3D de type différences finies et éléments finis.L'influence de la résolution d'image, de la taille et de la valeur du seuil sur les propriétésprédites du revêtement a été étudiée. En outre, les effets de la méthode numérique et du typede condition aux limites ont été étudiés. En particulier, la quantification de l'effet Knudsen(effet de raréfaction) sur le transfert de chaleur à travers une structure poreuse a été réaliséepar modélisation numérique en combinaison avec l'analyse d'image. Les conductivitéseffectives obtenues par modélisation 3D s'avèrent plus élevées que celles obtenues en 2D, etaussi en meilleur accord avec les résultats mesurés. Une corrélation 2D/3D a été trouvéepour la modélisation de la conductivité thermique : cette corrélation permet de prédire lesvaleurs 3D à partir des valeurs calculées en 2D.

En tribologie, la modélisation numérique est aujourd'hui un outil indispensable pour étudier un contact afin de pallier les limites expérimentales. Pour comprendre de mieux en mieux les phénomènes mis en jeu, les modèles ne se situent plus à une seule échelle, mais en font intervenir plusieurs, rendant plus que jamais le concept de triplet tribologique incontournable. Travaillant avec cette philosophie et en se basant sur l'approche Non Smooth Contact Dynamics, dont nous rappelons les grandes lignes, nous proposons de franchir deux cas: proposer des modèles offrant des résultats quantitatifs et mettre en place les premières pièces d'une homogénéisation au niveau du contact (VER). Dans le premier cas, l'étude du couplage éléments finis/éléments discrets au sein d'une même simulation a pour but de proposer des modèles plus "réalistes". Même si l'interface utilisée est déjà présente au coeur du contact et ne va pas évoluer, elle permet de mettre en évidence l'utilisation d'outil de mesure permettant de lier le mouvement des particules aux instabilités dynamiques et permet d'avoir des résultats qualitatifs mais aussi quantitatifs puisque la comparaison avec les taux de contraintes expérimentaux sont en très bonne adéquation. Dans le second cas, le VER sous sollicitations tribologiques est étudié afin d'étendre les techniques d'homogénéisation aux problèmes de contact afin de s'affranchir de la description des interfaces aux grandes échelles en trouvant un moyen d'homogénéiser le comportement hétérogène de l'interface et de le faire dialoguer avec le comportement continu des corps en contact en faisant remonter, dans un sens, des grandeurs moyennées à l'échelle microscopique à l'échelle macroscopique des premiers corps et dans l'autre sens, se servir des données locales à l'échelle macroscopique comme conditions limites à l'échelle microscopique.

Modélisation des équations électromagnétiques 3D couplées avec les circuits d'alimentation électrique et le déplacement de parties mobiles, dans le cas de machines électriques tournantes. La formulation choisie utilise, comme inconnue principale, le potentiel vecteur magnétique nonjaugé, int~rpolé sur des éléments finis d'arête. Le caractère compatible de la formulation est assuré par l'introduction d'un potentiel vecteur électrique pour représenter le courant source, tant pour les inducteurs filaires maillés que pour les inducteurs filaires non-maillés. Le circuit électrique d'alimentation est analysé par la méthode des potentiels électriques intégrés dans le temps et les systèmes électrique et magnétique sont résolus simultanément. La connexion entre les parties mobiles est les parties fixes se fait par une technique d'interpolation des maillages, adaptée aux éléments d'arête. Des méthodes numériques spécifiques ont été développées et implantées, notamment pour l'intégration en pas-à-pas dans le temps du système temporel du premier ordre et pour la résolution du système matriciel résultant. L'ensemble des développements effectués a été validé sur un moteur à réluctance variable.

L'insuffisance cardiaque est une maladie à grande incidence dont le traitement définitif est difficile. Les pompes d'assistance ventriculaire ont été proposées comme thérapie alternative à long terme, mais la technologie est relativement jeune et selon son design, axial ou centrifuge, le dispositif favorise soit l'hémolyse, soit la stagnation de l'écoulement sanguin. Les pompes à écoulement mixte, combinant certaines propriétés des deux types, ont été proposées comme solution intermédiaire. Pour évaluer leurs performances, nous avons effectué des comparaisons numériques entre huit pompes, deux axiales, deux centrifuges, et quatre mixtes, en employant un modèle Windkessel du système cardiovasculaire avec paramètres optimisés pour l'insuffisance cardiaque résolu avec une méthode Radau IIA3, une méthode de résolution de système d'équations différentielles ordinaires L-stable appartenant à la famille des méthodes Runge-Kutta implicites. Nos résultats semblent suggérer que les pompes d'assistance mixtes ne démontrent qu'un léger avantage comparativement aux autres types en terme de performance optimale dans le cas de l'insuffisance cardiaque, mais il faudrait effectuer plus d'essais numériques avec un modèle plus complet, entre autres avec contrôles nerveux implémentés. / Heart failure is a disease with a high incidence rate that is difficult to treat definitively. Ventricular assist devices have been proposed as alternative long-term therapeutic options, but the technology is fairly recent and depending on the design, axial or centrifugal, the device will favor hemolysis or blood flow stagnation. Mixed flow devices, which combine properties of the two previous types, have therefore been proposed as intermediate solutions. To evaluate their performance, we did a numerical comparison between eight pumps (two axial, two centrifugal, and four mixed) using a Windkessel model of the cardiovascular system with parameters optimised for heart failure and solved using a Radau IIA3 method, an L-stable method used to numerically solve systems of ordinary differential equations that belongs to the family of implicit Runge-Kutta methods. Our results suggest that in terms of optimal performance, mixed flow pumps are modest improvements over the other two types, and more tests would have to be performed using a more complete model with, among other modifications, neural control systems implemented.