Ce projet de maîtrise vise à faire une caractérisation préliminaire de l'écoulement d'une turbine Francis de moyenne chute au régime sans charge dans le cadre du projet Tr-Francis au Laboratoire de machines hydrauliques. Concrètement, le projet a validé la méthodologie numérique utilisée à l'aide de mesures expérimentales préliminaires. De plus, les simulations numérique sont permis de fournir les chargements fluides afin que des simulations structurelles puissent être réalisées. Finalement, les résultats ont permis de cibler les structures dominantes dans l'écoulement qui causent les plus importantes fluctuations de pression et ainsi aider l'équipe expérimentale à positionner les capteurs de pression et à choisir leurs plans de mesures. Pour y arriver, la méthodologie est basée sur des interpolations de mesures effectuées par Hydro-Québec au régime sans charge sur le prototype. Les données sont prises telles quelles pour les simulations effectuées à l'échelle prototype et elles sont mise à l'échelle grâce aux lois de similitude de la norme IEC 60193 pour les simulations à l'échelle modèle. L'effet des conditions d'entrée dans le domaine est étudié sur des domaines partiels avec des maillages de différentes densités ainsi que sur le domaine complet. Les pressions obtenues numériquement dans l'aspirateur sont comparées aux mesures expérimentales préliminaires et les résultats sont concordants. De plus, les poches de cavitations observées numériquement au bord de fuite remontant sur le côté pression des aubes sont, de plus, confirmées par l'expérimentation. Des analyses fréquentielles à partir des signaux de pression sont utilisées dans ce projet afin de caractériser les phénomènes. En plus d'utiliser différentes méthodes de visualisation numériques afin d'isoler et d'analyser les structures principales dans l'écoulement. Imposer un profil de vitesse uniforme comparativement à un profil de vitesse tiré de simulations de la bâche et de la conduite d'amenée ne modifie pas les résultats significativement à l'échelle modèle avec le maillage le plus fin. Lors des simulations avec le domaine complet, plusieurs phénomènes d'importances ont été observés. Notamment des tourbillons inter-aubes, qui ont une modulation à f/n = 1, une importante zone de recirculation dans l'aspirateur qui remonte jusque dans la roue, des tourbillons à l'interface entre la roue et l'aspirateur qui se développent dans la couche cisaillée, un débalancement de l'aspirateur causé par le coude ainsi qu'une zone cavitante au bord de fuite. / This thesis presents the preliminary characterization of the flow in a medium head Francis turbine at speed no-load within the scope of the Tr-Francis project at the Hydraulic Machinery. Laboratory. Concretely, the project aims at validating the numerical methodology used with preliminary experimental measurements. Moreover, the numerical simulations will provide the fluid load for the FEA simulations. Ultimately, the results will allow identifying the dominant structures in the flow causing important fluctuations. Thus helping the experimental team to find the optimal location for the pressure sensors and the measuring planes. To do so, the numerical methodology is based on the measurements interpolation performed by Hydro-Québec on the prototype turbine at speed no-load. The data are applied, as they are, as the initial conditions for the prototype scale simulations. For the model scale simulations, they are scaled down using the similitude laws from the IEC 60193 standard. The effects of the inlet conditions are studied on partial domains with different mesh densities as well as on the complete domain. To validate the simulations, the pressure measurements obtained numerically in the draft tube are compared with the preliminary measurements and the results are in good agreement. Furthermore, the trailing edge cavitation observed numerically is also visible in experimental flow visualizations. Spectral analyses of pressure signals are used to help to characterize the phenomena. Different numerical visualization techniques are also used to isolate and analyse the main flow structures. Imposing a uniform velocity profile compared to the one obtained from the penstock and spiral case simulation does not significantly affect the results at model scale with the finest grid. Many important phenomena such as modulated inter-blades vortices, an important backflow region in the draft tube coming up to the runner, vortices generated in the draft tube in the shear layer, static pressure imbalance in the draft tube caused by the elbow and trailing blade cavitation were observed on the complete domain model scale simulations.

Ce mémoire de maîtrise en milieu industriel présente les travaux réalisés chez Maya Heat Transfer Technologies afin d'implémenter la technique de la simulation des grandes échelles dans le solveur de dynamique des fluides parallèle NIECE. Les modèles de sous-maille [sic] algébriques de Smagorinsky, WALE (Wale Adapting Local Eddy-Viscosity) et le modèle de Vreman ont été retenus et validés.Ce projet de maîtrise se divise [entre] quatre volets. La première partie présente l'implémentation et l'analyse d'un nouveau schéma de discrétisation spatial et temporel dans le contexte d'un algorithme de résolution des équations de Navier-Stokes par une approche mixte [de] volumes éléments finis sur maillage non structuré. Une étude comparative sur la stabilité, la diffusion et [la] dipersion [sic] numérique des nouveaux schémas et ceux déjà existants dans le code permet d'identifier le schéma centré de second ordre en espace, couplé au schéma de Crank-Nicolson en temps comme un compromis adéquat pour la simulation des grandes échelles dans un contexte industriel. Le second volet détaille les résulats [i.e. résultats] de la simulation de la décroissance d'une turbulence homogène isotrope. La décroissance des spectres d'énergie est présentée pour chacun des modèles de sous-maille, puis comparée à la solution sans modèle. L'effet du raffinement du maillage est investigué.Ce volet détaille également l'implémentation et la validation de l'algorithme permettant l'initialisation d'une turbulence isotrope incompressible. Une discussion sur les formulations analytiques des spectres d'énergie de Passot-Pouquet et de von Karman-Pao est finalement présentée. Le problème de la répresentation [i.e. représentation] de l'écoulement à l'entrée du domaine de calcul est ensuite abordé.Ce troisième volet détaille la méthode de génération de turbulence, appelée méthode des tourbillons synthétiques, servant à générer un champ de vitesse cohérent et représentatif d'une turbulence réelle. Une étude de l'influence des paramètres numériques de la méthode permet de dériver un critère automatique de sélection pour l'utilisateur. Pour terminer, une étude comparative de la méthode des tourbillons synthétiques à une méthode d'injection aléatoire est présentée. En dernier lieu, la réalisation du cas de validation du canal plan périodique permet d'investiguer la capacité du code à simuler les propriétés d'une turbulence cisaillée. Les résultats des modèles de sous-maille sont comparés au cas sans modèle de turbulence. L'effet du raffinement du maillage est étudié sur trois différentes grilles. Les statistiques de premier ordre (profil de vitesse moyenne) et les satistiques [i.e. statistiques] de second ordre (profils des composantes du tenseur de Reynolds) sont détaillées et analysées. Finalement, la contribution de l'énergie cinétique résolue et de l'énergie de sous-maille permet de quantifier la contribution des modèles de turbulence sur la solution.

Le développement des hydroliennes se base souvent sur des conditions d'écoulement idéalisées qui ne reflètent pas entièrement le courant présent sur un site réel. Dans cette optique, ce mémoire investigue l'effet d'un écoulement non aligné ou cisaillé sur les performances et le chargement de deux types d'hydrolienne : les hydroliennes à rotor axial (HRA) et à aile oscillante (HAO). Cette étude est réalisée à l'aide de simulations numériques et les résultats démontrent que pour les deux types d'hydrolienne, le désalignement produit une réduction de la puissance et de la traînée, tandis que le cisaillement n'a qu'un effet de second ordre sur ces derniers. Le chargement additionnel sur la structure de support causé par des conditions d'écoulement non idéal est aussi similaire pour les deux technologies. Le désalignement et le cisaillement affectent cependant plus significativement le chargement en fatigue des pales de l'HRA que celui de l'aile de l'HAO. / The development of hydrokinetic turbines is often based on idealized flow conditions which do not fully refect river or tidal currents. In this regard, this thesis investigates the effect of non-aligned or sheared flows on the performances and loading of two turbine types: the axial rotor (ART) and oscillating foil (OFT) turbines. This study was conducted with unsteady numerical simulations and the obtained results show that for both turbine types, misalignment produces a reduction in power and thrust, while the shear has a limited effect on those quantities. The additional loading on the support structure caused by the non-idealized flow is also similar for both devices. However, misalignment and shear affect more severely the fatigue loading on ART blades than they do for OFT.

Cette thèse porte sur la simulation numérique des écoulements bi-fluides par l'approche à interface diffuse. La description mathématique d'un écoulement bi-fluide par l'approche à interface diffuse consiste en une équation de Navier-Stokes modifiée couplée à un modèle de capture de l'interface mobile entre les deux fluides. Dans cette thèse, pour la capture de l'interface mobile nous avons porté notre attention sur le modèle de Allen-Cahn. En premier lieu nous nous sommes intéressés de prime abord à la résolution numérique du système semi-discrétisé en temps et totalement implicite Navier-Stokes/Allen-Cahn (NSAC). Pour ce faire nous avons développé un algorithme itératif basé sur la méthode du point fixe. Nous avons montré qu'à chaque itération, le système d'équations (contenu dans l'algorithme) est bien défini ; de plus, la solution de l'équation de Allen-Cahn satisfait le principe du maximum. Par la suite nous avons montré que l'algorithme de point fixe converge et que sa limite est la solution du système NSAC semi-discret. Si cette première méthode itérative nous a donné une méthode de résolution, elle n'est pas satisfaisante quant à la performance. En second lieu nous proposons un nouveau schéma de discrétisation en temps à pas fractionnaire inconditionnellement stable. Utilisant une approche de type point fixe (une projection couplée) nous avons montré la convergence à chaque pas de temps et que la limite correspond à la solution du système semi-discret (donnant ainsi l'existence et l'unicité de la solution). Nous concluons enfin avec des applications numériques aux fins d'illustrer la pertinence et les potentielles limitations du modèle d'une part, puis les performances de notre méthode de résolution du système Navier-Stokes/Allen-Cahn d'autre part.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / Ce mémoire de maîtrise porte sur la modélisation numérique de la problématique d'ingestion de gaz chauds dans la cavité rotor-stator du premier étage d'une turbine à gaz moderne. Dans le contexte de développement d'une aviation plus écologique, ce sujet mérite une attention marquée car des gains importants en ce qui a trait à l'efficacité globale du moteur peuvent être atteints par de meilleurs designs de joints rotor-stator limitant l'ingestion. La modélisation physique des phénomènes affectant l'ingestion devient alors essentielle afin de développer une meilleure compréhension des écoulements et mécanismes en présence. C'est dans cette optique que Pratt & Whitney Canada a mandaté le Laboratoire de Mécanique des Fluides Numérique (LMFN) de l'Université Laval afin de conduire une étude approfondie du phénomène d'ingestion. Ce mémoire a donc pour objectif principal de décrire le développement, l'application et la validation de certaines techniques de simulation numérique visant à caractériser qualitativement et quantitativement le phénomène d'ingestion observé à travers le joint et dans la cavité du groupe rotor-stator d'une turbine à gaz. Une investigation des limites de fiabilité d'une telle modélisation numérique dans un contexte réel d'opération du moteur est également présentée. La paramétrisation adéquate du problème et la validité d'extrapoler des résultats numériques d'ingestion obtenus aux conditions "banc d'essai" afin de prédire l'ingestion réelle aux conditions "moteur" font également partie des sujets abordés dans ce mémoire. La méthodologie développée considère deux phases de calculs menant à des diagnostics d'ingestion par le biais d'un marqueur passif. Une série de validations numériques et des comparaisons avec des résultats expérimentaux sont également présentées. La méthodologie de modélisation numérique est considérée robuste et fiable, mais seulement partiellement validée.

Thèse ou mémoire avec insertion d'articles / Titre de l'écran-titre (visionné le 11 janvier 2024) / Cette thèse de doctorat étudie le transport de fluides viscoplastiques à travers des canaux dotés de parois superhydrophobes striées, en utilisant une approche de modélisation complète et des simulations numériques. L'objectif est de comprendre les interactions entre le fluide viscoplastique et la surface superhydrophobe, et d'identifier les facteurs qui influencent le comportement de l'écoulement. On suppose que les rainures emprisonnent des poches d'air à la surface superhydrophobe, maintenant l'interface liquide/air plate et fixée aux bords des rainures, ce qui entraîne une condition de glissement-collage. L'analyse est complète, couvrant à la fois les régimes de fluage et inertiel dans les limites des canaux épais et minces. En particulier, six paramètres, notamment le nombre de Reynolds (*R*), le nombre de Bingham (*B*), le nombre de glissement (*b*), la longueur de périodicité des rainures (ℓ), la fraction de zone de glissement (*φ*), et l'angle d'orientation des rainures (*θ*), sont pris en compte. Ensuite, les effets de ces paramètres sont analysés sur les variables d'écoulement d'intérêt, telles que le champ de perturbation, le champ de vitesse, la magnitude du taux de déformation, les zones de bouchon non cédées, la longueur de glissement effective, le facteur de frottement, l'asymétrie de l'écoulement, le mélange, et les classifications de régime. Dans les canaux épais (ℓ ≪ 1) avec des rainures transversales (*θ* = 90°), l'augmentation de *B* et de *b* entraîne une croissance des champs de perturbation et de glissement, qui deviennent asymétriques en augmentant *R*. Dans certaines conditions d'écoulement, le bouchon central se rompt ou une zone de bouchon non cédée peut se former à l'interface liquide/air de la paroi superhydrophobe, ce qui nous permet de classifier les régimes d'écoulement. Dans les canaux minces (ℓ ≫ 1), l'augmentation de *b* ou de *θ* entraîne la déformation du bouchon central non cédé, sauf pour les écoulements avec des rainures longitudinales (*θ* = 0). Pour les rainures obliques (0 < *θ* < 90°), un écoulement secondaire apparaît, déclenchant le mélange de l'écoulement. Enfin, l'analyse de stabilité linéaire de la condition de glissement homogène (*φ* = 1) révèle les effets stabilisants/déstabilisants des conditions de glissement dans le sens du courant/transversales. / This Ph.D. thesis investigates the transport of viscoplastic fluids through channels with grooved superhydrophobic (SH) walls, using a comprehensive modeling approach and numerical simulations. The goal is to understand the interactions between the viscoplastic fluid and the superhydrophobic surface and identify factors that influence the flow behaviour. It is assumed that the grooves trap air pockets on the SH surface, keeping the liquid/air interface flat and pinned at the groove edges, resulting in a slip-stick condition. The analysis is comprehensive, covering both creeping and inertial regimes in thick and thin channel limits. In particular, six parameters, including the Reynolds number (*R*), Bingham number (*B*), slip number (*b*), groove periodicity length (ℓ), slip area fraction (*φ*), and groove orientation angle (*θ*), are considered. Subsequently, the effects of these parameters are analyzed on the flow variables of interest, such as the perturbation field, velocity field, strain-rate magnitude, unyielded plug zones, effective slip length, friction factor, flow asymmetry, mixing, and regime classifications. In thick channels (ℓ ≪ 1) with transverse grooves (*θ* = 90°), increasing *B* and *b* leads to growing the perturbation and slip velocity fields, which become asymmetric by increasing *R*. Under certain flow conditions, the center plug breaks or an unyielded plug zone can form on the SH wall liquid/air interface, allowing us to classify the flow regimes. In thin channels (ℓ ≫ 1), increasing either *b* or *θ* causes the the unyielded center plug to deform and eventually break, except for flows with longitudinal grooves (*θ* = 0). For the oblique grooves (0 < *θ* < 90°), a secondary flow appears, triggering flow mixing. Finally, the linear stability analysis of the homogeneous slip condition (*φ* = 1) reveals stabilizing/destabilizing effects of the streamwise/spanwise slip conditions.

Cette étude s'intéresse à une méthode vortex lagrangienne 2-D permettant de simuler des écoulements incompressibles externes autour de multiples corps en mouvement arbitraire. L'implantation initiale de la méthode a été développée à 1' Université Catholique de Louvain puis des modifications ont été apportées en collaboration avec le Laboratoire de Mécanique des Fluides Numérique pour prendre en compte le mouvement et la présence simultanée de plusieurs corps. Cependant, ces modifications n'étaient pas validées et aucune méthode de calcul des forces et du moment sur les corps individuellement n'était disponible. Au sein de cette étude, plusieurs approches de calcul des forces et du moment sont présentées et détaillées, toutes ayant la particularité de ne pas nécessiter la connaissance du champ de pression de l'écoulement. La première approche utilise des surfaces de contrôle englobant chaque corps ainsi que la surface même des corps. La seconde est développée à partir d'une surface de contrôle à l'infini et ne peut fournir que les forces nettes totales sur l'ensemble des corps en présence. La troisième approche ne requiert que de l'information (vitesse et vorticité) à la surface des corps et elle est la seule pour l'instant qui permet de fournir également le moment de force. Toutes ces approches ainsi que les modifications effectuées pour prendre en compte le mouvement sont validées pour des écoulements à bas Reynolds autour de corps en translation et rotation combinés. Les cas d'un cylindre unique, de deux cylindres en tandem et d'une aile oscillante sont étudiés. Les comparaisons sont effectuées à l'aide de résultats produits par le code commercial eulérien de volumes finis Fluent 6.2. Les résultats mettent en évidence l'habilité de la méthode vortex à capturer précisément la dynamique des fluides, ainsi que l'efficacité de chacune des méthodes de calcul des forces et du moment. / This work is about a 2-D lagrangian vortex method allowing simulation of external incompressible flow around arbitrary moving bodies. The initial implementation of the method was developed at Université Catholique de Louvain and adapted for multiple moving bodies in collaboration with Laboratoire de Mécanique des Fluides Numérique. However, these modifications were not validated and no method of force and moment calculation on individual body was available. In this study, several approaches of force and moment calculation are presented and detailed, ail of which do not require information on the pressure field. The first approach uses control surfaces surrounding each body as well as the body surfaces themselves. The second one is developed from a control surface taken at infinity and can only provide total net forces on ail bodies taken together. The third one requires information (velocity and vorticity) only at body surfaces. The latter is the only approach considered up to now that can provide moment of force on the bodies. Ail three methods as well as ail the modifications that were implemented to adapt the solver to moving bodies are validated for low Reynolds number flows around bodies undergoing combined translation and angular motion. Cases of a single cylinder, two cylinders in tandem, and oscillating airfoil are studied. Comparisons are then made with results obtained with the eulerian finite volume commercial code Fluent 6.2. Results show the ability of the vortex method to correctly capture fluid dynamics, as well as the efficiency of each method for the force and moment calculations.

Les suspensivores ont la tâche importante de séparer les particules de l'eau. Bien qu'une grande gamme de morphologies existe pour les structures d'alimentation, elles sont pratiquement toutes constituées de rangées de cylindres qui interagissent avec leur environnement fluide. Le mécanisme de capture des particules utilisé dépend des contraintes morphologiques, des besoins énergétiques et des conditions d'écoulement. Comme nos objectifs étaient de comprendre ces relations, nous avons eu recours à des études de comparaison pour interpréter les tendances en nature et pour comprendre les conditions qui provoquent de nouveaux fonctionnements. Nous avons utilisé la dynamique des fluides numérique (computational fluid dynamics, CFD) pour créer des expériences contrôlées et pour simplifier les analyses. Notre première étude démontre que les coûts énergétiques associés au pompage dans les espaces petits sont élevés. De plus, le CFD suggère que les fentes branchiales des ptérobranches sont des structures rudimentaires, d'un ancêtre plus grande. Ce dernier point confirme l'hypothèse qu'un ver se nourrit par filtration tel que l'ancêtre des deuterostomes. Notre deuxième étude détermine la gamme du nombre de Reynolds number critique où la performance d'un filtre de balane change. Quand le Re est très bas, les différences morphologiques n'ont pas un grand effet sur le fonctionnement. Cependant, une pagaie devient une passoire lorsque le Re se trouve entre 1 et 3,5. Le CFD s’est dévoilé être un outil très utile qui a permis d’obtenir des détails sur les microfluides. Ces études montrent comment la morphologie et les dynamiques des fluides interagissent avec la mécanisme de capture ou de structures utilisées, ainsi que comment des petits changements de taille, de forme, ou de vitesse d'écoulement peuvent conduire à un nouveau fonctionnement. / Suspension feeders have the important task of separating particles from the liquid medium surrounding them. Although a wide range of morphologies exist for feeding structures, essentially all consist of arrays of cylinders interacting with their fluid environment. The particle capture mechanism employed depends on morphological constraints, energy requirements, and flow conditions, therefore our objectives were to understand these relationships through comparison studies to interpret trends in nature and to understand the conditions in which novel functioning arises. Our studies used computational fluid dynamics (CFD) to create controlled experiments and to simplify analyses. The first study demonstrates the unfeasibly high energetic costs of pumping in small and medium pharynx sizes of hemichordates, implying the gill pores of pterobranchs are likely vestigial structures from a larger ancestor. This last point further supports the hypothesis of an acorn worm as the ancestor to all deuterostomes. Our second study determined the critical Reynolds number (Re) range where barnacle filters transition from a paddle to a sieve. At very low Re, morphological differences have no major effect on functioning and filters behave as paddles, however, at Re 1 - 3.5, these differences dictate when a paddle becomes a sieve. CFD proved to be a very useful tool for simplifying studies and providing detailed microfluidics. These studies demonstrate how morphology and fluid dynamics interact to dictate the capture mechanisms or appendages employed and how simple changes in size, shape, or flow speed can lead to novel functioning.

Dans cette dissertation, nous considérons l’écoulement des liquides conducteurs d’électricité dans un champ magnétique externe. De tels écoulements sont décrits par les équations de la magnétohydrodynamique (MHD) quasi-statique, et sont fréquemment rencontrés dans des applications pratiques. Il suit qu’il y a un intérêt fort pour des outils numérques qui peuvent simuler ces écoulements dans des géometries complexes.<p>La première partie de cette thèse (chapitres 2 et 3) est dédiée à la présentation de la machinerie numérique qui a été utilisée et implémentée afin de résoudre les équations de la MHD quasi-statistique (incompressible). Plus précisément, nous avons contribué au développement d’un solveur volumes finis non-structuré parallèle. La discussion sur ces méthodes est accompagnée d’une analyse numérique qui est aussi valable pour des mailles non-structurées. Dans le chapitre 3, nous vérifions notre implémentation par la simulation d’un certain nombre de cas tests avec un accent sur des écoulements dans un champ magnétique intense.<p>Dans la deuxième partie de cette thèse (chapitres 4-6), nous avons utilsé ce solveur pour étudier des écoulements MHD de proche paroi .La première géometrie considérée (chapitre 4) est celle d’une conduite circulaire infini d’axe à haut nombre de Hartmann. Nous avons investitgué la sensitivité des résultats numériques au schéma de discrétisation et à la topologie de la maille. Nos résultats permettent de caractériser in extenso l’écoulement MHD dans une conduite avec des bords bien conducteurs par moyen des lois d’échelle.<p>Le sujet du cinquième chapitre est l’écoulement dans une conduite toroïdale à section carée. Une étude du régime laminaire confirme une analyse asymptotique pour ce qui concerne les couches de cisaillement. Nous avons aussi effectué des simulations des écoulements turbulents afin d’évaluer l’effet d’un champ magnétique externe sur l’état des couches limites limites.<p>Finalement, dans le chapitre 6, nous investiguons l’écoulement MHD et dans un U-bend et dans un coude arrière. Nous expliquons comment générer une maille qui permet de toutes les couches de cisaillement à un coût computationelle acceptable. Nous comparons nos résultats aux solutions asymptotiques. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Les suspensivores ont la tâche importante de séparer les particules de l'eau. Bien qu'une grande gamme de morphologies existe pour les structures d'alimentation, elles sont pratiquement toutes constituées de rangées de cylindres qui interagissent avec leur environnement fluide. Le mécanisme de capture des particules utilisé dépend des contraintes morphologiques, des besoins énergétiques et des conditions d'écoulement. Comme nos objectifs étaient de comprendre ces relations, nous avons eu recours à des études de comparaison pour interpréter les tendances en nature et pour comprendre les conditions qui provoquent de nouveaux fonctionnements. Nous avons utilisé la dynamique des fluides numérique (computational fluid dynamics, CFD) pour créer des expériences contrôlées et pour simplifier les analyses. Notre première étude démontre que les coûts énergétiques associés au pompage dans les espaces petits sont élevés. De plus, le CFD suggère que les fentes branchiales des ptérobranches sont des structures rudimentaires, d'un ancêtre plus grande. Ce dernier point confirme l'hypothèse qu'un ver se nourrit par filtration tel que l'ancêtre des deuterostomes. Notre deuxième étude détermine la gamme du nombre de Reynolds number critique où la performance d'un filtre de balane change. Quand le Re est très bas, les différences morphologiques n'ont pas un grand effet sur le fonctionnement. Cependant, une pagaie devient une passoire lorsque le Re se trouve entre 1 et 3,5. Le CFD s’est dévoilé être un outil très utile qui a permis d’obtenir des détails sur les microfluides. Ces études montrent comment la morphologie et les dynamiques des fluides interagissent avec la mécanisme de capture ou de structures utilisées, ainsi que comment des petits changements de taille, de forme, ou de vitesse d'écoulement peuvent conduire à un nouveau fonctionnement. / Suspension feeders have the important task of separating particles from the liquid medium surrounding them. Although a wide range of morphologies exist for feeding structures, essentially all consist of arrays of cylinders interacting with their fluid environment. The particle capture mechanism employed depends on morphological constraints, energy requirements, and flow conditions, therefore our objectives were to understand these relationships through comparison studies to interpret trends in nature and to understand the conditions in which novel functioning arises. Our studies used computational fluid dynamics (CFD) to create controlled experiments and to simplify analyses. The first study demonstrates the unfeasibly high energetic costs of pumping in small and medium pharynx sizes of hemichordates, implying the gill pores of pterobranchs are likely vestigial structures from a larger ancestor. This last point further supports the hypothesis of an acorn worm as the ancestor to all deuterostomes. Our second study determined the critical Reynolds number (Re) range where barnacle filters transition from a paddle to a sieve. At very low Re, morphological differences have no major effect on functioning and filters behave as paddles, however, at Re 1 - 3.5, these differences dictate when a paddle becomes a sieve. CFD proved to be a very useful tool for simplifying studies and providing detailed microfluidics. These studies demonstrate how morphology and fluid dynamics interact to dictate the capture mechanisms or appendages employed and how simple changes in size, shape, or flow speed can lead to novel functioning.