Dans ce travail, nous avons déterminé le diagramme de phases liquide-vapeur par la méthode des équations intégrales pour des modèles d interactions représentant non seulement des liquides simples, mais aussi des liquides complexes grâce au concept de ¡gros grains¡ . Notre choix s est porté sur une nouvelle formulation de l équation intégrale de Sarkisov qui a une bonne cohérence thermodynamique sans paramètre ajustable. La résolution du système d équations intégrales est opéré par la méthode de Labik et al., alors que le calcul des dérivées des fonctions de corrélation est réalisé de manière formellement exact grâce à un outil de différentiation automatique de code par la méthode linéaire tangent. Enfin, l exploration efficace du diagramme de phases est rendue possible par une technique adaptative garantissant une détermination précise des lignes binodales et spinodales ainsi que du point critique. Dans la première phase, nous avons mis au point et validé notre approche pour un fluide de sphères dures et pour celui de Lennard-Jones. Nous avons ensuite étudié les propriétés structurales, les grandeurs thermodynamiques et le diagramme de phases de deux familles de potentiels¡: la première regroupant les potentiels de type Yukawa, la seconde correspondant aux potentiels discrets. Dans chacun de ces cas, les résultats ont été validés par comparaison avec les résultats de simulations numériques disponibles dans la littérature. Nous avons également pu étudier l évolution du diagramme de phases avec les caractéristiques du potentiel, ainsi que l apparition d un transition liquide-liquide dans un cas particulier de potentiels discrets. Enfin, nous avons aussi entendu notre étude aux systèmes métalliques dont les potentiels effectifs dépendent de la densité. / N this work we determined the liquid-vapour phase diagram from the integral equations method in the case of interaction potentials representing not only simple liquids, but also complex ones through the coarse-grain concept. We chose a new formulation of the integral equation of Sarkisov characterised by high thermodynamic self-consistency without adjustable parameter. The set of integral equations is solved numerically according to Labik et al. algorithm, while the derivatives of the correlation functions are computed in a formally exact way by a tangent linear method. Finally, the phase diagram is scanned resorting to a self-adapting method producing accurate determinations of both the spinodal and binodal lines as well as of the critical point. In a first stage, we built up and tested our approach in the cases of hard-spheres and Lennard-Jones fluids. Then, we studied structural and thermodynamic properties, as well as phase diagrams of two families of fluids: the first one interacting with hard-core Yukawa-type potentials and the second one with discrete potentials. In each case, our results were compared with simulation predictions available in the literature. We were also able to study the evolution of the phase diagram with the features of the potential, and we observed a liquid-liquid phase transition in a specific discrete potential fluid. Finally, we also studied some metallic systems whose effective potentials depend on the density.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2005-2006 / Un modèle numérique tridimensionnel a été développé pour la simulation du système chimique quartz-eau couplé avec l’écoulement à densité et viscosité variable dans les milieux poreux discrètement fracturés. Le nouveau modèle simule aussi le transfert de chaleur dans les milieux poreux fracturés en supposant que l’expansion thermique du milieu est négligeable. Les propriétés du fluide, densité et viscosité, ainsi que les constantes chimiques (constant de taux de dissolution, constant d’équilibre, coefficient d’activité) sont calculées en fonction de la concentration des ions majeurs et de la température. Des paramètres de réaction et d’écoulement, comme la surface spécifique du minéral et la perméabilité sont mis jour à la fin de chaque pas de temps avec des taux de réaction explicitement calculés. Le modèle suppose que des changements de la porosite et des ouvertures de fractures n’ont pas d’impact sur l’emmagasinement spécifique. Des pas de temps adaptatifs sont utilisés pour accélérer et ralentir la simulation afin d’empêcher des résultats non physiques. Les nouveaux incréments de temps dépendent des changements maximum de la porosité et/ou de l’ouverture de fracture. Des taux de réaction au niveau temporel L+1 (schéma de pondération temporelle implicite) sont utilisés pour renouveler tous les paramètres du modèle afin de garantir la stabilité numérique. Le modèle a été vérifié avec des problèmes analytiques, numériques et physiques de l’écoulement à densité variable, transport réactif et transfert de chaleur dans les milieux poreux fracturés. La complexité de la formulation du modèle permet d’étudier des réactions chimiques et l’écoulement à densité variable d’une façon plus réaliste qu’auparavant possible. En premier lieu, cette étude adresse le phénomène de l’écoulement et du transport à densité variable dans les milieux poreux fracturés avec une seule fracture à inclinaison arbitraire. Une formulation mathématique générale du terme de flottabilité est dérivée qui tient compte de l’écoulement et du transport à densité variable dans des fractures de toute orientation. Des simulations de l’écoulement et du transport à densité variable dans une seule fracture implanté dans une matrice poreuse ont été effectuées. Les simulations montrent que l’écoulement à densité variable dans une fracture cause la convection dans la matrice poreuse et que la fracture à perméabilité élevée agit comme barrière pour la convection. Le nouveau modèle a été appliqué afin de simuler des exemples, comme le mouvement horizontal d’un panache de fluide chaud dans un milieu fracturé chimiquement réactif. Le transport thermohalin (double-diffusif) influence non seulement l’écoulement à densité variable mais aussi les réactions chimiques. L’écoulement à convection libre dépend du contraste de densité entre le fluide (panache chaud ou de l’eau salée froide) et le fluide de référence. Dans l’exemple, des contrastes de densité sont généralement faibles et des fractures n’agissent pas comme des chemins préférés mais contribuent à la dispersion transverse du panache. Des zones chaudes correspondent aux régions de dissolution de quartz tandis que dans les zones froides, la silice mobile précipite. La concentration de silice est inversement proportionnelle à la salinité dans les régions à salinité élevée et directement proportionnelle à la température dans les régions à salinité faible. Le système est le plus sensible aux inexactitudes de température. Ceci est parce que la température influence non seulement la cinétique de dissolution (équation d’Arrhenius), mais aussi la solubilité de quartz. / A three-dimensional numerical model is developed that couples the quartz-water chemical system with variable-density, variable-viscosity flow in fractured porous media. The new model also solves for heat transfer in fractured porous media, under the assumption of negligible thermal expansion of the rock. The fluid properties density and viscosity as well as chemistry constants (dissolution rate constant, equilibrium constant and activity coefficient) are calculated as a function of the concentrations of major ions and of temperature. Reaction and flow parameters, such as mineral surface area and permeability, are updated at the end of each time step with explicitly calculated reaction rates. The impact of porosity and aperture changes on specific storage is neglected. Adaptive time stepping is used to accelerate and slow down the simulation process in order to prevent physically unrealistic results. New time increments depend on maximum changes in matrix porosity and/or fracture aperture. Reaction rates at time level L+1 (implicit time weighting scheme) are used to renew all model parameters to ensure numerical stability. The model is verified against existing analytical, numerical and physical benchmark problems of variable-density flow, reactive solute transport and heat transfer in fractured porous media. The complexity of the model formulation allows chemical reactions and variable-density flow to be studied in a more realistic way than previously possible. The present study first addresses the phenomenon of variable-density flow and transport in fractured porous media, where a single fracture of an arbitrary incline can occur. A general mathematical formulation of the body force vector is derived, which accounts for variable-density flow and transport in fractures of any orientation. Simulations of variable-density flow and solute transport are conducted for a single fracture, embedded in a porous matrix. The simulations show that density-driven flow in the fracture causes convective flow within the porous matrix and that the highpermeability fracture acts as a barrier for convection. The new model was applied to simulate illustrative examples, such as the horizontal movement of a hot plume in a chemically reactive fractured medium. Thermohaline (double-diffusive) transport impacts both buoyancy-driven flow and chemical reactions. Free convective flow depends on the density contrast between the fluid (hot brine or cool saltwater) and the reference fluid. In the example, density contrasts are generally small and fractures do not act like preferential pathways but contribute to transverse dispersion of the plume. Hot zones correspond to areas of quartz dissolution while in cooler zones, precipitation of imported silica prevails. The silica concentration is inversely proportional to salinity in high-salinity regions and directly proportional to temperature in low-salinity regions. The system is the most sensitive to temperature inaccuracy. This is because temperature impacts both the dissolution kinetics (Arrhenius equation) and the quartz solubility.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2015-2016 / Ce mémoire investigue le phénomène de reconnexion visqueuse de deux tourbillons initialement placés de façon orthogonale ou antiparallèle. Les équations de Navier-Stokes pour un fluide incompressible sont résolues directement (DNS) à l’aide d’un algorithme pseudospectral utilisant des expansions périodiques dans les trois directions cartésiennes. La condition de circulation nulle inhérente à cette méthode numérique est contournée en résolvant les équations dans un repère tournant approprié. Une méthode simple utilisant des lignes de vorticité est proposée afin de calculer le pourcentage de reconnexion instantané η de deux tourbillons. Cette méthode est également employée pour séparer le champ de vorticité en ses composantes reconnectées et non-reconnectées, ce qui facilite l’identification visuelle des différentes étapes du processus de reconnexion. Finalement, l’échelle de temps de la reconnexion tourbillonnaire Trec est calculée pour différents nombres de Reynolds (500 ≤ Re ≤ 10000). Il est trouvé que l’ordre de grandeur de Trec varie de façon continue de [symbol] à mesure que Re augmente. / This work focuses on the viscous reconnection phenomenon of two vortex tubes that are initially antiparallel or orthogonal to each other. The incompressible Navier-stokes equations are solved directly (DNS) using a Fourier pseudospectral algorithm with triply periodic boundary conditions. The associated zero-circulation constraint is circumvented by solving the governing equations in a proper rotating frame of reference. A simple method using vortex lines is proposed to compute the instantaneous reconnection level η of two vortices. The proposed method is also used to split the vorticity field into its reconnected and non-reconnected parts, which allows for a clear and intuitive visual identification of the different reconnection phases. Finally, the Reynolds number dependence of the reconnection timescale Trec is investigated for 500 ≤ Re ≤ 10000. The scaling is found to vary continuously as Re is increased from [symbol].

Dans ce mémoire, une approche de calcul de sensibilité utilisée dans le domaine de l'optimisation, appelée la méthode adjointe, est explorée pour estimer et contrôler efficacement l'erreur en mécanique des fluides numérique (MFN). Les développements effectués ont été greffés à Fluent, un logiciel commercial. L'erreur est calculée en effectuant le produit scalaire du vecteur d'erreur de troncature liée à chaque équation de bilan et du vecteur des variables adjointes. Ces dernières sont intrinsèquement liées à une fonctionnelle sélectionnée par l'utilisateur (le rendement, la traînée, les pertes, etc.). En adaptant le maillage avec cette estimation d'erreur, l'adaptation devient très orientée vers le but spécifique de la simulation. Les résultats obtenus montrent que la méthode adjointe permet d'identifier les régions sensibles de la simulation. En adaptant la taille, l'étirement et l'orientation des cellules du maillage en conséquence, il est alors possible d'améliorer la précision numérique ou encore de limiter la taille de maillage requise pour une précision donnée. Quoique se limitant aux cas 2-D stationnaires, les résultats obtenus laissent entrevoir que les maillages adaptés représentent une amélioration substantielle comparativement aux maillages structurés utilisés actuellement. À terme, l'intégration des techniques d'adaptation dans les processus de simulation pourrait fortement influencer les façons de faire de l'Institut de Recherche d'Hydro-Québec (IREQ) ainsi que du Laboratoire de Machines Hydrauliques (LAMH) de l'Université Laval.

Les systèmes flottants de production, stockage et de déchargement (FPSO) ont été introduits dans les secteurs d'exploitation des hydrocarbures offshore en tant qu'outils facilement déplaçables pour l’exploitation de champs de pétrole et de gaz de petites ‘a moyenne tailles ou lorsque ceux-ci sont éloignés des côtes ou en eaux profondes. Ces systèmes sont de plus en plus envisagés pour les opérations de traitement et de raffinage des hydrocarbures à proximité des sites d'extraction des réservoirs sous-marins en utilisant des laveurs et des réacteurs à lit fixe embarqués. De nombreuses études dans la littérature pour découvrir l'hydrodynamique de l'écoulement polyphasiques dans des lits garnis ont révélé que la maîtrise de tels réacteurs continue d’être un défi quant à leur conception /mise à l'échelle ou à leur fonctionnement. De plus, lorsque de tels réacteurs sont soumis à des conditions fluctuantes propres au contexte marin, l'interaction des phases devient encore plus complexe, ce qui entraîne encore plus de défis dans leur conception. Les travaux de recherche proposés visent à fournir des informations cruciales sur les performances des réacteurs à lit fixes à deux phases dans le cadre d'applications industrielles flottantes. Pour atteindre cet objectif, un simulateur de mouvement de navire de type hexapode avec des mouvements à six degrés de liberté a été utilisé pour simuler les mouvements des FPSO tandis que des capteurs à maillage capacitif (WMS) et un tomographe à capacitance électrique (ECT) couplés avec le lit garni ont permis de suivre en ligne les caractéristiques dynamiques locales des écoulements diphasiques. L'effet des inclinaisons et des oscillations de la colonne sur le comportement hydrodynamique des lits garnis biphasiques a été étudié, puis les résultats ont été comparés à leurs analogues terrestres correspondants (colonne verticale immobile). De plus, des stratégies opérationnelles potentielles ont été proposées pour atténuer la maldistribution des fluides résultant des oscillations du lit ainsi que pour intensifier le processus de réactions dans les réacteurs à lit fixe. Parallèlement aux études expérimentales, un modèle Eulérien CFD transitoire 3D a été développé pour simuler le comportement hydrodynamique de lits garnis polyphasiques sous des inclinaisons et des oscillations de colonnes. Enfin, pour compléter le travail expérimental, une étude systématique a été réalisée pour étudier les performances de capture de CO2 à base d'amines d’un laveur à garnissage (en vrac et structuré) émulant une colonne à bord des ... / Floating production storage and offloading (FPSO) systems have been introduced to offshore hydrocarbon exploitation sectors as readily movable tools for development of small or remote oil and gas fields in deeper water. These systems are increasingly contemplated for onboard treatment and refining operations of hydrocarbons extracted from undersea reservoirs near extraction sites using embarked packed-bed scrubbers and reactors. Numerous efforts in the literature to uncover the hydrodynamics of multiphase flow in packed beds have disclosed that such reactors continue to challenge us either in their design/scale-up or their operation. Furthermore, when such reactors are subjected to marine conditions, the interaction of phases becomes even more complex, resulting in further challenges for design and scale-up. The proposed research aims at providing important insights into the performance of two-phase flow packed-bed reactors in the context of floating industrial applications. To achieve this aim, a hexapod ship motion simulator with six-degree-of-freedom motions was employed to emulate FPSO movements while capacitance wire mesh sensors (WMS) and electrical capacitance tomography (ECT) coupled with the packed bed scrutinized on-line and locally the two-phase flow dynamic features. The effect of column tilts and oscillations on the hydrodynamic behavior of multiphase packed beds was investigated and then the results were compared with their corresponding onshore analogs. Moreover, potential operational strategies were proposed to diminish fluid maldistribution resulting from bed oscillations as well as for process intensification of heterogeneous catalytic reactions in packed-bed reactors. In parallel with the experiment studies, a 3D transient Eulerian CFD model was developed to simulate the hydrodynamic behavior of multiphase packed beds under column tilts and oscillations. Ultimately, a systematic experimental study was performed to address the amine-based CO2 capture performance of packed-bed scrubbers on board offshore floating vessels/platforms. Apart from gaining a comprehensive knowledge on the influence of translational and rotational movements on multiphase flows in porous media, oil and gas sectors and ship industry would benefit from the results of this work for design and scale-up of industrial reactors and scrubbers. / Unité flottante de production, de stockage et de déchargement

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.

À l’heure actuelle, les écoulements à deux fluides non miscibles jouent un rôle très important dans plusieurs domaines, que ça soit en science ou en ingénierie. Leur complexité est tellement élevée que les modèles actuels ne permettent de résoudre que des cas particuliers ou simplifiés avec un degré de précision qui demeurent souvent plutôt modeste. Une nouvelle approche numérique parait être une nécessité pour capturer la complexité physique du phénomène. Pour ce faire nous avons besoin d’outils robustes. Au niveau de l’interface de séparation entre les deux fluides non miscibles, les variables physiques sont discontinues, ce qui pose un défi majeur dans la description des variables et des conditions aux limites à l’interface. Le fait que les densités et les viscosités de chaque fluide soient différentes de part et d’autre de l’interface donne naissance à des défauts et des impuretés dans le champ des vitesses, ce qu’on appelle une discontinuité faible. Pour sa part, l’existence de la force de tension superficielle au niveau de l’interface crée une discontinuité sur le champ de pression, ce qu’on appelle une discontinuité forte. Un autre grand problème se pose au niveau de l’étude numérique du problème, où les méthodes numériques classiques ont une précision assez limitée dans ce genre de situation. L’objectif de ce travail est de fournir une étude complète de la dynamique de l’interface entre deux fluides non miscibles à l’aide d’outils mathématiques, physiques et numériques robustes. D’abord, une étude analytique du problème a été faite où l’équation de Navier-Stokes et les conditions de saut sur les variables physiques au niveau de l’interface de séparation entre les fluides ont été prouvées en détail. Pour traiter les discontinuités, nous avons discrétisé nos variables à l’aide de la méthode XFEM. Dû aux larges distorsions rencontrées dans ce genre d’écoulement, nous avons utilisé l’approche Eulérienne, pour corriger les oscillations des solutions dues aux choix du système de coordonnées nous avons utilisé les techniques de stabilisation SUPG/PSPG. Le traitement de la courbure des interfaces K été fait à l’aide de l’opérateur Laplace Beltrami et le suivi d’interface à l’aide de la méthode ¨Level-set¨. Pour le traitement des conditions de saut au niveau de l’interface la méthode Nitsche est développée dans différents contextes. Après avoir développé un modèle physique et mathématique dans les premières parties de notre travail, nous avons fait une étude numérique à l’aide de la plateforme de calcul FEniCS, qui est une plateforme de développement en langage C++ avec une interface Python. Un code de calcul a été développé dans le cas des écoulements de deux fluides non miscibles avec les modèles physiques et les outils mathématiques développés dans les sections précédentes. / The two-phase flow problems have an important role in the multitude of domains in science and engineering. Their complexity is so high that the actual models can solve only particular or simplified cases with a certain degree of precision. A new approach is a necessity to understand the evolution of new ideas and the physical complexity in this kind of flow, to contribute to the study of this field. A good study requires solid and robust tools to have performing results and a maximum of efficacy. At the interface of separation between the two immiscible fluids, the physical parameters are discontinuous, which gives us difficulties for the description of the physical variables at the interface and boundary conditions. The fact that the density and the viscosity are discontinuous at the interface creates kinks in the velocity, which represent a weak discontinuity. The existence of the surface tension at the interface create a discontinuity for the pressure field, it represents a strong discontinuity. The main objective of this work is to make a complete study based on strong and robust physical, mathematical and numerical tools. A strong combination, capable of capturing the physical aspect of the interface between the two fluids with a very good precision. Building such a robust, cost effective and accurate numerical model is challenging and requires lots of efforts and a multidisciplinary knowledge in mathematics, physics and computer science. First, an analytical study was made where the one fluid model of the Navier-Stokes equation was proved from Newton’s laws and jump conditions at the interface was proved and detailed analytically. To treat the problem of discontinuity, we used the XFEM method to discretize our discontinuous variables. Due to the large distortion encountered in this kind of fluid mechanic problems, we are going to use the Eulerian approach, and to correct the oscillation of solutions we will use the SUPG/PSPG stabilization technic. The treatment of the interface curvature k was done with the Laplace Beltrami operator and the interface tracking with the Level-set method. To treat the jump conditions with a very sharp precision we used the Nitsche’s method, developed in different cases. After building a strong mathematical and physical model in the first parts of our work, we did a numerical study using the FEniCS computational platform, which is a platform of computational development based on C++ with a Python interface. A numerical code was developed in this study, in the case of two-phase flow problem, based on the previous mathematical and physical models detailed in previous sections.

Les aspirateurs de turbines hydrauliques jouent un rôle crucial dans l’extraction de l’énergie disponible. Dans ce projet, les écoulements dans l’aspirateur d’une turbine de basse chute ont été simulés à l'aide de différents modèles de turbulence dont le modèle DDES, un hybride LES/RANS, qui permet de résoudre une partie du spectre turbulent. Déterminer des conditions aux limites pour ce modèle à l’entrée de l’aspirateur est un défi. Des profils d’entrée 1D axisymétriques et 2D instationnaires tenant compte des sillages et vortex induits par les aubes de la roue ont notamment été testés. Une fluctuation artificielle a également été imposée, afin d’imiter la turbulence qui existe juste après la roue. Les simulations ont été effectuées pour deux configurations d’aspirateur du projet BulbT. Pour la deuxième, plusieurs comparaisons avec des données expérimentales ont été faites pour deux conditions d'opération, à charge partielle et dans la zone de baisse rapide du rendement après le point de meilleur rendement. Cela a permis d’évaluer l'efficacité et les lacunes de la modélisation turbulente et des conditions limites à travers leurs effets sur les quantités globales et locales. Les résultats ont montrés que les structures tourbillonnaires et sillages sortant de la roue sont adéquatement résolus par les simulations DDES de l’aspirateur, en appliquant les profils instationnaires bidimensionnels et un schéma de faible dissipation pour le terme convectif. En outre, les effets de la turbulence artificielle à l'entrée de l’aspirateur ont été explorés à l'aide de l’estimation de l’intermittence du décollement, de corrélations en deux points, du spectre d'énergie et du concept de structures cohérentes lagrangiennes. Ces analyses ont montré que les détails de la dynamique de l'écoulement et de la séparation sont modifiés, ainsi que les patrons des lignes de transport à divers endroits de l’aspirateur. Cependant, les quantités globales comme le coefficient de récupération de l’aspirateur ne sont pas influencées par ces spécificités locales. / Draft tubes play a crucial role in elevating the available energy extraction of hydroturbines. In this project, turbulent flows in the draft tube of a low-head bulb turbine were simulated using, among others, an advance hybrid LES/RANS turbulent model, called DDES, which can resolve portions of the turbulent spectrum. Providing appropriate inflow boundary conditions for such models is a challenging issue. In this regard, different inflow boundary conditions were tested, including axisymmetric 1D profiles, and unsteady 2D inflow profiles that take runner blade wakes and vortices into account. Artificial fluctuation at the inlet section of the draft tube was also included to mimic the turbulence existing after the runner. Simulations were conducted for two draft tube configurations of the BulbT project. For one of them, intensive comparisons with experimental data were done for two operating conditions, one at part load and another in the sharp drop-off portion of the efficiency hill after the best efficiency point. This allowed to assess the effectiveness and shortcomings of the adopted turbulence modeling and boundary conditions through their effects on the global and local quantities. The results showed that the runner-related vortical structures and wakes are appropriately resolved using stand-alone DDES simulation of the draft tube flows. This is achieved by applying unsteady 2D inflow profiles along with adopting low dissipation scheme for the convective term. Furthermore, the effects of applying artificial turbulence at inlet were explored using separation intermittency, two-point correlation, energy spectrum and Lagrangian coherent structure concepts. These analyses revealed that the type of inflow boundary conditions modifies the details of the flow and separation dynamics as well as patterns of the transport barriers in different regions of the draft tube. However, the global quantities such as recovery coefficient are not influenced by these local features.

La biomasse est une des sources importantes d'énergie primaire et renouvelable. Le développement d'un procédé basé sur la conversion de celle-ci en énergie tout en demeurant respectueux de l'environnement, fait l'objet de recherches intenses aussi bien dans les mondes académique qu'industriel. La gazéification pour produire un gaz de biosynthèse est considérée comme une des options les plus prometteuses via la valorisation des sources de résidus de biomasse. La thermodynamique et la cinétique intrinsèque imposent que les réactions de gazéification de la biomasse doivent être effectuées à des températures élevées, exigeant la fourniture et la récupération de chaleur de manière efficace. Le concept de gazéification allotherme (par opposition à son pendant autotherme) offre une solution attrayante pour la mise en oeuvre à haute température du couplage de réactions fortement endothermique avec des réactions exothermiques. Toutefois, la mise en oeuvre pratique du concept sous haute température n'est pas aisée. Dans ce travail, un nouveau concept pour la gazéification de résidus de la biomasse est proposé impliquant l'hybridation de réactions à hautes températures de la gazéification et de la combustion dans un réacteur monolithique structuré. Clairement, le design et l'optimisation de ce nouveau procédé hybride requiert la compréhension précise, non seulement des phénomènes physico-chimiques de la conversion thermochimique de la biomasse, mais aussi du comportement hydrodynamique, complexe, des deux phases mises en oeuvre dans un microréacteur monolithique à lit fluidisé. À cet égard, la caractéristique hydrodynamique de la distribution des écoulements des phases gaz-solide au sein du réacteur revêt une importance cruciale pour la prédiction du comportement des processus de gazéification/combustion et pour l'examen de stratégies d'opération du procédé. En particulier, en raison de la nature complexe de l'interaction entre le gaz et les particules solides ainsi que la phase stationnaire représentée par le microréacteur monolithique, un des défis dans le design et l'opération de ces réacteurs est la prévention de la maldistribution des phases. Dans ce travail, la mécanique des fluides numériques (MFN) est mise à profit comme outil de simulation permettant d'explorer les distributions des écoulements gaz-solide dans un réacteur monolithique. L'ensemble des sections structurée 111 (le monolithe) et les parties terminales non-structurées (lits fixes aléatoires permettant l'alimentation et l'évacuation de la suspension gaz-solide) est globalement considéré dans la simulation afin de capturer les tendances lourdes des mécanismes contribuant à la dynamique gaz-solide. Les résultats des simulations ont démontré la capacité de la MFN à capturer la caractéristique de non-uniformité de l'écoulement dans ce type de géométrie

No description available.