Cette thèse porte sur l'étude et la modélisation de l'orientation de fibres de verre dans des thermoplastiques renforcés. Dans une optique de maîtrise et de compréhension du comportement des composites renforcés de fibres longues, nous proposons dans ce travail d'étudier les phénomènes d'orientation et la rhéologie de ces suspensions. Une approche expérimentale a permis de caractériser l'orientation et la structure des fibres dans des pièces injectées. L'analyse de la distribution de longueur des fibres dans une pièce a mis en évidence que la casse des fibres se faisait préférentiellement durant la phase de la plastification de la matière. Les mesures d'orientation, réalisées avec des techniques de mesure 2D et 3D montrent que la concentration en fibres est un paramètre déterminant sur la formation et les caractéristiques de la structure coeur-peau. Une approche numérique de modélisation de l'orientation des fibres dans des pièces injectées, basée sur le calcul direct du mouvement d'orientation d'une population de fibres en écoulement, a permis d'accéder à la physique intrinsèque de la suspension. Ce modèle repose sur une formulation éléments finis multi- domaine développée dans le cadre du logiciel Rem3D. La principale originalité du modèle est qu'il n'est pas nécessaire d'exprimer de manière explicite l'ensemble des forces et interactions hydrodynamiques régissant le système. Lorsqu'on applique un cisaillement simple à un volume élémentaire représentatif, le calcul de simulation directe nous permet d'étudier les phénomènes d'interaction entre particules en fonction de la concentration et du rapport de forme des fibres. Ainsi, nous proposons une méthode d'identification du coefficient d'interaction Ci , paramètre déterminant du modèle statistique de Folgar & Tucker. Les valeurs de Ci sont fonction du rapport de forme et de la concentration des particules. On s'intéresse d'autre part à la validité des approximations de fermeture. On montre que l'approximation hybride est satisfaisante dans le cas d'un écoulement de cisaillement simple. Enfin, ce modèle est appliqué à l'étude rhéologique d'une suspension de fibres.

Le développement de nouvelles technologies pour le transport aérien moins polluant est de plus en plus basé sur la simulation numérique, qui nécessite alors une description fiable de la chimie. Pour la plupart des carburants, la description de la combustion nécessite des mécanismes détaillés mais leur utilisation dans une simulation numérique de combustion turbulente est limitée par le coût calcul. Des mécanismes cinétiques réduits et des méthodes de tabulation ont été proposés pour surmonter ce problème. Ces descriptions chimiques simplifiées ayant été développées dans le cadre de configurations laminaires, cette thèse propose de les évaluer dans des configurations turbulentes: une DNS de flamme prémélangée méthane/air de type Bunsen et une LES d'un brûleur expérimental. Les mécanismes sont analysés en termes de structure de flamme, paramètres de flamme globaux, longuer de flamme, prediction des concentrations en espèces majoritaires et des émissions polluantes. Une méthodologie pour évaluer a priori la capacité d'un mécanisme à prédire correctement des phénomènes chimiques tridimensionnels est proposée en se basant sur les résultats de flammes laminaires monodimensionnelles non étirées et étirées. Il ressort que, d'une part, pour construire un mécanisme réduit, il est nécessaire de faire un compromis entre coût calcul, robustesse et qualité des résultats. D'autre part, la qualité des résultats de DNS et LES de configurations tridimensionnelles turbulentes peut être anticipée par une analyse du comportement des schémas réduits dans des configurations simplifiées de flammes monodimensionnelles laminaires non étirées et étirées.

Les jets annulaires turbulents sont des écoulements soumis à de fortes instabilités qui peuvent devenir un handicap pour certaines applications industrielles (procédés de fibrage du verre, ventilation individuelle, ventilation opérationnelle, brûleurs, neige de culture ...). Ces écoulements particuliers de jets axisymétriques sont complexes puisqu'ils regroupent des phénomènes de sillage et de mélange. L'objectif de cette étude est de déterminer des méthodes de contrôle de la zone initiale des jets annulaires d'air de très grand rapport de diamètres (r=0,91), pour des nombres de Reynolds élevés (Ree=5130, Ree=2565, Ree=1368, e étant l'épaisseur de l'anneau de jet de 2,565mm), afin de répondre aux attentes des industriels mais aussi de mieux comprendre les mécanismes mis en jeu dans la formation de l'écoulement. Deux méthodes ont été étudiées. Une voie de contrôle passif par modification de géométrie de l'obstacle central a permis d'agir sur le sillage et de modifier la morphologie du jet. Les jets annulaires basique, conique et sphéroïdal ont été comparés pour obtenir des jets plus efficaces pour certaines applications. Un contrôle actif par forçage acoustique a été appliqué au jet annulaire basique afin de réduire le battement du jet autour de son axe sans modifier la géométrie de la buse. Différents paramètres de forçage (fréquence, position de l'excitation, phase de l'onde) ont été observés pour déterminer une configuration efficace de contrôle des instabilités de l'écoulement. Les méthodes de diagnostic utilisées sont intrusives (anémométrie par fil chaud) ou non intrusives (Tomographie rapide, Vélocimétrie par Images de Particules, Anémométrie Doppler Laser 2D et 3D) et sont complétées par un post-traitement de Décomposition Orthogonales en Modes Propres (P.O.D.) pour analyser les structures cohérentes de l'écoulement et déterminer l'influence du contrôle passif ou actif sur ces structures de différentes échelles. Une Simulation Numérique Directe (D.N.S.) est enfin proposée pour modéliser le contrôle passif des écoulements de jets annulaires turbulents de très grands rapports de diamètres.

La turbulence modifie sensiblement les propriétés des écoulements évoluant en zone de proche paroi en provoquant notamment une augmentation significative du frottement pariétal. Dans le but de réduire cette traînée turbulente, on propose dans cette étude un contrôle actif de type dual dont la principale caractéristique repose sur l'interaction de deux schémas de contrôle, un schéma en boucle ouverte par soufflage localisé en espace et périodique en temps, suivi d'un schéma en boucle fermée de type suboptimal. L'idée principale est alors d'organiser la turbulence par une excitation oscillante de manière à la rendre plus prédictible et de ce fait plus contrôlable.<br />On commencera par étudier séparément l'impact de ces deux types de contrôle sur la turbulence de proche paroi. Une étude expérimentale est menée afin d'investiguer en détail l'effet de l'instationnarité induite par une excitation périodique. Les méthodes suboptimales ne pouvant être à ce jour appliquées expérimentalement, une Simulation Numérique Directe s'avère être l'outil le mieux adapté pour étudier leurs comportements. On montre finalement qu'une stratégie de contrôle duale, c'est à dire combinant les deux manipulations (Périodique + Suboptimale), peut s'avérer être une méthode de contrôle efficace afin de réduire la contrainte turbulente.

L'approche de simulation numérique directe est utilisée pour la simulation d'un écoulement en canal pleinement turbulent afin d'étudier l'influence des grandes échelles de l'écoulement ainsi que la dynamique du transport des contraintes de Reynolds et de la vorticité. Les simulations sont réalisées sur un domaine de calcul de grande taille afin de pouvoir capturer l'intégralité des grandes structures de l'écoulement, et portent sur une gamme relativement étendue de nombres de Reynolds (Reτ =180, 395, 590 et 1100) allant des écoulements faiblement turbulents à des écoulements modérément turbulents. L'invariance remarquable des fluctuations de vorticité normale est expliquée à travers une analyse spectrale de la vorticité. L'étude des différents termes du transport de l'intensité turbulente de la vorticité révèle par ailleurs que le pic de production de la vorticité transverse est situé à proximité immédiate de la paroi et pourrait ouvrir la voie à des stratégies de réduction de la traînée basées sur la réduction de la vorticité transverse. Le transport des contraintes de Reynolds dans la couche interne et dans la couche de recouvrement est également étudié. A proximité des parois, la dépendance des termes de transport avec le nombre de Reynolds dans les échelles internes montre que ces dernières ne suffisent pas à caractériser la dynamique des contraintes de Reynolds dans cette zone. Cette insuffisance des échelles internes nous a amenés à nous intéresser plus particulièrement au processus de production à travers les statistiques de la production conditionnées par le passage par niveau des fluctuations de la vitesse normale ou longitudinale. Cette étude nous a permis d'identifier les fluctuations qui contribuent le plus à la production et celles qui sont à l'origine de la dépendance avec le nombre de Reynolds. / The direct numerical simulations of a fully turbulent channel flow are investigated to study the large scales effects on the flow quantities such as the Reynolds stresses and vorticity transport processes. Large computational domains are used so as to cover the largest scales of the flow. The simulations are performed in a wide range of Reynolds numbers (Reτ=180, 395, 590 and 1100) going from weakly to moderately high Reynolds number turbulent flows. The invariance of the wall-normal vorticity fluctuations scaled in wall variables in the inner layer versus the Reynolds number is analyzed using a spectral analysis. The vorticity transport equations are investigated in detail, presumably for the first time. The transport mechanism of the Reynolds shear stresses are subsequently analyzed in the inner layer and the overlapping zone. In the wall layer, different terms of the Reynolds stresses transport expressed in inner scales depend on the Reynolds number. This scaling failure lead us to focus on the statistics of the production when the streamwise or normal velocity fluctuations cross a given level, through the conditional Palm statistics. The main aim is to identify those amplitudes of the fluctuations that contribute more to the production and those which are responsible for the production Reynolds dependence.

Les travaux menés, depuis plusieurs années, au CORIA ont abouti à la construction d’un outil numérique (ARCHER) permettant la simulation numérique directe d’écoulements diphasiques et notamment l’atomisation d’un jet liquide à haute vitesse. Ce type de simulation permet de capturer les phénomènes d’atomisation au voisinage de l’injecteur difficilement caractérisables par les outils expérimentaux actuels. Ces simulations requièrent des conditions d’injection délicates à évaluer a priori car elles dépendent des caractéristiques de l’écoulement au sein de l’injecteur. Or, certains jets présentent une grande sensibilité à ces conditions d’injection. Dès lors, il est nécessaire de simuler l’écoulement au sein de l’injecteur afin d’appréhender la nature de cette sensibilité. L’utilisation d’un maillage cartésien par le code ARCHER conjuguée à la volonté de simuler le système d’atomisation dans son ensemble ont orienté ces travaux vers l’utilisation d’une méthode de frontière immergée. Ces travaux ont ainsi permis de reproduire des écoulements au sein d’injecteurs de forme quelconque tout en conservant le maillage cartésien d’origine, précieux tant pour l’efficacité du solveur que pour sa précision. Dans un premier temps, l’implantation dans le code ARCHER d’une méthode de frontière immergée a été réalisée et testée sur des configurations de canal et de conduite et de l’écoulement autour d’un cylindre. L’application de cette méthode a porté sur la simulation de l’écoulement au sein d’un injecteur triple disque mono-trou et a notamment permis de caractériser l’origine de l’écoulement secondaire formé dans l’orifice de décharge. Afin d’évoluer vers la construction d’un outil numérique capable de simuler le système d’atomisation dans son ensemble, un couplage entre la méthode de frontière immergée et la méthode Ghost fluid a été nécessaire. La version bi-dimensionnelle développée a été testée sur la relaxation d’une goutte posée sur une paroi. Cette version a permis de simuler des écoulements au sein de canaux à différents rapports de longueur sur diamètre et l’écoulement au sein d’une buse convergente. La simulation simultanée de l’écoulement interne et externe a permis de lier les fluctuations de vitesses des écoulements internes à la création de surface engendrée sur les écoulements externes. / Since several years, the research conducted at the CORIA laboratory led to the development of a numerical tool (ARCHER) alllowing direct numerical simulations of two phase flows. In particular, the simulations of high speed liquid jet primary break-up have been strongly investigated. These simulations are able to capture primary break-up phenomena near the nozzle exit where experimental characterisations are difficult to conduct. These simulations need injection conditions tricky to gauge a priori, since they depend on the flow characteristics inside the nozzle. Moreover, some jets are highly sensitive to these injection conditions. Therefore, it becomes necessary to simulate the flow inside the nozzle to better understand this sensitive nature. The objective to simulate the whole atomization system guided the present work dedicated to the use of an immersed boundary method (IBM). Such an approach allows reproducing flows inside nozzles of arbitrary shape while keeping the original cartesian mesh valuable for numerical efficiency and accuracy. As a first step, the implementation of an IBM in ARCHER was carried out and tested on channels, pipes and uniform flows past a circular cylinder. An industrial application focused on the flow inside a triple disk compound injector. This work led to a refined description of the secondary flow origin in the discharge hole. In order to move towards the design of a numerical tool able to simulate the whole injection system, a coupling between IBM and the Ghost Fluid Method (GFM) has been found necessary. This allows accounting for two phase flows inside the nozzle where the dynamics of the triple line has to be considered. The bidimensional developments have been tested on drops released on walls. This version enabled to simulate flows inside channels with different ratios of length over diameter and the flow inside a convergent nozzle. The simultaneous computation of flows inside and outside nozzle has enabled to link the velocity fluctuations of internals flows to the surface setting-up gene-rated on external flows.

La stabilisation et la propagation d’une flamme laminaire pré-mélangée méthane/air dans un canal étroit, sont revisitées à partir de simulations numériques. La combustion est modélisée à l’aide d’une chimie et de propriétés de transport complexes, ainsi que du couplage des transferts thermiques à l’interface et dans les parois. Premièrement, une procédure de réduction des mécanismes chimiques adaptée à cette application est appliquée. Deuxièmement, la réponse de la forme de flamme lorsque soumise à diverses conditions thermiques est analysée en termes de vitesse de propagation et de topologie de l’écoulement au voisinage du front de réaction. Troisièmement, le mécanisme de transfert thermique déclenchant la propagation de flamme lorsque celle-ci est soumise à un préchauffage est montré être principalement convectif. Pour finir, le rôle prépondérant de la gravité, via l’action du moment barocline, sur la symétrie des flammes se propageant dans des canaux étroits, est démontré. / The flow physics controlling the stabilization and propagation of a methane/air laminar premixed flame in a narrow channel is revisited from numerical simulations. Combustion is described with complex chemistry and transport properties, along with a coupled simulation of heat transfer at and within the wall. First, a chemistry mechanism reduction procedure fitted to this application is applied. Second, the response of the premixed flame shape to various heat transfer conditions is analyzed in terms of flame propagation velocity and flow topology in the vicinity of the reactive front. Third, the heat transfer mechanism triggering the flame movement when this last is submitted to an upstream wall preheating is revealed to be mainly convective. To finish, the preponderant role of gravity, via an impact on the baroclinic torque, in the symmetry breaking of small-scale channel flames is demonstrated.

Dans ce travail de thèse, les mécanismes à l'origine de la formation des chocs associés à la perception de crackle proche de jets supersoniques axisymétriques sont étudiés à l'aide de simulations numériques. Dans ces simulations, les équations de Navier-Stokes instationnaires et compressibles sont résolues en coordonnées cylindriques à l'aide de différences finies d'ordre élevé peu dissipatives et peu dispersives. Quatre jets temporels à des nombres de Mach de 2 et~3 et à des nombres de Reynolds compris entre 3125 et 50000 sont simulés dans un premier temps. Des ondes acoustiques de forte amplitude présentant d'importants gradients de pression sont mises en évidence à proximité des jets. Elles se forment par un mécanisme de raidissement à la source qui est étudié par le calcul de moyennes conditionnelles synchronisées autour des pics de pression en champ proche. Ces moyennes montrent un lien direct entre ces ondes non linéaires et la convection de structures cohérentes à desvitesses supersoniques dans les couches de~mélange. L'influence de la température sur la formation de ces ondes est examinée dans un second temps par le calcul de cinq jets temporels à des rapports de température de 1, 2 et 4, et à des nombres de Mach acoustique compris entre 2 et 4. À vitesse d'éjection constante, les niveaux de bruit produits par les jets chauds sont moins élevés que ceux du jet isotherme, mais les ondes non linéaires qu'ils rayonnent sont peu affectées par une hausse de température. À nombre de Mach constant, les niveaux augmentent avec la température, de même que l'asymétrie des fluctuations de pression, traduisant un renforcement du caractère non linéaire des ondes rayonnées. Ces variations pourraient être dues à celles de la vitesse de convection des structures cohérentes, qui augmente de façon significative avec la température lorsque le nombre de Mach est constant, mais diminue légèrement à vitesse~constante. Finalement, trois simulations de jets spatiaux isothermes et chauds à un nombre de Mach acoustique de 2 et à des nombres de Reynolds de 12500 et 50000 sont mises en \oe uvre. Des ondes de Mach présentant d'importants gradients de pression sont visibles au voisinage direct des jets. La formation de ces ondes est liée, comme dans le cas des jets temporels, à la convection supersonique de structures cohérentes dans les couches de mélange. Le champ lointain acoustique est enfin déterminé par des méthodes d'extrapolation linéaire et non linéaire. Lorsque la propagation est non linéaire, un raidissement additionnel des fronts d'onde est constaté en champ lointain. / Numerical simulations are carried out with the aim of investigating the formation of nonlinear steepened waves at the origin of crackle in the near acoustic field of supersonic jets. In these simulations, the compressible Navier-Stokes equations are solved in cylindrical coordinates using high-order low-dissipative and low-dispersive finite difference schemes.Four temporally-developing isothermal round jets are first simulated at Mach numbers of~2 and~3 and at Reynolds numbers ranging from 3,125 to 50,000. Strong acoustic waves containing sharp pressure variations are observed in the vicinity of the jets. Their formation process is described by the computation of conditional averages which are triggered by the detection of strong pressure peaks in the near field. Such steepened waves are then shown to be produced by the supersonic motion of coherent structures inside the jet shear layers.Temperature effects are then investigated by considering five temporal round jets at temperature ratios of 1, 2 and~4 and at acoustic Mach numbers of 2, 2.8 and 4. For a given jet speed, the sound levels produced by the hot jets are lower than those of the isothermal one. However, the properties of the steepened waves they generate are not significantly affected by a rise of temperature. On the contrary, when the Mach number is held constant, pressure levels are higher at high temperature. The skewness and kurtosis factors of pressure fluctuations are also increased, which indicates a strengthening of the asymmetry and the intermittency of the pressure fluctuations. It is likely that the influence of temperature on these waves results from the variations of the convection speed, which is found to significantly increase with temperature at constant Mach number, but to slightly decrease at constant jet speed.Finally, three simulations of spatially-developing axisymmetric, isothermal and hot jets at a Mach number of~2 and at Reynolds numbers of 12500 and 50000 are performed. Strong Mach waves possessing the distinctive features of crackle are visible in the near vicinity of the jets. As observed for temporal simulations, their formations are associated with the supersonic motion of large-scale coherent structures inside the jet shear layers. The far acoustic field is determined using linear as well as nonlinear extrapolation methods. When nonlinear propagation effects are taken into account, a further steepening of the wavefronts is observed with increasing propagation distance.

Lors d'un hypothétique accident grave de réacteur à eau sous pression, un mélange de matériaux fondus, appelé corium, issu de la fusion du cœur peut se relocaliser dans le puits de cuve constitué par un radier en béton. Les codes d'évaluation réacteur pour simuler la phénoménologie de l'interaction corium-béton sont basés sur une description à grande échelle des échanges qui soulève de nombreuses questions, tant sur la prise en compte des phénomènes multi-échelles mis en jeu que sur la structure adoptée de la couche limite au voisinage du front d'ablation. Dans ce contexte, l'objectif principal de ce travail consiste à aborder le problème de la structure de la couche limite par simulation numérique directe. Ce travail s'inscrit dans le cadre plus général d'une description et d'une modélisation multi-échelle des échanges, c'est-à-dire de l'échelle locale associée au voisinage du front d'ablation jusqu'à l'échelle du code d'évaluation réacteur. Une telle description multi-échelle des échanges soulève le problème de la description locale de l'écoulement multiphasique multiconstituant mais aussi le problème du changement d'échelle et en particulier le passage de l'échelle locale à l'échelle de description supérieure dite macroscopique associée aux mouvements convectifs dans le bain de corium. Parmi les difficultés associées au changement d'échelle, nous nous intéressons à la problématique de la construction de conditions aux limites effectives ou lois de parois pour les modèles macroscopiques. Devant la complexité du problème multiphasique multiconstituant posé au voisinage du front, cette contribution a été abordée sur un problème modèle. Des conditions aux limites dites effectives ont été construites dans le cadre d'une méthode de décomposition de domaine puis testées pour un problème d'écoulement laminaire de convection naturelle sur parois rugueuses. Mˆeme si le problème traité reste encore éloigné des applications visées, cette contribution offre de nombreuses perspectives et constitue une première étape d'une modélisation multiéchelle des échanges pour la problématique de l'interaction corium-béton. Dans le cas plus complexe des écoulements multiphasiques multiconstituants et devant les difficultés expérimentales associées, le développement de lois de parois pour les outils existants aux échelles de description supérieures nécessite, au préalable, de disposer d'un outil de simulation numérique directe de l'écoulement au voisinage du front d'ablation. L'outil développé dans ce travail correspond à un modèle de Cahn-Hilliard/Navier-Stokes pour un mélange diphasique (liquide-gaz) compositionnel (corium-béton fondu) s'appuyant sur une description du système selon trois paramètres d'ordre associés respectivement aux fractions volumiques du gaz et aux deux espèces miscibles de la phase liquide ainsi que sur une décomposition de l'énergie libre selon une contribution diphasique et compositionnelle. Les équations de transport sont dérivées dans le cadre de la thermodynamique des processus irréversibles et résolues sur la base d'une application éléments finis de la plate-forme PELICANS. Plusieurs expériences numériques illustrent la validité et les potentialités d'application de cet outil sur des problèmes diphasiques et/ou compositionnels. Enfin, à partir de l'outil développé, nous abordons par simulation numérique directe une étude de la structure de la couche limite au voisinage du front d'ablation pour des bétons siliceux et silico-calcaire. / In the late phases of some scenario of hypothetical severe accident in Pressurized Water Reactors, a molten mixture of core and vessel structures, called corium, comes to interact with the concrete basemat. The safety numerical tools are lumped parameter codes. They are based on a large averaged description of heat and mass transfers which raises some uncertainties about the multi-scale description of the exchanges but also about the adopted boundary layer structure in the vicinity of the ablation front. In this context, the aim of this work is to tackle the problem of the boundary layer structure by means of direct numerical simulation. This work joins within the more general framework of a multi-scale description and a multi-scale modeling, namely from the local scale associated with the vicinity of the ablation front to the scale associated with the lumped parameter codes. Such a multi-scale description raises not only the problem of the local description of the multiphase multicomponent flow but also the problem of the upscaling between the local- and the macro-scale which is associated with the convective structures within the pool of corium. Here, we are particularly interested in the building of effective boundary conditions or wall laws for macro-scale models. The difficulty of the multiphase multicomponent problem at the local scale leads us to consider a relatively simplified problem. Effective boundary conditions are built in the frame of a domain decomposition method and numerical experiments are performed for a natural convection problem in a stamp shaped cavity to assess the validity of the proposed wall laws. Even if the treated problem is still far from the target applications, this contribution can be viewed as a first step of a multi-scale modeling of the exchanges for the molten core concrete issue. In the more complicated case of multiphase multicomponent flows, it is necessary to have a direct numerical simulation tool of the flow at the local scale to build wall laws for macro-scale models. Here, the developed tool corresponds to a Cahn-Hilliard/Navier-Stokes model for a two-phase compositional system. It relies on a description of the system by three volume fractions and on a free energy composed by a two-phase part and a compositional part. The governing equations are derived in the frame of the thermodynamic of irreversible processes. They are solved on the basis of a finite element application of the object-oriented software component library PELICANS. Several numerical experiments illustrate the validity and the potentialities of application of this tool on two-phase compositional problems. Finally, using the developed tool, we tackle by means of direct numerical simulation the problem boundary layer structure in the vicinity of the ablation front for limestone-sand and siliceous concretes.

Cette 'Habilitation à Diriger des Recherches' a pour objectif de retracer le parcours et les activités du candidat sur les dix dernières années. Elle donne un aperçu de l'expérience acquise dans les domaines de l'enseignement, de la recherche et de l'implication dans la vie de l'établissement. La première partie récapitule le parcours académique et les activités du candidat ainsi que sa production scientifique. La deuxième partie présente une description synthétique des activités de recherche portant sur le thème principal qui traite de la modélisation analytique des couches limites turbulentes et la dispersion de particules solides en suspension. L'idée principale de cette approche consiste à proposer à l'ingénieur des outils à la fois pratiques et précis. Afin d'améliorer la description des écoulements turbulents proche parois, nous avons commencé par travailler à la résolution analytique de l'équation de l'énergie cinétique turbulente k. Nous avons proposé une solution générale qui a permis d'obtenir une fonction au voisinage immédiat de la paroi. Cette fonction a été d'abord validée par des données de simulation numérique directe (DNS) puis évaluée au niveau du code Fluent. Ces résultats ont permis de développer des modèles analytiques de viscosité turbulente qui ont permis de reconstituer avec précision le champ des vitesses en canal turbulent. A partir d'une analyse des équations de Navier-Stokes moyennées au sens de Reynolds (RANS) appliquées aux écoulements à surface libre, nous avons obtenu une équation différentielle ordinaire (EDO) pour les vitesses ainsi que sa solution semi-analytique. L'application aux canaux étroits (rapport d'aspect < 5) montre des profils de vitesses avec un maximum au-dessous de la surface libre. Nous avons développé des formulations analytiques du nombre de Schmidt turbulent nécessaire à l'évaluation du transfert de masse et des concentrations au sein des milieux fluides en mouvement avec particules. Nous avons considéré le cas de particules de sable en suspension sous l'effet d'un écoulement oscillant sur un fond de rides. Le modèle analytique de viscosité turbulente a été généralisé à ce type d'écoulements et validé par le modèle à deux équations k-ω. Cet exemple a permis en particulier de distinguer l'effet de sédimentation du phénomène de détachement tourbillonnaire. En conclusion, cette présentation sera complétée par des perspectives à court et à moyen terme.