Cette thèse a pour objectif de détecter les structures turbulentes aux grandes échelles présentes dans une couche de mélange de type Rayleigh-Taylor incompressible à faible nombre d'Atwood. Diverses grandeurs statistiques conditionnées par la présence de ces structures ont été obtenues, et il est désormais possible de les comparer avec les résultats des modèles statistiques turbulents dits bi-structure, tel le modèle 2SFK développé au CEA. Afin de réaliser les simulations numériques directes du mélange turbulent, un code numérique tridimensionnel incompressible à densité variable a été développé. Ce code a été parallélisé dans les trois directions. Plusieurs méthodes de détection de structure ont été conçues et testées. Bien que toutes ces méthodes présentent différents intérêts, seule la plus efficace vis-à-vis de nos critères de détection a été gardée pour faire des simulations à forte résolution (plus d'un milliard de mailles, 1024^3). Un filtrage temporel de la vitesse verticale est utilisé dans cette méthode de détection afin de : 1) corriger les distorsions dues aux points d'arrêt et zones de recirculation dans l'écoulement, 2) minimiser l'effet de la turbulence aux petites échelles et mieux mettre en évidence les grandes échelles, 3) introduire un effet mémoire permettant de prolonger la bimodalité du champ de détection depuis les zones laminaires extérieures jusqu'au centre de la zone de mélange turbulent. Plusieurs simulations numériques directes 1024^3 ont été effectuées. Les résultats viennent conforter ceux obtenus avec le modèle bi-structure 2SFK et justifient une étude plus poussée des grandeurs statistiques en vue de sa validation. / This thesis aims at detecting large-scale turbulent structures in incompressible Rayleigh-Taylor mixing layers at low Atwood number. Various statistical quantities conditioned by structure presence have been obtained and it is now possible to compare them with results from two-structure statistical turbulent models such as the 2SFK model developed at CEA. In order to produce direct numerical simulations of the turbulent mixing, a three-dimensional, incompressible, variable-density numerical code was developed. This code is parallelized in the three directions. Several structure detection methods have been designed and tested. Although all these methods are of interest, only the most efficient with respect to our detection criteria has been retained for simulations at high resolution (over a billion cells, 1024^3). A time filtering of vertical velocity is used in this method to: 1) correct distortions due to stagnation points and recirculation zones in the flow, 2) minimize small-scale turbulence effects and better highlight large-scales, 3) introduce a memory effect in order to extend bimodality of the detection field from the external laminar zones up to the centre of the turbulent mixing zone. Several direct numerical simulations at 1024^3 have been achieved. Results support those obtained with two-structure 2SFK model and justify further studies for its validation.

Dans le contexte péri-urbain à dominance agricole de l’Ile de France, les pressions anthropiques exercées sur les rivières jouent autant sur la qualité que sur le débit des cours d'eau. Située au sud-ouest de Paris, l’Yvette draine un bassin versant de 202 km²et à la géologie homogène. Elle est alimentée directement par des stations d'épuration (STEP) et par la nappe des sables de Fontainebleau qui représente la principale source d’eau. Sur ce type de bassin, la gestion durable du système riverain repose sur la connaissance de la distribution des flux nappe – rivière et sur l’impact de cette répartition sur la qualité des cours d’eau.La dynamique des flux d’eau a été suivie grâce à l’implantation de stations hydrométriques aux points clés du réseau hydrographique (i.e. exutoires des principaux affluents, aval des confluences sur l’Yvette). Les flux chimiques ont été étudiés par analyses d’échantillons d’eau recueillis lors de campagnes de terrain effectuées en période d’étiage sur les cours d’eau. Un modèle conceptuel global, calibré à partir des chroniques de débit, a permis (i) d’estimer des variations spatiales de la recharge de la nappe (60 – 160 mm.an-1), et (ii) d’établir une répartition journalière des parts d’eau de nappe, de STEP et de ruissellement en rivière, pour la période 2001-2014. Par ailleurs, l’analyse des traceurs géochimiques (e.g. Cl-, NO3-, SO42-) et isotopiques (222Rn, δ18Oeau, δ2Heau) soutient la prédominance, aussi bien quantitative que qualitative, de la nappe sur les rivières.D’un point de vue hydrodynamique, le contexte géomorphologique homogène procure aux rivières une dynamique similaire avec des épisodes de crues très courts (de l’ordre de quelques heures) et des périodes d’étiages marqués, quelle que soit la saison. L’étude par sous-bassin a mis à jour une différence entre les bassins topographiques et les bassins d’écoulement souterrains, créant des déficits hydriques sur l’amont de certains cours d’eau (e.g. Mérantaise, Ru des Vaux) au profit d’autres (e.g. Rhodon). Le parallèle entre bassins topographique et souterrain n’a pu se faire qu’au niveau du cours principal de l’Yvette. Sur la période 2001-2014, le débit de l’Yvette provient en moyenne à 55 % de la nappe, à 38 % du ruissellement et à 8 % des STEP. En période de basses eaux, la contribution des STEP reste sensiblement identique tandis que la nappe constitue la principale alimentation des rivières (90 %), contrôlant donc leur qualité. Mais la composition chimique de ce soutien souterrain n’est pas homogène. Pour déterminer l’origine de ces disparités, un travail à plus petite échelle a été conduit sur un affluent majeur de l’Yvette (le Rhodon). La décharge de la nappe en rivière y est bimodale : 15 % arrive par voie souterraine et 85 % transite par les milieux humides en surface. Au sein des zones humides sont observées de fortes hétérogénéités dans les écoulements et leur chimie, avec des flux de subsurface totalement épurés en nitrates et du ruissellement riche en nutriments. La prédominance de des flux de surface réduit fortement le pouvoir épurateur des milieux humides, avec seulement 6 % des eaux de nappe épurées avant leur décharge en rivière. Cette faible efficacité renforce l’impact direct de la qualité de la nappe sur les rivières. Sur le bassin, la nappe est stratifiée par l’intervention de nombreux processus de recharge (infiltrations d’eau d’étang, de STEP, de zones humides de plateau). Les petits sous-bassins (< 50 km²) semblent dominés par cette stratification, ainsi que par les apports de STEP qui contribuent aux flux chimiques des rivières à hauteur de 30 à 50 %. Le cours principal de l’Yvette, dans sa partie aval, dépend des mélanges issus des confluences avec ses affluents. Les outils utilisés ici ont donc mis en avant les différences de fonctionnement des bassins selon l’échelle concernée, tout en reliant l’importance de la compréhension du système à petite échelle pour la gestion à grande échelle. / In the peri-urban context with large farm practices of the Paris region, anthropogenic pressures on streams and rivers impact both their flow and quality. Located in the southwestern of Paris, the Yvette stream drains a watershed of 202 km², in a homogeneous geological context. The stream is sustained by waste water treatment plant (WWTP) and by the Fontainebleau sands aquifer which represents the main source of water.In this context, the sustainable management of the stream and its riparian zones needs the assessment of groundwater discharge locations and chemistry, and their impact on the stream water quality.Water fluxes in the stream and in its main tributaries have been recorded at 11 stations. Dissolved elements distributions have been studied by field works in low-flows conditions. A lumped model, calibrated with stream flow hydrographs, has allowed (i) the estimation of groundwater recharge rates (60 – 160 mm.y-1) and (ii) the distribution of stream flow contributions between 2001 and 2014 (i.e. 55 % from groundwater, 38 % from runoff and 8 % from WWTP). Moreover, the use of geochemical and isotopic tracers (e.g. Cl-, NO3-, SO42-, 222Rn, δ18Owater, δ2Hwater) confirmed the large impact of groundwater on stream flow and quality.The hydrological dynamics of the Yvette stream and its tributaries are similar because of the homogeneity of the geomorphological context of the catchment. The flood events are really fast (on the order of hours) and low-flow conditions can be observed in every season. The discretization of the Yvette catchment in 11 sub-watersheds highlighted a divergence between the topographical and hydrogeological catchments. This generates a water deficit in some tributaries heads (e.g. Mérantaise, Ru des Vaux) in favor of some others (e.g. Rhodon). The matching limits of these catchments can only be made on the Yvette channel. During the 2001 – 2014 period, considering a fairly constant WWTP contribution, the groundwater discharge provides more than 90 % of the stream flow during low-flow conditions. This highlights a stream quality dominated by groundwater discharge. However, groundwater chemistry largely varies over the catchment. To study this phenomenon, a small-scale field work has been directed on the Rhodon stream, a main tributary of the Yvette River. At this scale, the groundwater discharge appeared to have two main pathways: 15 % come from the streambed while 85 % transit in wetlands before going to the stream. Within these wetlands, a large part of groundwater flows at the surface and has high levels of nutrients, whereas a small part (6 %) goes through sediments and is nutrients free. reduces The purifying power of these wetlands, significantly reduced by this distribution, is likely related to historical drainage and reinforces the direct connection between the stream and groundwater quality. On the Yvette catchment, the stratification of groundwater chemistry is due to several recharge processes (e.g. pond water infiltration, WWTP effluents infiltration, wetlands water infiltration). Small sub-watersheds (< 50 km²) seem to be dominated by (i) this stratification, and (ii) the WWTP discharges which can contribute from 30 % to 50 % to dissolved fluxes in the streams. The main channel of the Yvette stream primarily depends of the mixing of its tributaries waters. The methods developed here allowed highlighting (i) differences of catchment functioning according to the scale used, and (ii) the utmost importance of the small-scale assessment to understand/decipher and manage streams at larger scales.

Mémoire ou thèse en cotutelle / L'électrodialyse (ÉD) est largement utilisée dans diverses industries, y compris l'industrie alimentaire (pour la désacidification des jus de fruits et du vin, la déminéralisation du lactosérum, etc.). Cependant, les membranes utilisées dans le traitement des solutions alimentaires sont souvent sujettes au colmatage. Le colmatage, dépôt de différentes espèces sur la surface ou à l'intérieur du matériau, dégrade considérablement les performances et la durée de vie des membranes, ce qui conduit à la nécessité de leur remplacement et de leur nettoyage régulier. À cet égard, la formation de dépôts à la surface des membranes entraîne une augmentation du coût de l'ensemble du procédé d'électrodialyse (ÉD). La détermination des moyens pour réduire le colmatage est une problématique réelle, en particulier dans l'industrie alimentaire, qui traite des solutions complexes. Cette thèse consistait à déterminer la contribution du débit de la solution dans la cellule d'ÉD à l'échelle laboratoire, ainsi que des modes de courant électrique sur la cinétique de colmatage protéique (formation de dépôts protéiques) des membranes bipolaires (MB) lors du traitement d'électroacidification par ÉD de solutions laitières. La nouveauté scientifique et l'originalité du travail tiennent au fait qu'à notre connaissance, il n'y a jamais eu d'évaluation de l'effet combiné du débit de solution et d'utilisation de régimes de champ électrique pulsé (CÉP), alternant des impulsions de courant et des pauses (courant égal à zero), sur cette cinétique de colmatage . De plus, dans le cadre de cette étude, un modèle mathématique a été développé, pour la première fois, afin de tenter de décrire la cinétique de formation et de dissolution des dépôts de protéines au cours du processus d'ÉDMB. L'objectif principal de la première étape de la thèse était de déterminer l'effet du débit de solution et du mode de courant électrique sur la cinétique de colmatage protéique des membranes lors du procédé d'électrodialyse à membranes bipolaires du lait écrémé. Deux modes de courant électrique ont été utilisés : le mode de courant traditionnel, à savoir le courant continu (CC) et le CÉP. Un seul régime d'impulsion-pause de CÉP a été étudié, dans lequel 10 s d'impulsion alternaient avec 50 s de pause. De plus, cinq débits différents de solution ont été testés et correspondant aux nombres de Reynolds (Re) suivants : 162, 242, 323, 404 et 485. Les résultats obtenus en mode CÉP ont été comparés aux données obtenues en mode CC sur une même base de nombre de charges transportées pendant l'expérience. Des dépôts de protéines sur la couche échangeuse de cations de la membrane bipolaire (MB) en mode CC ont été trouvés, quel que soit le débit de solution appliqué. Les résultats ont montré que l'utilisation de CÉP empêchait presque complètement la formation de colmatage protéique sur la surface de la MB, quel que soit le débit de la solution dans le canal de la cellule d'ÉD. Il a ainsi été démontré que la minimisation du colmatage est principalement liée à l'utilisation de CÉP en comparaison avec l'augmentation du débit de solution. L'objectif de la deuxième étape du projet était d'optimiser le régime actuel de CÉP afin de minimiser davantage le colmatage des protéines sur la surface de la MB. Dans cette partie du travail, l'effet combiné du débit de la solution et de différents régimes d'impulsion-pause de CÉP a été étudié sur la quantité de protéines déposée durant le processus d'ÉDMB d'une solution modèle de caséinate de sodium. Ainsi, trois débits (correspondant à des nombres de Reynolds de 187, 374 et 560) pour la solution d'alimentation et quatre régimes d'impulsion-pause (avec des rapports Ton/Toff égaux à 10 s-10 s, 10 s-20 s, 10 s-33 s et 10 s-50 s) ont été évalués à une densité de courant de 5 mA/cm². Un colmatage protéique a été observé sur la couche échangeuse de cations de la MB dans toutes les conditions expérimentales considérées, et cela quels que soient le régime de courant et le débit de solution. Il a cependant été constaté que l'augmentation du temps de pause et du débit de la solution de caséinate avaient un impact positif sur la minimisation du colmatage protéique à la surface cationique de la MB. Une pause longue et un débit élevé ont contribué à une diminution plus efficace de la concentration de protons et d'anions caséinates à la surface de la MBP: Une diminution de 87 % de la quantité de colmatage de la membrane a été obtenue avec 50 s de temps de pause en CÉP et un débit correspondant à Re = 374. Une augmentation supplémentaire du temps de pause en CÉP (au-dessus de 50 s) ou du débit (au-dessus de Re = 374) n'a pas entraîné de diminution significative de la quantité de colmatage. La dernière étape du travail a été la modélisation mathématique du processus à l'étude afin de mieux comprendre le mécanisme d'attachement et de croissance du colmatage au cours de l'ÉDBM. Dans cette partie du projet, un modèle unidimensionnel d'ÉDBM d'une solution de caséinate de sodium a été développé. Il est basé sur l'équation de Nernst-Planck sous la condition d'électroneutralité locale. Le modèle prend en compte le transfert de 6 espèces ioniques (Na⁺, K⁺, H⁺, Cl⁻, OH⁻ et anions caséinate (Cas⁻)), ainsi que des molécules non chargées de caséine, qui précipitent (HCas⁰). La géométrie du modèle cible la surface de la couche échangeuse de cations de la MB d'un côté et le volume de la solution agitée de l'autre. Les paramètres hydrodynamiques du système sont pris en compte à travers l'épaisseur de la couche de diffusion. Le modèle développé permet le calcul des concentrations ioniques dans la solution traitée et décrit qualitativement la cinétique de colmatage des protéines à la surface de la MB. Il a ainsi été montré qu'une augmentation du temps de pause en CÉP entraîne une diminution de l'épaisseur des dépôts protéiques en surface. Les meilleurs résultats en termes de minimisation des dépôts ont été obtenus pour le régime de CÉP de 10 s-50 s, ce qui est cohérent avec les résultats expérimentaux obtenus précédemment. Le plus grand effet anticolmatant a été obtenu avec une durée d'impulsion de courant de 10 s et une pause de 50 s en combinaison avec un nombre de Reynolds de 560. La combinaison optimale du mode de courant et paramètre hydrodynamique a réduit la masse du dépôt de caséine à la surface de la MB de 78 % par rapport au régime de CÉP de 10 s-10 s combiné avec un nombre de Reynolds de 187. Les résultats obtenus dans cette thèse de doctorat ont démontré pour la première fois les impacts du mode hydrodynamique dans la cellule d'ÉD, du mode de courant électrique ainsi que du régime de courant en CÉP sur la cinétique de colmatage dans la cellule d'ÉDBM lors de l'électroacidification de solutions laitières. De plus, le modèle développé a permis de prédire qualitativement la quantité de dépôt à la surface de la membrane en fonction du débit (nombre de Reynolds) et de la combinaison impulsion/pause de CÉP appliqués. Les perspectives futures du travail sont la transition vers un système d'ÉD semi-industriel, l'utilisation d'autres types d'espaceurs, l'utilisation de membranes modifiées avec une hétérogénéité géométrique, le test d'autres modes impulsion-pause de CÉP. De plus, concernant la modélisation, le modèle développé pourrait être amélioré en passant par un système bidimensionnel, prenant en compte les molécules protéiques chargées positivement, ainsi que la présence d'un espaceur dans le canal d'acidification/déminéralisation. / Electrodialysis (ED) is widely used in various industries, including food industry (for deacidification of fruit juices and wine, demineralization of whey, etc.). However, membranes used in the processing of food solutions are often prone to fouling due to the solution complexity. Fouling, deposition of different species on the surface or inside the material, significantly degrades the performance and lifetime of membranes, leading to their regular replacement and cleaning. Hence, the formation of deposits on the surface of membranes leads to an increase in the cost of the entire electrodialysis (ED) process. The determination of ways to decrease fouling is a main actual concern, especially in the food industry, which treats complex solutions. This PhD thesis consisted in determining the contribution of the flow rate of the solution in an ED cell at laboratory scale, as well as of the electric current modes on the protein fouling kinetics (formation of protein deposits) on bipolar membranes (BPM) during electroacidification by ED of milk solutions. The scientific novelty and originality of the work are based on the fact that, as far as we know, there have been no assessments of the combined effect of both solution flow rate and use of pulsed electric field (PEF) regimes, in which CC pulses alternate with pauses (zero current), on such protein fouling kinetics. Also, within the framework of the study, a tentative mathematical model was developed, for the first time, to allow to describe the formation and dissolution kinetics of protein deposits during EDBM process. The main objective of the study at the first stage of the thesis was to determine the effect of solution flow rate and electrical current mode on the protein fouling kinetics on membranes during the process of electrodialysis with bipolar membranes of skim milk. Two electrical modes of current were used: the traditional mode of current, namely continuous current (CC) and PEF. One PEF regime was studied, in which 10 s of pulse lapse alternated with 50 s of pause. Additionally, five different flow rates of solution were tested corresponding to the following Reynolds numbers (Re): 162, 242, 323, 404, and 485. The results obtained in the PEF mode were compared with the data obtained in the CC mode with respect to the same number of charges transported during the experiment. Protein deposits on the cation exchange layer of the BPM using CC were found, regardless of the solution flow rate applied. However, the use of PEF allows to prevent the protein fouling formation on the BPM surface almost completely, regardless of the solution flow rate in the channel of the of ED cell. Furthermore, it has been shown that the use of PEF has the main effect in terms of fouling minimization compared to the effect of increasing solution flow rate. The objective of the thesis at the second stage of the project was to optimize the PEF current mode to further minimize protein fouling on the BPM surface. During this stage, the combined effect of solution flow rate and different PEF regimes on the amount of protein deposit during EDBM of a model sodium caseinate solution was investigated. To achieve this objective three different flow rates of the feed solution (corresponding to Reynolds numbers of 187, 374 and 560) and four different pulse-pause regimes (Ton/Toff ratios equal to 10 s-10 s, 10 s- 20 s, 10 s-33 s and 10 s-50 s) during EDBM treatment were evaluated at a current density of 5 mA/cm². Protein fouling was observed on the cation-exchange layer of BPM for all the experimental conditions considered, regardless of the PEF regime and the solution flow rate. It was found that increasing both pause duration and caseinate solution flow rate had a positive impact on protein fouling minimization occurring on the BPM cationic surface. A long pause associated with a high flow rate contribute to a more effective decrease in the concentration of protons and caseinate anions at the BPM surface: a very good membrane performance was achieved with 50 s of pause duration of PEF and a flow rate corresponding to Re = 374 (decrease in the amount of fouling by 87%). A further increase in PEF pause duration (above 50 s) or flow rate (above Re = 374) did not lead to a significant decrease in the amount of fouling. The next step of the work was the mathematical modelling of the process under study for better understanding the mechanism of fouling attachment and growth during EDBM. In this part of the project, a one-dimensional model of EDBM for sodium caseinate solution was developed. It is based on the Nernst-Planck equation under the local electroneutrality condition. The model takes into account the transfer of 6 ionic species (Na⁺, K⁺, H⁺, Cl⁻ , OH⁻ and caseinate anions (Cas⁻)), as well as uncharged molecules of casein, which precipitate (HCas⁰). The geometry of the model is the surface of the BPM cation-exchange layer on the one side and the volume of the stirred solution on the other. The hydrodynamic parameters in the system are taken into account through the diffusion layer thickness. The developed model allows the calculation of the ionic concentrations in the system and to qualitatively describe the kinetics of BPM protein fouling. It has been shown that an increase in the pause duration of PEF leads to a decrease in the thickness of the protein deposits at the BPM surface. The best results in terms of deposit minimization were achieved under the 10 s-50 s PEF regime, which is consistent with the experimental results obtained earlier. The greatest positive antifouling effect was achieved with a current pulse duration of 10 s and a pause of 50 s in combination with Reynolds number of 560. The synergistic combination of optimal current and hydrodynamic modes reduced the mass of casein deposits on the surface of the BPM by 78% compared to the 10 s-10 s PEF regime in combination with Reynolds number of 187. The results obtained in this Ph.D. thesis demonstrated for the first time the impact of the hydrodynamic parameter in the ED cell, the mode of electric current and the PEF regime on fouling kinetics in the process of EDBM of dairy solutions. Additionally, the developed model makes it possible to qualitatively predict the amount of deposit on the membrane surface depending on the flow rate (the Reynolds number) in the desalination compartment and PEF regime applied. The main perspectives of the work are the transition to a semi-industrial ED system, use of other types of spacers, use of modified membranes with geometric heterogeneity, testing of other pulse-pause modes of PEF. Regarding the modeling, the developed model could be improved, namely to transit to a two-dimensional system, taking into account positively charged protein molecules, as well as the presence of a spacer in the acidification/desalting channel.

Thèse ou mémoire avec insertion d'articles. / Le compactage des sols agricoles affecte négativement les performances des cultures et devrait être évité grâce à des stratégies de gestion des sols appropriées. Les recherches antérieures montrent un manque d'études complètes sur l'impact des différents niveaux de compactage du sol sur les propriétés hydrodynamiques. Cette étude s'est concentrée sur l'effet du compactage sur les propriétés hydrodynamiques du sol dans des conditions non saturées et saturées en utilisant le système Hydraulic Property Analyzer (HYPROP). Nous avons étudié l'impact de cinq niveaux de compactage sur des échantillons de sol de sable limoneux prélevés dans un champ de pommes de terre dans les régions agricoles du nord du Québec, au Canada, en utilisant quatre répétitions par niveau. Tout d'abord, la conductivité hydraulique saturée de chaque échantillon a été mesurée en utilisant la méthode de la charge constante. Les données relatives à la courbe de rétention d'eau du sol (CRS) et à la conductivité hydraulique non saturée ont été obtenues à l'aide de l'appareil de mesure de l'évaporation HYPROP et d'un potentiomètre à point de rosée WP4C. Seize modèles hydrauliques de sol du logiciel HYPROP-Fit ont été simulés. Nos résultats soutiennent l'application de la variante Peters-Durner-Iden (PDI) du modèle bimodal sans contrainte de van Genuchten (VGm-b-PDI) pour l'estimation complète du CES sur la base de la métrique de l'erreur quadratique moyenne (EQM). Le modèle VGm-b-PDI a été implémenté dans le logiciel HYDRUS-1D pour calculer la quantité d'irrigation pour différents niveaux de compactage et fournir des estimations plus précises des propriétés hydrauliques du sol. Nous avons simulé des scénarios de compactage et d'irrigation du sol avec le modèle de double porosité et les résultats ont indiqué que le compactage du sol peut fortement influencer les propriétés hydrauliques du sol de différentes manières. La conductivité saturée et non saturée diminue avec l'augmentation du compactage du sol. À un niveau modéré de compactage, la quantité d'irrigation pour la culture de pommes de terre était optimale. Dans l'ensemble, la méthodologie combinée d'HYPROP et d'HYDRUS 1D a donné d'excellents résultats en termes de comportement hydraulique des sols compactés. / Agricultural soil compaction adversely affects crop performance and should be avoided through appropriate soil management strategies. Previous research shows a lack of comprehensive studies on the impact of different levels of soil compaction on hydrodynamic properties. This study focused on the effect of compaction on soil hydrodynamic properties under unsaturated and saturated conditions using the Hydraulic Property Analyzer (HYPROP) system. We studied the impact of five compaction levels on silty sand soil samples taken from a potato field in the northern agricultural regions of Quebec, Canada, using four replicates per level. First, the saturated hydraulic conductivity of each sample was measured using the constant-load method. Soil water retention curve (SWRC) and unsaturated hydraulic conductivity data were obtained using the HYPROP evaporation measurement device and a WP4C dew point potentiometer. Sixteen hydraulic soil models in the HYPROP-Fit software were simulated. Our results support the application of the Peters-Durner-Iden (PDI) variant of van Genuchten's bimodal unconstrained model (VGm-b-PDI) for full SWRC estimation based on the root mean square error (RMSE) metric. The VGm-b-PDI model was implemented in HYDRUS-1D software to calculate the amount of irrigation for different compaction levels and provide more accurate estimates of soil hydraulic properties. We simulated soil compaction and irrigation scenarios with the dual-porosity model and the results indicated that soil compaction can strongly influence soil hydraulic properties in different ways. Saturated and unsaturated conductivity decreases with increasing soil compaction. At a moderate level of compaction, the amount of irrigation for potato cultivation was optimal. Overall, the combined methodology of HYPROP and HYDRUS 1D gave excellent results in terms of the hydraulic behavior of compacted soils.

Ce mémoire s'inscrit dans le cadre plus général d 'un projet multidisciplinaire mené à l'Université Laval sur la technologie des aéro/hydrogénérateurs à ailes oscillantes. Cette partie du projet s'intéresse plus particulièrement à une mesure expérimentale des forces s'exerçant sur une aile oscillante. Pour y arriver, cette étude se découpe en deux phases bien distinctes. Tout d 'abord, il s'est révélé primordial de caractériser et d 'optimiser au mieux l'environnement de travail. Au cours de cette première phase, la mise en place de pertes de charge accompagnée par des mesures de vitesse dans le canal furent effectuées. De même, le système mécanique devait être piloté afin de reproduire le plus fidèlement possible les mouvements demandés. Par la suite, la deuxième phase impliquait la création puis l'utilisation d 'un instrument de mesure. Une comparaison avec des résultats issus de simulation numérique a alors permis de valider le concept de la balance hydrodynamique. Toutefois, plusieurs problèmes furent identifiés et des recommandations faites en vue de travaux futurs.

Thèse ou mémoire avec insertion d'articles. / Ce mémoire présente le développement, à l'aide de simulations numériques d'écoulement, d'une méthode de contrôle des tourbillons en colonne dans une turbine hydraulique axiale de vitesse spécifique élevée opérée en marche à vide. Les observations sur le développement des tourbillons en colonne ont été généralisées à la géométrie d'une turbine Francis de vitesse spécifique moyenne. La roue, qui en raison de son emplacement et de sa géométrie vient interagir avec les tourbillons en colonne, a été retirée de la simulation instationnaire. La zone de recirculation obtenue sans la roue est pratiquement inchangée et des tourbillons en colonne ont été observés. Cela suggère que les tourbillons en colonne dans les turbines axiales et les tourbillons inter-aubes dans les turbines Francis sont liés à la même instabilité hydrodynamique. Selon la géométrie et les paramètres d'opération, différents modes de l'instabilité seraient excités. L'existence des tourbillons en colonne ne semble donc pas sensible à la forme du canal méridien, mais plutôt à la quantité de swirl injectée. Pour contrôler les tourbillons en colonne, le swirl injecté dans la roue est réduit en augmentant l'angle d'ouverture des directrices. Deux mécanismes générant des pertes de charge supplémentaires dans l'écoulement ont été développés pour recréer les conditions d'opération en marche à vide à grande ouverture des directrices. Les mécanismes, similaires à des volets sur des ailes d'avion, sont installés sur les directrices. Des simulations numériques instationnaires SAS-SST ont montré que les deux mécanismes ont permis d'éliminer les tourbillons en colonne dans l'entrefer et dans la roue. Au lieu de cela, six tourbillons inter-aubes co-rotatifs ont été obtenus dans la roue. Aucune structure dommageable n'a été identifiée dans l'aspirateur. Les fluctuations de pression dans la roue et les fluctuations de couple sur les aubes ont été significativement réduites, prouvant ainsi la validité du concept de volets sur les directrices. / This thesis presents the development, using numerical flow simulations, of a mitigation method against columnar vortices in a high-specific-speed axial hydraulic turbine operated at speed-no-load. Observations regarding the development of columnar vortices were generalized to the geometry of a Francis turbine with a medium specific speed. The runner, due to its location and geometry, interacts with the columnar vortices. It was removed from the unsteady simulation. The backflow obtained without the runner is practically unaffected, and columnar vortices were observed in the meridional channel. This suggests that columnar vortices in axial turbines and inter-blade vortices in Francis turbines are linked to the same hydrodynamic instability. Depending on the geometry and operating parameters, different modes of instability would be excited. Therefore, the existence of columnar vortices does not seem to be sensitive to the shape of the meridional channel, but rather to the amount of swirl injected. To mitigate the columnar vortices, the swirl injected into the runner is reduced by increasing the guide vanes opening angle. Two mechanisms that generate additional pressure losses in the flow were developed to recreate the operating conditions at speed-no-load with a larger guide vane opening angle. The mechanisms, similar to spoilers on airplane wings, are installed on the guide vanes. Unsteady SAS-SST numerical simulations showed that the two mechanisms eliminated columnar vortices in the vaneless space and the runner. Instead, six co-rotating inter-blade vortices are obtained in the runner. No damaging structures were identified in the draft tube. The pressure fluctuations in the runner and the blade torque fluctuations were significantly reduced, thus proving the validity of the concept of spoilers on the guide vanes.

Cette maîtrise s'inscrit dans la continuité d'un projet multidisciplinaire de grande envergure se déroulant à l'Université LAVAL, qui a pour but d'analyser les différentes caractéristiques et d'étudier la faisabilité d'une nouvelle technologie de production d'énergie : l'hydrogénérateur à ailes oscillantes. On s'intéresse ici à l'étude expérimentale des forces hydrodynamiques créées par un écoulement d'eau autour de ces ailes. M. Jean-Frederick Faure ayant déterminé dans un précédent mémoire les différentes technologies permettant la mesure de ces forces, ce mémoire présente essentiellement les travaux concernant l'amélioration du montage (installation et étalonnage des nouveaux instruments de mesure) et les résultats et analyses découlant des différents essais réalisés. Des efficacités énergétiques supérieures à 30% ont pu être trouvées et comparées favorablement avec les valeurs de simulations numériques effectuées à l'Université Laval. Une campagne d'essai a en outre permis d'évaluer l'influence des différents paramètres du mouvement (fréquence de cycle, amplitudes de pilonnement et de tangage) sur cette efficacité.

La flottation industrielle est un procédé continu qui se déroule souvent en plusieurs étapes et dans lequel des particules d'une espèce de minéral donnée (généralement celles d'intérêt), présentes en différentes tailles, rencontrent une grande quantité de bulles de gaz (normalement de l'air) pour produire des agrégats bulles-particules minérales, qui sont extraits du dispositif de flottation (colonne ou cellule) en tant que produit de valeur (concentré). Le contenu en bulles est décrit par les conditions hydrodynamiques régnant dans le réacteur de flottation. Celles-ci sont reconnues pour leur influence sur la performance de la flottation. Ce projet de recherche porte sur deux sujets majeurs. Le premier est l'analyse de l’impact des particules solides sur les variables hydrodynamiques et l’effet de ces variables hydrodynamiques sur la récupération d’eau au concentré. Pour ce faire, l'effet du solide sur la distribution de la taille des bulles et le taux de rétention de l’air, ainsi que la corrélation entre la distribution de taille des bulles et le taux de rétention de l’air dans une colonne de flottation ont été étudiés. L'effet du taux de rétention de l’air, de la dimension des bulles et du taux surfacique de bulles (Sb) sur la quantité d’eau extraite au concentré a ensuite été analysé. Le second sujet traite de l'utilisation des variables hydrodynamiques pour la modélisation de la cinétique du procédé de flottation selon distribution granulométrique des particules introduites. La surface inter-faciale de bulle (Ib) est introduite à cet égard comme une variable hydrodynamique fournissant plus d'informations sur la distribution de taille de bulle que le taux surfacique de bulles qui est plus couramment utilisé. De plus, la corrélation entre la constante cinétique, la taille des particules et certaines variables hydrodynamiques a été analysée en utilisant une projection de structures latentes (PSL). Les résultats indiquent que l'importance relative des variables hydrodynamiques pour la modélisation de la cinétique de flottation dépend de la distribution granulométrique des particules. Finalement, les variables hydrodynamiques suggérées pour chaque classe granulométrique considérée ont été utilisées pour produire des modèles de régression mono-variable de la constante cinétique. / Industrial flotation is a continuous and often multistage process, where particles of a given mineral species (usually the targeted one), present in different sizes, encounter a large amount of gas bubbles (normally air) to produce mineral–bubble aggregates, which are removed from the flotation device (cell or column) as a valuable product (concentrate). The bubble content inside the cell is characterized by the prevailing hydrodynamic conditions (known as gas dispersion variables), which in turn are known to influence the flotation performance. This research project deals with two major topics. The first one is identifying the effect of mineral particles on hydrodynamic variables, and the effects of hydrodynamic variables on the final water recovery. For this purpose, the effect of solid particles on the bubble size distribution and gas hold-up, as well as the correlation between bubble size distribution and gas hold-up in column flotation were studied. It is followed by an assessment of the effect of the gas hold-up, bubble size and bubble surface area flux (Sb) on the amount of water reporting to the concentrate. The second topic deals with applying appropriate hydrodynamic variables for flotation modeling based on a given introduced particle size distribution. The interfacial area of bubbles (Ib) is introduced to address this issue as a hydrodynamic variable providing more information about the size distribution of bubbles than the commonly used bubble surface area flux. The correlation between the flotation rate constant and particle size as well as given hydrodynamic variables using a Projection to Latent Structures (PLS) has been analyzed. Results suggest that the relative importance of hydrodynamic variables for flotation rate modeling depends on the particle size distribution. Finally the suggested hydrodynamic variables for each of the various particle size-classes considered were used to produce single variable models for the flotation rate constant.

Ce mémoire présente l’étude d’un prototype novateur d’hydrolienne basé sur l’utilisation d’ailes oscillantes. L’Hydrolienne à Ailes Oscillantes (HAO) est en développement depuis une dizaine d’années à l’Université Laval et le potentiel de cette technologie a d’ailleurs été vérifié numériquement et expérimentalement. Il est maintenant nécessaire de développer des outils permettant de prédire le comportement de l’hydrolienne lorsqu’installée en rivière ou en courant de marée. Pour ce faire, la dynamique des fluides numérique (CFD) est utilisée afin d’étudier l’impact de différents paramètres sur les performances de l’hydrolienne. L’étude présentée dans ce mémoire décrit notamment l’influence de la modélisation de la turbulence, de la présence de la structure de l’hydrolienne et des effets de confinement. Dans un premier temps, une étude sur l’aile oscillante comparant deux niveaux de modélisation de la turbulence est présentée. Cette étude a permis de montrer que malgré la présence de structures turbulentes plus fines dans le sillage de l’aile avec le modèle Scale-Adaptive Simulation, les signaux de forces instantanées ainsi que les paramètres moyens de performance sont très similaires à ceux obtenus avec le modèle Spalart-Allmaras qui est de fait utilisé pour simuler l’hydrolienne HAO complète. Ensuite, l’hydrolienne HAO constituée d’une paire d’ailes oscillantes à l’intérieur de sa structure est simulée. La technique de maillage par grilles superposées est utilisée afin de simuler le mouvement relatif des différents corps. Cette représentation de l’hydrolienne a permis d’étudier l’impact de la structure de celle-ci sur ses performances et ainsi d’optimiser sa forme afin de maximiser l’extraction d’énergie. En plus d’améliorer les performances, le carénage optimisé permet d’atténuer la sensibilité de l’hydrolienne à un écoulement amont désaligné. Le troisième principal aspect étudié est l’impact du confinement sur les performances de l’aile oscillante. En effet, les parois d’un canal d’essais, la topologie des fonds marins ainsi que la proximité de la surface de l’eau sont susceptibles d’avoir un impact sur les performances hydrodynamiques d’une hydrolienne. Les simulations réalisées sur une aile oscillante à différents niveaux de confinement ont montré que la puissance extraite augmente avec le niveau de blocage, mais en plus, que cette relation est linéaire pour un confinement inférieur à 40%. Finalement, une technique est suggérée afin de corréler les performances de l’aile oscillante dans différents environnements confinés. / This master’s thesis focuses on a novel prototype of hydrokinetic turbine based on oscillating foils. This concept known as HAO, which stands for “Hydrolienne à Ailes Oscillantes”, has been under development for about 10 years at Laval University and its potential in power extraction has been confirmed through numerical and experimental studies. Efforts are now focused on developing tools to predict the turbines behavior prior to its deployment in rivers or tidal streams. To achieve this goal, computational fluid dynamics (CFD) is used to investigate the impact of different parameters on the power-extraction performance of the HAO turbine. This study describes, among other things, the influence of the turbulence modeling, the presence of the frame structure and the blockage effects. First of all, a methodological study performed on a single oscillating foil is presented which compares two different turbulence modeling approaches. This work has shown that even if the Scale-Adaptive Simulation model presents finer structures in the wake of the foil, instantaneous forces and mean performance parameters closely match the results obtained with the Spalart-Allmaras model which is thus used to simulate the complete HAO hydrokinetic turbine prototype. In a second study, the HAO hydrokinetic turbine is simulated considering two hydrofoils oscillating within the frame structure. The overset mesh technique is used to represent the relative motions of the different bodies. This methodology allows to study the impact of the frame structure on the turbine performance and to optimize its shape in order to increase the power extracted. In addition to the enhanced performances, the optimized frame shape provides an improved robustness to misaligned upstream flows. The third principal aspect addressed in this thesis is the impact of flow confinement on the performance of oscillating-foils. In fact, towing tank walls, sea and river bed topology and free surface proximity are likely to have an impact on the turbine hydrodynamic performance. Simulations of a single oscillating foil for different blockage levels have shown that the power extracted increases with the blockage ratio, but more precisely that this relation is linear for confinement of less than 40%. Finally, a technique is suggested to correlate the performance of the oscillating-foils turbines in different confined environments.

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.