Premièrement, nous allons explorer la modélisation des cavités diélectriques bidimensionnelles. Plus spécifiquement, nous allons développer différentes méthodes de modélisation valides pour des cavités diélectriques à géométrie et profil d’indice de réfraction arbitraires. Ce degré de liberté supplémentaire pourra être utilisé dans le design de microcavités pour des applications spécifiques. Un formalisme de diffusion permettra de définir les modes caractéristiques de ce type de structure et d’en calculer les résonances. Une analyse numérique des équations résultantes montrera que les méthodes intégrales sont possiblement meilleures que les méthodes différentielles. Deuxièmement, nous discuterons de la modélisation de structures radiatives. Nous utiliserons les méthodes développées dans la section précédente pour modéliser les propriétés lasers des microcavités bidimensionnelles prédites par la théorie SALT. Nous aborderons aussi la modélisation de fibres-antennes RF, plus particulièrement les câbles coaxiaux à perte radiative, dans le but d’intégrer des fonctionnalités radio dans un textile de manière transparente à l’utilisateur. / In this essay, we will develop different modelization techniques valid for bidimensional dielectric cavities having arbitrary geometries and refractive index profiles and provide a way to accurately compute the resonances of such structures. The refractive index thus becomes an additional design variable for dielectric cavities. A numerical analysis of of the underlying equations of the theory will reveal that perhaps it is best to forego differential equations in favour of integral ones for the scattering problem. In the second part, we will discuss the modelization of radiating structures. Using the formalism developed in the previous section, we will study the lasing properties of bidimensional cavities using the newly developed self-consistent ab initio laser theory (SALT). We will also touch on the modelization of the class of antenna known as leaky coax

Cette thèse porte sur la modélisation numérique des problèmes de changement de phase. Le but est de proposer un algorithme robuste, basé essentiellement sur la méthode CutFEM, pour résoudre le problème avec la plus grande précision possible. Dans notre démarche, la méthode level-set est employée pour une représentation implicite de l'interface entre différentes phases, tandis que l'imposition des conditions essentielles sur cette dernière se fait à l'aide de la variante symétrique de la méthode de Nitsche. Le premier résultat de la thèse porte sur le développement d'un post-traitement fiable, facile à implémenter, pour une meilleure estimation du flux à travers une frontière immergée. La solution proposée est inspirée de la méthode de l'intégrale de domaine, à la différence que nous introduisons une stabilisation de type Ghost Penalty pour éliminer la sensibilité de la qualité de l'approximation à l'interface par rapport aux nœuds du maillage. L'idée est validée à travers une série d'exemples 2D et 3D.Le deuxième résultat concerne le développement d'un solveur CutFEM pour la simulation des problèmes de Stefan à deux phases sans convection. Dans notre approche, nous nous servons du premier résultat pour une meilleure évaluation de la vitesse de l'interface. Quant à la stabilisation des formulations d'extension de la vitesse et du transport de la fonction level-set, nous avons opté pour la méthode CIP (Continuous Interior Penalty) qui a déjà fait preuve d'une grande efficacité pour des problèmes fortement dominés par la convection. Des exemples numériques 2D permettent de confirmer l'efficacité de notre algorithme puisque des taux de convergence optimaux sont obtenus. Une simulation 3D est également proposée pour démontrer la validité du solveur dans une configuration plus réaliste. Pour le dernier résultat, nous reconsidérons le problème de Stefan, mais cette fois en prenant en compte la convection naturelle dans la phase liquide. L'équation de conduction de la chaleur est alors modifiée par l'ajout d'un terme de convection. De plus, cette équation est couplée aux équations de Navier-Stokes, introduisant de fait des non-linéarités supplémentaires. Un schéma itératif basé sur la méthode de Newton est alors proposé pour une résolution efficace du système. Les autres étapes de résolution sont assez similaires à celles de l'algorithme développé pour les problèmes de Stefan sans convection. Une étude numérique approfondie du solveur est présentée. Premièrement, l'efficacité de la discrétisation en temps et en espace est étudiée pour des problèmes (relativement simples) avec des solutions manufacturées. Ensuite, pour un problème beaucoup plus complexe modélisant la solidification de l'eau, une comparaison des résultats obtenus numériquement avec les résultats expérimentaux connus dans la littérature est présentée. Enfin, d'autres problèmes de changement de phase encore plus complexes tels que la fusion du gallium dans une cavité rectangulaire, ou encore celle du n-octadécane dans une cavité cubique sont considérés, et une comparaison des résultats de simulation avec les données expérimentales est également présentée. / This thesis focuses on the numerical modeling of phase change problems. The aim is to proposea robust CutFEM-based algorithm to solve the problem with the highest possible accuracy.In our approach, the level-set technique is used for an implicit representation of the interface,whereas the imposition of the essential conditions on the interface is done using the symmetric variant of the Nitsche method.The first result of the thesis is the development of a reliable, easy-to-implement post-processingfor an accurate estimate of the flux across an interface. The proposed solution is inspired by the domain integral method, with the difference that we introduce Ghost Penalty stabilization to eliminate the sensitivity to the interface location relative to the mesh nodes. The accuracy of the method is validated through a series of two- and three-dimensional examples. The second result concerns the development of a CutFEM method for the numerical simulation of two-phase Stefan's problems without convection. In our approach, we use the idea of the first result for an accurate evaluation of the normal speed of the interface. Regarding the stabilization of the velocity extension and front transportation problems, we opted for the Continuous Interior Penalty method which has shown great effectiveness for convection-dominated problems. Two-dimensional numerical examples are used to confirm the effectiveness of our algorithm since optimal convergence rates are achieved. A three-dimensional simulation is also provided to demonstrate the validity of the solver in a more realistic configuration. For the last result, we reconsider the Stefan problem, but this time taking into account natural convection in the liquid phase. In this phase, the heat conduction equation is then modified by the addition of a convection term. Furthermore, this equation is coupled to the Navier-Stokes equations, there by introducing additional nonlinearities. An iterative scheme based on Newton's method is then proposed for an efficient solution of the system. The other steps of the proposed solver are quite similar to those of the solver developed for Stefan problems without convection. An in-depth numerical study of the solver is presented. First, the efficiency of discretization in time and space is studied for (relatively simple) problems with manufactured solutions. Then, for a much more complex problem modeling the solidification of water, acomparison of the results obtained numerically with the experimental results known in the literature is presented. Finally, other even more complex phase change problems such as the fusion of gallium in a rectangular cavity, or that of n-octadecane in a cubical test cellare considered, and a comparison of the simulation results with experimental data is alsoprovided.

Les propriétés physiques, mécaniques et chimiques des matières premières ont un effet majeur sur la qualité des anodes en carbone pour le procédé de production d’aluminium. Ce travail tente d’étudier la faisabilité de l’application de simulation de la Méthode des Élément Discrets (DEM) à la technologie de production d’anodes. L’effet de la forme des particules et de la distribution de leurs tailles sur la densité apparente vibrée (VBD) d’échantillons de coke sec est étudié. Les particules de coke sont numérisées en utilisant des techniques d’imagerie à deux et trois dimensions. Ces images donnent les formes et les aspects réels des particules qui sont utilisées pour les modèles DEM pour les tests VBD pour le codage de flux de particules (PFC). Le coefficient de friction interne des particules de coke est estimé par la méthode de mesure d’angle au repos. Les résultats ont montrés comme attendu, que la VBD des échantillons de coke est affectée par la forme et la distribution de taille des particules. Les simulations à deux dimensions ont confirmé qu’en général, les échantillons formés de particules de tailles poly-dispersées ont une VBD plus haute que ceux dont la taille des particules est mono-dispersée. De plus, la VBD des échantillons augmente lorsque la fraction de grosses particules augmente. Cependant, la présence de 10 % massique de particules fines est nécessaire pour remplir les pores entre les grosses particules. De même pour la simulation 3D, le modèle suit la tendance des données expérimentales montrant que dans une éprouvette de 2,9 cm de diamètre, l’augmentation de la quantité de particules de - 4+6 mesh (de 3,36 à 4,76 mm) engendre une augmentation de la VBD. En conclusion, un modèle DEM approprié est capable de prédire le réarrangement des particules et l’évolution de la densité pendant le processus de vibration. / Physical, mechanical and chemical properties of raw materials have considerable effects on quality of carbon anodes for aluminium smelting process. The present work attempts to investigate the feasibility of application of Discrete Element Method (DEM) simulations in anode production technology. Effects of coke particles shape and size distribution on vibrated bulk density (VBD) of dry coke samples are studied. Coke particles are digitized using two-dimensional and three-dimensional imaging techniques and real-shape particles are used in DEM models of VBD test in Particle Flow Code (PFC). Internal friction coefficient of coke particles were estimated by means of angle of repose tests. Results showed that, as expected, VBD of coke samples is affected by shape and size distribution of the particles. Two-dimensional simulations confirmed that in general, mixed-sized samples have higher VBD than mono-sized cokes and as the fraction of coarse particles increases vibrated bulk density increases. However, existence of 10 wt.% of fine particles to fill the pores between coarse particles is essential. For 3D simulations also, the model follows the trend of experimental data showing that in the container of 2.9 mm diameter, as the content of -4+6 mesh (3.36-4.76 mm) particles increase, VBD increases. It can be concluded that a well-tailored DEM model is capable of predicting the particle rearrangement and density evolution during the vibration process.

Cette recherche introduit un modèle de décision multiobjectif qui aide le gestionnaire de projets à allouer les ressources à son projet. Le modèle comporte trois axes de décision représentant les compromis à faire entre la durée du projet, le coût du projet et les quantités de ressources allouées. L'allocation des ressources est déterminée en considérant le coût de chaque ressource. Le modèle incorpore le processus d'assignation des ressources, une problématique typique de l'affectation des ressources en gestion projets, à celui de !'ordonnancement des activités. Pour le résoudre, il a fallu concevoir et expérimenter de nouveaux algorithmes. La recherche comporte trois volets où chacun des volets alimente le volet suivant. Le premier volet propose de nouvelles procédures pour minimiser la durée du projet sous contraintes de disponibilités variables de ressources, dans le cas de problèmes ayant des activités avec un seul mode de réalisation ou de problèmes ayant des activités avec plusieurs modes de réalisation. Ces procédures comprennent une nouvelle règle de priorité selon un schème que nous qualifions d'ensembliste et une adaptation de la méthode Tabou. Les résultats de l'expérimentation montrent l'efficacité de la nouvelle règle de priorité et de l'adaptation de la méthode Tabou. Le deuxième volet propose des procédures pour minimiser le coût de disponibilité et le coût d'assignation des ressources allouées sous contrainte d'une date d'échéance du projet. La procédure incorpore l'adaptation de la méthode Tabou élaborée au volet précédent. Les résultats de l'expérimentation montrent la bonne performance des procédures développées. Le troisième volet propose une approche interactive par laquelle le gestionnaire de projets évalue des choix d'allocation de ressources afin de parvenir à une solution de compromis entre les trois axes de décision. Pour y parvenir, on évalue des solutions dans une région délimitée par les buts du gestionnaire de projets au lieu d'effectuer une approximation de la frontière efficace. L'expérimentation montre la faisabilité de l'approche interactive ainsi que l'efficacité de l'adaptation proposée de la méthode Tabou. Ces trois volets permettent d'élaborer l'architecture d'un Système Interactif d'Aide à la Décision pour l'allocation multiobjective des ressources en gestion de projets.

Ce travail de recherche a pour objectif de caractériser les pièces de la roue et de la trompette d’une turbine bulbe à partir d’un modèle réel. Suite à la reconstruction du modèle numérisé de la roue, cette géométrie sert de base pour le développement d’une aube instrumentée pour des mesures de pression. Un appareil de mesure laser sans contact est utilisé conjointement avec un bras de mesure portable de précision afin de récupérer les géométries du modèle sous forme de nuages de points. Ils sont ensuite traités dans un modeleur 3D afin d’obtenir la géométrie numérisée. Des analyses de déviations sont effectuées pour valider les modèles informatiques, quantifier les déviations et le respect de la norme CEI60193 par la roue. Une aube de roue instrumentée avec capteurs de pression embarqués est finalement développée pour investiguer les fluctuations de pression sur la roue en régime permanent et transitoire. / This research work deals with the geometry recovery of the runner and the draft tube of a bulb turbine, based on the real scale model. Once the geometry recovery of the runner completed, this computerized geometry allows to design an instrumented blade with embedded pressure sensors. A non-contact laser scanner on a precision measuring arm is used to create precise and detailed point clouds of the runner and draft tube geometries. Then, these point clouds are processed in a computer-aided design (CAD) software, in order to achieve a precise computerized geometry of the model. Deviation analyses allow to quantify the deviations so as to validate the CAD geometries and the IEC60193 standard agreement for the runner. Finally, a new instrumented blade with thirty embedded pressure sensors is developed so as to investigate pressure fluctuations in the runner in constant and transient operating conditions.

Les neurones thalamiques de relai ont un rôle exclusif dans la transformation et de transfert de presque toute l'information sensorielle dans le cortex. L'intégration synaptique et la réponse électrophysiologique des neurones thalamiques de relai sont déterminées non seulement par l’état du réseau impliqué, mais ils sont également contrôlés par leurs propriétés intrinsèques tels les divers canaux ioniques voltage-dépendants ainsi que l’arborisation dendritique élaboré. Par conséquent, investiguer sur le profil complexe de morphologie dendritique et sur les propriétés dendritiques actives révèle des renseignements importants sur la fonction d'entrée-sortie de neurones thalamiques de relai. Dans cette étude, nous avons reconstruit huit neurones thalamocorticaux (TC) du noyau VPL de chat adulte. En se basant sur ces données morphologiques complètes, nous avons développé plusieurs modèles multicompartimentaux afin de trouver un rôle potentiellement important des arbres dendritiques des neurones de TC dans l'intégration synaptique et l’intégration neuronale. L'analyse des caractéristiques morphologiques des neurones TC accordent des valeurs précises à des paramètres géométriques semblables ou différents de ceux publiés antérieurement. En outre, cette analyse fait ressortir de tous nouveaux renseignements concernant le patron de connectivité entre les sections dendritiques telles que l'index de l'asymétrie et la longueur de parcours moyen (c'est-à-dire, les paramètres topologiques). Nous avons confirmé l’étendue des valeurs rapportée antérieurement pour plusieurs paramètres géométriques tels que la zone somatique (2956.24±918.89 m2), la longueur dendritique totale (168017.49±4364.64 m) et le nombre de sous-arbres (8.3±1.5) pour huit neurones TC. Cependant, contrairement aux données rapportées antérieurement, le patron de ramification dendritique (avec des cas de bifurcation 98 %) ne suit pas la règle de puissance de Rall 3/2 pour le ratio géométrique (GR), et la valeur moyenne de GR pour un signal de propagation est 2,5 fois plus grande que pour un signal rétropropagé. Nous avons également démontré une variabilité significative dans l'index de symétrie entre les sous-arbres de neurones TC, mais la longueur du parcours moyen n'a pas montré une grande variation à travers les ramifications dendritiques des différents neurones. Nous avons examiné la conséquence d’une distribution non-uniforme des canaux T le long de l'arbre dendritique sur la réponse électrophysiologique émergeante, soit le potentiel Ca 2+ à seuil bas (low-threshold calcium spike, LTS) des neurones TC. En appliquant l'hypothèse du «coût minimal métabolique», nous avons constaté que le neurone modélisé nécessite un nombre minimal de canaux-T pour générer un LTS, lorsque les canaux-T sont situés dans les dendrites proximales. Dans la prochaine étude, notre modèle informatique a illustré l'étendue d'une rétropropagation du potentiel d'action et de l'efficacité de la propagation vers des PPSEs générés aux branches dendritiques distales. Nous avons démontré que la propagation dendritique des signaux électriques est fortement contrôlée par les paramètres morphologiques comme illustré par les différents paliers de polarisation obtenus par un neurone à équidistance de soma pendant la propagation et la rétropropagation des signaux électriques. Nos résultats ont révélé que les propriétés géométriques (c.-à-d. diamètre, GR) ont un impact plus fort sur la propagation du signal électrique que les propriétés topologiques. Nous concluons que (1) la diversité dans les propriétés morphologiques entre les sous-arbres d'un seul neurone TC donne une capacité spécifique pour l'intégration synaptique et l’intégration neuronale des différents dendrites, (2) le paramètre géométrique d'un arbre dendritique fournissent une influence plus élevée sur le contrôle de l'efficacité synaptique et l'étendue du potentiel d'action rétropropagé que les propriétés topologiques, (3) neurones TC suivent le principe d’optimisation pour la distribution de la conductance voltage-dépendant sur les arbres dendritiques. / Thalamic relay neurons have an exclusive role in processing and transferring nearly all sensory information into the cortex. The synaptic integration and the electrophysiological response of thalamic relay neurons are determined not only by a state of the involved network, but they are also controlled by their intrinsic properties; such as diverse voltage-dependent ionic channels as well as by elaborated dendritic arborization. Therefore, investigating the complex pattern of dendritic morphology and dendritic active properties reveals important information on the input-output function of thalamic relay neurons. In this study, we reconstructed eight thalamocortical (TC) neurons from the VPL nucleus of adult cats. Based on these complete morphological data, we developed several multi-compartment models in order to find a potentially important role for dendritic trees of TC neurons in the synaptic integration and neuronal computation. The analysis of morphological features of TC neurons yield precise values of geometrical parameters either similar or different from those previously reported. In addition, this analysis extracted new information regarding the pattern of connectivity between dendritic sections such as asymmetry index and mean path length (i.e., topological parameters). We confirmed the same range of previously reported value for several geometric parameters such as the somatic area (2956.24±918.89 m2), the total dendritic length (168017.49±4364.64 m) and the number of subtrees (8.3±1.5) for eight TC neurons. However, contrary to previously reported data, the dendritic branching pattern (with 98% bifurcation cases) does not follow Rall’s 3/2 power rule for the geometrical ratio (GR), and the average GR value for a forward propagation signal was 2.5 times bigger than for a backward propagating signal. We also demonstrated a significant variability in the symmetry index between subtrees of TC neurons, but the mean path length did not show a large variation through the dendritic arborizations of different neurons. We examined the consequence of non-uniform distribution of T-channels along the dendritic tree on the prominent electrophysiological response, the low-threshold Ca2+ spike (LTS) of TC neurons. By applying the hypothesis of “minimizing metabolic cost”, we found that the modeled neuron needed a minimum number of T-channels to generate low-threshold Ca2+ spike (LTS), when T-channels were located in proximal dendrites. In the next study, our computational model illustrated the extent of an action potential back propagation and the efficacy of forward propagation of EPSPs arriving at the distal dendritic branches. We demonstrated that dendritic propagation of electrical signals is strongly controlled by morphological parameters as shown by different levels of polarization achieved by a neuron at equidistance from the soma during back and forward propagation of electrical signals. Our results revealed that geometrical properties (i.e. diameter, GR) have a stronger impact on the electrical signal propagation than topological properties. We conclude that (1) diversity in the morphological properties between subtrees of a single TC neuron lead to a specific ability for synaptic integration and neuronal computation of different dendrites, (2) geometrical parameter of a dendritic tree provide higher influence on the control of synaptic efficacy and the extent of the back propagating action potential than topological properties, (3) TC neurons follow the optimized principle for distribution of voltage-dependent conductance on dendritic trees.

Cette thèse traite du problème d’optimisation de formes dans le contexte de la mécanique des solides en contact. Le modèle physique considéré est celui de solides linéaires élastiques en petites déformations, en contact (glissant ou avec frottement de Tresca) avec un corps rigide. Les formulations mathématiques étudiées sont deux versions régularisées de l’inéquation variationnelle décrivant le système d’origine : la formulation pénalisée et la formulation Lagrangien augmenté. Comme ces deux formulations présentent des non différentiabilités, nous proposons une approche par dérivées directionnelles pour obtenir les dérivées de forme associées. En particulier, pour chacune des formulations, nous exprimons des conditions suffisantes pour que la solution soit dérivable par rapport à la forme. Ceci nous permet de construire un algorithme d’optimisation topologique de type gradient, s’appuyant sur les dérivées obtenues et une représentation des formes par des ensembles de niveau (level-sets). L’algorithme bénéficie en outre d’une technique de découpage de maillage qui permet d’obtenir une représentation explicite de la forme à chaque itération, et ainsi d’appliquer fortement les conditions aux limites sur la zone de contact. Après avoir détaillé les différentes étapes de la méthode, nous présentons des résultats numériques en deux et trois dimensions pour en tester la validité. / This thesis deals with shape optimization for contact mechanics. More specifically, the linear elasticity model is considered under the small deformations hypothesis, and the elastic bodyis assumed to be in contact (sliding or with Tresca friction) with a rigid foundation. The mathematical formulations studied are two regularized versions of the original variational inequality: the penalty formulation and the augmented Lagrangian formulation. In order to get the shape derivatives associated to those two non-differentiable formulations, we suggest an approach based on directional derivatives. Especially, we derive sufficient conditions for the solution to be shape differentiable. This allows to develop a gradient-based topology optimization algorithm, built on these derivatives and a level-set representation of shapes. The algorithm also benefits from a mesh-cutting technique, which gives an explicit representation of the shape at each iteration, and enables to apply the boundary conditions strongly on the contact zone. The different steps of the method are detailed. Then, to validate the approach,some numerical results on two-dimensional and three-dimensional benchmarks are presented.

Les systèmes réels subissant des perturbations par l’interaction avec leur environnement sont susceptibles d’être entraînés vers des transitions irréversibles de leur principal état d’activité. Avec la croissance de l’empreinte humaine mondiale sur les écosystèmes, la caractérisation de la résilience de ces systèmes complexes est un enjeu majeur du 21e siècle. Cette thèse s’intéresse aux systèmes complexes pour lesquels il existe un réseau d’interactions et où les composantes sont des variables dynamiques. L’étude de leur résilience exige la description de leurs états dynamiques qui peuvent avoir jusqu’à plusieurs milliers de dimensions. Cette thèse propose trois nouvelles méthodes permettant de faire des mesures de la dynamique en fonction de la structure du réseau. L’originalité de ce travail vient de la diversité des approches présentées pour traiter la résilience, en débutant avec des outils basés sur des modèles dynamiques définis et en terminant avec d’autres n’exploitant que des données récoltées. D’abord, une solution exacte à une dynamique de cascade (modèle de feu de forêt) est développée et accompagnée d’un algorithme optimisé. Comme sa portée pratique s’arrête aux petits réseaux, cette méthode signale les limitations d’une approche avec un grand nombre de dimensions. Ensuite, une méthode de réduction dimensionnelle est introduite pour établir les bifurcations dynamiques d’un système. Cette contribution renforce les fondements théoriques et élargit le domaine d’applications de méthodes existantes. Enfin, le problème de retracer l’origine structurelle d’une perturbation est traité au moyen de l’apprentissage automatique. La validité de l’outil est supportée par une analyse numérique sur des dynamiques de propagation, de populations d’espèces et de neurones. Les principaux résultats indiquent que de fines anomalies observées dans la dynamique d’un système peuvent être détectées et suffisent pour retracer la cause de la perturbation. L’analyse témoigne également du rôle que l’apprentissage automatique pourrait jouer dans l’étude de la résilience de systèmes réels. / Real complex systems are often driven by external perturbations toward irreversible transitions of their dynamical state. With the rise of the human footprint on ecosystems, these perturbations will likely become more persistent so that characterizing resilience of complex systems has become a major issue of the 21st century. This thesis focuses on complex systems that exhibit networked interactions where the components present dynamical states. Studying the resilience of these networks demands depicting their dynamical portraits which may feature thousands of dimensions. In this thesis, three contrasting methods are introduced for studying the dynamical properties as a function of the network structure. Apart from the methods themselves, the originality of the thesis lies in the wide vision of resilience analysis, opening with model-based approaches and concluding with data-driven tools. We begin by developing an exact solution to binary cascades on networks (forest fire type) and follow with an optimized algorithm. Because its practical range is restricted to small networks, this method highlights the limitations of using model-based and highly dimensional tools. Wethen introduce a dimension reduction method to predict dynamical bifurcations of networked systems. This contribution builds up on theoretical foundations and expands possible applications of existing frameworks. Finally, we examine the task of extracting the structural causesof perturbations using machine learning. The validity of the developed tool is supported by an extended numerical analysis of spreading, population, and neural dynamics. The results indicate that subtle dynamical anomalies may suffice to infer the causes of perturbations. It also shows the leading role that machine learning may have to play in the future of resilience of real complex systems.

Ce mémoire présente les travaux de recherche ayant trait à la modélisation hydrodynamique du ruisseau Pratt et à l’évaluation du risque lié aux inondations à Coaticook. Quatre chapitres composent ce document. Le contexte de la recherche y est d’abord présenté. L’état des connaissances concernant le cours d’eau complexe qu’est le ruisseau Pratt est synthétisé, notamment son historique d’inondations, ses aspects hydrologiques et hydrauliques ainsi que les solutions envisagées jusqu’ici pour régler les problématiques d’inondations. La construction d’un bassin de rétention et l’érection de murs au centre-ville sont favorisées. Par la suite, la méthodologie employée pour développer des modèles hydrodynamiques 1D et 2D sur HECRAS 5.0 est présentée. Les différentes simulations effectuées pour évaluer le risque lié aux inondations de même que la façon dont ces résultats ont été exploités sont ensuite détaillées. La démarche d’évaluation du dommage moyen annuel (DMA) à l’aide de l’aléa des inondations, de l’exposition des bâtiments et de leur vulnérabilité y est notamment décrite. Le troisième chapitre du mémoire contient les résultats de l’étude. Les étendues d’inondation de même que la cartographie des zones inondables pour trois situations et pour les récurrences de 2, 20 et 100 ans sont présentées. Des séries de cartes exploitant des indicateurs de danger potentiel se trouvent également dans cette section. Enfin, les résultats du calcul du DMA pour différents scénarios sont montrés. Selon ces résultats, le DMA lié aux inondations de la ville de Coaticook se trouverait actuellement entre 145 k $-et 228 k $-an. La construction d’un bassin de rétention pourrait toutefois abaisser ces valeurs de plus de 80 %. Finalement, différentes analyses sont exposées dans le dernier chapitre du mémoire. Un jugement sur la légitimité des résultats qui comprennent des incertitudes importantes est d’abord avancé. Ceux-ci sont estimés réalistes, malgré certaines ambiguïtés. L’écoulement du ruisseau Pratt lors d’inondation est ensuite apprécié à l’aide d’une caractérisation de la propagation des inondations et d’une analyse des résultats de modélisation. Les risques liés aux inondations et la pertinence des mesures de mitigations sont finalement analysés. Techniquement, les risques pourraient être grandement diminués à l’aide du bassin de rétention et de l’érection des murs. La première solution pourrait d’ailleurs être rentable si son coût de construction demeure inférieur à 2,4 M $ -voire 3,7 M $ -en considérant les impacts potentiels des changements climatiques sur le régime hydrologique du ruisseau Pratt. Des risques supplémentaires peuvent toutefois accompagner la mise en place de ses mesures. Ces derniers sont décrits à la fin du quatrième chapitre de ce mémoire.

La sauvagine canadienne bénéficie des initiatives de suivi de la faune parmi les plus anciennes et les plus vastes au monde. En particulier, initié en 1955 afin d’obtenir une estimation annuelle des effectifs continentaux et fixer les quotas de chasse, l’Inventaire des Populations Nicheuses et de l'Habitat de la Sauvagine couvre aujourd’hui plus de 3 millions de kilomètres carrés d’habitats, du nord des États-Unis à l’Arctique et de l’Atlantique à l’Alaska. Les dimensions hors-normes de ces données ont contribué à la bonne atteinte des objectifs de gestion des populations de sauvagine et à une meilleure connaissance de l’écologie des espèces. Cependant, les observations issues des inventaires sont spatialement-discrètes (segments aériens, quadrats d’hélicoptères ou points d’eau) et leur couverture géographique reste relativement faible une fois rapportée à la superficie du Canada (≈1% de 10 millions de kilomètres carrés). Or, pour être efficace, la planification de la conservation requiert l’accès à une information quantitative spatialement-continue sur la distribution et l’abondance des espèces. La valorisation de données d’occurrences d’espèces spatialement-discrètes en couches d’information géographique continue est en revanche une application commune des modèles de distribution d’espèces dont les prédictions peuvent être étendues à l’échelle à laquelle les covariables environnementales sont disponibles. Si de tels modèles existent déjà pour guider la conservation de la sauvagine canadienne, des lacunes identifiées dans les approches méthodologiques, la couverture géographique et les thématiques environnementales explorées par les tentatives existantes suggèrent que des améliorations ou des extensions importantes pourraient être apportées. C’est l’objectif général de cette thèse : développer une nouvelle génération de modèles prédictifs de la distribution et de l’abondance de la sauvagine au Canada afin d’assister la planification spatiale des mesures de conservation des espèces et de leurs habitats. Dans le premier chapitre, une revue de la littérature a été effectuée avec pour objectif d’identifier les covariables environnementales d’intérêt pour modéliser la sauvagine. La principale contribution de cette étude a été la création d’une base de données de 533 associations canard-habitat attribuables à 133 covariables qui aidera le développement des futurs modèles. Dans le second chapitre, l’objectif a été d’améliorer les modèles nationaux d’abondance de sauvagine existants en développant de nouveaux modèles permettant de considérer explicitement les variations spatiotemporelles et testant des associations sauvagine habitat inexplorées grâce à l’utilisation d’un jeu initial de 232 covariables candidates. Pour cela, nous avons élaboré une méthode combinant techniques d’apprentissage automatique, procédures de sélection des covariables et approches hiérarchiques bayésiennes. La principale contribution de cette étude a été la mise à disposition de cartes annuelles de l’abondance de 18 espèces de sauvagines pour la période 1990-2015. Dans le troisième chapitre, l’objectif a été d’évaluer le potentiel d’intégration des données d’inventaires aériens standardisés et de science citoyenne (eBird) afin de modéliser la distribution de la sauvagine à l’échelle de la forêt boréale de l’Ouest canadien. La principale contribution de ce chapitre a été la formalisation d’une approche de modélisation intégrée de la distribution des espèces utilisant un processus par points dans une représentation état-espace. Nous avons démontré la capacité de cette approche à efficacement combiner des jeux de données d’occurrence hétérogènes afin de bénéficier de la complémentarité de leurs observations et de leurs couvertures spatiales. Dans le quatrième chapitre, l’objectif a été d’évaluer les effets potentiels des changements climatiques sur la distribution et l’abondance de la sauvagine dans l’est du Canada. Les résultats ont mis en évidence que les changements climatiques pourraient avoir un effet positif sur l’abondance de 7/12 des espèces évaluées, alors que 5/12 pourraient décliner. Une contribution majeure de cette étude a été le développement d’indices spatialement-explicites de la compatibilité climatique future de chacune des 12 espèces. Dans leur ensemble, les résultats issus de cette thèse ont permis (i) d’approfondir et de synthétiser l’état des connaissances sur l’écologie de la sauvagine canadienne, (ii) de fournir de nouvelles données pour guider les mesures de conservation et (iii) de développer des méthodes innovantes et efficaces pour modéliser la distribution et l’abondance des espèces à large échelle. Les plus de 1000 cartes et couches de données raster rendues publiques constituent une contribution d’envergure majeure pour l'élaboration d'indicateurs de biodiversité, l'évaluation et l'exécution de stratégies de conservation des espèces et de leur habitat, ou encore le suivi des services écologiques. / Canadian waterfowl benefit from some of the longest and most spatially extensive wildlife monitoring initiatives worldwide. The annual Waterfowl Breeding Population and Habitat Survey in particular, initiated in 1955 for estimating annual continental population sizes and setting hunting quotas, now covers more than 3 million square kilometres of breeding habitats, from the northern United States to the Arctic and from the Atlantic coast to Alaska. The exceptional dimensions of these data have historically contributed to the successful achievement of waterfowl population management goals and to a better understanding of the ecology of these species. However, inventories’ records are spatially discrete (aerial segments, helicopter plots or water bodies) and their geographic coverage remains small when compared to the area of Canada (≈1% of 10 million square kilometres). In order to be effective, conservation planning requires spatially-continuous, quantitative information on species distribution and abundance. This is a common application of species distribution models which, by predicting the distribution and abundance of individuals from occurrence data and environmental covariates, can provide continuous geospatial information. While such models already exist for Canadian waterfowl, identified deficiencies in the methodological approaches, geographic coverage orenvironmental themes explored by existing attempts suggest that significant improvements or extensions could be made. This is the overall objective of this thesis: to develop a new generation of models to predict waterfowl distribution and abundance across Canada in order to assist the spatial planning of conservation measures for species and their habitats. In the first chapter, a literature review was conducted with the specific objective of identifying environmental covariates of interest for waterfowl modeling. The main contribution of this study was the creation of a database gathering 533 duck-habitat associations attributable to 133 covariates, which will assist in the development of future models. In the second chapter, the objective was to refine existing national waterfowl models by developing new, more interpretable models that explicitly account for spatiotemporal variations in abundances, while testing for unexplored waterfowl-habitat associations by using a set of 232 newly-available candidate covariates. For this purpose, we developed a method combining machine learning techniques, covariate selection procedures and hierarchical Bayesian approaches. The main contribution of this study was the provision of annual abundance maps of 18 waterfowl species for the period 1990-2015. In the third chapter, the objective was to assess the potential of integrating standardized inventory and citizen science data to model waterfowl distribution across the Canadian western boreal forest. The main contribution of this chapter was the formalisation of an integrated species distribution modelling approach, which was based on a state-space point process framework. We demonstrated the ability of this approach to efficiently combine heterogeneous occurrence datasets in order to benefit from the complementarity of their records and spatial coverages. In the fourth chapter, the objective was to assess the potential effects of climate change on the distribution and abundance of waterfowl in Eastern Canada. Results showed that climate change could have a positive effect on the abundance of 7/12 of the species evaluated, while 5/12 could decline. A major contribution of this study was the development of spatially-explicit future climate suitability indices for each of the 12 species. Overall, results from this thesis have allowed to (i) deepen and synthesize the state of knowledge on the ecology of Canadian waterfowl, (ii) provide new data to guide conservation measures, and (iii) develop innovative and efficient methods to model large-scale species distribution and abundance. The more than 1,000 maps and raster layers made publicly available constitute a major contribution for the development of biodiversity indicators, the evaluation and execution of conservation planning strategies, and ecosystem services monitoring.