Tableau d'honneur de la Faculté des études supérieures et postdorales, 2014-2015 / Cette thèse présente une étude de la turbulence et du scalaire passif dans la zone autosimilaire d’un jet rond libre faiblement chauffé, caractérisé par un haut nombre de Reynolds ([symbol]). D’abord, une analyse des limitations des sondes à fils chauds croisés stationnaires basée sur des mesures LDV est présentée. Celle-ci met en évidence l’importance des événements de vitesse à grands angles d’incidence dans les valeurs mesurées. Ensuite, la turbulence et la chaleur employée comme scalaire passif dans le jet sont étudiées à l’aide du LDV, de la thermométrie à fil froid et d’une combinaison de ces deux méthodes. On mesure ainsi les distributions radiales des valeurs moyennes et des moments de vitesse, de température et de leurs corrélations croisées. Ces statistiques incluent notamment [symbol]. Ces mesures sont caractérisées par un contrôle précis des conditions aux limites et une grande convergence statistique. On remarque un bon accord entres les valeurs mesurées de [symbol] , [symbol] ainsi que [symbol] et leur expressions analytiques, dérivées des équations de continuité, du mouvement moyen et d’enthalpie, respectivement. Les bilans de transports de l’énergie cinétique turbulente (k), du scalaire passif ([symbol] ) ainsi que des flux de chaleur turbulents [symbol] sont ensuite établis, dans le même écoulement de jet rond faiblement chauffé. Pour le bilan de k, le modèle proposé par Lumley (1978) est employé pour représenter les termes de diffusion de pression. Pour les deux bilans, l’intégration à zéro sur un volume de contrôle des termes de diffusion est vérifiée et la portée de ce critère de validité est analysée. Aussi, les restes des deux bilans constituent leurs termes de dissipation respectifs, ek et eq . Dans le cas de ek, le modèle prédictif mis de l’avant par Thiesset et al (2014) ainsi qu’une analyse originale des bases de données de Panchapakesan & Lumley (1993a) et de Hussein et al (1994) laisse entrevoir une possibilité de réconciliation en ce qui a trait aux résultats produits par ces deux expériences réputées. Les budgets des composantes du tenseur de Reynolds ([symbol] ) sont aussi établis et analysés. La valeur de eq sur l’axe du jet est employée comme référence dans une étude exhaustive des techniques de mesures directes de ce taux de dissipation. Cette étude expose clairement les difficultés de ces mesures, notamment par des analyses spectrales. Les données amassées permettent d’établir les distributions radiales de paramètres d’intérêt pour la modélisation, tels que la viscosité tourbillonnaire ([symbol] ), la diffusivité tourbillonnaire ([symbol] ) et le nombre de Prandtl turbulent ([symbol] ). L’assomption de [symbol] constant est examinée. Le ratio entre les échelles de temps de dissipation de température et d’énergie cinétique turbulente (R) est aussi calculé et analysé. / This thesis presents a study of the turbulence generated by a slightly heated turbulent round jet, evolving in still air. The study focuses on the self-preserving region of the high Reynolds number jet ([symbol] ). The jet is heated to only 20°C above ambient, so that temperature is treated as a passive scalar. In the first chapter of the thesis, a new and novel analysis of the stationary crossed hot wire probes based on unbiased LDV data is exposed, thus revealing the importance of the wide incidence angle events on the measured values. LDV, cold wire thermometry and a combination of both methods is used to study the slightly heated jet. Radial distributions of mean velocities and heating level, second and third order moments of fluctuating velocity and temperature ([symbol]) along with cross correlations, such as [symbol] are made available. During data acquisition, great attention was devoted to the close control of boundary conditions and statistical convergence. The measured values of V, uv as well as vq are in good agreement with their analytical expressions, provided by continuity, mean momentum and enthalpy equations, respectively. The transport budgets for the turbulent kinetic energy, ([symbol] ), passive scalar ([symbol] ), and temperature turbulent fluxes ([symbol] ) are then established using the data acquired in the very same, slightly heated, jet flow. The model introduced by Lumley (1978) for the pressure diffusion, is used as a surrogate for this term, in the case of the k budget. The integration to zero of the diffusion terms of both budgets over a control volume is verified. The significance of this validity criterion is analysed. The remainders of both budgets, are their respective dissipation terms, ek and eq . In the case of ek, the predictive model recently put forward by Thiesset et al (2014) along with a novel treatment of the databases of Panchapakesan & Lumley (1993a) and Hussein et al (1994), may help bring the outcome of the two famous experiments to a reconciliation. The Reynolds stress tensor component budgets were also established, leading to an analysis of the importance of the redistribution terms. The [symbol] value provided by the [symbol] budget is used as a reference in an exhaustive study of the means of dissipation measurement, on the centerline of the jet. This study reveals the shortcoming of the direct measurements of eq , namely thought spectral analysis. Radial distributions of turbulence modeling parameters, such as turbulent momentum viscosity ([symbol] ), turbulent thermal diffusivity ([symbol] ) and turbulent Prandtl number ([symbol] ) were established. The assumption of a constant PrT is scrutinized. The radial distribution of the thermal-to-mechanical timescale ratio (R) is calculated and analysed.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2017-2018 / Ce mémoire vise à contribuer à l'amélioration de la prédiction numérique des performances d'une turbine hydraulique par la quantification de l'impact de certains paramètres d'influence, autant physiques que méthodologiques. Pour ce faire, l'écoulement dans une turbine caractérisée par une chute de rendement près du meilleur point de fonctionnement a été étudié. Cette dernière est associée à une baisse abrupte du coefficient de récupération de l'aspirateur, une composante connue pour jouer un rôle capital dans les performances globales de la machine. L'accent a particulièrement été mis sur l'étude du débalancement de l'écoulement sous la roue, observé à la suite d'une campagne de mesure expérimentale effectuée sur le banc d'essai du CREMHyG à Grenoble. En effet, l'écoulement en sortie de roue est souvent considéré axisymétrique dans les simulations numériques, une simplification qui est en désaccord avec les données expérimentales du cas étudié dans ce projet de recherche. Afin d'identifier la nature du débalancement, une analyse critique a été menée sur les données expérimentales. Toutefois, aucune conclusion satisfaisante n'a pu être tirée puisqu'une inconsistance entre les données obtenues par LDV et celles obtenues par PIV a été décelée. Considéré statique, il a été démontré que la présence d'un tel débalancement à l'entrée de l'aspirateur altère considérablement la topologie de l'écoulement en aval ainsi que la prédiction des performances. De plus, certains aspects des simulations obtenus avec ce type de débalancement ne concordent pas avec les observations expérimentales, c'est notamment le cas pour la zone d'amorçage du décollement mais aussi pour l'écart de débit entre les deux pertuis. Il a également été montré que les simulations des composantes fixes en amont de l'aspirateur, incluant une importante partie de la conduite d'amenée, ne permettent pas de retrouver l'asymétrie de l'écoulement observée expérimentalement en sortie de roue. En effet, l'écoulement secondaire, associé à la présence d'une conduite coudée en amont de la bâche spirale, s'est avéré avoir très peu d'influence sur la distribution de l'écoulement dans l'avant-distributeur et dans le distributeur. La conduite menant jusqu'à l'entrée de la bâche spirale semble être assez longue pour que les modifications apportées au profil de vitesse deviennent insuffisantes pour perturber l'écoulement à cet endroit. / This master degree thesis aims at enhancing the numerical prediction of the performance of a hydraulic turbine by quantifying the impact of physical and methodological influencing parameters. To this end, the flow in a turbine characterised by an efficiency drop near the best efficiency point was studied. The latter is associated with a sharp drop in the recovery coefficient of the draft tube, a component known to play a key role in the global performance of the machine. Particular emphasis was placed on the study of the imbalance of the flow at the runner exit, observed following an experimental measurement campaign carried out on the CREMHyG test bench in Grenoble. Indeed, the flow at the runner exit is often considered axisymetric in numerical simulations, a simplification that is not consistent with the experimental data of the case studied in this research project. In order to identify the nature of the imbalance, a critical analysis was conducted on the experimental data. However, no satisfactory conclusion could be drawn since inconsistency between the data obtained by LDV and those obtained by PIV was identified. Considered static, it has been shown that the presence of such an imbalance at the inlet of the draft tube alters considerably the topology of the downstream flow as well as the prediction of the performances. Moreover, some aspects of the simulations obtained with this type of imbalance do not agree with the experimental observations, it is notably the case for the initiation zone of the flow detachment but also for the flow rate imbalance between the two draft tube channels. It has also been shown that the numerical simulations of fixed components upstream of the draft tube, including a large part of the penstock, do not allow to recover the asymmetry of the flow experimentally observed at the exit of the runner. Indeed, the secondary flow, associated with the presence of a bend upstream of the spiral case, proved to have a very little influence on the distribution of the flow in the distributor. The pipe leading to the entrance of the spiral case appears to be long enough for changes made in the velocity profile to become insufficient to disrupt the flow at this place.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2006-2007 / L'objectif de ce projet de maîtrise est de contribuer à une meilleure compréhension des écoulements dans les conduits rectangulaires en rotation. Une attention particulière est apportée aux effets de bout et au développement de l'écoulement. Des simulations RANS en régime stationnaire ont démontré que la présence d'un écoulement secondaire de faible intensité est suffisante pour modifier le profil de vitesse axiale de façon à que ce dernier diffère significativement de la solution théorique 2D. En conséquence, les comparaisons entre les résultats expérimentaux et les simulations numériques 2D devraient être faites avec prudence. Les calculs stationnaires ont également montré que le développement de l'écoulement dans un canal tournant est affecté par l'apparition de cellules de type Taylor-Görtler le long du mur pression. Puisque ces tourbillons sont légèrement inclinés par rapport au plan horizontal, une solution stationnaire pleinement développée n'est jamais atteinte. Des simulations subséquentes en régime instationnaire dans un domaine périodique ont prouvé qu'une solution pleinement développée peut cependant exister dans un canal tournant, mais l'écoulement est alors formellement instationnaire car les cellules longitudinales sont convectées vers les bouts de la conduite par l'action de l'écoulement secondaire. / The objective of the present work is to provide a basic description and to contribute toward a better understanding of rotating high aspect ratio duct flows. Particular attention is given to end-wall effects and flow development. Steady RANS simulations have shown that the presence of a weak, but unavoidable, end-wall generated secondary flow produces sufficient transverse advection to alter the streamwise velocity profile in such a way that it differs noticeably from the theoretical 2D channel case. Thus, comparisons between experimental results and 2D numerical simulations should be performed very cautiously. Steady state calculations have also revealed that the flow development in a rotating duct is subtly affected by the appearance of longitudinal Taylor-Görtler type cells near the pressure wall. Since the roll cells are tilted toward the end-wall, no steady fully developed solution is ever reached. Subsequent unsteady flow calculation in a streamwise periodic domain have shown that a fully developed solution can exist in a rotating duct, but the resulting flow field is formally unsteady due to the convection of the roll cells toward the end-wall.

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.

Des écoulements instationnaires au niveau du distributeur des turbines à réaction ont été expérimentalement étudiés dans ce projet doctoral. Le but est de développer une expertise intégrée pour mesurer des écoulements instationnaires au niveau des distributeurs des turbines à réaction en utilisant une vélocimétrie par image de particules (PIV) et d'autres techniques de mesure. Les expériences ont été effectuées en deux phases appelées configuration 1 et configuration 2. La configuration 1 est composée d'une directrice, installée dans un canal bidimensionnel d'eau, examinée aux nombres de Reynolds comparable à ceux dans les modèles standards de turbines Francis. La configuration 2 consiste en un distributeur d'un modèle réduit d'une turbine Francis. Dans la configuration 1, les champs d'écoulement ont été mesurés avec un PIV autour d'une directrice asymétrique dans un canal d'eau, à différents angles d'attaque entre 0° et 20°. Les données obtenues comprennent les champs de vitesse instantanée et moyenne. Pour les mêmes conditions, les pertes de charge dues à la directrice ont été mesurées et comparées. Dans la configuration 2, les champs d'écoulement entre deux directrices voisines ont été mesurés avec le PIV au travers d'un accès optique spécifique. Des résultats instantanés tels que les champs de vitesse et les résultats moyens tels que la fluctuation de vitesse et les iso-contours de vitesse ont été présentés. La distribution de pression du côté intrados et du côté extrados de deux directrices voisines a été mesurée avec des capteurs de pression miniatures montés sur les surfaces des directrices. 'Basé sur les mesures, les effets des points d'opération sur le comportement des écoulements sont analysés.

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.

Parmi les réacteurs triphasiques gaz-liquide-solide utilisés dans la pratique industrielle, les réacteurs catalytiques à lit fixe arrosé à cocourant de gaz et de liquide vers le bas, i.e., trickle bed reactors (TBR), sont très répandus en particulier dans divers processus de transformation à hautes température et pression. Les travaux expérimentaux se poursuivent depuis plus de quatre décennies sur la quantification des paramètres hydrodynamiques (transition des régimes d'écoulement, perte de pression biphasique, rétention liquide, efficacité de mouillage, etc.) pour cette configuration de réacteurs. Différentes approches ont été mises en œuvre par un grand nombre d’équipes de recherche pour mesurer ces paramètres hydrodynamiques dans le but de construire des outils de prédiction et de description par rapport aux conditions réelles d’opération des processus à l’échelle industrielle. La présente contribution se propose de répondre à la question suivante : Dans quelle mesure les connaissances accumulées à partir d’observations à l’échelle laboratoire dans les conditions ambiantes sont-elles fiables pour opérer un TBR à pression et température élevées? Une question sous-jacente à la précédente concerne le comportement hydrodynamique avec la température lorsque le réacteur est alimenté par un liquide non-newtonien. L'intensification des procédés est une approche en vogue et prometteuse pour continuer à apporter des perfectionnements (gains en économie et en efficacité) au réacteur TBR. Aussi, l’induction artificielle d’impulsions est-elle envisagée dans cette étude en tant que méthode d'intensification de processus pour des températures et pressions non-ambiantes. Le présent travail tentera de démontrer les avantages de plusieurs variantes de l'opération périodique sur l'hydrodynamique des TBR pour des systèmes coalescent, non-newtonien et moussant à des températures et pressions augmentées. / Trickle bed reactor (TBR) is one of the most widely used three-phase reactors in various processes mostly operated at high temperature and high pressure. The ongoing experimental work on the hydrodynamic parameters (flow regime transition, pressure drop, liquid holdup, wetting efficiency etc.) of this packed bed reactor configuration goes to early 1960’s. Different techniques were applied by different researchers for the measurement of these hydrodynamic parameters which let the comparison and the decision of more convenient method by means of doing investigations at conditions near to that of industrial processes. Process intensification is considered to be a leading approach for the ongoing research on the economic reduction and reactor efficiency enhancement. Artificial induction of pulses is pronounced as one of the methods for the process intensification in TBRs. As trickle bed reactor is also used in biochemical processes, and the initial liquid behaving like a Newtonian fluid could turn into a non-Newtonian fluid after various biochemical processes; it is emphatic to study TBR hydrodynamics with non-Newtonian systems. Despite large amount of work exists in the literature for steady state hydrodynamics of TBR operating at high pressure; the hydrodynamic behavior of TBR at high temperature has been left as a concealed issue. Additionally none of the experimental work performed to demonstrate the advantages of periodic operation on TBR hydrodynamics dealt with the effects of increased temperature and pressure. This study illustrates the hydrodynamics of TBR at increased temperature and pressure under constant throughput flow and cyclic operation.

Ce mémoire présente l’étude expérimentale de l’écoulement d’entrée d’un aspirateur d’une turbine bulbe présentant une chute abrupte de performance. Des mesures par vélocimétrie laser à effet Doppler (LDV) ont été réalisées sur deux axes soit en aval des pales de la roue et en aval du moyeu de la roue. Une particularité de cette étude est la conception d’un montage permettant de mesurer la vitesse axiale proche de la paroi du cône. De plus, une méthode d’estimation de la vitesse radiale moyenne a été développée. Ces mesures ont permis de caractériser l’écoulement primaire et les écoulements secondaires et d’analyser leur évolution entre les deux axes. De plus, l’évolution de ces écoulements est analysée en fonction de la chute de performance de la turbine. Les principales particularités de l’écoulement sont la présence d’une recirculation sous le moyeu, d’une zone contrarotative, des sillages des directrices et des tourbillons de bout de pale. / This thesis presents the experimental study of the draft tube inflow of a bulb turbine with a steep drop in performance. Measurements by laser Doppler velocimetry (LDV) were performed on two axes: downstream of the blades of the wheel and downstream of the wheel hub. A feature of this study is the design of a setup that allows the measurement of the axial velocity near the wall of the cone. In addition, a method of estimating the average radial velocity was developed. These measurements were used to characterize the primary flow and the secondary flows and analyze their evolution between the two axes. Moreover, the evolution of these flows is analyzed according to the turbine performance drop. The main features of the flow are the presence of recirculation downstream of the hub, a counter-rotating area, the guide vane wakes and of the blade tip vortices.

Ce travail de maitrise s’inscrit dans le cadre des activités de recherche du Laboratoire de Machines Hydraulique de l’Université Laval et a pour objectif la caractérisation de l’écoulement dans le canal d’entrée d’un modèle de turbine hydraulique de type bulbe. La représentation des champs de vitesses moyennes et des fluctuations sous différentes conditions d’opérations ont été obtenues en utilisant un système de mesure LDV. Un débalancement du débit et des structures tourbillonnaires dans le canal d’entrée ont été mises en évidence. La conception d’une géométrie d’obstacle provoquant une non-uniformité dans le canal d’entrée a été développée en tenant compte des prédictions de simulations numériques. Des simulations numériques de la machine complète en régime stationnaire et instationnaire selon deux configurations géométriques ont été menées pour déterminer l’influence des conditions d’entrée de l’écoulement sur les performances de la machine. Des comparaisons entre les quantités expérimentales et numériques ont été réalisées. / This work is part of the research activities of the Hydraulic Machines Laboratory of the Laval University and its objective is to contribute to the characterisation of the intake flow in a model of a bulb turbine. The representation of the mean velocity fields and the turbulent fluctuations under predefined operating conditions were obtained by the use of a LDV measurement system. Mass-flow imbalance and vortices in the intake channel were identified. The conception of an obstacle geometry causing a non-uniformity in the intake channel has been developed taking in consideration the predictions of numerical simulation. Numerical simulations of the complete machine for both steady and unsteady case were performed with and without obstacle in the intake channel. The objective of this process was to evaluate the influence of the intake flow condition on the turbine performances. Moreover, comparisons between experimental and numerical quantities were made.

Un colmatage des réacteurs à lit fixe arrosé, i.e., Trickle Bed Reactor (TBR), peut se produire lorsque des suspensions liquide-solide diluées y sont traitées ou lorsque des processus chimiques induisent la formation de solides, par exemple du coke, à l'instar de certaines opérations de raffinage du pétrole. Suite à l'accumulation excessive de ces fines solides, le gradient de pression peut atteindre un niveau limite mettant en jeu la tolérance physique du catalyseur. Ainsi, indépendamment de l'activité chimique résiduelle de catalyseur, l'arrêt du réacteur peut s'imposer pour un remplacement du lit catalytique par un lot de catalyseurs frais avec les conséquences économiques associées. Les activités de recherche dans la littérature ouverte concernant ce domaine ont été limitées à quelques travaux expérimentaux et de modélisation. La nécessité d'une recherche plus approfondies et plus systématique s'est donc imposée à nous afin de déchiffrer les mécanismes sous-jacents de la filtration en mode biphasique et de proposer éventuellement de nouvelles solutions d'opération de ce type de réacteurs. L'objectif de la recherche proposée était de fournir une meilleure compréhension des mécanismes et de la structure des dépôts dans les réacteurs à lit fixe arrosé traitant des suspensions d'hydrocarbures afin de proposer de nouvelles stratégies de décolmatage. En raison de la nature opaque du lit, les observations pariétales, les mesures de perte de pression et de dépôts spécifiques, ont été couplés à des mesures non-invasives basées sur l'imagerie par la tomographic électrique capacitive afin de dévoiler la structure interne du lit lors de la déposition de fines particules. Pour surmonter les forces qui concourent à la déposition sur les collecteurs, de nouvelles stratégies, incluant des configurations différentes de réacteur ou des modes d'alimentation cyclique ont été abordées dans notre étude. Les réacteurs de type monolithe, consistant en des blocs à multiples canaux droits, ont été étudiés avec l'objectif de minimisation du dépôt spécifique et de la perte de pression biphasique. La filtration biphasique de suspensions a été étudiée aussi bien en régime d'écoulement de Taylor (petits canaux) qu'en régime de film (grands canaux). Les forces inter-particulaires fine-fine et fine-collecteur sont déterminantes pour que la déposition de fines puisse avoir lieu. Par conséquent, différentes opérations cycliques ont été examinées pour étudier l'effet positif des pulsations hydrodynamiques liquides pour surmonter les forces d'adhésion entre dépôts et lit. La mise en oeuvre de forces répulsives fortes via la modulation de la conductivité électrique dans des milieux non-aqueux a également été étudiée pour montrer l'efficacité de cette approche pour la prévention de colmatage.