Les jets sont des phénomènes d’éjection collimatée de plasma magnétisé. Ces ph́énomènes liés à l’accrétion d’un disque sur un objet central, sont relativement répandus dans l’univers : les environnement des étoiles jeunes (objets Herbig-Haro, étoiles T Tauri), des binaires X, des sursauts gamma et les noyaux actifs de galaxies... Les jets extra-galactiques sont issus des trous noirs super-massifs au centre de galaxies telles que les quasars ou les radiogalaxies. Ils sont caractérisés par leur taille, leur puissance et la vitesse du plasma.Les jets extragalactiques sont étudiés dans de ce travail de thèse, même si les outils et méthodes développés peuvent être utilisés pour les binaires X et les micro-quasars. Nous poserons en particulier les questions des mécanismes de lancement, d’accélération et de collimation de ces écoulements. Nous traiterons également de la source énergétique à l’origine de l’écoulement qui peut atteindre une puissance de l’ordre de 10^47 erg.s−1.Le liens avec l’accrétion, la proximité de la base des jets avec le trou noir central, les vitesses d’écoulement observées dans certains jets, montrent que le traitement de ces questions doit inclure les effets de la relativité générale. Nous étudierons donc des solutions de la décomposition 3+1 des équations de la magnéto-hydrodynamique en métrique de Kerr. Nous nous appliquerons au développement d’un modèle d’écoulement méridional auto-similaire avec un traitement consistant du cylindre de lumière. Ce modèle pouvant s’appliquer à la fois au jet et à l’accrétion. Nous explorons les mécanismes d’accélération et de collimation des solutions produites. Nous calculerons des solutions de l’écoulement entrant dans l’horizon et de l’écoulement sortant à l’infini incluant des termes d’injection de paires. Le rôle du mécanisme de création de paires et des processus d’extraction de l’énergie du trou noir sera exploré. / Jets are collimated ejection phenomena of magnetized plasma. These phenomena related to the accretion of a disk on a central object, are relatively common in the universe: the environment of young stars (Herbig- Haro Objects, T Tauri stars...), X-ray binaries, Gamma-ray-bursts, and active galactic nuclei... Extragalactic jets come from super-massive black holes in the center of galaxies such as quasars or radiogalaxies. They are characterized by their size, their power and velecity of the plasma.Extragalactic jets will be the subject of studies in this thesis work, although the tools and methods developed can be used for X-ray binaries and microquasars. In particular, we will ask questions about the mechanisms of launching, accelerating and collimating these flows, but also about the energy source at the origin of the flow that can reach a power in the order of 10^47erg.s−1.The links with the accretion, the proximity of the jet base to the central black hole, flow velocities observed in some jets, show that the treatment of these issues must include the effects of general relativity. We will therefore study solutions of the 3+1 decomposition of magneto-hydrodynamic equations in Kerr metric. We will apply ourselves the development of a meridional self-similar magnetized flow model with a consistent treatment of the light cylinder effect. This model can be applied to both spine jet and accretion. We explore the mechanisms of acceleration and collimation of the obtained solutions. We will calculate solutions of the incoming flow in the horizon and the outgoing flow reaching infinity including injection terms. The role of the pair creation mechanism and the processes of extracting energy from the black hole are explored.

Une zone de constriction dans un micro-canal fluidique peut générer, si le débit est suffisant, un écoulement bi-phasique. Ceci est l’origine de la cavitation hydrodynamique. Les échanges de chaleur latente générés par l’apparition et l’implosion des bulles impliquent une variabilité importante de la température dans les zones au-delà de la constriction. En ajoutant des sondes de température nanométriques dans le fluide et en utilisant un microscope confocal on peut déterminer la température en un point. Ainsi on a pu établir des cartographies thermiques en 2 et 3 dimensions à l’intérieur d’un écoulement stationnaire bi-phasique. La technique permet en outre d’avoir accès à la quantité de gaz ce qui permet de corréler les gradients de température avec les zones de transitions de phases. Des zones de très forts refroidissements sont observées après la constriction, là où les bulles apparaissent. Par contre on n’observe pas les zones d’échauffement attendu à cause de la condensation. Une méthode complémentaire, moins sensible, utilisant la spectroscopie Raman a aussi été utilisée pour confirmer ce résultat. Par ailleurs une nouvelle classe de matériaux luminescents sensible à la température et la pression a été étudiée. Enfin une étude de la production de radicaux lors de l’implosion des bulles a été menée en utilisant la chimiluminescence du luminol. La technique utilisée par comptage de photons a permis de quantifier cette production et une cartographie de l’émission du luminol a permis d’associer celle-ci avec la zone d’implosion des bulles / A constriction in the microchannel can be used to establish a two-phase flow, when a sufficient liquid flux is introduced. This is known as hydrodynamic cavitation. The latent heat resulting from the growing and collapsing vapor bubbles makes it interesting to observe the temperature conditions in the flow downstream of the constriction. Using fluorescence microscopy, with the addition of temperature sensitive nano probes into the working fluid, we can determine the temperature at a single point, averaged over the integration time. Coupled with a confocal microscope, we were able to produce two and three dimensional temperature maps of the steady state flow in the microchannel by the use of ratiometric intensity measurements. This technic allows us to observe temperature gradients in two-phase flow as well yielding the void fraction information. Areas of substantial cooling are observed downstream the constriction in the two-phase flow, linked to the bubble growth, while heating regions due to condensations are missing. A complementary, yet less sensitive probe-less technique using the inherent Raman scattering signal of the liquid, was used to confirm the findings. A separate study evaluating a new group of luminescent materials for optical temperature and pressure probes is performed and discussed herein. Finally, the luminol chemiluminescent reaction with radicals produced by the cavitating flow, is used to obtain a corresponding photon yield. By counting the photons produced, an estimate on the radical yield can be obtained. Additionally, rudimentary mapping of the chemiluminescence signal allows the localization of the bubble collapse regions

Les Moteurs à Propergol Solide sont le siège d'instabilités de combustion qui se manifestent par des fluctuations de la pression interne et, en conséquence, de la poussée induite. Les phénomènes oscillatoires observés résultent d'un couplage entre l'acoustique de la chambre de combustion et une ou plusieurs instabilités hydrodynamiques de l'écoulement interne.Dans cette thèse, on s'intéresse à l'étude de l'écoulement induit par injection pariétale en géométrie axisymétrique. Plus particulièrement, on se concentre sur l'interaction qui peut avoir lieu entre cet écoulement et une structure : les Protections Thermiques de Face.Une étude expérimentale reposant sur le montage gaz froid VALDO de l'ONERA, modélisation représentative des MPS P230, a permis d'examiner l'influence d'obstacles, rigides et en élastomère, sur les instationnarités de l'écoulement. Le caractère instable de l'écoulement induit par injection pariétale est alors confirmé et le rôle primordial joué par les PTF sur les fréquences émergeant au sein de l'écoulement a été identifié.L'analyse de stabilité linéaire de l'écoulement, conduite en approche globale, permet l'étude de sa dynamique en tant qu'amplificateur de bruit. Les m¶mécanismes mis en jeu dans l'instabilité de l'écoulement induit par injection pariétale ainsi que son interaction complexe avec la couche de cisaillement issue de l'obstacle sont analysés. La réponse de l'écoulement à un forçage harmonique est alors examinée à la lumière des résultats expérimentaux.Cette étude a été complétée par une approche numérique de l'Interaction Fluide-Structure sur une configuration jugée critique. Une étape préliminaire à la simulation multi-physique est l'élaboration d'un modèle réaliste du comportement des PTF en élastomère. L'identification et la prise en considération des mécanismes propres aux élastomères dans la loi de comportement est ainsi détaillée. Les résultats issus des simulations sont confrontés aux résultats expérimentaux et rendent possible la validation d'un scénario de couplage dédié à l'Interaction Fluide-Structure au sein des MPS. / Fluid Structure Interaction of an inhibitor with the internal flow induced by wall injection was studied in an axisymmetric cold flow apparatus. Experiments were carried out over a wide range of injection velocities in order to underline how the obstacle not only modifies the mean flowfield but also its entire dynamic behaviour.The resulting instability (from the interaction of the unstable shear layer with the Taylor-Culick flow) exhibits a significant shift with respect to the Taylor-Culick instability and therefore emphasizes the strong impact of the inhibitor on hydrodynamics.The mecanisms responsible of such behaviour are studied in this thesis with the help of global linear stability analysis and multi-physics numerical computations.

De nombreux projets visant à réduire l’impact environnemental global des avions sont lancés au niveau européen. L’un des moteurs étudié pour les avions moyens et longs courriers est le moteur Ultra High Bypass Ratio (UHBR) : un moteur simple corps, double flux, à flux externe fortement augmenté. Le moteur UHBR doit être équipé d’un réducteur épicycloïdal, qui est un composant jamais utilisé dans un turboréacteur jusqu’à présent. L’optimisation d’un réducteur épicycloïdal a conduit à l’utilisation de paliers hydrodynamiques pour supporter les pignons satellites du porte-satellites.Pour une telle application, le palier hydrodynamique subit une déformation très élevée due aux charges de l’engrènement sur le pignon satellite et à l’effet centrifuge engendré par la rotation du porte-satellites. La géométrie optimisée des composants du palier varie avec le comportement thermique et mécanique des pièces, nécessitant la prise en compte d’une modélisation thermoélastohydrodynamique (TEHD).Afin de modéliser précisément ces phénomènes, un modèle conservatif dans la zone inactive en régime TEHD a été développé et validé à l’aide des résultats d’essais de la littérature et des bancs d’essais de Safran Transmission Systems. Les résultats obtenus montrent un double champ de pression dans le palier hydrodynamique du réducteur épicycloïdal, engendré par la forte déformation élastique du pignon agissant ainsi de façon significative sur le comportement dynamique du palier. De plus, l’influence de l’effet centrifuge sur l’huile dans le palier a été également examinée. / Several projects aimed at reducing the overall environmental impact of aircrafts are launched at European level. One of the engines studied for medium and long-haul aircraft is the Ultra High Bypass Ratio (UHBR) engine: a single-body, dual-flow gas turbine, with a greatly increased external flow rate. The UHBR engine must be equipped with an epicyclic reduction gearbox, which is a component never used before in a turbofan engine. The optimization of an epicyclic gearbox has resulted to the use of hydrodynamic bearings for supporting the satellite gears on the planet carrier.The hydrodynamic bearing of such application undergoes high deformations due to the gear contact loads on the satellite gear and to the centrifugal effect generated by the rotation of the planet carrier. The optimized geometry of the bearing components is influenced by the thermal and mechanical behavior of the bearing components, which requires a thermoelastohydrodynamic (TEHD) modeling taking into account the real film thickness under operation.In order to precisely simulate these phenomena, a conservative modeling in the inactive zone, under a TEHD regime has been developed and validated by comparing the predictions to the test results of the literature and of the test benches of Safran Transmission Systems. The results obtained show a double pressure field in the hydrodynamic bearing of the epicyclic reduction gearbox, due to the strongly deformed film thickness, which greatly affect the dynamic behavior of the bearing. In addition, the influence of the centrifugal force on the oil pressure in the oil film was also examined.

L’engouement pour l’originalité des propriétés physiques des nanoparticules s’est accompagné d’un développement de nombreuses méthodes de synthèse depuis un demi siècle. Parmi elles, l’ablation laser en liquide permet de produire des nanoparticules avec des surfaces libres de tout contaminant et ce pour une multitude de combinaisons de matériaux et de solvants. Cependant, la simplicité apparente de cette technique dissimule la complexité des mécanismes physico-chimiques, ce qui entraîne actuellement un manque de contrôle des objets synthétisés. Tout d’abord, nous nous sommes intéressés à la cinétique des bulles pour laquelle les conditionsextrêmes d’ablation laser en liquide présentent des cas originaux de cinétique dans le domaine de la mécanique des fluides. Puis, ce travail de thèse vise à donner de plus amples perspectives quant à une meilleure maîtrise de la structure cristalline des nanoparticules et de la stabilité colloïdale. Une manière plus directe de contrôler la taille, la phase cristalline et la stabilité colloïdale des solutions contenant des nanoparticules est d’ajouter des ligands. Nous avons donc étudié les mécanismes de stabilisation de ces solutions en utilisant des ions qui se complexent aux nanoparticules d’or. Nous avons aussi réussi à synthétiser des nanoparticules de rubis (alumine dopée chrome). La stabilisation de ces nanoparticules dans une phase métastable en utilisant des ligands organiques a été expliquée par une étude théorique / Laser generation of nanoparticles in liquids : new insights on crystal structure control and colloidal stability The great interest of nanoparticles for their original physical and an chemical properties has been supported by the development of numerous methods of synthesis. In the nineties, laser generation of nanoparticles in liquids appeared, including Pulsed Laser Ablation in Liquids (PLAL). The PLAL technique enables to produce surface free particles for plenty of material and solvent combinations. However, the apparent simplicity of its implementation hides complex physico-chemical mechanisms resulting in a lack of control of the final products. We firstly investigated the dynamics of the laser-generated bubbles for which the PLAL extreme conditions present new studied cases of bubbles dynamics not encountered in the field of fluid mechanics. Then, we aim to bring new insights into better control of the nanoparticles morphology and their colloidal stability. A straight way to tune sizes, crystal structures and the colloidal stability consists in the addition of stabilizing agents. Hence, we investigated the mechanisms of stabilization of colloidal gold using complexing ions. We also succeed to synthesis nano-rubies, i.e. chromium doped corundum alumina nanoparticles, unexpected at nanoscale. The stabilization of the metastable crystal structure using ligands is explained thanks to a comprehensive theoretical approach

Ce travail de thèse met en évidence deux phénomènes microscopiques dissipatifs au voisinage de la ligne de contact dans les phénomènes de capillarité et d'adhésion. L'étude s'appuie sur des expériences dynamiques et des modélisations théoriques. La mesure expérimentale de quantités macroscopiques et la théorie hydrodynamique permettent d'extraire les informations dynamiques localisées au voisinage la ligne de contact. Les différents phénomènes dissipatifs localisés au voisinage de la ligne de contact ont pour origines les propriétés des substrats sur lesquels se déplace la ligne de contact. Pour une surface rigide hétérogène, nous avons développé un modèle rhéologique de la ligne de contact fondé sur l'hydrodynamique, permettant d'établir théoriquement l'évolution temporelle de la ligne de contact et de ses déformations. Une décomposition modale fondée sur la réduction de l'énergie par la théorie du chemin de réaction fournit une prédiction quantitative de la dynamique thermiquement activée de la ligne de contact, en accord avec l'expérience réalisée. Pour un substrat déformable, à partir de l'analyse de deux expériences différentes de mouillage dynamique et de l'estimation de la dissipation dans le substrat, fondée sur sa viscoélasticité, nous avons développé une compréhension générale du comportement dynamique d'une ligne de contact sur un substrat viscoélastique. Pour finir, ce modèle de dissipation viscoélastique est appliqué au cas de l'adhésion réversible, où expérimentalement nous mesurons la dynamique de pelage et de recollement sur un substrat viscoélastique. Cette extension à l'adhésion permet de relier les phénomènes interfaciaux en une compréhension générale / This thesis reveals two dissipative microscopic phenomena close to the contact line in the fields of capillarity and adhesion. The study is based on dynamic experiments and theoretical predictions. Experimental measurement of macroscopic quantities and the hydrodynamic theory give access to dynamic information located close to the contact line. The different dissipative phenomena, located close to the contact line, originate from the properties of the substrates on which the contact line moves. For a heterogeneous rigid surface, we have developed a rheological model of the contact line based on hydrodynamics, in order to theoretically establish the temporal evolution of the contact line and its deformations. A modal decomposition based on the reduction of the energy by the reaction path theory allows a quantitative prediction of the thermally activated dynamics of the contact line, in agreements with the experiment carried out. For a deformable substrate, based on the analysis of two different experiments of wetting dynamics and on the estimation of the dissipation in the substrate founded on its viscoelasticity, a general understanding of the dynamical behavior of contact lines on viscoelastic substrates is achieve. Finally, this viscoelastic model is applied to the case of weak adhesion, where experimentally we measured the peeling dynamics from a viscoelastic substrate. This extension to adhesion bridge the gap between different interfacial phenomena into a general understanding

La dynamique des écoulements géophysiques planétaires est fortement influencée par des processus physiques souvent non résolus par les modèles numériques de circulation générale. Il est essentiel de comprendre les mécanismes physiques sous-jacents pour paramétrer l’effet des petites échelles sur les grandes. Dans cette thèse, nous étudions un problème emblématique d’interactions entre ondes et écoulements moyens : la dynamique des écoulements zonaux forces par des ondes internes de gravite. Une manifestation remarquable de ces interactions est l’oscillation quasi-biennale (QBO) des vents équatoriaux dans l’atmosphère terrestre. Dans un premier temps, nous décrivons une transition vers le chaos dans un modèle quasi-linéaire classique du QBO. Nous montrons que ces bifurcations persistent dans des simulations numériques directes. A l’aune de ces résultats, nous proposons une interprétation de l’observation d’une rupture inattendue de la périodicité du QBO en 2016. Le mécanisme de génération d’écoulements moyens par les ondes dans les fluides stratifies nécessite la prise en compte d’effets dissipatifs. Il s’agit d’un phénomène analogue au "streaming" acoustique. Dans un second temps nous exploitons cette analogie en étudiant la génération d’écoulements moyens par les ondes internes proche d’une paroi, a l’aide d’approches asymptotiques multi échelles. Finalement, nous proposons une approche inertielle pour décrire l’émergence spontanée d’écoulements vorticaux en présence d’ondes : nous appliquons les outils de mécanique statistique pour calculer la partition d’énergie entre petites et grandes échelles dans le modèle d’eau peu profonde. / The dynamics of planetary-scale geophysical flows is strongly influenced by physicalprocesses, mostly unresolved by general circulation numerical models. To parametrisethe coupling between small and large scales, it is essential to understand the underlying physical mechanisms. In this thesis, we study an emblematic problem of interactions between waves and mean flows: the dynamics of zonal flows forced by internal gravity waves. A striking manifestation of these interactions is the quasi-biennial oscillation (QBO) of equatorial winds in the Earth’s atmosphere. First, we describe a transition to chaos in a classical quasilinear model of the QBO and show that these bifurcations persist in direct numerical simulations. Based on these results, we suggest an interpretation for the observation of the unexpected periodicity disruption of the QBO in 2016. The mechanism by which mean flows are generated by waves in stratified fluids requires the consideration of dissipative effects. This phenomenon is analogous to acoustic "streaming". In a second time, we exploit this analogy to study the generation of mean flows by internal gravity waves close to a wall, using multi-scale asymptotic approaches. Finally, we propose an inertial approach to describe the spontaneous emergence of vortical flows in the presence of waves: we apply the tools of statistical mechanics to calculate the partition of energy between small and large scales in the shallow-water model.

L'environnement côtier de l'île de Mayotte (nord du canal du Mozambique, océan Indien) se caractérise par l'existence d'un vaste complexe récifo-lagonaire au sein duquel se développe un littoral particulièrement diversifié et marqué par la présence de nombreuses mangroves. Plusieurs campagnes de mesures hydrodynamiques et d'observations morphologiques ont été effectuées afin d'apporter des éléments de compréhension sur l'évolution et le fonctionnement morphodynamique de ce type de littoral mixte à mangroves et récifs coralliens. Les observations de terrain conjuguées à une analyse diachronique de photographies aériennes ont permis de montrer le caractère très variable de la dynamique des mangroves avec une progression ou une stabilité au nord et à l'est de l'île, et inversement, une évolution nettement régressive au sud et à l'ouest. Les mesures hydrodynamiques ont donné lieu à une analyse des paramètres de houle et de courant en se focalisant sur leurs variabilités spatiales et temporelles. Les résultats obtenus montrent l'importance de l'architecture récifale sur les processus hydrosédimentaires et le développement des formes littorales. L'évolution et la dynamique des mangroves de Mayotte apparaissent ainsi fortement influencées par les caractéristiques morphologiques des structures récifales amont qui engendrent une plus forte vulnérabilité des mangroves du sud et de l'ouest face à une déstabilisation naturelle ou anthropique. La disparition des formations végétales des mangroves du sud et de l'ouest de l'île entraine une réorganisation des répartitions sédimentaires et la mise en place de nouveaux systèmes de fonctionnement morphodynamique. L'ampleur des phénomènes constatés a montré la nécessité de mettre en place un observatoire opérationnel de surveillance de la dynamique côtière. Cet outil, basé sur un réseau de mesures et d'observations de terrain et sur la réalisation d'un Système d'Information Géographique (SIG), a pour but de favoriser le suivi et l'évaluation de l'état du littoral. Dans un contexte local de fort développement et dans celui, plus global, de mutations environnementales, cet outil de surveillance et d'analyse permet une meilleure compréhension de la morphodynamique des côtes mahoraises et apporte une aide à la décision et à la gouvernance afin d'orienter les politiques de gestion de ce milieu fragile, à forts enjeux et en constante évolution.

Cette étude définit le fonctionnement du système hydrothermal de la Léchère, grâce à une approche interdisciplinaire associant (1) géologie, (2) hydrodynamique, (3) hydrogéochimie et (4) modélisation hydrodynamique/thermique.<br />(1) L'accident tectonique, marquant le contact entre la couverture carbonatée et les micaschistes du socle, est jalonné par des écailles discontinues d'évaporites triasiques. A l'intersection entre cet accident et la vallée de l'Isère, les eaux thermales se diffusent dans l'ombilic glaciaire par l'intermédiaire d'un écroulement ancien, identifié par la réinterprétation des descriptions de forages.<br />(2) Le système thermal est caractérisé par la présence d'une zone décomprimée dans laquelle ont lieu les circulations superficielles et qui sert également de trop-plein au système profond. L'ensemble est alimenté de manière diffuse par l'infiltration des eaux de pluie et de manière ponctuelle par deux pertes sur le réseau hydrographique. Dans la zone d'émergence, les circulations resteront ascendantes tant que le niveau de l'aquifère thermal restera supérieur à celui de l'aquifère superficiel.<br />(3) Au début de leur parcours, les eaux thermales interagissent avec les dolomies, cargneules et gypses (Ca, SO4, Mg), au même titre que les eaux triasiques superficielles. Elles se différencient ensuite par un parcours profond de plusieurs milliers d'années, leur permettant d'interagir avec la halite (Cl), les roches métamorphiques (Si, K, échange Ca/Na) et les roches permo-carbonifères (matière organique, H2S), avant<br />de se mélanger avec les eaux superficielles de la zone d'émergence.<br />(4) La modélisation numérique couplée hydrodynamique / thermique a mis en évidence la chasse des eaux chaudes vers la vallée lors de la réactivation post-glaciaire des écoulements. La mise en exploitation de la ressource thermale accélère, toutefois de façon imperceptible à l'échelle d'observation, le refroidissement naturel de l'aquifère thermal.

Les travaux du PIREN-Seine, dans lesquels s'inscrit mon travail de thèse, se sont dès le début orientés vers la recherche d'une description du fonctionnement écologique du système aussi déterministe que possible. Dans ce but, l'approche adoptée consiste à décrire les organismes et constituants ainsi que les processus qui ont une influence sur la qualité de l'eau dans le réseau hydrographique. La modélisation a été l'un des outils utilisés pour améliorer notre connaissance des processus et ainsi aider les différents acteurs de l'eau à la prise de décision. Pour étudier l'impact des rejets urbains dans les zones fortement influencées par l'activité humaine, il a été décidé de mettre en oeuvre un modèle bi-dimensionnel plan capable de simuler de façon précise le devenir des différents polluants rejetés en Seine. Grâce à la collaboration entre l'École des Mines et la Direction Études et Recherches d'Électricité de France, nous avons pu utiliser le code de calcul hydrodynamique aux éléments finis, TELEMAC 2D, d'E.D.F. Une première étape a consisté à mettre en oeuvre ce modèle dans le bief de Seine caractérisé par les rejets de la station d'épuration d'Achères pour estimer les vitesses de l'écoulement et les hauteurs d'eau à l'étiage et de connaître, entre autres, le temps de résidence de l'eau pour quelques situations permanentes de bas débit. Les travaux antérieurs du PIREN-Seine ont largement montré que la pollution d'origine organique était étroitement liée aux particules. Bien que l'essentiel du flux de matériel particulaire soit transporté au cours des crues, l'étude du transport en période d'étiage ne doit pas être négligée. La dilution des rejets urbains par l'eau de Seine est plus faible en période d'étiage. Le dépôt est plus intense, le matériel qui le constitue est fortement pollué et les zones de stockage sont plus nombreuses. Il nous est donc apparu important de simuler convenablement le transport solide pour la Seine à l'étiage. Nous avons alors modifié le modèle de transport bi-dimensionnel d'EDF, SUBIEF 2D. La navigation, intense dans cette partie du fleuve, joue un rôle considérable sur les temps de résidence des particules et ainsi sur le devenir des produits contaminants peu ou non solubles. Nous avons donc développé un module de remise en suspension qui prend en compte de façon précise la circulation fluviale : la modélisation du courant de retour, du champ de vitesses généré sur le lit de la rivière par les jets d'hélice, l'utilisation d'une formulation énergétique des processus de dépôt/érosion basée sur la capacité de transport d'un écoulement, nous permet de reproduire les concentrations en MES mesurées dans la colonne d'eau sur l'ensemble du domaine lors des différentes campagnes organisées. Nous en tirons des informations sur la quantité de particules sédimentées dans le bief, la localisation des dépôts frais et le temps de transfert moyen des particules en périodes de basses eaux. Nous appuyant sur le modèle hydraulique et de transport précédemment défini, nous avons mis en place un module biologique reprenant le schéma conceptuel du modèle RIVE utilisé dans le PIREN-Seine. Nous simulons le cycle de l'oxygène en relation avec celui du carbone. Les processus modélisés sont le développement phytoplanctonique, la dégradation de la matière organique par les bactéries hétérotrophes et la production/consommation d'oxygène associée. Une ré-aération à l'interface air-eau est aussi prise en compte. Le modèle simule les évolutions de 24 variables d'état relatives à la colonne d'eau et à la vase non consolidée déposée sur le lit du fleuve. Les processus d'échange entre la colonne d'eau et les sédiments, qu'ils soient diffusifs ou contraints par le dépôt et l'érosion, ont été introduits explicitement dans le modèle. La couche de vase est également le siège de réactions biologiques. Les résultats des simulations pour différents étiages sont confrontés aux données expérimentales dont nous disposons. Ces résultats montrent que la conceptualisation retenue pour le modèle de simulation de la qualité de l'eau est adaptée à la simulation des processus qui dominent dans ce secteur. Une comparaison entre le modèle ProSe, développé à l'École des Mines de Paris, dans ses versions mono-dimensionnelles et à tube de courant, et le modèle bi-dimensionnel TELEMAC-SUBIEF a été réalisée. Cette intercomparaison permet, d'une part, de valider les développements réalisés au sein de SUBIEF 2D, le modèle ProSe ayant été validé dans de nombreuses situations et, d'autre part, de montrer les différences de comportement des modèles mono et bi-dimensionnels dans le secteur d'étude.