La thèse s'intéresse aux problèmes d'investissement en moyens de production électriques qui sont apparues ces dernières années dans les pays ayant déréglementé cette industrie. Dans la première partie, les différentes crises ayant affecté les marchés de gros sont recensées et analysées, afin de montrer qu'il existe un problème récurrent de sous-investissement en moyens de pointe. Les prix sur les marchés de gros sont tels que les coûts fixes de nouveaux moyens de pointe ne sont pas couverts, phénomène qui a trouvé l'appelation de « revenus manquants » dans la récente littérature théorique. Dans la deuxième partie, un modèle décrivant l'investissement d'un duopole en situation de demande aléatoire est comparé à un modèle de construction similaire décrivant l'investissement d'un monopole public. La comparaison permet de montrer qu'une hausse des prix de l'électricité sur les marchés de gros ne peut palier au problème des « revenus manquants ». Le résultat est étendu au cas d'un n-oligopole. La troisième partie s'intéresse donc aux différents moyens de mettre en place une rémunération complémentaire de la capacité des moyens de pointe. Différents designs de marchés de capacités sont analysés à la lumière des conclusions de la deuxième partie. Enfin, la nécessité de maintenir des capacités comme « réserves opérationnelles » est étudiée. La conclusion, à la lumière d'un article important de Joskow et Tirole, est que celles-ci doivent être considérées comme un bien public. Cet aspect semble empêcher la mise en place d'un marché de capacité, et, ce faisant pose un problème important à la déréglementation du secteur électrique.

L'objectif de cette étude visait à développer des modèles mathématiques capables de prédire la croissance d'une levure désacidifiante, la consommation des substrats carbonés (lactose, lactate) et la perte de poids au cours de l'affinage de fromages de type Munster. Afin de constituer une base de données servant à l'établissement et à la validation des modèles, des fabrications de fromages à pâte molle et à croûte lavée de type Munster ont été réalisées à échelle pilote et dans des conditions aseptiques. Cette base comporte les résultats des analyses microbiologiques et biochimiques ainsi que les conditions opératoires d'affinage. Les fromages ont été ensemencés, en plus de la flore lactique, par une flore d'affinage composée d'une levure (Debaryomyces hansenii 304, GMPA) et d'une bactérie de surface (Brevibacterium aurantiacum ATCC 9175). L'affinage de ces fromages est limité à 14 jours et il est conduit selon un plan d'expérience à deux facteurs (température et humidité relative) et à trois niveaux (respectivement 8, 12, 16° C et 85, 93, 99 %). Trois modèles mathématiques ont été construits selon une approche mécanistique dynamique. Le premier modèle (dit "microbiologique") a réussi à prédire, en fonction de la température et de l'humidité relative du hâloir, la croissance de D. hansenii et l'évolution des concentrations en lactose et en lactate durant l'affinage. Les erreurs résiduelles standard (RSE) de ces prédictions sont satisfaisantes par comparaison avec l'écart type des trois essais répétés sous les conditions du point central du plan d'expérience. Les résultats de ce modèle permettent de soutenir l'hypothèse selon laquelle D. hansenii consomme le lactose pour sa croissance et le lactate pour sa maintenance. Le deuxième modèle (dit "perte de poids") est parvenu à prédire la perte de poids et l'évolution de la teneur en matière sèche, avec des erreurs résiduelles standard moyennes proches de la précision des mesures expérimentales. Les résultats obtenus avec ce modèle montrent que les hypothèses avancées, selon lesquelles le lactose et le lactate sont complètement oxydés en CO2 et en O2 et que l'oxygène atmosphérique ne contribue pas à la perte de poids, sont justifiées. Le troisième modèle (dit "généralisé") a été construit en combinant les deux modèles précédents. Il est capable de prédire, dès le début de l'affinage, en fonction de la température et de l'humidité relative de la chambre d'affinage, la croissance de D. hansenii, la perte de poids et les évolutions des concentrations en lactose et lactate ainsi que la teneur en matière sèche du fromage.

Doctorat en sciences sociales, politiques et économiques / info:eu-repo/semantics/nonPublished

Un système multiagents (SMA) est constitué de plusieurs agents qui interagissent dans un environnement pour atteindre leurs objectifs. Une interaction peut être une coopération si les agents coopèrent. L'interaction peut aussi être une négociation si les agents sont en compétition. Il est toutefois possible que les agents soient situés dans un environnement semi compétitif où ils peuvent, à la fois, coopérer et être en compétition. Dans ce mémoire, nous présentons un modèle dynamique de pilotage d'un SMA semicompétitif qui permet à l'agent de choisir ses partenaires d'interaction et d'élaborer son comportement vis-à-vis des agents. Le modèle ainsi développé est basé sur l'évaluation de la réputation des agents. Si l'agent est bien réputé, alors il constituera un bon partenaire d'interaction. L'évaluation de la réputation est effectuée à partir de plusieurs sources d'informations comme l'historique des expériences directes avec l'agent à évaluer, et le témoignage des autres agents de la société classés selon des groupes. Le modèle dynamique de pilotage a été validé en considérant les trois scénarios les plus représentatifs du contexte étudié. Le premier scénario suppose qu'il n'existe aucune collaboration entre les agents (1). Le deuxième scénario implique les agents dans une relation de compétition (2). Enfin, le troisième scénario utilise le modèle dynamique de pilotage pour permettre à l'agent de choisir son comportement vis-à-vis ses partenaires d'interaction (3). Les résultats de la simulation ont montré que les agents ayant des bonnes réputations ont gagné en termes d'utilité. / A multi-agent system (MAS) consists of several agents interacting in an environment to achieve their goals. An interaction can be cooperation if the agents cooperate. It can be also a negotiation if the agents are in competition. However, it is possible that the agents are located in a semi-competitive environment. They can both be in competition and cooperate. In this memory, we present a dynamic model of piloting a semicompetitive MAS which allows the agent to choose its interaction partners and to define its behaviour towards the agents. This developed model is based on the evaluation of the reputation of the agents. If the agent has a good reputation, then it will be considered as a good interaction partner. The evaluation of the reputation is carried out starting from several sources of information like the direct experience with the evaluated agent and the witness of the other agents. The dynamic model was validated by a simulation of three scenarios. The first scenario supposes that there is not any collaboration between the agents. The second scenario implies the agents in competition. Lastly, the third scenario uses the dynamic model of piloting to allow agent to choose its behaviour toward its interaction partners. The simulation results showed that the agents having good reputations improved their utility.

Ce mémoire décrit le modèle numérique Smoothed Particle Ice Dynamics Equations for Rivers (SPIDER) développé par Nolin, Roubtsova et Morse afin de modéliser le transport et l'accumulation de glace à la surface d'une rivière. Dans SPIDER, la méthode Smoothed Particle Hydrodynamics (SPH) est utilisée pour résoudre les équations de dynamique de la glace en deux dimensions. Le modèle SPIDER a été utilisé afin de simuler la formation d'embâcle fluvial. Les profils d'embâcle obtenus sans considérer le frottement des berges sont identiques à ceux présentés dans les études antérieures. Les profils obtenus en considérant le frottement des berges diffèrent de ceux présentés dans les études antérieures mais sont conformes à la méthode utilisée pour modéliser le frottement. Deux études paramétriques ont été réalisées afin de caractériser l'influence de l'angle de frottement interne de la glace ([phi]) et du coefficient de traînée eau-glace (cw) sur les profils d'embâcle. Ces études montrent que plus l'angle ʹ est grand, plus l'embâcle est long et moins la glace au pied de l'embâcle est épaisse. De même, plus cw est grand, plus l'embâcle est court et plus l'épaisseur de glace au pied de l'embâcle est grande. Le relâchement d'un embâcle fluvial a été simulé en couplant SPIDER au modèle hydrodynamique KIMBULA développé par Roubtsova V. et Kahawita R. à l'École Polytechnique de Montréal. Les résultats préliminaires simulés par ce modèle hybride sont similaires à ceux présentés dans les études antérieures.

Les travaux de cette thèse portent sur la conception optimale des systèmes de production constitués d'équipements ayant une fiabilité donnée. Les systèmes étudiés utilisent la redondance et/ou des stocks tampons comme technique d'amélioration de la performance. La méthodologie proposée développe des modèles de conception optimale, des méthodes d'évaluation et des algorithmes de résolution robustes et efficaces à base de métaheuristiques pour ces systèmes. Cette méthodologie est validée en l'appliquant à trois problèmes de conception optimale des systèmes de production caractérisés comme N-P difficiles : 1. Le problème d'allocation de la redondance (PAR) : Un algorithme hybride très efficace est proposé pour résoudre ce problème d'optimisation. Cet algorithme est basé sur une hybridation de l'algorithme à colonies de fourmis avec l'algorithme du grand déluge étendu. Les résultats numériques obtenus pour 33 problèmes tests, existant dans la littérature du PAR, montrent que l'algorithme proposé est très compétitif par rapport aux autres approches présentées dans la littérature. 2. Le problème d'allocation optimale des stocks : Deux types de lignes sont traitées, à savoir les lignes homogènes et les lignes non-homogènes. Une nouvelle approche basée sur l'algorithme du grand déluge étendu est proposée. Une étude comparative avec les meilleurs résultats publiés dans la littérature (recuit simulé) montre l'efficacité de cette approche tant au niveau du temps de convergence qu'au niveau de la qualité des solutions. 3. Le problème d'allocation optimale des stocks et de la redondance : Une méthode d'évaluation du taux de production d'une ligne de production série-parallèle est proposée et validée par simulation. Deux algorithmes utilisant des métaheuristiques, à savoir l'algorithme à colonies de fourmis et le recuit simulé, sont développés et comparés. Les résultats numériques obtenus, pour différentes instances générées aléatoirement, illustrent l'efficacité de l'algorithme à colonies de fourmis lorsqu'il est couplé à une procédure d'amélioration.

Le défi relevé dans ce mémoire est d'élaborer un nouveau modèle d'écoulement « Q3D » permettant le calcul des débits dans les cours d'eau naturels, en présence et absence d'un couvert de glace, à partir des relevés de vitesses et de profondeurs enregistrés par un instrument installé dans le fond du cours d'eau. Ce modèle, dans sa composante d'éléments finis, résout l'équation de Poisson pour générer la forme de la distribution des isotâches (contours d'égale vitesse). Les vitesses sont ensuite obtenues à l'aide de l'équation proposée par l'auteur du modèle POME (Principle Of Maximum Entropy) reliant ces contours au degré de turbulence caractérisant l'écoulement. Q3D est testé sur les données provenant des instruments SONTEK Argonaut Shallow water (SW) installés à deux sites voisins sur la rivière Châteauguay, QC. Les premiers résultats obtenus indiquent que Q3D reste stable dans l'estimation des paramètres hydrauliques, malgré le changement des conditions d'écoulement. La comparaison des débits du modèle Q3D avec les jaugeages réalisés indique que Q3D est sensible à la forme géométrique de la section d'écoulement. En ayant l'information d'un seul jaugeage, il est possible d'introduire un facteur de correction des débits. Une fois appliqué, Q3D présente des résultats excellents : Pour le site aval de Châteauguay, été comme hiver, 12 fois sur 13 l'erreur sur le débit est inférieure à 7%. Cela est très prometteur pour de futures investigations. / This study introduces a new model "Q3D" to calculate flow in the natural rivers during summer and winter conditions using velocity and depth data obtained by an instrument installed in the bottom of the river. The model uses a finite element mesh to solve the Poisson's equation to generate the form of the distribution of the isovels (contours equal velocity) in the gauged cross-section. As a function of turbulence intensity, numerical values are then assigned to the isovels using the equation suggested by the authors of the POME model (Principle Of Maximum Entropy). Q3D is tested on the data of the Châteauguay River, QC. Results indicate that Q3D remains stable in the estimate of the hydraulic parameters for a wide variety of flow conditions. The comparison of Q3D with measured discharges indicates that Q3D is sensitive to the geometrical form of the cross-section of stream. However, when one single gauging is used to calibrate Q3D, it yields excellent results: the error in discharge was less than 7% in 12 times out of 13 measurements. Therefore, winter or summer, Q3D is a very promising candidate for use in conjunction with the velocity index methods to monitor streamflow.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2015-2016 / Ce mémoire investigue le phénomène de reconnexion visqueuse de deux tourbillons initialement placés de façon orthogonale ou antiparallèle. Les équations de Navier-Stokes pour un fluide incompressible sont résolues directement (DNS) à l’aide d’un algorithme pseudospectral utilisant des expansions périodiques dans les trois directions cartésiennes. La condition de circulation nulle inhérente à cette méthode numérique est contournée en résolvant les équations dans un repère tournant approprié. Une méthode simple utilisant des lignes de vorticité est proposée afin de calculer le pourcentage de reconnexion instantané η de deux tourbillons. Cette méthode est également employée pour séparer le champ de vorticité en ses composantes reconnectées et non-reconnectées, ce qui facilite l’identification visuelle des différentes étapes du processus de reconnexion. Finalement, l’échelle de temps de la reconnexion tourbillonnaire Trec est calculée pour différents nombres de Reynolds (500 ≤ Re ≤ 10000). Il est trouvé que l’ordre de grandeur de Trec varie de façon continue de [symbol] à mesure que Re augmente. / This work focuses on the viscous reconnection phenomenon of two vortex tubes that are initially antiparallel or orthogonal to each other. The incompressible Navier-stokes equations are solved directly (DNS) using a Fourier pseudospectral algorithm with triply periodic boundary conditions. The associated zero-circulation constraint is circumvented by solving the governing equations in a proper rotating frame of reference. A simple method using vortex lines is proposed to compute the instantaneous reconnection level η of two vortices. The proposed method is also used to split the vorticity field into its reconnected and non-reconnected parts, which allows for a clear and intuitive visual identification of the different reconnection phases. Finally, the Reynolds number dependence of the reconnection timescale Trec is investigated for 500 ≤ Re ≤ 10000. The scaling is found to vary continuously as Re is increased from [symbol].

La géosimulation est une approche qui permet la modélisation et la simulation de phénomènes dynamiques. Elle est appliquée dans un grand nombre de domaines tels que l’urbanisme, la gestion de l’environnement et la santé publique. La géosimulation se base sur des données géo-référencées qui fournissent le contexte spatial virtuel de la simulation. La dynamique du phénomène simulé est généralement fondée sur des modèles, des algorithmes et/ou des équations mathématiques. Dans le domaine de la santé publique qui nous concerne ici, plusieurs approches de simulation ont été proposées pour simuler la propagation des maladies vectorielles ou zoonoses. Une zoonose est une maladie propagée par un virus ou une bactérie transmise par des insectes (ex. moustiques, tiques) à des animaux (rongeurs, oiseaux, mammifères) qui peuvent être infectés, et à leur tour transmettre le virus ou la bactérie à des insectes sains qui les piquent. Éventuellement les insectes peuvent aussi piquer les humains et leur transmettre la maladie. Dans ces approches, on représente habituellement la dynamique de la zoonose par des modèles mathématiques (dits ‘à base de compartiments’) en faisant l’hypothèse que l’espace est homogène. Ainsi elles ne tiennent pas compte de l’influence du ‘paysage’ (en particulier de la couverture du sol). Pourtant prendre en compte les caractéristiques spatiales est important si on veut simuler de façon plausible la dynamique des zoonoses en tenant compte des zones favorables à la survie ou à la prolifération de certaines espèces (ex. moustiques, tiques) de leur dispersion par les espèces transportrices (ex. oiseaux, chevreuils) et de leurs comportements de mobilité. Dans le projet ZoonosisMAGS, dans lequel s’insère ce travail de maîtrise, on vise à développer une approche générique de simulation de la propagation des zoonoses en utilisant des données géo-référencées et en créant un environnement géographique virtuel efficace pour tenir compte de l’influence de la couverture du sol. Pour la dynamique du phénomène qui reflète l’évolution des espèces (des stades œuf, larve, nymphe à adulte), leurs interactions et leur statut épidémiologique (sain, susceptible, infecté) on utilise un modèle enrichi à base de compartiments qui représente les stades d’évolution biologiques et épidémiologiques des différentes espèces impliquées dans la zoonose. Les populations de chaque espèce en interaction et à divers stades de maturité sont associées aux cellules qui représentent les unités de l’environnement géographique virtuel. La transition des individus d’une espèce d’un compartiment à un autre, ainsi que l’interaction entre les espèces se fait à travers des calculs d’expressions mathématiques à chaque pas de simulation, et ceci pour chaque cellule. Par exemple, dans le cadre d’une simulation de la propagation de la maladie de Lyme sur 365 jours d’une année, sur une région d’environ 300km2, on doit évaluer environ 150 expressions mathématiques pour une trentaine de compartiments et environ 20000 cellules à chaque pas de temps. On comprend qu’un défi majeur est celui de l’efficacité de la spécification et de l’évaluation des expressions mathématiques. Aussi, nous traitons dans ce mémoire une première problématique qui concerne la spécification et l’évaluation efficace des expressions mathématiques et booléennes pour de telles simulations. C’est dans ce contexte que nous avons développé un système qui permet à l’usager de spécifier les expressions mathématiques requises par de telles simulations, ainsi qu’un système qui les évalue par la suite de manière très efficace. La spécification des modèles de compartiments des espèces, de leurs interactions et de l’infection peut être une tâche complexe pour certaines zoonoses comme la maladie de Lyme. À cet effet, nous avons traité une deuxième problématique visant à simplifier la tâche de spécification des modèles en créant un module graphique couplé au système d’expressions mathématiques. En nous basant sur la librairie de dessin Visio, nous avons développé un tel système qui est complètement intégré au simulateur ZoonosisMAGS (version C++). Dans ce mémoire nous présentons le contexte de la recherche, les étapes d’analyse et de conception des systèmes développés. Nous montrons à travers des cas d’utilisation comment ces systèmes sont pratiquement utilisés grâce aux interfaces utilisateurs pratiques que nous avons développées. Nous rendons compte de tests d’efficacité qui ont été conduits pour montrer la contribution de notre travail à travers les résultats d’exécution des expressions mathématiques dans une simulation utilisant un modèle à base de compartiments de la maladie de Lyme.

Les matériaux composites constituent de nos jours, de par leurs propriétés très prisées, des matériaux incontournables dans plusieurs secteurs industriels telles les industries aéronautique, automobile et navale. Bien que les applications industrielles aient été axées durant des décennies sur les matériaux composites à matrices thermodurcissables, les composites thermoplastiques à renforts continus connaissent actuellement un grand essor. Dans ce sens, des recherches considérables ont été menées sur la modélisation du comportement des matériaux thermoplastiques chargés de fibres continues à température ambiante. Cependant, la modélisation du comportement de ces matériaux à hautes températures et dans les conditions de mise en forme demeure une problématique d'actualité et relativement peu étudiée. Ce travail de thèse se veut une contribution à la modélisation du comportement des composites thermoplastiques à renforts continus dans les procédés de mise en forme. Deux directions de recherches ont été explorées. La première direction consiste à développer une loi constitutive viscoélastique orthotrope pour les milieux incompressibles, et sous l'hypothèse d'un processus isotherme. Dans cette optique, on se propose d'exploiter la théorie du modèle K-BKZ, destiné à la rhéologie des fluides viscoélastiques isotropes, pour une extension aux cas orthotrope et transversalement isotrope. Les tenseurs d'orientation des fibres sont introduits comme arguments de la fonctionnelle densité d'énergie libre pour tenir compte de la symétrie matérielle. La réponse du milieu est formulée en considérant les apports des constituants du milieu, pris séparément, augmentés de la réponse due aux interactions entre la résine et les fibres. Ce faisant, la densité d'énergie libre sera décomposée en un produit d'une fonction mémoire, responsable de la relaxation, et d'une fonction des déformations, laquelle a été scindée en la somme de quatre composantes : i) une fonction rendant compte du comportement de la matrice, ii) une composante responsable de la réponse de la première famille de fibres, iii) une composante responsable de la seconde famille de fibres et iv) une fonction rendant compte des interactions matrice-fibres. La seconde direction vise l'étude du phénomène de l'endommagement susceptible de toucher les composites thermoplastiques à renforts continus dans des conditions proches de celles de la mise en forme. Dans cette perspective, une loi constitutive hyperélastique orthotrope avec endommagement sera proposée. On fait appel à une approche thermodynamique classique, basée sur l'existence d'un potentiel d'énergie libre duquel découle le tenseur de contraintes, et au principe de la contrainte effective, bien connu de la théorie de l'endommagement. À l'instar du modèle viscoélastique orthotrope, la fonctionnelle densité d'énergie libre sera décomposée en la somme de quatre composantes responsables des réponses de la résine, des deux familles de fibres et des interactions matrice-fibres. Pour introduire le phénomène de l'endommagement, on associe à chaque constituant du milieu (i.e. les fibres et la résine) une variable d'endommagement isotrope. Chacune de ces variables est associée à la composante de l'énergie libre rendant compte du comportement du constituant qu'elle caractérise. Une loi d'évolution de l'endommagement isotrope est également associée à la composante de l'énergie libre responsable des interactions matrice-fibres. Les hypothèses d'un milieu incompressible et d'un processus isotherme ont été notamment retenues. Une campagne d'essais expérimentaux, axée sur des essais de traction dans les directions 0°/90° et ±45°, a été réalisée pour caractériser l'endommagement d'une classe de composite thermoplastique à renforts continus. Les résultats de ces essais ainsi que ceux rencontrés dans la littérature ont été exploités afin de réaliser une procédure d'identification inverse visant à tester l'aptitude des modèles de comportement proposés à reproduire des résultats expérimentaux. Des simulations par éléments finis de la mise en forme d'un composite thermoplastique à renforts continus ont été réalisées. Les résultats obtenus montrent que les modèles utilisés reproduisent des résultats physiquement admissibles et en concordance avec des observations expérimentales.