Titre de l'écran-titre (visionné le 6 novembre 2023) / Le but est d'étudier l'athérosclérose comme une application de la mécanique des fluides qui porte sur l'écoulement sanguin dans une artère sténosée. Le sang et la forme de l'accumulation des graisses sur la paroi artérielle jouent un rôle fondamentale dans le comportement de l'écoulement. Pour se mettre dans le cadre mathématique, un survol sur la mécanique des fluides permet d'obtenir la formulation de Navier-Stokes que l'on résout numériquement à l'aide des schémas de projection. L'adaptation de maillage est un outil numérique performant pour un tel type de calcules. Les conditions aux limites en pression à l'entrée et à la sortie de l'artère sont plus recommandés pour ce type de problèmes, mais leurs insertion exige des espaces d'interpolation particuliers que l'on explore en parcourant la littérature correspondante.

Ce mémoire de maîtrise présente l'étude réalisée an d'obtenir une mesure de référence du frottement pariétal dans le canal hydrodynamique conçu par Coulaud (2013). Ce canal est destin é à étalonner des capteurs de frottement pariétal afin de pouvoir obtenir un meilleur aperçu de l'écoulement dans les turbines hydrauliques. Le canal hydrodynamique du laboratoire doit alors représenter un écoulement de Poiseuille plan dont le frottement pariétal moyen est connu. Ainsi, le frottement pariétal moyen a été obtenu par une mesure du gradient longitudinal de pression. Utilisée lors des essais préliminaires de Coulaud (2013), cette méthode a été améliorée par l'emploi d'un manomètre à 13 colonnes d'eau. À cause du rapport de forme du canal, le frottement pariétal moyen obtenu par le gradient longitudinal de pression doit être corrigé (Vinuesa et al., 2015) pour devenir une valeur de référence. Afin de vérifier ces résultats, une autre méthode de mesure du frottement pariétal basée sur les fluctuations des composantes débitante et transversale de la vitesse a été mise en place. Un système de vélocimétrie laser à effet Doppler (LDV) a été utilisé pour ces expériences. Avec les deux méthodes de mesure, la tendance d'évolution du frottement pariétal moyen en fonction du nombre de Reynolds est identique. Par contre, un écart significatif entre les valeurs obtenues est mesuré. En effet, par la méthode utilisant le LDV, la valeur du frottement pariétal moyen est toujours supérieure à celle issue du gradient longitudinal de pression. Dans l'objectif de confirmer que l'écoulement dans la section centrale du canal était conforme à l'écoulement de Poiseuille plan, le même système LDV a été utilisé. En effet, les profils transversaux des grandeurs moyennes et des niveaux de turbulence des composantes débitante et transversale de la vitesse ont montré que l'écoulement est pleinement développé et bidimensionnel. / This Master's thesis presents the study conducted to obtain a reference measurement of the wall shear stress in the hydrodynamic channel designed by Coulaud (2013). This channel is destined to calibrate wall shear stress sensors in order to get a better understanding of the flow in hydraulic turbines. The laboratory's hydrodynamic channel therefore needs to recreate a planar Poiseuille flow in which the mean wall shear stress is known. Thus, mean wall shear stress was obtained by measuring longitudinal pressure gradient. Used during preliminary tests of Coulaud (2013), this method was improved with a 13 water columns manometer. Because of the aspect ratio of the channel, wall shear stress obtained by this method needs a correction (Vinuesa et al., 2015) to provide a reference value. In order to verify those results, another wall shear stress measurement method based on fluctuations of longitudinal and transverse velocity components was performed. A laser Doppler velocimetry system was used for these experiments. Both methods present the same trends. However, a significant deviation between obtained values is measured. Indeed, with the LDV method, mean wall shear stress values are always greater than those from longitudinal pressure gradient method. To conform that the flow in the central core of the test section was equivalent to a planar Poiseuille flow, the LDV system was used. Indeed, transverse profiles of mean and turbulent levels of longitudinal and transverse velocity components demonstrate that the flow is fully developped and two-dimensional.

Titre de l'écran-titre (visionné le 24 juillet 2023) / Ce mémoire vise à identifier et caractériser les phénomènes d'écoulement dans le cône de la turbine modèle à moyenne chute opérée en marche à vide à la vitesse de synchronisation et installée au laboratoire Heki de l'Université Laval. Avec peu de recherches réalisées sur le régime en marche à vide et des machines fonctionnant de plus en plus dans ce régime, le mémoire aidera l'industrie à comprendre par quels mécanismes la machine est endommagée par le régime en marche à vide. Les mesures expérimentales ont utilisé la vélocimétrie par image de particules stéréoscopique dans un plan radial-axial du diffuseur conique. Les mesures étaient calées sur la position de la roue de turbine de façon verrouillé en phase). Cette méthode de mesure a été appliquée à deux hauteurs dans le diffuseur conique pour couvrir la majeure partie de sa hauteur. Les phases mesurées sont équidistantes d'un intervalle de 6° et couvrent 360°. Pour la section de mesure supérieure, une résolution plus fine de 3° a été utilisée entre les phases 300° et 360°. La structure d'écoulement globale discutée dans la littérature s'est avérée présente dans le modèle étudié et a été mis en évidence. Un débit à haute vitesse avec une grande vitesse tangentielle et axiale a été mesurée près de la paroi interne du cône. Une zone de recirculation avec des vitesses beaucoup plus faibles a été mesurée au centre du diffuseur. La couche de cisaillement mentionnée par quelques auteurs a été caractérisée par le calcul des composantes de cisaillement. L'écoulement s'est avéré ne pas être axisymétrique. Des traces de tourbillons inter-aubes dans le diffuseur conique ont été trouvées dans les mesures à moyenne de phase, comme cela avait été prédit par les simulations numériques. / This thesis aims to characterize the flow and identify flow phenomena in the draft tube of the Tr-Francis medium head model turbine operated at speed-no load and installed at the Heki laboratory at Laval University. With little research accomplished on the speed no-load regime and machines being operated increasingly in this regime, the thesis will help the industry understand by what mechanisms the machine is damaged by the no load regime. The experiments used stereoscopic particle image velocimetry in one radial-axial plane of the conical diffuser. The measurements were locked on the position of the turbine runner (phase-locked). This measurement method was applied at two heights in the conical diffuser to cover most of its height. The measured phases are equally spaced at an interval of 6° and cover 360°. For the upper measurement section, a finer resolution of 3° was used between phases 300° and 360°. The global flow structure discussed in the literature was found to be present in the studied model. A high velocity throughflow with a large tangential and axial velocity was measured near the diffuser's inner wall. A recirculation zone with much lower velocities was measured at the center of the diffuser. The shear layer mentioned by a few authors was characterised by calculating shear strain components. The flow was found not to be axisymmetric. Traces of inter-blade vortices in the draft tube were found in the phase-averaged measurements as was predicted by the numerical simulations.

Nous présentons la généralisation d'un modèle de dépôt de particules aux écoulements non isothermes. Ce travail comporte deux parties. La première partie présente l'approche de modélisation lagrangienne stochastique adoptée pour simuler le dépôt de particules colloïdales et traite du raccord entre le Modèle de Langevin Généralisé qui calcule le transport des particules dans le cœur de l'écoulement et le Modèle de Proche Paroi destiné au transport des particules au voisinage de la paroi. Le raccord est assuré par l'équilibre du flux de particules traversant l'interface entre les deux modèles. La validation du couplage s'effectue par le passage en particules fluides, il s'agit de vérifier que le raccord est valable dans le cas limite du fluide. En effet, des accumulations non physiques de particules fluides peuvent se produire autour de l'interface. Le travail de couplage a donné lieu à des développements théoriques et numériques sur les statistiques conditionnées d'une particule traversant une interface. Des calculs de dépôt de particules montrent que le modèle reproduit correctement les résultats expérimentaux pour des écoulements isothermes. Cependant, dans un écoulement non isotherme, les particules peuvent faire l'expérience du phénomène de thermophorèse (qui est le mécanisme de déplacement des particules induit par un gradient de température) qui peut considérablement influencer le dépôt de particules. Nous avons donc choisi d'intégrer la modélisation de la thermophorèse dans le Modèle de Proche Paroi afin d'évaluer au mieux le taux de dépôt de particules dans un écoulement non isotherme. Nous commençons par étudier la thermophorèse en milieu gazeux car elle est très bien expliquée par la théorie cinétique des gaz. L'étude de la thermophorèse en milieu liquide est plus ardue car elle ne dispose pas encore de théorie unifié ni de travaux expérimentaux qui étudient ce phénomène. La dernière partie de ce document est donc consacrée à la prise en compte du phénomène de thermophorèse en liquide dans le Modèle de Proche Paroi. Cette étude fournit de nouveaux résultats numériques de dépôt de particules en présence d'un gradient de température dans un milieu liquide.

Blood is a complex suspension of deformable particles, red blood cells, which exhibits a sophisticated dynamics when flowing in the microvasculature. Most of these complex phenomena, non-linear rheology, structuration of the suspension, heterogeneities of the hematocrit distribution, are directly connected to the rich microscopic dynamics of individual red blood cells, and their hydrodynamics interactions. We investigate a few aspects of the dynamics of red blood cells and giant vesicles - a simple model for RBCs. A study on the dynamics of very deflated vesicles, with shapes similar to those of red blood cells, shows that these objects which haven't received a lot of attention so far can exhibit richer than expected dynamics. We then mainly focus on the still unexplored problem of hydrodynamic interactions between vesicles or red blood cells and their consequences at the scale of the suspension. An experimental study of the interaction of two identical vesicles in shear flow shows that there is a net repulsion between the cells that leads to an increase of the distance between vesicles in a suspension. Scaling arguments are proposed for this interaction and a comparison with numerical results is performed and a quantitative estimation of a shear induced diffusion coefficient obtained by averaging the results for pair interactions is found. Finally, we investigate the diffusion of a cloud of red blood cells in Poiseuille flow in order to directly determine diffusion coefficients. The experiment shows that the cloud widens when traveling along the channel with a power law behaviour indicating sub-diffusion. This effect is confirmed by a theoretical analysis of a few limit cases.

Ce mémoire est consacré à l'étude des effets des interactions hydrodynamiques sur le mouvement et le transport des particules sphériques et cylindriques dans les fluides non newtoniens. Des résultats importants ont été obtenus pour expliquer la physique des suspensions en fluide non newtonien. Le modèle non newtonien choisi est celui d'Ostwald car il décrit bien les effets de rhéofluidification et de rhéoépaississement qui caractérisent la plupart des fluides non newtoniens. Le premier résultat a été de montrer comment l'indice de fluidité affecte le comportement des suspensions par le biais de la longueur d'écran hydrodynamique autour de chaque particule. Ceci nous a permis de donner les bonnes valeurs de la traînée subie par chaque particule en milieu infini, et de montrer l'apparition du paradoxe de Whitehead dans le cas de la sphère à partir de n=2 et de Stokes dans le cas du cylindre à partir de n=1. Lorsque n est voisin de ces valeurs critiques, la détermination des forces hydrodynamiques devient très sensible à l'inertie. Le deuxième résultat important pour l'industrie de l'injection des matériaux composites, a été de montrer par une méthode inverse que les interactions hydrodynamiques pouvaient induire un retard plus ou moins important par rapport au cas newtonien dans le transport de particules. Le troisième résultat important pour l'analyse des processus d'agrégation de particules sphériques ou cylindriques a été obtenu dans le cas d'une particule sphérique ou cylindrique en sédimentation vers un plan fixe. Un calcul asymptotique dans le cas non newtonien en régime de lubrification, comparé avec succès à celui obtenu numériquement par la méthode de maillage dynamique, nous a permis d'obtenir les lois d'évolution de la force subie par ces particules entrant en contact avec un plan. Ces derniers résultats trouvent une application dans l'utilisation des machines dynamiques à force de surface dans le cadre de la nanorhéologie.

Dans ces travaux , nous avons deux principaux objectifs physique et numérique. Le problème physique consiste à trouver la solution de Rayleigh-Bénard pour des fluides newtoniens et non-newtoniens. Dans la présente étude, une présentation générale des résultats de la convection de Rayleigh-Bénard (RBC) est donnée dans le cas des fluides newtoniens et non-newtoniens tels que des fluides rhéofluidifiants modélisés par la loi puissance et des fluides viscoplastiques (fluides de Bingham, Herschel-Bulkley et Casson), en régime permanent et transitoire. Dans le cas des fluides viscoplastiques, les modèles macroscopiques ne prenant pas bien en compte la réalité physique de la contrainte seuil ont fait l'objet d'une modélisation. Un modèle mesoscopique proposé par Hébraud et Lequeux a été utilisé. Le problème numérique consiste à développer la méthode de résolution PGD (Proper Generalized Decomposition) pour résoudre les modèles non linéaires couplés transitoires, dans le cas du problème de Rayleigh-Bénard. Cette méthode est également utilisée pour résoudre le problème RBC paramétrique en y ajoutant quelques variables physiques comme coordonnées supplémentaires. Par ailleurs, dans le cas des fluides non-newtoniens, nous avons utilisé la PGD pour résoudre les équations mesoscopiques et macroscopiques couplées.

Le prototypage virtuel des circuits pneumatiques de puissance, par exemple les circuits de freinage des véhicules industriels, constitue un enjeu important en raison de la complexité des écoulements en régime transsonique et des couplages entre les échelles locales et macroscopiques. Ces problèmes sont rencontrés lors de la conception, de la synthèse des commandes et de l'analyse des performances statiques et dynamiques de ces circuits et l'analyse. La mise au point des modèles numériques de ces systèmes induit des coûts et des temps importants par rapport à d'autres systèmes. La démarche proposée dans cette thèse repose sur la construction numérique de bases de données permettant de caractériser le comportement local et macroscopique d'un composant de circuit en fonction de la variation de certains paramètres physiques ou géométriques par rapport à un point de fonctionnement de référence. Les bases de données résultent de l'extrapolation de la solution des équations de Navier Stokes moyennées (RANS) pour le point de référence considéré obtenu à l'aide d’un logiciel de paramétrisation en mécanique des fluides (Turb’Opty). La contribution de cette thèse repose pour l'essentiel dans un travail d'analyse des solutions issues de la paramétrisation dans deux contextes différents: la tuyère De Laval et un élément "coude", des composants élémentaires de circuit. Nous avons montré que ces exemples "simples" conduisent déjà à des difficultés importantes en termes de paramétrisation du problème et du calcul des dérivées des champs aérodynamiques en raison de la taille du problème. Pour pallier cette difficulté, nous avons proposé de déraffiner le maillage et nous avons alors montré que cette démarche conduit parfois à déplacer ou à atténuer certains phénomènes (chocs). La deuxième contribution de ce travail repose sur l'évaluation de la qualité des solutions extrapolées, de leur domaine de validité et la construction des liens entre grandeurs locales et macroscopiques. Nous avons enfin proposé une démarche permettant de reconstruire la caractéristique en débit d'un composant à partir de la détermination de la solution extrapolée pour un nombre limité de points de référence. / Virtual prototyping pneumatic circuits for power transmission, for example braking circuits of trucks, is still a difficulty because of the complexity of the flow behavior in transonic conditions and of the coupling between local and macroscopic scales. These problems are met during system design, control synthesis and for static and dynamic performance analysis. Tuning accurate numerical models requires important costs and time when compared to other systems. The methodology proposed in this PhD thesis relies on numerically determining a data base that characterizes the local and macroscopic behavior of a circuit component according the variation from a reference point of some physical or geometrical parameters. The data bases are obtained from the extrapolation of the Mean Navier Stokes solution (RANS) for a given reference point with the help of a parametrization software dedicated to fluid mechanics (Turb'Flow). The main contribution of this thesis relies io the analysis of the solution obtained from the parametrization in two different cases: the De Laval nozzle and un "elbow" connecting element, which are elementary component in a circuit. We have shown that these two "simple" cases lead already to important difficulties in term of problem parametrization and calculation of the derivatives of the aerodynamic fields because of the problem dimension. In order to tackle this, we proposed to reduce the spatial discretization (mesh derefining) and we showed that this approach could sometimes lead to damp or move some phenomena (shocks). The second contribution of this work relies on evaluating the quality of the extrapolated solution and their validity domain, and on building links between local and macroscopic behavior. Finally, we proposed a method that allows the mass flow rate characteristic of a component to be determined from the calculation of the extrapolated solution issued from a limited number of reference points.

La présente étude s'intéresse à un écoulement d'acoustic streaming, c'est-à-dire un écoulement généré par la propagation d'une onde acoustique dans un fluide. Le travail consiste à comparer deux approches : expérimentale et numérique. Les ultrasons sont émis à 2MHz par un transducteur piézo-électrique de 28.5mm de diamètre. Ce dernier est plongé dans une cuve d'eau équipée de deux parois absorbantes : l'une sert à séparer le champ proche du champ lointain et l'autre est placée à l'extrémité du domaine fluide afin d'éviter toute réflexion. On réalise ainsi une étude en champ proche et une étude en champ lointain. Les mesures sont de deux types : champ de pression acoustique (hydrophone) et champ de vitesse (PIV). En parallèle, on effectue des simulations numériques directes avec le logiciel StarCCM+TM. Il s'agit de résoudre les équations de Navier-Stokes en fluide incompressible complétées d'un terme source de force acoustique. L'expression de ce dernier est obtenue par séparation des échelles de temps, ce qui consiste à négliger à l'échelle de temps acoustique les variations temporelles lentes, de l'écoulement généré. La démarche est ensuite analogue à celle utilisée en turbulence pour le calcul des tenseurs de Reynolds. On obtient finalement un bon accord entre les résultats expérimentaux et ceux de la modélisation numérique.

Les phénomènes d'érosion sont une cause majeure dans la rupture des ouvrages hydrauliques tels que les digues ou les barrages. Le territoire compte 10000 km de digues, d'où un enjeu important en terme de risque d'inondation. Une approche plus générale est essentielle à une compréhension fine des phénomènes. La modélisation des mécanismes de l'érosion doit permettre l'élaboration de solutions pour sécuriser les ouvrages et les rendre moins sensibles à ces phénomènes. Ce travail est consacré à la modélisation d'une interface entre un milieu de sol et un milieu d'écoulement. Le mouvement de cette interface est régi par l'érosion du sol sous l'effet de l'écoulement. L'écoulement est calculé par les équations de Stokes munies d'un terme de pénalisation permettant la représentation unifiée des comportements des deux milieux. Dans une première partie, nous aborderons la problématique de la représentation de l'érosion. La seconde partie traite plus précisément de la représentation unifiée des comportements des milieux, notamment grâce à l'utilisation des méthodes de pénalisation. La troisième partie constitue une revue des principales méthodes permettant le suivit d'une interface. Une description précise de la méthode Level Set est donnée. Dans la quatrième partie, nous abordons la résolution numérique des équations d'écoulement. La cinquième partie rassemble l'ensemble des tests de la méthode. Enfin, la dernière partie est consacrée à l'application du modèle sur trois cas : le " Hole Erosion Test ", le phénomène de suffusion et l'érosion d'un lit de rivière autour d'une pile de pont.