Le self-induced ion plating (S.I.I.P.) est un procédé hybride entre l’évaporation sous vide et la pulvérisation cathodique magnétron. L’ionisation d’un gaz rare (argon) est obtenue grâce à la différence de potentiel appliquée entre une cathode d’étain (matériau à évaporer) et un substrat (tôle d’acier à revêtir) relié à la masse. Les ions positifs Ar+ sont accélérés vers la cible par le champ électrique et bombardent celle-ci avec une énergie importante. Ce bombardement ionique de la surface de la cible et l’isolation thermique du système induisent l’échauffement de la cible, sa fusion et enfin son évaporation. Afin de garantir un bon rendement au procédé, un magnétron d’aimants permanents est placé sous la cible à évaporer (phénomène de confinement magnétique). Le gaz métallique évaporé se solidifie au contact du substrat formant la couche de dépôt souhaitée. L’objectif de cette thèse revêt deux aspects : la mise en évidence des phénomènes physiques présents au sein du système et la simulation numérique de ceux-ci afin de déterminer l’épaisseur de dépôt obtenue sur la largeur du substrat. Le modèle de simulation du S.I.I.P. est scindé en quatre étapes. La première étape consiste à calculer le champ magnétique produit par le magnétron. Ce champ magnétique dicte le comportement des électrons qui sont à l’origine de l’ionisation de l’argon et de ce fait, il est une donnée indispensable à la deuxième étape à savoir, la modélisation du bombardement ionique et la détermination de la distribution du flux de chaleur transmis à la cible par celui-ci. Dans cette deuxième étape nous utilisons un modèle de suivi de particules basé sur une méthode statistique connue sous le nom de méthode de Monte-Carlo. Le phénomène du bombardement ionique modélisé, au cours de la troisième étape, nous calculons les échanges thermiques qui ont lieu dans le S.I.I.P. Les modes de transfert considérés sont : la conduction, le rayonnement et la convection (convection naturelle, convection électromagnétique et convection Marangoni). Enfin, le champ de température obtenu par la modélisation thermique, couplé à la théorie de l'évaporation, nous permet de déterminer le flux de matière évaporé et déposé sur le substrat suite à son passage au-dessus du système de dépôt. L’étude du S.I.I.P. a mis en évidence la complexité et la multidisciplinarité des phénomènes physiques mis en jeu dans les procédés PVD. De plus les résultats obtenus, validés à partir des mesures, nous ont permis de mettre en évidence l’importance significative du phénomène de convection Marangoni.

Une scierie présente souvent une grande variabilité du niveau de production journalier. Une partie de cette fluctuation est causée par la variation de la dimension de la matière première. Cette variation engendre des problèmes importants de flux du bois. Pensons notamment aux surcharges de pièces à certains centres-machines, alors que d'autres manquent de matériel. Cela occasionne une sous-utilisation du potentiel des équipements de sciage et résulte en une perte de productivité. Certaines scieries ont implanté des méthodes de triage en forêt selon la dimension des tiges. Les tiges sont séparées en deux ou trois catégories de diamètres. Elles sont empilées séparément dans la cour à bois de la scierie. La scierie est ensuite alimentée en utilisant les empilements en alternance. Cela permet d'exercer un certain contrôle sur la façon dont le bois entre dans la scierie. Toutefois, les gains éventuels de cette façon d'opérer n'ont jamais été documentés ou évalués clairement. Comme le tri des tiges engendre des coûts supplémentaires, il apparaît important d'établir s'il y a un gain de productivité pour la scierie, afin d'évaluer la rentabilité globale de cette opération de triage. L'objectif principal de ce projet est donc d'évaluer les éventuels bénéfices associés à l'utilisation d'une matière première triée selon ses dimensions. Le logiciel QUEST, un simulateur de procédés manufacturiers par événements discrets, a été configuré pour représenter les aspects concernant le flux des pièces d'une scierie. Auparavant, le logiciel Optitek, un simulateur de débitage, a été mis à contribution pour obtenir des solutions de débitage, propres à la scierie étudiée, pour 255 billes échantillons scannées. Le design de la scierie a été recréé dans le logiciel QUEST. Les vitesses d'avance des convoyeurs et les vitesses d'alimentation des équipements ont été incluses ainsi que les différentes options de cheminement des pièces dans la scierie. L'efficacité des équipements a été considérée en incorporant la fréquence et la durée des arrêts par machine. De plus, des efforts ont été mis afin de bien cerner et de reproduire, dans le modèle, les décisions humaines ayant une influence sur le flux des pièces. Un modèle de flux représentant la situation actuelle de la scierie a ainsi été obtenu. Le modèle a été ensuite modifié afin d'y inclure une étape de tri des tiges selon trois classes de diamètre à hauteur de souche (DHS). Pour gérer la façon dont le bois des différentes classes de dimension est amenée à la scierie, différents indicateurs de l'état du système de production ont été utilisés à tour de rôle. Ces indicateurs sont : le volume moyen des billes écorcées des cinq dernières minutes, le niveau d'occupation des bacs de classement des billes et le volume moyen des tiges tronçonnées. Pour chacun de ces indicateurs, des bornes ont été établies, par essais et erreurs, en vue d'appliquer des recettes d'alimentation à partir des trois catégories de dimension des tiges. Le modèle de référence, c'est-à-dire celui où il n'y a aucun tri de la matière première, a été utilisé afin de simuler 65 factions de travail. Les résultats ont ensuite fait l'objet d'analyses statistiques afin de vérifier l'existence d'une relation entre le volume de sciage produit dans une faction de travail et certaines variables explicatives. Des équations, obtenues par analyses de régression multiple, ont été établies afin d'exprimer la productivité de la scierie sous l'angle de la production de sciage brut en volume nominal. Les variables explicatives qui ont été retenues dans ces équations sont le volume des billes sciées et les arrêts au poste d'éboutage. De la même façon, le modèle incluant le tri des tiges en trois classes de dimension a été simulé. Des séries de simulations ont été réalisées pour chacune des trois méthodes de contrôle de l'alimentation. Les résultats ont été analysés par régression multiple et des équations exprimant la production en fonction du volume des billes sciées et des arrêts au poste d'éboutage ont été établies. En comparant les différentes équations de productivité, et ce au même niveau des variables explicatives, il a été possible de vérifier les gains possibles associés au triage des tiges. Il ressort des résultats que le tri des tiges selon leurs dimensions peut améliorer de façon appréciable la productivité de la scierie. Cependant, la méthode utilisée pour contrôler l'alimentation influence de façon très importante les gains.

Titre de l'écran-titre (visionné le 30 novembre 2023) / Dans cette thèse, la dynamique d'excitation électronique et d'ionisation de H₂ dans une impulsion laser intense est étudiée en utilisant une approche ab initio pour résoudre numériquement l'équation de Schrödinger dépendante du temps (TDSE) pour ce système. Nous avons développé une nouvelle méthodologie utilisant des fonctions B-Splines comme base de développement des fonctions d'onde multi-configurationnelles du système. Afin de décrire et d'analyser la dynamique d'ionisation et d'excitation électronique, nous faisons appel au programme MEDYS (Many-Electron-Dynamics System), conçu à notre laboratoire en interne, et dont l'adaptation en base B-Spline pour donner la version MEDYS-BSpline est un des objectifs de la thèse. Ce programme utilise une méthode d'interaction de configuration dépendante du temps (TDCI) pour décrire la dynamique temporelle de l'ionisation sur les voies de l'espace lié et celui des cations. En application de la méthodologie, le travail continue avec la détermination du régime d'ionisation, tunnel ou multiphotonique, quand la molécule H₂ est soumise à un rayonnement intense dans l'infra-rouge proche (de longueur d'onde λ = 800 nm). Le travail entreprend également une évaluation numérique de l'approximation du Champ Fort en comparant les résultats de calculs de la dynamique électronique utilisant l'approximation de Strong Field Approximation (SFA) avec ceux utilisant une représentation complète et non-SFA du propagateur de l'électron ionisé. / In this thesis, the dynamics of electronic excitation and ionization of H₂ in an intense laser pulse are studied using an ab initio approach to numerically solve the time-dependent Schrödinger equation (TDSE) for this system. We have developed a new methodology using B-Spline functions as a basis for developing the multiconfigurational wave functions of the system. To describe and analyze the dynamics of ionization and electronic excitation, we employ the in-house program called Many-Electron-Dynamics System (MEDYS), and one of the objectives of the thesis is to adapt it to the B-Spline basis, resulting in the MEDYS-BSpline version. This program utilizes a time-dependent configuration interaction (TDCI) method to describe the temporal dynamics of ionization in both bound and cationic states. Applying this methodology, the work continues with the determination of the ionization regime, either tunneling or multiphoton, when H₂ is subjected to intense radiation in the near-infrared (wavelength λ = 800 nm). The work also undertakes a numerical evaluation of the Strong Field Approximation (SFA) by comparing the results of electronic dynamics calculations using the SFA approximation with those using a full and non-SFA representation of the ionized electron propagator.

L’objectif de l’étude est d’évaluer des stratégies en flux tiré et flux poussé pour la planification de la production en utilisant comme cas d’étude une usine de sciage située au Québec. Un outil de planification de la production (APS, Advanced Planning and Scheduling System), basé sur une architecture logicielle distribuée, simule les divers processus de production et planification des opérations de la scierie (approvisionnement, sciage, séchage, rabotage, entreposage et distribution) en les représentant comme des agents logiciels spécifiques. Pour les expériences, des configurations de pilotage définies par deux facteurs contrôlables sont utilisées, à savoir : la position du point de découplage et les niveaux des contrats pour une famille de produits. Ensuite, des scénarios ont été définis par deux facteurs non contrôlables : la qualité de l’approvisionnement et le différentiel des prix pour les produits visés par les contrats. Ces configurations et ces scénarios génèrent un plan d’expériences de niveaux mixtes avec cinquante-quatre expériences. Pour chacun des 54 plans de production généré par l’APS trois indicateurs de performance sont calculés : le taux de service, le niveau de stock de produits en-cours et le profit potentiel de la production (Potential Monetary Throughput). Les résultats montrent une relation directe entre le niveau de pénétration du point de découplage, le flux tiré et le taux de service, pour chaque niveau de demande des produits tirés par les contrats. Ainsi, une planification qui incorpore le flux tiré génère en moyenne une augmentation de performance de 100% dans le niveau de service, par rapport à une planification en flux poussé pour tous les niveaux de contrats. Cette augmentation de la performance du niveau de service a un coût financier d’approximativement 7% du profit potentiel, lequel devrait être compensé par de meilleures conditions des contrats (un premium payé sur le prix marché) et des coûts de stockage plus bas. Ce « compromis » s’avère un résultat direct de la divergence dans la production de bois d’œuvre. Ainsi, dans un contexte d’affaires qui privilégie la qualité du service et où les clients sont prêts à payer pour ce service, l’utilisation de ces stratégies de planification de la production, où le point de découplage est situé en amont de la chaîne, est une opportunité possible d’avantage concurrentiel. / The objective of this study is the evaluation of pull and push strategies in lumber production planning using a Quebec sawmill as case study. An Advanced Planning and Scheduling System (APS), based on a distributed software architecture, simulates the main operations planning and production processes of the sawmill (sourcing, sawing, drying, finishing, warehousing and delivery) representing them as autonomous software agents. Push and pull strategies are simulated using different penetration positions of the demand information decoupling point over the value chain. To set experiments, configurations are defined by two controllable factors, namely: the decoupling point position and the level of contracts for a product family. Following, a set of scenarios are generated by two uncontrollable factors: the quality of supply and market prices differential for products under contracts. These configurations and scenarios leads to a mixed levels experimental design with fifty four runs. Three performance indicators: orders fill rate, work in process, and potential monetary throughput; are calculated for every one of the 54 production plans generated by the APS. Results show a direct relation between the orders fill rate and the position of the decoupling point, pull strategy, for the three levels of demand on products under contract. Accordingly, at every demand level, production plans under pull strategies generate improvements of 100% compared with equivalent plans under push strategy. This service level performance improvement has a financial cost of about 7% of the Potential Monetary Throughput which should be compensated externally with better contract conditions and internally by lower costs of inventory management. This trade-off seems to be a direct consequence of the divergent nature of lumber production. Consequently, in a business context that privileges service quality and where customers are willing to pay for it, the use of this kind of demand driven strategies in production planning represents a source of competitive advantage.

Le train d’atterrissage est une structure assurant le mouvement de l’avion en vol-sol et sol-vol pour laquelle beaucoup de travaux restent encore à compléter afin d’avoir une structure performante, durable et qui demande moins de maintenance. Ce travail de recherche est effectué dans le but de contribuer à la modélisation avancée et au développement d’un outil de simulation dynamique d’un train d’atterrissage à amortissement oléopneumatique. Le développement d’un modèle mathématique est réalisé à cet égard, et ce dernier essaie de reproduire le fonctionnement physique d’un système train d’atterrissage en particulier celui de l’amortisseur. L’amortisseur d’un train d’atterrissage est une structure complexe composée de plusieurs sous-ensembles ayant des rôles complémentaires et qui aide l’avion à amortir les chocs pendant l’atterrissage, le roulage et le décollage de l’avion. La structure d’un amortisseur oléopneumatique est constituée des plusieurs chambres contenant du gaz ou un fluide hydraulique. Cependant la structure interne de l’amortisseur demande beaucoup de connaissances physique et technique afin de bien comprendre sa conception interne et externe. Les pneumatiques sont aussi pris en compte dans le modèle et elles sont caractérisées par deux modèles à savoir un basé sur ADAMS/Aircraft et un autre modèle développé en éléments finis sur ABAQUS. Ces aspects sont réunis dans un code numérique développé dans VBA dans un but de caractériser, modéliser et simuler un train d’atterrissage. / Although an aircraft landing gear has a well-known role of allowing the aircraft to complete its movement during landing, taxiing and takeoff, some works are still required today in order to obtain an efficient, more reliable and environmentally friendly landing gear structure. The scope of this project is to develop a numerical simulator of an aircraft landing taking into account the dynamic of the landing gear and the internal design of the shock strut. Landing gear design can be affected by phenomena such as sudden temperature variations as well as extreme temperatures. This can change the behavior of the hydraulic fluid of the shock strut. This work will be conducted by developing a mathematical model computer code in VBA that describes the landing gear physical function. The role of a landing gear shock strut is to absorb the kinetic energy during the landing, taxiing and takeoff phases of the aircraft. The landing gear shock strut is an assembly of different chambers that have a well-defined role. These chambers contain a gas and a hydraulic fluid and work together to provide the force required to maintain the system. These aspects are combined in a numerical code developed with VBA. In this project the effect of the tires was taken into account and characterized using two different models (finite elements and analytical model) using two different numerical tools (Abaqus and Adams). The scope is to compare and see how the two models are able to describe the energy absorption process.

Cet ouvrage a trait à la modélisation, à la conception de simulateurs et au contrôle de procédés. La simulation de procédés est employée largement pour l'enseignement des notions relatives au contrôle et aussi dans l'industrie pour, par exemple, tester au préalable de nouvelles stratégies de commande. Afin d'évaluer la pertinence de proposer une approche de conception de simulateurs, cet ouvrage passe d'abord en revue les principaux logiciels commerciaux de simulation disponibles sur le marché dans les domaines du génie chimique, du génie métallurgique et du traitement des eaux. Ensuite, différentes approches de conception de simulateurs de procédés dont la description du fonctionnement est basée sur des équations différentielles et algébriques sont présentées. À cet effet, quelques méthodes d'implantation d'un système d'équations de cette nature à l'aide de l'environnement de développement causal Matlab/Simulink sont exposées. Puis, étant donné l'importance croissante de la nouvelle approche de modélisation acausale utilisant l'approche orienté-objet et orienté-équation, deux environnements (Maple/MapleSim et Dymola) comportant une passerelle vers Matlab/Simulink et supportant cette approche via le langage Modelica sont présentés. Les simulateurs réalisés au cours de ce projet sont un réservoir avec élément chauffant, un ensemble de quatre procédés simples caractérisés par des réactions chimiques à l'intérieur de réacteurs, un bioprocédé de culture de microalgues en photobioréacteur et, finalement, un procédé minéralurgique de broyage couplé à la séparation par flottation. Pour les deux derniers procédés, des exemples de stratégies de contrôle basées sur la commande algébrique en continu sont présentés afin d'illustrer l'utilisation des simulateurs dans un contexte pratique.

La physique quantique est un domaine qui est encore aujourd'hui difficile à approcher numériquement ; les calculs s'y rattachant souffrent du grand nombre de degrés de liberté que comportent les systèmes quantiques intéressants. Nous explorons deux alternatives aux solutions traditionnelles en les appliquant à des problèmes de la physique quantique. Nous étudions d'abord l'effet tunnel dans un potentiel double puits à l'aide du concept de l'action quantique. Nous montrons que l'action quantique reproduit avec précision la physique du système dans le régime quantique et qu'elle présente de nombreuses caractéristiques intéressantes. Nos résultats suggèrent aussi que l'action quantique peut être utilisée pour des potentiels asymétriques. Nous présentons ensuite la construction d'un Hamiltonien Monte Carlo pour la théorie de jauge sur réseau U(l). Nous calculons explicitement la mesure exacte de l'intégrale' de chemin pour U(l). Enfin, nous présentons comment ces résultats peuvent être généralisés pour des théories de jauge non-abéliennes SU(N).

Les systèmes de communication sans fil de nouvelle génération nécessitent des techniques avancées afin de supporter un nombre élevé d'usagers et d'augmenter le débit sans avoir besoin de largeur de bande additionnelle. Une de ces techniques est l'utilisation des systèmes MIMO (Multi-Input, Multi-Output). De tels systèmes augmentent le débit, la portée et la robustesse aux évanouissements du canal sans consommer de largeur de bande additionnelle. Dans cette thèse, nous modélisons la norme 802.11 pour les réseaux ad-hoc à accès MIMO. Nous dégageons des scénarios de simulation permettant de mettre en lumière certaines situations problématiques des réseaux sans-fil. Ce mémoire introduit dans un premier temps le canal sans fil et les spécifications de la norme 802.11. Dans un deuxième temps, nous introduisons le simulateur réseau NS-2 ainsi que certaines simulations des réseaux ad-hoc utilisant la norme 802.11. Et pour finir on a testé notre extension MIMO par des scénarios de simulation de réseaux étendus.

Nous avons effectué 14 simulations de collisions majeures de galaxies spirales riches en gaz. Ces simulations ont été réalisées grâce à GCD+, un algorithme qui inclue la gravité, l'hydrodynamique, la formation stellaire et un traitement détaillé de l'enrichissement en métaux. Nous avons analysé les propriétés cinématiques, structurelles et chimiques des étoiles formées avant, pendant et après la collision. Ces collisions forment une galaxie ayant un disque pouvant être divisé en deux composantes. Ces deux composantes peuvent correspondre au disque mince et au disque épais d'une galaxie, et leur profil de luminosité peut être ajustées par une loi exponentielle. Les étoiles formées avant et pendant la collision ont une longueur d'échelle plus grande que les étoiles formées après la collision par 20% en moyenne. Du point de vue de la cinématique, les étoiles vieilles ont des dispersions en vitesse plus élevées et sont en retard sur les étoiles jeunes pour ce qui est de la vitesse de rotation. Le sursaut de formation d'étoiles associé à la collision permet d'enrichir rapidement le gaz en différent métaux. Les explosions de supernovae de type II qui ont lieu rapidement après la collision, étant donné la courte durée de vie des étoiles qui les produisent, enrichissent le milieu intergalactique en éléments a. Les supernovae de type la, ayant une distribution plus étendue dans le temps, permettent l'enrichissement en fer des deux populations associées aux composantes du disque, ce qui permet d'obtenir une population stellaire vieille ayant un rapport [a/Fe] supérieur à celui des étoiles jeunes et ce, même à des métallicités relativement élevées ([Fe/H] = —0.5). Ce résultat pourrait expliquer le rapport [a/Fe] élevé que l'on retrouve chez les étoiles du disque épais de la Voie lactée. / We employ GCD+, a N-body, smoothed particle hydrodynamic simulation, including star formation and a detailed treatment of chemical enrichment, to follow 14 gas-rich mergers that resuit almost ail in a galaxy with disk morphology. We trace the kinematic, structural, and chemical properties of stars formed before, during, and after the merger. We show that such merger produces two exponential disk components, with the older, hotter component having a scale length 20% larger than the later forming, cold disk. On a kinematical point of view the old stellar population clearly lags the rotation velocity of the young disk and hâve a higher rotational velocity dispersion. Rapid star formation during the mergers quickly enriches the protogalactic gas réservoir, resulting in high metallicities of the forming stars. Thèse stars form from gas largely polluted by Type II supernovae, which form rapidly in the merger-induced starburst. After the mergers, a thin disk forms from gas that has had time to be polluted by type la supernovae. This fact lead to an old stellar population with a higher [a/Fe] ratio than the young population at quite high metalicity ([Fe/H] = —0.5). We examine the proposai that increased star formation during gas-rich mergers may explain the high a-to-iron abundance ratios that exist in the relatively high-metallicity, thick-disk component of the Milky way.

L'étude porte sur les écoulements et les transferts en convection naturelle au sein d'empilements de produits alimentaires et plus particulièrement dans des réfrigérateurs ménagers. L'approche est à la fois expérimentale et numérique. Une première étude a été effectuée dans un empilement ordonné de sphères soumis à une convection à faible vitesse d'air (u<0.2 ms-1). L'intensité des échanges convectifs entre les objets et l'air ainsi qu'entre la paroi de la cavité et l'air a été caractérisée en fonction de la position et de la vitesse. Par ailleurs, les résultats ont montré que les autres modes de transfert (conduction et rayonnement entre les objets et avec les parois) sont du même ordre de grandeur que la convection. Des mesures de vitesse d'air (par vélocimétrie par images de particules) et de température ont été effectuées dans une enceinte fermée représentant une maquette d'un réfrigérateur ménager (à l'échelle 1) dans différentes configurations (maquette vide, maquette chargée d'empilement de produits inertes thermiquement et de produits échangeant de la chaleur). Des instationnarités dans le bas de l'enceinte vide ont été mis en évidence (Ra" 3´108). La présence d'empilements d'objets (même inerte thermiquement) modifie de façon importante les écoulements et les transferts. Les simulations numériques effectuées avec un logiciel de mécanique de fluide pour les différentes configurations étudiées sur la maquette et dans le cas d'un réfrigérateur réel ont donné des prédictions satisfaisantes en comparaison avec les mesures expérimentales. Mots clés: milieu macro-poreux, convection naturelle, coefficient de transfert convectif, conduction, rayonnement, réfrigérateur ménager, simulations numérique (CFD).