Transition to turbulence in suspension flows : numerical simulations and physical analysis
Transition à la turbulence des écoulements de suspension : simulations numériques et analyse physique
Résumé
This PhD addresses the influence of macroscopic and neutrally buoyant particles on the channel flows close to the laminar-turbulent transition regime. The suspension flow is moderately concentrated (solid volumetric concentration _ = 5%). The hydrodynamic coupling between the dispersed and carrier fluid is numerically resolved using the Force-Coupling Method approach. Particle trajectories are obtained by lagrangian tracking. In laminar wall-bounded flows, we show that the stress induced by the solid phase increases with inertia, and that the effect of the concentration is weaker than in the Stokes regime. The particles lag the flow and they migrate across the streamlines (Segré-Silberberg effect in Poiseuille flow). The migration and slip velocities are of the same order of magnitude for Rep = O(1). When the particles are freely suspended in a Poiseuille flow below the transition threshold, the initially homogeneous suspension (_ = 5%) becomes stratified after several ten time units (channel height/average flow velocity). After a hundred time units, the different rheological properties of the suspension segregated parts induce an instability yielding the formation of dune-like patterns, sustained by weak but finite secondary flows. In the fluctuating flow regime, we studied the profiles of the continuous and dispersed phase and realized 3D visualizations in order to understand why finite size particles delay the relaminarization threshold. The particles induce a significant increase of the velocity fluctuations in the transverse directions and they modify the rotational flow structures, which become smaller, more numerous and more energetic (larger rotation velocity). The wall-friction coefficient of the suspension flow in the transition regime is larger than the single-phase flow case. When the Reynolds number is decreased and the flow becomes laminar, the friction coefficient recovers the laminar law of a single phase flow provided that the fluid viscosity is replaced by the effective suspension viscosity in the Reynolds number definition. Our results clearly show that the two-phase channel flow turbulence is conserved down to a threshold well below the single phase flow limit, in agreement with the observations of Matas, Morris et Guazzelli (PRL, 2003) for a cylindrical geometry. In addition, we show that according to the transition direction, i.e. laminar 7! turbulent or turbulent 7! laminar, the switch from a regime to another does not occur at the same critical Reynolds number. Finally, in the limit of moderately concentrated (_ = 2.5−5%) suspension flow in a channel, the concentration has no significant influence on the critical Reynolds number.
Le travail de cette thèse aborde le sujet de l’influence des particules non-pesantes et de taille macroscopique sur les écoulements en canal dans des conditions proches du seuil de la transition laminaire-turbulent. Les suspensions sont faiblement concentrées (fraction volumique φ = 5%). Le couplage hydrodynamique existant entre la phase dispersée et la phase continue est résolu numériquement par la Force-Coupling Method, et les particules sont suivies d’une façon lagrangienne. Dans un écoulement laminaire de Couette ou de Poiseuille plan, nous montrons que les contraintes induites par la phase solide augmentent avec l’inertie, et l’influence de la concentration est plus faible qu’en régime de Stokes. Les particules avancent avec un retard dans la direction de l’écoulement et migrent à travers les lignes de courant (effet Segré-Silberberg en Poiseuille). Les vitesses de migration et de glissement s’amplifient avec l’inertie et sont du même ordre de grandeur quand Rep = O(1). Quand les particules sont lâchées librement dans un écoulement de Poiseuille plan en-deça du seuil critique de transition à la turbulence, la suspension initiale- ment homogène (φ = 5%) devient stratifiée, après un temps d’écoulement de plusieurs dizaines d’unités de temps (rapport de la hauteur du canal sur la vitesse moyenne de l’écoulement). Après une centaine d’unités de temps, nous observons le développement d’une instabilité à l’interface entre la zone chargée en particules et la zone de fluide pur. Des motifs dunaires prennent place dans la direction de la vorticité. Ces motifs sont soutenus par des écoulements secondaires d’intensités faibles mais non-nulles. Dans un écoulement au-dessus du seuil de transition, nous avons étudié les profils des phases continues et dispersées et réalisé des visualisations 3D afin de comprendre pourquoi les particules macroscopiques diminuent le nombre de Reynolds critique de relaminarisation de l’écoulement. Nous observons que les particules provoquent une augmentation significative des fluctuations de vitesses dans les directions transverses et qu’elles modifient les structures rotationnelles de l’écoulement, qui deviennent plus petites, plus nombreuses et plus énergétiques (plus grandes vitesses de rotation). Le coefficient de frottement pariétal de l’écoulement de suspension en régime de transition est supérieur à celui de l’écoulement monophasique. Quand le nombre de Reynolds est diminué et que l’écoulement devient finalement laminaire, le coefficient de frottement pariétal rejoint la loi laminaire d’un écoulement monophasique, à condition de substituer la viscosité effective de la suspension à la viscosité du fluide dans l’expression du nombre de Reynolds. D’après nos résultats, la turbulence de l’écoulement de suspension est conservée jusqu’à des nombres de Reynolds bien inférieurs à celui de l’écoulement monophasique en canal, en accord avec les observations ex- périmentales de Matas, Morris et Guazzelli (PRL, 2003) pour une géométrie cylindrique. Par ailleurs, nous montrons que selon le sens de la transition, laminaire → turbulent ou turbulent → laminaire, le nombre de Reynolds critique de transition d’un régime à l’autre n’est pas le même. Nous n’avons pas observé d’influence significative de la concentration en ce qui concerne la valeur du nombre de Reynolds critique de relaminarisation pour les deux concentrations étudiées (φ = 2.5% et 5%).
Origine : Version validée par le jury (STAR)