Les moteurs à combustion interne utilisés pour propulser les véhicules routiers produisent des gaz à effet de serre, l’une des causes probables des changements climatiques actuellement observés partout sur la planète. L’utilisation de véhicules électriques représente une alternative intéressante à cette problématique. Grâce aux nouvelles générations de piles, il est désormais possible de recharger les véhicules électriques rapidement, ce qui peut favoriser leur adoption. Afin de réduire l’impact de cette recharge rapide sur le réseau électrique, une solution est d’emmagasiner de l’énergie lors de périodes creuses de consommation énergétique pour recharger les véhicules électrique rapidement lors des périodes de forte demande. La présente thèse fait état de travaux visant le développement d’un système de stockage inertiel d’énergie à moyen terme destiné à cette problématique des véhicules électriques. La conception de ce type de système nécessite l’étude de plusieurs domaines du génie mécanique et du génie électrique. Les recherches présentées dans cette thèse s’intéressent principalement à la dynamique des rotors, à la mécanique des matériaux et aux roulements magnétiques passifs. Plusieurs outils théoriques ont été développés pour prévoir le comportement dynamique du système et pour prédire les contraintes internes induites par la rotation. Plusieurs parmi ceux-ci ont été validés théoriquement grâce à des cas présentés dans la littérature, ainsi qu’expérimentalement à l’aide d’un système de stockage expérimental original. Ce dernier a aussi été utilisé pour valider plusieurs concepts théoriques de la dynamique des rotors qui sont peu discutés dans la littérature, c’est-à-dire l’influence des imperfections des roulements magnétiques passifs radiaux et du moteur à aimants permanents sur le comportement du système. Pour ce faire, des roulements magnétiques passifs radiaux réduisant les pertes mécaniques dans le système ont été conçus, fabriqués et caractérisés en utilisant une méthodologie novatrice. Il a été démontré expérimentalement qu’il est possible d’utiliser le principal défaut de ce type de roulements, c’est-à-dire leur instabilité, afin de réduire les pertes induites par une butée axiale à bille. Finalement, les premières étapes du développement d’un composite lourd utilisant l’acier comme renfort sont présentées. Mots-clés : Dynamique de rotor; Roulement magnétique passif; composite lourd; renfort d’acier; validation expérimentale; Moteur synchrone à aimants permanents; Stockage inertiel. / Internal combustion engines used to propel road vehicles produce greenhouse gases, one of the probable causes of global warming, a phenomenon currently observed everywhere on the planet. The use of electric vehicles is an interesting alternative to reduce the production of those gas emissions. New battery technologies enable rapid electric vehicles recharge, a feature which might increase their adoption. In order to reduce the impact of this rapid charging on the grid, an attractive solution is to store energy during periods of low consumption and then use it rapidly to recharge the onboard batteries during high demand periods. This thesis reports the work aimed at the development of a flywheel energy storage system that would solve this problem for electric vehicles. The design of such a system requires the study of multiple disciplines related to both mechanical and electrical engineering. The areas of research presented in this thesis are rotor dynamics, materials strength and behavior and passive magnetic bearings. Several theoretical tools were developed in order to predict the behavior of the system during operation and to estimate the internal stresses induced by the rotation. These tools were theoretically validated using reference cases presented in the literature and experimentally with the help of an actual storage system developed for this research. This system was also used to validate several rotor dynamics theoretical concepts that are rarely discussed in the literature, namely the influence of the imperfections in the passive magnetic bearings on the behavior of the system. To achieve this study, radial passive magnetic bearings were developed, fabricated and characterized using a novel technique. It has been demonstrated that it is possible to use the main flaw of this type of bearing, its instability, to reduce the mechanical losses induced by the mechanical thrust bearing. Finally, the first steps of the development of a new heavy composite that uses steel wires as reinforcement are presented. Keywords: Passive magnetic bearing; Rotor dynamics; Heavy composite; Steel reinforced plastics; Experimental validation; Permanent magnet synchronous motor; Flywheel energy storage.

À l’heure actuelle, les écoulements à deux fluides non miscibles jouent un rôle très important dans plusieurs domaines, que ça soit en science ou en ingénierie. Leur complexité est tellement élevée que les modèles actuels ne permettent de résoudre que des cas particuliers ou simplifiés avec un degré de précision qui demeurent souvent plutôt modeste. Une nouvelle approche numérique parait être une nécessité pour capturer la complexité physique du phénomène. Pour ce faire nous avons besoin d’outils robustes. Au niveau de l’interface de séparation entre les deux fluides non miscibles, les variables physiques sont discontinues, ce qui pose un défi majeur dans la description des variables et des conditions aux limites à l’interface. Le fait que les densités et les viscosités de chaque fluide soient différentes de part et d’autre de l’interface donne naissance à des défauts et des impuretés dans le champ des vitesses, ce qu’on appelle une discontinuité faible. Pour sa part, l’existence de la force de tension superficielle au niveau de l’interface crée une discontinuité sur le champ de pression, ce qu’on appelle une discontinuité forte. Un autre grand problème se pose au niveau de l’étude numérique du problème, où les méthodes numériques classiques ont une précision assez limitée dans ce genre de situation. L’objectif de ce travail est de fournir une étude complète de la dynamique de l’interface entre deux fluides non miscibles à l’aide d’outils mathématiques, physiques et numériques robustes. D’abord, une étude analytique du problème a été faite où l’équation de Navier-Stokes et les conditions de saut sur les variables physiques au niveau de l’interface de séparation entre les fluides ont été prouvées en détail. Pour traiter les discontinuités, nous avons discrétisé nos variables à l’aide de la méthode XFEM. Dû aux larges distorsions rencontrées dans ce genre d’écoulement, nous avons utilisé l’approche Eulérienne, pour corriger les oscillations des solutions dues aux choix du système de coordonnées nous avons utilisé les techniques de stabilisation SUPG/PSPG. Le traitement de la courbure des interfaces K été fait à l’aide de l’opérateur Laplace Beltrami et le suivi d’interface à l’aide de la méthode ¨Level-set¨. Pour le traitement des conditions de saut au niveau de l’interface la méthode Nitsche est développée dans différents contextes. Après avoir développé un modèle physique et mathématique dans les premières parties de notre travail, nous avons fait une étude numérique à l’aide de la plateforme de calcul FEniCS, qui est une plateforme de développement en langage C++ avec une interface Python. Un code de calcul a été développé dans le cas des écoulements de deux fluides non miscibles avec les modèles physiques et les outils mathématiques développés dans les sections précédentes. / The two-phase flow problems have an important role in the multitude of domains in science and engineering. Their complexity is so high that the actual models can solve only particular or simplified cases with a certain degree of precision. A new approach is a necessity to understand the evolution of new ideas and the physical complexity in this kind of flow, to contribute to the study of this field. A good study requires solid and robust tools to have performing results and a maximum of efficacy. At the interface of separation between the two immiscible fluids, the physical parameters are discontinuous, which gives us difficulties for the description of the physical variables at the interface and boundary conditions. The fact that the density and the viscosity are discontinuous at the interface creates kinks in the velocity, which represent a weak discontinuity. The existence of the surface tension at the interface create a discontinuity for the pressure field, it represents a strong discontinuity. The main objective of this work is to make a complete study based on strong and robust physical, mathematical and numerical tools. A strong combination, capable of capturing the physical aspect of the interface between the two fluids with a very good precision. Building such a robust, cost effective and accurate numerical model is challenging and requires lots of efforts and a multidisciplinary knowledge in mathematics, physics and computer science. First, an analytical study was made where the one fluid model of the Navier-Stokes equation was proved from Newton’s laws and jump conditions at the interface was proved and detailed analytically. To treat the problem of discontinuity, we used the XFEM method to discretize our discontinuous variables. Due to the large distortion encountered in this kind of fluid mechanic problems, we are going to use the Eulerian approach, and to correct the oscillation of solutions we will use the SUPG/PSPG stabilization technic. The treatment of the interface curvature k was done with the Laplace Beltrami operator and the interface tracking with the Level-set method. To treat the jump conditions with a very sharp precision we used the Nitsche’s method, developed in different cases. After building a strong mathematical and physical model in the first parts of our work, we did a numerical study using the FEniCS computational platform, which is a platform of computational development based on C++ with a Python interface. A numerical code was developed in this study, in the case of two-phase flow problem, based on the previous mathematical and physical models detailed in previous sections.

Suite au déluge qui a marqué la région du Saguenay en 1996, une couche de nouveaux sédiments a recouvert les sédiments du secteur amont du fjord du Saguenay. Avant le déluge cette zone était caractérisée par la présence de sédiments contaminés, associés au déversement incontrôlé d’effluents liquides, pratiqué dans le passé par les industries de la région. La couche de sédiments déposés en 1996 a enfoui les sédiments contaminés et constitue donc une barrière isolante composée de matériaux plus propres. Afin d'estimer l'efficacité à long terme de la couche, un nouveau modèle numérique (TRANSCAP-1D) a été développé. Ce modèle simule la migration de composantes dissoutes dans une colonne de sédiments et considère l'advection, la diffusion/dispersion, et l'effet de la bio-irrigation. La formulation mathématique représente un milieu à double porosité, composé de pores et de trous ou tubes de vers. Le modèle a été calé avec les profils de concentration de l'arsenic dissous qui ont été mesurés à deux stations du fjord du Saguenay, après la déposition de la couche de recouvrement naturelle. Par la suite, une analyse de sensibilité a été réalisée afin d'évaluer l'impact de certains paramètres caractérisés par une variabilité ou incertitude importante aux sites à l'étude. Les résultats montrent que les paramètres associés à la bio-irrigation ont un impact significatif sur la migration des contaminants dissous vers la colonne d'eau. Par la suite, nous avons réalisé une analyse d'incertitude en utilisant la méthode Monte Carlo et les résultats ont été intégrés dans une analyse de décision pour le design d'une couche de recouvrement potentielle. Le cas présenté est hypothétique et concerne la réhabilitation d'un site contaminé qui est fréquenté régulièrement par la population de bélugas du Saint Laurent. L'exemple illustre les avantages de l'application de l'analyse de décision pour trouver l'alternative correspondante au coût total plus bas, en considérant les coûts et les risques d'échec associés au projet de réhabilitation. / In 1996 two days of intense rainfalls caused severe flooding in the Saguenay region and a new sediment layer was deposited on the upstream area of the Saguenay fjord. In the past, this area was exposed to the uncontrolled discharge of industrial effluents and the sediments were contaminated. The new sediment layer buried the contaminated sediments and constitutes a natural barrier of cleaner material, which isolates the contaminants from the water column. A new numerical model (TRANSCAP-1D) was developed in order to estimate the effectiveness of the natural barrier in isolating the contaminated sediments from the overlying water column. The model simulates the migration of dissolved compounds in a sediment column and includes advection, diffusion/dispersion and the effect of bio-irrigation. The mathematical formulation represents a double porosity medium, composed of sediment pores and tubes or burrows dug by worms. The model was calibrated using the concentration profiles of dissolved arsenic measured at two stations of the Saguenay fjord, after the capping event. Thereafter, a sensitivity analysis was performed in order to evaluate the impact of certain parameters showing a great variability and uncertainty at the studied sites. The results indicate that the parameters associated to bio-irrigation have a significant impact on the migration of dissolved contaminants towards the water column. Thereafter, we performed an uncertainty analysis, using the Monte Carlo method. The results were integrated in a decision analysis for the design of a capping layer in a hypothetical case of sediment remediation. The presented case considers a contaminated site regularly frequented by the beluga population of the St. Lawrence Estuary. The example illustrates the advantages of the application of the decision analysis method, which is used to find the least cost option, considering the costs and failure risks associated to the remediation project.

Ce mémoire utilise une approche basée sur les méthodes d’identification paramétrique afin de représenter la relation entre signaux électromyographiques du biceps et du triceps versus le moment de force au coude variant dynamiquement lors de contractions à angle constant. La performance de deux modèles linéaires, le modèle «Auto Régressive with eXogenous inputs (ARX)» et le modèle «Output Error(OE)», de plusieurs ordres sont comparées. La performance des modèles est évaluée lors de l’utilisation de différents pré-processeurs, basée sur des canaux uniques ou multiples et l’utilisation d’un filtre «adaptive whitening». Le pré-processeur «multiple channels with adaptive whitening» améliore les performances d’identification. Le modèle OE permet une meilleure identification mais est sujet à des problèmes de convergence. / This thesis analyzes different parametric System identification (ID) methods to détermine a linear relationship between electromyographic signais from the biceps and the triceps and a dynamically varying torque about the elbow for constant angles contractions. The performance of two different linear models, Auto Régressive with eXogenous inputs (ARX) and Output Error (OE) of different orders are compared. The performance of the models are estimated with different pre-processors, single and multiple, with / without adaptive whitening filter. The multiple channels with adaptive whitening pre-processor increases the performance of the two models tested. The OE model leads to a better identification performance but is subject to convergence problems. / Québec Université Laval, Bibliothèque 2019

La géométrie et la flexibilité d'un canyon sont les paramètres qui affectent grandement la valeur des périodes naturelles dans les barrages en terre. Le canyon entourant des barrages peut être considéré comme un domaine illimité. Pour prendre en compte ces deux effets, le canyon a été modélisé par SBFEM et le barrage en terre, à géométrie limitée, par FEM. La technique hybride SBFEM-FEM pour l'analyse tridimensionnelle dynamique de l'interaction sol-barrage a été validée avec les résultats disponibles dans la littérature. Comme la matrice de rigidité dynamique du domaine non borné est complexe et dépendante de la fréquence, la méthode classique de superposition de modes n'est pas simple pour le système d'interaction sol-structure. Ainsi, pour obtenir la fréquence propre fondamentale, le barrage a été excité en direction amont-aval. Les périodes naturelles du barrage de terre pour des canyons de formes géométriques et de coefficient de impédance différents ont été obtenues. Ils se sont avérés avoir des effets significatifs sur la période naturelle. Les résultats ont été comparés aux données enregistrées réelles. Il a été constaté que les graphiques proposés dans cette étude peuvent être utilisés par des concepteurs de barrages pour l'estimation des périodes naturelles des barrages en terre dans des canyons de formes et de propriétés matérielles différentes. Plusieurs fonctions d'amplification correspondant à différentes conditions de canyon ont été obtenues en appliquant un déplacement uniforme à la limite du canyon. Une étude approfondie a été réalisée pour examiner les effets de la géométrie et de la flexibilité du canyon sur la réponse en régime permanent du barrage. Ces deux effets ont influencé de manière importante les fonctions d'amplification. Alors que la flexibilité du canyon affecte de manière significative la valeur de la fonction d'amplification maximale, cette valeur ne change pas pour les barrages en terre dans lesquels les canyons ont des formes différentes et la même longueur. De plus, la réponse latérale du barrage de terre dans le domaine temporel a été calculée pour analyser les effets susmentionnés lors d'un tremblement de terre réel. Les fonctions d'amplification proposées ont été utilisées pour comparer les spectres de réponse enregistrés du barrage d'El Infiernillo lors des tremblements de terre de 1966 avec la fonction d'amplification calculée. Un accord raisonnable a été observé entre eux. La méthode linéaire équivalente (EQL) a été implémentée dans le FEM. La technique FEMSBFEM a été étendue pour prendre en compte l'effet du comportement non linéaire des barrages en terre. Il a été observé que le comportement non linéaire affecte grandement la fréquence naturelle, la fonction d'amplification et l'accélération de crête maximale du barrage de terre situé dans les canyons. Les effets de la géométrie et de la flexibilité du canyon sur le comportement non linéaire ont été examinés, et on a vu qu'en augmentant la flexibilité du canyon, l'effet de la non-linéarité était diminué. Le barrage d'El Infiernillo a été modélisé par FEM-SBFEM non linéaire 3D, et une comparaison de la fonction d'amplification de crête obtenue par la méthode proposée avec les données enregistrées montre la précision du FEM-SBFEM non linéaire. / The canyon surrounding a dam can be assumed as an unbounded domain, and the geometry and flexibility of a canyon are parameters that greatly affect the values of natural periods in earth dams. In this thesis, in order to take into account these two effects, canyons are modeled by SBFEM, and earth dams, which have limited geometries, are modeled by FEM. The hybrid FEM-SBFEM technique used for the dynamic three-dimensional analysis of soil-earth dam interactions is validated with results available in the literature. Because the dynamic-stiffness matrix of the unbounded domain is complex and frequency-dependent, the classical mode-superposition method is not straightforward for a soil-structure interaction system, and thus, to obtain their fundamental natural frequencies, the modeled dams were excited in the upstream-downstream direction. The natural periods of earth dams in canyons with different geometries shapes and impedance ratios are obtained, and are found to have significant effects on the dams' natural periods. The results are compared with actual recorded data, and it is found that the graphs put forward in this study may be used by practical engineers for the estimation of natural periods of earth dams in canyons with different shapes and material properties. Several amplification functions corresponding to different canyon conditions are obtained by applying a uniform displacement at the canyons' boundaries. A comprehensive study is performed to examine the effects of canyon geometry and flexibility on the steady-state responses of the dams, and it is found that these two effects significantly influence the amplification functions. While the flexibility of the canyon does affect the maximum amplification function value, this value does not change for earth dams in canyons that have different shapes but the same length. In addition, the lateral responses of earth dams in the time domain are computed in order to analyze the aforementioned effects under an actual earthquake. The proposed amplification functions are used to compare the recorded response spectra of the El Infiernillo dam under the two 1966 earthquakes with the calculated amplification function, and a reasonable agreement is observed between them. The equivalent linear method (EQL) is implemented into the FEM, and the FEM-SBFEM technique is extended in order to take into consideration the effect of earth dams' nonlinear behavior. It is observed that such nonlinear behavior greatly affects the natural frequency, the amplification function, and peak crest acceleration of earth dams located in canyons. The effects of canyon geometry and flexibility on the nonlinear behavior are examined, and it is found that by increasing canyon flexibility, the effect of nonlinearity is decreased. The El Infiernillo dam is modeled by the 3D nonlinear FEM-SBFEM, and comparison of the crest amplification function obtained by the proposed method with the recorded data shows the accuracy of the nonlinear FEM-SBFEM.

Ce mémoire a pour objectif d'étudier les cycles électoraux au Canada. Le cadre théorique mobilisé pour cette recherche s'appuie sur deux théories de ce champ d'étude soit la théorie des élections de premier et de deuxième ordre ainsi que la théorie du contrepoids électoral (« electoral balancing »). Dans un système politique à multi-niveaux, ces deux théories postulent que les résultats électoraux des différents paliers électoraux sont en fait reliés entre eux. D'après ces théories, les élections nationales auraient plus d'importance pour les électeurs, ce pourquoi les résultats d'une élection régionale seraient en fait déterminés par des facteurs provenant à l'arène politique nationale. Au Canada, cela implique que les résultats d'une élection provinciale seraient alors influencés par l'élection fédérale précédente. Pour étudier ce phénomène, l'ensemble des résultats électoraux des élections générales fédérales et provinciales de 1988 jusqu'à aujourd'hui ont été compilé dans une base de données. Des régressions linéaires et polynomiales ont ensuite été réalisées en prenant les résultats fédéraux comme variable explicative et les résultats provinciaux comme variable expliquée dans un premier temps, puis en étudiant la relation inverse. Les résultats électoraux ont été divisés par province puis en cinq familles de partis politiques : libérale, conservatrice, travailliste, écologiste et régionaliste. Les résultats de cette recherche varient selon la province considérée : les élections fédérales ont parfois plus d'influence sur l'élection provinciale suivante dans certains cas alors que dans d'autres il s'agit de la relation inverse qui prévaut. Ce mémoire conclut également qu'il y a une dynamique de balancement entre le vote fédéral et provincial dans certaines provinces, mais il n'y a aucune tendance générale qui se dégage à l'échelle du Canada. Cette recherche apporte une contribution significative à la littérature sur les cycles électoraux au Canada. / This research project is studying electoral cycles in Canada. Two main theories from this field of research are mobilized in this paper. The first one is the theory of first and second order election and the second is the theory of electoral balancing. In a multi-level electoral system, these two theories suggest that electoral results from different electoral levels are related to each other. These theories assumed that national elections are more important to voters. Therefore, regional electoral results are explained by factors coming from the national electoral level. In Canada, this theory implies that provincial electoral results are influenced by the previous federal election. To study this phenomenon, all provincial and federal electoral results from 1988 until today were coded into a dataset. A first set of linear and polynomial regressions was performed using federal electoral results as the independent variable and provincial electoral results as the dependent variable. Then, another set of regressions was performed but to study the opposite relation. Electoral results were split by province and into five families of political parties: liberal, conservative, labor, ecologist, and regionalist. Results vary between provinces: federal elections have more influence on provincial results in some provinces while the opposite relation is more important in others. This paper also concludes that there is a dynamic of electoral balancing between federal and provincial results in certain provinces, but that no general tendency can be seen across Canada. This research project is a significant contribution to the literature on electoral cycles in Canada.

Les structures de commande prédictive (MPC) utilisant des simulateurs dynamiques comme modèles de procédé ne font pas légion. Cette rareté s'explique en grande partie par la difficulté de résolution des problèmes d'optimisation résultants. En effet, les algorithmes de programmation non linéaire ne sont pas toujours adaptés pour atteindre efficacement l'optimum des fonctions basées sur des modèles physiques, à plus forte raison encore si les équations utilisées sont inconnues (modèles de type boîte noire). C'est en s'appuyant sur ce constat qu'une nouvelle approche pour la MPC est proposée. Le principe consiste à substituer la minimisation explicite de la fonction objectif par une simulation du système en boucle fermée pour solutionner de façon généralement sous-optimale le problème de contrôle en boucle ouverte. Cette nouvelle méthode permet d'exploiter plusieurs des avantages de la commande prédictive sans être limitée par la complexité des modèles. Deux algorithmes sont présentés: décentralisé et découplé. Bien qu'avantageuse sur le plan du réglage, la structure découplée ne permet pas la même liberté que la structure décentralisée pour le choix de l'horizon de prédiction qui doit généralement être du même ordre de grandeur que le temps de réponse en boucle ouverte. Le développement d'une structure de simulation du comportement dynamique de la flottation en colonne, un procédé de séparation des minéraux, représente la seconde contribution de cette thèse. Une grande lacune des simulateurs proposés à ce jour demeure que même lorsque les équations différentielles de conservation sont utilisées, les variations dynamiques du niveau de pulpe ne sont jamais considérées. La structure présentée ici s'intéresse à la simulation des mouvements des phases présentes et de leurs effets sur le niveau de pulpe et sur les débits de sortie. Par ailleurs, comme bien d'autres procédés minéralurgiques, la flottation en colonne a peu bénéficié des avancées en contrôle de procédés. C'est donc sur ce terrain que les deux thèmes étudiés se rencontrent lors d'une mise à l'épreuve conjointe. L'étude de cas proposée s'intéresse à l'asservissement de trois variables d'opération critiques pour le bon fonctionnement du procédé à savoir, la concentration d'air dans la zone de pulpe, le flux net d'eau à l'interface et le niveau de pulpe à l'aide d'un contrôleur prédictif basé sur la simulation. / Applications of dynamic simulators for model predictive controllers design are rather scarce in the litterature. The complexity of solving the resulting optimization problems may explain this lack of popularity. In fact, nonlinear programming algorithms are not always well suited to efficiently reach the optimum of a fundamentaly-based cost function. The situation is even worse when the equations used in the model are unknown by the control designers (black box models). The simulation-based model predictive controller is an alternative formulation to perform model predictive control (MPC) without making use of any explicit optimization solver, but rather based on an easy-to-compute closed-loop simulation. The resulting scheme generally provides a sub-optimal solution and benefits from many interesting features of conventional MPC without being restricted by the model complexity. Two algorithms are proposed: decentralized and decoupled. The decentralized simulation structure allows a flexible setting of the prediction horizon (Hp) that is not possible in the decoupled case, easier to tune, but where Hp must generally be in the same order of magnitude that the system settling time. A second contribution of this thesis is the development of a framework for the dynamic simulation of a mineral separation process: column flotation. Until now, most of the proposed models or simulators were restricted to the steady-state behavior. When dynamic mass-balance equations were considered, a constant pulp level during the simulation was always assumed. The presented framework aims to simulate water, solids and gas motion and their effect on the pulp level and output flow rates. As it often happens in mineral processing, the column flotation process has not benefited from advanced control techniques. This is where the two previous subjects merge. The proposed simulation framework is used to design a simulation-based model predictive controller for process variables having a strong influence on metallurgical results (grade and recovery). A case study is presented where the pulp level, bias and air hold-up in the pulp zone are kept within an acceptable operating region.

Cette thèse considère des modèles de volatilité où la distribution conditionnelle des données est un cas particulier de la loi "Generalized Hyperbolic" de Barndorff-Nielsen (1977). Ces modèles permettent de capter les principales caractéristiques des séries financières à haute fréquence, à savoir le groupement de volatilité (volatility clustering), l'excès de kurtosis et de skewness ainsi que l'effet de levier qui s'applique au rendements des marchés boursiers. Etant donnée la forme fortement non linéaire de cette densité, nous utilisons l'approche Bayesienne basée sur les méthodes Markov Chain Monte Carlo pour l'estimation et l'inférence Cette approche est relativement simple à mettre en oeuvre et permet une inférence exacte et valable en échantillon fini ainsi que la comparaison de modèles qui ne sont pas forcément emboîtés. A titre illustratif, nous proposons des applications empiriques en employons des données journalières de l'indice boursier S&P500. D'abord, nous considérons un modèle de volatilité stochastique basé sur un mélange des lois normale et inverse-Gaussien où la variance conditionnelle est considérée comme un processus stochastique latent généré par la loi inverse-Gaussian. Conditionnellement à la volatilité, la loi des données est une normale. Il en résulte la loi normal inverse Gaussian (NIG) de Barndorff-Nielsen (1997) pour les données qui présente beaucoup de flexibilité pour capter les excès de kurtosis et de skewness. Dans ce modèle la volatilité est traitée de façon similaire aux paramètres du modèle et elle est simulée par l'échantillonneur de Gibbs. Ce modèle s'avère plus performant que les modèles GARCH asymétriques de Ding et al (1993). Par ailleurs, nous proposons les lois NIG de Barndorff-Nielsen (1997) et GH-skew student de de Barndorff-Nielsen et Shepard (2001) comme densités alternatives aux modèles GARCH asymétriques. Formellement, nous considérons deux modèles GARCH asymétriques à la Ding et al (1993), l'un avec une loi NIG et l'autre avec une loi GH-skew student. Dans ce contexte la volatilité est calculée de façon récursive sur la base de données passées. Les résultats sont quelque peu décevants pour la loi GH-skew student, puisque la performance de ce modèle est comparable à celle d'un modèle GARCH asymétrique standard

Les aspirateurs de turbines hydrauliques jouent un rôle crucial dans l’extraction de l’énergie disponible. Dans ce projet, les écoulements dans l’aspirateur d’une turbine de basse chute ont été simulés à l'aide de différents modèles de turbulence dont le modèle DDES, un hybride LES/RANS, qui permet de résoudre une partie du spectre turbulent. Déterminer des conditions aux limites pour ce modèle à l’entrée de l’aspirateur est un défi. Des profils d’entrée 1D axisymétriques et 2D instationnaires tenant compte des sillages et vortex induits par les aubes de la roue ont notamment été testés. Une fluctuation artificielle a également été imposée, afin d’imiter la turbulence qui existe juste après la roue. Les simulations ont été effectuées pour deux configurations d’aspirateur du projet BulbT. Pour la deuxième, plusieurs comparaisons avec des données expérimentales ont été faites pour deux conditions d'opération, à charge partielle et dans la zone de baisse rapide du rendement après le point de meilleur rendement. Cela a permis d’évaluer l'efficacité et les lacunes de la modélisation turbulente et des conditions limites à travers leurs effets sur les quantités globales et locales. Les résultats ont montrés que les structures tourbillonnaires et sillages sortant de la roue sont adéquatement résolus par les simulations DDES de l’aspirateur, en appliquant les profils instationnaires bidimensionnels et un schéma de faible dissipation pour le terme convectif. En outre, les effets de la turbulence artificielle à l'entrée de l’aspirateur ont été explorés à l'aide de l’estimation de l’intermittence du décollement, de corrélations en deux points, du spectre d'énergie et du concept de structures cohérentes lagrangiennes. Ces analyses ont montré que les détails de la dynamique de l'écoulement et de la séparation sont modifiés, ainsi que les patrons des lignes de transport à divers endroits de l’aspirateur. Cependant, les quantités globales comme le coefficient de récupération de l’aspirateur ne sont pas influencées par ces spécificités locales. / Draft tubes play a crucial role in elevating the available energy extraction of hydroturbines. In this project, turbulent flows in the draft tube of a low-head bulb turbine were simulated using, among others, an advance hybrid LES/RANS turbulent model, called DDES, which can resolve portions of the turbulent spectrum. Providing appropriate inflow boundary conditions for such models is a challenging issue. In this regard, different inflow boundary conditions were tested, including axisymmetric 1D profiles, and unsteady 2D inflow profiles that take runner blade wakes and vortices into account. Artificial fluctuation at the inlet section of the draft tube was also included to mimic the turbulence existing after the runner. Simulations were conducted for two draft tube configurations of the BulbT project. For one of them, intensive comparisons with experimental data were done for two operating conditions, one at part load and another in the sharp drop-off portion of the efficiency hill after the best efficiency point. This allowed to assess the effectiveness and shortcomings of the adopted turbulence modeling and boundary conditions through their effects on the global and local quantities. The results showed that the runner-related vortical structures and wakes are appropriately resolved using stand-alone DDES simulation of the draft tube flows. This is achieved by applying unsteady 2D inflow profiles along with adopting low dissipation scheme for the convective term. Furthermore, the effects of applying artificial turbulence at inlet were explored using separation intermittency, two-point correlation, energy spectrum and Lagrangian coherent structure concepts. These analyses revealed that the type of inflow boundary conditions modifies the details of the flow and separation dynamics as well as patterns of the transport barriers in different regions of the draft tube. However, the global quantities such as recovery coefficient are not influenced by these local features.

Le modèle de Skyrme est une théorie des champs classique non-linéaire qui possède des solitons topologiques comme solutions. De façon innovatrice, ces solutions sont utilisées par Skyrme pour décrire les baryons. Ceci est justifié puisque le modèle est une théorie efficace de la chromodynamique quantique à basse énergie et puisque les solutions possèdent un nombre topologique conservé qu'on associe au nombre baryonique. Le modèle original permet de décrire les propriétés des nucléons avec une précision de 30% ou mieux, mais a de la difficulté à reproduire les caractéristiques des baryons de masse élevé. En effet, l'énergie de liaison des multi-skyrmions est très grande et les densités baryonique possèdent des structures spatiales complexes et de type coquille ce qui s'éloignent des observations expérimentales. On propose alors une extension du modèle afin d'améliorer précisément les deux problèmes précédents. Ce travail introduit le modèle de Skyrme quasi-BPS dont les solutions analytiques saturent presque la borne de Bogomol'nyi améliorant les prédictions de l'énergie de liaison des multi-skyrmions. Le choix d'un potentiel permettant d'obtenir des densités baryonique de configurations pleines et l'ajout de l'énergie de Coulomb et de brisure de symétrie d'isospin à la masse du skyrmion sont les contributions originales de ce travail. Ceci permet d'obtenir des meilleurs prédictions physiques pour les multi-skyrmions de nombre baryonique très élevé que dans le modèle original et l'analyticité des solutions en simplifie le traitement. Pour ces raisons, le modèle de Skyrme quasi-BPS s'avère une bonne théorie efficace pour la descriptions des noyaux atomiques.