Cette étude s'intéresse à une méthode vortex lagrangienne 2-D permettant de simuler des écoulements incompressibles externes autour de multiples corps en mouvement arbitraire. L'implantation initiale de la méthode a été développée à 1' Université Catholique de Louvain puis des modifications ont été apportées en collaboration avec le Laboratoire de Mécanique des Fluides Numérique pour prendre en compte le mouvement et la présence simultanée de plusieurs corps. Cependant, ces modifications n'étaient pas validées et aucune méthode de calcul des forces et du moment sur les corps individuellement n'était disponible. Au sein de cette étude, plusieurs approches de calcul des forces et du moment sont présentées et détaillées, toutes ayant la particularité de ne pas nécessiter la connaissance du champ de pression de l'écoulement. La première approche utilise des surfaces de contrôle englobant chaque corps ainsi que la surface même des corps. La seconde est développée à partir d'une surface de contrôle à l'infini et ne peut fournir que les forces nettes totales sur l'ensemble des corps en présence. La troisième approche ne requiert que de l'information (vitesse et vorticité) à la surface des corps et elle est la seule pour l'instant qui permet de fournir également le moment de force. Toutes ces approches ainsi que les modifications effectuées pour prendre en compte le mouvement sont validées pour des écoulements à bas Reynolds autour de corps en translation et rotation combinés. Les cas d'un cylindre unique, de deux cylindres en tandem et d'une aile oscillante sont étudiés. Les comparaisons sont effectuées à l'aide de résultats produits par le code commercial eulérien de volumes finis Fluent 6.2. Les résultats mettent en évidence l'habilité de la méthode vortex à capturer précisément la dynamique des fluides, ainsi que l'efficacité de chacune des méthodes de calcul des forces et du moment. / This work is about a 2-D lagrangian vortex method allowing simulation of external incompressible flow around arbitrary moving bodies. The initial implementation of the method was developed at Université Catholique de Louvain and adapted for multiple moving bodies in collaboration with Laboratoire de Mécanique des Fluides Numérique. However, these modifications were not validated and no method of force and moment calculation on individual body was available. In this study, several approaches of force and moment calculation are presented and detailed, ail of which do not require information on the pressure field. The first approach uses control surfaces surrounding each body as well as the body surfaces themselves. The second one is developed from a control surface taken at infinity and can only provide total net forces on ail bodies taken together. The third one requires information (velocity and vorticity) only at body surfaces. The latter is the only approach considered up to now that can provide moment of force on the bodies. Ail three methods as well as ail the modifications that were implemented to adapt the solver to moving bodies are validated for low Reynolds number flows around bodies undergoing combined translation and angular motion. Cases of a single cylinder, two cylinders in tandem, and oscillating airfoil are studied. Comparisons are then made with results obtained with the eulerian finite volume commercial code Fluent 6.2. Results show the ability of the vortex method to correctly capture fluid dynamics, as well as the efficiency of each method for the force and moment calculations.

L'objectif de cette étude est de mettre au point un modèle mathématique reliant la hauteur dominante et l'âge des peuplements dominés par l'épinette noire (Picea mariana (Mill.) B.S.P.) au Québec, ce qui inclut les pessières irrégulières. Ce modèle, qui tient compte de la structure diamétrale du peuplement, a été étalonné à partir des données provenant des inventaires du réseau de placettes échantillons permanentes, recueillies entre 1970 et 2003 dans la province de Québec. Plusieurs indices de régularité ont été testés selon l'approche méthodologique de Raulier et al. (2003). Le modèle hauteur-âge utilisant l'indice de Shannon comme indice de régularité s'est avéré le plus performant pour prédire la hauteur dominante de l'épinette noire. Pour un âge donné, le modèle prédit une hauteur dominante plus élevée dans le cas des peuplements irréguliers que dans celui des peuplements réguliers. Par ailleurs, l'asymptote de hauteur dominante des peuplements irréguliers est atteinte beaucoup plus tard que dans le cas des peuplements réguliers. Ces effets sont plus marqués pour les stations caractérisées par une forte valeur d'indice de qualité de station. Cette différence est probablement liée au comportement adopté par les arbres des peuplements irréguliers qui, au stade juvénile, auraient subi une période d'oppression. Une fois que ces arbres sont dominants, ils adoptent une croissance en hauteur semblable à celle des arbres plus jeunes. Finalement, l'équation hauteur-âge ainsi mise au point est indépendante de l'âge de référence, ce qui permet la prédiction de la hauteur dominante à n'importe quel âge désiré si la hauteur dominante à un autre âge est connue.

L'objectif de nos travaux est de faire la synthèse de mouvements plausibles et réalistes de marche humaine dans divers environnements de synthèse. Bien que la solution proposée puisse également s'appliquer aux autres mouvements de locomotion humains ou animaux, nos travaux traitent uniquement du problème de la marche humaine. Afin de résoudre ce problème, nous avons développé une approche permettant de générer une multitude de variations d'une animation issue de capture de mouvement. Ces variations sont obtenues en adaptant le mouvement original à un environnement de synthèse dont les paramètres, tels que l'inclinaison du sol ou la courbure de la trajectoire, sont variés. Nous sommes donc en mesure de produire un mouvement de marche courbe ou de marche sur un plan incliné à partir d'un mouvement de marche en ligne droite sur un sol horizontal, ce que nous qualifions d'extrapolation de mouvement. Une animation initiale, obtenue par capture de mouvement, est essentielle à la solution proposée. Adapter ce mouvement à un nouvel environnement de synthèse consiste essentiellement à ajuster les caractéristiques globales du mouvement, telles que l'orientation du personnage et sa vitesse de déplacement. Ce faisant, nous sommes en mesure de conserver les détails plus fins du mouvement qui lui confèrent son aspect humain, tels que le mouvement des bras ou la vitesse avec laquelle un pied entre en contact avec le sol. En conservant les détails fins du mouvement d'origine, la solution proposée assure un certain réalisme dans les mouvements synthétisés. Dans la solution proposée, l'adaptation du mouvement initial est basée sur le paradigme des contraintes spatio-temporelles, où la synthèse du mouvement est posée comme un problème d'optimisation numérique. En plus d'être une formulation élégante du problème, ce paradigme est tout indiqué pour faire la synthèse de mouvements physiquement plausibles. En combinant ce paradigme avec l'utilisation d'une animation initiale issue de capture de mouvement, nous sommes en mesure de produire des animations de mouvements humains plausibles et réalistes. En pratique, le problème d'optimisation sous-tendu par l'adaptation d'un mouvement par contraintes spatio-temporelles est fortement non linéaire et opère dans un espace à très grande dimensionnalité. Cette complexité peut fortement ralentir le processus d'optimisation et aller jusqu'à en empêcher la convergence. La solution proposée fait donc appel à plusieurs mécanismes afin de réduire cette complexité. Notons qu'aucun de ces mécanismes ne vient compromettre la polyvalence de l'approche, en limitant la complexité du modèle biomécanique du personnage par exemple. Parmi ces mécanismes, deux sont des contributions originales : une technique d'estimation rapide des forces de réaction du sol et une approche d'optimisation incrémentale. Ces deux mécanismes visent à simplifier le processus d'optimisation en fournissant une solution initiale très proche de la solution optimale. La technique d'estimation des forces de réaction du sol sert à donner à ces paramètres une valeur initiale qui est relativement proche de leur valeur optimale, ce qui simplifie significativement la tâche d'optimisation subséquente. Cette technique consiste à trouver, pour les phases de support double, les forces de réaction du sol minimisant l'effort interne du personnage. Ce problème peut être exprimé comme une séquence de sous-problèmes de programmation quadratiques. Cette formulation est un aspect central de notre contribution et elle permet d'atteindre la solution très efficacement. L'approche d'optimisation incrémentale proposée s'inspire des méthodes de continuation. Le mouvement original est considéré comme une solution, un mouvement optimal, pour l'environnement de capture. L'environnement de synthèse est ensuite modifié graduellement, en augmentant l'inclinaison du sol par petits incréments par exemple. À chaque incrément, un nouveau mouvement optimal est trouvé en utilisant la solution de l'incrément précédent comme point de départ. On procède de la sorte jusqu'à l'obtention du mouvement désiré pour l'environnement de synthèse considéré. Si les incréments sont suffisamment petits, la différence entre deux problèmes d'optimisation consécutifs sera petite et il en sera de même pour leur optimum respectif.

La thèse présentée a été effectuée en cotutelle entre l'Institut Polytechnique de Grenoble et l'université Laval à Québec. Les travaux ont impliqué les laboratoires universitaire GIPSA-lab à Grenoble et le Laboratoire de Vision et Systèmes Numériques de l'Université Laval (LVSN). Notre étude porte sur la segmentation des contours internes et externes des lèvres. L'objectif visé dans notre étude est de proposer un ensemble de méthodes permettant de modéliser précisément la zone de la bouche avec la meilleure robustesse possible. Par robustesse, notre [sic] entendons obtenir une méthode fiable qui ne nécessite pas de réglage de paramètres et qui permette une segmentation fidèle des contours externes et internes de la bouche. Dans un premier temps, une approche combinée région-contour est introduite dans le but d'obtenir une segmentation multi-locuteur de la bouche sur des images de visage en couleurs. Nous décrivons une approche par décorrélation permettant d'augmenter le contraste entre la peau et les lèvres sur des images en couleurs ainsi qu'une étude sur les gradients multi-échelles pour améliorer la robustesse de la modélisation des contours de la bouche. Ensuite nous présentons notre méthode de localisation et de segmentation région-contour de la bouche sur des images de visage en couleurs. Dans un second temps nous nous sommes intéressés à la détection de l'état de la bouche. L'étape de détection de l'état de la bouche est nécessaire à la modélisation de la région interne qui présente une grande variabilité de forme et de texture. Une approche bioinspirée, basée sur un modèle de rétine et de cortex visuel conduisant au calcul d'un spectre Log-polaire, a été développée pour modéliser la zone de la bouche. Ces spectres sont, ensuite, utilisés pour entraîner un réseau SVM destiné à identifier l'état de la bouche. D'autre part, de nombreux auteurs ont suggéré d'utiliser la modalité infrarouge en analyse faciale. Le LVSN possédant une solide expertise dans le domaine de la vision infrarouge, une étude sur la pertinence de la modalité infrarouge dans le cadre de la segmentation des lèvres est proposée. Après avoir localisé la bouche et identifié l'état de la bouche, nous nous intéressons alors à la segmentation des contours externes et internes de la bouche. Un modèle polynomial de contour externe, dont la complexité sera automatiquement adaptée en fonction de la bouche traitée, est présenté. L'aspect de la zone interne de la bouche pouvant varier très rapidement, l'extraction du contour interne est particulièrement difficile. Nous proposons de traiter ce problème par une méthode de classification non-supervisée pour sélectionner les régions internes de la bouche. La méthode de modélisation de contour par un modèle polynomial est par la suite appliquée pour extraire le contour interne de la bouche. Enfin, une analyse quantitative de la performance globale de l'extraction des contours internes et externes est réalisée par comparaison avec des vérités-terrain.

Ce mémoire développe un modèle d’équilibre général dynamique pour analyser les impacts des facteurs démographiques sur les soldes extérieurs des pays. Le modèle de Gertler [1999] est généralisé à une économie internationale regroupant deux zones économiques, l’une représentant les pays du Nord, développés et dont la population est vieillissante, et l’autre regroupant les pays du Sud, en développement et ayant une population jeune. Les résultats des simulations révèlent que le taux d’épargne d’un pays dépend de la structure d’âge de sa population. Les pays du Sud ont un taux d’épargne élevé alors que les pays du Nord ont un faible taux d’épargne. En conséquence, le modèle prédit des soldes externes excédentaires couplés à l’endettement pour les pays du Sud et des soldes externes déficitaires pour les pays du Nord. / This thesis develops a dynamic general equilibrium model to analyze the impacts of the demographic factors on the external balance of countries. The model of Gertler [1999] is generalized in an international economy including two economic zones, the first one representing the developed countries, with an ageing population, and the other one including the developing countries, having a young population. The results of the simulations reveal that the rate of savings of a country depends on the structure of age of his population. The developing countries have a high rate of savings while developed countries have a low rate of savings. As a consequence, the model predicts that developing countries have an external balance surplus coupled with the debts and developed countries have an external balance deficit.

Dans l’optique de réhabiliter les eaux souterraines d’un site contaminé par des hydrocarbures incluant le benzène, le toluène et le xylène (BTX), une méthode efficace qui est souvent utilisée afin de stimuler la biodégradation des BTX est l’injection de persulfate de sodium (Na2S2O8). Cette oxydation chimique in situ (ISCO) permet d’éliminer l’essentiel de la masse contaminée, par réaction avec le persulfate (S2O8 2-). Le reste de la masse est ensuite éliminée par augmentation de la biorestauration (EBR), puisque le BTX restant réagit avec le sulfate produit de la réaction d’oxydation chimique. Cette méthode de traitement séquentiel a notamment été utilisée lors d’un essai-terrain à Borden, en Ontario, dans le cadre d’un projet RDC-CRSNG entre 2010-2015. Sur le terrain, la méthode de suivi habituelle de l’efficacité d’une technique de réhabilitation est la mesure des concentrations. Toutefois, la diminution des concentrations peut être attribuée à plusieurs effets tels que la complexité d’une zone source (phase libre), l’écoulement transitoire de l’eau souterraine, les hétérogénéités du milieu poreux ainsi que la dispersion hydrodynamique. La séparation des processus responsables de la diminution de masse n’est pas facile à évaluer sur le terrain. Pour analyser les changements de zones oxydo-réductrices, l’analyse isotopique δ13C et δ2H – méthode précise pour situer les réactions – peut être utilisée. Afin de mieux comprendre le comportement des signatures sous conditions variables, et d’identifier les limitations de l’approche, la méthodologie doit être testée en conditions de contrôle. À Borden, ce fractionnement isotopique n’avait pas encore été simulé. Les simulations numériques à l’aide du modèle BIONAPL/3D sont réalisées afin de reproduire l’essai-terrain de Borden, l’objectif étant de mieux comprendre le cheminement de dégradation des hydrocarbures en eau souterraine. Les compositions isotopiques permettent de conclure quant aux différents processus responsables de la dégradation des BTX. Toutefois, ces compositions ne permettent pas de quantifier la masse dégradée. / An emerging strategy for the remediation of contaminated sites is the integration of different treatment technologies. One example of a synergistic sequential treatment system is to use persulfate, a strong chemical oxidant (ChemOx), to target the bulk of the contaminant mass in the high concentration zones and then allow the produced sulfate to enhance biodegradation of the remaining mass. The design and subsequent performance of this combined remedy depends on the development of an enhanced bio-remediation (EBR) zone from the initial ChemOx zone. To provide insight into the development of these two dynamic mass removal zones, the isotopic fractionation signatures of representative petroleum hydrocarbons were simulated using the BIONAPL/3D model, which includes groundwater flow and multi–component reactive transport. This model was used to simulate a pilot-scale experiment conducted at the Canadian Forces Base Borden where persulfate was injected into a controlled-release plume of dissolved BTX components (benzene, toluene and xylene). Isotopic signatures of δ13C and δ2H were used to identify redox reactions, and to provide insight into hydrocarbon degradation pathways including aerobic biodegradation, chemical oxidation by persulfate, and anaerobic biodegradation by microbial sulfate reduction. Isotopic fractionation of C and H was included by associating the fractionation factors with the maximum substrate utilization rates. Simulated isotopic signatures agree reasonably well with the observed isotope ratios, showing increasing shifts of δ13C and δ2H over time in the remaining dissolved BTX. Spatially distinct redox zones, which are also affected by advection and dispersion, were identifiable from the simulated isotope ratios. The simulation results show that mass loss was dominated by chemical oxidation followed by sulfate EBR. The modelling tool and approach will be useful for application at other sites to support the design of persulfate/EBR sequential treatment systems, and to investigate the role of engineering controls on system behavior.

Le but de ce travail était de modéliser le comportement de micropolluants en plus des polluants traditionnels dans une station d'épuration. Deux études de cas ont ainsi été réalisées. La première concerne une usine pilote de la station d'épuration de Lynetten à Copenhague, au Danemark. Un modèle tenant compte à la fois de la dégradation aérobie par une bactérie spécifique et de la sorption sur les matières en suspension a été développé afin d'y décrire le comportement de composés organiques xénobiotiques. Les expériences ont été effectuées avec un perturbateur endocrinien, le Bisphénol A. Les résultats démontrent que le modèle parvient à représenter adéquatement la dynamique d'enlèvement du Bisphénol A. Le processus utilisé pour évaluer l'adaptation de la biomasse afin de valider le modèle n'épouse pas parfaitement les données mesurées, mais il permet toutefois de prédire la structure générale du comportement du Bisphénol A. La seconde étude concerne la station de Norwich, au Royaume-Uni. Dans ce cas-ci, le modèle, qui considère la sorption des métaux lourds sur les matières en suspension comme processus réactionnel, a été évalué à partir de données recueillies à même la station. Les résultats démontrent la grande sensibilité du modèle face aux variations instantanées des concentrations de métaux lourds dans l'affluent. Il est néanmoins capable de prévoir l'évolution des concentrations des métaux lourds à l'intérieur de la station ainsi qu'à l'effluent.

Ce mémoire traite un problème de sélection des transporteurs à un niveau opérationnel de la hiérarchie de décision, d'un système de trafic constitué de trois flottes de transporteurs (interne, externe et spot), plusieurs entrepôts, plusieurs centres de distributions et plusieurs produits. Les demandes de produits sont considérées aléatoires et les transporteurs des flottes, interne et externe, sont sujets à des périodes de non-disponibilités. Pour permettre de suivre l'évolution du système dans le temps, la politique de sélection devrait inclure une stratégie d'affectation des quantités transférées en fonction de l'état du système. Les transporteurs externes proviennent d'un processus de sélection au niveau stratégique. Les transporteurs spots sont appelés du marché au comptant. Dans la première partie de ce mémoire, un cadre conceptuel est proposé permettant de raccorder les décisions de sélection prises au niveau stratégique et opérationnel. La sélection stratégique utilise un mécanisme d'enchères combinatoires dans lequel les transporteurs se mettent en concurrence en soumettant des mises en paquet sur les offres proposées par les expéditeurs. Un modèle mathématique mixte en nombres entiers est exposé pour déterminer les transporteurs gagnants de l'enchère avec lesquels l'expéditeur va s'engager sur une longue période. Les transporteurs gagnants de cette phase seront regroupés dans une seule flotte, appelée flotte externe. Dans la deuxième partie du projet, le problème de sélection des transporteurs, qui traite le cas d'une chaîne logistique simple, constituée d'un seul entrepôt, un seul centre de distribution, un seul produit et les trois flottes des transporteurs, est abordé. De ce fait, une formulation en programmation dynamique stochastique du problème est développée. Une résolution numérique est menée par la suite pour résoudre les équations d'optimalité obtenues. Pour illustrer l'utilité pratique et le comportement robuste de la politique développée, plusieurs analyses de sensibilité sont effectuées. Cette politique combine les deux aspects : la sélection des transporteurs et le contrôle des inventaires dans les centres de distribution. Par la suite, une approche expérimentale basée sur la simulation et des techniques d'optimisation statistiques est adoptée, en vue de valider, quantifier et optimiser la structure de la politique obtenue par modélisation mathématique. Une comparaison statistique est par la suite menée permettant de montrer la supériorité de la politique paramétrée développée par rapport à une autre politique basée principalement sur les coûts du transport comme critères de sélection. Finalement, une politique générale est proposée permettant de résoudre le problème de sélection des transporteurs à la phase opérationnelle de la hiérarchie. Cette politique est basée sur la structure globale de la politique développée du cas simple dans la deuxième partie du projet. Elle est appliquée dans un cadre plus représentatif de la réalité. Une comparaison statistique est par la suite menée montrant les gains que la politique générale proposée occasionnent par rapport à une autre politique générale visant un minimum de stock dans les centres de distribution et un coût de transport réduit. Plusieurs interprétations ont été dégagées à travers cette étude.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2012-2013. / Ce mémoire présente une étude numérique du phénomène d'oscillations auto-induites d'une aile rigide montée sur un support élastique. Ces oscillations ont été rapportées expérimentalement au Collège Militaire Royal du Canada par l'équipe du professeur Poirel. Ils ont montré que le phénomène a lieu dans une plage de nombres de Reynolds spécifique où la transition de la couche limite peut survenir : 5 x 10⁴ < Rec < 1.3 x 10⁵. Des oscillations en tangage seulement ainsi qu'en tangage et pilonnement ont été observées. Les oscillations en tangage seulement ont une amplitude d'environ 5 degrés et une fréquence aux alentours de 3 Hz. Les oscillations en tangage et pilonnement ont des amplitudes de tangage pouvant atteindre 65 degrés selon la rigidité structurale et des fréquences allant de 3 à 5 Hz. Le phénomène a été étudié ici par la mécanique des fluides numérique. Le code libre OpenFOAM utilisant la méthode des volumes finis a été utilisé pour simuler le problème aéroélastique. Dans le cas des oscillations en tangage, une très bonne comparaison entre les résultats numériques et expérimentaux a été obtenue. L'utilisation d'un modèle de transition a entraîné une amélioration par rapport aux simulations numériques réalisées dans le passé et a contribué à mieux élucider la physique en jeu. La séparation de la couche limite laminaire étant le mécanisme déclencheur du phénomène, ces oscillations sont appelées flottement de séparation laminaire. L'impact de 1' intensité turbulente de 1' écoulement sur les oscillations a été étudié et s'est révélé jouer un rôle très important: un haut niveau empêchant l'apparition des oscillations. Le caractère secondaire du rôle joué par les structures d'écoulement à haute fréquence a été démontré ainsi que les différents mécanismes de dissipation d'énergie en jeu. Les oscillations auto-induites en tangage et pilonnement combinés ont également été simulées. La comparaison entre les résultats numériques et expérimentaux n'est pas aussi bonne que dans le cas de oscillations en tangage, mais des tendances similaires sont tout de même observées. Lorsque la rigidité structurale en pilonnement est petite, des oscillations de faibles amplitudes en tangage et pilonnement sont obtenues, tel que dans le cas en tangage pure. Lorsque la rigidité structurale est grande, d' importantes amplitudes de tangage sont obtenues qui s'avèrent du même ordre de grandeur que celles observées en expérimental. Ces oscillations diffèrent du cas en tangage puisqu'elles sont caractérisées par un flottement de coalescence plutôt qu'un flottement de séparation laminaire.

Des données de terrain et des études de modélisation ont montré que la recharge d'eau de fonte sous les calottes glaciaires peut avoir un impact important et durable sur l'écoulement des eaux souterraines. En Estonie, au nord du Bassin Artésien Balte (BAB), ce mécanisme de recharge est invoqué pour expliquer la présence d’importants volumes d'eaux souterraines marquées par un signal isotopique et géochimique glaciaire caractéristique, étant donné que la région a connu plusieurs glaciations durant le Pléistocène et a été entièrement recouverte par la calotte Fennoscandienne au cours du Dernier Maximum Glaciaire (DMG), il y a 20000 ans environ. Cette étude vise à tester cette hypothèse à l’aide de simulations numériques. En premier lieu, une étude conceptuelle a été effectuée pour déterminer quels processus sous-glaciaires doivent être représentés dans un modèle numérique qui reproduise adéquatement les écoulements souterrains et le transport de solutés. Les processus suivants ont été étudiés: la recharge sous-glaciaire d'eau de fonte, la déformation poroélastique du milieu poreux sous le poids de la glace, l’isostasie, l’évolution du drainage en surface, le permafrost et les écoulements densitaires impliquant des eaux douces de fonte et des saumures profondes. Ces processus ont été simulés dans un modèle représentant un bassin sédimentaire conceptuel, au cours d'un épisode glaciaire suivi d'une période postglaciaire. Le transport de trois traceurs d’eau glaciaire a été simulé: δ18O, solides dissouts et âge de l’eau. Les résultats montrent que la simulation de la recharge sous-glaciaire avec une condition-limite de type 1 (Dirichlet) n'est pertinente que pour des flux de faible amplitude, ce qui pourrait être le cas sous des calottes glaciaires dont la base n’est que partiellement en fusion. La compression de la matrice rocheuse diminue les surpressions, qui apparaissent uniquement dans les couches à faible diffusivité hydraulique et épaisses. Si la recharge sous-glaciaire est faible, la compression de la matrice rocheuse peut entraîner des sous-pressions après le retrait de la calotte glaciaire. L’isostasie réduit considérablement l'infiltration d'eau de fonte et les écoulements d'eau souterraine. Sous la couche de pergélisol, l'écoulement des eaux souterraines est réduit en-dessous de la calotte glaciaire mais augmente en région périglaciaire. Tenir compte des variations de densité en lien avec la salinité diminue l'infiltration d'eau de fonte en profondeur. Cette étude montre que chaque processus sous-glaciaire est potentiellement important et devrait être pris en compte dans des modèles d’écoulement des eaux souterraines et de transport de solutés en milieu sous-glaciaire. Cependant, il est raisonnable de ne représenter que la recharge sous-glaciaire si les informations manquent pour décrire correctement les autres processus. Par conséquent, ce seul processus a été simulé pour reproduire les écoulements d'eau souterraine sous la calotte Fennoscandienne dans le BAB. Les simulations ont été réalisées dans deux modèles 2D verticaux, afin de vérifier si la recharge sous-glaciaire d’eau de fonte peut expliquer la distribution particulière de δ18O (un traceur d’eau de fonte) dans les eaux souterraines de la région. L’un recoupe l’Estonie, l’autre la Lettonie et les îles estoniennes dans le Golfe de Riga. L'écoulement des eaux souterraines est simulé durant 28000 ans, depuis le DGM jusqu’à aujourd’hui, de même que le transport de δ18O pour tracer l'eau de fonte et confronter les résultats des simulations avec les données de terrain. L'espace d’incertitude de certains paramètres a été exploré, comme l’intensité et la durée de la recharge sous-glaciaire, ainsi que la composition isotopique initiale de l'eau de fonte. Les simulations fournissent un ajustement satisfaisant entre les valeurs observées et calculées de δ18O, confirmant l’hypothèse que le BAB a subi une phase de recharge sous-glaciaire durant le DMG. Elles montrent que la recharge sous-glaciaire a créé une inversion de l'écoulement des eaux souterraines dans le bassin. L’eau de fonte a infiltré tous les aquifères, en particulier les aquifères non confinés. Après le retrait de la calotte Fennoscandienne, l'eau de fonte a été entièrement remplacée par de l'eau météorique moderne, excepté dans les aquifères confinés où de l’eau de fonte a été préservée à proximité des zones de décharge. Par ailleurs, d’importants volumes d'eau de fonte sont probablement préservés sous la mer Baltique. Les simulations indiquent enfin que des épisodes de recharge sous-glaciaire antérieurs au DGM doivent être considérés afin d'expliquer les valeurs de δ18O dans la partie plus profonde du bassin. / Field evidence and modelling studies have shown that subglacial recharge of meltwater under wet-based ice-sheets can have a significant and long-lasting impact on groundwater flow. In the northern Baltic Artesian Basin (BAB), in Estonia, this mechanism of recharge is thought to be responsible of the presence of large volumes of groundwater with a characteristic glacial isotopic and geochemical signal, because the region experienced several glaciations during the Pleistocene and was entirely covered by the Fennoscandian ice-sheet during the Last Glacial Maximum (LGM), some 20 ky BP. The present study aims at testing this hypothesis by means of numerical simulations. First, a conceptual numerical study was performed to determine which glacial and subglacial processes need to be represented in numerical models for adequately capturing subglacial groundwater flow dynamics and solute transport. The relevance of the following processes was studied: subglacial recharge of meltwater, poroelastic deformation of the porous medium under ice-sheet loading, isostasy, evolution of surface drainage, permafrost, and density-dependent flow involving fresh glacial meltwater and deep brines. Simulations of these processes were conducted in a generic sedimentary basin during a single glacial event followed by a postglacial period. The transport of three common tracers of subglacial recharge was simulated: δ18O, TDS, and groundwater age. Results show that simulating subglacial recharge with a fixed flux boundary condition is relevant only for low fluxes, which could be the case under partially wet-based ice-sheets. Glacial loading decreases overpressures, which appear only in thick and low hydraulic diffusivity layers. If subglacial recharge is low, glacial loading can lead to underpressures after the retreat of the ice-sheet. Isostasy considerably reduces the infiltration of meltwater and the groundwater flow rates. Below permafrost, groundwater flow is reduced under the ice-sheet but is enhanced beyond the ice-sheet front. Accounting for salinity-dependent density reduces the infiltration of meltwater at depth. This study shows that each glacial process is potentially relevant in models of subglacial groundwater flow and solute transport. However, representing only subglacial recharge can be a reasonable assumption if information is missing to describe the other processes properly. Therefore, this single process is simulated to reproduce groundwater flow beneath the Fennoscandian ice-sheet in the northern BAB. Simulations are performed in two cross-sectional models, in order to check whether subglacial recharge of meltwater can explain the unusual distribution of δ18O in groundwater in the region, which serves as a tracer of glacial meltwater. One model crosses Estonia, the other crosses Latvia and Estonian islands in the Gulf of Riga. Groundwater flow is simulated over 28 ky, from the Last Glacial Maximum (LGM) to present-day, along with δ18O transport for tracing meltwater and to compare the results of the simulations with field data. Parameter space exploration of subglacial recharge conditions is used to tackle the uncertainty in the intensity and duration of subglacial recharge in the northern BAB, as well as in the isotopic composition of meltwater. Simulations provide a satisfying fit between the observed and the computed values of δ18O, supporting the idea that subglacial recharge happened in the northern BAB during the LGM. Simulations show that subglacial recharge created a flow reversal in the basin. Meltwater infiltrated into all aquifers, especially the shallow ones. After the retreat of the Fennoscandian ice-sheet, meltwater was entirely replaced by modern meteoric water, excepted in confined aquifers where some meltwater has been preserved close to the discharge areas. Large volumes of meltwater are also probably preserved beneath the Baltic Sea. Simulations also indicate that episodes of subglacial recharge prior to the LGM must be considered in order to explain the values of δ18O in the deeper basin.