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.

Thèse ou mémoire avec insertion d’articles / Les puits de pétrole et de gaz ont leur propre cycle de vie: le forage, la production, et finalement le bouchage et l'abandon (P&A) des puits. Les opérations de P&A impliquent généralement la mise en place d'un bouchon de ciment dans des sections spécifiques du puits de forage, afin de fournir des isolations zonales et de prévenir toute fuite. L'une des techniques courantes de mise en place d'un bouchon de ciment est la méthode du dump bailing, dans laquelle un coulis de ciment est injecté via un conteneur excentrique dans le puits de forage rempli d'un fluide in situ. Le processus est affecté par divers paramètres d'écoulement, tels que la différence de densité (flottabilité), le taux d'injection, les considérations géométriques, l'angle d'inclinaison, le rapport de viscosité et la contrainte d'écoulement. Motivée par cette application industrielle, cette thèse de doctorat étudie l'écoulement par injection miscible flottante d'un fluide lourd (dense) dans une conduite inclinée à extrémité fermée remplie d'un fluide léger (moins dense). Pour ce faire, des expériences accompagnées de méthodes semi-analytiques et informatiques complémentaires sont utilisées pour analyser les effets des paramètres susmentionnés sur la dynamique de l'écoulement. Pour généraliser nos résultats, nous les présentons en utilisant les nombres adimensionnels pertinents de l'écoulement, y compris le nombre de Reynolds (Re) qui représente les contraintes inertielles par rapport aux contraintes visqueuses, le nombre de Froude (Fr) qui montre le rapport des forces d'inertie aux forces de flottabilité, le nombre d'Archimède (Ar) et le nombre de flottabilité (χ) qui sont tous deux le rapport des contraintes de flottabilité aux contraintes visqueuses, le nombre de Bingham (B) qui représente la limite d'élasticité par rapport aux contraintes visqueuses, le rapport de viscosité (M), l'angle d'inclinaison (β), le rapport de rayon (Rc) qui est le rapport du rayon intérieur du tube intérieur au rayon intérieur du tube extérieur, le rapport d'aspect (δ) qui représente le rapport de la hauteur de dégagement au rayon intérieur du tube extérieur, et l'excentricité (E). Pour les flux d'injection miscibles flottants iso-visqueux, on observe quatre étapes distinctes au cours du développement du flux: le jet initial, la région de mélange, la région d'affaissement et le front de fluide lourd atteignant l'extrémité de la conduite et retournant vers la région de mélange. Les caractéristiques de l'écoulement sont quantifiées pour chaque phase d'écoulement, sur la base de Re, Ar, Fr et β. Cette étude est ensuite étendue aux effets des rapports de viscosité faibles à modérés (0.04 ≤ M ≤1), à divers rapports d'aspect (δ) et à la dynamique de l'écoulement dans la région annulaire entre la conduite d'injection (intérieure) et la conduite extérieure. On constate que le placement de fluide le plus efficace se produit lorsque δ est petit. En outre, la diminution de M améliore l'élimination du fluide in situ, en augmentant la vitesse du front du fluide lourd dans la section circulaire du domaine d'écoulement. Les résultats des flux d'injection de contrastes de viscosité élevés (0.0025 ≤ M ≤ 0.10) montrent qu'une séparation du front du fluide lourd se produit lorsque M est petit, Re est grand, et β ≫ 30, ce qui conduit à classer les flux en régimes de séparation et de non-séparation, dans un plan de groupe adimensionnel basé sur une combinaison de Ar, β, Re, et M. Pour l'injection miscible flottante d'un fluide à forte limite d'élasticité dans un fluide newtonien léger, des expériences complémentaires et une modélisation semi-analytique ont permis d'identifier trois régimes d'écoulement distincts, à savoir les régimes de rupture, d'enroulement et de flambage (bulging). Dans cette étude, un modèle d'approximation de la lubrification est développé en utilisant l'équation constitutive de Herschel-Bulkley. Sur la base d'une prédiction raisonnable du début de la cédulation, les régimes d'écoulement sont classés par rapport à une combinaison élégante des paramètres d'écoulement sans dimension, y compris B, χ, M, et δ. Un bon accord est obtenu entre la carte des régimes d'écoulement fournie par le modèle et les expériences. Le modèle permet d'étudier les variations de la limite de transition du régime en fonction du rapport des rayons (Rc) et de l'excentricité (E). Les trois principaux régimes d'écoulement observés sont étudiés plus en détail à l'aide de techniques d'imagerie par caméra à grande vitesse, de fluorescence induite par laser et de vélocimétrie Doppler à ultrasons. En conséquence, les caractéristiques détaillées de l'écoulement et les sous-régimes à plus long terme sont analysés. En particulier, le régime de rupture est étudié en termes de processus de cédulation, d'encolure et de pincement, de longueur des filaments et de leurs comportements de chute ultérieurs. Le régime d'enroulement est évalué en termes de comportements d'enroulement régulier, libre et irrégulier. Enfin, le régime de flambage (bulging) est analysé en termes d'apparition d'une partie tombante détachée et de profondeur de pénétration maximale. / Oil and gas wells have their own life cycle: drilling, production, and eventually plug and abandonment (P&A) of the wells. P&A operations usually involve the placement of a cement plug in specified sections of the wellbore, to provide zonal isolations and prevent any leakages. One of the common cement plug placement techniques is the dump bailing method, in which cement slurry is injected via an eccentric container into the wellbore filled with an in-situ fluid. The process is affected by various flow parameters, such as the density difference (buoyancy), the injection rate, the geometrical considerations, the inclination angle, the viscosity ratio, and the yield stress. Motivated by this industrial application, in this Ph.D. thesis, buoyant miscible injection flow of a heavy (dense) fluid into an inclined closed-end pipe filled with a light (less dense) fluid is studied. To do so, experiments accompanied by complementary semi-analytical and computational methods are used to analyze the effects of the aforementioned parameters on the flow dynamics. To generalize our results, we present them using the relevant dimensionless numbers of the flow, including the Reynolds number (Re) which represents the inertial stresses to the viscous stresses, the Froude number (Fr) which shows the ratio of the inertial forces to the buoyant forces, the Archimedes number (Ar) and the buoyancy number (χ) which both are the ratio of the buoyant stresses to the viscous stresses, the Bingham number (B) which represents the yield stress to viscous stresses, the viscosity ratio (M), the inclination angle (β), the radius ratio (R[subscript c]) which is the ratio of the inner radius of the inner pipe to the inner radius of the outer pipe, the aspect ratio (δ) which represents the ratio of the releasing height to the inner radius of the outer pipe, and the eccentricity (E). For iso-viscous buoyant miscible injection flows, four distinct flow stages during the flow development are observed: the initial jet, the mixing region, the slumping region, and the heavy fluid front reaching the pipe end and returning toward the mixing region. The flow features are quantified for each flow stage, based on Re, Ar, Fr and β. Then, this study is extended by considering the effects of low-to-moderate viscosity ratios (0.04 ≤ M ≤1), various aspect ratios (δ), and taking into account the flow dynamics in the annular region between the injection (inner) pipe and the outer pipe. It is found that the most efficient fluid placement occurs when δ is small. Also, decreasing M enhances the removal of the in-situ fluid, via increasing the heavy fluid front velocity in the circular section of the flow domain. The results of high viscosity contrasts (0.0025 ≤ M ≤ 0.10) injection flows show that a heavy fuid front separation occurs when M is small, Re is large, and β ≫ 30, leading to classify the flows into separation and non-separation regimes, in a dimensionless group plane based on a combination of Ar, β, Re, and M. For the buoyant miscible injection of a heavy yield stress fluid into a light Newtonian fluid, using complementary experiments and semi-analytical modelling, three distinct flow regimes, namely, the breakup, coiling, and buckling (bulging) regimes are identified. In this study, a lubrication approximation model is developed using the Herschel-Bulkley constitutive equation. Based on a reasonable prediction to the yielding onset, the flow regimes are classified versus an elegant combination of the dimensionless flow parameters, including B, χ, M, and δ. A good agreement is achieved between the flow regime map provided by the model and experiments. Using the model, the regime transition boundary variations versus the radius ratio (R[subscript c]) and the eccentricity (E) are investigated. The three main observed flow regimes are further investigated based on high-speed camera imaging, laser-induced fluorescence, and ultrasound Doppler velocimetry techniques. Accordingly, detailed flow features and sub-regimes at longer times are analyzed. In particular, the breakup regime is studied in terms of the yielding, necking and pinch-off process, the length of filaments, and their subsequent falling behaviours. The coiling regime is evaluated in terms of the regular, free and irregular coiling behaviours. Finally, the buckling (bulging) regime is analyzed in terms of the appearance of a detached falling part and the maximum penetration depth.

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.

Trente-quatre essais d'écoulement provenant de quatre groupes d'auteurs ont été sélectionnés dans la littérature afin d'étudier l'instabilité interne de matériaux de faible plasticité soumis à l'écoulement de l'eau et initiée par suffosion. Les essais sélectionnés ont tous été effectués dans un perméamètre en appliquant les contraintes hydrodynamiques de manière ascendante et par palier. L'objectif était de déterminer la vitesse critique (vcr(l))p qui initie la mobilisation des particules fines dans les pores de la structure prédominante des matériaux, soit celle qui forme les constrictions. La méthode d'analyse présentée par Kovacs et Ujfaludi (1983), mettant en relation la résistance à l'écoulement ([Delta]h/Q) en fonction de la charge totale appliquée ([Delta]h), a été utilisée afin de déterminer les vitesses critiques (vcr(¹))p. Les essais ont tous été réinterprétés avec cette méthode et celle-ci a permis d'extraire des conditions critiques à des stades inférieurs durant les essais comparativement aux méthodes d'analyses utilisées par leurs auteurs respectifs. Ainsi, par cette approche, un nouveau critère de stabilité interne, basé sur celui développé par Kovacs et Ujfaludi (1983) pour l'application sur les puits, a été proposé pour l'étude des barrages en remblai. Puis, toujours suite à l'analyse des essais avec la méthode [Delta]h/Q - [Delta]h, les matériaux utilisés lors des essais ont été classés en trois catégories selon leur type de comportement sous écoulement après la vitesse critique (vcr(¹))p dans le repère [Delta]h/Q - [Delta]h. Il a ensuite été discuté que la sélection de matériaux stables pour l'étude des barrages doit tenir compte du type de comportement sous écoulement des matériaux et doit inclure un critère admissible de perte de fines, tel que présenté par l'approche russe. Enfin, via l'étude de certains paramètres qui influencent la vitesse critique, il a été démontré que la seule façon d'obtenir une tendance dans les résultats était de représenter dans un repère cartésien logarithmique les vitesses d'écoulement critiques (vcr(¹))p par rapport à la conductivité hydraulique (K) des matériaux.

Ce mémoire présente les résultats de la première étape d'une étude sur un jet circulaire turbulent faiblement chauffé et confiné. Ces travaux, réalisés d'août 2004 à mai 2006, ont mené au développement d'une nouvelle approche expérimentale pour l'étude des jets circulaires confinés. L'étude a été réalisée dans trois enceintes de surfaces différentes, permettant ainsi d'étudier le comportement d'un jet faiblement confiné (pour (An)05 = 42 D) à celui d'un jet moyennement confiné (pour (AR)0,5 = 30Det20D). Les profils de la vitesse et de la température moyennes du jet ont été mesurés et comparés avec des résultats expérimentaux issus de la littérature. Ces essais ont permis d'établir que le jet est autosimilaire à une distance d'environ x/D = 20 en aval de la sortie. La plupart des mesures ont été effectuées dans la région auto-similaire, c'est-à-dire à des distances supérieures à x/D = 20. L'objectif était de vérifier l'effet du nombre de Reynolds sur la diffusion du jet dans la région où la turbulence est pleinement développée. Le mémoire présente aussi une étude expérimentale dans la zone la plus proche de la sortie du jet dénommée "cône potentiel", ainsi que la mesure de l'ordre de grandeur des échelles de Kolmogorov dans la région auto-similaire du jet.

Thèse ou mémoire avec insertion d'articles / Les gisements aurifères orogéniques, lesquels représentent la principale source d'or au Canada, sont relativement bien documentés, mais la source des fluides hydrothermaux aurifères dans ces gisements reste équivoque. Le district aurifère de Meliadine (MGD), situé dans la province du Churchill occidental au Nunavut, renferme approximativement 9 Moz (~280 tonnes) d'or (production historique, réserves prouvée et probables, et ressources mesurées, indiquées, et inférées), sous la forme d'une série de gisements situés dans la ceinture de roches vertes de Rankin Inlet (Tiriginaq, Wesmeg, Normeg et Pump). Ces gisements sont spatialement associés à la faille Pyke et ses subsidiaires, formant des corridors de veines de quartz-carbonate subparallèles à travers des unités géologiques d'âge archéen. Ces formations incluent des formations de fer, des roches métasédimentaires clastiques et des roches métavolcaniques. Les zones minéralisées sont de style veines et de style remplacement des roches hôtes en bordure des veines. Les zones de remplacement sont encaissées dans les lithologies riches en fer sous la forme d'or associé à de multiples générations de sulfures. La présente étude documente la signature multi-isotopique du soufre (δ$^{34}$SS, Δ$^{33}$SS et Δ$^{36}$S) des sulfures associés à l'or dans les veines d'or orogéniques et les zones de remplacement associées, afin de caractériser la source et tracer l'évolution des fluides aurifères et de mieux comprendre les contrôles sur la précipitation de l'or. La séquence paragénésique établie décrit la génération 1 de sulfures comme correspondant à la pyrrhotite±(pyrite-chalcopyrite) précoce qui n'est pas associée à l'or, alors que la génération 2 est associée à la sulfuration des unités riches en fer et se compose d'arsénopyrite-pyrrhotite±(pyrite-galène-chalcopyrite) associée à l'or, subdivisée en 2a et 2b. L'aspy2a est riche en inclusions, alors que la génération 2b est pauvre en inclusions et est associée à la déformation et à la recristallisation tardive des sulfures 2a le long de fractures riches en Au libre, Pb, Te et Bi, et aux halos de recristallisation enrichis en As et Se et appauvris en Sb. Des analyses multi-isotopiques de soufre ($^{32}$SS, $^{33}$SS, $^{34}$SS et $^{36}$S) montrent que le Δ$^{33}$SS=0.3 ‰ ± 0.2 et le Δ$^{36}$S=-0.7 ± 0.6 ‰ sont homogènes à travers les générations de sulfures. Le δ$^{34}$SS de la génération 1 varie de 1.4 ‰ à 2.4 ‰, de 1.7 ‰ à 5.5 ‰ pour la génération 2a et de -0.5 ‰ à 5.6 ‰ pour la génération 2b. En se basant sur la concentration de l'or le long de fractures de l'aspy2a, leur abondance d'inclusions, ainsi que sur les enrichissements et appauvrissements systématiques en Sb, Se et As, il apparait que les sulfures de génération 2a ont été partiellement dissous par réactions couplées de dissolution-reprécipitation, libérant et transportant l'Au, le Pb, le Bi et le Te sous forme de produits de fusion (« melts ») d'éléments chalcophiles à bas point de fusion et/ou sous forme de colloïdes, pour reprécipiter les sulfures de génération 2b avec l'or libre et les métaux associés le long de fractures et d'auréoles de recristallisation. Le fractionnement indépendant de la masse du soufre (S-MIF) est présent et constant dans tous les échantillons et indique clairement une source unique et une contribution de soufre métasédimentaire archéen aux fluides aurifères responsables de la formation des gisements du MGD au Paléoprotérozoïque, soit environ 800 millions d'années après la formation des roches hôtes. Ce nouvel ensemble de données montre que les réactions fluides-roches, représentées par les multiples épisodes de remobilisation et de reprécipitation ont été importantes pour la formation des zones minéralisées, et que la dévolatilisation des métasédiments et des métavolcanites en profondeur pendant le métamorphisme prograde peut jouer un rôle important dans l'approvisionnement en soufre, élément clé pour le transport de l'or à l'échelle de la croûte terrestre.

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.

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.

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.