Nous nous intéressons dans cette thèse, à l’étude des méthodes itératives pour la résolutiond’équations matricielles de grande taille : Lyapunov, Sylvester, Riccati et Riccatinon symétrique.L’objectif est de chercher des méthodes itératives plus efficaces et plus rapides pour résoudreles équations matricielles de grande taille. Nous proposons des méthodes itérativesde type projection sur des sous espaces de Krylov par blocs Km(A, V ) = Image{V,AV, . . . ,Am−1V }, ou des sous espaces de Krylov étendus par blocs Kem(A, V ) = Image{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V } . Ces méthodes sont généralement plus efficaces et rapides pour les problèmes de grande dimension. Nous avons traité d'abord la résolution numérique des équations matricielles linéaires : Lyapunov, Sylvester, Stein. Nous avons proposé une nouvelle méthode itérative basée sur la minimisation de résidu MR et la projection sur des sous espaces de Krylov étendus par blocs Kem(A, V ). L'algorithme d'Arnoldi étendu par blocs permet de donner un problème de minimisation projeté de petite taille. Le problème de minimisation de taille réduit est résolu par différentes méthodes directes ou itératives. Nous avons présenté ainsi la méthode de minimisation de résidu basée sur l'approche global à la place de l'approche bloc. Nous projetons sur des sous espaces de Krylov étendus Global Kem(A, V ) = sev{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V }. Nous nous sommes intéressés en deuxième lieu à des équations matricielles non linéaires, et tout particulièrement l'équation matricielle de Riccati dans le cas continu et dans le cas non symétrique appliquée dans les problèmes de transport. Nous avons utilisé la méthode de Newtown et l'algorithme MINRES pour résoudre le problème de minimisation projeté. Enfin, nous avons proposé deux nouvelles méthodes itératives pour résoudre les équations de Riccati non symétriques de grande taille : la première basée sur l'algorithme d'Arnoldi étendu par bloc et la condition d'orthogonalité de Galerkin, la deuxième est de type Newton-Krylov, basée sur la méthode de Newton et la résolution d'une équation de Sylvester de grande taille par une méthode de type Krylov par blocs. Pour toutes ces méthodes, les approximations sont données sous la forme factorisée, ce qui nous permet d'économiser la place mémoire en programmation. Nous avons donné des exemples numériques qui montrent bien l'efficacité des méthodes proposées dans le cas de grandes tailles. / In this thesis, we focus in the studying of some iterative methods for solving large matrix equations such as Lyapunov, Sylvester, Riccati and nonsymmetric algebraic Riccati equation. We look for the most efficient and faster iterative methods for solving large matrix equations. We propose iterative methods such as projection on block Krylov subspaces Km(A, V ) = Range{V,AV, . . . ,Am−1V }, or block extended Krylov subspaces Kem(A, V ) = Range{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V }. These methods are generally most efficient and faster for large problems. We first treat the numerical solution of the following linear matrix equations : Lyapunov, Sylvester and Stein matrix equations. We have proposed a new iterative method based on Minimal Residual MR and projection on block extended Krylov subspaces Kem(A, V ). The extended block Arnoldi algorithm gives a projected minimization problem of small size. The reduced size of the minimization problem is solved by direct or iterative methods. We also introduced the Minimal Residual method based on the global approach instead of the block approach. We projected on the global extended Krylov subspace Kem(A, V ) = Span{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V }. Secondly, we focus on nonlinear matrix equations, especially the matrix Riccati equation in the continuous case and the nonsymmetric case applied in transportation problems. We used the Newton method and MINRES algorithm to solve the projected minimization problem. Finally, we proposed two new iterative methods for solving large nonsymmetric Riccati equation : the first based on the algorithm of extended block Arnoldi and Galerkin condition, the second type is Newton-Krylov, based on Newton’s method and the resolution of the large matrix Sylvester equation by using block Krylov method. For all these methods, approximations are given in low rank form, wich allow us to save memory space. We have given numerical examples that show the effectiveness of the methods proposed in the case of large sizes.

Les procédés de gazéification de bois à lits fixes étagés sont adaptés à la production d'électricité de petites puissances. Dans ces procédés, la pyrolyse est opérée dans un réacteur continu à lit fixe descendant. La particularité de ce type de réacteur est son fonctionnement autothermique. L'énergie nécessaire au chauffage, au séchage et à la pyrolyse est apportée par la combustion partielle du bois : on parle de “pyrolyse oxydante”. L'injection d'air par le haut du réacteur provoque la propagation d'une zone d'oxydation dans le milieu poreux à contre-courant des écoulements des solides et des gaz. Les travaux présentés dans ce manuscrit visent une meilleure description de cette étape du procédé. Le problème posé est de type multi-échelles. Ainsi, nous avons préalablement mené une étude à l'échelle de la particule isolée avant de s'intéresser au comportement global du lit fixe. A l'échelle de la particule, nous avons quantifié l'effet de l'oxygène et de la taille des particules sur la cinétique de la pyrolyse oxydante. Cette étude nous a guidés pour la mise en place d'un modèle cinétique de cette transformation. A l'échelle du lit fixe, la propagation de la zone d'oxydation a été caractérisée par des approches expérimentale et numérique, offrant ainsi une meilleure connaissance de cette étape du procédé étagé. / Wood Multi-staged gasification in a fixed bed reactor is suitable for small-scale electricity generation. In these processes, the pyrolysis is performed in a continuous downward fixed bed reactor. The main feature of this reactor is the autothermal operation. Energy for heating, drying and pyrolysis is supplied by partial combustion of wood, known as “oxidative pyrolysis”. The air introduced from the top of the reactor induces a combustion front that propagates countercurrent with the solids and gazes flows. The work presented in this document aimed to achieve a better description of this process. A multi-scale approach was considered. Therefore, we have firstly studied the behavior of an isolated particle before focusing on the overall fixed bed. At the particle scale, we have quantified the effect of oxygen and of particle size on the oxidative pyrolysis kinetics. This led us to the setup of a kinetic model for this transformation.At the fixed bed scale, the propagation of the combustion front was studied considering the experimental and numerical approaches, which provides a better understanding of this step of the wood staged gasifiers.

Les régions d’interface du Soleil de la photosphère à la chromosphère et au delà de la basse couronne ont été étudiées à partir des spectres éclairs. Les éclipses sont les plus adaptées à ce type d’observation, car l’occultation a lieu en dehors de l’atmosphère terrestre et sont exemptes de lumière parasite. Les images Extrême-UV des régions du limbe obtenues récemment dans l’espace sont analysés avec des modèles hydrostatiques à une dimension, comme les modèles VAL, mais cette méthode ne tient pas compte du phénomène d’émergence du champ magnétique, associé au réseau chromosphérique qui est responsable de: i) les spicules et le milieu interspiculaire, ii) les jets coronaux et macrospicules. Un saut de température de 0.01 à 1 MK est observé autour de 2 Mm d’altitude plus loin, et produit plus loin le flot du vent solaire permanent. Le processus de chauffage responsable du saut de température et la source du vent solaire ne sont pas encore compris. Dans cette thèse, nous traitons ces problèmes à partir de spectres éclairs récents réalisés avec les technologies actuelles de détecteurs CCD rapides, images d’éclipse en lumière blanche et des images EUV obtenues avec des instruments de missions spatiales. Nous illustrons les mécanismes des émissions des raies à faible potentiel de première ionisation (FIP) présents dans les basses couches de l’atmosphère solaire. Nous identifions plus précisément les raies à bas FIP à la fois dans les interfaces, à l’intérieur et en dehors des protubérances. Nous caractérisons en détail les enveloppes d’hélium et la région de l’interface solaire. Nous discutons de l'enrichissement de la couronne en éléments low FIP. / The solar interface region from the photosphere to the chromosphere and to the lower corona has been studied using flash spectra obtained during solar total eclipses. Eclipses are very favourable for this type of observation as the occultation takes place outside the Earth atmosphere and are free of parasitic scattered light. Independently, EUV filtergrams of the limb region obtained in space were analyzed using one dimensional hydrostatic models like the VAL models but this method ignores the ubiquitous magnetic field emergence phenomenon associated with the chromospheric network and responsible for: i) spicules and interspicular regions, ii) coronal jets and macrospicules. The components of the solar interface region are dynamical and different type of waves and magnetic reconnections are suggested to be at work. A jump of temperature from 0.01 to 1 MK is observed near the 2 Mm heights and higher, further producing a permanent solar wind flow. The heating processes responsible for this temperature jump and for the flow are not yet fully understood. In this thesis, we reconsider these problems on the basis of original, superior flash spectra which benefit from present technology such as CCD detectors, white light (W-L) eclipse images and new EUV images obtained with space-borne instruments. We illustrate the mechanisms of low First Ionisation Potential (FIP) emission lines present in the low layers of the solar atmosphere and interfaces. We identify more precisely low FIP lines both inside and nearby prominences. We characterize in detail the He shells and the solar interface region. We discuss the enrichment of low FIP elements in the corona.

Les convertisseurs de forte puissance sont des éléments critiques des futurs réseaux HVDC. À ce titre, leur fiabilité et leur endurance sont primordiales. La défaillance d’un composant se produit soit en circuit ouvert, ou en court-circuit. Le composant défaillant en circuit ouvert est inadmissible pour les convertisseurs utilisant une topologie de mise en série. En particulier, dans certaines applications HVDC, les modules doivent être conçus de telle sorte que lorsqu'une défaillance se produit, le module défaillant doit se comporter comme un court-circuit et supporter ainsi le courant nominal qui le traverse. Un tel comportement est appelé “défaillance en court-circuit” ou “failure-to-short-circuit”. Actuellement, tous les modules de puissance ayant un mode de défaillance en court-circuit disponibles dans le commerce utilisent des semi-conducteurs en silicium. Les potentialités des semi-conducteurs en carbure de silicium (SiC) poussent, aujourd’hui, les industriels et les chercheurs à mener des investigations pour développer des modules Fail-to-short à base des puces SiC. C’est dans ce contexte que se situe ce travail de thèse, visant à concevoir un module à base de puces SiC offrant un mode de défaillance de court-circuit. Pour cela nous présentons d’abord une étude de l’énergie de défaillance des puces SiC, afin de définir les plages d’activation du mécanisme Fail-to-short. Ensuite, nous démontrons la nécessité de remplacer les interconnexions classiques (fils de bonding) par des contacts massifs sur la puce. Enfin, une mise en œuvre est présentée au travers d’un module “demi pont” à deux transistors MOSFET. / The reliability and endurance of high power converters are paramount for future HVDC networks. Generally, module’s failure behavior can be classified as open-circuit failure and short-circuit failure. A module which fails to an open circuit is considered as fatal for applications requiring series connection. Especially, in some HVDC application, modules must be designed such that when a failure occurs, the failed module still able to carry the load current by the formation of a stable short circuit. Such operation is referred to as short circuit failure mode operation. Currently, all commercially available power modules which offer a short circuit failure mode use silicon semiconductors. The benefits of SiC semiconductors prompts today the manufacturers and researchers to carry out investigations to develop power modules with Fail-to-short-circuit capability based on SiC dies. This represents a real challenge to replace silicon power module for high voltage applications in the future. The work presented in this thesis aims to design a SiC power module with failure to short-circuit failure mode capability. The first challenge of the research work is to define the energy leading to the failure of the SiC dies in order to define the activation range of the Fail-to-short mechanism. Then, we demonstrate the need of replacing the conventional interconnections (wire bonds) by massive contacts. Finally, an implementation is presented through a "half bridge" module with two MOSFETs.

Le développement de l’énergie renouvelable loin des zones urbaines demande le transport d'une grande quantité d’énergie sur de longues distances. Le transport d’électricité en courant continu haute tension (HVDC) présente beaucoup d’avantages par rapport à celui en courant alternatif. Dans ce contexte il est nécessaire de développer des convertisseurs de puissance constitués par des composants électroniques très haute tension, 10 à 30 kV. Si les composants en silicium ne peuvent pas atteindre ces objectifs, le carbure de silicium (SiC) se positionne comme un matériau semiconducteur alternatif prometteur. Pour supporter des tensions élevées, une région de "drift", relativement large et peu dopée constitue le cœur du composant de puissance. En pratique l’obtention d’une tension de blocage effective dépend de plusieurs facteurs et surtout de la conception d'une terminaison de jonction adaptée. Cette thèse présente une méthode pour améliorer la tenue en tension des composants en SiC basée sur l’utilisation des terminaisons de jonctions : Deep Trench Termination. Cette méthode utilise une tranchée gravée profonde en périphérie du composant, remplie avec un matériau diélectrique pour supporter l'étalement des lignes équipotentielles. La conception de la diode avec cette terminaison a été faite par simulation TCAD, avec deux niveaux de tension 3 et 20 kV. Les travaux ont pris en compte les caractéristiques du matériau, les charges à l’interface de la tranchée et les limites technologiques pour la fabrication. Ce travail a abouti sur la fabrication de démonstrateurs et leur caractérisation pour valider notre conception. Lors de la réalisation de ces structures, la gravure plasma du SiC a été optimisée dans un bâti ICP de manière à obtenir une vitesse de gravure élevée et en conservant une qualité électronique de l'état des surfaces gravées. Cette qualité est confirmée par les résultats de caractérisation obtenus avec des tenues en tension proches de celle idéale. / The development of renewable energy away from urban areas requires the transmission of a large amount of energy over long distances. High Voltage Direct Current (HVDC) power transmission has many advantages over AC power transmission. In this context, it is necessary to develop power converters based on high voltage power electronic components, 10 to 30 kV. If silicon components cannot achieve these objectives, silicon carbide (SiC) is positioned as a promising alternative semiconductor material. To support high voltages, a drift region, relatively wide and lightly doped is the heart of the power component. In practice obtaining an effective blocking voltage depends on several factors and especially the design of a suitable junction termination. This thesis presents a method to improve the voltage withstand of SiC components based on the use of junction terminations: Deep Trench Termination. This method uses a trench deep etching around the periphery of the component, filled with a dielectric material to support the spreading of the equipotential lines. The design of the diode with this termination was done by TCAD simulation, with two voltage levels 3 and 20 kV. The work took into account the characteristics of the material, the interface charge of the trench and the technological limits for the fabrication. This work resulted in the fabrication of demonstrators and their characterization to validate the design. During the production of these structures, plasma etching of SiC has been optimized in an ICP reactor so as to obtain a high etching rate and maintaining an electronic quality of the state of etched surfaces. This quality is confirmed by the results of characterization obtained with blocking voltage close to the ideal one.

La géométrie moléculaire, également connue sous le nom de conformation, est la représentation la plus intrinsèque et la plus informative des molécules. Cependant, prédire des conformations stables à partir de graphes moléculaires reste un problème difficile et fondamental en chimie et en biologie computationnelles. Les méthodes expérimentales et computationelles traditionnelles sont généralement coûteuses et chronophages. Récemment, nous avons assisté à des progrès considérables dans l'utilisation de l'apprentissage automatique, en particulier des modèles génératifs, pour accélérer cette procédure. Cependant, les approches actuelles basées sur les données n'ont généralement pas la capacité de modéliser des distributions complexes et ne tiennent pas compte de caractéristiques géométriques importantes. Dans cette thèse, nous cherchons à construire des modèles génératifs basés sur des principes pour la génération de conformation moléculaire qui peuvent surmonter les problèmes ci-dessus. Plus précisément, nous avons proposé des modèles de diffusion basés sur les flux, sur l'énergie et de débruitage pour la génération de structures moléculaires. Cependant, il n'est pas trivial d'appliquer ces modèles à cette tâche où la vraisemblance des géométries devrait avoir la propriété importante d'invariance par rotation par de translation. Inspirés par les progrès récents de l'apprentissage des représentations géométriques, nous fournissons à la fois une justification théorique et une mise en œuvre pratique sur la manière d'imposer cette propriété aux modèles. Des expériences approfondies sur des jeux de données de référence démontrent l'efficacité de nos approches proposées par rapport aux méthodes de référence existantes. / Molecular geometry, also known as conformation, is the most intrinsic and informative representation of molecules. However, predicting stable conformations from molecular graphs remains a challenging and fundamental problem in computational chemistry and biology. Traditional experimental and computational methods are usually expensive and time-consuming. Recently, we have witnessed considerable progress in using machine learning, especially generative models, to accelerate this procedure. However, current data-driven approaches usually lack the capacity for modeling complex distributions and fail to take important geometric features into account. In this thesis, we seek to build principled generative models for molecular conformation generation that can overcome the above problems. Specifically, we proposed flow-based, energy-based, and denoising diffusion models for molecular structure generation. However, it's nontrivial to apply these models to this task where the likelihood of the geometries should have the important property of rotational and translation invariance. Inspired by the recent progress of geometric representation learning, we provide both theoretical justification and practical implementation about how to impose this property into the models. Extensive experiments on common benchmark datasets demonstrate the effectiveness of our proposed approaches over existing baseline methods.

Cette thèse expérimentale s’inscrit dans le domaine des sciences séparatives et se base sur la technique de SPLITT (SPLIT-flow Thin fractionation). Son objectif consiste en l’étude des mécanismes qui sont à l’origine de la séparation, en continu et sans membrane, d’objets de taille micrométrique dans des mini-séparateurs fluidiques (Step-SPLITT). Les expériences menées, en laboratoire et lors de vols paraboliques, ont révélé le couplage complexe comme l’influence des effets hydrodynamiques et du champ gravitationnel sur la migration transverse des espèces en écoulement. Des visualisations tridimensionnelles par holographie digitale ont corroboré nos résultats et dévoilé des comportements inattendus. Les capacités séparatives des Step-SPLITT ont rendu possible l’analyse et la séparation d’objets biologiques et biomimétiques. Enfin, cette étude complétée par une modélisation tridimensionnelle de l’écoulement nous a permis de mettre au point un nouveau prototype de séparateur. This experimental thesis belongs to the field of separative sciences and is based on the SPLITT technique (SPLIT-flow Thin fractionation). The objective is to study the mechanisms that are at the origin of continuous and membraneless separation of micron-size species in mini fluidic separators (Step-SPLITT). Experiments undertaken in laboratory and during parabolic flights revealed the complex coupling of the hydrodynamic effects and the gravitational field influencing the transverse migration of the flowing species. Three-dimensional visualizations performed by digital holography confirmed our results and disclosed unexpected behaviours. The separation capacities of Step-SPLITT made the analysis and the separation of biological and biomimetic species possible. In addition this study in conjunction with a three-dimensional flow modelling enabled us to develop a new prototype of separator.

Les argiles raides sont généralement rencontrées à grande profondeur (plus de 300 mètres). Ces matériaux se caractérisent par leur teneur en eau relativement faible. De nombreux projets industriels justifient l'intérêt récent porté à ces matériaux. La thèse traite en particulier d'argiles raides susceptibles de répondre à des sollicitations par des déformations élastiques, plastiques et viscoplastiques. C'est la prise en compte simultanée d'une rupture engendrant des déformations plastiques et des déformations différées de types viscoplastiques qui font l'originalité de cette thèse. Cette nouvelle approche constitue une amélioration indispensable dans la description du comportement mécanique d'argiles raides. Ce travail est présenté en trois parties : la première partie est consacrée à l'étude expérimentale et à la modélisation du comportement mécanique des argiles raides. Dans cette partie, nous décrivons les matériaux considérés, et les essais réalisés. Nous discutons les résultats et élaborons un modèle viscoplastique avec rupture. La deuxième partie, « Application du modèle VPR au calcul d'un tunnel de section circulaire », est consacrée à l'élaboration d'une solution originale semi analytique et d'un algorithme implanté dans GEOMEC91. La troisième partie, « Application du modèle VPR au calcul 2d axisymétrique des tunnels soutenus et étude paramétrique », montre très clairement l'influence du modèle VPR sur la convergence du tunnel au moment de la pose du soutènement, et par conséquent sur les sollicitations de ce dernier. Les résultats des calculs montrent une forte influence de la cohésion à court terme sur la convergence du tunnel.

Cette thèse expérimentale s’inscrit dans le domaine des sciences séparatives et se base sur la technique de SPLITT (SPLIT-flow Thin fractionation). Son objectif consiste en l’étude des mécanismes qui sont à l’origine de la séparation, en continu et sans membrane, d’objets de taille micrométrique dans des mini-séparateurs fluidiques (Step-SPLITT). Les expériences menées, en laboratoire et lors de vols paraboliques, ont révélé le couplage complexe comme l’influence des effets hydrodynamiques et du champ gravitationnel sur la migration transverse des espèces en écoulement. Des visualisations tridimensionnelles par holographie digitale ont corroboré nos résultats et dévoilé des comportements inattendus. Les capacités séparatives des Step-SPLITT ont rendu possible l’analyse et la séparation d’objets biologiques et biomimétiques. Enfin, cette étude complétée par une modélisation tridimensionnelle de l’écoulement nous a permis de mettre au point un nouveau prototype de séparateur.<p><p>This experimental thesis belongs to the field of separative sciences and is based on the SPLITT technique (SPLIT-flow Thin fractionation). The objective is to study the mechanisms that are at the origin of continuous and membraneless separation of micron-size species in mini fluidic separators (Step-SPLITT). Experiments undertaken in laboratory and during parabolic flights revealed the complex coupling of the hydrodynamic effects and the gravitational field influencing the transverse migration of the flowing species. Three-dimensional visualizations performed by digital holography confirmed our results and disclosed unexpected behaviours. The separation capacities of Step-SPLITT made the analysis and the separation of biological and biomimetic species possible. In addition this study in conjunction with a three-dimensional flow modelling enabled us to develop a new prototype of separator. / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

Cette recherche démontre que, compte tenu des tendances passées dans les comportements individuels face à la participation au marché du travail, l’augmentation de l’âge normal de la retraite n’est ni urgente ni nécessaire au Canada si son seul objectif est de pallier une augmentation de la durée de la retraite. Par ailleurs, si l’accueil d’un nombre toujours grandissant d’immigrants alimente la croissance de la main-d’œuvre, il entraîne aussi un déficit fiscal non négligeable. Celui-ci est toutefois en bonne partie lié au sous-emploi de cette main-d’œuvre potentielle; sous-emploi qui affecte aussi la population native, alimentant ainsi les risques de pénurie de main-d’œuvre. D’une part, et contrairement à l’opinion répandue selon laquelle l’allongement de la vie a été accompagné d’un rétrécissement du nombre d’années dédiées au travail, la durée de vie au travail en 2016 a augmenté de 4,96 ans et représente une part plus importante de l’espérance de vie, soit 3,55 points de pourcentage, comparé à son niveau en 1981. En outre, bien que l’évolution de la structure par âge de la population ait exercé une pression à la baisse (−11,7%) sur la croissance de l’offre agrégée de travail, les changements de comportement individuel quant à la participation au travail (19,3%) et les heures travaillées (5,8%) ont plus que compensé cet effet, entre 1981 et 2016. Au cours de cette période, les comportements individuels, notamment la participation des femmes au marché du travail, ont connu de profonds changements. Si ces changements étaient déjà connus et pleinement documentés, la quantification de leur contribution à l’offre individuelle et agrégée de travail le sont moins. Cette étude vient combler ce vide. Par exemple, et comme on pouvait s’y attendre, la participation des femmes au marché du travail a été la principale source de l’offre additionnelle de travail entre 1981 et 2016, contribuant 9,6 ans (contre 0,25 an pour les hommes) au changement de la durée de vie au travail. Au vu de ces résultats, le vieillissement de la population n’a pas encore entraîné une diminution de la durée de la vie active qui justifierait une augmentation de l’âge normal de la retraite au Canada. D’autre part, si l’immigration explique à elle seule 32,9% de l’augmentation de l’offre de travail entre 1981 et 2016, cette contribution n’est pas à coût nul. Par exemple, entre 1997 et 2015, les immigrants, comparés aux natifs du même âge, ont reçu 110$ de plus et contribué 3 520$ de moins en transferts publics. Ainsi, au même âge, l’immigrant moyen a reçu 3 640 en transferts nets de plus que le natif moyen. Toutefois, 85% de ce surplus provient des déséquilibres sur le marché du travail, dont le sous-emploi de cette main-d’œuvre potentielle. En effet, les résultats de cette recherche montrent que le sous-emploi représente un défi important dans un contexte de vieillissement de la population et de pénurie de main-d’œuvre appréhendée. S’il touche plus durement la population immigrante, il affecte aussi les natifs. Ainsi, en moyenne entre 1981 et 2016, 20,1% des travailleurs canadiens ont été en situation de sous-emploi, alors que le plein emploi aurait contribué à l’ajout de 1,5 million (±0,2) de travailleurs équivalent temps plein. Ces résultats découlent de l’application de plusieurs méthodes notamment la méthode de Sullivan (Sullivan, 1971), le modèle de changement continu (Horiuchi et al., 2008) ainsi que l’optimal matching et l’analyse des clusters sur un large éventail de données. Les sources de données incluent les recensements, les estimations de la population, les enquêtes sur la main-d’œuvre, la santé, et la consommation, au Canada et aux États-Unis. Cette recherche vient nuancer certaines des conséquences prétendues du vieillissement de la population sur la main-d’œuvre canadienne, tout en apportant un nouvel éclairage sur les solutions à mettre de l’avant afin de faire face à certains défis au cours des prochaines décennies. Si l’immigration peut ici jouer un rôle non négligeable, un meilleur arrimage entre l’offre et la demande de travail réduirait le sous-emploi parmi l’ensemble de la population canadienne, tout en améliorant leur niveau de vie. Les discussions sur l’augmentation de l’âge normal de la retraite ou des quotas d’immigration ne peuvent ignorer un tel constat. / This research demonstrates that, given past trends in individual behaviour in the labour market, increasing the normal retirement age is neither urgent nor necessary in Canada if its sole purpose is to offset an increase in the duration of retirement. Moreover, while the inflow of ever-increasing numbers of immigrants fuels labour force growth, it also creates a significant fiscal deficit. However, this deficit is largely related to the underemployment of this potential workforce, which also affects the native population, thus fuelling the risk of labour shortages. On the one hand, and contrary to the widespread opinion that population ageing has led to shrinking worklife for financing longer lifespan, Worklife Duration increased by 4,96 years and 3,55 percentage points of life expectancy over the last four decades. Furthermore, although the change in the age structure of the population has put a slight downward pressure (−11,7%) on the growth of aggregated labour supply, changes in individual behaviour regarding labour participation (19,3%) and worked hours (5,8%) have more than compensated for this effect between 1981 and 2016. During this period, individual behaviours, including women’s participation in the labour market, have undergone profound changes. While these changes were already known and fully documented, their contribution to individual and aggregate labour supply is less so. This study fills this gap. For example, and not surprisingly, women’s participation in the labour market was the main source of additional labour supply between 1981 and 2016, contributing 9,6 years (compared to 0,25 year for men) to the change in worklife duration. Given these results, population aging has not yet resulted in a decrease in working life which would justify an increase in the normal retirement age in Canada. On the other hand, while immigration alone accounts for 32,9% of the increased labour supply between 1981 and 2016, this contribution is at not zero cost. For example, between 1997 and 2015, immigrants, compared to natives at the same age, received $110 more and contributed $ 3 520 less, in public transfer. As a result, the average immigrant has received $3 640 in net transfer more than the average native. However, 85% of this deficit arises from the labour market imbalances such as the underemployment of this potential labour supply. Indeed, the results of this research show that underemployment represents a major challenge in the context of an aging population and a perceived labour shortage. While it affects the immigrant population most severely, it is also prevalent among native workers. Thus, on average, between 1981 and 2016, 20,1% of Canadian workers were underemployed, while full employment would have contributed an additional 1,5 million (±0,2) full-time equivalent workers. This study uses various methods, including the Sullivan method (Sullivan, 1971) and the model of continuous change (Horiuchi et al., 2008), as well as optimal matching and cluster analysis on a wide range of data. The data sources include censuses, population estimates, labour force, and health and consumer surveys in Canada and the United States. This research brings new evidence to the debates around the consequences of population aging on the Canadian workforce while shedding new light on the solutions to be put for facing the challenges in the coming decades. While immigration can play an important role, a better match between labour supply and demand would reduce underemployment among the Canadian population while improving their standard of living. Discussions about increasing the normal retirement age or immigration quotas cannot ignore such a finding.