Les écoulements turbulents à phases séparées sont présents dans de très nombreuses applications. Cependant, la simulation de tels écoulements avec une interface déformable constitue l'un des problèmes les plus complexes de la mécanique des fluides numérique. La prise en compte du bilan des contraintes normales est au cœur du problème de déformation de l'interface. Dans le travail présenté ici, nous développons un algorithme permettant de simuler des écoulements diphasiques incompressibles et turbulents en suivant le déplacement de l'interface par une approche lagrangienne. Les équations de Navier-Stokes instationnaires écrites en variables vitesse-pression sont résolues dans les deux phases en utilisant des maillages curvilignes orthogonaux. Dans un premier temps, nous introduisons un schéma de raccordement des vitesses tangentielles et des cisaillements. Ce schéma est appliqué afin de simuler l'interaction de deux écoulements turbulents séparés par une interface plane. La turbulence est traitée par une approche de simulation des grandes échelles utilisant un modèle dynamique. Un algorithme original est ensuite développé dans le but de satisfaire de façon non-itérative à la fois la continuité des vitesses normales et des contraintes normales sur l'interface et l'incompressibilité dans les deux phases. Différentes simulations d'écoulements diphasiques avec interface déformable sont réalisées afin de valider ces développements. / Turbulent incompressible two-phase separated flows are present in many applications. However, simulation of such flows with a moving interface is one of the most challenging problems in todays computational fluid dynamics. Taking properly into account the normal stress budget accross the interface is the main difficulty of moving interface problems. This work deals with the development of a boundary-fitted method for computing turbulent incompressible two-phase flows. The interface displacement is achieved through a Lagrangian approach. The unsteady Navier-Stokes equations written in a velocity- ressure formulation are solved within the two phases using an orthogonal curvilinear grid. In a first step, we introduce a scheme allowing tangential velocities and shear stresses to match across the interface. We apply this technique to compute the countercurrent flow generated by two streams separated by a plane interface. This scheme is then applied to compute various situations involving the interaction between two turbulent flows separated by a flat interface. The turbulence is treated by using the Large Eddy Simulation approach with a dynamic model. An original algorithm is then developed to satisfy without any internal iteration the continuity of normal velocities and stresses across the interface and the incompressibility condition within both phases. Several simulations of two-phase flows with a moving interface are carried out to validate these developments.

La perspective lagrangienne, décrivant un écoulement selon les trajectoires de traceurs fluides, est une approche naturelle pour étudier les phénomènes de dispersion dans les écoulements turbulents. En turbulence de paroi, le mouvement des traceurs est influencé par le cisaillement moyen et par une forte inhomogénéité et anisotropie en proche paroi. On étudie les propriétés lagrangiennes d’un écoulement de canal turbulent par simulation numérique directe à un nombre de Reynolds modéré. Les statistiques d’accélération lagrangienne sont comparées aux expériences de suivi de particules réalisées en parallèle à ce travail. Comme en turbulence homogène isotrope (THI), les composantes d’accélération le long des trajectoires lagrangiennes se décorrèlent sur des temps comparables aux plus petites échelles de l’écoulement, tandis que la norme de l’accélération reste corrélée plus longtemps. La persistance d’anisotropie à petite échelle loin de la paroi est constatée par l’existence d’une corrélation croisée non nulle entredeux composantes de l’accélération. On montre que, en conséquence des flux moyens d’énergie cinétique en turbulence de paroi, près des parois les traceurs se déplacent et s’étalent sur des plus grandes distances quand ils sont suivis en arrière dans le temps qu’en avant. La dispersion relative de paires de traceurs est aussi étudiée. Aux temps courts, la séparation des paires est balistique pour toutes les distances à la paroi. Comme en THI, les traceurs se séparent plus rapidement lorsqu’ils sont suivis en arrière dans le temps. Aux temps plus longs, le cisaillement moyen accélère la séparation dans la direction de l’écoulement moyen. Un modèle de cascade balistique initialement proposé pour la THI est adapté aux écoulements inhomogènes / The Lagrangian perspective, describing a flow from the trajectories of fluid tracers, isa natural framework for studying dispersion phenomena in turbulent flows. In wall-boundedturbulence, the motion of fluid tracers is affected by mean shear and by strong inhomogeneityand anisotropy near walls. We investigate the Lagrangian properties of a turbulent channel flowusing direct numerical simulations at a moderate Reynolds number. Lagrangian accelerationstatistics are compared to particle tracking experiments performed in parallel to this work. Asin homogeneous isotropic turbulence (HIT), the acceleration components along Lagrangianpaths decorrelate over time scales representative of the smallest scales of the flow, while theacceleration norm stays correlated for much longer. The persistence of small-scale anisotropy farfrom the wall is demonstrated in the form of a non-zero cross-correlation between accelerationcomponents. As a result of the average fluxes of kinetic energy in wall turbulence, tracers initiallylocated close to the wall travel and spread over longer distances when tracked backwardsin time than forwards. The relative dispersion of tracer pairs is finally investigated. At shorttimes, pair separation is ballistic for all wall distances. As in HIT, relative dispersion is timeasymmetric, with tracers separating faster when tracked backwards in time. At longer times,mean shear dominates leading to rapid separation in the mean flow direction. A ballisticcascade model previously proposed for HIT is adapted to inhomogeneous flows

Cette thèse s'inscrit dans un projet de modélisation numérique de la dispersion atmosphérique de polluants à travers le code de mécanique des fluides numérique Code_Saturne. L'objectif est de pouvoir simuler la dispersion atmosphérique de polluants en environnement complexe, c'est-à-dire autour de centrales, sites industriels ou en milieu urbain. Dans ce contexte, nous nous concentrons sur la modélisation de la dispersion des polluants à micro-échelle, c'est-à-dire pour des distances de l'ordre de quelques mètres à quelques kilomètres et correspondant à des échelles de temps de l'ordre de quelques dizaines de secondes à quelques dizaines de minutes : on parle de modélisation en champ proche. L’approche suivie dans ces travaux de recherche suit une formulation hybride eulérienne/lagrangienne, où les champs dynamiques moyens relatifs au fluide porteur (pression, vitesse, température, turbulence) sont calculés via une approche eulérienne et sont ensuite fournis au solveur lagrangien. Ce type de formulation est couramment utilisé dans la littérature atmosphérique pour son efficacité numérique. Le modèle lagrangien stochastique considéré dans nos travaux est le Simplified Langevin Model (SLM), développé par Pope (1985,2000). Ce modèle appartient aux méthodes communément appelées méthodes PDF (Probability Density Function), et, à notre connaissance, n'a pas été exploité auparavant dans le contexte de la dispersion atmosphérique. Premièrement, nous montrons que le SLM respecte le critère dit de mélange homogène (Thomson, 1987). Ce critère, essentiel pour juger de la bonne qualité d'un modèle lagrangien stochastique, correspond au fait que si des particules sont initialement uniformément réparties dans un fluide incompressible, alors elles doivent le rester. Nous vérifions le bon respect du critère de mélange homogène pour trois cas de turbulence inhomogène représentatifs d'une large gamme d'applications pratiques : une couche de mélange, un canal plan infini, ainsi qu'un cas de type atmosphérique mettant en jeu un obstacle au sein d'une couche limite neutre. Nous montrons que le bon respect du critère de mélange homogène réside simplement en la bonne introduction du terme de gradient de pression en tant que terme de dérive moyen dans le modèle de Langevin (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Nous discutons parallèlement de l'importance de la consistance entre champs eulériens et lagrangiens dans le cadre de telles formulations hybrides eulériennes/lagrangiennes. Ensuite, nous validons le modèle dans le cas d'un rejet de polluant ponctuel et continu, en conditions de vent uniforme et turbulence homogène. Dans ces conditions, nous disposons en effet d'une solution analytique nous permettant une vérification précise. Nous observons que dans ce cas, le modèle lagrangien discrimine bien les deux différents régimes de diffusion de champ proche et champ lointain, ce qui n'est pas le cas d'un modèle eulérien à viscosité turbulente (Bahlali et al., 2018b).Enfin, nous travaillons sur la validation du modèle sur plusieurs campagnes expérimentales en atmosphère réelle, en tenant compte de la stratification thermique de l'atmosphère et de la présence de bâtiments. Le premier programme expérimental considéré dans nos travaux concerne le site du SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), dans la banlieue sud de Paris, et met en jeu une stratification stable de la couche limite atmosphérique. La seconde campagne étudiée est l'expérience MUST (Mock Urban Setting Test). Réalisée aux Etats-Unis, dans le désert de l'Utah, cette expérience a pour but de représenter une ville idéalisée, au travers d'un ensemble de lignées de conteneurs. Deux rejets ont été simulés et analysés, respectivement en conditions d'atmosphère neutre et stable (Bahlali et al., 2018a) / This Ph.D. thesis is part of a project that aims at modeling pollutant atmospheric dispersion with the Computational Fluid Dynamics code Code_Saturne. The objective is to simulate atmospheric dispersion of pollutants in a complex environment, that is to say around power plants, industrial sites or in urban areas. In this context, the focus is on modeling the dispersion at micro-scale, that is for distances of the order of a few meters to a few kilometers and corresponding to time scales of the order of a few tens of seconds to a few tens of minutes: this is also called the near field area. The approach followed in this thesis follows a hybrid Eulerian/Lagrangian formulation, where the mean dynamical fields relative to the carrier fluid (pressure, velocity, temperature, turbulence) are calculated through an Eulerian approach and are then provided to the Lagrangian solver. This type of formulation is commonly used in the atmospheric literature for its numerical efficiency. The Lagrangian stochastic model considered in our work is the Simplified Langevin Model (SLM), developed by Pope (1985,2000). This model belongs to the methods commonly referred to as PDF (Probability Density Function) methods, and, to our knowledge, has not been used before in the context of atmospheric dispersion. First, we show that the SLM meets the so-called well-mixed criterion (Thomson, 1987). This criterion, essential for any Lagrangian stochastic model to be regarded as acceptable, corresponds to the fact that if particles are initially uniformly distributed in an incompressible fluid, then they must remain so. We check the good respect of the well-mixed criterion for three cases of inhomogeneous turbulence representative of a wide range of practical applications: a mixing layer, an infinite plane channel, and an atmospheric-like case involving an obstacle within a neutral boundary layer. We show that the good respect of the well-mixed criterion lies simply in the good introduction of the pressure gradient term as the mean drift term in the Langevin model (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Also, we discuss the importance of consistency between Eulerian and Lagrangian fields in the framework of such Eulerian/Lagrangian hybrid formulations. Then, we validate the model in the case of continuous point source pollutant dispersion, under uniform wind and homogeneous turbulence. In these conditions, there is an analytical solution allowing a precise verification. We observe that in this case, the Lagrangian model discriminates well the two different near- and far-field diffusion regimes, which is not the case for an Eulerian model based on the eddy-viscosity hypothesis (Bahlali et al., 2018b).Finally, we work on the validation of the model on several experimental campaigns in real atmosphere, taking into account atmospheric thermal stratification and the presence of buildings. The first experimental program considered in our work has been conducted on the `SIRTA' site (Site Instrumental de Recherche par Télédétection Atmosphérique), in the southern suburb of Paris, and involves a stably stratified surface layer. The second campaign studied is the MUST (Mock Urban Setting Test) experiment. Conducted in the United States, in Utah's desert, this experiment aims at representing an idealized city, through several ranges of containers. Two cases are simulated and analyzed, respectively corresponding to neutral and stable atmospheric stratifications (Bahlali et al., 2018a)

Le travail présenté dans ce mémoire concerne l'étude du mouvement d'inclusions solides en suspension dans un écoulement turbulent de canal isotherme et anisotherme par simulation numérique directe. Grâce à cet outil de simulation, nous avons pu analyser l'influence de l'inertie et du croisement de trajectoires sur les caractéristiques dynamiques et thermiques de la phase dispersée ainsi que sur celles du fluide vu au sein d'une turbulence non homogène. Cela pourra notamment faciliter le développement futur des modélisations euléro-lagrangienne et euléro-eulérienne dans ce type de turbulence. Nous avons essayé d'examiner au mieux la conséquence de ces effets sur les statistiques de la phase dispersée telles que la concentration, la moyenne et l'écart type de la vitesse, les covariances fluide-particules, les corrélations triples de la vitesse des particules. Parallèlement à cela, les caractéristiques du fluide vu par les particules ont été étudiées, nous nous sommes intéressés à la vitesse de dérive, aux tensions de Reynolds du fluide vu, et la décorrélation temporelle des fluctuations de la vitesse du fluide vu. En ce qui concerne la partie thermique, nous présentons les statistiques thermiques de la phase dispersée et du fluide vu par les particules au sein de l'écoulement anisotherme vertical descendant. Tous ces aspects sont développés en gardant à l'esprit le cadre général de cette étude, c'est-à-dire comprendre le comportement thermique de la phase dispersée et fournir des informations concernant des grandeurs importantes intervenant dans la modélisation de tels écoulements.

Ce manuscrit comporte trois parties distinctes, la première concerne la simulation numérique d'écoulements de fluides. La deuxième partie porte sur les équations de l'électromagnétisme et leur couplage avec les équations cinétiques de Vlasov dans le cadre de simulations numériques en physique des plasmas. La dernière partie évoque rapidement des travaux qui ont donné lieu à des publications mais qui ne rentrent pas complètement dans un des deux cadres abordés précédemment. Dans la première partie, l'objectif est d'étendre aux maillages triangulaires non structurés une méthode numérique éprouvée pour résoudre les équations de Stokes bidimensionnelles : la méthode Marker And Cell qui a été développée sur des maillages en quadrilatères quasi-réguliers dans les années 60. L'idée proposée pour cela est de résoudre le problème de Stokes avec pour variables le tourbillon, la vitesse et la pression. Alors que les résultats numériques obtenus sur des maillages réguliers sont satisfaisants, ceux sur des maillages non structurés ne le sont pas. Il s'est avéré lors de l'étude théorique que ce problème est un problème de stabilité. On montre théoriquement et numériquement que la formulation tourbillon-vitesse-pression est une généralisation de la formulation fonction courant-tourbillon permettant la prise en compte de conditions limites plus générales. L'instabilité, due à des fonctions harmoniques discrètes, peut être levée en utilisant de véritables fonctions harmoniques, calculées à l'aide de leur représentation intégrale, dans le schéma numérique. On résout ainsi l'instabilité de la formulation fonction courant-tourbillon et on améliore les précédents résultats de convergence connus de cette formulation. En particulier, on démontre alors une convergence en moyenne quadratique du tourbillon de l'ordre de 3/2 à 2 dans les cas les plus réguliers (contre un-demi avant). Puis on utilise le fait que la formulation tourbillon-vitesse-pression est équivalente à la formulation fonction courant-tourbillon pour redéfinir une nouvelle formulation tourbillon-vitesse-pression. En effet, il est bien connu que la formulation classique en fonction courant-tourbillon (n'utilisée qu'en deux dimensions d'espace) est mal posée lorsque l'on cherche le tourbillon dans l'espace de Sobolev H1 car son gradient n'est alors pas contrôlé, mais bien posée dans un autre espace de fonctions moins régulières. On étend alors ce résultat au cas tri-dimensionnel et l'on obtient une nouvelle formulation en tourbillon-vitesse-pression bien posée dans un nouvel espace. On démontre aussi théoriquement que ce nouvel espace est bien celui introduit en 2D, ce qu'on confirme par des résultats numériques. La deuxième partie concerne la résolution d'équations cinétiques intervenant dans la simulation directe des plasmas et des faisceaux de particules chargées (modèles de Vlasov-Poisson ou Vlasov-Maxwell). Grâce à l'augmentation de la puissance de calculs des ordinateurs, la simulation de l'évolution des plasmas et des faisceaux de particules basée sur une résolution directe de l'équation de Vlasov sur un maillage de l'espace des phases devient une alternative aux méthodes particulaires (Particle In Cell) habituellement employées. La force de ces simulations directes réside dans le fait qu'elles ne sont pas bruitées (contrairement aux méthodes PIC) et que l'approximation est de même résolution sur tout l'espace des phases, en particulier dans les régions à faible densité de particules où des phénomènes physiques importants ont lieu. L'inconvénient principal est que beaucoup de points sont inutiles car la fonction de distribution des particules y est nulle, ce qui rend ces méthodes directes coûteuses en temps de calcul. On introduit alors une méthode de résolution directe de l'équation de Vlasov sur un maillage {\bf{mobile}} de l'espace des phases. Ce qui permet de ne mailler que la partie de l'espace des phases sur laquelle la fonction de distribution des particules est {\emph{a priori}} non nulle. Nous avons utilisé avec succès cette méthode de maillage mobile en 1D afin de simuler un problème d'interaction laser-plasma. Nous introduisons donc un maillage mobile de l'espace des phases qui suit parfaitement le développement des instabilités et permet de réduire drastiquement le temps de calcul. Nous avons aussi obtenu les premiers résultats d'une méthode de maillage mobile en 4D en couplant le maillage mobile et une méthode de décomposition de domaines. \\ En ce qui concerne la résolution des équations de Vlasov-Maxwell : nous travaillons sur le développement d'un solveur Maxwell éléments finis d'arêtes d'ordre élevé couplé à un code PIC en trois dimensions d'espaces (6D de l'espace des phases). Un point important afin que le couplage fonctionne est que l'équation de conservation de la charge doit être vérifiée au niveau numérique à chaque pas de temps. Ce qui implique que le courant obtenu à partir de l'évolution de l'équation de Vlasov doit être calculé d'une façon bien particulière. Les premiers résultats obtenus en 2D confirment que la méthode de calcul du courant proposée conserve bien la charge comme attendu.

Une méthode est proposée pour determiner la conductivité hydraulique, en fonction de la concentration en eau, de milieux poreux partiellement saturés susceptibles de se déformer. Elle est fondée sur la description eulerienne des transferts d'eau et de particules solides. elle s'appuie sur l'analyse en régime transitoire des profils d'humidité et de masse volumique séche obtenus par spectrométrie gamma double-sources, specialement developpée à cet effet, couplée à des mesures de pression capillaire par tensiométrie. les résultats présentés sont relatifs à des expériences d'infiltration verticale d'eau dans des matériaux poreux compactes et gonflant librement. le bon accord, entre les valeurs déterminées par l'approche lagrangienne, également mise en oeuvre, et celles obtenues dans le cadre de la description eulerienne confirme la pertinence des hypothèses sous-tendant cette dernière. de plus, elle apparait séduisante dans la mesure où elle peut constituer une formulation générale des écoulements d'eau, incluant les milieux déformables et rigides.

Le bloc opératoire constitue l'un des secteurs les plus coûteux et les plus importants dans un établissement hospitalier. Afin, d'utiliser de manière efficace et rationnelle les ressources (humaines et matérielles) disponibles tout en assurant une bonne qualité de service vis-à-vis des patients, la planification du bloc opératoire est devenue l'une des premières préoccupations des établissements hospitaliers. Le problème de planification du bloc opératoire consiste à déterminer pour un horizon de plusieurs jours, l'ensemble d'interventions qui seront réalisées dans chaque salle opératoire. Ce problème a été traité dans la littérature et plusieurs approches de planification ont été proposées. Toutefois, les approches existantes sont essentiellement basées sur des modèles déterministes qui font abstraction de toutes sortes d'aléas. Or, le bloc opératoire est sujet à nombreuses formes d'aléas qui concernent essentiellement la chirurgie d'urgence et les durées d'interventions. À ce titre, l'objectif de notre travail de recherche est de développer des modèles et méthodes pour la planification des activités chirurgicales dans le bloc opératoire tout en tenant compte des aléas importants. Dans le cadre de cette thèse, nous avons proposé des modèles stochastiques qui (i) capturent les éléments essentiels à considérer lors de la planification des activités chirurgicales, (ii) permettent de modéliser de manière explicite différentes formes d'aléas tels que les incertitudes relatives aux durées des interventions et à la chirurgie d'urgence, et (iii) peuvent être facilement étendus pour modéliser d'autres contraintes dues aux pratiques du terrain. Nous avons aussi développé plusieurs approches et méthodes d'optimisation complètes et complémentaires pour la planification stochastique du bloc opératoire. Moyennant des expérimentations numériques, nous avons évalué les performances de ces différentes approches et nous avons montré leurs efficacités. En particulier, nous avons montré que des gains considérables peuvent être réalisés moyennant une modélisation stochastique du problème de planification du bloc opératoire.

L'émergence des transformateurs piézoélectriques coïncident avec le développement dans les années 1950 des céramiques ferroélectriques appartenant à la famille cristalline des pérovskites qui n'ont cessé de s'améliorer depuis. Outre la compacité dont bénéficie ces structures, les transformateurs piézoélectriques offrent des performances remarquables en terme de gain en tension et rendement utiles pour des applications nécessitant une adaptation de tension ou une isolation galvanique, parfaitement dédiés aux applications de faibles puissances à haut rendement. Toutefois, les transformateurs piézoélectriques peuvent être déviés de leurs applications premières. En effet, la dernière décennie a été marquée par l'apparition de générateur de plasma par effet piézoélectrique utilisant principalement des architectures de type transformateur. Pourtant, si quelques applications usuelles illustrent parfaitement cette interaction, la compréhension des phénomènes physiques qui en sont à l'origine reste à approfondir. L'objectif de cette thèse est d'en expliquer les fondements par une approche méthodique. Ce travail s'articule autour de plusieurs étapes comprenant la mise en oeuvre d'une méthode systématique de la modélisation analytique d'un transformateur piézoélectrique, de l'étude de la carte de champ produit par un transformateur ainsi qu'une étude expérimentale vue des bornes en guise de premières investigations. La modélisation analytique est basée sur l'exploitation du Principe de Moindre Action (PMA). A partir de la théorie linéaire de la piézoélectricité, un modèle général applicable à toutes les géométries de transformateur, exploitant des modes de couplage piézoélectrique multiples, est proposé. Son caractère multimodal est par ailleurs mis en exergue. Cette modélisation est appliquée à une structure classique de transformateur piézoélectrique de type Rosen et les résultats obtenus sont validés d'une part par une identification numérique, issue d'un logiciel de calcul par éléments finis, et d'autre part par une caractérisation expérimentale. La modélisation analytique précédente ne tient pas compte dans sa mise en oeuvre de l'influence de l'environnement dans lequel évolue le transformateur piézoélectrique. Afin de caractériser le potentiel électrique produit, un modèle numérique 2D du champ électrique environnant est proposé selon la méthode des différences finies. Ce modèle est basé sur une extension du modèle analytique précédemment développé incluant les pertes mécaniques afin de quantifier le potentiel électrique de surface. Même si l'influence du plasma est négligée en première approximation, la modélisation permet de mettre en lumière les zones de fort champ correspondant aux zones de décharges luminescentes observées expérimentalement. Finalement, afin de valider le concept de générateur de plasma piézoélectrique, une caractérisation vue des bornes du transformateur piézoélectrique de type Rosen a été entreprise. Une étude systématique du déclenchement de la décharge plasma en fonction du niveau de tension et de la pression environnante a été menée. Cette part expérimentale de l'étude constitue une approche pionnière pour qualifier le comportement électromécanique du transformateur et a ainsi permis de mettre en évidence des comportements non linéaires issus de ce mode de fonctionnement atypique qu'est la génération de décharges de surface par effet piézoélectrique.

L'objectif de cette thèse sur la modélisation des décanteurs automobiles est de mettre au point un outil de simulation numérique permettant de se substituer à un banc organe pour déterminer l'efficacité des décanteurs à l'étape de conception. Après avoir étudié la composition du blow-by sur banc moteur et cela pour différents points de fonctionnement du moteur, nous avons pu comparer la pertinence des différents outils de mesure granulométrique choisis (PCS, ELPI, TEOM) en fonction du spectre granulométrique émis par le moteur. Le PCS a été retenu comme l'outil le plus adéquat en fonction des tailles de gouttes sur lesquelles nous souhaitions avoir des informations. En parallèle, nous avons mené des investigations pour mettre au point une méthodologie de calcul permettant de simuler l'efficacité du décanteur. Le suivi des particules d'huile s'effectue grâce à une simulation lagrangienne. Les forces à prendre en compte ont été déterminées grâce à des expériences numériques simples permettant de verifier si les différentes expressions de ces forces trouvées dans la littérature étaient compatibles des conditions de calculs rencontrées dans les décanteurs. Une autre question soulevée par l'utilisation d'une modélisation lagrangienne dans la littérature, est celle de la prise en compte ou non de la dispersion turbulente. L'étude proposée dans cette thèse nous a permis de montrer que l'implémentation d'un modèle de dispersion turbulente rudimentaire ne permet pas d'obtenir des résultats satisfaisants. La faisabilité de la mise au point d'une méthodologie de calcul fiable de l'efficacité du décanteur repose finalement sur la modification du type de maillage utilisé ainsi que le traitement du calcul de l'écoulement à la paroi. Nous avons montré dans cette thèse que la mise en place d'un maillage raffiné aux parois reposant sur des prismes avec un Y+ =1, associé à une loi de parois à deux couches permettait d'obtenir des résultats satisfaisants en terme de prédiction de l'efficacité du décanteur.

C'est sans doute en Arctique que le changement climatique est le plus visible, et semble affecter toutes les composantes du système Arctique, et notamment ses bilans d'eau douce et de chaleur. Alors que l'on s'attend à ce que le signal d'un changement local en Arctique ait son impact climatique le plus important lorsqu'il est exporté au sud d'un coté ou de l'autre du Groenland vers les mers subarctiques, où il peut moduler l'intensité de la circulation thermohaline, l'objectif de cette thèse est donc d'étudier les échanges de volume, de chaleur, d'eau douce et de glace de l'Océan Arctique vers l'Atlantique Nord. Tout d'abord, une simulation réaliste des années 1958 à 2002 basée sur un modèle global couplé glace/océan est utilisée pour étudier la variabilité du bilan d'eau douce en Arctique, afin de comprendre quelle composante de ce bilan contrôle les variations du contenu halin du bassin. On s'intéresse également à la variabilité des exports d'eau douce vers l'Atlantique Nord, et on montre que les exports d'eau douce vers l'Atlantique sont contrôlés par des mécanismes différents de part et d'autre du Groenland: dans les détroits canadiens, le transport de volume domine la variabilité, alors que salinité et courants contribuent à la variabilité dans le détroit de Fram. Par la suite, une réanalyse océanique des années récentes nous permet d'explorer les conséquences pour le contenu halin de l'Arctique du minimum record de l'extension de glace de l'été 2007. Une méthode numérique originale est ensuite utilisée pour comprendre l'origine des masses d'eau qui sont exportées de l'Arctique vers l'Atlantique Nord. On effectue ainsi une analyse lagrangienne qualitative à partir des sorties 3D mensuelles climatologiques d'un modèle global couplé glace/océan à haute résolution, qui permet de quantifier les contributions relatives des différentes branches de circulation à ces exports, ainsi que les échelles de temps et les transformations de masses d'eau associées. Un schéma complet de la circulation dans le bassin Arctique est ainsi proposé, et nous soulignons le rôle clé de la mer de Barents pour les transformations des eaux d'origine Atlantique. Enfin, nous examinons l'influence relative des différents forçages atmosphériques (vent, flux de chaleur et halins) sur les variations du volume de glace en Arctique. Des expériences de sensibilité sont réalisées à l'aide d'un modèle régional de l'Arctique et l'Atlantique Nord, permettant de mieux comprendre les distributions spatiales et temporelles des contributions des différents forçages atmosphériques.