Ce projet de maîtrise s’insère dans l’effort de modélisation entourant le phénomène de focalisation extrême amorcé dans les groupes de recherche des professeurs Thibault et Piché. Le mémoire présente d’abord une comparaison entre le formalisme de Richards-Wolf étendu (ERWT) et les solutions numériques de propagation d’onde électromagnétique obtenues à l’aide d’un algorithme FDTD. Les résultats montrent que l’utilisation du formalisme étendu permet de traiter le processus de focalisation non paraxiale pour une gamme étendue de surfaces réfléchissantes avec peu de variations entre les deux techniques. Les limitations intrinsèques à l’algorithme utilisé et la divergence de la fonction d’illumination imparfaitement traitée pourraient expliquer les déviations observées. Ensuite, cherchant à améliorer le traitement de cette illumination, le processus de focalisation inverse a été développé. En effet, il est possible de formuler une inversion du formalisme de Richards-Wolf (RWT) pour des systèmes à symétrie de révolution en définissant l’illumination à polarisation radiale en fonction d’une distribution recherchée au foyer selon un seul axe (radial ou sur l’axe optique). En utilisant un seul axe, le problème de surdéfinir le champ ne se pose pas et un critère est fourni afin de conclure de la validité physique du patron d’illumination calculé. Avec la méthode proposée, il est souvent possible d’obtenir des solutions analytiques qui sont essentielles à l’obtention d’une meilleure compréhension des différences entre les modèles paraxiaux et non paraxiaux. De plus, il est facile d’adapter cette dernière afin d’obtenir des solutions numériques pour des problèmes plus complexes. Ainsi, les figures d’illumination calculées peuvent directement être utilisées et sont particulièrement utiles lorsque la dimension optimale de la tâche focale pour une application donnée est connue. Un article a été soumis basé sur ces travaux. Enfin, une généralisation culminant avec le traitement de surfaces asphériques avec aberrations est présentée en adaptant le formalisme étendu (ERWT). Cet ajout au formalisme permet d’envisager de modéliser des surfaces complexes et possiblement hors de portée des corrections seulement en phase usuellement utilisées dans la littérature. Cependant, ces ajouts au formalisme complexifient de manière appréciable les intégrales de diffraction de Richards-Wolf. Finalement, une démonstration de la procédure a été effectuée pour un miroir parabolique comportant un léger tilt. / This project is part of the modeling effort around extreme focusing phenomenon taking place in the research groups of professor Thibault and Piché. At first, this thesis shows the comparison between the modeling by the extended Richards-Wolf formalism (ERWT) and by FDTD simulations of the field propagation. Both methods used for computing the electromagnetic distribution resulting from the reflection on a mirror in a non paraxial setup show satisfactory agreement. In fact, the extended formalism is suited to accurately model a large spectrum of reflecting surfaces such as elliptic mirrors, for which the description is not possible with the classical formalism of Richards-Wolf. Some intrinsic limitations of the used FDTD algorithm and the divergence of the illumination could explained the observed variations. Then, looking to solve the illumination problem, the inversion formalism has been developed. The integral of the Richards-Wolf formalism (RWT) of an axisymetric optical system can be inverted to define the radially polarized illumination pattern as a function of the electromagnetic distribution at the focus over a given axis (radial or the optical axis). Using only one axis at the focus, the field distribution is not over defined and a criterion is given to check the physical validity of the obtained illumination pattern. The method gives numerical or, in some cases, analytic solutions that can be used to obtain the optimized focal pattern fora given application. The analytical solutions are relevant as they can intuitively show the differences between the paraxial and non paraxial regimes. An article has been submitted on this particular subject. Finally, this thesis describes the generalization of the extended formalism to cover aspheric surfaces with aberrations. The procedure gives the possibility to accurately model complex reflecting surfaces that would otherwise be out of reach of the formalism. However, the complexity of the formalism increases compared to the initial diffraction integral. The proposed technique is demonstrated with a slightly tilted parabolic mirror.

L’anticlinal de San Corneli qui se situe dans les Pyrénées centrales espagnoles, est une structure plissée et faillée en contexte d’avant-pays de chaînes de montagnes. Cette structure compressive s’est mise en place lors du développement de la ceinture de chevauchement au-dessus de la plaque ibérique au Méso-Cénozoïque. Géométriquement, l’anticlinal de San Corneli correspond à un pli de propagation (8 km x 10 km) développé au-dessus du chevauchement de Boixòls, soit la structure la plus frontale de la nappe de Boixòls. Des structures extensives pré-compression ont été reconnues au sein de l’anticlinal ; il s’agit d’un réseau de failles normales et de fractures précoces ayant joué lors du dépôt de sédiments dans des bassins extensifs au Mésozoïque. Notre problématique de recherche est de développer une méthodologie de modélisation 3D à partir de techniques existantes pour permettre de vérifier l’influence de ces failles normales sur le développement de l’anticlinal de San Corneli et leur effet sur la géométrie particulière de ce pli de propagation. D’autre part, dans cette région le couvert végétal est minimal et l’anticlinal de San Corneli affleure bien, permettant l’acquisition de nombreuses données de surface et facilitant la modélisation de la géométrie du pli. Afin d’avoir une meilleure compréhension des relations qui existent entre la géométrie de formation d’un pli et les failles précoces, la modélisation 3D s’avère indispensable. La méthodologie utilisée consiste, dans un premier temps, en une modélisation inverse 3D. Le modèle géologique volumique est restauré à l’aide du module Restauration 3D proposé par gOcad, basée sur l’équilibrage de coupe. Cette technique permet d’obtenir un modèle paléogéographique 3D, et ainsi d’y déduire l’emplacement et la géométrie des failles normales préexistantes au moment du dépôt des sédiments dans le bassin. Dans le but de tester différentes hypothèses concernant le mode de développement de ce pli-faille nous avons effectué une série de modélisations directes 3D à partir du modèle paléogéographique. La modélisation directe 3D consiste à revenir à l’état actuel de déformation du modèle. Nous avons pu ainsi mieux apprécier la chronologie des évènements qui ont affecté le pli et voir dans quelle mesure les failles précoces ont influencé la formation du pli de San Corneli. En ce sens, nous avons pu faire ressortir que l’orientation des failles préexistantes par rapport à la contrainte régionale compressive avait une importance sur le devenir de ces structures. / The San Corneli anticline located within the Spanish central Pyrenees, is a foreland fault-related fold. This compressive structure developed within the Boixòls thrust sheet during the emplacement of the fold and thrust belt above the Iberian plate during the Meso-Cenozoic. Geometrically, the San Corneli anticline corresponds to a fault propagation fold (8 km by 10 km). Pre-compression extensive structures have been recognized within the anticline. These syntectonic normal faults and fractures were formed during deposition of sediments in the mesozoic rifted basin. The aim of our research is to determine the influence of these normal faults on the development of the San Corneli anticline and their effect on the particular geometry of this fault propagation fold. We used 3D geometrical modeling to gain a better understanding of the relationship between the pre-existing structures and the geometry and kinematics of the fault propagation fold. Furthermore, in this area the vegetative cover is minimal and the San Corneli anticline is very well exposed giving access to numerous field data and facilitating surface modeling of the present geometry of the fold. Our methodology consists, as a first step, in reverse 3D modeling. The geological model volume is restored using the 3D Restoration plugin proposed by gOcad based on cross section balancing. This technique allows to obtain a 3D paleogeographic model, and thus to deduce the location and geometry of preexisting normal faults with depositing sediments in the basin. In order to test different assumptions about the mode of development of this fault propagation fold, we have conducted a series of 3D geometrical models directly from the paleogeographic model. The aim of 3D forward modeling is to revert to the current state of deformation of the model. We were able to better appreciate the chronology of events that have affected the fold and evaluate to what extent the early faults have influenced the kinematics of the San Corneli fold. In this sense, we were able to show that the orientation of pre-existing faults with respect to the regional compressive stress influences the way in which these faults will be reactivated, either in a vertical (reverse) or horizontal sense (strike-slip).

Ce mémoire a comme objectif de prédire les performances d’une turbine modèle de type bulbe à l’aide de simulations numériques RANS k-ε. Tout d’abord, par des simulations numériques, on valide les éléments de la méthodologie les plus susceptibles d’influencer la qualité des résultats. Par la suite, on obtient la colline de rendement numérique et on analyse l’effet du jeu de bout d’aube ainsi que des jeux au moyeu sur la prédiction de performance. Les résultats indiquent que la colline de rendement numérique se compare bien à la littérature portant sur les turbines bulbes. On montre, de plus, que tenir compte du jeu de bout d’aube dans la simulation numérique diminue de manière significative le rendement prédit. L’impact des jeux au moyeu est cependant négligeable sur la prédiction de performance de la turbine. Finalement, une modification de l’aspirateur est proposée pour la suite du projet entrepris par le consortium de recherche.

La parallélisation est une solution possible pour améliorer un simulateur séquentiel s'il devient trop lent dû à une surcharge de calculs et qu'on ne désire pas concevoir à neuf un nouveau simulateur parallèle. Ce mémoire présente la conception UML d'une architecture pour un simulateur parallèle flexible et extensible capable de gérer différents environnements de déploiement par une configuration au temps d'exécution. Ce projet a vu le jour dans le but d'améliorer les performances de l'architecture KARMA, une architecture séquentielle pour la simulation d'engagements d'armes, pour réussir à atteindre le temps-réel dans une simulation haute-fidélité avec matériel dans la boucle. L'approche retenue propose une architecture non commerciale et développable à faibles coûts. L'implémentation et les tests préliminaires ont été basés sur un logiciel source libre et portable sur plusieurs plates-formes nommé ACE.

Titre de l'écran-titre (visionné le 26 juillet 2023) / Ce mémoire présente la modélisation et la simulation dynamique 3D de l'effecteur d'un robot humanoïde à l'aide du logiciel Simcenter Amesim. Ici, l'intérêt de modélisation et simulation porte sur un effecteur rigide pouvant entrer en contact multipoint avec son environnement et tenant lieu de mains et de pieds d'un robot humanoïde. On passe d'abord en revue les outils natifs disponibles dans Simcenter Amesim pour la gestion des contacts qui peuvent être intermittents. Les outils diffèrent selon que le modèle est en 1D, 2D ou 3D. Les outils 3D sont ensuite utilisés pour gérer les contacts multipoints entre un modèle d'effecteur simplifié et un environnement composé de plusieurs objets. Les résultats de simulation sont concluants, mais on constate des difficultés à effectuer les simulations lorsque le nombre de paires de contact est élevé. Ensuite, une extension des outils pour gérer le contact entre deux arêtes quelconques, toujours en utilisant les outils de contact natifs 3D et la programmation graphique est présentée. La technique fonctionne bien et permet, par exemple, de simuler un contact quelconque entre un effecteur et un barreau d'échelle ou entre un effecteur et une marche d'escalier.

La présente thèse intitulée «Développement d'outils intégrés pour la fabrication virtuelle de produits en aluminium et en acier par hydroformage de tubes pour l'industrie aéronautique » porte essentiellement sur l'étude de comportement du matériau métallique en grandes déformations. En utilisant la réponse de comportement, nous avons construit des outils de fabrication virtuelle qui détectent et corrigent « en ligne» les défauts qui peuvent survenir en cours de simulation de différents procédés de mise en forme, notamment, le procédé d'hydroformage de tubes. Pour être en mesure de modéliser efficacement le mouvement nonlinéaire de la matière et contrôler sa cinématique et sa dynamique, nous avons développé une loi hyperélastoplastique (écrite dans un repère corotationnel) en utilisant la théorie de grandes déformations (Simo) et nous l'avons couplé avec différents modèles d'écrouissage nonlinéaire (Voce et Simo) et d'endommagement ductile (Lemaitre et Chaboche). Nous avons calibré la loi de comportement en utilisant la courbe contrainte-déformation vraie (obtenue à partir des essais expérimentaux (uniaxiaux ou biaxiaux)) par l'utilisation de la méthode statistique de moindres carrés. À l'aide d'un utilitaire matériau (VUMAT) disponible dans l'environnement du logiciel Abaqus/Explicit, nous l'avons implémenté dans le but de vérifier l'évolution des différentes variables spécifiques au comportement réel du matériau (le champ de déformation et des contraintes, l'écrouissage, l'endommagement, etc.). En ce qui concerne la construction des outils adaptatifs pour les différents procédés de mise en forme et le contrôle du comportement du matériau, nous avons couplé les utilitaires de l'usager VUMAT (comportement du matériau), VDLOAD ( chargement) et avec un algorithme adaptatif qui utilise les résultats de VUMAT pour évaluer le chargement requis pour éviter des défauts et le transférer à l'utilitaire VDLOAD qui applique le chargement. Le premier sous programme adaptatif a été développé pour le procédé de gonflement de tubes avec les extrémités fixées. Il a été conçu pour évaluer la pression interne nécessaire pour gonfler le tube jusqu'à ce qu'il atteigne la circonférence maximale. Nous avons proposé une relation de calcul de l'incrément de pression fictive dans chaque point matériel. Cette relation contient des informations liées à l'écoulement (la limite d'écoulement) et à l'évolution de la plasticité équivalente. Après avoir trouvé la valeur maximale de la pression fictive dans chaque point matériel (pour un incrément de temps donné), nous avons évalué finalement la valeur maximale à partir de tous les points d'intégrations. Nous utilisons ultérieurement cette valeur maximale pour réaliser la mise à jour de la pression à l'intérieur du tube, afin de mieux contrôler le déroulement de l'opération dans le prochain pas de temps. Dans le but de couvrir une large gamme de produits fabriqués par hydroformage de tubes, nous avons proposé deux sous-programmes de contrôle pour évaluer la pression fictive dans chaque point matériel du tube. Le premier sous programme utilise le contrôle circonférentiel du champ de contraintes tandis que le second est basé sur le contrôle en compression axiale. Afin de vérifier la validité et la précision de nos travaux de recherche, nous avons évalué la qualité de nos résultats numériques en les comparant à des résultats expérimentaux disponibles dans la littérature scientifique. Dans ce contexte, nous avons utilisé les cas d'hydroformage dans une matrice conique et dans une matrice en « Y ». Parallèlement, nous avons tenté d'explorer la possibilité de déterminer la formabilité de tubes en alliages en aluminium en utilisant les outils de fabrication virtuelle que nous avons développés pour le gonflement de tubes avec extrémités mobiles. Dans ce cadre, nous avons proposé une méthode originale de traçage du diagramme limite de formage en utilisant les données du test de gonflement de tubes (avec extrémités fixées) combinés avec un calcul analytique qui évalue la plastification matérielle causée par le déplacement axial du tube. En dernier lieu, nous avons exploré l'hydroformage en double cavité et l'hydroformage des pièces qui sont d'abord cintrées et ensuite hydroformées. Dans ce contexte, nous avons réalisé une étude sur l'évolution de l'épaisseur par rapport aux changements de la géométrie (les rayons de raccordement à la profondeur des cavités) et par rapport au changement dans les conditions du contact entre le tube et la matrice.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2010-2011 / Le présent mémoire a comme objectif principal de procéder à la simulation d'un aspirateur de turbine existant en utilisant le DES, soit une nouvelle approche de modélisation de la turbulence. Il est postulé que la meilleure représentation de l'écoulement instationnaire qui existe dans cette partie de la machine pourrait faire la lumière sur la discordance entre la courbe de rendement applatie mesurée en centrale et la courbe prédite par les simulations RANS. Afin de procéder à ces simulations avancées, l'approche DES est d'abord testée avec un cas d'écoulement tourbillonnant dans un diffuseur conique pour lequel des données expérimentales sont disponibles. Par la suite, le cas de l'aspirateur réel est traité et plusieurs comparaisons entre les approches RANS, URANS et DES sont dégagées. Les résultats démontrent l'importance d'établir une méthodologie rigoureuse pour appliquer l'approche DES à un écoulement tourbillonnant dans un diffuseur. Principalement, le phénomène néfaste propre au DES de la séparation induite par le maillage (GIS) se manifeste très facilement pour ce type d'écoulement. Il est montré que les modèles qui protègent la couche limite, nommément les modèles DDES SA et DDES SST, sont beaucoup plus polyvalents et mieux adaptés à ces écoulements. Aussi, l'analyse des conditions d'entrée a montré, tant pour le diffuseur que pour l'aspirateur, que le niveau de viscosité tourbillonnaire qui y est prescrit doit être significativement plus faible que celui imposé en RANS. Les résultats finaux montrent une claire supériorité de l'approche DDES sur l'approche URANS pour prédire les instationnarités présentes pour les cas étudiés. Plus particulièrement pour l'aspirateur, le DDES permet de faire une analyse du tourbillon de sortie de roue et de ses caractéristiques associées à la charge. Par rapport à la prédiction de la performance, l'utilisation du DDES ne modifie pas la forme de la courbe de rendement. Il semble que la prédiction de la performance de l'aspirateur soit majoritairement gouvernée par la modélisation RANS en proche paroi et que le plateau dans la courbe de rendement mesurée ne soit pas entièrement attribuable à un phénomène instationnaire prenant place au centre de l'aspirateur.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2013-2014. / Le présent mémoire a pour objectif d’améliorer la prédiction des phénomènes de perte dans les aspirateurs de turbines hydrauliques. Pour ce faire, l’écoulement dans un aspirateur caractérisé par une diminution abrupte du coefficient de récupération près du meilleur point de fonctionnement a été étudié. Une méthode avancée de modélisation de la turbulence, le DES, a été mise à l’essai, afin de déterminer les gains associés à une représentation plus fine des mouvements turbulents dans cette composante. Les requis méthodologiques liés à cette approche, notamment par rapport à la condition d’entrée, ont été explorés, dans le but de développer une meilleure expertise d’utilisation du DES. Il a été démontré que le profil de vitesse radiale imposé dans le plan d’entrée du domaine de calcul altère de façon significative l’écoulement en aval et les performances prédites. Avec le profil de vitesse radiale mesuré expérimentalement, l’allure de la courbe de performance de l’aspirateur a pu être assez bien reproduite avec l’approche de modélisation de la turbulence URANS. Toutefois, certains aspects des simulations ne concordent pas bien avec les mesures expérimentales, c’est notamment le cas de l’écart de débit entre les deux canaux de sortie. Il a été établi que les structures d’écoulement en rotation sous les aubes de la roue nécessitent une discrétisation spatiale et temporelle extrêmement fine pour éviter qu’elles se diffusent prématurément sous le plan d’entrée. Toutefois, au point d’opération considéré, leur influence sur les performances de l’aspirateur s’est avérée très faible. Les simulations DES et URANS de l’aspirateur où des conditions d’entrée axisymétriques ont été imposées ont prédit des performances similaires. Cependant, le DES permet de simuler une dynamique tourbillonnaire beaucoup plus riche, avec un maillage et un pas de temps similaire au URANS, tout en étant largement moins dépendant des quantités turbulentes modélisées imposées dans le plan d’entrée. / The work carried throughout this thesis has for objective to enhance losses predictions in hydraulic turbines draft tube. In order to acheive this, the flow in a draft tube charaterized by a sharp drop in the pressure recovery coefficient near the best efficiency point was studied. Detached Eddy Simulation (DES), an advanced turbulence modeling approach, was put to the test, in order to asses the gain attributable to a finer and more precise description of turbulent motions in this component. The numerical methods required associated to this approach, especially regarding the inlet boundary condition, were investigated. It was shown that the radial velocity profile specified at the inlet of the computational domain alters significantly the flow downstream and the predicted performance. With the measured radial velocity profile specified at the inlet of the draft tube, reasonnable agreement was found between URANS numerical results and experimental measurements of pressure recovery. However, some aspects of the numerical simulations does not agree well with experimental data. It is the case for flow imbalance between the two outlet bays. It was established that rotating flow structures underneath the runner blades require extremely fine grid and time step resolution to avoid their premature diffusion underneath the inlet plane. Nevertheless, at the studied operating point, their influence on draft tube performance was found to be very limited. DES and URANS simulations of the draft tube where axisymmetric inlet boundary conditions were imposed predicted similar pressure recovery. However, DES enables to simulate much more complex and rich turbulent motions, at a computational cost similar to the one of a URANS simulatation and with much less influence from the modeled turbulent quantities specified at the inlet plane.

Dans des conditions de climat humide et lorsque les sols sont peu perméables, les systèmes de drainage souterrains sont généralement utilisés pour contrôler le niveau de la nappe phréatique et améliorer la production agricole. Cependant, les drains souterrains modifient à la fois les voies d'écoulement hydrologique et les taux de transport des nutriments des terres cultivées vers les eaux de surface, pouvant détériorer la qualité des eaux souterraines et de surface. De plus, des macropores sont souvent présents dans les sols composés de till argileux, ce qui génère un flux d'eau rapide et riche en nutriments de la surface du sol vers les drains souterrains. Une approche rentable pour réduire le lessivage des nutriments provenant de l'agriculture consiste à imposer des restrictions uniquement dans les zones vulnérables à la contamination de l'eau. Ces zones peuvent être identifiées à l'aide de modèles hydrologiques distribués. Les résultats obtenus sur de petits bassins versants expérimentaux doivent être simplifiés pour être appliqués à des échelles plus grandes, généralement requises pour l'élaboration de politiques. L'objectif de cette étude était d'examiner les avancées et les limitations de l'inclusion des drains souterrains dans les modèles d'écoulement de surface et souterrain. Les objectifs spécifiques étaient de i) démontrer l'utilisation des estimations de conductivité électrique spécifique (CE), pour améliorer les simulations hydrologiques dans un champ drainé, ii) étudier l'efficacité d'un modèle hydrologique et de transport de soluté tridimensionnel pour simuler un test de traçage de bromure (Br) dans un champ drainé et iii) évaluer différents modèles conceptuels de drains souterrains et d'hétérogénéité du sol pour la simulation numérique du drainage dans un bassin versant agricole au Danemark. Les résultats suggèrent que la simulation de la profondeur de la nappe phréatique peut être améliorée par l'inclusion d'hétérogénéités basées sur des estimations de la CE. L'approche des seepage nodes était appropriée pour simuler les débits de drainage, cependant la précision des simulations était meilleure pour les modèles à l'échelle du terrain. À l'échelle du bassin versant, le fait de ne représenter que les drains principaux est approprié pour pouvoir utiliser des maillages plus grossiers et pour simuler le débit des cours d'eau et les faibles profondeurs des eaux de surface dans les zones drainées. Des résultats similaires ont été obtenus lorsque les seepage nodes ont été appliqués sur l'ensemble des zones agricoles, sans tenir compte de l'emplacement spécifique des drains souterrains. Cette dernière approche peut être appliquée lorsque les drains souterrains ne sont pas cartographiés, ce qui est généralement le cas. Une représentation simplifiée de l'hétérogénéité et de la macroporosité peut expliquer les différences entre ls valeurs observées et simulées des charges hydrauliques, débits de drainage et processus de transport de solutés. Les approches de modélisation étudiées dans cette thèse peuvent améliorer la représentation de la dynamique de l'écoulement souterrain et les simulations du transport de substances agrochimiques lessivées des champs cultivés, telles que le nitrate et phosphate. / Under humid climate conditions and for low-permeability soils, subsurface tile drains are usually employed to lower the water table and enhance agricultural production. However, tile drains alter both the hydrologic flow pathways in agricultural catchments and the rates of nutrient transport from cropland to surface water bodies, potentially impairing the groundwater and surface water quality. Furthermore, macropores are often present in clayey till soils, generating rapid and nutrient-rich water flow from the ground surface to the tile drains. A cost-effective approach to reduce nutrient leaching from agriculture is to impose restrictions only in vulnerable areas to water contamination, which can be identified using distributed hydrological models. Results on small experimental catchments need to be simplified for application on larger scales, usually required for policy-making purposes. The objective of this study was to investigate the outcomes and limitations of including tile drains in surface and subsurface flow models. Specific objectives were to i) demonstrate the use of electrical conductivity (EC) estimates to improve hydrological simulations in a tile-drained field, ii) investigate the efficiency of a three-dimensional hydrological and solute transport model to simulate a bromide (Br) tracer test in a tile-drained field and iii) assess different conceptual models for tile drains and soil heterogeneity for the numerical simulation of tile drainage in an agricultural catchment in Denmark. The results suggest that the simulation of the water table depth can be improved by the inclusion of heterogeneities based on EC estimates. The seepage nodes approach was suitable to simulate drain discharge, however the accuracy of the simulations was better for the field-scale models. At the catchment scale, representing only the main drains was suitable to reduce the mesh refinement and simulate stream flow and low surface water depths in drained areas. Similar results were obtained when seepage nodes were applied all over the agricultural areas, without considering the specific location of tile drains. The later approach can be applied when tile drains are not mapped, which is usually the case. The misrepresentation of heterogeneity and macroporosity may explain the differences between observed and simulated hydraulic heads, drain discharge and solute transport processes. The modeling approaches investigated in this dissertation can improve subsequent simulations of tile drainage and the transport and fate of leached agrochemicals such as nitrate or phosphate.

Les opérations courantes de la centrale nucléaire Gentilly-2 nécessitent l'utilisation d'une grande quantité d'eau pour le refroidissement. En effet, les activités de la centrale nucléaire utilisent environ trente mètres cubes d'eau brute du fleuve par seconde. De plus, il faut considérer le caloporteur et le modérateur. Ces deux systèmes utilisent de l'eau lourde (D₂O). La conception de ces systèmes prévoit des systèmes de récupération en cas de fuites qui peuvent provenir par exemple des garnitures des pompes. Les deux systèmes ont aussi des systèmes de purification pour maintenir le caloporteur et le modérateur à des niveaux très bas d'impuretés et limiter ainsi la formation de produits d'activation. Ces systèmes de purification ont des filtres pour les particules solides et des résines cationiques et anioniques pour éliminer les impuretés solubles. Une fois saturées, ces résines sont dédeutérisées et envoyées dans des fosses d'entreposage avec une certaine quantité d'eau légère (pour faciliter la manipulation des résines). Le surplus est envoyé périodiquement dans les réservoirs de rejets radioactifs qui sont vidangés au fleuve et mélangés aux trente mètres cubes d'eau rejetée au canal de rejet de la centrale nucléaire. Cette étude s'intéresse au devenir des contaminants à partir de leurs rejets dans le fleuve Saint-Laurent et plus particulièrement aux trois radionucléides les plus rejetés à l'effluent liquide par la centrale nucléaire Gentilly 2 c'est-à-dire le carbone-14, le tritium et le cobalt-60. Des modèles informatiques ont été choisis pour étudier tout le phénomène de dispersion, de dilution, de sédimentation et de concentration dans les espèces de ces trois principaux contaminants rejetés dans le fleuve Saint-Laurent. Ces modèles permettent de traiter un grand nombre de contaminants différents. De plus, à partir des résultats obtenus, il a été possible de calculer les doses aux organismes qui se situent dans le panache de rejet de la centrale. Ces calculs de doses peuvent êtres ensuite extrapolés pour un grand nombre de rejets différents. En parallèle avec ce projet, plusieurs échantillons ont été recueillis dans la région sur une espèce implantée artificiellement dans la région, l'Eliptio Complanata, pour valider les valeurs obtenues à l'aide des modèles.