L'objectif principal de ce mémoire est de présenter le processus de conception qui a mené à la réalisation d'un prototype expérimental d'hydrogénérateur à ailes oscillantes de 1 kW. Ce dernier a été réalisé de manière à démontrer le potentiel d'extraction d'énergie d'un hydrogénérateur à ailes oscillantes et pour valider certaines prédictions numériques effectuées par le Laboratoire de Mécanique des Fluides Numériques (LMFN) de l'Université Laval. Une revue de littérature présentant tous les hydrogénérateurs à aile(s) oscillante(s) répertoriés dans la littérature et leur principe de fonctionnement est d'abord présentée. Ensuite, la conception et l'expérimentation d'un premier prototype de petite taille est montrée. Des observations pertinentes quant à la conception du prototype de 1 kW sont aussi incluses. Suit alors la description explicite de la conception de l'embarcation motorisée permettant de faire fonctionner le prototype de 1 kW, la présentation du choix de la topologie de mécanisme et son processus d'optimisation et une description détaillée du prototype de 1 kW lui-même. Les résultats expérimentaux obtenus avec le prototype de 1 kW sont présentés de façon sommaire. Une analyse approfondie de ces derniers sera effectuée par Kinsey ([1]).

L'objectif principal de ce mémoire est de mettre à jour les étapes parcourues dans la synthèse de mécanismes pour une génératrice éolienne/hydrolienne à ailes oscillantes. En effet, l'utilisation d'ailes oscillant dans un écoulement d'air ou d'eau afin d'en extraire de l'énergie nécessite un contrôle par un mécanisme contraignant cette dite oscillation. Une architecture doit donc être conçue afin d'imposer le mouvement adéquat de tangage et de pilonnement aux ailes, en plus de transférer aux alternateurs l'énergie cinétique extraite. De ce fait, l'analyse de l'aile oscillante à auto-compensation, un mécanisme proposé dans la littérature, est d'abord conduite. Deux lacunes sont ainsi identifiées : le manque d'adaptabilité du système à l'environnement et l'impossibilité d'atteindre des amplitudes de tangage supérieures à 62 °. S'ensuit alors une recherche de solutions pouvant procurer les mouvements établis qui se solde en l'élaboration de deux architectures à deux degrés de liberté : Valkyrie 2 et AEGIR. L'examen de celles-ci révèle ensuite que Valkyrie 2 est plus difficile à contrôler qu'AEGIR, puisque ce dernier possède une commande découplée. De son côté, Valkyrie 2 utilise seulement des barres et des liaisons rotoïdes, ce qui est un avantage du point de vue de l'efficacité énergétique. En étudiant par la suite la possibilité d'utiliser ces mécanismes dans des systèmes tandem, l'avantage du contrôle découplé d'AEGIR se fait plus important, puisqu'il est plus facile de diriger l'orientation des ailes par le pivot central à l'aide de courroies plutôt que par de nombreuses membrures. Finalement, la proposition d'ajouter des masses décentrées au système AEGIR tandem est apportée afin d'améliorer l'allure de la courbe de puissance en sortie du système.

La présente thèse propose une étude expérimentale du décollement dans le diffuseur d’un modèle de turbine hydroélectrique bulbe. Le décollement se produit quand la turbine est opérée à forte charge et il réduit la section effective de récupération du diffuseur. La diminution de la performance du diffuseur à forte charge engendre une baisse brusque de l’efficacité de la turbine et de la puissance extraite. Le modèle réduit de bulbe est fidèle aux machines modernes avec un diffuseur particulièrement divergent. Les performances de la turbine sont mesurées sur une large gamme de points d’opération pour déterminer les conditions les plus intéressantes pour l’étude du décollement et pour étudier la distribution paramétrique de ce phénomène. La pression est mesurée le long de l’aspirateur par des capteurs dynamiques affleurants alors que les champs de vitesse dans la zone de décollement sont mesurés avec une méthode PIV à deux composantes. Les observations à la paroi sont pour leur part faites à l’aide de brins de laine. Pour un débit suffisant, le gradient de pression adverse induit par la géométrie du diffuseur affaiblit suffisamment la couche limite, entraînant ainsi l’éjection de fluide de la paroi le long d’une large enveloppe tridimensionelle. Le décollement instationnaire tridimensionnel se situe dans la même zone du diffuseur indépendamment du point d’opération. L’augmentation du débit provoque à la fois une extension de la zone de décollement et une augmentation de l’occurrence de ses manifestations. La position et la forme du front de décollement fluctue significativement sans périodicité. L’analyse topologique et celle des tourbillons des champs de vitesse instantanés montrent une topologie du front de décollement complexe qui diffère beaucoup d’une réalisation à l’autre. Bien que l’écoulement soit turbulent, les tourbillons associés aux foyers du front sont clairement plus gros et plus intenses que ceux de la turbulence. Cela suggère que le mécanisme d’enroulement menant aux tourbillons du décollement est clairement distinct des mécanismes de la turbulence. / This thesis presents an experimental investigation of flow separation inside the diffuser of a small scale model of a bulb turbine. The flow separation occurs when the turbine is operated at high discharge and it reduces the diffuser effective area. In the case of bulb turbines, the kinetic energy recovered by the diffuser represents an important part of the total net head available for the runner energy extraction. The decrease of the diffuser efficiency leads to a sudden drop in the turbine efficiency and in the power extraction. The small scale model is faithful to modern turbines with a particularly divergent diffuser. The turbine performances are measured in a large range of operating conditions in order to select the most interesting ones and to investigate the parametric range of the phenomena. The pressure is measured along the diffuser by flush mounted dynamic sensors while the velocity fields inside the separation zone are obtained by a two-component PIV method. Separation observations on the wall are done using tufts. For a sufficient flow rate, the adverse pressure gradient induced by the diffuser geometry sufficiently weakens the boundary layer, thus leading to fluid ejection from the wall along a large three-dimensional envelope. The three-dimensional unsteady flow separation zone is located in the same area independently of the operating points. The flow rate increase leads to a wider separation zone occurring more frequently. The separation front fluctuates significantly both in location and in shape with no periodicity. Topological and vortex analyses on instantaneous velocity fields show a complex separation front topology which differs greatly from one realisation to another. Despite the highly turbulent flow, the separation front vortices are definitely bigger and more intense than turbulent vortices. It suggests that the roll-up mechanisms leading to separation surface vortices appear to be distinct from those of turbulent vortices.

Le développement d'une sonde de traînée pour la mesure instantanée de la vitesse de l'aluminium en fusion dans une cellule d'électrolyse est l'objectif principal de ce projet de maitrise. L'effet de l'hydrodynamique des fluides sur l'efficacité énergétique d'une cuve d'électrolyse a incité les chercheurs à étudier ce phénomène par le biais de simulations numériques depuis plusieurs décennies. Cependant, il existe peu de données expérimentales disponibles pour valider les profils de vitesse obtenus à partir de ces simulations numériques. Une revue de la littérature sur les méthodes possibles de mesure de la vitesse dans une cellule d'électrolyse en aluminium a été réalisée. Cependant, la plupart des méthodes suggérées dans la littérature ne permettent pas de résoudre les problèmes de vélocimétrie à l'intérieur de la cellule d'électrolyse. La force du champ magnétique généré par le passage d’un fort courant dans les cuves d’électrolyse affecterait fortement les instruments électroniques. Une température élevée de l’aluminium allant jusqu'à 960°C est la raison pour laquelle de nombreuses techniques suggérées sont écartées. En outre, un environnement agressif de la cuve limite le choix des matériaux pour les méthodes avec des pièces immergées dans la cellule. Compte tenu de ces difficultés, une nouvelle tentative est entreprise pour mesurer la vitesse de l'aluminium en fusion avec une sonde de traînée inspirée par un dispositif éprouvé initialement utilisé pour la mesure de la vitesse du vent. Ce nouveau dispositif est conçu pour minimiser les effets de vortex, pour augmenter le coefficient de traînée et pour être applicable dans les cuves d’électrolyse. Ce document présente la sonde adaptée à la cuve d'électrolyse, son procédé de calibration et de validation dans un canal d'eau à température ambiante et les résultats obtenus en cuve d'électrolyse. Enfin, une conclusion sur les résultats obtenus et une suggestion de travaux futurs possibles est fournie / The development of a drag probe for instantaneous velocity measurement of molten aluminium in an electrolysis cell is performed in this master project. Because of the strong effects of flow hydrodynamics on energy efficiency of electrolysis cell, researchers studied this phenomenon through numerical simulations for decades. However, there is little experimental data to validate velocity profiles obtained from these numerical simulations. A literature review over possible methods for velocity measurement in aluminium electrolysis cell has been done. However, the challenges of velocimetry inside the electrolysis cell are beyond capabilities of most of methods suggested in literature. The high magnitude of the magnetic field from the passage of high DC current through the process would strongly affect electronic instruments. The high temperature of molten aluminium (up to 960°C) is reason to discard many suggested techniques. Besides, the aggressive environment of reduction cell limits the choice of materials for measurement methods with parts immersing inside the cell. Considering all those difficulties, a new attempt is undertaken to measure molten aluminium velocity with a drag probe inspired from a proven device first used for wind velocity measurement. The new device is designed to minimize vortex shedding, to increase the drag coefficient and to be applicable in the harsh environment of aluminium reduction cell. This document presents the probe developed for the electrolysis cell, its calibration in a water channel at room temperature, and the velocity results obtained within electrolysis cell.

L'écoulement de von Karman est l'écoulement produit par la rotation de deux disques coaxiaux placés en vis-à-vis. Dans la première partie de cette thèse, nous présentons un certain nombre de résultats expérimentaux visant à la caractérisation des fluctuations à grande échelle d'un écoulement d'allongement sensiblement unitaire. Dans certaines circonstances, nous mettons notamment à jour une "bifurcation globale" ainsi qu'un comportement "intermittent" de la structure de l'écoulement. Nous présentons par ailleurs une relation qui relie les statistiques des fluctuations de l'écoulement aux valeurs des couples fournis par les moteurs. Enfin, nous ébauchons l'étude de l'influence de la force de Coriolis sur l'écoulement. La deuxième partie de cette thèse présente des résultats numériques et expérimentaux portant sur l'étude des caractéristiques magnéto-hydrodynamiques de l'écoulement, et notamment sur sa capacité à produire un champ magnétique par "effet dynamo". Les résultats expérimentaux ont été obtenus à haut nombre de Reynolds magnétique dans le montage "VKS", qui utilise du sodium liquide.

In this thesis, we construct a new hydrodynamical algorithm able of handling general compressible reactive flow problems, based on a finite-volume method inspired by the original MUSCL scheme of van Leer (1979). The algorithm is of second-order in the smooth part of the flow and avoids dimensional splitting. It uses MPI to achieve parallelism, and includes an astrophysical equation of state and a nuclear reaction network. It proves to be robust to tests cases. In particular it reproduces quite well the reactive and non-reactive results obtained with two different numerical methods (Fryxell & al. 1989, Busegnies & al. 2007). Moreover the time-dependent results are in agreement with the corresponding steady state solution. This gives us confidence in applying it to an astrophysical situation which has never been studied, the propagation of a detonation in conditions relevant to superbursts. The algorithm is described in (Noel & al. 2007). In a firt step we obtain the detonation profiles in pure carbon and in a mixture of carbon and iron. In both cases we underline the large difference between the total reaction length and the length on which some species burn. This difference leads to enormous numerical difficulties because all the length scales cannot be resolved at the same time in a single simulation. We show that the carbon detonation might be studied in a partial resolution approach like the one of Gamezo & al. (1999). In a second step we construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegrations of heavy elements, like ruthenium, which are thought to occur during superbursts in mixed H/He accreting systems. Using this new nuclear network we simulate detonations in mixture of carbon and ruthenium. An interesting feature is that, in this case, all the reaction lengths can be resolved in the same simulation. This makes the C/Ru detonations easier to study in future multi-dimensional simulations than the pure carbon ones (Noel & al. 2007b). Finally we perform some numerical experiments which show that our algorithm is able to deal with initially inhomogeneous medium, and that the multi-dimensional simulations are attainable even if they are quite computational time consuming. - B. Van Leer, J. Comp. Phys., 21, 101, 1979 - Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989 - Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007 - Gamezo, V.N., Wheeler, J.C., Khokhlov, A.M., and Oran, E.S., ApJ, 512, 827, 1999 - Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007 - Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b / Un algorithme parallèle basé sur une méthode aux volumes finis inspirée du schéma MUSCL de Van Leer (1979) a été construit. Il a été développé sur base de la méthode de Lappas & al. (1999) qui permet de résoudre simultanément toutes les dimensions spatiales. Cette méthode se base sur la construction de surfaces appropriées dans l'espace-temps, le long desquelles les équations de bilan se découplent en équations plus simples à intégrer. Cet algorithme est actuellement le seul à éviter le "splitting" des dimensions spatiales. Dans les modèles conventionnels (PPM, FCT, etc.), l'intégration spatiale des équations est réalisée de manière unidimensionnelle pour chaque direction. Un réseau de réactions nucléaires ainsi qu'une équation d'état astrophysique ont été inclus dans l'algorithme et celui-ci a ensuite été soumis à une grande variété de cas tests réactifs et non réactifs. Il a été comparé à d'autres codes généralement utilisés en astrophysique (Fryxell & al. 1989, Fryxell & al. 2000, Busegnies & al. 2007) et il reproduit correctement leurs résultats. L'algorithme est décrit dans Noël & al. (2007). Sur base de cet algorithme, les premières simulations de détonation dans des conditions thermodynamiques représentatives des Superbursts ont été réalisées. Différentes compositions du milieu ont été envisagées (carbone pur, mélange de carbone et de fer, mélange de carbone et de cendres du processus rp). Dans la plupart des systèmes où des Superbursts ont été observés, la matière accrétée est un mélange d'hydrogène et d'hélium. Dans ce cas, des phases de combustion précédant le Superburst produisent des nucléides plus lourd que le fer (Schatz & al. 2003). Ces nucléides peuvent être photodésintégrés durant le Superburst. Pour prendre en compte ces réactions endothermiques de photodésintégration, nous avons construit un nouveau réseau réduit de réactions nucléaires qui a été incorporé dans l'algorithme hydrodynamique (Noël & al. 2007b). Ce réseau réduit reproduit globalement l'énergétique d'un réseau complet et a permis de faire la première simulation numérique de détonation dans des conditions caractéristiques de systèmes accréteurs d'un mélange hydrogène-hélium. Finallement quelques simulations multidimensionelles préliminaires ont éte réalisées. - Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007 - Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989 - Fryxell, B.A., Olson, K., Ricker, P. & al., ApJS, 131, 273, 2000 - Lappas, T., Leonard, A. and Dimotakis, P.E., SIAM J. Sci. Comput., 20, 1064, 1999 - Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007 - Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b - Röpke, F. K. PhD thesis, Technischen Universitat Munchen, 2003 - Schatz, H., Bildsten, L., Cumming, A. and Ouellette, M., Nuclear Physics A, 718, 247, 2003 - Van Leer, B. Comp. Phys., 21, 101, 1979 - Weinberg, N.N. and Bildsten, L., ArXiv e-prints, 0706.3062, 2007

Ce travail porte sur l'étude du comportement rhéologique de fluide en milieu confiné. Pour cela le levier d'un microscope à force atomique (AFM) est utilisé pour sonder les propriétés rhéologiques d'un fluide confiné entre deux surfaces : la surface d'une sphère collée à l'extrémité du levier et une surface plane sur lequel le fluide est déposé. Le dispositif expérimental est constitué du système de mesure d'un AFM et d'un piézoélectrique permettant d'approcher ou d'éloigner de la sphère la surface plane. Un modèle analytique permet d'extraire les propriétés rhéologiques du fluide confiné à partir de la déflexion du levier induite par le pincement du fluide. Cette méthode a été validée pour les fluides newtoniens. Par contre pour les fluides non-newtoniens comme par exemple la solution de polyacrylamide nous avons trouvé que la viscosité dépend de la distance D et que le cisaillement n'est pas le seul paramètre pertinent pour interpréter les propriétés rhéologiques.

La stabilité linéaire des jets axisymétriques est étudiée dans un contexte global, c'est à dire en prenant en compte sans approximation la géométrie et le non-parallélisme de l'écoulement. Afin de caractériser le "mode préféré" du jet observé expérimentalement et numériquement, différentes analyses ont été mises en oeuvre. Une conjecture étant que ces structures à une résonance entre des perturbations extérieures et le mode propre le moins atténué, une analyse modale de l'écoulement a tout d'abord été effectuée. Cette étude a demandé la mise au point d'une méthode numérique spécifique pour pouvoir traiter les écoulements compressibles de manière efficace. Cependant, les résultats ont montré qu'une représentation modale n'est pas adaptée pour décrire la dynamique des écoulements stable dominés par l'advection. Des modèles simplifiés permettant de mieux comprendre les limites de cette approche sont présentés. Cette dynamique peut cependant être caractérisée par le calcul des perturbations optimales et du forçage harmonique optimal. Cette dernière approche reproduit de manière robuste les observations expérimentales concernant la fréquence et la structure spatiale du "mode préféré". La structure de ce paquet d'onde global est interprétée comme provenant de la coopération entre différentes familles de modes locaux. L'analyse présentée dans cette thèse montre que le "mode préféré" du jet ne vient pas de la résonance d'un mode propre, mais qu'il s'agit en fait d'une pseudo-résonance.

Situé au SE de l’Espagne, le Campo de Cartagena est un cas emblématique extrême des changements hydrologiques et environnementaux causés par l’utilisation intensive des eaux souterraines pour l’agriculture dans les zones semi-arides du pourtour méditerranéen. Le développement agricole a entraîné la surexploitation des horizons les plus profonds de l’aquifère multicouche, tandis que l’augmentation de la recharge liée au retour d’irrigation causait une remontée des niveaux dans la nappe superficielle libre. Un grand nombre de forages multi-crépinés permettent une connexion artificielle entre les différents aquifères, et donc un possible transfert de contaminants d’origine agricole depuis la nappe superficielle vers les nappes profondes. Suite à la montée des niveaux piézométriques dans la nappe superficielle, un débit permanent est maintenant observé dans le réseau hydrographique. S’y ajoutent les rejets incontrôlés de saumures issues d’usines privées de désalinisation d’eau souterraine, aboutissant eux aussi à l'exutoire du système, la lagune de la Mer Mineure. Dans ce contexte, comprendre la complexité de l’évolution du bilan hydrique est un enjeu scientifique majeur. Le premier axe de recherche a consisté en une large tâche de collecte et d’interprétation d’informations sur l’évolution du système aquifère multicouche au cours du siècle passé. Elle a permis de mettre en évidence l’inversion des gradients hydrauliques entre les trois premiers aquifères et des baisses de niveau piézométrique atteignant 500m. Malgré 40 ans de chroniques piézométriques et géochimiques, les possibilités d’interprétation étaient limitées par les rares descriptions techniques complètes de forages. Une méthode reposant sur la méthode d’apprentissage automatique Random Forest a donc été mise en œuvre et a permis d’identifier l’aquifère d’origine de 107 échantillons d’eau souterraine, avec des résultats bien meilleurs à ceux de l’Analyse Discriminante Linéaire (LDA) ou des arbres de décision (CART). Le deuxième axe de recherche a été motivé par la difficulté d’actualiser le bilan hydrique du système aquifère, où le retour d’irrigation s’ajoute à l’infiltration des précipitations comme source de recharge de la nappe superficielle. Des traceurs environnementaux (14C, 13C, 2H, 18O, 3H) ont été employés pour caractériser l’évolution à long terme de la recharge des nappes ainsi que les principaux processus de mélanges, discriminés entre échelle régionale (entre aquifères) et locale (au sein des forages). L’identification des échantillons correspondant à une recharge ancienne (pré-anthropisation) ou à une recharge moderne (post-anthropisation) à permis de quantifier une augmentation des taux de recharge dans la plaine de plus d’un ordre de grandeur en conséquence de la mise en place de l’activité agricole, soit de 10 mm.an-1 à 210 mm.an-1. Le troisième axe de recherche a été la quantification de la décharge sous-marine d’eau souterraine (SGD) vers la lagune de la Mer Mineure à l’aide des isotopes du radon (222Rn) et du radium (223Ra, 224Ra), combinés à une modélisation hydrodynamique de la lagune pour comprendre l’impact des apports anthropiques sur le bilan de radionucléides. Les principales zones de SGD ont été localisées, de même que la zone d’influence du panache de radionucléides issus des eaux de surface et des rejets de saumure. Finalement, les bilans de masse en été et en hiver aboutissent à des flux de SGD totaux (incluant les recirculation salines) de 7.2 à 15.9 108 m3.an-1 (222Rn), 21.9-44.7 108 m3.an-1 (224Ra) et 6.9 108 m3.an-1 (223Ra, hiver). Cette étude d'un cas extrême d’anthropisation présente un large intérêt bien au-delà du système aquifère multicouche du Campo de Cartagena ou de la région de Murcie. Les méthodes mises en place dans cette thèse pourront être appliquées dans d’autres sites méditerranéens où l’exploitation de l’eau souterraine suit une évolution comparable, continue et apparemment inexorable. / The Campo de Cartagena area in the Murcia region (SE Spain) is an emblematic case of the hydrological and environmental changes caused by the intensive use of groundwater for agriculture in semiarid Mediterranean areas. Agricultural development supported by the multi-layer aquifer system has led to the overexploitation of the deep layers, while irrigation return flow and the subsequent increased recharge rates increased water table levels in the unconfined layer. In addition, a large number of boreholes of the area are screened in several aquifers and allow an artificial connection between different groundwater masses. Moreover, as a consequence of the water table increase in the shallow aquifer, a permanent flow appeared in the last kilometres of the surface watershed. Together with the uncontrolled release of brines from private groundwater desalination, it induced a permanent surface flow of water to the main outlet of the system, the Mar Menor lagoon. In this context, understanding the complex evolution of the whole system and how the water balance is affected is a hard task.The first research focus aimed at collecting and reviewing all kind of existing data in order to reconstruct the evolution along one century of the multi-layer aquifer system. It highlighted an inversion of the vertical hydraulic gradient between aquifers and the decrease of water table levels up to 500 m for the deepest layer. Facing the difficult identification of the origin of groundwater samples, a method based on the Random Forest (RF) machine learning technique was developed. Despite the difficulty of an unequal dataset, accuracy over 90% was reached and 107 groundwater samples of unknown origin could be classified. Results were much better than using Linear Discriminant Analysis (LDA) or Decision Trees (CART).The second research axis was motivated by the difficulty to update the water balance of the multi-layer aquifer where irrigation return flow represents an additional source of recharge, added to the limited rainfall infiltration. Environmental tracers (14C, 13C, 2H, 18O, 3H) were combined to high-resolution temperature loggings to investigate the long-term evolution of recharge in the Campo de Cartagena aquifer system and discriminate local mixing processes (infra-boreholes) from regional mixing processes (between aquifers). Both pre-anthropization and post-anthropization recharge regime could be identified and quantified. Before the development of agriculture, recharge varied from 17 mm.a-1 in the mountain ranges to 6 mm.a-1 in the plain. In response to the increase of agricultural activity, recharge fluxes were amplified up to 210 mm.a-1 in irrigated areas.The third research axis consisted in quantifying submarine groundwater discharge (SGD). In order to decipher the influence of the different water sources on the Mar Menor, a radon (222Rn) and radium (223Ra, 224Ra) survey was combined with the hydrodynamic modeling of the lagoon. The areas of influence of a plume of radionuclides from the river were identified, the main areas of SGD were located and a location for a submarine emissary was proposed. Mass balances in winter and summer seasons provided total yearly SGD fluxes of water of 7.2-15.9 108 m3.a-1 (222Rn), 21.9-44.7 108 m3.a-1 (224Ra) and 6.9 108 m3.a-1 (223Ra, in winter). Water level effect, rather than tidal pumping, was identified as the main driver for recirculated saline groundwater, while fresh submarine groundwater discharge from the aquifer was about 1% of total SGD. As the case study is an extreme case of anthropization, these results present a wide interest that is not limited to the Murcia region or the Campo de Cartagena aquifer. The methods developed in this thesis might be used in other Mediterranean sites where groundwater exploitation seems to follow a continuous and inevitable increase.

L’objectif de ce travail est d’étudier l’influence de jets d’impact libres utilisés comme mode de mise en contact des phases dans les procédés de recristallisation utilisant un fluide supercritique comme anti-solvant. L’influence des variations de plusieurs paramètres opératoires sur l’hydrodynamique des jets et sur les caractéristiques des poudres de Sulfathiazole a été étudiée. Les paramètres sont la vitesse des jets (de 0,25 m.s-1 à 25,92 m,s-1), le rapport molaire solvant/CO2 (de 2,5 % à 20 %), la température (de 313 K à 343 K), la pression (10 MPa à 20 MPa) et la concentration du soluté dans la solution (de 0,5 % à 1,8 %). Les conditions de mélange ont été caractérisées par l’estimation des puissances dissipées par les jets d’impact, variant de 0,1 à 158 W.kg-1 dans les conditions étudiées. Les résultats ont montré que la vitesse des jets et la sursaturation sont deux paramètres-clés contrôlant la cristallisation avec un effet prépondérant de la sursaturation. Les cristaux obtenus ont des tailles, distributions de taille, faciès et nature polymorphique différents selon les conditions. La forme polymorphique la plus stable a été obtenue pure ou en mélange. Une comparaison avec le procédé SAS classique a montré que les particules sont significativement plus petites lorsqu’elles sont cristallisées avec les jets d’impact, ce qui confirme que ce dispositif créant un mélange plus intense, permet d’accélérer la cinétique de nucléation.Ce travail devrait contribuer à une meilleure maîtrise des procédés de cristallisation en milieu supercritique. / The aim of this work is to study the influence of free impinging jets used for the fluids’introduction in supercritical anti-solvent (SAS) processes. The influence of the variations of several operating parameters upon jets’ hydrodynamics and upon the powder characteristics is studied. Parameters are jets velocity (0.25 m.s-1 to 25.92 m.s-1), molar ratio solvent / CO2 (2.5 % to 20 %), temperature (313 K to 343 K), pressure (10 MPa to 20 MPa) and solute concentration in the organic solution (0.5 % to 1.8 %). Mixing conditions have been characterized estimating the dissipated powers of the impinging jets, varying from 0,1 à 158 W.kg-1 in the studied conditions. The results showed that jets velocity and supersaturation are two key-parameters controlling the crystallization with a preponderant effect of supersaturation. The obtained crystals have different size, particle size distribution, habit and polymorphic nature depending on the conditions. The most stable polymorphic form has been obtained pure or in mixture.A comparison with the classical SAS process showed that mean particle sizes are significantly smaller with impinging jets device proving that this device, creating a more efficient mixing, enhances the nucleation kinetics.This work may contribute to a better control of processes of crystallization in supercritical media.