Mesurer les poussées des glaces sur les barrages hydroélectriques nécessite la mesure des contraintes dans les couverts de glace (CdG). Ces derniers sont des couches de glace qui se forment en hiver sur des réservoirs d'eau des barrages dans des zones hivernales (exemple, Québec). Le CdG se dilate lors d'une hausse de température, et est déstabilisé par la variation du niveau d'eau (VNE) dans le réservoir. Si de plus, le CdG est en contact avec le mur du barrage, des forces naissent sur le mur. Pour mesurer ces forces, dans un cadre de sécurité des barrages, des panneaux capteur de pression sont placés sur le mur du barrage, et des jauges sont positionnées dans le CdG. Parmi les jauges figurent les jauges biaxiales (JB) de 12 pouces (JB12) ou 30.48 cm de hauteur, et celles de 4 pouces (JB4) ou 10.16cm de hauteur. La petite hauteur de la JB4 lui permet de mesurer les contraintes selon la profondeur à laquelle elle est placée. Vu sa forte rigidité, la JB4 attire les forces sur le plan horizontal et sur le plan vertical. Une théorie mathématique permet de traiter le problème d'attraction des forces à l'horizontale (Savin, 1961), mais le problème d'attraction des forces en trois dimensions (effets 3D du CdG) n'est pas résolu. Or, les effets 3D causeraient le fait que la JB4 surestime les mesures. L'objectif de ce mémoire est de quantifier les effets 3D et évaluer l'impact du fluage de la glace, sur les contraintes estimées par la JB4 et sur l'amplification des contraintes (AC) autour d'elle. Pour se faire, sont développés, à l'aide du logiciel ANSYS, deux modèles numériques basés sur la méthode des éléments finis. Le premier modèle, en 2D, permet de valider l'approche numérique et évaluer l'impact du fluage. Le modèle 3D, sert à quantifier les effets 3D. Glace et JB4 sont homogènes et isotropes dans tout le projet. Le modèle 2D de base s'est avéré robuste, car les champs de contrainte et de déplacement ont au plus 0.5 % d'écart par rapport à la théorie de Savin. À partir du modèle robuste, est développé un modèle 2D (véritable outil) pour étudier l'impact du fluage sur PAC autour de la JB4 et sur la mesure de la JB4. Les résultats qui en découlent sont indicatifs en raison de la limitation au fluage secondaire. Un modèle de fluage plus complet permettrait de valider l'impact du fluage secondaire. Cet impact s'est avéré plus important pour des contraintes plus élevées. Par exemple, pour un chargement uniaxial de 250kPa, la diminution horaire moyenne de la contrainte principale maximale est de 0.92kPa/h, contre 1.79kPa/h pour une charge de 500kPa. Les effets 3D ont été évalués (matériaux uniquement élastiques). L'AC en 3D est au maximum de 5% plus élevé que TAC en 2D. Cet impact maximal correspond aux épaisseurs du CdG supérieures à 30 cm. Mais, la force moyenne (suivant la hauteur) reçue par le capteur est difficile à évaluer. En effet, le champ de contrainte autour de la JB4 et dans la direction verticale est très complexe et doit être profondément examiné. Cette complexité est due à la brusque variation de l'épaisseur du capteur le long de la hauteur. Ainsi, le modèle 3D doit faire l'objet d'une validation.

Titre de l'écran-titre (visionné le 14 août 2023) / La maladie d'Alzheimer (MA) est la cause la plus commune de démence chez les personnes de 65 ans et plus. Qu'est-ce qui déclenche réellement cette maladie ? La réponse est probablement une interaction dynamique de diverses voies, y compris l'amyloïde-bêta, les protéines tau hyperphosphorylées et les cytokines inflammatoires, modulées par des facteurs de risque tels que la génétique et le sexe. Cependant, tester une telle hypothèse multifactorielle chez l'humain est logistiquement presque impossible. Les modèles mathématiques pourraient résoudre ce problème en offrant un moyen de démêler les relations causales. Dans le cadre de ce mémoire, nous avons construit un modèle mathématique décrivant l'évolution d'un cerveau normal vers un état pathologique. Il est composé de dix-neuf équations différentielles ordinaires utilisant des paramètres issus de la littérature. Nos variables incluent l'amyloïde-bêta, les protéines tau, les neurones, les astrocytes activés, les microglies et macrophages, et certaines cytokines. Nous obtenons des résultats pour ces variables pour le vieillissement entre 30 et 80 ans d'âge. En raison de son importance dans la MA, nous avons séparé l'accumulation d'amyloïde-bêta en plusieurs composantes, soit les monomères intracellulaires et extracellulaires, les oligomères et les plaques. L'effet du sexe et la présence de l'allèle APOE4 ont été pris en compte via des paramètres sélectionnés. Le modèle est également sensible à l'insuline, étant donné que le diabète est un facteur de risque bien connu de la MA. Notre modèle relie avec succès la plupart des variables par rapport à l'âge. Nous présentons aussi deux variations du modèle, soit sans la variation de la concentration d'insuline, ou avec une diminution du taux d'activation des microglies. Les résultats obtenus sont cohérents avec les observations expérimentales. Malgré cela, notre modèle est encore sujet à plusieurs améliorations et une validation avec des données expérimentales est à faire. / Alzheimer's disease (AD) is the most common cause of dementia in people aged 65 and over. What actually triggers this disease? The solution is likely a dynamic interplay of various pathways, including amyloid-beta, hyperphosphorylated tau proteins, and inflammatory cytokines, modulated by risk factors such as genetics and gender. However, testing such a multifactorial hypothesis is logistically almost impossible to achieve in humans in practice. Mathematical models could solve this problem by offering a way to untangle causal relationships. As part of this thesis, we have built a mathematical model presenting the evolution of a normal brain towards a pathological state. It is composed of nineteen ordinary differential equations using parameters taken from the literature. Our variables include amyloid-beta, tau proteins, neurons, activated astrocytes, microglia and macrophages, and some cytokines. We obtain results for these variables for aging between 30 and 80 years of age. Because of its importance in AD, we have separated amyloid-beta accumulation into several components, namely intracellular and extracellular monomers, oligomers and plaques. The effect of sex and the presence of the APOE4 allele were taken into account via selected parameters. The model is also insulin sensitive, as diabetes is a well-known risk factor for AD. Our model successfully links most of the variables with age. We also present two variations of the model, either without the variation in insulin concentration, or with a decrease in the rate of activation of microglia. The results obtained are consistent with the experimental observations. Despite this, our model is still subject to several improvements and a validation with experimental data is to be done.

La démocratisation des véhicules hybrides branchables ainsi que des véhicules purement électriques implique un surplus de demande sur les réseaux de distribution. Le Vehicle-to-Grid (V2G) ou le Vehicle-to-Building (V2B) visent à répondre à cette demande accrue en utilisant les véhicules non plus comme de simples charges pour le réseau électrique mais comme des acteurs effectuant des échanges bidirectionnels. Les travaux présentés dans ce mémoire montrent, avec des données réelles du campus de l’Université Laval, une modélisation de flottes de véhicules et l’application d’un modèle d’optimisation linéaire que le V2B peut permettre de réaliser des gains financiers partagés entre les acteurs tout en rechargeant efficacement les véhicules participants. / The democratization of plug-in hybrid electric vehicles along with purely electric vehicles causes an increased electric demand on the power grid. Vehicle-to-Grid (V2G) or Vehicle-to-Building (V2B) aim to bring an appropriate response to this increased demand, by not simply considering vehicles as loads for the grid but as actors making bidirectionnal exchanges. The works presented in this master’s thesis show, with real data on the Université Laval campus, a modelling of vehicle fleets and the application of a linear optimization model, that V2B can provide financial gain shared between the actors of the system, while charging the vehicles efficiently.

La grande densité de couple qui fait la réputation des machines à flux transverse (MFT) est attribuable à l’agencement particulier de ses circuits magnétiques et électriques. Cependant, cet agencement complique aussi l’étape de modélisation nécessaire à la conception de la machine. Des facteurs tels que le flux de fuite, la saturation et l’aspect tridimensionnel du flux expliquent en grande partie la difficulté à élaborer un modèle magnétique analytique précis pour la machine. Par ailleurs, l’emploi d’un modèle imprécis pour l’optimisation des caractéristiques de la machine peut conduire à des résultats erronés. Pour cette raison, la plupart des concepteurs fait appel à une modélisation numérique avec des outils de calcul des champs 3-D, coûteux en ressources informatiques, en temps de calcul et malheureusement peu adaptée au processus de conception itératif d’une machine électrique. La thèse défendue dans ce document propose une alternative : elle présente une approche de modélisation hybride, faisant appel à un modèle magnétique analytique assisté par le calcul des champs. L’approche est appliquée pour l’optimisation des performances d’une structure de MFT particulière : la MFT à griffes (ou "claw-pole") à stator hybride (MFTCPSH). Dans un premier temps, on montre comment des facteurs de corrections, dont les valeurs sont déterminées à partir de simulations par éléments finis 3-D, peuvent permettre de palier aux problèmes de précision d’un modèle analytique au cours de l’optimisation de la machine. Un exemple de maximisation du flux à vide d’une MFTCPSH démontre la rapidité, l’efficacité et la précision de l’approche proposée. Dans un second temps, le modèle hybride et la méthode d’optimisation sont développés pour maximiser le couple d’une MFTCPSH. En se basant sur le cahier des charges d’une application de traction électrique de type moteur-roue, il est montré comment le modèle et la méthode proposés permettent, en quelques heures seulement, de déterminer les dimensions géométriques optimales qui maximisent la densité de couple de la machine. / The high power and torque capabilities of Transverse Flux Machines (TFM) are well known. The particular arrangement of the TFM’s electric and magnetic circuits, which mainly explains its exceptional features, also greatly complicates the TFM design process. Building an accurate analytical magnetic model for a TFM is a hard task: factors as the leakage flux, magnetic saturation or 3-D flux paths greatly affect the model accuracy. As the use of an inaccurate model can lead to erroneous results, TFM designers extensively use numerical models, as 3-D finite element analysis (FEA). Despite their accuracy, such numerical tools are still time and resource consuming but also tricky to use for optimization purpose. This thesis proposes an alternative approach for the design of TFM: a hybrid model, based on analytical calculations, corrected with finite element results inside the optimization process. Such an approach is used in order to optimize the performances of a claw-pole TFM with hybrid stator (CTFMHS). This work first presents how FEA-derived correction factors can be inserted in the analytical model to improve the model accuracy, in the course of the design process. The principles and the effectiveness of the proposed approach are shown through an example of CTFMHS no-load flux maximization. The approach is then derived in order to maximize the torque density of the CTFMHS. Considering the specifications of an in-wheel motor application, it is shown how the developed approach can determine, in a few hours only, the optimal machine geometry maximizing its torque density.

Le calcul des surfaces libres a fait l'objet de plusieurs ouvrages. L'objectif de ces travaux est de calculer avec précision les interfaces entre fluides immiscibles en minimisant la diffusion numérique et les pertes de masse. Plusieurs techniques existent pour effectuer de telles simulations et dans ce travail, on présente une variante améliorée de la méthode des surfaces de niveau (« level-set method »). Cette méthode requiert la résolution d'une équation de transport pour l'interface, couplée à la résolution des équations de NavierStokes. Dans la méthode proposée, l'interface est représentée par la surface de niveau 0 de la fonction distance signée. Cette interface présente une zone de transition de largeur 2e des paramètres rhéologiques des deux fluides comme la viscosité et la densité. Ces paramètres sont régularisés à l'interface pour éviter de gérer des discontinuités. La tension superficielle caractérise la résistance à la déformation de l'interface et joue un rôle très important. Il est donc primordial de bien l'imposer comme terme source dans l'équation de Navier-Stokes, ce qui est réalisé par l'introduction d'un tenseur qui agit uniquement à l'interface. Pour remédier au problème de conservation de la masse, nous avons développé une technique qui consiste à modifier la surface de niveau 0 de manière à conserver le volume. Cette correction est calculée avec la méthode de la sécante. Toutes ces stratégies sont facilitées par l'introduction d'une méthode d'adaptation de maillage instationnaire permettant de concentrer les éléments au voisinage de l'interface, là où le besoin s'en fait sentir. Cela assure une imposition très précise des forces surfaciques ainsi que de la transition des caractéristiques rhéologiques des différents fluides. Ce remaillage adaptatif est basé sur un estimateur d'erreur hiérarchique en dimension 2 et sur un estimateur basé sur une métrique en dimension 3. La résolution des équations de Navier-Stokes est faite par une méthode directe en dimension 2 et par une méthode itérative en dimension 3. L'élément de Taylor-Hood et un schéma implicite de différences arrières d'ordre 2 sont utilisés pour les discrétisations en espace et en temps. On utilise une méthode SUPG pour les équations de transport et de régularisation de l'interface. Enfin, plusieurs problèmes sont traités afin de valider la méthode et de comparer nos résultats avec ceux obtenus par d'autres chercheurs et différentes méthodes. On soulève enfin un certain nombre de questions sur les méthodes employées dans la littérature.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2008-2009 / Les objectifs de ce mémoire peuvent se scinder en deux parties distinctes. On retrouve tout d'abord un objectif de nature méthodologique qui consiste à réduire les coûts de calcul d'une méthode lagrangienne par particules vortex. Le second objectif est pour sa part davantage relié aux phénomènes physiques et à l'optimisation d'un concept de turbine à aile oscillante. L'objectif méthodologique concerne une méthode vortex lagrangienne 2-D permettant de résoudre des écoulements externes, visqueux et incompressibles, autour de multiples corps en mouvement arbitraire. Avec une telle méthode, le contrôle de la relaxation de la résolution spatiale dans le sillage des corps est plus complexe que dans le cas d'une méthode eulérienne (méthode de grille), puisque la discrétisation des champs est réalisée via un ensemble de particules vortex. Étant donné le nombre sans cesse croissant de ces particules dans le sillage, les coûts de calcul deviennent alors de plus en plus importants. Afin de mieux contrôler cette augmentation du nombre de particules, une approche de redistribution par boîtes est donc proposée et permet à l'utilisateur de définir des boîtes de dimensions et de résolution variables. Les résultats montrent que dans plusieurs simulations, le nombre de particules peut être réduit jusqu'à 50% sans toutefois affecter significativement les prédictions des forces agissant sur les corps. L'atteinte de l'objectif méthodologique s'est avérée nécessaire pour réaliser l'optimisation du rendement d'un concept de turbine à aile oscillante. Cette optimisation a été étudiée via la modification des fonctions de mouvement de l'aile et l'utilisation d'un volet au bord de fuite. Les résultats montrent que la modification des fonctions de mouvement ne permet qu'une faible amélioration du rendement. Toutefois, le déploiement d'un volet permet quant à lui d'aller jusqu'à doubler le rendement pour certaines conditions d'opération.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2013-2014. / Ce mémoire présente la méthodologie utilisée et les résultats d’optimisation de champs d’échangeurs géothermiques verticaux. Le contenu est présenté en deux articles de journaux scientifiques. Le premier article traite d’une étude sur l’impact de la stratigraphie du sol sur l’évaluation de performance d’échangeur géothermique. L’étude dévoile que les résultats des simulations prédisent une consommation d’énergie équivalente du système de thermopompe géothermique, que le sol soit stratifié ou homogène avec une conductivité thermique moyenne. Le premier article confirme donc l’hypothèse selon laquelle un système de pompe à chaleur couplée au sol peut être modélisé en utilisant une conductivité thermique moyenne. Le deuxième article présente une méthode d’optimisation de la conception d’échangeurs géothermiques afin de minimiser le coût du système global. Les résultats montrent entre autres que le nombre de puits dans une grille d’échangeurs, de même que leur profondeur, sont des paramètres d’importance dans la variation du coût. / This master thesis presents the method used and the optimization results for a ground coupled heat exchangers borefield. The content is divided into two papers submitted for publication to scientific journals. The first paper reports a study on soil stratigraphy impact on performance evaluation for a geothermal heat pump. The study reveals that the results obtained from simulation predict the same energy consumption of heat pump system for stratified soil and for homogenous soil defined by an average thermal conductivity. It confirms the assumption that a ground coupled heat exchanger system can be modeled using an average value of thermal conductivity. The second paper presents an optimization method of geothermal borefield design to minimize the cost of the system. Results show that number of boreholes in the grid, and their depth, are two influential parameters having strong effect on total cost.

Cette thèse porte sur la planification des ventes et des opérations (S±&OP) dans une chaîne d'approvisionnements axée sur la demande. L'objectif de la S±&OP, dans un tel contexte, est de tirer profit de l'alignement de la demande des clients avec la capacité de la chaîne d'approvisionnement par la coordination de la planification des ventes, de la production, de la distribution et de l'approvisionnement. Un tel processus de planification exige une collaboration multifonctionnelle profonde ainsi que l'intégration de la planification. Le but étant d'anticiper l'impact des décisions de vente sur les performances de la chaîne logistique , alors que l'influence de la dynamique des marchés est prise en compte pour les décisions concernant la production, la distribution et l'approvisionnement. La recherche a été menée dans un environnement logistique manufacturier multi-site et multi-produit, avec un approvisionnement et des ventes régis par des contrats ou le marché. Cette thèse examine deux approches de S±&OP et fournit un support à la décision pour l'implantation de ces méthodes dans une chaîne logistique multi-site de fabrication sur commande. Dans cette thèse, une planification traditionnelle des ventes et de la production basée sur la S±feOP et une planification S±fcOP plus avancée de la chaîne logistique sont tout d'abord caractérisées. Dans le système de chaîne logistique manufacturière multi-site, nous définissons la S±&OP traditionnelle comme un système dans lequel la planification des ventes et de la production est effectuée conjointement et centralement, tandis que la planification de la distribution et de l'approvisionnement est effectuée séparément et localement à chaque emplacement. D'autre part, la S±fcOP avancée de la chaîne logistique consiste en la planification des ventes, de la production, de la distribution et de l'approvisionnement d'une chaîne d'approvisionnement effectuée conjointement et centralement. Basés sur cette classification, des modèles de programmation en nombres entiers et des modèles de simulation sur un horizon roulant sont développés, représentant, respectivement, les approches de S±&OP traditionnelle et avancée, et également, une planification découplée traditionnelle, dans laquelle la planification des ventes est effectuée centralement et la planification de la production, la distribution et l'approvisionnement est effectuée séparément et localement par les unités d'affaires. La validation des modèles et l'évaluation pré-implantation sont effectuées à l'aide d'un cas industriel réel utilisant les données d'une compagnie de panneaux de lamelles orientées. Les résultats obtenus démontrent que les deux méthodes de S±feOP (traditionnelle et avancée) offrent une performance significativement supérieure à celle de la planification découplée, avec des bénéfices prévus supérieurs de 3,5% et 4,5%, respectivement. Les résultats sont très sensibles aux conditions de marché. Lorsque les prix du marché descendent ou que la demande augmente, de plus grands bénéfices peuvent être réalisés. Dans le cadre de cette recherche, les décisions de vente impliquent des ventes régies par des contrats et le marché. Les décisions de contrat non optimales affectent non seulement les revenus, mais également la performance manufacturière et logistique et les décisions de contrats d'approvisionnement en matière première. Le grand défi est de concevoir et d'offrir les bonnes politiques de contrat aux bons clients de sorte que la satisfaction des clients soit garantie et que l'attribution de la capacité de la compagnie soit optimisée. Également, il faut choisir les bons contrats des bons fournisseurs, de sorte que les approvisionnements en matière première soient garantis et que les objectifs financiers de la compagnie soient atteints. Dans cette thèse, un modèle coordonné d'aide à la décision pour les contrats e développé afin de fournir une aide à l'intégration de la conception de contrats, de l'attribution de capacité et des décisions de contrats d'approvisionnement pour une chaîne logistique multi-site à trois niveaux. En utilisant la programmation stochastique à deux étapes avec recours, les incertitudes liées à l'environnement et au système sont anticipées et des décisions robustes peuvent être obtenues. Les résultats informatiques montrent que l'approche de modélisation proposée fournit des solutions de contrats plus réalistes et plus robustes, avec une performance prévue supérieure d'environ 12% aux solutions fournies par un modèle déterministe.

Le bris d'un barrage est un événement exceptionnel et rare. Toutefois, lorsqu'il se produit, il peut engendrer de nombreuses pertes aussi bien matérielles qu'humaines. La modélisation de la propagation d'ondes de crue intumescentes issues d'un bris de barrage a fait l'objet d'inombrables études durant les quarante dernières années. À ce jour, il existe de nombreux modèles commerciaux capables de prédire avec satisfaction la propagation d'une onde de crue dans un canal à surface libre. Cependant, ces modèles ne sont pas conçus pour décrire la propagation dynamique des ondes de crue dans le cas des conditions hivernales caractérisées par la présence d'un couvert de glace à la surface des rivières. Pour les analyses de sécurité ou les études d'impact, les gestionnaires d'ouvrages de retenues et d'autres aussi qui s'intéressent aux écoulements fluviaux ont donc besoin de nouveaux outils de prédiction mieux adaptés à ces conditions particulières. Cette thèse se situe dans le cadre des efforts menés actuellement pour palier ce manque. L'objectif spécifique de la présente recherche est d'étudier la propagation des ondes hydrodynamiques dans un canal recouvert d'un couvert de glace. Pour y arriver, l'analyse s'est faite en deux parties. Premièrement, une analyse statique de la réponse du couvert de glace aux vagues de différentes formes est présentée. Cette dernière, faite à l'aide de modèles numériques uni- et bi-dimendionnels, a pour but de formuler des critères simples de bris du couvert de glace en rivière. Dans la seconde partie, un nouveau modèle numérique couplé eau-glace intitulé HYDROBEAM est bâti pour simuler la propagation dynamique des ondes dans un canal recouvert de glace. Pour ce nouveau modèle couplé, les équations de Saint-Venant, sous leurs formes conservatives et en présence d'un couvert flottant, sont reformulées. Les équations de Saint-Venant sont ensuite couplées aux équations décrivant le mouvement d'un couvert de glace semi-rigide. Les résultats du modèle couplé démontrent que la présence du couvert a pour effet d'atténuer et de ralentir les ondes et que l'effet se fait ressentir surtout sur des ondes très courtes (typiquement inférieures à 100 m).

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2012-2013. / Le comportement des bétons fibrés à ultra-hautes performances (BFUP) dépend fortement de l’orientation des fibres vis-à-vis de la direction des sollicitations. L’orientation des fibres étant principalement due à la mise en oeuvre de la structure, la ductilité d’un ouvrage s’en retrouve fortement dépendant. Dans ce mémoire de maîtrise, un modèle numérique est développé pour capturer l’effet de l’orientation des fibres via leurs mécanismes d’extraction. Il s’agit d’un modèle d’endommagement micromécanique, utilisé dans le cadre de la mécanique de la rupture linéaire élastique permettant de modéliser la propagation des fissures, la dissipation d’énergie et donc la ductilité de la structure. Cette approche permet également de rendre le modèle indépendant de la finesse du maillage et de modéliser l’effet d’échelle. Dans ce travail, les BFUP sont composés de fibres courtes en acier dont le pourcentage volumique est inférieur à 3%. La résistance en compression du BFUP est si grande que le matériau est admis être purement linéaire élastique en compression. Le modèle est implanté dans le logiciel Code_Aster pour calibration, validation et application sur des essais de traction et sur des essais de flexion quatre points : le modèle développé reproduit très bien ces résultats expérimentaux. / The behavior of ultra-high performance fiber reinforced concrete (UHPFRC) strongly depends on fibers’ orientations relatively to the direction of stresses. Fibers’ orientations are mostly due to concrete casting of the structure and therefore the ductility of a structure is highly dependent of those orientations. In this MSc thesis, a numerical model is developed to capture the effect of fiber orientation considering the fiber pull-out mechanism. A micromechanical damage model based on linear elastic fracture mechanics to model crack propagation, energy dissipation and thus structure’s ductility. This approach also corrects spurious mesh sensitivity and captures size effect. In this work, UHPFRC are made of short steel fibers with fiber content inferior to 3%. The compressive strength of UHPFRC is so important that the material is assumed to be purely linear elastic in compression. The model is implemented in Code_Aster software for calibration, validation and application on tensile tests and four-point bending tests : the developed model can reproduce these experimental results.